[Pw_forum] questions on vc-relax
Guoying Gao
gaoguoying at gmail.com
Sun Dec 3 07:44:56 CET 2006
Dear Axel,
I'm sorry for my ambiguous esperssions. I mean that in the process
of relaxing, either "etot_conv_thr" or "forc_con_thr" ( or both )
does not converge. I have attached my input and output files(in the accessory).
Another question:
>>first of all, you have to make sure that your stress tensor is
>>reasonably converged with respect to the basis set
How to judge that stress tensor converges with respect to the basis set.
The following is my input files:
&control
calculation='vc-relax',
restart_mode='from_scratch',
dt=30,
prefix='n2',
pseudo_dir = '/home/cpmd/yanming/ggy/pseudo/',
outdir='/home/cpmd/yanming/ggy/tmp/tmp0/',
tstress=.t.,
tprnfor=.t.,
nstep= 100
/
&system
ibrav= 8, celldm(1) =7.85937,celldm(2)=0.66482327482567924981,
celldm(3)=1.21158932435681654243,
nat= 8, ntyp= 1,ecutwfc =30, ecutrho=300, nosym = .t.
/
&electrons
mixing_beta = 0.7
conv_thr = 1.0d-8
/
&IONS
ion_dynamics = 'damp'
upscale = 20
/
&CELL
cell_dynamics = 'damp-pr'
press =800
/
ATOMIC_SPECIES
N 14.0067 N.pbe-van_ak.UPF
ATOMIC_POSITIONS {crystal}
N 0.39750000 0.86660000 0.19890000
N 0.92010000 0.64090000 0.05920000
N 0.60250000 0.13340000 0.69890000
N 0.07990000 0.35910000 0.55920000
N 0.60250000 0.86660000 0.30110000
N 0.07990000 0.64090000 0.44080000
N 0.39750000 0.13340000 0.80110000
N 0.92010000 0.35910000 -0.05920000
K_POINTS { automatic }
5 8 4 0 0 0
Thank you again!
gao
More information about the Pw_forum
mailing list