<div dir="ltr"><pre>Dear <span id="m_-8361387898297949902m_5422436294486064676:53f.4">QuantumEspresso</span> user group,<br><br></pre><pre>i have a problem with my calculations. I investigate tungsten carbide (WC) layer with 14 A vacuum and adsorb atomic oxygen on the surface.<br></pre><pre>I use PWSCF v.6.3MaX <br></pre><pre>If try a relax Job with a starting magnetization for the oxygen, then the following happens while the SCF-Cycles:<br></pre><pre>The first optimization step converge in a range of 17 iterations. <br></pre><pre>Step 1:</pre><pre><span id="m_-8361387898297949902m_5422436294486064676:53f.6"></span>     negative rho (up, down):  1.163E-01 7.746E-02<br>     total energy              =  -18470.12492847 Ry<br>     Harris-Foulkes estimate   =  -18470.12491708 Ry<br>     estimated scf accuracy    <       0.00000263 Ry<br>     magnetization             =    -1.69 Bohr mag/cell<br>     absolute magnetization    =     2.25 Bohr mag/cell<br>--<br>     negative rho (up, down):  1.168E-01 7.723E-02<br>     total energy              =  -18470.12493846 Ry<br>     Harris-Foulkes estimate   =  -18470.12493246 Ry<br>     estimated scf accuracy    <       0.00000146 Ry<br>     total magnetization       =    -1.72 Bohr mag/cell<br>     absolute magnetization    =     2.30 Bohr mag/cell<br>--<br>     negative rho (up, down):  1.176E-01 7.701E-02<br>     total cpu time spent up to now is     3103.3 secs<br>     total energy              =  -18470.12494014 Ry<br>     Harris-Foulkes estimate   =  -18470.12494049 Ry<br>     estimated scf accuracy    <       0.00000116 Ry<br>     total magnetization       =    -1.75 Bohr mag/cell<br>     absolute magnetization    =     2.34 Bohr mag/cell</pre><pre>All seems in a region to be okay for the first step, but then some parameter going crazy,<br>Step 2<br>     negative rho (up, down):  1.103E+02 1.177E+02<br>     total energy              =  -18691.45057043 Ry<br>     Harris-Foulkes estimate   =  -18587.46248750 Ry<br>     estimated scf accuracy    <   45389.31161404 Ry<br>     total magnetization       =    -1.39 Bohr mag/cell<br>     absolute magnetization    =    36.78 Bohr mag/cell<br>--<br>     negative rho (up, down):  1.211E+02 1.247E+02<br>     total energy              =  -26401.86319230 Ry<br>     Harris-Foulkes estimate   =  -20544.68384009 Ry <br>     estimated scf accuracy    <   91328.60181762 Ry<br>     total magnetization       =    -0.44 Bohr mag/cell<br>     absolute magnetization    =    23.60 Bohr mag/cell<br><br></pre><pre>I try a lot of standard things like:<br><br></pre><pre>testing: k-points 221 - 881 ,  ecutwfc ecutrho 60-100 600-1000, different smearing methods and degauss values 0.005 - 0.04, <br></pre><pre>change mixing_mode, mixing_beta between 0.7-0.1, different diagonalization themes, nbnd between 196-300, some pseudo potentials PAW/US,<br>verify my input via xcrysden, using nspin=1 instead of 2 </pre><pre> <br></pre><pre>but nothing works. Can someone explain to me how the absolute magnetization can rise up to 600 and higher? While the scf accuracy<br></pre><pre>goes to infinity?<br></pre><pre>And what can I do to get rid of this errors? I have no ideas left ... maybe I overlook something easy.<br><br></pre><div class="m_-8361387898297949902gmail-oSioSc"><div id="m_-8361387898297949902gmail-tw-target"><div id="m_-8361387898297949902gmail-kAz1tf" class="m_-8361387898297949902gmail-g9WsWb"><div class="m_-8361387898297949902gmail-tw-ta-container m_-8361387898297949902gmail-tw-nfl" id="m_-8361387898297949902gmail-tw-target-text-container"><pre class="m_-8361387898297949902gmail-tw-data-text m_-8361387898297949902gmail-tw-ta m_-8361387898297949902gmail-tw-text-small" id="m_-8361387898297949902gmail-tw-target-text" style="text-align:left" dir="ltr"><span lang="en">I will print out one of my millions of input data at the end.<br><br></span></pre><pre><span lang="en">Many thanks and <br>Kind Regards,</span><br><br></pre><pre>Wigbert Polet<br></pre><pre><br></pre><pre>MY FIRST INPUT:<br><br>// n_make_cells.sh<br>&control<br> calculation = 'relax',<br> pseudo_dir = 'NODES_QE_PP',<br> prefix = 'NODES_QE_PREFIX',<br> outdir = 'NODES_QE_OUT',<br> nstep = 200,<br> etot_conv_thr = 0.000001,<br> forc_conv_thr = 0.0003,<br> /<br><br>&system<br> ibrav = 0,<br> nat = 49, ntyp = 3,<br> ecutwfc = 60, ecutrho = 600, nbnd = 280,<br> nspin = 2,<br> starting_magnetization(1)=0.0,<br> starting_magnetization(2)=0.0,<br> starting_magnetization(3)=1.0,<br> occupations = 'smearing',<br> smearing = 'methfessel-paxton',<br> degauss = 0.005,<br> input_dft = 'vdw-df2-b86r'<br> /<br><br>&ELECTRONS<br> conv_thr = 1.D-6,<br> mixing_mode='local-TF',<br> mixing_beta = 0.5,<br> diagonalization = 'david',<br> electron_maxstep = 250<br> /<br><br>&IONS<br> ion_dynamics='bfgs'<br> /<br><br>ATOMIC_SPECIES<br>C   12.00000000    C.pbe-n-rrkjus_psl.1.0.0.UPF<br>W  183.9509312     W.pbe-spn-rrkjus_psl.1.0.1.UPF (This doesn't work see <a href="https://dalcorso.github.io/pslibrary/PP_list.html">https://dalcorso.github.io/pslibrary/PP_list.html</a>) We used later other<br>O   15.99491461956 O.pbe-n-rrkjus_psl.1.0.0.UPF<br><br>CELL_PARAMETERS (angstrom)<br>8.502074  0.000000   0.000000<br>0.000000 10.084398   0.000005<br>0.000000 -0.000010 -19.823287<br><br>ATOMIC_POSITIONS (angstrom)<br> C 7.085057 8.407257 0.000011 0 0 0<br> C 4.251032 8.407257 0.000011 0 0 0<br> C 1.417007 8.407257 0.000011 0 0 0<br> C 7.085057 3.365058 0.000008 0 0 0<br> C 4.251032 3.365058 0.000008 0 0 0<br> C 1.417007 3.365058 0.000008 0 0 0<br> W 5.668050 5.038705 0.000002 0 0 0<br> W 2.834025 5.038705 0.000002 0 0 0<br> W 0.000000 5.038705 0.000002 0 0 0<br> W 5.668050 -0.003494 -0.000001 0 0 0<br> W 2.834025 -0.003494 -0.000001 0 0 0<br> W 0.000000 -0.003494 -0.000001 0 0 0<br> W 5.668050 7.559787 -1.455798 1 1 1<br> W 2.834025 7.559787 -1.455798 1 1 1<br> W 0.000000 7.559787 -1.455798 1 1 1<br> W 5.668050 2.517588 -1.455800 1 1 1<br> W 2.834025 2.517588 -1.455800 1 1 1<br> W 0.000000 2.517588 -1.455800 1 1 1<br> C 7.085057 5.886154 -1.455824 1 1 1<br> C 4.251032 5.886154 -1.455824 1 1 1<br> C 1.417007 5.886154 -1.455824 1 1 1<br> C 7.085057 0.843955 -1.455827 1 1 1<br> C 4.251032 0.843955 -1.455827 1 1 1<br> C 1.417007 0.843955 -1.455827 1 1 1<br> C 7.085057 8.407256 -2.911633 1 1 1<br> C 4.251032 8.407256 -2.911633 1 1 1<br> C 1.417007 8.407256 -2.911633 1 1 1<br> C 7.085057 3.365057 -2.911635 1 1 1<br> C 4.251032 3.365057 -2.911635 1 1 1<br> C 1.417007 3.365057 -2.911635 1 1 1<br> W 5.668050 5.038704 -2.911642 1 1 1<br> W 2.834025 5.038704 -2.911642 1 1 1<br> W 0.000000 5.038704 -2.911642 1 1 1<br> W 5.668050 -0.003495 -2.911644 1 1 1<br> W 2.834025 -0.003495 -2.911644 1 1 1<br> W 0.000000 -0.003495 -2.911644 1 1 1<br> W 5.668050 7.559786 -4.367441 1 1 1<br> W 2.834025 7.559786 -4.367441 1 1 1<br> W 0.000000 7.559786 -4.367441 1 1 1<br> W 5.668050 2.517587 -4.367444 1 1 1<br> W 2.834025 2.517587 -4.367444 1 1 1<br> W 0.000000 2.517587 -4.367444 1 1 1<br> C 7.085057 5.886152 -4.367467 1 1 1<br> C 4.251032 5.886152 -4.367467 1 1 1<br> C 1.417007 5.886152 -4.367467 1 1 1<br> C 7.085057 0.843953 -4.367470 1 1 1<br> C 4.251032 0.843953 -4.367470 1 1 1<br> C 1.417007 0.843953 -4.367470 1 1 1<br> O 2.834025 2.517587 -6.367444 1 1 1<br><br>K_POINTS {automatic}<br>2 2 1 0 0 0<br><br><br><br></pre></div></div></div></div><pre> <br></pre><pre><br><br><br><br></pre></div>