[CPMD-list] convergence problems with linear response EPR calculation

Silviu Zilberman silviu at Princeton.EDU
Wed May 17 10:37:41 CEST 2006


Hi,

I've been unsuccessful in converging an EPR linear response calculation. 
Below is the input file. The calculation stops with an error message

 *** TOO MANY RESETS. EXITING.
I tried tweaking the parameters in the &resp section:

&RESP
  EPR
  CG-ANALYTIC
  20
  CG-FACTOR
  0.5
  POLAK
&END


I also tried CG-ANALYTIC=10 and CG-FACTOR=0.7.

Are there any other ideas/tricks to make the computation converge?

Thanks, Silviu.

--------------------------------------------------------------------------------------------------

INPUT FILE:

&CPMD
 LINEAR RESPONSE
 RESTART ALL LATEST
 CONJUGATE GRADIENT ELECTRONS
 CONVERGENCE ORBITALS
  3.0d-4
  MAXSTEP
   50000
  TRAJECTORY XYZ
  REAL SPACE WF KEEP
  FILEPATH
  /pvfs/silviu/EPR_PDT_active
  LSD
&END

&SYSTEM
  ANGSTROM
  SYMMETRY
  1
  CELL
  13.229425  1.0  1.0  0  0  0
  CUTOFF
  120.
  CHARGE
  -2
  MULTIPLICITY
  2
&END

&RESP
  EPR
  CG-ANALYTIC
  20
  CG-FACTOR
  0.5
  POLAK
&END

&ATOMS

*Fe-q16
    LMAX=D
    2
        4.59378 6.01722 6.64393
        7.25407 5.60785 6.76516

*S-q6
   LMAX=P
   3
        6.02759 6.9518 8.24564
        6.1248 6.84295 5.07668
        9.17715 7.06325 6.72235
*O-q6
   LMAX=P
   4
        8.4227 4.13464 8.97638
        2.99429 4.85191 8.78731
        3.14056 4.55461 4.58191
        5.64473 3.18001 6.76378
*N-q5
   LMAX=P
   2
        8.58601 3.92803 4.56286
        2.46491 8.27062 6.4422
*C-q4
   LMAX=P
   10
        5.85625 9.29547 6.58413
        6.44051 8.63744 5.3385
        6.3696 8.7282 7.90236
        8.0854 4.58616 5.40449
        7.98264 4.74203 8.07122
        3.66034 5.30188 7.93443
        3.31794 7.46029 6.51862
        10.6484 5.98372 6.66049
        5.92112 4.34117 6.73045
        3.75089 5.14112 5.39089

*H-q1
   LMAX=S
   9
        7.4678 8.82514 7.97528
        5.9162 9.27653 8.74375
        6.12808 10.3697 6.55365
        7.54167 8.73005 5.32463
        6.0404 9.1328 4.43948
        11.5255 6.6054 6.42739
        10.8232 5.48561 7.62698
        10.5277 5.22004 5.87937
        4.76028 9.23949 6.55651
&END

&DFT
  NEWCODE
  FUNCTIONAL PBE
&END




More information about the CPMD-list mailing list