Cassandra 1.2 > Restarting Simulations

  • piskuliche
  • piskuliche's Avatar Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
7 years 9 months ago #168 by piskuliche
Ryan,

Thanks for your help a few weeks ago. Is there some change we need to make with the MCF file now that we are using 1.2 (versus 1.1) i.e. can we use a 1.1 MCF for a 1.2 simulation or do we need to convert them? I ask this because in my restarts of 1.1 simulations using 1.2 I am seeing the vapor box collapsing (i.e. the volume drops by two orders of magnitude immediately upon restarting a 1.1 simulation with 1.2 using the same input parameters etc). This is happening whenever I use the read_config or checkpoint commands. Do you have any idea what would cause this? (I can send you my files if that would be helpful.)

Thanks again for all of your help,
Zeke

Please Log in to join the conversation.

More
7 years 9 months ago #169 by ryangmullen
Hi Zeke,

1.2 can read Mie parameters from the MCF, while 1.1 will only read LJ parameters. But any MCF that worked with 1.1 should also work with 1.2.

I can't think what would be causing the vapor box to collapse like that. Please upload your files and we'll try to straighten this out.

Ryan

Please Log in to join the conversation.

  • piskuliche
  • piskuliche's Avatar Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
7 years 9 months ago - 7 years 9 months ago #171 by piskuliche
Hello Ryan,

Google Drive Link : drive.google.com/open?id=0B35lw-A81eT_Ti1ncElCMUFWY1E

Thanks for the update. I have included a zip folder with this email that should contain all the information needed (I hope).

The top level folder contains:
* The MCF files for my system
* A graph indicating the problem (Volume vs. Sweeps) - red line is the original simulation, green line is the restarted simulation from the checkpoint.
* The two fragment libraries
* Three Sub directories.

The Subdirectory 'Original_Run' includes all of the Cassandra files from the simulation being restarted.
The Subdirectory 'Checkpoint_Cont' includes all of the Cassandra files from the simulation restarted from the checkpoint
The Subdirectory 'Read_Config_Cont' includes all of the Cassandra files from the simulation restarted from the last configuration of the original system (i.e. by reading in the xyz file)

Thanks!
Zeke
Last Edit: 7 years 9 months ago by piskuliche.

Please Log in to join the conversation.

More
7 years 9 months ago #181 by ryangmullen
Hi Zeke,

There's definitely something buggy in the v1.2 GEMC_NPT routines. It looks like both boxes are converging to the liquid composition / volume. We are looking for the root cause now.

Ryan

Please Log in to join the conversation.

  • piskuliche
  • piskuliche's Avatar Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
7 years 9 months ago #182 by piskuliche
Ryan-

Good thing to know I am not crazy. I appreciate the response and look forward to hearing from you. On a side note- I believe I am seeing something similar restarting my single component GEMC simulations as well (GEMC_NVT) but in this case the liquid is collapsing and the vapor box is expanding. Thanks for your help!

Thanks,
Zeke

Please Log in to join the conversation.

More
7 years 9 months ago #185 by ryangmullen
Hi Zeke,

We found the bug. A variable in the acceptance criteria was being erroneously reset in the particle transfer routine. This effected both NVT and NPT GEMC simulations, so I hope it solves both the problems you were seeing. I tested it on the methanol/CO2 mixture simulations you shared and now 1.2 gives the same composition and vapor density as you were seeing with 1.1.

You can download the corrected code now.

Ryan
The following user(s) said Thank You: piskuliche

Please Log in to join the conversation.

Time to create page: 0.131 seconds