[Pw_forum] Fw: relax is too slow
vega
vegalew at hotmail.com
Fri Sep 12 11:26:01 CEST 2008
>I think it is the too much k-points that caused the problem .when I set the k-point to Gamma ,it iterated very quickly.but I don't know whether it is right or wrong that I set it to Gamma.
I don't know either. I think you could compare the Gamma results and 2X2X2 or 4X4X4 if necessary. If the energy different and coordinate is nearly the same. you may could use a small K-mesh to save computional cost.
As far as I tested, the energy would be within 1% for large supercells, but the coordinates of several atoms may differ 200%. So I have the same question like you.
Could anyone help us, please~ what is the benchmark for judging whether the k-mesh is adequate? the most important thing is the energy or the coordinates of atoms? If the energy difference is with 1% and 80% of atoms or more in the system is nearly at the same position, but coordinate some of atoms very different for the two case. Could we trust the results calculated by relative small number of k points?
vega
=================================================================================
Vega Lew (weijia liu)
PH.D Candidate in Chemical Engineering
State Key Laboratory of Materials-oriented Chemical Engineering
College of Chemistry and Chemical Engineering
Nanjing University of Technology, 210009, Nanjing, Jiangsu, China
From: wangqj1
Sent: Friday, September 12, 2008 3:34 PM
To: vegalew at hotmail.com
Subject: relax is too slow
Thank you very much for your advice !
My input file as follows:
&CONTROL
title = ZnOrelax ,
calculation = 'relax' ,
restart_mode = 'from_scratch' ,
pseudo_dir = '/root/tool/espresso-3.2.3/pseudo/' ,
prefix = ZnO ,
forc_conv_thr = 1.D-4 ,
tstress = .true. ,
tprnfor = .true. ,
dt = 100 ,
/
&SYSTEM
ibrav = 14,
celldm(1) = 18.6488,
celldm(2) = 0.666667,
celldm(3) = 1.0691,
celldm(4) = 0,
celldm(5) = 0,
celldm(6) = -0.5,
nat = 48,
ntyp = 2,
ecutwfc = 40 ,
ecutrho = 200 ,
nosym = .true. ,
nbnd = 220,
occupations = 'smearing' ,
degauss = 0.002 ,
smearing = 'methfessel-paxton' ,
nspin = 1 ,
lda_plus_u = .false. ,
/
&ELECTRONS
startingpot = 'atomic' ,
startingwfc = 'atomic' ,
mixing_mode = 'plain' ,
mixing_beta = 0.7 ,
diagonalization = 'cg' ,
/
&IONS
ion_dynamics = 'damp' ,
phase_space = 'full' ,
pot_extrapolation = 'second_order' ,
wfc_extrapolation = 'second_order' ,
/
ATOMIC_SPECIES
Zn 65.00000 Zn.pbe-van.UPF
O 16.00000 O.pbe-rrkjus.UPF
ATOMIC_POSITIONS angstrom
Zn 0.000000000 1.899193710 0.000000000 1 1 1
O 0.000000000 1.899193710 2.018253260 1 1 1
Zn 1.644750000 0.949596855 2.637550000 1 1 1
O 1.644750000 0.949596855 4.655803260 1 1 1
Zn 3.289500000 1.899193710 0.000000000 1 1 1
O 3.289500000 1.899193710 2.018253260 1 1 1
Zn 4.934250000 0.949596855 2.637550000 0 0 0
O 4.934250000 0.949596855 4.655803260 0 0 0
Zn 6.579000000 1.899193710 0.000000000 0 0 0
O 6.579000000 1.899193710 2.018253260 0 0 0
Zn 8.223750000 0.949596855 2.637550000 0 0 0
O 8.223750000 0.949596855 4.655803260 0 0 0
Zn -1.644750000 4.747984276 0.000000000 1 1 1
O -1.644750000 4.747984276 2.018253260 1 1 1
Zn 0.000000000 3.798387421 2.637550000 1 1 1
O 0.000000000 3.798387421 4.655803260 1 1 1
Zn 1.644750000 4.747984276 0.000000000 1 1 1
O 1.644750000 4.747984276 2.018253260 1 1 1
Zn 3.289500000 3.798387421 2.637550000 0 0 0
O 3.289500000 3.798387421 4.655803260 0 0 0
Zn 4.934250000 4.747984276 0.000000000 0 0 0
O 4.934250000 4.747984276 2.018253260 0 0 0
Zn 6.579000000 3.798387421 2.637550000 0 0 0
O 6.579000000 3.798387421 4.655803260 0 0 0
Zn 0.000000000 1.899193710 5.275100000 1 1 1
O 0.000000000 1.899193710 7.293353260 1 1 1
Zn 1.644750000 0.949596855 7.912650000 1 1 1
O 1.644750000 0.949596855 9.930903260 1 1 1
Zn 3.289500000 1.899193710 5.275100000 1 1 1
O 3.289500000 1.899193710 7.293353260 1 1 1
Zn 4.934250000 0.949596855 7.912650000 0 0 0
O 4.934250000 0.949596855 9.930903260 0 0 0
Zn 6.579000000 1.899193710 5.275100000 0 0 0
O 6.579000000 1.899193710 7.293353260 0 0 0
Zn 8.223750000 0.949596855 7.912650000 0 0 0
O 8.223750000 0.949596855 9.930903260 0 0 0
Zn -1.644750000 4.747984276 5.275100000 1 1 1
O -1.644750000 4.747984276 7.293353260 1 1 1
Zn 0.000000000 3.798387421 7.912650000 1 1 1
O 0.000000000 3.798387421 9.930903260 1 1 1
Zn 1.644750000 4.747984276 5.275100000 1 1 1
O 1.644750000 4.747984276 7.293353260 1 1 1
Zn 3.289500000 3.798387421 7.912650000 0 0 0
O 3.289500000 3.798387421 9.930903260 0 0 0
Zn 4.934250000 4.747984276 5.275100000 0 0 0
O 4.934250000 4.747984276 7.293353260 0 0 0
Zn 6.579000000 3.798387421 7.912650000 0 0 0
O 6.579000000 3.798387421 9.930903260 0 0 0
K_POINTS automatic
4 4 4 0 0 0
According your advice I found the CPU user% ratio is about 98% and the swap is about 16%.I think it is the too much k-points that caused the problem .when I set the k-point to Gamma ,it iterated very quickly.but I don't know whether it is right or wrong that I set it to Gamma.
Sincerely
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080912/e29e69ea/attachment-0001.htm
More information about the Pw_forum
mailing list