<div dir="ltr"><div><font size="4" face="georgia, serif">Dear Experts,</font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">I want to simulate MoS2 and similar 2D systems with ferromagnetic calculation with and without SOC in  spin-polarised situation and a collinear ferromagnetic case.</font></div><div><font size="4" face="georgia, serif">Below is my research that I have done before posting this query here:</font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">example07 tells me to include </font></div><div><font size="4" face="georgia, serif">  Â  lspinorb=.true.,<br>  Â  noncolin=.true.,  to perform SOC with fully relativistic PPs.<br></font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">from example11,  I came to know that, if I insert </font></div><div><font size="4" face="georgia, serif">  Â  starting_magnetization(1)=0.5, in SYSTEM control then it will be ferromagnetic.<br></font></div><div><font size="4" face="georgia, serif">Additionally, we need to define "nspin=2" to do a spin-polarised calculation.</font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">Therefore, for a  ferromagnetic+SOC (spin-polarised)  calculation, the &SYSTEM control will be</font></div><div><font size="4" face="georgia, serif"> &SYSTEM<br></font></div><div><div><font size="4" face="georgia, serif">lspinorb=.true.,<br>noncolin=.true., </font></div><div><font size="4" face="georgia, serif">starting_magnetization(1)=0.5,</font></div><div><font size="4" face="georgia, serif">nspin=2</font></div></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">But then it becomes non-collinear (example11: QE-6.4.1), NO?? </font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">From here I have few basic question that needs answer:</font></div><div><font size="4" face="georgia, serif">1. Is, my below QE.in file is correct for  ferromagnetic+SOC+Spin-polarised case?</font></div><div><font size="4" face="georgia, serif">2. Can we do vc-relax with ferromagnetic+SOC+Spin-polarised?</font></div><div><font size="4" face="georgia, serif">3. What else I need to do if I want to perform a collinear calculation?</font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">I want to test spin-polarised ferromagnetic, spin-polarised ferromagnetic+SOC and collinear ferromagnetic with PBE and PBE+HSE on 2D systems one which is under mentioned:</font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">I have prepared below structure for ferromagneitc+SOC+spin-polarised calculation. </font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">&CONTROL<br>  Â  Â  Â  Â  Â  Â  Â  Â calculation = 'scf'<br>  Â  Â  Â  Â  Â  Â  Â  restart_mode = 'from_scratch'<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  outdir = './tmp'<br>  Â  Â  Â  Â  Â  Â  Â  Â  pseudo_dir = '.'<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  prefix = 'pwscf'<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â disk_io = 'none'<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â verbosity = 'default'<br>  Â  Â  Â  Â  Â  Â  Â etot_conv_thr = 0.0001<br>  Â  Â  Â  Â  Â  Â  Â forc_conv_thr = 0.001<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â nstep = 400<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â tstress = .true.<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â tprnfor = .true.<br> /<br> &SYSTEM<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ibrav = 0<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â nat = 4<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ntyp = 2<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ecutwfc = 60<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ecutrho = 480<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  vdw_corr = 'Grimme-D3  ! layered system<br> <br>  Â  Â  Â  Â  Â  Â  Â  Â occupations = 'smearing',  Â  ! Spin-polarised <br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  smearing = 'mv',<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â degauss = 0.01,<br><br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â nspin = 2,  Â  Â  Â  Â  Â  Â  Â  ! Spin-polarised<br>  Â starting_magnetization(1) = 0.6  ! Ferromagneitc</font></div><div><font size="4" face="georgia, serif">  Â starting_magnetization(2) = 0.6  Â  ! Ferromagneitc<br><br>  Â  Â  Â  Â  Â  Â  Â  Â  Â lspinorb = .true.,  Â  Â  Â  ! SOC<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â noncolin = .true.,  Â  Â  Â  !SOC<br><br><br> /<br> &ELECTRONS<br>  Â  Â  Â  Â  Â  electron_maxstep = 200<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  conv_thr = 1.0D-7<br>  Â  Â  Â  Â  Â  Â  diago_thr_init = 1e-4<br> /<br>&IONS<br>  Â  Â  Â  Â  Â  Â  Â  ion_dynamics = 'bfgs'<br> /<br><br>ATOMIC_SPECIES<br>  Â Mo Â  95.9400000000 Â Mo.rel_pz.upf<br>  Â  S Â  32.0650000000 Â S.upf<br><br>ATOMIC_POSITIONS crystal<br>Mo Â  Â  Â 0.99980000000000 Â  Â  Â 0.99967000000000 Â  Â  Â 0.00140000000000<br>Mo Â  Â  Â 0.41010000000000 Â  Â  Â 0.50033000000000 Â  Â  Â 0.99860000000000<br> S Â  Â  Â 0.76020000000000 Â  Â  Â 0.47161000000000 Â  Â  Â 0.12240000000000<br> S Â  Â  Â 0.27390000000000 Â  Â  Â 0.96106000000000 Â  Â  Â 0.16787000000000<br> S Â  Â  Â 0.64960000000000 Â  Â  Â 0.02839000000000 Â  Â  Â 0.87760000000000<br> S Â  Â  Â 0.13600000000000 Â  Â  Â 0.53894000000000 Â  Â  Â 0.83213000000000<br><br>K_POINTS automatic<br>3 Â 4 Â 1 Â  0 0 0<br><br>CELL_PARAMETERS {bohr}<br>  Â  10.77143895744817 Â  Â  Â 0.00000000000000 Â  Â  Â 0.00000000000000<br>  Â  Â 0.00000000000000 Â  Â  Â 6.06413116043003 Â  Â  Â 0.00000000000000<br>  Â  Â 1.28290955440224 Â  Â  Â 1.40801058817318 Â  Â  18.81430995418407<br></font></div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">I have gone thorough mailing list and example06,07,11,12 before writing this email.</font></div><div><div><div><div><font size="4" face="georgia, serif"><br></font></div><div><font size="4" face="georgia, serif">[1]. <a href="https://www.mail-archive.com/users@lists.quantum-espresso.org/msg17101.html" target="_blank">https://www.mail-archive.com/users@lists.quantum-espresso.org/msg17101.html</a></font></div><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Thank you very much for any help.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Best wishes</font></div><div><font size="4" face="georgia, serif"><br></font></div><div><div><font face="georgia, serif" color="#000000" size="4">Ms. Rekha</font></div><div><span style="color:rgb(0,0,0)"><font size="4" face="georgia, serif">Teaching Assistant,</font></span></div><div><font color="#000000" face="georgia, serif" size="4">Department of Physics</font></div><div><span style="color:rgb(0,0,0)"><font size="4" face="georgia, serif">RK college, Jaipur, </font></span><span style="color:rgb(0,0,0)"><font size="4" style="" face="georgia, serif">India</font><br></span></div></div></div></div></div></div></div></div></div>