Dear Dr. Anupriya Nyayban,

so maybe also 'cell_dofree' should be relevant in your case.


Best regards
Chiara Cignarella
PhD THEOS - EPFL
Switzerland

________________________________
Da: users <users-boun...@lists.quantum-espresso.org> per conto di Anupriya 
Nyayban <anupriya...@phy.nits.ac.in>
Inviato: venerdì 26 gennaio 2024 11:47:11
A: users@lists.quantum-espresso.org
Oggetto: Re: [QE-users] Issues with vc-relax for a slab

Dear Chiara Cignarella,

Yes, this has a vacuum of 20 Angstrom along the z direction.
 I will check by reducing the mixing_beta first and then check with 
mixing_mode/ mixing_ndim.

Thank you!

with regards
Anupriya Nyayban
Postdoc, IIT Madras

On Fri, Jan 26, 2024 at 3:58 PM Anupriya Nyayban 
<anupriya...@phy.nits.ac.in<mailto:anupriya...@phy.nits.ac.in>> wrote:
Dear Dr. Chirantan Pramanik,

Thank you for the reply!
No, I haven't done vc-relax for a surface earlier but with relax also I am 
facing a similar problem. However, I will try 'relax calculation' by increasing 
the electron_maxstep.

with regards
Anupriya Nyayban
Postdoc, IIT Madras


On Fri, Jan 26, 2024 at 12:44 PM Anupriya Nyayban 
<anupriya...@phy.nits.ac.in<mailto:anupriya...@phy.nits.ac.in>> wrote:
Dear users and experts,

I have considered a (311) plane of 2 layers of ZnFe2O4 structure. The attached 
input file is for the vc-relax calculation.

&CONTROL
  calculation='vc-relax'
  restart_mode='from_scratch',
  prefix='zn'
  pseudo_dir = '.'
  outdir='tmp'
  verbosity='high'
  etot_conv_thr = 1.0D-5
  forc_conv_thr = 1.0D-4
/
&SYSTEM
  ibrav = 0
  A =   21.62370
  nat = 112
  ntyp = 3
  ecutwfc = 30.0,
  ecutrho = 240.0,
  occupations = 'smearing',
  smearing = 'mv',
  degauss = 0.02,
  vdw_corr  = 'DFT-D3'
/
&ELECTRONS
   conv_thr =  1.d-9
   mixing_beta = 0.3
  electron_maxstep = 300
  startingpot = "atomic"
  startingwfc = "atomic+random"
/
&IONS
  ion_dynamics = "bfgs"
/
&CELL
  cell_dynamics = "bfgs"
  press_conv_thr = 5.00000e-01
/
CELL_PARAMETERS {alat}
  1.000000000000000   0.000000000000000   0.000000000000000
  0.299999977738132   0.331662474918549   0.000000000000000
  0.000000000000000   0.000000000000000   2.019068156482969
ATOMIC_SPECIES
  Zn   65.38000  Zn_ONCV_PBE-1.0.upf
  Fe   55.84500  Fe_ONCV_PBE-1.0.upf
   O   15.99900  O_ONCV_PBE-1.0.upf
