[Pw_forum] phonon calculation stops at a particular mode
lung Fermin
ferminlung at gmail.com
Mon Dec 7 04:57:59 CET 2015
Dear QE community,
Recently, I am trying to calculate the phonon dispersion in LiOsO3 (with
R-3c space group). The program always stop at a mode which the scf cycle
does not converge after 100 iteration. I have tried reducing alpha_mix and
increasing the q-point density but the problem persists. The following are
the inputs and output. Any help or suggestions would be much appreciated.
::::::::::::::
scf.in <http://lioso3.scf.in/>
::::::::::::::
&control
calculation='scf',
restart_mode='from_scratch',
prefix='LiOsO3',
outdir = './tmp/',
pseudo_dir = '/pslibrary.1.0.0/pbe/PSEUDOPOTENTIALS/',
max_seconds = 172000
/
&system
nat = 10
ntyp = 3
ibrav = 5
celldm(1) = 10.0031
celldm(4) = 0.530117
ecutwfc =60
ecutrho=480
/
&electrons
! mixing_beta=0.3D0,
! mixing_ndim=10
conv_thr=1.D-10
/
&ions
/
&cell
/
ATOMIC_SPECIES
Li 6.941 Li.pbe-sl-rrkjus_psl.1.0.0.UPF
Os 190.23 Os.pbe-spn-rrkjus_psl.1.0.0.UPF
O 15.999 O.pbe-n-rrkjus_psl.1.0.0.UPF
ATOMIC_POSITIONS crystal
Li 0.250000000 0.250000000 0.250000000
Li 0.750000000 0.750000000 0.750000000
Os 0.000000000 0.000000000 0.000000000
Os 0.500000000 0.500000000 0.500000000
O 0.623133124 0.876866876 0.250000000
O 0.376866876 0.123133124 0.750000000
O 0.876866876 0.250000000 0.623133124
O 0.123133124 0.750000000 0.376866876
O 0.250000000 0.623133124 0.876866876
O 0.750000000 0.376866876 0.123133124
K_POINTS AUTOMATIC
15 15 15 1 1 1
::::::::::::::
phGrid.in <http://lioso3.phgrid.in/>
::::::::::::::
&inputph
prefix='LiOsO3',
outdir='./tmp/'
ldisp=.true.
fildyn='LiOsO3.dyn',
nq1=4,
nq2=4,
nq3=4,
alpha_mix(1)=0.05,
tr2_ph=1.0d-10,
/
::::::::::::::
(Some lines in the output)
iter # 96 total cpu time : 75159.6 secs av.it.: 38.6
thresh= 1.755E-06 alpha_mix = 0.050 |ddv_scf|^2 = 2.806E-10
iter # 97 total cpu time : 75686.1 secs av.it.: 40.8
thresh= 1.675E-06 alpha_mix = 0.050 |ddv_scf|^2 = 1.889E-10
iter # 98 total cpu time : 76134.3 secs av.it.: 33.7
thresh= 1.374E-06 alpha_mix = 0.050 |ddv_scf|^2 = 1.803E-10
iter # 99 total cpu time : 76695.2 secs av.it.: 41.8
thresh= 1.343E-06 alpha_mix = 0.050 |ddv_scf|^2 = 2.824E-10
iter # 100 total cpu time : 77231.0 secs av.it.: 41.2
thresh= 1.680E-06 alpha_mix = 0.050 |ddv_scf|^2 = 3.136E-10
End of self-consistent calculation
No convergence has been achieved
=========
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20151207/3d114f3e/attachment.html>
More information about the users
mailing list