<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=utf-8">
<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:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.E-MailFormatvorlage19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
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="DE" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US">Hi Torstein<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US">Did you select the thermostat TYPE CSVR? I do not see such a line in your input file. The default TYPE is NOSE and thus your settings in the &CSVR section are possibly not considered,
but the default settings of the &NOSE section, e.g. TIMECON 1000.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US">HTH<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US">Matthias<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><b>Von:</b> cp...@googlegroups.com <...@googlegroups.com> <b>
Im Auftrag von </b>Torstein Fjermestad<br>
<b>Gesendet:</b> Mittwoch, 3. April 2019 18:56<br>
<b>An:</b> cp2k <...@googlegroups.com><br>
<b>Betreff:</b> [CP2K:11528] The temperature of an NVT simulation is decreasing from 673 K to approximately 350 K<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Dear cp2k users,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I have done a NVT simulation with the temperature set to 673 K<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">The idea was to keep the temperature around this value by the use of the CSVR thermostat.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">However, the temperature quickly decreases to around 350 K as shown in the attached figure.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">The MD section of the input file is as follows:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"> &MD<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ENSEMBLE NVT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> TEMPERATURE [K] 673<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> TIMESTEP [fs] 0.30<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &THERMOSTAT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &CSVR<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> TIMECON [fs] 10<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END CSVR<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END THERMOSTAT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> STEPS 100000<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &PRINT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &PROGRAM_RUN_INFO<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &EACH<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> MD 20<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END EACH<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END PROGRAM_RUN_INFO<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END PRINT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END MD<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Why is the temperature decreasing?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Shouldn't the thermostat be able to bring the temperature up again to the initial value?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Is there something I do incorrectly in the set-up of this simulation?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks a lot for your help. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Regards,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Torstein<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<p class="MsoNormal">-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to
<a href="mailto:cp2k+unsu...@googlegroups.com">cp2k+unsu...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k">https://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout">https://groups.google.com/d/optout</a>.<o:p></o:p></p>
</div>
</body>
</html>