Reaction rates are a complicated function of molecular interactions, which can be selected from vast chemical design spaces. Seeking the design that optimizes a rate is a particularly challenging problem since the rate calculation for any one design is itself a difficult computation. Toward this end, we demonstrate a strategy based on transition path sampling to generate an ensemble of designs and reactive trajectories with a preference for fast reaction rates. Each step of the Monte Carlo procedure requires a measure of how a design constrains molecular configurations, expressed via the reciprocal of the partition function for the design. Although the reciprocal of the partition function would be prohibitively expensive to compute, we apply Booth’s method for generating unbiased estimates of a reciprocal of an integral to sample designs without bias. A generalization with multiple trajectories introduces a stronger preference for fast rates, pushing the sampled designs closer to the optimal design. We illustrate the methodology on two toy models of increasing complexity: escape of a single particle from a Lennard-Jones potential well of tunable depth and escape from a metastable tetrahedral cluster with tunable pair potentials.
Skip Nav Destination
Article navigation
28 November 2020
Research Article|
November 23 2020
Estimating reciprocal partition functions to enable design space sampling
Special Collection:
2020 JCP Emerging Investigators Special Collection
Alex Albaugh
;
Alex Albaugh
Department of Chemistry, Northwestern University
, 2145 Sheridan Road, Evanston, Illinois 60208, USA
Search for other works by this author on:
Todd R. Gingrich
Todd R. Gingrich
a)
Department of Chemistry, Northwestern University
, 2145 Sheridan Road, Evanston, Illinois 60208, USA
a)Author to whom correspondence should be addressed: todd.gingrich@northwestern.edu
Search for other works by this author on:
a)Author to whom correspondence should be addressed: todd.gingrich@northwestern.edu
Note: This paper is part of the JCP Emerging Investigators Special Collection.
J. Chem. Phys. 153, 204102 (2020)
Article history
Received:
August 14 2020
Accepted:
October 30 2020
Citation
Alex Albaugh, Todd R. Gingrich; Estimating reciprocal partition functions to enable design space sampling. J. Chem. Phys. 28 November 2020; 153 (20): 204102. https://doi.org/10.1063/5.0025358
Download citation file:
Sign in
Don't already have an account? Register
Sign In
You could not be signed in. Please check your credentials and make sure you have an active account and try again.
Sign in via your Institution
Sign in via your InstitutionPay-Per-View Access
$40.00