News:

SMF - Just Installed!

Main Menu

Anisotropic volume change with Monte Carlo

Started by MichaelH, April 18, 2019, 08:31:52 PM

Previous topic - Next topic

MichaelH

Dear David,

Im looking at the advanced examples 1 and 2, trying to adapt them to study swelling (in the z-direction) of multi-layered clay(ff) system. Since my clay layers are periodic in x and y, I would like to avoid isotropic volume changes and use some sort of anisotropic-like barostatting instead.

So my question is: When running MonteCarlo, does the volume move (VolumeChangeProbability) only do isotropic volume scaling? or is there some way to activate anisotropic volume trial moves? or simply just do volume moves along the z-direction? I have seen the BoxShapeChangeProbability option, but Im not sure if its effective.

Best regards,

Michael

David Dubbeldam

You could change the source code to achieve that.
Look at the function 'int VolumeMove(void)' in mc_moves.c and change the isotropic scaling into only scaling in the z-direction.