Welcome To Support Community

Materials Studio

Advanced Search
Ask Search:
farzinfarzin 

Why dosn't the cell length change under pressure?

Hi all,

I simulated a periodic box including a graphene sheet at the middle (Fig. is attached). I did some NPT dynamics under the atmosphere pressure, but the length of the box didn't reduced in C direction (normal to graphene). This also happens when some polymer chains are packed into the empty space of the box with low density (e.g. 0.5 g/cm3). What is the reason?


User-added image

Thank you in advance.
Farzin
Best Answer chosen by farzin
CarstenCarsten (Accelrys) 
Please use Souza-Martin or Parrinello as barostat.

Your system is quite anisotropic and Anderson and Berendsen barostats are not a good choice then.

From the help for Berendsen:
As implemented, this method changes the cell uniformly, so that the size of the cell is changed, but not its shape. Therefore, it is inappropriate for simulations such as crystal phase transitions, where both the cell size and shape are expected to change.

And for Andersen:
However, this method is not very useful for studying materials under anisotropic stress or undergoing phase transitions, which involve changes in both cell lengths and cell angles

All Answers

CarstenCarsten (Accelrys) 
Please use Souza-Martin or Parrinello as barostat.

Your system is quite anisotropic and Anderson and Berendsen barostats are not a good choice then.

From the help for Berendsen:
As implemented, this method changes the cell uniformly, so that the size of the cell is changed, but not its shape. Therefore, it is inappropriate for simulations such as crystal phase transitions, where both the cell size and shape are expected to change.

And for Andersen:
However, this method is not very useful for studying materials under anisotropic stress or undergoing phase transitions, which involve changes in both cell lengths and cell angles
This was selected as the best answer
farzinfarzin
Dear Carsten,
Thank you for your helpful answer.

Farzin