[CP2K-user] Cell Optimization for Copper.
Travis
polla... at gmail.com
Fri Aug 16 15:25:45 UTC 2019
Hi,
It's not too bad but can probably be better. You can check for convergence
of the lattice parameter by increasing k-points, changing basis, changing
cutoffs, tightening convergence criteria, etc. Different potentials will
give different results - and will vary in quality across the periodic
table. The warnings are due to features missing in CP2K when using Fermi
smearing. They are safe to ignore. Converting a CIF file to XYZ loses
symmetry operations, so you need to fill the unit cell. I usually convert
from CIF to XYZ because sometimes CP2K struggles to read in CIFs and quits
out.
-T
On Friday, August 16, 2019 at 7:48:13 AM UTC-3, Hasan Al-Mahayni wrote:
>
> I ran the files you posted. I get 3.651 Angstroms as a result in few
> iterations, which is still far from 3.621. I don't know if that's ok, but I
> get 11 warnings in my output file that I attached. Also, what is the
> meaning to the fact that the XYZ file has 4 atoms of copper only?
>
> Thanks for your help,
> Hasan.
>
> On Wednesday, August 14, 2019 at 8:25:15 PM UTC-4, Hasan Al-Mahayni wrote:
>>
>> Hello all,
>>
>> I am fairly new to cp2k so I believe this issue is easily resolvable. My
>> end goal is calculating adsorption energy of different molecules on copper
>> surfaces (slab). However, the first step I need to do is optimize the cell
>> parameters of copper(cell size). I know from materials project site that it
>> is 3.62126 Angstroms. That is the configuration I start with and is
>> attached by the name of "copper.xyz". The problem I am having is that
>> when I run cell optimization on this (input file is "copper.inp"), I get
>> the wrong results as the cell size shrinks (I get 3.4A). I think this is
>> because of the vacuum between the copper cell and the simulation cell. My
>> question is, what is the right way to set up the input file and especially
>> the &CELL part. I started off by making the A B C cell size 10 A bigger
>> than 3.62A, but that gives me the wrong result. If the A B C cell size is
>> too small, I get an error saying geometry is wrong. Also, maybe GEO_OPT
>> would be a more appropriate run type... I tried making supercell 2x2x2,
>> 3x3x3, and 4x4x4 while keeping ABC cell size 10 A bigger every time, and
>> the result gets closer to 3.62 but this method is not ideal as the system
>> gets bigger and bigger and the calculation is more computationally
>> expensive. Any ideas?
>>
>> Thanks,
>>
>> Hasan.
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20190816/3a708fce/attachment.htm>
More information about the CP2K-user
mailing list