mikimo: microkinetic modeling and microkinetic volcano plots for homogeneous catalytic reactions

☄️ Why use microkinetic modeling? ☄️

Elegant way to deal with complex reaction pathway thermodynamics and kinetics.

Accounts for reaction conditions: temperature effects, concentration effects, reaction time, etc.

☄️ What are microkinetic volcano plots? ☄️

Volcano plot: diagrams that show the activity (or selectivity) of catalysts plotted against a descriptor variable that identifies a specific catalyst. Based on linear free energy scaling relationships.

Microkinetic volcano plot: volcano plots in which the activity/selectivity is expressed as the final product concentration, or a ratio of concentrations, after a given time.

Dependencies
Install
Usages
Examples
Limitations
Citation

Dependencies ↑

The code runs on pure python with the following dependencies:

numpy
scipy
autograd
matplotlib
pandas
scipy
h5py
fire
navicat_volcanic
openpyxl
tqdm

Install ↑

pip install .

python setup.py install

Usages ↑

The code requires two essential inputs for the analysis: energy data and reaction network (along with initial concentrations to be incorporated in the reaction network). All these files must be in the working directory or in the directory targeted with the -d option..

energy data: reaction_data (in csv or xlsx format)
reaction network: rxn_network (in csv or xlsx format)
(optional) kinetic data: kinetic_data (in csv or xlsx format)

The energy data must be named "reaction_data" and share similar format as in navicat_volcanic found here. In cases involving multiple reaction pathways, the final column of each pathway (including those pathways leading to the resting state) should contain "Prod" (case-insensitive).

The reaction network must be provided as a pandas-compatible csv or xlsx file named "rxn_network". Each row in the network represents an elementary step, while the columns represent the chemical species involved in the mechanism, excluding transition states. When filling in the reaction network, it is crucial to ensure that the species names match those in the energy data.

For each step (denoted as i) in the reaction network, assign a value of -n to the species on the left side of the elementary step equation and +n (or simply n) to the species on the right side, where n represents the stoichiometric coefficient. If a chemical species is not involved in step i, leave the corresponding cell empty or fill it with 0.

The initial concentrations should be specified in the last row of the reaction network file. This row can be named "c0", "initial_conc," or "initial conc."

The code offers three modes of operation:

mkm: This mode is for a single MKM run. If there are multiple profiles in the reaction data file, the top-most row is read.
cond: Use this mode for screening over reaction time and/or temperature.
vp: This mode is for screening over all energy profiles in the reaction data file. Note that it's only applicable when the reaction data contains more than one energy profile.

Once all input files are ready and mikimo is installed, several run options are available:

Call just microkinetic solver:

python -m navicat_mikimo mkm -d [DIR]

Microkinetic modeling for all reaction profiles:

python -m navicat_mikimo vp -d [DIR] -nd 0

To construct microkinetic volcano plot:

python -m navicat_mikimo vp -d [DIR] -nd 1

To construct microkinetic activity/selectivity map with a descriptor variable representing catalyst and temperature (or reaction time) (T1 and T2 are the lower and upper bounds of the temperature range, respectively):

python -m navicat_mikimo vp -d [DIR] -nd 1 -t T1 T2

To construct microkinetic activity/selectivity map:

python -m navicat_mikimo vp -d [DIR] -nd 2

To screen over reaction time and temperature:

python -m navicat_mikimo cond -d [DIR]

To smooth the volcano plots generated:

python replot.py [i]

Examples ↑

Microkinetic modeling for Pd-catalyzed carbocyclization-borylation of enallene in the presence of chiral phosphoric acid (298.15 K, 1 min):

python -m navicat_mikimo mkm -d test_cases/pd_carbocylic_borylation/ -t 298.15 -Tf 60

Microkinetic modeling for all profiles of the catalytic hydrosilylation of carbon dioxide with metal pincer complexes (323.15 K, 2 h):

python -m navicat_mikimo vp -d test_cases/pincer_CO2/ -t 323.15 -Tf 7200 -nd 0

Constructing microkinetic volcano plot for the catalytic hydrosilylation of carbon dioxide with metal pincer complexes (323.15 K, 2 h):

python -m navicat_mikimo vp -d test_cases/pincer_CO2/ -t 323.15 -Tf 7200 -nd 1 -ncore 24

Constructing microkinetic activity and selectivity maps for the catalytic hydrosilylation of carbon dioxide with metal pincer complexes (323.15 K, 2 h):

python -m navicat_mikimo vp -d test_cases/pincer_CO2/ -t 323.15 -Tf 7200 -nd 2 -ncore 24

Constructing microkinetic activity and selectivity maps with a descriptor variable representing catalyst and temperature [273.15-423.15 K] for the catalytic hydrosilylation of carbon dioxide with metal pincer complexes (2h):

python -m navicat_mikimo vp -d test_cases/pincer_CO2/ -t 273.15 423.15 -Tf 7200 -nd 1 -ncore 24

Constructing microkinetic activity and selectivity maps with reaction time [2-24 hr] and temperature [273.15-423.15 K] as descriptors for the catalytic hydrosilylation of carbon dioxide with the Co pincer complex:

python -m navicat_mikimo cond -d test_cases/pincer_CO2_jacs/ -tt -Tf 7200 86400 -t 273.15 423.15 -ncore 24

Using filtering to smooth the plot:

python -m navicat_mikimo replot examples/data/vp/data_a.h5 -p 3 3 3 -w 20 20 20

You can find examples demonstrating how to read h5 files and regenerate plots in the "examples" folder.

If the kinetic profile is detected in the directory, the code will prompt the user with an option to choose between using the kinetic profile instead of the energy profile. However, it's important to note that selecting the kinetic profile will restrict the user from screening over a range of temperatures or utilizing different temperature settings. Additionally, using the kinetic profile will limit the information about species names (for ax labelling purposes) and may worsen the quality of the linear scaling relationships in volcano plot generation.

Limitations ↑

Overlapping states of different pathways before the reference state (starting point).
Bridging states between otherwise separate pathways.
Different TSs connecting the same 2 intermediates: just choose the lowest one or compute an effective TS energy that corresponds to the sum of the two rate constants.

To overcome these limitations and offer more flexibility, users have the option to input a kinetic profile named "kinetic_profile" in either csv or xlsx format, replacing the conventional energy profile. However, this choice comes with the trade-off of disabling the ability to screen over a range of temperatures or use different temperature settings.

Citation ↑

If you use navicat_mikimo in your work, please cite our work and the forthcoming publication.

lcmd-epfl/mikimo