Tutorial simple-adsorption-workflow
Prerequisites
- Install Apptainer
- Download the Apptainer image:
1 - Choose Simulation Parameters
A window with two tabs opens:
Select the desired parameters or fill them in manually
Modify advanced parameters
Save the parameter file
The structures come from the MOFX-DB database based on a server access pointing to the original structural database (CoRE MOF 2019). These structures are cleaned (no solvent, no disorder, etc.) from X-ray diffraction resolved structures from the CSD (Crystallographic Structural Database), where the identifier is a six-letter key.
Currently, partial charges can be determined by two methods:
- without partial charges:
None
- by charge equilibration
EQeq
(see this article)
The usable force fields are:
ExampleMOFsForceField
: a generic force field for MOFs, using a parameter set combining parameters from the Dreiding and UFF force fields.
2 - Run the Simulations
The basic script runs as many GCMC simulations using one CPU core each as there are combinations of input parameters. Example: 3 structures x 2 Temperatures x 5 Pressure points x 2 charge methods = 60 simulations.
Note: In its containerized version, the workflow cannot be used with a scheduler, so the user must ensure to launch no more simulations than there are accessible CPU cores to ensure acceptable performance.
The generated file structure is as follows:
The adsorption property database is located in the gcmc/run<index>.json
file, where <index>
is the experiment identifier.
3 - Update the Database
To update a database already generated by a past experiment (e.g., run<index1>.json
), a new database can be generated with the command:
This results in two new files:
run_merged.json
: the entire databaseisotherms.json
: the file containing the isotherms
The
isotherms.json
file does not contain all the metadata for each single-CPU simulation, unlike therun_merged.json
file, but it is very useful for grouping data and displaying it simply (see the next section).
4 - Visualize Results
Select the file with the isotherm data
Visualize the isotherms by selecting the desired parameters