<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi Kejiang,</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
Structural relaxation in cp.x (step 2 as you mentioned) can involve the cell degrees of freedom (in addition to the ionic ones) and thus is in principle equivalent to doing "vc-relax" in pw.x. As such, it is unnecessary to do a "vc-relax" with pw.x.</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
In practice, proper parameter setting relies on some intuition of the targeted physical system. I think the difficulty in force convergence probably is related to "emass" and "dt" (e.g., when the emass is set too large, ionic motion during relaxation could
fictitiously heat up the electronic degrees of freedom and lead to fluctuations in forces). Nevertheless, detailed input/output information would be needed for more definitive comments.<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
Best,</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hsin-Yu<br>
</div>
<div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="Signature">
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<font color="#888888">--
<div dir="ltr">
<div dir="ltr">
<div><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Hsin-Yu Ko</span></span></div>
<div><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Postdoctoral Research Fellow<br>
</span></span></div>
<span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Department of Chemistry and Chemical Biology<br>
</span></span><span><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Cornell University</span></span></span></div>
</div>
</font><br>
</div>
</div>
</div>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of likejiang@ustb.edu.cn <likejiang@ustb.edu.cn><br>
<b>Sent:</b> Friday, April 8, 2022 2:44 AM<br>
<b>To:</b> users@lists.quantum-espresso.org <users@lists.quantum-espresso.org>; users-request@lists.quantum-espresso.org <users-request@lists.quantum-espresso.org><br>
<b>Subject:</b> [QE-users] Is it necessary to do "vc-relax" with pw.x before doing CPMD with cp.x ?</font>
<div> </div>
</div>
<style>
<!--
@font-face
{font-family:SimSun}
@font-face
{font-family:"Cambria Math"}
@font-face
{font-family:DengXian}
@font-face
{font-family:Calibri}
@font-face
{font-family:SimSun}
@font-face
{font-family:DengXian}
p.x_MsoNormal, li.x_MsoNormal, div.x_MsoNormal
{margin:0cm;
text-align:justify;
text-justify:inter-ideograph;
font-size:10.5pt;
font-family:DengXian}
a:link, span.x_MsoHyperlink
{color:#0563C1;
text-decoration:underline}
span.x_EmailStyle21
{font-family:DengXian;
color:windowtext}
.x_MsoChpDefault
{font-size:10.0pt}
@page WordSection1
{margin:72.0pt 90.0pt 72.0pt 90.0pt}
div.x_WordSection1
{}
-->
</style>
<div lang="ZH-CN" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="x_WordSection1">
<p class="x_MsoNormal" style=""><span lang="EN-US">Dear QE-Users,<br>
<br>
I am doing CPMD with cp.x for iron oxide system following the official tutorial which introduced two steps (Step1: Reaching the electronic ground state; Step 2: Relax the system) before doing cpmd.<br>
<br>
However, in both the electronic minimization and ionic minimization processes, I found that the force cannot converge (forc_conv_thr=1.0D-3 cannot be satisfied), while it is relative easy for electronic energy and total energy to converge (ekin_conv_thr=1.0D-6,
etot_conv_thr= 1.0D-4). <br>
<br>
I am wondering whether this is caused by that the initial crystal structure is not fully relaxed with the current pseudopotential. Is it necessary to do "vc-relax" with pw.x before doing CPMD with cp.x ? If yes, which parameters should be kept to be the same
in both pw.x and cp.x?<br>
<br>
Any comment or suggestion will be much appreciated.<br>
<br>
Thanks,<br>
Kejiang</span><b><span lang="EN-US" style="font-size:11.0pt"><br>
<br>
</span></b><span lang="EN-US">-------------------------------------------------<br>
</span><span lang="EN-US" style="font-size:9.0pt">Dr. Kejiang Li<br>
School of Metallurgical and Ecological Engineering,<br>
University of Science and Technology Beijing<br>
Xueyuan Rd., Haidian District, Beijing 100083, P. R. China</span><span lang="EN-US" style="font-size:7.5pt"><br>
<a href="https://kejiangli.com/" target="_blank">https://kejiangli.com/</a></span><span lang="EN-US" style="font-size:9.0pt"></span></p>
<p class="x_MsoNormal"><span lang="EN-US"> </span></p>
</div>
</div>
</body>
</html>