Author Topic: MIL-53(Cr) flexible MOF MD problem  (Read 796 times)

merdan

  • Newbie
  • *
  • Posts: 2
    • View Profile
MIL-53(Cr) flexible MOF MD problem
« on: April 20, 2020, 01:32:50 PM »
Hi everyone, I have a problem with MD simulation of water diffusion in MIL-53(Cr) flexible MOF. I want to reproduce some simulation results that i read in the some article.  In the article, it is mentioned that diffusion of 5 water molecule (Tip4p2005 model) per unit cell of MIL-53(Cr) flexible MOF will decrease volume of MOF from 1480 A^3 to around 1050 A^3. Simulation conditions are given in the article as follows,

simulation type MD
ensemble NPT
temperature 300K
timestep 0.001
total time 2ns

I used forcefield mixing rules and forcefield definition files which i found in the raspa examples files for MIL-53Cr MOF. But after simulation finished, there is no observable reduction in volume, it changed very little around initial volume 1480 A^3. I don't know what is the problem, i will share my input files. Can you give some idea? thanks.

link for files : https://drive.google.com/open?id=1hsrNft_pGlJzLnRTQcEGRceDNCYNLBtC
« Last Edit: April 20, 2020, 01:36:41 PM by merdan »

David Dubbeldam

  • Full Member
  • ***
  • Posts: 141
    • View Profile
Re: MIL-53(Cr) flexible MOF MD problem
« Reply #1 on: April 24, 2020, 07:26:57 PM »
I would do this as a function of loading. Around 2 waters per unit cell many models start to collapse (the same for CO2).
So five water molecule might be so high that the unit cell starts to become larger again.