1
|
Cannon WR, Britton S, Banwarth-Kuhn M, Alber M. Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models. J Chem Phys 2024; 160:214123. [PMID: 38842085 DOI: 10.1063/5.0180417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 05/13/2024] [Indexed: 06/07/2024] Open
Abstract
We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.
Collapse
Affiliation(s)
- William R Cannon
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
- Department of Mathematics, University of California, Riverside, California 92505, USA
- Center for Quantitative Modeling in Biology, University of California Riverside, Riverside, California 92505, USA
| | - Samuel Britton
- Department of Mathematics, University of California, Riverside, California 92505, USA
- Center for Quantitative Modeling in Biology, University of California Riverside, Riverside, California 92505, USA
| | - Mikahl Banwarth-Kuhn
- Center for Quantitative Modeling in Biology, University of California Riverside, Riverside, California 92505, USA
- Department of Mathematics, California State University East Bay, Hayward, California 94542, USA
| | - Mark Alber
- Department of Mathematics, University of California, Riverside, California 92505, USA
- Center for Quantitative Modeling in Biology, University of California Riverside, Riverside, California 92505, USA
| |
Collapse
|
2
|
Schöller A, Woodcock HL, Boresch S. Exploring Routes to Enhance the Calculation of Free Energy Differences via Non-Equilibrium Work SQM/MM Switching Simulations Using Hybrid Charge Intermediates between MM and SQM Levels of Theory or Non-Linear Switching Schemes. Molecules 2023; 28:4006. [PMID: 37241747 PMCID: PMC10222338 DOI: 10.3390/molecules28104006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Revised: 04/25/2023] [Accepted: 04/26/2023] [Indexed: 05/28/2023] Open
Abstract
Non-equilibrium work switching simulations and Jarzynski's equation are a reliable method for computing free energy differences, ΔAlow→high, between two levels of theory, such as a pure molecular mechanical (MM) and a quantum mechanical/molecular mechanical (QM/MM) description of a system of interest. Despite the inherent parallelism, the computational cost of this approach can quickly become very high. This is particularly true for systems where the core region, the part of the system to be described at different levels of theory, is embedded in an environment such as explicit solvent water. We find that even for relatively simple solute-water systems, switching lengths of at least 5 ps are necessary to compute ΔAlow→high reliably. In this study, we investigate two approaches towards an affordable protocol, with an emphasis on keeping the switching length well below 5 ps. Inserting a hybrid charge intermediate state with modified partial charges, which resembles the charge distribution of the desired high level, makes it possible to obtain reliable calculations with 2 ps switches. Attempts using step-wise linear switching paths, on the other hand, did not lead to improvement, i.e., a faster convergence for all systems. To understand these findings, we analyzed the solutes' properties as a function of the partial charges used and the number of water molecules in direct contact with the solute, and studied the time needed for water molecules to reorient themselves upon a change in the solute's charge distribution.
Collapse
Affiliation(s)
- Andreas Schöller
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Währingerstr. 17, A-1090 Vienna, Austria
- Vienna Doctoral School in Chemistry (DoSChem), University of Vienna, Währingerstr. 42, A-1090 Vienna, Austria
| | - H. Lee Woodcock
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA;
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Währingerstr. 17, A-1090 Vienna, Austria
| |
Collapse
|
3
|
Blaber S, Sivak DA. Skewed thermodynamic geometry and optimal free energy estimation. J Chem Phys 2020; 153:244119. [PMID: 33380076 DOI: 10.1063/5.0033405] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
Free energy differences are a central quantity of interest in physics, chemistry, and biology. We develop design principles that improve the precision and accuracy of free energy estimators, which have potential applications to screening for targeted drug discovery. Specifically, by exploiting the connection between the work statistics of time-reversed protocol pairs, we develop near-equilibrium approximations for moments of the excess work and analyze the dominant contributions to the precision and accuracy of standard nonequilibrium free-energy estimators. Within linear response, minimum-dissipation protocols follow the geodesics of the Riemannian metric induced by the Stokes friction tensor. We find that the next-order contribution arises from the rank-3 supra-Stokes tensor that skews the geometric structure such that minimum-dissipation protocols follow the geodesics of a generalized cubic Finsler metric. Thus, near equilibrium, the supra-Stokes tensor determines the leading-order contribution to the bias of bidirectional free-energy estimators.
Collapse
Affiliation(s)
- Steven Blaber
- Department of Physics, Simon Fraser University, Burnaby, British Columbia V5A 1S6, Canada
| | - David A Sivak
- Department of Physics, Simon Fraser University, Burnaby, British Columbia V5A 1S6, Canada
| |
Collapse
|
4
|
Sun Z, Wang X, Song J. Extensive Assessment of Various Computational Methods for Aspartate's pK a Shift. J Chem Inf Model 2017. [PMID: 28644624 DOI: 10.1021/acs.jcim.7b00177] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
A series of computational methods for pKa shift prediction are extensively tested on a set of benchmark protein systems, aiming at identifying pitfalls and evaluating their performance on high variants. Including 19 ASP residues in 10 protein systems, the benchmark set consists of both residues with highly shifted pKa values as well as those varying little from the reference value, with an experimental RMS free energy differences of 2.49 kcal/mol with respect to blocked amino acid, namely the RMS pKa shift being 1.82 pKa units. The constant pH molecular dynamics (MD), alchemical methods, PROPKA3.1, and multiconformation continuum electrostatics give RMSDs of 1.52, 2.58, 1.37, and 3.52 pKa units, respectively, on the benchmark set. The empirical scoring method is the most accurate one with extremely low computational cost, and the pH-dependent model is also able to provide accurate results, while the accuracy of MD sampling incorporating alchemical free energy simulation is prohibited by convergence achievement and the performance of conformational search incorporating multiconformation continuum electrostatics is bad. Former research works did not define statistical uncertainty with care and yielded the questionable conclusion that alchemical methods perform well in most benchmarks. In this work the traditional alchemical methods are thoroughly tested for high variants. We also performed the first application of nonequilibrium alchemical methods to the pKa cases.
Collapse
Affiliation(s)
- Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Physics and Material Science, East China Normal University , Shanghai 200062, China
| | - Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Material Science, East China Normal University , Shanghai 200062, China
| | - Jianing Song
- NYU-ECNU Center for Computational Chemistry, NYU Shanghai , Shanghai 200062, China.,School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China
| |
Collapse
|
5
|
Sun ZX, Wang XH, Zhang JZH. BAR-based optimum adaptive sampling regime for variance minimization in alchemical transformation. Phys Chem Chem Phys 2017; 19:15005-15020. [DOI: 10.1039/c7cp01561e] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The efficiency of alchemical free energy simulations with a staging strategy is improved by adaptively manipulating the significance of each ensemble followed by importance sampling.
Collapse
Affiliation(s)
- Zhao X. Sun
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- East China Normal University
- Shanghai 200062
- China
| | - Xiao H. Wang
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- East China Normal University
- Shanghai 200062
- China
| | - John Z. H. Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai
- Shanghai 200062
- China
- School of Chemistry and Molecular Engineering
- East China Normal University
| |
Collapse
|
6
|
Gong Z, Lan Y, Quan HT. Stochastic Thermodynamics of a Particle in a Box. PHYSICAL REVIEW LETTERS 2016; 117:180603. [PMID: 27835006 DOI: 10.1103/physrevlett.117.180603] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2016] [Indexed: 06/06/2023]
Abstract
The piston system (particles in a box) is the simplest paradigmatic model in traditional thermodynamics. However, the recently established framework of stochastic thermodynamics (ST) fails to apply to this model system due to the embedded singularity in the potential. In this Letter, we study the ST of a particle in a box by adopting a novel coordinate transformation technique. Through comparing with the exact solution of a breathing harmonic oscillator, we obtain analytical results of work distribution for an arbitrary protocol in the linear response regime and verify various predictions of the fluctuation-dissipation relation. When applying to the Brownian Szilard engine model, we obtain the optimal protocol λ_{t}=λ_{0}2^{t/τ} for a given sufficiently long total time τ. Our study not only establishes a paradigm for studying ST of a particle in a box but also bridges the long-standing gap in the development of ST.
Collapse
Affiliation(s)
- Zongping Gong
- School of Physics, Peking University, Beijing 100871, China
- Department of Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Yueheng Lan
- Department of Physics, Beijing University of Posts and Telecommunications, Beijing 100876, China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
| | - H T Quan
- School of Physics, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
| |
Collapse
|
7
|
Rotskoff GM, Crooks GE. Optimal control in nonequilibrium systems: Dynamic Riemannian geometry of the Ising model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:060102. [PMID: 26764609 DOI: 10.1103/physreve.92.060102] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2015] [Indexed: 06/05/2023]
Abstract
A general understanding of optimal control in nonequilibrium systems would illuminate the operational principles of biological and artificial nanoscale machines. Recent work has shown that a system driven out of equilibrium by a linear response protocol is endowed with a Riemannian metric related to generalized susceptibilities, and that geodesics on this manifold are the nonequilibrium control protocols with the lowest achievable dissipation. While this elegant mathematical framework has inspired numerous studies of exactly solvable systems, no description of the thermodynamic geometry yet exists when the metric cannot be derived analytically. Herein, we numerically construct the dynamic metric of the two-dimensional Ising model in order to study optimal protocols for reversing the net magnetization.
Collapse
Affiliation(s)
- Grant M Rotskoff
- Biophysics Graduate Group, University of California, Berkeley, California 94720, USA
| | - Gavin E Crooks
- Molecular Biophysics Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Kavli Energy NanoSciences Institute, Berkeley, California 94720, USA
| |
Collapse
|
8
|
Acconcia TV, Bonança MVS. Degenerate optimal paths in thermally isolated systems. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:042141. [PMID: 25974472 DOI: 10.1103/physreve.91.042141] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/13/2015] [Indexed: 06/04/2023]
Abstract
We present an analysis of the work performed on a system of interest that is kept thermally isolated during the switching of a control parameter. We show that there exists, for a certain class of systems, a finite-time family of switching protocols for which the work is equal to the quasistatic value. These optimal paths are obtained within linear response for systems initially prepared in a canonical distribution. According to our approach, such protocols are composed of a linear part plus a function which is odd with respect to time reversal. For systems with one degree of freedom, we claim that these optimal paths may also lead to the conservation of the corresponding adiabatic invariant. This points to an interesting connection between work and the conservation of the volume enclosed by the energy shell. To illustrate our findings, we solve analytically the harmonic oscillator and present numerical results for certain anharmonic examples.
Collapse
Affiliation(s)
- Thiago V Acconcia
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, 13083-859 Campinas, São Paulo, Brazil
| | - Marcus V S Bonança
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, 13083-859 Campinas, São Paulo, Brazil
| |
Collapse
|
9
|
Bonança MVS, Deffner S. Optimal driving of isothermal processes close to equilibrium. J Chem Phys 2014; 140:244119. [DOI: 10.1063/1.4885277] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Marcus V. S. Bonança
- Instituto de Física “Gleb Wataghin,” Universidade Estadual de Campinas, 13083-859, Campinas, São Paulo, Brazil
- Department of Chemistry and Biochemistry and Institute of Physical Sciences and Technology, University of Maryland, College Park, Maryland 20742, USA
| | - Sebastian Deffner
- Department of Chemistry and Biochemistry and Institute of Physical Sciences and Technology, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
10
|
Andresen B. Aktuelle Trends in der Thermodynamik in endlicher Zeit. Angew Chem Int Ed Engl 2011. [DOI: 10.1002/ange.201001411] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Bjarne Andresen
- Niels‐Bohr‐Institut, Universität Kopenhagen, Universitetsparken 5, DK‐2100 Kopenhagen (Dänemark) http://www.fys.ku.dk/∼andresen/BAhome/welcome.html
| |
Collapse
|
11
|
Andresen B. Current trends in finite-time thermodynamics. Angew Chem Int Ed Engl 2011; 50:2690-704. [PMID: 21374763 DOI: 10.1002/anie.201001411] [Citation(s) in RCA: 98] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2010] [Indexed: 11/08/2022]
Abstract
The cornerstone of finite-time thermodynamics is all about the price of haste and how to minimize it. Reversible processes may be ultimately efficient, but they are unrealistically slow. In all situations-chemical, mechanical, economical-we pay extra to get the job done quickly. Finite-time thermodynamics can be used to develop methods to limit that extra expenditure, be it in energy, entropy production, money, or something entirely different. Finite-time thermodynamics also includes methods to calculate the optimal path or mode of operation to achieve this minimal expenditure. The concept is to place the system of interest in contact with a time-varying environment which will coax the system along the desired path, much like guiding a horse along by waving a carrot in front of it.
Collapse
|
12
|
|
13
|
Minh DDL, Chodera JD. Optimal estimators and asymptotic variances for nonequilibrium path-ensemble averages. J Chem Phys 2009; 131:134110. [PMID: 19814546 PMCID: PMC2771048 DOI: 10.1063/1.3242285] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2009] [Accepted: 09/13/2009] [Indexed: 11/14/2022] Open
Abstract
Existing optimal estimators of nonequilibrium path-ensemble averages are shown to fall within the framework of extended bridge sampling. Using this framework, we derive a general minimal-variance estimator that can combine nonequilibrium trajectory data sampled from multiple path-ensembles to estimate arbitrary functions of nonequilibrium expectations. The framework is also applied to obtain asymptotic variance estimates, which are a useful measure of statistical uncertainty. In particular, we develop asymptotic variance estimates pertaining to Jarzynski's equality for free energies and the Hummer-Szabo expressions for the potential of mean force, calculated from uni- or bidirectional path samples. These estimators are demonstrated on a model single-molecule pulling experiment. In these simulations, the asymptotic variance expression is found to accurately characterize the confidence intervals around estimators when the bias is small. Hence, the confidence intervals are inaccurately described for unidirectional estimates with large bias, but for this model it largely reflects the true error in a bidirectional estimator derived by Minh and Adib.
Collapse
Affiliation(s)
- David D L Minh
- Laboratory of Chemical Physics, NIDDK, National Institutes of Health, Bethesda, Maryland 20892, USA.
| | | |
Collapse
|
14
|
Minh DDL. Density-dependent analysis of nonequilibrium paths improves free energy estimates. J Chem Phys 2009; 130:204102. [PMID: 19485432 PMCID: PMC2832053 DOI: 10.1063/1.3139189] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2009] [Accepted: 04/29/2009] [Indexed: 11/15/2022] Open
Abstract
When a system is driven out of equilibrium by a time-dependent protocol that modifies the Hamiltonian, it follows a nonequilibrium path. Samples of these paths can be used in nonequilibrium work theorems to estimate equilibrium quantities such as free energy differences. Here, we consider analyzing paths generated with one protocol using another one. It is posited that analysis protocols which minimize the lag, the difference between the nonequilibrium and the instantaneous equilibrium densities, will reduce the dissipation of reprocessed trajectories and lead to better free energy estimates. Indeed, when minimal lag analysis protocols based on exactly soluble propagators or relative entropies are applied to several test cases, substantial gains in the accuracy and precision of estimated free energy differences are observed.
Collapse
Affiliation(s)
- David D L Minh
- Laboratory of Chemical Physics, NIDDK, National Institutes of Health, Bethesda, Maryland 20892-0520, USA.
| |
Collapse
|
15
|
Adaptive Path-Integral Monte Carlo Methods for Accurate Computation of Molecular Thermodynamic Properties. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470141649.ch5] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
16
|
Calculation of the Free Energy and the Entropy of Macromolecular Systems by Computer Simulation. REVIEWS IN COMPUTATIONAL CHEMISTRY 2007. [DOI: 10.1002/9780470125892.ch1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
17
|
de Koning M. Optimizing the driving function for nonequilibrium free-energy calculations in the linear regime: a variational approach. J Chem Phys 2006; 122:104106. [PMID: 15836308 DOI: 10.1063/1.1860556] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We consider the issue of optimizing linear-regime nonequilibrium simulations to estimate free-energy differences. In particular, we focus on the problem of finding the best-possible driving function lambda(t) that, for a given thermodynamic path, simulation algorithm, and amount of computational effort, minimizes dissipation. From the fluctuation-dissipation theorem it follows that, in the linear-response regime, the dissipation is controlled by the magnitude and characteristic correlation time of the equilibrium fluctuations in the driving force. As a result, the problem of finding the optimal switching scheme involves the solution of a standard problem in variational calculus: the minimization of a functional with respect to the switching function. In practice, the minimization involves solving the associated Euler-Lagrange equation subject to a set of boundary conditions. As a demonstration we apply the approach to the simple, yet illustrative problem of computing the free-energy difference between two classical harmonic oscillators with very different characteristic frequencies.
Collapse
Affiliation(s)
- Maurice de Koning
- Instituto de Física, Universidade de São Paulo, Caixa Postal 66318, 05315-970 São Paulo, São Paulo, Brazil
| |
Collapse
|
18
|
|
19
|
Djikaev Y. Histogram analysis as a method for determining the line tension of a three-phase contact region by Monte Carlo simulations. J Chem Phys 2005; 123:184704. [PMID: 16292918 DOI: 10.1063/1.2056548] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method is proposed for determining the line tension, which is the main physical characteristic of a three-phase contact region, by Monte Carlo (MC) simulations. The key idea of the proposed method is that if a three-phase equilibrium involves a three-phase contact region, the probability distribution of states of a system as a function of two order parameters depends not only on the surface tension, but also on the line tension. This probability distribution can be obtained as a normalized histogram by appropriate MC simulations, so one can use the combination of histogram analysis and finite-size scaling to study the properties of a three phase contact region. Every histogram and results extracted therefrom will depend on the size of the simulated system. Carrying out MC simulations for a series of system sizes and extrapolating the results, obtained from the corresponding series of histograms, to infinite size, one can determine the line tension of the three phase contact region and the interfacial tensions of all three interfaces (and hence the contact angles) in an infinite system. To illustrate the proposed method, it is applied to the three-dimensional ternary fluid mixture, in which molecular pairs of like species do not interact whereas those of unlike species interact as hard spheres. The simulated results are in agreement with expectations.
Collapse
Affiliation(s)
- Yuri Djikaev
- Department of Chemistry, Baker Laboratory, Cornell University, Ithaca, New York 14853-1301, USA.
| |
Collapse
|
20
|
Gloor GJ, Jackson G, Blas FJ, de Miguel E. Test-area simulation method for the direct determination of the interfacial tension of systems with continuous or discontinuous potentials. J Chem Phys 2005; 123:134703. [PMID: 16223322 DOI: 10.1063/1.2038827] [Citation(s) in RCA: 244] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A novel test-area (TA) technique for the direct simulation of the interfacial tension of systems interacting through arbitrary intermolecular potentials is presented in this paper. The most commonly used method invokes the mechanical relation for the interfacial tension in terms of the tangential and normal components of the pressure tensor relative to the interface (the relation of Kirkwood and Buff [J. Chem. Phys. 17, 338 (1949)]). For particles interacting through discontinuous intermolecular potentials (e.g., hard-core fluids) this involves the determination of delta functions which are impractical to evaluate, particularly in the case of nonspherical molecules. By contrast we employ a thermodynamic route to determine the surface tension from a free-energy perturbation due to a test change in the surface area. There are important distinctions between our test-area approach and the computation of a free-energy difference of two (or more) systems with different interfacial areas (the method of Bennett [J. Comput. Phys. 22, 245 (1976)]), which can also be used to determine the surface tension. In order to demonstrate the adequacy of the method, the surface tension computed from test-area Monte Carlo (TAMC) simulations are compared with the data obtained with other techniques (e.g., mechanical and free-energy differences) for the vapor-liquid interface of Lennard-Jones and square-well fluids; the latter corresponds to a discontinuous potential which is difficult to treat with standard methods. Our thermodynamic test-area approach offers advantages over existing techniques of computational efficiency, ease of implementation, and generality. The TA method can easily be implemented within either Monte Carlo (TAMC) or molecular-dynamics (TAMD) algorithms for different types of interfaces (vapor-liquid, liquid-liquid, fluid-solid, etc.) of pure systems and mixtures consisting of complex polyatomic molecules.
Collapse
Affiliation(s)
- Guy J Gloor
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | | | | | | |
Collapse
|
21
|
Hoy RS, Robbins MO. Fcc-bcc transition for Yukawa interactions determined by applied strain deformation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 69:056103. [PMID: 15244879 DOI: 10.1103/physreve.69.056103] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/31/2003] [Indexed: 05/24/2023]
Abstract
Calculations of the work required to transform between bcc and fcc phases yield a high-precision bcc-fcc transition line for monodisperse point Yukawa (screened-Coulomb) systems. Our results agree qualitatively but not quantitatively with recently published simulations and phenomenological criteria for the bcc-fcc transition. In particular, the bcc-fcc-fluid triple point lies at a higher inverse screening length than previously reported.
Collapse
Affiliation(s)
- Robert S Hoy
- Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | | |
Collapse
|
22
|
Ytreberg FM, Zuckerman DM. Efficient use of nonequilibrium measurement to estimate free energy differences for molecular systems. J Comput Chem 2004; 25:1749-59. [PMID: 15362132 DOI: 10.1002/jcc.20103] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
A promising method for calculating free energy differences DeltaF is to generate nonequilibrium data via "fast-growth" simulations or by experiments--and then use Jarzynski's equality. However, a difficulty with using Jarzynski's equality is that DeltaF estimates converge very slowly and unreliably due to the nonlinear nature of the calculation--thus requiring large, costly data sets. The purpose of the work presented here is to determine the best estimate for DeltaF given a (finite) set of work values previously generated by simulation or experiment. Exploiting statistical properties of Jarzynski's equality, we present two fully automated analyses of nonequilibrium data from a toy model, and various simulated molecular systems. Both schemes remove at least several k(B)T of bias from DeltaF estimates, compared to direct application of Jarzynski's equality, for modest sized data sets (100 work values), in all tested systems. Results from one of the new methods suggest that good estimates of DeltaF can be obtained using 5-40-fold less data than was previously possible. Extending previous work, the new results exploit the systematic behavior of bias due to finite sample size. A key innovation is better use of the more statistically reliable information available from the raw data.
Collapse
Affiliation(s)
- F Marty Ytreberg
- Center for Computational Biology and Bioinformatics, University of Pittsburgh, 200 Lothrop St., Pittsburgh, Pennsylvania 15261, USA.
| | | |
Collapse
|
23
|
Lu N, Adhikari J, Kofke DA. Variational formula for the free energy based on incomplete sampling in a molecular simulation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 68:026122. [PMID: 14525064 DOI: 10.1103/physreve.68.026122] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2003] [Indexed: 05/24/2023]
Abstract
Finite sampling in free-energy perturbation (FEP) calculations by molecular simulation leads to reproducible systematic errors, with sign shown to depend (in a known way) only on which system governs sampling in the simulation. Thus the result of a FEP calculation can be used as a bound on the true free energy. This inequality is of a wholly different nature from established forms such as the Gibbs-Bogoliubov inequality or the second law, in that its origins relate to the performance of a molecular simulation. If one can identify a suitable reference system having a free energy known as a function of some defining parameter, variational schemes based on the finite-sampling inequalities can be implemented. This idea is demonstrated through calculation of the free energy of a hard-sphere solid by perturbing from harmonic references and of a hard-sphere fluid by perturbing from infinitely polydisperse references. The tightness of the bounds increases with the amount of sampling in the simulation and correlates with the entropy difference between the target and reference systems. The bounds are tightest near the point where the entropy difference is least.
Collapse
Affiliation(s)
- Nandou Lu
- Department of Chemical Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, USA
| | | | | |
Collapse
|
24
|
|
25
|
Zuckerman DM, Woolf TB. Overcoming finite-sampling errors in fast-switching free-energy estimates: extrapolative analysis of a molecular system. Chem Phys Lett 2002. [DOI: 10.1016/s0009-2614(01)01397-5] [Citation(s) in RCA: 56] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
26
|
de Koning M, Antonelli A, Yip S. Single-simulation determination of phase boundaries: A dynamic Clausius–Clapeyron integration method. J Chem Phys 2001. [DOI: 10.1063/1.1420486] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
|
28
|
Hendrix DA, Jarzynski C. A “fast growth” method of computing free energy differences. J Chem Phys 2001. [DOI: 10.1063/1.1353552] [Citation(s) in RCA: 136] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
29
|
|
30
|
Miller MA, Reinhardt WP. Efficient free energy calculations by variationally optimized metric scaling: Concepts and applications to the volume dependence of cluster free energies and to solid–solid phase transitions. J Chem Phys 2000. [DOI: 10.1063/1.1313537] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
Amon LM, Reinhardt WP. Development of reference states for use in absolute free energy calculations of atomic clusters with application to 55-atom Lennard-Jones clusters in the solid and liquid states. J Chem Phys 2000. [DOI: 10.1063/1.1286808] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
de Koning M, Wei Cai, Antonelli A, Yip S. Efficient free-energy calculations by the simulation of nonequilibrium processes. Comput Sci Eng 2000. [DOI: 10.1109/5992.841802] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
33
|
Cummins PL, Gready JE. Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules. J Comput Chem 1999. [DOI: 10.1002/(sici)1096-987x(19990730)20:10<1028::aid-jcc5>3.0.co;2-7] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
34
|
Schön JC. A thermodynamic distance criterion of optimality for the calculation of free energy changes from computer simulations. J Chem Phys 1996. [DOI: 10.1063/1.472836] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Ebeling M, Nadler W. Reconstruction of a first-order phase transition from computer simulations of individual phases and subphases. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1996; 53:3365-3368. [PMID: 9964644 DOI: 10.1103/physreve.53.3365] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
36
|
Hunter JE, Reinhardt WP. Finite‐size scaling behavior of the free energy barrier between coexisting phases: Determination of the critical temperature and interfacial tension of the Lennard‐Jones fluid. J Chem Phys 1995. [DOI: 10.1063/1.470121] [Citation(s) in RCA: 50] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Hogenson GJ, Reinhardt WP. Variational upper and lower bounds on quantum free energy and energy differences via path integral Monte Carlo. J Chem Phys 1995. [DOI: 10.1063/1.468543] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
Tsao L, Sheu S, Mou C. Absolute entropy of simple point charge model water by adiabatic switching processes. J Chem Phys 1994. [DOI: 10.1063/1.467670] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|