<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-2">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-ligatures:standardcontextual;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hello,<o:p></o:p></p>
<p class="MsoNormal">I am trying to construct the proper input for vc-relax calculations for ZrSnS3, using following crystal data and qe input file, see below.<o:p></o:p></p>
<p class="MsoNormal">I’m getting "invalid np," error which typically occurs when the specified parameters for FFT are incorrect, if I read correctly.<o:p></o:p></p>
<p class="MsoNormal">I would appreciate any insights in correcting my input file.<o:p></o:p></p>
<p class="MsoNormal">Regards<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">CIF:<o:p></o:p></p>
<p class="MsoNormal">_chemical_name_systematic 'Zirconium Tin Sulfide'<o:p></o:p></p>
<p class="MsoNormal">_chemical_formula_moiety 'Zr Sn S3'<o:p></o:p></p>
<p class="MsoNormal">_chemical_formula_sum 'S3 Sn Zr'<o:p></o:p></p>
<p class="MsoNormal">_chemical_formula_weight 306.13<o:p></o:p></p>
<p class="MsoNormal">_cell_length_a 9.188(1)<o:p></o:p></p>
<p class="MsoNormal">_cell_length_b 3.717(1)<o:p></o:p></p>
<p class="MsoNormal">_cell_length_c 13.839(1)<o:p></o:p></p>
<p class="MsoNormal">_cell_angle_alpha 90<o:p></o:p></p>
<p class="MsoNormal">_cell_angle_beta 90<o:p></o:p></p>
<p class="MsoNormal">_cell_angle_gamma 90<o:p></o:p></p>
<p class="MsoNormal">_cell_formula_units_Z 4<o:p></o:p></p>
<p class="MsoNormal">_cell_volume 472.63<o:p></o:p></p>
<p class="MsoNormal">_symmetry_cell_setting orthorhombic<o:p></o:p></p>
<p class="MsoNormal">_symmetry_space_group_name_H-M Pnma<o:p></o:p></p>
<p class="MsoNormal">_symmetry_Int_Tables_number 62<o:p></o:p></p>
<p class="MsoNormal">loop_<o:p></o:p></p>
<p class="MsoNormal">_symmetry_equiv_pos_site_id<o:p></o:p></p>
<p class="MsoNormal">_symmetry_equiv_pos_as_xyz<o:p></o:p></p>
<p class="MsoNormal">1 x,y,z<o:p></o:p></p>
<p class="MsoNormal">2 -x,-y,-z<o:p></o:p></p>
<p class="MsoNormal">3 x+1/2,-y+1/2,-z+1/2<o:p></o:p></p>
<p class="MsoNormal">4 -x+1/2,y+1/2,z+1/2<o:p></o:p></p>
<p class="MsoNormal">5 -x,y+1/2,-z<o:p></o:p></p>
<p class="MsoNormal">6 x,-y+1/2,z<o:p></o:p></p>
<p class="MsoNormal">7 -x+1/2,-y,z+1/2<o:p></o:p></p>
<p class="MsoNormal">8 x+1/2,y,-z+1/2<o:p></o:p></p>
<p class="MsoNormal">loop_<o:p></o:p></p>
<p class="MsoNormal">_atom_site_label<o:p></o:p></p>
<p class="MsoNormal">_atom_site_type_symbol<o:p></o:p></p>
<p class="MsoNormal">_atom_site_fract_x<o:p></o:p></p>
<p class="MsoNormal">_atom_site_fract_y<o:p></o:p></p>
<p class="MsoNormal">_atom_site_fract_z<o:p></o:p></p>
<p class="MsoNormal">_atom_site_U_iso_or_equiv<o:p></o:p></p>
<p class="MsoNormal">_atom_site_symmetry_multiplicity<o:p></o:p></p>
<p class="MsoNormal">_atom_site_Wyckoff_symbol<o:p></o:p></p>
<p class="MsoNormal">_atom_site_occupancy<o:p></o:p></p>
<p class="MsoNormal">Sn1 Sn 0.46328 -0.25 0.1715 0.0177 4 c 1.0<o:p></o:p></p>
<p class="MsoNormal">Zr2 Zr 0.16377 0.25 0.04954 0.0106 4 c 1.0<o:p></o:p></p>
<p class="MsoNormal">S3 S 0.2742 0.25 0.21656 0.0143 4 c 1.0<o:p></o:p></p>
<p class="MsoNormal">S4 S 0.33181 -0.25 -0.0082 0.0111 4 c 1.0<o:p></o:p></p>
<p class="MsoNormal">S5 S 0.01227 0.25 -0.10761 0.0095 4 c 1.0<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Input file for qe:<o:p></o:p></p>
<p class="MsoNormal">&control<o:p></o:p></p>
<p class="MsoNormal"> calculation = 'vc-relax',<o:p></o:p></p>
<p class="MsoNormal"> restart_mode = 'from_scratch',<o:p></o:p></p>
<p class="MsoNormal"> prefix = '$NAME',<o:p></o:p></p>
<p class="MsoNormal"> verbosity = 'high',<o:p></o:p></p>
<p class="MsoNormal"> wf_collect = .true.,<o:p></o:p></p>
<p class="MsoNormal"> pseudo_dir = '$PSEUDO_DIR/',<o:p></o:p></p>
<p class="MsoNormal"> outdir = '$TMP_DIR/'<o:p></o:p></p>
<p class="MsoNormal">/<o:p></o:p></p>
<p class="MsoNormal">&system<o:p></o:p></p>
<p class="MsoNormal"> ecutwfc = 45.0,<o:p></o:p></p>
<p class="MsoNormal"> ibrav = 8,<o:p></o:p></p>
<p class="MsoNormal"> celldm(1) = 9.188,<o:p></o:p></p>
<p class="MsoNormal"> celldm(2) = 0.0404549,<o:p></o:p></p>
<p class="MsoNormal"> celldm(3) = 1.507, <o:p></o:p></p>
<p class="MsoNormal"> nat = 5, <o:p></o:p></p>
<p class="MsoNormal"> ntyp = 3, <o:p></o:p></p>
<p class="MsoNormal"> ecutrho = 360,<o:p></o:p></p>
<p class="MsoNormal"> nbnd=68, <o:p></o:p></p>
<p class="MsoNormal"> occupations = 'smearing', <o:p></o:p></p>
<p class="MsoNormal"> smearing = 'mv', <o:p></o:p></p>
<p class="MsoNormal"> degauss = 0.02, <o:p></o:p></p>
<p class="MsoNormal"> !space_group = 62, <o:p></o:p></p>
<p class="MsoNormal"> nosym = .false., <o:p></o:p></p>
<p class="MsoNormal"> noinv = .false.,<o:p></o:p></p>
<p class="MsoNormal">/<o:p></o:p></p>
<p class="MsoNormal">&electrons<o:p></o:p></p>
<p class="MsoNormal"> diagonalization = 'david',<o:p></o:p></p>
<p class="MsoNormal"> mixing_mode = 'plain',<o:p></o:p></p>
<p class="MsoNormal"> mixing_beta = 0.6,<o:p></o:p></p>
<p class="MsoNormal"> conv_thr = 1.0d-8<o:p></o:p></p>
<p class="MsoNormal">/<o:p></o:p></p>
<p class="MsoNormal">&ions<o:p></o:p></p>
<p class="MsoNormal"> ion_dynamics = 'bfgs'<o:p></o:p></p>
<p class="MsoNormal">/<o:p></o:p></p>
<p class="MsoNormal">&cell<o:p></o:p></p>
<p class="MsoNormal"> cell_dynamics = 'bfgs',<o:p></o:p></p>
<p class="MsoNormal"> press_conv_thr = 0.5D0<o:p></o:p></p>
<p class="MsoNormal">/<o:p></o:p></p>
<p class="MsoNormal">ATOMIC_SPECIES<o:p></o:p></p>
<p class="MsoNormal">Sn 118.710 Sn_pbe_v1.uspp.F.UPF<o:p></o:p></p>
<p class="MsoNormal">Zr 91.224 Zr_pbe_v1.uspp.F.UPF<o:p></o:p></p>
<p class="MsoNormal">S 32.065 s_pbe_v1.4.uspp.F.UPF<o:p></o:p></p>
<p class="MsoNormal">ATOMIC_POSITIONS {crystal}<o:p></o:p></p>
<p class="MsoNormal">Sn 0.46328 -0.25 0.1715<o:p></o:p></p>
<p class="MsoNormal">Zr 0.16377 0.25 0.04954<o:p></o:p></p>
<p class="MsoNormal">S 0.2742 0.25 0.21656<o:p></o:p></p>
<p class="MsoNormal">S 0.33181 -0.25 -0.0082<o:p></o:p></p>
<p class="MsoNormal">S 0.01227 0.25 -0.10761<o:p></o:p></p>
<p class="MsoNormal">K_POINTS {automatic}<o:p></o:p></p>
<p class="MsoNormal">6 6 4 0 0 0<o:p></o:p></p>
</div>
</body>
</html>