Cassandra 1.2 > Restarting Simulations
- piskuliche
- Topic Author
- Offline
- Junior Member
Less
More
- Posts: 30
- Thank you received: 1
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
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.
- ryangmullen
- Offline
- Administrator
Less
More
- Posts: 124
- Karma: 4
- Thank you received: 24
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
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
- Topic Author
- Offline
- Junior Member
Less
More
- Posts: 30
- Thank you received: 1
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
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.
- ryangmullen
- Offline
- Administrator
Less
More
- Posts: 124
- Karma: 4
- Thank you received: 24
- piskuliche
- Topic Author
- Offline
- Junior Member
Less
More
- Posts: 30
- Thank you received: 1
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
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.
- ryangmullen
- Offline
- Administrator
Less
More
- Posts: 124
- Karma: 4
- Thank you received: 24
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
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