[CP2K:5603] Re: Bulk Nickel (and possibly other newbie questions)

Ralph Koitz ralph... at gmail.com
Fri Aug 15 23:04:25 UTC 2014


I used the following
      BASIS_SET  DZVP-MOLOPT-SR-GTH
      POTENTIAL  GTH-PBE-q18

and set the multiplicity to n(atom)+1, basically s=1/2 for each atom in the
supercell.


On Fri, Aug 15, 2014 at 3:14 PM, Alex <nedo... at gmail.com> wrote:

> Hi Ralph,
>
> This is all very useful, thanks a bunch. I added cell multiplicity
> explicitly, as you suggested.
> Since I am not the owner of the installation (I'm running it on a remote
> cluster), I won't be able to compile anything myself, so I'll stick with
> what's there now (I think we have Debian binaries). I am looking at the
> reference you provided... For PBE/revPBE XC, what potential and basis did
> you use, exactly (assuming they are part of the package)? Also, did you use
> s=1/2 for the multiplicity in this case, i.e. the MULTIPLICITY param would
> be 2?..
>
> Thank you,
>
> Alex
>
>
> On Friday, August 15, 2014 3:36:47 PM UTC-6, Ralph wrote:
>
>> Hi,
>>
>> To use spin polarization, you should specify the keyword UKS in the DFT
>> section and define a MULTIPLICITY.
>> You may also want to turn on RELAX_MULTIPLICITY in that section.
>>
>> To use a supercell, you can either manually repeat your coordinates and
>> load the entire supercell as input coordinates, or specify
>> MULTIPLE_UNIT_CELL in both the CELL and TOPOLOGY sections of the SUBSYS.
>> A 6x6x6 supercell should be okay to start out with, but you may want to
>> check convergence.
>>
>> We have had decent experience with revPBE as an XC Functional for Ni [1].
>> BEEFvdW, a more recent addition to cp2k also gives very good bulk nickel,
>> but it's currently somewhat of a beta-version and you need to compile with
>> libxc to get it to work.
>>
>> Best,
>> Ralph
>>
>> [1] Gomez-Diaz et al, Theor Chim Acta, 2013
>>
>> On Fri, Aug 15, 2014 at 2:24 PM, Alex <ned... at gmail.com> wrote:
>>
>>>  Also, I just redid the calculation with
>>>
>>>     BASIS_SET DZV-GTH-PADE-q18
>>>    POTENTIAL GTH-PADE-q18
>>>
>>> Same energy within I believe five decimal places. As far as the correct
>>> choice of the XC-potential-basis combination, what would be better in my
>>> case? As I've said before, I am really new to this... Any reference to that
>>> effect would be great.
>>>
>>> Thank you.
>>>
>>> On Friday, August 15, 2014 3:05:37 PM UTC-6, Marcella Iannuzzi wrote:
>>>>
>>>> Hi
>>>>
>>>> Just few remarks.
>>>>
>>>> For the fcc bulk Ni energy, you need to construct a supercell, since
>>>> there is no k-point sampling,
>>>> and check the convergence of the result with system size.
>>>>
>>>> PADE is probably not an optimal choice for the XC functional, anyway,
>>>> you should use consistent potential and basis set,
>>>> i.e. for the same number of valence electrons (in your input, the PP is
>>>> for 10 v.e. and the BS for 18)
>>>>
>>>> If you don't specify in in the input, no spin polarisation is
>>>> considered.
>>>>
>>>> Regards
>>>>
>>>> Marcella
>>>>
>>>>
>>>> On Friday, August 15, 2014 10:19:59 PM UTC+2, Alex wrote:
>>>>>
>>>>> Hi all,
>>>>>
>>>>> I am very new to DFT calculations, let alone CP2k, so some level of
>>>>> idiocy on my part should be expected.
>>>>> As a simple test, I am trying to calculate the binding energy of a Ni
>>>>> atom in a bulk crystal. The relevant portion of the input shown below:
>>>>>
>>>>> ***
>>>>> &GLOBAL
>>>>>   PROJECT Ni_inp_test
>>>>>   RUN_TYPE ENERGY_FORCE
>>>>>   PRINT_LEVEL LOW
>>>>> &END GLOBAL
>>>>> &FORCE_EVAL
>>>>>   METHOD Quickstep
>>>>>   &SUBSYS
>>>>>     &KIND Ni
>>>>>       ELEMENT Ni
>>>>>       BASIS_SET DZV-GTH-PADE-q18
>>>>>       POTENTIAL GTH-PADE-q10
>>>>>     &END KIND
>>>>>     &CELL
>>>>>       A     1.765000    1.765000    0.000000
>>>>>       B     0.000000    1.765000    1.765000
>>>>>       C    1.765000    0.000000    1.765000
>>>>>       PERIODIC XYZ
>>>>>     &END CELL
>>>>>     &COORD
>>>>>       Ni    0.000000000    0.000000000    0.000000000
>>>>>     &END COORD
>>>>>   &END SUBSYS
>>>>>   &DFT
>>>>>     BASIS_SET_FILE_NAME  BASIS_SET
>>>>>     POTENTIAL_FILE_NAME  GTH_POTENTIALS
>>>>>     &QS
>>>>>       EPS_DEFAULT 1.0E-10
>>>>>     &END QS
>>>>>     &MGRID
>>>>>       NGRIDS 4
>>>>>       CUTOFF 300
>>>>>       REL_CUTOFF 60
>>>>>     &END MGRID
>>>>>     &XC
>>>>>       &XC_FUNCTIONAL PADE
>>>>>       &END XC_FUNCTIONAL
>>>>>     &END XC
>>>>>     &SCF
>>>>>       SCF_GUESS ATOMIC
>>>>>       EPS_SCF 1.0E-7
>>>>>       MAX_SCF 300
>>>>>       ADDED_MOS 10
>>>>>       &DIAGONALIZATION  ON
>>>>>         ALGORITHM STANDARD
>>>>>       &END DIAGONALIZATION
>>>>>       &MIXING  T
>>>>>         METHOD BROYDEN_MIXING
>>>>>         ALPHA 0.4
>>>>>         NBROYDEN 8
>>>>>       &END MIXING
>>>>>       &SMEAR ON
>>>>>         METHOD FERMI_DIRAC
>>>>>         ELECTRONIC_TEMPERATURE [K] 300
>>>>>       &END SMEAR
>>>>>     &END SCF
>>>>>   &END DFT
>>>>>   &PRINT
>>>>>     &FORCES ON
>>>>>     &END FORCES
>>>>>   &END PRINT
>>>>> &END FORCE_EVAL
>>>>>
>>>>> ***
>>>>>
>>>>> This yields a total energy of E1=-35.155 a.u. after convergence.
>>>>>
>>>>> Then I decided to calculate the "vacuum" energy of an isolated atom,
>>>>> input below:
>>>>>
>>>>> &GLOBAL
>>>>>   PROJECT Ni_inp_test
>>>>>   RUN_TYPE ENERGY_FORCE
>>>>>   PRINT_LEVEL LOW
>>>>> &END GLOBAL
>>>>> &FORCE_EVAL
>>>>>   METHOD Quickstep
>>>>>   &SUBSYS
>>>>>     &KIND Ni
>>>>>       ELEMENT Ni
>>>>>       BASIS_SET DZV-GTH-PADE-q18
>>>>>       POTENTIAL GTH-PADE-q10
>>>>>     &END KIND
>>>>>     &CELL
>>>>>       A     30.00000    0.000000    0.000000
>>>>>       B     0.000000    30.00000    0.000000
>>>>>       C    0.000000    0.000000    30.00000
>>>>>     &END CELL
>>>>>     &COORD
>>>>>       Ni    0.000000000    0.000000000    0.000000000
>>>>>     &END COORD
>>>>>   &END SUBSYS
>>>>>   &DFT
>>>>>     BASIS_SET_FILE_NAME  BASIS_SET
>>>>>     POTENTIAL_FILE_NAME  GTH_POTENTIALS
>>>>>     &QS
>>>>>       EPS_DEFAULT 1.0E-10
>>>>>     &END QS
>>>>>     &MGRID
>>>>>       NGRIDS 4
>>>>>       CUTOFF 300
>>>>>       REL_CUTOFF 60
>>>>>     &END MGRID
>>>>>     &XC
>>>>>       &XC_FUNCTIONAL PADE
>>>>>       &END XC_FUNCTIONAL
>>>>>     &END XC
>>>>>     &SCF
>>>>>       SCF_GUESS ATOMIC
>>>>>       EPS_SCF 1.0E-7
>>>>>       MAX_SCF 300
>>>>>       ADDED_MOS 10
>>>>>       &DIAGONALIZATION  ON
>>>>>         ALGORITHM STANDARD
>>>>>       &END DIAGONALIZATION
>>>>>       &MIXING  T
>>>>>         METHOD BROYDEN_MIXING
>>>>>         ALPHA 0.4
>>>>>         NBROYDEN 8
>>>>>       &END MIXING
>>>>>       &SMEAR ON
>>>>>         METHOD FERMI_DIRAC
>>>>>         ELECTRONIC_TEMPERATURE [K] 300
>>>>>       &END SMEAR
>>>>>     &END SCF
>>>>>   &END DFT
>>>>>   &PRINT
>>>>>     &FORCES ON
>>>>>     &END FORCES
>>>>>   &END PRINT
>>>>> &END FORCE_EVAL
>>>>>
>>>>> ***
>>>>>
>>>>> This also converges and yields a total energy E2=-34.555 a.u.
>>>>>
>>>>>
>>>>> Hence, my questions:
>>>>>
>>>>> 1. Is this even the correct way of calculating what I want, including
>>>>> the energy calculations, XC functional, and basis?
>>>>> 2. Should the spin properties be explicitly set in the input? There
>>>>> are none now.
>>>>> 3. Am I setting up the FCC lattice correctly (first input file)? My
>>>>> translation vectors are set by the ABC values, but I have no idea whether
>>>>> this is right.
>>>>> 4. If the first simulation yields the total energy of the system and
>>>>> the FCC lattice implies 12 nearest neighbors, then removing the center
>>>>> would change the total energy by (E1-E2)/6, which isn't the experimental
>>>>> -4.4 eV. Am I completely off track here? :)
>>>>>
>>>>> Thanks a lot!
>>>>>
>>>>  --
>>> You received this message because you are subscribed to the Google
>>> Groups "cp2k" group.
>>> To unsubscribe from this group and stop receiving emails from it, send
>>> an email to cp2k+... at googlegroups.com.
>>> To post to this group, send email to cp... at googlegroups.com.
>>>
>>> Visit this group at http://groups.google.com/group/cp2k.
>>> For more options, visit https://groups.google.com/d/optout.
>>>
>>
>>  --
> You received this message because you are subscribed to the Google Groups
> "cp2k" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to cp2k+uns... at googlegroups.com.
> To post to this group, send email to cp... at googlegroups.com.
> Visit this group at http://groups.google.com/group/cp2k.
> For more options, visit https://groups.google.com/d/optout.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20140815/f3e335b5/attachment.htm>


More information about the CP2K-user mailing list