[Pw_forum] something wrong in relax, BFGS
xywu
xywu at imr.ac.cn
Mon Apr 18 10:34:19 CEST 2005
Dear All,
There're something wrong in relax.
The follows occurs:
---------------
number of scf cycles = 117
number of bfgs steps = 62
energy old = -161.4284071964 ryd
energy new = -161.4284071569 ryd
CASE: energy_new > energy_old
new trust radius = 0.0000032530 bohr
WARNING : something is going wrong
---------------------
Then the scf cycle continues ,but the bfgs step does not.
Please give me some advice.Thank you in advance!
The input file is as follows:
&CONTROL
title = 'A--0.01-' ,
calculation = 'relax' ,
restart_mode = 'from_scratch' ,
outdir = '/root/tmp/' ,
pseudo_dir = '/root/software/pwscf/espresso/pseudo/' ,
prefix = 'A--0.01-' ,
nstep = 150 ,
tstress = .true. ,
tprnfor = .true. ,
tefield = .true. ,
/
&SYSTEM
ibrav = 6,
celldm(1) = 10.2,
celldm(3) = 5,
nat = 19,
ntyp = 2,
ecutwfc = 25 ,
ecutrho = 120 ,
occupations = 'smearing' ,
degauss = 0.005 ,
smearing = 'gaussian' ,
edir = 3 ,
emaxpos = 0.71 ,
eopreg = 0.01 ,
eamp = 0.01 ,
/
&ELECTRONS
electron_maxstep = 100,
conv_thr = 1.0d-8 ,
mixing_mode = 'plain' ,
mixing_beta = 0.1 ,
mixing_ndim = 16,
diagonalization = 'cg' ,
/
&IONS
ion_dynamics = 'bfgs' ,
upscale = 10 ,
potential_extrapolation = 'wfc2' ,
/
ATOMIC_SPECIES
Si 28.08600 Si.pz-vbc.UPF
N 14.00670 N.pz-vbc.UPF
ATOMIC_POSITIONS angstrom
Si 0.000000000 0.000000000 5.397600000
Si 2.698800000 2.698800000 5.397600000
Si 1.349400000 4.048200000 4.048200000
Si 4.048200000 1.349400000 4.048200000
Si 1.349400000 4.048200000 9.445800000
Si 4.048200000 1.349400000 9.445800000
Si 0.000000000 2.698800000 2.698800000 0 0 0
Si 2.698800000 0.000000000 2.698800000 0 0 0
Si 0.000000000 2.698800000 8.096400000
Si 2.698800000 0.000000000 8.096400000
Si 1.349400000 1.349400000 1.349400000 0 0 0
Si 4.048200000 4.048200000 1.349400000 0 0 0
Si 1.349400000 1.349400000 6.747000000
Si 4.048200000 4.048200000 6.747000000
Si 0.000000000 0.000000000 0.000000000 0 0 0
Si 2.698800000 2.698800000 0.000000000 0 0 0
Si 0.000000000 0.000000000 10.795200000
Si 2.698800000 2.698800000 10.795200000
N 2.698800000 2.698800000 16.795200000 1 1 0
K_POINTS automatic
4 4 1 1 1 1
xywu
xywu at imr.ac.cn
More information about the Pw_forum
mailing list