Optimal estimators and asymptotic variances for nonequilibrium path-ensemble averages.

David D L Minh, John D Chodera
Author Information
  1. David D L Minh: Laboratory of Chemical Physics, NIDDK, National Institutes of Health, Bethesda, Maryland 20892, USA. daveminh@gmail.com

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.

References

  1. Biophys J. 2007 Nov 15;93(10):3373-81 [PMID: 17704183]
  2. J Chem Phys. 2008 Jul 14;129(2):024102 [PMID: 18624511]
  3. J Phys Chem B. 2005 Apr 14;109(14):6902-15 [PMID: 16851777]
  4. J Chem Phys. 2005 Apr 8;122(14):144107 [PMID: 15847516]
  5. J Chem Phys. 2008 Jul 14;129(2):024114 [PMID: 18624523]
  6. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1999 Sep;60(3):2721-6 [PMID: 11970075]
  7. J Chem Phys. 2007 Jun 28;126(24):244101 [PMID: 17614531]
  8. J Chem Phys. 2005 Nov 22;123(20):204909 [PMID: 16351319]
  9. J Chem Theory Comput. 2007 Jan;3(1):26-41 [PMID: 26627148]
  10. J Phys Chem B. 2008 May 15;112(19):5892-7 [PMID: 18088108]
  11. J Comput Chem. 2009 Aug;30(11):1726-36 [PMID: 19408280]
  12. Phys Rev Lett. 2007 Mar 9;98(10):108301 [PMID: 17358574]
  13. Phys Rev E Stat Nonlin Soft Matter Phys. 2006 Dec;74(6 Pt 1):061120 [PMID: 17280051]
  14. Phys Rev E Stat Nonlin Soft Matter Phys. 2009 Jan;79(1 Pt 1):012104 [PMID: 19257090]
  15. Proc Natl Acad Sci U S A. 2003 Oct 28;100(22):12564-9 [PMID: 14528008]
  16. J Phys Chem B. 2007 Apr 26;111(16):4137-40 [PMID: 17407342]
  17. Phys Rev E Stat Nonlin Soft Matter Phys. 2006 Apr;73(4 Pt 2):046105 [PMID: 16711876]
  18. J Chem Phys. 2004 Dec 1;121(21):10392-400 [PMID: 15549919]
  19. Phys Rev Lett. 2003 Oct 3;91(14):140601 [PMID: 14611511]
  20. J Chem Phys. 2009 Apr 14;130(14):144908 [PMID: 19368472]
  21. Phys Rev Lett. 2009 Feb 20;102(7):070602 [PMID: 19257655]
  22. J Chem Phys. 2009 May 28;130(20):204102 [PMID: 19485432]
  23. Nature. 2005 Sep 8;437(7056):231-4 [PMID: 16148928]
  24. Science. 2002 Jun 7;296(5574):1832-5 [PMID: 12052949]
  25. Phys Rev Lett. 2002 Oct 28;89(18):180602 [PMID: 12398588]
  26. Phys Rev Lett. 2006 Mar 17;96(10):100602 [PMID: 16605720]
  27. Acc Chem Res. 2005 Jul;38(7):504-13 [PMID: 16028884]
  28. Phys Rev Lett. 1989 Sep 18;63(12):1195-1198 [PMID: 10040500]
  29. J Chem Phys. 2004 Jun 15;120(23):10876-9 [PMID: 15268117]
  30. J Chem Phys. 2007 Aug 14;127(6):065101 [PMID: 17705625]
  31. J Chem Phys. 2008 Jan 14;128(2):024104 [PMID: 18205440]
  32. J Chem Phys. 2008 Sep 28;129(12):124105 [PMID: 19045004]
  33. Proc Natl Acad Sci U S A. 2001 Mar 27;98(7):3658-61 [PMID: 11274384]
  34. Phys Rev Lett. 2008 May 9;100(18):180602 [PMID: 18518359]

Grants

  1. /Intramural NIH HHS

MeSH Term

Models, Chemical
Reproducibility of Results
Thermodynamics

Word Cloud

Created with Highcharts 10.0.0estimatorsnonequilibriumasymptoticframeworkvarianceestimatespath-ensembleaveragesestimatorbidirectionalmodelconfidenceintervalsbiasExistingoptimalshownfallwithinextendedbridgesamplingUsingderivegeneralminimal-variancecancombinetrajectorydatasampledmultiplepath-ensemblesestimatearbitraryfunctionsexpectationsalsoappliedobtainusefulmeasurestatisticaluncertaintyparticulardeveloppertainingJarzynski'sequalityfreeenergiesHummer-Szaboexpressionspotentialmeanforcecalculateduni-pathsamplesdemonstratedsingle-moleculepullingexperimentsimulationsexpressionfoundaccuratelycharacterizearoundsmallHenceinaccuratelydescribedunidirectionallargelargelyreflectstrueerrorderivedMinhAdibOptimalvariances

Similar Articles

Cited By