[CP2K:6953] Use the LSD option for an odd number of electrons
Zhendong Guo
zhendon... at gmail.com
Sat Sep 5 12:16:45 UTC 2015
Dear Alfred
You make a very very small mistake. Nowadays, the total number of
electrons in your system is even. So should add "LSD" command into your
input file. Please Put "LSD" after "CHARGE -1". Everything will be fine.
Best regards.
zhendong
On Sat, Sep 5, 2015 at 12:53 PM, Kit Tang <scott... at gmail.com> wrote:
> Dear Expert,
>
> I'm a newbie to CP2K, just started to learn about MD last week.
>
> I'm constructing an input for radical anion CO2. (single point calculation)
> however, i got this error "Use the LSD option for an odd number of
> electrons"
>
> &GLOBAL
> PROJECT co2-6w
> RUN_TYPE ENERGY_FORCE
> PRINT_LEVEL LOW
> &END GLOBAL
> &FORCE_EVAL
> METHOD Quickstep
> &SUBSYS
> &KIND O
> ELEMENT O
> BASIS_SET TZV2P-MOLOPT-GTH
> POTENTIAL GTH-PBE-q6
> &END KIND
> &KIND H
> ELEMENT H
> BASIS_SET TZV2P-MOLOPT-GTH
> POTENTIAL GTH-PBE-q1
> &END KIND
> &KIND C
> ELEMENT C
> BASIS_SET TZV2P-MOLOPT-GTH
> POTENTIAL GTH-PBE-q4
> &END KIND
> &CELL
> ABC 8.0 8.0 8.0
> &END CELL
> &COORD
> O 2.29129800 -1.18830800 -0.47113100
> O 0.36340100 -0.54189000 -2.11624300
> O 0.11560000 2.06897400 -1.26425100
> O 0.18463600 -1.95825500 1.32803500
> O 2.04036400 1.44194800 0.70523000
> O -0.12580300 0.67826100 2.29573900
> O -1.88449400 -1.21917500 -0.52771200
> O -2.10558400 0.92308900 0.18316900
> H 1.60817400 -0.97433800 -1.16834000
> H 2.98547000 -1.68509900 -0.90611900
> H -0.44547000 -0.87949700 -1.68917300
> H 0.29416800 0.42859200 -1.99445000
> H -0.68834500 1.86302000 -0.75927800
> H 0.83356900 2.02229400 -0.60641600
> H 0.98596400 -1.83704500 0.79754000
> H -0.54600000 -1.86021900 0.69414700
> H 1.36730700 1.25002700 1.38629600
> H 2.31400700 0.56982900 0.38640600
> H -0.85915200 0.90354400 1.70135400
> H -0.02622800 -0.28151100 2.16433600
> C -2.47646900 -0.19279100 -0.17816500
> &END COORD
> &END SUBSYS
> &DFT
> BASIS_SET_FILE_NAME /share/apps/cp2k/cp2k-2.6.1/data/BASIS_MOLOPT
> POTENTIAL_FILE_NAME /share/apps/cp2k/cp2k-2.6.1/data/GTH_POTENTIALS
> MULTIPLICITY 2
> CHARGE -1
> &QS
> EPS_DEFAULT 1.0E-10
> &END QS
> &MGRID
> NGRIDS 4
> CUTOFF 280
> REL_CUTOFF 60
> &END MGRID
> &XC
> &XC_FUNCTIONAL PBE
> &END XC_FUNCTIONAL
> &END XC
> &SCF
> SCF_GUESS ATOMIC
> EPS_SCF 1.0E-7
> MAX_SCF 300
> &DIAGONALIZATION ON
> ALGORITHM STANDARD
> &END DIAGONALIZATION
> &MIXING T
> METHOD BROYDEN_MIXING
> ALPHA 0.4
> NBROYDEN 8
> &END MIXING
> &END SCF
> &END DFT
> &PRINT
> &FORCES ON
> &END FORCES
> &END PRINT
> &END FORCE_EVAL
>
> Im not sure if i did it correctly. Please correct me if i were wrong.
>
> Rgds
> Alfred
>
> --
> 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.
>
--
Zhendong Guo
Phone: +41786943316
Email: zhendon... at gmail.com
zhendo... at epfl.ch
-------------------------------------------------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150905/9d0b2dc9/attachment.htm>
More information about the CP2K-user
mailing list