ATOMIC_POSITIONS {crystal}
Zn   0.863636000000000   0.954545000000000   0.475797000000000
Zn   0.863636000000000   0.454545000000000   0.475797000000000
Zn   0.363636000000000   0.454545000000000   0.475797000000000
Zn   0.363636000000000   0.954545000000000   0.475797000000000
Zn   0.545454000000000   0.181819000000000   0.487603000000000
Zn   0.045455000000000   0.181818000000000   0.487603000000000
Zn   0.045455000000000   0.681818000000000   0.487603000000000
Zn   0.545455000000000   0.681818000000000   0.487602000000000
Zn   0.590909000000000   0.863636000000000   0.523020000000000
Zn   0.590909000000000   0.363636000000000   0.523020000000000
Zn   0.090909000000000   0.363636000000000   0.523020000000000
Zn   0.090909000000000   0.863636000000000   0.523020000000000
Zn   0.272727000000000   0.090909000000000   0.534825000000000
Zn   0.772727000000000   0.090909000000000   0.534826000000000
Zn   0.772727000000000   0.590909000000000   0.534826000000000
Zn   0.272727000000000   0.590909000000000   0.534825000000000
Fe   0.954546000000000   0.068182000000000   0.481700000000000
Fe   0.704545000000000   0.318182000000000   0.481700000000000
Fe   0.204546000000000   0.568182000000000   0.481700000000000
Fe   0.590909000000000   0.863636000000000   0.458088000000000
Fe   0.204546000000000   0.818181000000000   0.481700000000000
Fe   0.704545000000000   0.568182000000000   0.481700000000000
Fe   0.590909000000000   0.363636000000000   0.458088000000000
Fe   0.954545000000000   0.568182000000000   0.481700000000000
Fe   0.704546000000000   0.818182000000000   0.481700000000000
Fe   0.204545000000000   0.068182000000000   0.481700000000000
Fe   0.204545000000000   0.318182000000000   0.481700000000000
Fe   0.704546000000000   0.068181000000000   0.481700000000000
Fe   0.454545000000000   0.068182000000000   0.481700000000000
Fe   0.090909000000000   0.863637000000000   0.458088000000000
Fe   0.090909000000000   0.363636000000000   0.458088000000000
Fe   0.454545000000000   0.568182000000000   0.481700000000000
Fe   0.681818000000000   0.977272000000000   0.528922000000000
Fe   0.431818000000000   0.227273000000000   0.528922000000000
Fe   0.931818000000000   0.477273000000000   0.528923000000000
Fe   0.318182000000000   0.772727000000000   0.505311000000000
Fe   0.931818000000000   0.727272000000000   0.528923000000000
Fe   0.431818000000000   0.477273000000000   0.528922000000000
Fe   0.318182000000000   0.272727000000000   0.505311000000000
Fe   0.681818000000000   0.477273000000000   0.528923000000000
Fe   0.431818000000000   0.727273000000000   0.528923000000000
Fe   0.931818000000000   0.977273000000000   0.528923000000000
Fe   0.931818000000000   0.227273000000000   0.528923000000000
Fe   0.431818000000000   0.977272000000000   0.528923000000000
Fe   0.181818000000000   0.977273000000000   0.528923000000000
Fe   0.818182000000000   0.772728000000000   0.505311000000000
Fe   0.818182000000000   0.272727000000000   0.505311000000000
Fe   0.181818000000000   0.477273000000000   0.528923000000000
 O   0.021208000000000   0.351541000000000   0.468710000000000
 O   0.642439000000000   0.036341000000000   0.493900000000000
 O   0.269690000000000   0.845466000000000   0.471866000000000
 O   0.521208000000000   0.584834000000000   0.468711000000000
 O   0.521208000000000   0.084834000000000   0.468710000000000
 O   0.021208000000000   0.851541000000000   0.468710000000000
 O   0.269690000000000   0.345466000000000   0.471866000000000
 O   0.642439000000000   0.536342000000000   0.493900000000000
 O   0.142439000000000   0.536342000000000   0.493900000000000
 O   0.021208000000000   0.084834000000000   0.468710000000000
 O   0.769690000000000   0.345466000000000   0.471866000000000
 O   0.521208000000000   0.851541000000000   0.468711000000000
 O   0.887883000000000   0.284822000000000   0.494689000000000
 O   0.766652000000000   0.100022000000000   0.469499000000000
 O   0.387883000000000   0.551530000000000   0.494689000000000
 O   0.139401000000000   0.790898000000000   0.491533000000000
 O   0.639401000000000   0.790898000000000   0.491533000000000
 O   0.387883000000000   0.284822000000000   0.494689000000000
 O   0.887883000000000   0.551530000000000   0.494689000000000
 O   0.266652000000000   0.100022000000000   0.469499000000000
 O   0.266652000000000   0.600022000000000   0.469499000000000
 O   0.639401000000000   0.290898000000000   0.491533000000000
 O   0.887882000000000   0.051530000000000   0.494689000000000
 O   0.387883000000000   0.784822000000000   0.494689000000000
 O   0.887883000000000   0.784822000000000   0.494689000000000
 O   0.387883000000000   0.051530000000000   0.494689000000000
 O   0.766652000000000   0.600022000000000   0.469499000000000
 O   0.142439000000000   0.036341000000000   0.493900000000000
 O   0.021208000000000   0.584834000000000   0.468710000000000
 O   0.521209000000000   0.351541000000000   0.468710000000000
 O   0.769690000000000   0.845466000000000   0.471866000000000
 O   0.139401000000000   0.290898000000000   0.491533000000000
 O   0.748481000000000   0.260632000000000   0.515933000000000
 O   0.369712000000000   0.945432000000000   0.541123000000000
 O   0.996962000000000   0.754557000000000   0.519089000000000
 O   0.248481000000000   0.493925000000000   0.515933000000000
 O   0.248481000000000   0.993924000000000   0.515933000000000
 O   0.748481000000000   0.760632000000000   0.515933000000000
 O   0.996962000000000   0.254556000000000   0.519089000000000
 O   0.369712000000000   0.445432000000000   0.541123000000000
 O   0.869712000000000   0.445433000000000   0.541123000000000
 O   0.748481000000000   0.993924000000000   0.515933000000000
 O   0.496962000000000   0.254556000000000   0.519089000000000
 O   0.248481000000000   0.760632000000000   0.515933000000000
 O   0.615155000000000   0.193913000000000   0.541912000000000
 O   0.493924000000000   0.009113000000000   0.516722000000000
 O   0.115155000000000   0.460621000000000   0.541912000000000
 O   0.866674000000000   0.699989000000000   0.538756000000000
 O   0.366674000000000   0.699989000000000   0.538756000000000
 O   0.115155000000000   0.193913000000000   0.541912000000000
 O   0.615155000000000   0.460621000000000   0.541912000000000
 O   0.993924000000000   0.009113000000000   0.516722000000000
 O   0.993924000000000   0.509113000000000   0.516722000000000
 O   0.366674000000000   0.199989000000000   0.538756000000000
 O   0.615155000000000   0.960621000000000   0.541912000000000
 O   0.115155000000000   0.693913000000000   0.541912000000000
 O   0.615155000000000   0.693913000000000   0.541912000000000
 O   0.115155000000000   0.960621000000000   0.541912000000000
 O   0.493924000000000   0.509113000000000   0.516722000000000
 O   0.869712000000000   0.945432000000000   0.541123000000000
 O   0.748481000000000   0.493925000000000   0.515933000000000
 O   0.248481000000000   0.260632000000000   0.515933000000000
 O   0.496962000000000   0.754557000000000   0.519089000000000
 O   0.866674000000000   0.199989000000000   0.538756000000000
