<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;">Dear QE Fellows,</span><br style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><br style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;">I am trying to apply LDA+U+SOC Scheme on Mn-doped silicon system (Magnetism is also incorporated). </span><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13px;">I am using version espresso-5.0.2.</span><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"> After running till first scf cycle (couple of iterations), when it comes to write the forces acting and coordinates of relaxed system, the calculation is crashed with the following error:</span><br style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><br style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><font color="#ff0000" style="font-family: arial, sans-serif; font-size: 13.333333969116211px;"> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%<wbr>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<wbr>%%%%%%%%%%%%%%%%%%%<br>     Error in routine force_hub (1):<br>      forces in full LDA+U scheme are not yet implemented<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%<wbr>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<wbr>%%%%%%%%%%%%%%%%%%%<br><br>     stopping ...<br><br>.......</font><br style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><br></span><div><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><br></span></div><div><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;">Can somebody check the input and suggest some solution? </span><div style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><br><br>Thank you in advance. Below is the input:<br><br><br>------------------------------<wbr>------------------------------<wbr>-----------------<br>&control<br>    calculation='relax',<br>    restart_mode='from_scratch',<br>    prefix='sil',<br>   ! verbosity='high',<br>   ! tprnfor=.true.<br>    pseudo_dir = '/home/sattar',<br>    etot_conv_thr=1.0d-4,           ! Total Energy conf. for ionic min.<br>    forc_conv_thr=1.0d-3,           ! Force conv for ion min.<br>   ! tstress=.true.<br>    nstep=150,<br>    outdir='tmp/',<br> /<br>&system<br>    ibrav = 4,<br>    a=15.44, b=15.44, c=25,<br>    cosab=-0.5, cosbc=0.0, cosac=0.0<br>    nat=34, ntyp= 3,<br>    ecutwfc = 40,ecutrho=350,<br>    occupations='smearing', smearing='fd', degauss=0.01<br>    nspin =4,<br>    starting_magnetization(1)= 0.0<br>    starting_magnetization(2)= 0.5<br>    lda_plus_u = .true.<br>    lda_plus_u_kind = 1<br>    lspinorb=.true.,<br>    noncolin=.true.,<br>    Hubbard_U(2) = 4.0,<br>    Hubbard_J=0.0<br>    nosym = .FALSE.<br> /<br> &electrons<br>    mixing_beta = 0.3,<br>    conv_thr =  1.0d-7,<br>    electron_maxstep=300,<br> /<br>&CELL<br>    cell_dynamics='bfgs',<br> /<br>&IONS<br>    ion_dynamics='bfgs',<br>/<br>ATOMIC_SPECIES<br>Si  28.086  Si.pbe-n-rrkjus_psl.0.1.UPF<br>Mn 54.938  Mn.rel-pbe-spn-rrkjus_psl.0.<wbr>2.3.UPF<br>ATOMIC_POSITIONS (angstrom)<br>Si ....<br>Mn ....<br>;;;<br>;;;<br>;;;<br>;;;<br>;;;<br><br>K_POINTS AUTOMATIC<br>4 4 1 0 0 0<br>------------------------------<wbr>------------------------------<wbr>--------------------</div><div class="yj6qo ajU" style="cursor: pointer; outline: none; padding: 10px 0px; width: 22px; color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><div id=":2ln" class="ajR ajS" role="button" tabindex="0" data-tooltip="Hide expanded content" aria-label="Hide expanded content" style="background-color: rgb(216, 216, 216); border: 1px solid rgb(205, 205, 205); clear: both; line-height: 6px; outline: none; position: relative; width: 20px;"><img class="ajT" src="https://mail.google.com/mail/u/0/images/cleardot.gif" style="background-image: url(https://ssl.gstatic.com/ui/v1/icons/mail/ellipsis.png); height: 8px; opacity: 0.9; width: 20px; background-position: initial initial; background-repeat: no-repeat no-repeat;"></div></div><span class="HOEnZb adL" style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13.333333969116211px;"><font color="#888888"><div><br></div><div><br><br><br>-- <br>Kind Regards,<br>*Shahid Sattar*<br>Ph. D. Student<br>Material Science and Engineering (MSE)<br></div><div>KAUST</div></font></span></div>                                        </div></body>
</html>