<div dir="ltr"><div dir="ltr"><div>Dear Experts,</div><div><br></div><div>I am trying 
to include 
spin-orbit coupling in energy calculations of a 3x3x3 supercell of a 
material. I looked it up on the internet and wrote the input code as 
given below (with atomic positions cropped but the complete input file 
is attached herewith) for vc-relax calculations. However, the generated output file does not contain any new "Atomic positions" or "Cell
 parameters". It looks more like an scf output file. Hence can you 
kindly let me know where I went wrong?</div><div><br></div><div>Input code:</div><div><span style="font-family:monospace">&CONTROL<br>  calculation = 'vc-relax'<br>  etot_conv_thr =   1.3500000000d-03<br>  forc_conv_thr =   1.0000000000d-04<br>  outdir = './outdir'<br>  prefix = 'cs'<br>  pseudo_dir = './'<br>  nstep = 500<br>  tprnfor = .true.<br>  tstress = .true.<br>  verbosity = 'high'<br>/<br>&SYSTEM<br>  degauss =   1.4699723600d-02<br>  ecutrho =   4.4000000000d+02<br>  ecutwfc =   5.5000000000d+01<br>  ibrav = 0<br>  nat = 135<br>  noncolin = .true.<br>  lspinorb = .true.<br>  !nosym = .false.<br>  !nspin = 2<br>  ntyp = 4<br>  occupations = 'smearing'<br>  smearing = 'cold'<br>  starting_magnetization(1) =   1.0000000000d-01<br>  starting_magnetization(2) =   1.0000000000d-01<br>  starting_magnetization(3) =   1.0000000000d-01<br>  starting_magnetization(4) =   1.0000000000d-01<br>/<br>&ELECTRONS<br>  conv_thr =   2.7000000000d-08<br>  electron_maxstep = 80<br>  mixing_beta =   4.0000000000d-01<br>/<br>&IONS<br>  ion_dynamics = 'bfgs'<br>/<br>&CELL<br>  cell_dynamics = 'bfgs'<br>/<br>ATOMIC_SPECIES<br>Cs     132.9054519 Cs.rel-pbesol-spn-kjpaw_psl.1.0.0.UPF <br>Cu     63.546 Cu.rel-pbesol-dn-kjpaw_psl.1.0.0.UPF<br>I      126.90447 I.rel-pbesol-n-kjpaw_psl.1.0.0.UPF<br>Pb     207.2 Pb.rel-pbesol-dn-kjpaw_psl.1.0.0.UPF<br>ATOMIC_POSITIONS crystal<br>Cs           0.1666669990       0.1666669990       0.1666669990 <br>Cs           0.5000000000       0.1666669990       0.1666669990 <br>Cs           0.8333330150       0.1666669990       0.1666669990 <br>Cs           0.1666669990       0.5000000000       0.1666669990 <br>Cs           0.5000000000       0.5000000000       0.1666669990 <br><br>K_POINTS automatic<br>2 2 2 0 0 0<br>CELL_PARAMETERS angstrom<br>     18.7120380402       0.0000000000       0.0000000000<br>      0.0000000000      18.7120380402       0.0000000000<br>      0.0000000000       0.0000000000      18.7120380402</span></div><div><br></div><div>Thanking you,</div><div>Tarik</div><div>Indian Institute of Technology</div><div>Indore</div><div>India</div></div></div>