<div dir="ltr"><div>Dear all,</div><div>I am doing AIMD simulations using a cuboidal box with an initial edge length of 14.31 angstroms. After running 50 picoseconds in NPT ensemble (at thermostat 333 K temp. and barostat 1 bar pressure) the average pressure comes out to be around 1.2 bar but the instantaneous fluctuations in pressure at each time step are in the order of -10000 to 10000 bar. Please find the attached file having the plot of pressure at each time step. I have the following 3 questions-  <br></div><div><br></div><div>1. I want to know up to which order the instantaneous fluctuations in pressure are ok?</div><div>2. What can be done to decrease these instantaneous fluctuations at each time step?</div><div>3. Why there are 2 separate temperatures in the out file like one is thermostat temperature other is barostat temperature while the temperature part is only mentioned once in the input file?<br></div><div><br></div><div>TEMPERATURE [K]              =                     310.171              333.009</div><div>BAROSTAT TEMP[K]             =          0.652301869614E+03   0.334642361097E+03</div><div><br></div><div>Following is the cp2k code section of the input file:</div><div>______________________________________________<br></div><div>&MOTION<br>  &MD<br>   ENSEMBLE NPT_I<br>   TIMESTEP [fs] 0.5<br>   TEMPERATURE [K] 333<br>   STEPS 100000<br>   &BAROSTAT<br>      TIMECON [wavenumber_t] 1000.0<br>      PRESSURE  1<br>   &END BAROSTAT<br>   &THERMOSTAT<br>    TYPE NOSE<br>    REGION MASSIVE<br>    &NOSE<br>     LENGTH 3<br>     YOSHIDA 3<br>     TIMECON [wavenumber_t] 1000.0<br>     MTS 2<br>    &END NOSE<br>   &END THERMOSTAT</div><div>______________________________________________<br></div><div>Following is the outfile : *******************************************************************************<br> ENSEMBLE TYPE                =                                            NPT_I<br> STEP NUMBER                  =                                           107644<br> TIME [fs]                    =                                     53822.000000<br> CONSERVED QUANTITY [hartree] =                              -0.135935656268E+04<br><br>                                              INSTANTANEOUS             AVERAGES<br> CPU TIME [s]                 =                       45.30                14.39<br> ENERGY DRIFT PER ATOM [K]    =         -0.331821283044E+02  -0.299197457243E+02<br> POTENTIAL ENERGY[hartree]    =         -0.136035799755E+04  -0.136034247538E+04<br> KINETIC ENERGY [hartree]     =          0.452327402751E+00   0.485632690237E+00<br> TEMPERATURE [K]              =                     310.171              333.009<br> PRESSURE [bar]               =          0.256385237325E+04   0.129828680436E+01<br> BAROSTAT TEMP[K]             =          0.652301869614E+03   0.334642361097E+03<br> VOLUME[bohr^3]               =          0.188345152694E+05   0.189171690534E+05<br> CELL LNTHS[bohr]             =    0.2660632E+02   0.2660632E+02   0.2660632E+02<br> AVE. CELL LNTHS[bohr]        =    0.2664482E+02   0.2664482E+02   0.2664482E+02<br> *******************************************************************************<br> </div></div>