<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=gb2312">
<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:等线;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"\@等线";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
text-align:justify;
text-justify:inter-ideograph;
font-size:10.5pt;
font-family:等线;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:等线;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
mso-ligatures:none;}
/* Page Definitions */
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;}
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="ZH-CN" link="#0563C1" vlink="#954F72" style="word-wrap:break-word;text-justify-trim:punctuation">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Dear QE users,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I am trying to run cpmd simulation with NPT ensemble to find the final lattice constant of ZrS<sup></sup><sub>3</sub> at specific temperature. I am wondering if Ekinc shows zero in evp file, and eigen values are all
zero at first 100 iterations being a problem or not.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The input is:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">###################################################################################<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&control<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">calculation= 'vc-cp'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">prefix= '300'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pseudo_dir= '../'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">nstep = 200000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">dt = 20.D0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">iprint = 100<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">saverho = .TRUE.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">restart_mode = 'from_scratch'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&SYSTEM<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> a = 1.01942e+01<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> b = 7.23461e+00<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> c = 2.00000e+01<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> cosab = -3.08225e-03<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> cosac = -1.26753e-01<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> cosbc = 1.95477e-04<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ibrav = 14<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> nat = 32<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ntyp = 2<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ecutwfc= 80.0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ecutrho= 320.0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> input_dft= 'pbe'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> vdw_corr= 'Grimme-D2'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&electrons<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> electron_maxstep = 500<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> conv_thr = 1e-06<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> electron_dynamics = 'cg'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> emass = 400.D0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> emass_cutoff = 3.D0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&ions<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ion_dynamics = 'verlet',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> ion_temperature = 'nose'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> fnosep = 6.D0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> tempw = 300<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&cell<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> cell_dynamics= 'pr'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> cell_dofree= '2Dxy'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> press = 0.001<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ATOMIC_SPECIES<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 91.22400 Zr_ONCV_PBE-1.2.upf<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 32.06600 S_ONCV_PBE-1.2.upf<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ATOMIC_POSITIONS {angstrom}<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 0.642320 0.919656 5.802935<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 0.631170 4.536946 5.802935<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 5.739425 0.919656 5.802935<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 5.728275 4.536946 5.802935<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 3.251654 2.707991 3.145483<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 3.240504 6.325280 3.145483<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 8.348759 2.707991 3.145483<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Zr 8.337609 6.325280 3.145483<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 3.167220 0.988938 4.906504<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 3.156070 4.606228 4.906504<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 8.264325 0.988938 4.906504<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 8.253175 4.606228 4.906504<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 0.712834 2.909611 3.950696<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 0.701684 6.526901 3.950696<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 5.809939 2.909611 3.950696<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 5.798789 6.526901 3.950696<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 2.347000 0.882416 1.538609<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 2.335851 4.499706 1.538609<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 7.444105 0.882416 1.538609<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 7.432956 4.499706 1.538609<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 1.726982 2.611313 7.291868<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 1.715832 6.228603 7.291868<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 6.824087 2.611313 7.291868<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 6.812937 6.228603 7.291868<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 4.389913 0.811184 1.393192<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 4.378764 4.428474 1.393192<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 9.487018 0.811184 1.393192<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 9.475869 4.428474 1.393192<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S -0.348116 2.608655 7.527313<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S -0.359266 6.225945 7.527313<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 4.748989 2.608655 7.527313<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">S 4.737839 6.225945 7.527313<o:p></o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-bottom:double windowtext 6.75pt;padding:0cm 0cm 1.0pt 0cm">
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The evp file:<o:p></o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-bottom:double windowtext 6.75pt;padding:0cm 0cm 1.0pt 0cm">
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span lang="EN-US"># nfi time(ps) ekinc T_cell(K) Tion(K) etot enthal econs econt Volume Pressure(GPa) EXX EVDW<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> 100 4.837769E-02 0.000000E+00 2.951422E+03 3.436941E+02 -626.40209263 -18.56732187 -626.34984874 -626.37459988 1.033521E+04 -0.02214<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US">200 9.675537E-02 0.000000E+00 6.723219E+02 3.693105E+02 -626.37151330 -35.54672732 -626.31537554 -626.36051078 1.094949E+04 -0.70109<o:p></o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-bottom:double windowtext 6.75pt;padding:0cm 0cm 1.0pt 0cm">
<p class="MsoNormal" style="text-indent:21.0pt;border:none;padding:0cm"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US">The eig file:<o:p></o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-bottom:double windowtext 6.75pt;padding:0cm 0cm 1.0pt 0cm">
<p class="MsoNormal" style="text-indent:21.0pt;border:none;padding:0cm"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> STEP: 100 0.05<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> Eigenvalues (eV), kp = 1 , spin = 1<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> STEP: 200 0.10<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> Eigenvalues (eV), kp = 1 , spin = 1<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-bottom:double windowtext 6.75pt;padding:0cm 0cm 1.0pt 0cm">
<p class="MsoNormal" style="text-indent:21.0pt;border:none;padding:0cm"><span lang="EN-US"> 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00<o:p></o:p></span></p>
</div>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US">The output file is attached.<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:21.0pt"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Best regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Jiaming Li<o:p></o:p></span></p>
</div>
</body>
</html>