<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="arial, sans-serif">Dear QE Developers,</font><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">I am running a constrained optimization of a slab, with two constrained atoms, using DAMP, but the constraint does not consistently hold.  </font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">In my example (attached), the constraint on one of the two atoms holds during many geometry steps.  The constraint on the other atom holds for 9 iterations, and then moves on the 10th.  I have run similar jobs on multiple architectures now, and the pattern remains- one atom stays fixed by the constraint, and the other does not always.  </font></div><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif"><br></font></span></div><div><font face="arial, sans-serif">Here is a grep of the position of the copper atom that moves:<span style="color:rgb(0,0,0)"></span><br></font></div><div><font face="arial, sans-serif"><br></font></div><div>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.625170579    0   0   0</font></span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif">NCu      0.500000000   0.500000000   0.624738221    0   0   0</font></span></p><p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif"><br></font></span></p><p style="margin:0px;font-stretch:normal;line-height:normal"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif"><font color="#000000">(I originally thought that there was a symmetry constraint that was causing the problem, <span style="caret-color: rgb(0, 0, 0);">which</span> is why I created atoms called "NCu", the Cu atoms I'm constraining.  )</font></font></span></p><p style="margin:0px;font-stretch:normal;line-height:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font face="arial, sans-serif"><br></font></span></p></div><div><font face="arial, sans-serif"> I have attached the input file + pseudopotentials, and the output file.<br></font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">As you can see from the output file, this is for <span style="color:rgb(0,0,0)">PWSCF v.6.3, and a single node job with 64 processors, but this error also occurs on other clusters, and parallelized across multiple nodes.</span></font></div><div><font face="arial, sans-serif"><span style="color:rgb(0,0,0)"><br></span></font></div><div><font face="arial, sans-serif"><font color="#000000">Thanks, and please <span style="caret-color: rgb(0, 0, 0);">let</span> me know if additional information would be helpful.</font></font></div><div><div><br></div><div>Best,</div><div>Katie </div><div><br></div><div>Katie Schwarz</div><div>Research Chemist</div></div><div>National Institute of Standards and Technology</div></div></div></div></div>