Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Messages - isaaraujo

Pages: [1]
1
General / NPT Simulations
« on: December 01, 2021, 05:25:16 PM »
Dear Dr. Dubbeldam,

Is it possible to run NPT simulations in a rigid framework ( in my case it is a clay slit pore)? When I run the NPT simulations keeping framework rigid, I see that I have two simulation boxes: One for the framework(which is a slit pore) and one for the the adsorbate(water molecules). Only the simulation cell of the water molecules change volume, the simulation cell of the adsorbent is kept fixed, therefore when there is change in volume in the water molecules simulation cell, the framework and molecules get very close and the simulation crashes. Is there a way to keep framework rigid but allow change in the framework simulation box (specially in the z direction)?

Thank you
Isa

2
General / Re: NPT-Flexible framework
« on: December 01, 2021, 04:29:31 PM »
Thank you very much for your response!
Do you know if it is possible to run NPT simulations with a rigid framework? When I run the NPT simulations keeping framework rigid, I see that I have two simulation boxes: One for the framework(which is a slit pore) and one for the the adsorbate(water molecules). Only the simulation cell of the water molecules change volume, the simulation cell of the adsorbent is kept fixed, therefore when there is change in volume in the water molecules simulation cell, the framework and molecules get very close and the simulation crashes. Is there a way to keep framework rigid but allow change in the framework simulation box?
Thank you

3
General / NPT-Flexible framework
« on: November 11, 2021, 08:37:01 AM »
Hello,
I am performing MD in the NPT ensemble. The system is an illite slit pore filled with water and an electrolyte. I am considering the framework flexible (the only bond considered in the framework is between hydrogen and oxygen). When I run NPT simulations the structure moves out of the simulation box and looks pretty deformed.
If I consider the framework to be rigid in the NPT ensemble the simulation crashes even at very small time steps and the output file for the framework (.cif) shows NAN values.
What could be a reason for that, and is there any way to address this problem in RASPA? Is there a way to restrain all hydrogen bonds in the framework?

#CoreShells bond  BondDipoles UreyBradley bend  inv  tors improper-torsion bond/bond bond/bend bend/bend stretch/torsion bend/torsion
          0    1            0           0    0    0     0                0         0         0         0               0            0
#bond stretch atom n1-n2, equilibrium distance, bondforce-constant
OH  HO  HARMONIC_BOND     557706.438693      1.00       
Thank you

Pages: [1]