I have a quick question. Is there a way to use langevin dynamics (thermostat), keeping volume as independent and adjusting parameter (i.e., volume change during the course of the simulation)? In the doc file, it is suggested to use nvt or nve with fix langevin in which the volume remain constant. Or before running such simulation, I need to equilibrate the volume using some other fixes??