[DFTB-Plus-User] Does Barostat really work?

范芫钟 fyanzh at mail2.sysu.edu.cn
Wed May 6 18:42:10 CEST 2015


Dear Sirs,
       I'm try to do some MD research of NPT ensemble, with the detailed information of Drive listed below. (DFTB+ Version 1.2)
      Driver = {
  VelocityVerlet = {
  OutputPrefix = "XXX"
  Steps = 100000
  TimeStep[fs]= 0.5
  KeepStationary = Yes
  MDRestartFrequency = 10
  Thermostat = Berendsen {
  Temperature[K] = 298.15
  CouplingStrength[ps] = 0.01
}
  Barostat = {
  Pressure[Pa] = 101325.0
  TimeScale[ps] = 0.1
}
}
}


       The system is periodic, with the LatticeVector is,
11.16658929  -4.884981308e-15  12.61723333
-5.583294646  9.67055  12.61723333
-5.583294646  -9.67055  12.61723333

      However, the LatticeVector doesn't change in .gen, and the output in .dftb shows that the Mode is NVT.
Mode:                        MD with scaling of velocities according to temperature
                             (a.k.a. not NVT ensemble using Berendsen thermostating)

     I know that if I employ Nose-Hoover, it'll be a standard NVT, but it should be NPT instead.
    The similar output appears when I want to do the research of NPH, and it becomes NVE.
    Therefore, I doubt that whether I have some mistakes in Drive or the Barostat doesn't work at all?
    I'm very grateful to receive any response from you!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20150507/752c5efa/attachment.htm>


More information about the DFTB-Plus-User mailing list