<html><head><link type="text/css" rel="stylesheet" href="/gw_resource/css/namo_basic.css"><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body class=""><style type="text/css" id="NamoSE__GeneralStyle"> body{font-family :돋움; color : #000000; font-size : 10pt; margin : 7px 0 0 7px;} p,li{line-height:1.2; word-wrap: break-word; margin-top:0; margin-bottom:0;} body{overflow:auto;}.NamoSE_layoutlock_show { word-break: break-all;}</style><p>Dear QE users,</p><p><br></p><p>I would like to get your opinion about few points which are not clear to me:</p><p><br></p><p>- If I would like to calculate the Magnetic Anisotropy Energy (MAE) using the total energy difference (Non-collinear + SOC), should I include the Hubbard U correction during the calculation or it is not needed since it is an energy difference calculation ?</p><p><br></p><p>- If I suspect that my system is frustrated (it is Mn doped DMS) - because Mn itself has this non-collinear magnetic property - and I run a calculation of non-collinear scf calculation in-plane (IP) and out-of-plane (OOP) without SOC. Then I do the same with including SOC this time. What information we can get from the difference between {E(IP)-E(OOP)}without SOC and {E(IP)-E(OOP)}with SOC ?</p><p>Is it for example that there are another causes of the non-collinearity besides the SOC, like a geometrical frustration ?</p><p><br></p><p>- Last question is regarding the relaxation of a large supercell: I used 0.05 gaussian smearing. I didn't check the convergence w.r.t. to the smearing term, because of the expensive calculation. Do you have any idea if this smearing value is okay or may be it is very large ? (just for relaxing the system)<br></p><p>some details from the input file:</p><p><br></p><p>&control<br> calculation ='relax'<br> restart_mode='restart'<br> prefix = 'mn451_relax'<br> outdir = '/home/houcine/outdir'<br> pseudo_dir = '/home/houcine/pseudo'<br> nstep=200<br> etot_conv_thr = 1e-6<br> forc_conv_thr = 1e-4<br> verbosity = 'high'<br></p><p>&system<br> ibrav = 0,<br> nat = 60,<br> ntyp = 3,<br> ecutwfc=50<br> ecutrho=400<br> vdw_corr = 'dft-d3'<br> occupations='smearing',smearing='gaussian',degauss=0.05<br> nspin = 2<br> starting_magnetization(3)=0.5<br>/<br>&electrons<br> electron_maxstep=120<br> conv_thr=1e-6<br> mixing_beta = 0.2</p><p>ATOMIC_SPECIES<br>Sn 118.71000 Sn.pbe-dn-kjpaw_psl.1.0.0.UPF<br>S 32.06600 S.pbe-n-kjpaw_psl.1.0.0.UPF<br>Mn 54.93804 Mn.pbe-spn-kjpaw_psl.0.3.1.UPF<br></p><p>K_POINTS (automatic)<br>5 4 12 0 0 0<br><br>it is Mn-doped SnS2, semiconductor, in this calculation among the 20 Sn atoms only one was replaced by Mn</p><p><br></p><p>Thank you</p><p>Best Regards<br></p><p><br></p><div name="gwSign"><link type="text/css" rel="stylesheet" href="/gw_resource/css/namo_basic.css"><p><span style="color: rgb(0, 37, 222); font-family: Calibri, sans-serif; font-size: 10pt; background-color: rgb(255, 255, 255);">Houcine BOUZID </span><span style="color: rgb(0, 37, 222); font-family: Calibri, sans-serif; font-size: 10pt; background-color: rgb(255, 255, 255);"><br>Sungkyunkwan University, <br>Korea </span></p><span style="color: rgb(0, 173, 254); font-family: Calibri, sans-serif; font-size: 9pt; background-color: rgb(255, 255, 255);"></span><br></div><img src="http://ptmsg.skku.edu/emate_app/ematemdn.nsf/mdnform?OpenForm&sender=bouzid@skku.edu&receiver=users@lists.quantum-espresso.org&key=5CF4BE48AC12D74B492584E1004D82B6" width=0 height=0 border=0></body></html>