RPMDMonteCarloBarostat not changing box volume?

The functionality of OpenMM will (eventually) include everything that one would need to run modern molecular simulation.
POST REPLY
User avatar
Nabil Faruk
Posts: 5
Joined: Mon Jul 09, 2012 11:33 am

RPMDMonteCarloBarostat not changing box volume?

Post by Nabil Faruk » Thu Jun 18, 2015 1:42 pm

Hi, I tried out the new RPMDMonteCarloBarostat with a ~(2 nm)^3 box of (q-SPC/Fw) water at 1 bar and various update freqs (5-50 steps) and platforms (Reference, OpenCL, CUDA). However, the volume reported when I call getState() on the RPMD integrator and then getPeriodicBoxVolume() on the state does not change during the course of the simulation. In contrast, the regular MonteCarloBarostat with the classical Langevin integrator produces volume changes with similar settings.

User avatar
Peter Eastman
Posts: 2607
Joined: Thu Aug 09, 2007 1:25 pm

Re: RPMDMonteCarloBarostat not changing box volume?

Post by Peter Eastman » Mon Jun 22, 2015 1:11 pm

Hi Nabil,

How long did you run your simulation for?

It's hard to guess at what might be going on without more information. I don't suppose you can post your script and input files?

Peter

POST REPLY