<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p></p>
<div dir="ltr" style="font-family: Calibri, Helvetica, sans-serif; font-size: 16px;">
<div>Dear DFTB+ programmers</div>
<div><br>
</div>
<div>I'm a new user of DFTB+ (using version 18.2) and I'm trying to use NPT molecular dynamics on some polimeryzed C60 phases but I'm having some problems. <br>
</div>
<div><br>
</div>
<div>Using some particular cells the atoms start getting out of the cell  even if the tag "KeepStationary = Yes" is present. A lot of times when this happens the final temperature is higher than the given threshold. I  tried different thermostats however the
 outcome was always the same (so it shouldn't be the ice cube effect). Although if I stop the calculation at some step, then translate the atoms back to the cell (making all fractional positions to be between 0 and 1), and then continue the calculation (without
 giving the velocities and charges) it will go as it should (making the temperature profile as required). Furthermore I tried to restart the calculation without changing the positions to be sure that it was the problem and indeed the threshold of temperature
 was again overshoot (I'm attaching this input). Normally this overshoot of temperature happens with a huge deformations of the cell. If you could elucidate me on why the keep stationary tag is not working properly and on whats happening afterwards it would
 be much appreciated.<br>
</div>
<div><br>
</div>
<div>Besides this when I try different description of the system (all of them with the same number of atoms) the transformation is deferent. Should this happen since it is exactly the same structure just described in a different cell?</div>
<div><br>
</div>
<div>Finally I would like to ask how is the pressure applied to the cell. It seems that if we give XXgpa of pressure it applies XXgpa in each lattice vector although in the provided literature it isn't clear how the program works.  Furthermore I would like
 to request the inclusion of a pressure matrix in future versions.   <br>
</div>
<div><br>
</div>
<div> Best Regards,</div>
<div><br>
</div>
<div>Jorge Laranjeira</div>
<div><br>
</div>
</div>
<br>
<p></p>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Jorge Laranjeira<br>
<b>Sent:</b> Monday, March 4, 2019 6:16 PM<br>
<b>To:</b> dftb-plus-user@mailman.zfn.uni-bremen.de<br>
<b>Subject:</b> DFTB+ 18.2 NPT atoms going out of cell</font>
<div> </div>
</div>
<div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p></p>
<div dir="ltr">
<div>Dear DFTB+ programmers</div>
<div><br>
</div>
<div>I'm a new user of DFTB+ (using version 18.2) and I'm trying to use NPT molecular dynamics on some polimeryzed C60 phases but I'm having some problems.
<br>
</div>
<div><br>
</div>
<div>Using some particular cells the atoms start getting out of the cell  even if the tag "KeepStationary = Yes" is present. A lot of times when this happens the final temperature is higher than the given threshold. I  tried different thermostats however the
 outcome was always the same (so it shouldn't be the ice cube effect). Although if I stop the calculation at some step, then translate the atoms back to the cell (making all fractional positions to be between 0 and 1), and then continue the calculation (without
 giving the velocities and charges) it will go as it should (making the temperature profile as required). Furthermore I tried to restart the calculation without changing the positions to be sure that it was the problem and indeed the threshold of temperature
 was again overshoot (I'm attaching this input). Normally this overshoot of temperature happens with a huge deformations of the cell. If you could elucidate me on why the keep stationary tag is not working properly and on whats happening afterwards it would
 be much appreciated.<br>
</div>
<div><br>
</div>
<div>Besides this when I try different description of the system (all of them with the same number of atoms) the transformation is deferent. Should this happen since it is exactly the same structure just described in a different cell?</div>
<div><br>
</div>
<div>Finally I would like to ask how is the pressure applied to the cell. It seems that if we give XXgpa of pressure it applies XXgpa in each lattice vector although in the provided literature it isn't clear how the program works.  Furthermore I would like
 to request the inclusion of a pressure matrix in future versions.   <br>
</div>
<div><br>
</div>
<div> Best Regards,</div>
<div><br>
</div>
<div>Jorge Laranjeira</div>
</div>
<br>
<p></p>
</div>
</div>
</div>
</div>
</body>
</html>