[QE-users] DFT relaxation steps
Giuseppe Mattioli
giuseppe.mattioli at ism.cnr.it
Wed Dec 1 08:55:31 CET 2021
Dear Sunil Kumar
Check the structure first. I opened your file with xcrysden and even
if it is not totally weird, it is a bit strange (Ni-Y distances of 2.2
A are likely wrong...). Moreover, are you sure you want to use the LDA
(PZ) functional? The electronic structure of a TM oxide will likely be
strongly affected by delocalization errors. Finally, I expect that you
norm-conserving (so no need to specify ecutrho, it has a fixed value
of 4*ecutwfc) mt-fhi pseudopotentials are quite hard, with ecutwfc
maybe as high as 120-140 Ry. In difficult cases (but check the
structure first!) you may lower mixing_beta down to 0.01.
HTH
Giuseppe
Quoting "Dr. SUNIL KUMAR" <suniliitd14 at gmail.com>:
> Dear QE developers and Users.
>
> We are doing DFT simulation of Monoclinic YNiO3 with and without oxygen
> vacancy. During the vc-relax and relax relaxation step, energy does not
> converge and simulation ends with error msg " convergence NOT achieved
> after 200 iterations: stopping" . I have tried many different values of
> ecutwfc and mixing_beta but unable to find convergence. Your suggestions
> will be highly appreciated.
> Following is the script for QE.
>
> Thank you
>
> &CONTROL
> calculation = "relax"
> prefix = "espresso"
> wf_collect = .TRUE.
> forc_conv_thr = 1.00000e-03
> max_seconds = 1.72800e+05
> nstep = 100
> pseudo_dir = "/home/pseudopot"
> outdir = "./outdir"
> wfcdir = "./wfcdir"
> /
>
> &SYSTEM
> a = 5.17932e+00
> b = 5.51529e+00
> c = 7.41657e+00
> cosac = -1.34914e-03
> degauss = 1.00000e-02
> ecutrho = 2.25000e+02
> ecutwfc = 2.50000e+01
> ibrav = -12
> nat = 19
> nspin = 2
> ntyp = 3
> occupations = "smearing"
> smearing = "gaussian"
> starting_magnetization(1) = 2.00000e-01
> starting_magnetization(2) = 0.00000e+00
> starting_magnetization(3) = 2.00000e-01
> /
>
> &ELECTRONS
> conv_thr = 1.00000e-06
> electron_maxstep = 200
> mixing_beta = 4.00000e-01
> startingpot = "atomic"
> startingwfc = "atomic+random"
> /
>
> &IONS
> ion_dynamics = "bfgs"
> /
>
> &CELL
> /
>
> K_POINTS {automatic}
> 2 2 2 0 0 0
>
> ATOMIC_SPECIES
> Ni 58.69340 Ni.pz-mt_fhi.UPF
> O 15.99940 O.pz-mt_fhi.UPF
> Y 88.90585 Y.pz-mt_fhi.UPF
>
> ATOMIC_POSITIONS {angstrom}
> Ni 2.589660 0.000000 0.000000
> Ni 0.000000 2.757645 0.000000
> Ni -0.005000 0.000000 3.708280
> Ni 2.584660 2.757645 3.708280
> O 1.796655 1.696862 0.356885
> O 1.238791 4.384435 0.375842
> O 0.756099 2.606316 1.812340
> O 4.418221 5.363961 1.895940
> O 3.935529 1.626790 3.332439
> O 3.377665 4.454507 3.351395
> O 1.233790 3.888500 4.084122
> O 0.751099 0.151329 5.520620
> O 4.413221 2.908974 5.604220
> O 3.930529 1.130855 7.040719
> O 3.372665 3.818428 7.059675
> Y 3.782084 3.157134 1.851618
> Y 1.392236 0.399489 1.856662
> Y 3.777084 5.115801 5.559899
> Y 1.387236 2.358156 5.564942
>
>
> Dr. Sunil Kumar
> Ph.D (Chemical Engg. IIT Delhi)
> M.Tech (Chemical Engg. IIT Delhi)
> B.Tech (Chemical Engg. IET-CSJMU Kanpur)
> Senior Scientist (E1) and Assistant Professor
> CSIR-National Metallurgical Laboratory Jamshedpur-831007
> http://www.nmlindia.org/
> https://scholar.google.co.in/citations?user=OchYqugAAAAJ&hl=en&oi=sra
GIUSEPPE MATTIOLI
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA
Via Salaria Km 29,300 - C.P. 10
I-00015 - Monterotondo Scalo (RM)
Mob (*preferred*) +39 373 7305625
Tel + 39 06 90672342 - Fax +39 06 90672316
E-mail: <giuseppe.mattioli at ism.cnr.it>
More information about the users
mailing list