<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear QE users,<div><br></div><div>I'm trying to do a DFT+U calculation on the material monolayer 1T' WTe2, which is known to have strong SOC effects. It is mentioned in <a href="https://www.quantum-espresso.org/Doc/user_guide_PDF/Hubbard_input.pdf">Hubbard_input.pdf (quantum-espresso.org)</a> that "however, Liechtenstein’s DFT+U supports the noncollinear spin-polarized case, and so if you use this case then the code will automatically switch to the Liechtenstein’s DFT+U". The documentation however doesn't explicitly mention how to carry out a noncollinear DFT+U calculation. Currently I'm using the following input file:</div><div><br></div><div><div>&CONTROL</div><div>  calculation = 'scf'</div><div>  etot_conv_thr =   1.2000000000d-04</div><div>  forc_conv_thr =   1.0000000000d-04</div><div>  outdir = './'</div><div>  prefix = 'WTe2'</div><div>  pseudo_dir = '../pseudo/'</div><div>  tprnfor = .true.</div><div>  tstress = .true.</div><div>  verbosity = 'high'</div><div>/</div><div>&SYSTEM</div><div>  degauss =   7.3498618000d-03</div><div>  ecutrho =   640 </div><div>  ecutwfc =   80 </div><div>  ibrav = 0</div><div>  nat = 6 </div><div>  nosym = .false.</div><div>  ntyp = 2</div><div>  occupations = 'smearing'</div><div>  smearing = 'gaussian'</div><div>  noncolin = .true.</div><div>  lspinorb = .true.</div><div>/</div><div>&ELECTRONS</div><div>  conv_thr =   2.4000000000d-09</div><div>  electron_maxstep = 80</div><div>  mixing_beta =   4.0000000000d-01</div><div>/</div><div>ATOMIC_SPECIES</div><div>Te     127.6    Te_pbe_v1.uspp.F.UPF</div><div>W      183.84   W_pbe_v1.2.uspp.F.UPF</div><div>CELL_PARAMETERS (angstrom)</div><div>   6.336599029   0.000000000   0.000000000</div><div>   0.000000000   3.520370517   0.000000000</div><div>   0.000000000   0.000000000  20.160000000</div><div>ATOMIC_POSITIONS (crystal)</div><div>W             0.8214704044        0.7500000000        0.4951630000    1   1   0</div><div>W             0.1785295956        0.2500000000        0.5048380000    1   1   0</div><div>Te            0.0812386096        0.7500000000        0.6035080000    1   1   0</div><div>Te            0.5726099734        0.2500000000        0.5730930000    1   1   0</div><div>Te            0.4273900266        0.7500000000        0.4269070000    1   1   0</div><div>Te            0.9187613904        0.2500000000        0.3964920000    1   1   0</div><div>HUBBARD (ortho-atomic)</div><div>U W-5d 4</div><div>J W-5d 4</div><div>K_POINTS automatic</div><div>40 40 1 0 0 0</div></div><div><br></div><div>This input file however results in "Noncollinear case is not supported". I'm wondering how I should modify it to let it work.</div><div><br></div><div>Best regards,</div><div>Jinyuan</div></div></div></div></div></div>