[CPMD-list] question about CUTOFF

liu liuxiandong xiandongliu at gmail.com
Fri Feb 22 17:13:43 CET 2008


Dear CPMD users,

  I want to optimize the geometry of a silicate material. The
following is part of the input,
&CPMD	
 RESTART WAVEFUNCTIONS COORDINATES HESSIAN	
 OPTIMIZE GEOMETRY XYZ
 HESSIAN UNIT

 CONVERGENCE ORBITALS
 5.0D-6
 CONVERGENCE GEOMETRY
 5.0D-4

 PCG MINIMIZE
 TIMESTEP
 10
 LBFGS NREM
 40
 LSD

&END
...
&SYSTEM
  CUTOFF
   70.0
....
&END

&ATOMS

ISOTOPES
...
*O_MT_BLYP.psp KLEINMAN-BYLANDER
 LMAX=P
*H_MT_BLYP.psp KLEINMAN-BYLANDER
 LMAX=S
*Si_MT_BLYP.psp KLEINMAN-BYLANDER
 LMAX=D
*Al_MT_GIA_BLYP.psp KLEINMAN-BYLANDER
 LMAX=D LOC=0
*Fe_MT_BLYP_NLCC.psp NLCC KLEINMAN-BYLANDER
 LMAX=D LOC=S
*Na_SEMI_BLYP.psp KLEINMAN-BYLANDER
 LMAX=P
...............
The total energy of the relaxed structure is  -983.949852 at CUTOFF 70Ry.
Then I increased CUTOFF to 75 Ry and performed single-point
calculation with this optimized structure,  the total energy decreased
to about -984.845447.

Since such LSD calculation is very time demanding, is it necessary to
optimize the structure with higher CUTOFF ?
Thanks.

Xiandong Liu


More information about the CPMD-list mailing list