A process-level perspective of the impact of molecular force fields on the computational screening of MOFs for carbon capture

21 March 2023, Version 1
This content is a preprint and has not undergone peer review at the time of posting.


The question we pose in this study is to what extent the ranking of metal organic frameworks (MOFs) for adsorption-based carbon capture, and the selection of top performers identified in PSA process modelling, depends on the choice of the commonly available forcefields. To answer this question, we first generated distributions of CO2 and N2 adsorption isotherms via molecular simulation in 690 MOFs using six typical forcefields: the UFF or Dreiding sets of Lennard-Jones parameters, in combination with partial charges derived from ab initio calculations (DDEC scheme) or by charge equilibration (Qeq an EQeq schemes). We then conducted a systematic uncertainty quantification study using PSA process-level modelling. We observe that: (i) the ranking of MOFs significantly depends on the choice of forcefield; (ii) partial charge assignment is the prevailing source of uncertainty, and that charge equilibration schemes produce results which are inconsistent with ab initio-derived charges; iii) the choice of Lennard-Jones parameters is a considerable source of uncertainty. Our work highlights that is not really possible to obtain material rankings with high resolution using a single molecular modelling approach and that, as a minimum, some uncertainty should be estimated for the performance of MOFs shortlisted using high throughput computational screening workflows. Future prospects for computational screening studies are also discussed.


metal organic framework
carbon capture
computational screening
molecular simulation

Supplementary materials

Supplementary information
Forcefield parameters; mixture adsorption predictions and validation; additional PSA simulation details; process sensitivity analysis; machine learning surrogate model details; additional notes and figures to support the results and discussion; comparison of CALF-20 process-level performance with CRAFTED-u MOFs within the CL1 class.


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.