Dear all,
I am currently working on an optimization for the ecutwfc parameter in Quantum
Espresso for HgTe, including spin-orbit coupling. However, I suspect there
might be an error in my input file.
Has anyone previously performed a Quantum Espresso calculation for HgTe
considering spin-orbit coupling? If so, could you kindly share your experience
or offer guidance on any potential mistakes in the input? Any help would be
greatly appreciated!
Thank you in advance for your time and assistance.
Best regards,
Sally
FYI:
PP:
Hg 200.59 Hg.rel-pz-dn-kjpaw_psl.0.2.2.UPF
Te 127.60 Te.rel-pz-dn-kjpaw_psl.0.2.2.UPF
------
for ECUT in 30 40 50 60 70; do
ECUTRHO=$((8 * ECUT))
here my input :>>>>
&CONTROL
calculation = 'vc-relax'
prefix = 'HgTe'
pseudo_dir = "$PSEUDO_DIR"
outdir = "$OUT_DIR"
nstep = 100 ! Limit relaxation steps
verbosity = 'high'
tprnfor = .true.
tstress = .true.
etot_conv_thr = 1.0d-6
forc_conv_thr = 1.0d-3
/
&SYSTEM
ibrav = 2
celldm(1) = 12.22 !bohr (Experimental lattice constant for HgTe = 6.46 Å)
nat = 2
ntyp = 2
ecutwfc = $ECUT
ecutrho = $ECUTRHO
lspinorb = .true. ! Include spin-orbit interactions
noncolin = .true. ! Non-collinear magnetism for SOC
/
&ELECTRONS
electron_maxstep = 500
conv_thr = 1.0d-9
mixing_beta = 0.7
diagonalization = 'david' ! cg: Conjugate Gradient method for better SCF
stability
/
&IONS
ion_dynamics = 'bfgs'
/
&CELL
cell_dynamics = 'bfgs'
cell_dofree = 'all'
/
$ATOMIC_SPECIES
$ATOMIC_POSITIONS
K_POINTS automatic
6 6 6 0 0 0
_______________________________________________________________________________
The Quantum ESPRESSO Foundation stands in solidarity with all civilians
worldwide who are victims of terrorism, military aggression, and indiscriminate
warfare.
--------------------------------------------------------------------------------
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
users mailing list [email protected]
https://lists.quantum-espresso.org/mailman/listinfo/users
