<p dir="ltr">Dear kanak <br>
Ecutwfc is too small. Find the optimal value. Set ecutrho = 5*ecutwfc.<br>
Set mixing mode local TF because of inhomogenity.<br>
Set conv thr of order 10^-9set dt= 20</p>
<p dir="ltr">Kind regards<br>
Ashkan Shekaari<br>
Tell: +98 933 459 7122; +98 921 346 7384</p>
<div class="gmail_quote">On Jul 25, 2015 5:59 PM, "Kanak Datta" <<a href="mailto:kanakeee08@gmail.com">kanakeee08@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear researchers<div><br></div><div>For BFGS relaxation of bilayer WS2..I have been using the following code...<br><div><font size="1">&CONTROL</font></div><div><font size="1">      calculation = 'relax',</font></div><div><font size="1">     restart_mode = 'from_scratch',</font></div><div><font size="1">           outdir = 'D:\QuantumEspresso\Quantum ESPRESSO 64-bit 5.1.2-mpich2\WS2Bilayer\outdir',</font></div><div><font size="1">           prefix = 'WS2Bilayer',</font></div><div><font size="1">/</font></div><div><font size="1">&SYSTEM</font></div><div><span style="font-size:x-small"> ibrav = 0,</span></div><div><span style="font-size:x-small">a = 1.86,</span><font size="1"><br></font></div><div><span style="font-size:x-small">nat = 6,</span><font size="1">                </font></div><div><div><font size="1">ntyp = 2,</font></div></div><div><span style="font-size:x-small">occupations = 'smearing',</span><font size="1"><br></font></div><div><span style="font-size:x-small">smearing = 'methfessel-paxton',</span><span style="font-size:x-small">            </span></div><div><span style="font-size:x-small">degauss = 0.001,</span><font size="1"><br></font></div><div><span style="font-size:x-small">nspin=2,</span><span style="font-size:x-small">          </span></div><div><span style="font-size:x-small">ecutwfc = 20,</span><font size="1">      </font></div><div><span style="font-size:x-small">ecutrho = 160,</span><font size="1"><br></font></div><div><span style="font-size:x-small">nbnd = 30,</span><font size="1">       </font></div><div><span style="font-size:x-small">starting_magnetization=0.6, </span><font size="1"><br></font></div><div><font size="1"> </font><span style="font-size:x-small">/</span></div><div><font size="1"><br></font></div><div><font size="1">&ELECTRONS</font></div><div><span style="font-size:x-small">conv_thr = 1.0d-5,</span><font size="1"><br></font></div><div><span style="font-size:x-small">mixing_mode = 'plain',</span><font size="1">         </font></div><div><span style="font-size:x-small">mixing_beta = 0.7,</span><font size="1"><br></font></div><div><span style="font-size:x-small">diagonalization = 'david',</span><font size="1">      </font></div><div><span style="font-size:x-small">diago_full_acc = .true.</span><font size="1"><br></font></div><div><font size="1"> </font><span style="font-size:x-small">/</span><span style="font-size:x-small">     </span></div><div><div><font size="1">&IONS </font></div><div><font size="1">ion_dynamics      = 'bfgs', </font></div><div><font size="1">pot_extrapolation = 'second_order', </font></div><div><font size="1">wfc_extrapolation = 'second_order', </font></div><div><font size="1">upscale           = 100,  </font></div><div><font size="1">/</font></div></div><div><font size="1">CELL_PARAMETERS </font></div><div><font size="1">1.5 0.8660254038 0.00000</font></div><div><font size="1">1.5 -0.8660254038 0.00000</font></div><div><font size="1">0.0000 0.00000 15.0  </font></div><div><span style="font-size:x-small">ATOMIC_SPECIES</span><br></div><div><font size="1">W 183.84 W.pbe-hgh.UPF</font></div><div><font size="1">S 32.066 S.pbe-hgh.UPF</font></div><div><span style="font-size:x-small">ATOMIC_POSITIONS (angstrom)</span><br></div><div><font size="1">S       0.000000000  0.000000000   5.15</font></div><div><font size="1">S       0.000000000  0.000000000   2</font></div><div><font size="1">W       1.818653348  0.000000000   3.575</font></div><div><font size="1">S       0.000000000  0.000000000   8.6</font></div><div><font size="1">S       0.000000000  0.000000000   10.175</font></div><div><font size="1">W       1.818653348  0.000000000   11.75</font></div><div><span style="font-size:x-small">K_POINTS {automatic}</span><br></div><div><span style="font-size:x-small">24 24 1   0 0 0</span><font size="1">  </font></div><div><font size="1">  </font></div><div><br></div>1. I just wanted to know whether the cutoff values are ok or not.......As I am from EEE background..I have very little knowledge on these things.....<br>2. Besides should the van der waals corrections be applied in relaxation part from the beginning?</div><div><br></div><div>Thanks in advance.</div><div><br></div><div>Sincerely yours </div><div>kanak</div><div>BUET</div></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div>