K_POINTS {automatic}
 2 2 1 1 1 1

The following is how the energy changes after each optimization step. At the 
last optimization step ( Energy error = 1.5E+00 Ry), scf didn't converge even 
after 300 iterations.
    Energy error            =      0.0E+00 Ry
     Energy error            =      1.9E+01 Ry
     Energy error            =      2.5E+01 Ry
     Energy error            =      2.7E+01 Ry
     Energy error            =      2.7E+01 Ry
     Energy error            =      2.7E+01 Ry
     Energy error            =      2.4E+01 Ry
     Energy error            =      2.0E+01 Ry
     Energy error            =      1.5E+01 Ry
     Energy error            =      1.1E+01 Ry
     Energy error            =      8.3E+00 Ry
     Energy error            =      6.5E+00 Ry
     Energy error            =      5.1E+00 Ry
     Energy error            =      4.3E+00 Ry
     Energy error            =      3.7E+00 Ry
     Energy error            =      3.2E+00 Ry
     Energy error            =      2.8E+00 Ry
     Energy error            =      2.5E+00 Ry
     Energy error            =      2.2E+00 Ry
     Energy error            =      2.0E+00 Ry
     Energy error            =      1.7E+00 Ry
     Energy error            =      1.5E+00 Ry
Please suggest how I can overcome the issue and what may be the reason! It will 
be really helpful!


Thank you!
Anupriya Nyayban
Postdoc, IIT Madras

On Wed, Jan 24, 2024 at 6:02 PM Anupriya Nyayban 
<anupriya...@phy.nits.ac.in<mailto:anupriya...@phy.nits.ac.in>> wrote:
Dear users and experts,

I have considered a (311) plane of 2 layers of ZnFe2O4 structure. The attached 
input file is for the vc-relax calculation. The following is how the energy 
changes after each optimization step. At the last optimization step ( Energy 
error = 1.5E+00 Ry), scf didn't converge even after 300 iterations.
    Energy error            =      0.0E+00 Ry
     Energy error            =      1.9E+01 Ry
     Energy error            =      2.5E+01 Ry
     Energy error            =      2.7E+01 Ry
     Energy error            =      2.7E+01 Ry
     Energy error            =      2.7E+01 Ry
     Energy error            =      2.4E+01 Ry
     Energy error            =      2.0E+01 Ry
     Energy error            =      1.5E+01 Ry
     Energy error            =      1.1E+01 Ry
     Energy error            =      8.3E+00 Ry
     Energy error            =      6.5E+00 Ry
     Energy error            =      5.1E+00 Ry
     Energy error            =      4.3E+00 Ry
     Energy error            =      3.7E+00 Ry
     Energy error            =      3.2E+00 Ry
     Energy error            =      2.8E+00 Ry
     Energy error            =      2.5E+00 Ry
     Energy error            =      2.2E+00 Ry
     Energy error            =      2.0E+00 Ry
     Energy error            =      1.7E+00 Ry
     Energy error            =      1.5E+00 Ry
Please suggest how I can overcome the issue and what may be the reason!

Thank you!

With regards
Anupriya Nyayban
IIT Madras, India
_______________________________________________
The Quantum ESPRESSO community stands by the Ukrainian
people and expresses its concerns about the devastating
effects that the Russian military offensive has on their
country and on the free and peaceful scientific, cultural,
and economic cooperation amongst peoples
_______________________________________________
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
users mailing list users@lists.quantum-espresso.org
https://lists.quantum-espresso.org/mailman/listinfo/users

Reply via email to