1
|
Hatch HW, Siderius DW, Shen VK. Monte Carlo molecular simulations with FEASST version 0.25.1. J Chem Phys 2024; 161:092501. [PMID: 39234968 DOI: 10.1063/5.0224283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2024] [Accepted: 08/08/2024] [Indexed: 09/06/2024] Open
Abstract
FEASST is an open-source Monte Carlo software package for particle-based simulations. This software, which was released in 2017, has been used to study phase equilibrium, self-assembly, aggregation or gelation in biological materials, colloids, polymers, ionic liquids, and adsorption in porous networks. We highlight some of the unique features available in FEASST, such as flat-histogram grand canonical ensemble, Gibbs ensemble, and Mayer-sampling simulations with support for anisotropic models and parallelization with flat-histogram and prefetching. We also discuss how the challenges of supporting a variety of Monte Carlo algorithms were overcome by an object-oriented design. This also allows others to extend classes, which improves software interoperability, as inspired by LAMMPS classes and user packages. This article describes version 0.25.1 with benchmarks, compilation instructions, and introductory tutorials for running, restarting, and testing simulations, user guidelines, software design strategies, alternative interfaces, and the test-driven development strategy.
Collapse
Affiliation(s)
- Harold W Hatch
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| | - Daniel W Siderius
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| | - Vincent K Shen
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| |
Collapse
|
2
|
Siderius DW, Hatch HW, Shen VK. Flat-Histogram Monte Carlo Simulation of Water Adsorption in Metal-Organic Frameworks. J Phys Chem B 2024; 128:4830-4845. [PMID: 38676704 PMCID: PMC11175621 DOI: 10.1021/acs.jpcb.4c00753] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2024]
Abstract
Molecular simulations of water adsorption in porous materials often converge slowly due to sampling bottlenecks that follow from hydrogen bonding and, in many cases, the formation of water clusters. These effects may be exacerbated in metal-organic framework (MOF) adsorbents, due to the presence of pore spaces (cages) that promote the formation of discrete-size clusters and hydrophobic effects (if present), among other reasons. In Grand Canonical Monte Carlo (MC) simulations, these sampling challenges are typically manifested by low MC acceptance ratios, a tendency for the simulation to become stuck in a particular loading state (i.e., macrostates), and the persistence of specific clusters for long periods of the simulation. We present simulation strategies to address these sampling challenges, by applying flat-histogram MC (FHMC) methods and specialized MC move types to simulations of water adsorption. FHMC, in both Transition-matrix and Wang-Landau forms, drives the simulation to sample relevant macrostates by incorporating weights that are self-consistently adjusted throughout the simulation and generate the macrostate probability distribution (MPD). Specialized MC moves, based on aggregation-volume bias and configurational bias methods, separately address low acceptance ratios for basic MC trial moves and specifically target water molecules in clusters; in turn, the specialized MC moves improve the efficiency of generating new configurations which is ultimately reflected in improved statistics collected by FHMC. The combined strategies are applied to study the adsorption of water in CuBTC and ZIF-8 at 300 K, through examination of the MPD and the adsorption isotherm generated by histogram reweighting. A key result is the appearance of nontrivial oscillations in the MPD, which we show to be associated with water clusters in the adsorption system. Additionally, we show that the probabilities of certain clusters become similar in value near the boundaries of the isotherm hysteresis loop, indicating a strong connection between cluster formation/destruction and the thermodynamic limits of stability. For a hydrophobic MOF, the FHMC results show that the phase transition from low density to high density is suppressed to water pressure far above the bulk-fluid saturation pressure; this is consistent with results presented elsewhere. We also compare our FHMC simulation isotherm to one measured by a different technique but with ostensibly the same molecular interactions and comment on observed differences and the need for follow-up work. The simulation strategies presented here can be applied to the simulation of water in other MOFs using heuristic guidelines laid out in our text, which should facilitate the more consistent and efficient simulation of water adsorption in porous materials in future applications.
Collapse
Affiliation(s)
- Daniel W. Siderius
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, United States
| | - Harold W. Hatch
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, United States
| | - Vincent K. Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, United States
| |
Collapse
|
3
|
Elsässer P, Kuhnhold A. Phase behavior of the generalized chiral Lebwohl-Lasher model in bulk and confinement. Phys Rev E 2022; 105:054704. [PMID: 35706156 DOI: 10.1103/physreve.105.054704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Accepted: 05/06/2022] [Indexed: 06/15/2023]
Abstract
One of the simplest models that is used to study the isotropic-nematic tranition in liquid crystalline systems is the Lebwohl-Lasher model. Several extensions of this model further enhanced its applicability. We combine two of these extensions (a generalization and the inclusion of a chiral term) and study the phase behavior and the nature of phase transitions of the resulting generalized chiral Lebwohl-Lasher model using Monte Carlo simulations. We find that the type (and even the existence) of the transition depends on the combination of the width of the interaction potential, the strength of the chiral part of the interaction, and the geometry of the system. As well, the pitch of the cholesteric bulk phase changes on approaching the phase transition if the interaction width differs from the one of the original Lebwohl-Lasher model. Thus, we identify parameter combinations that allow one to tune the properties of the ordered phase and the nature of the order-disorder transition.
Collapse
Affiliation(s)
- P Elsässer
- Institute of Physics, University of Freiburg, 79104 Freiburg (Breisgau), Germany
| | - A Kuhnhold
- Institute of Physics, University of Freiburg, 79104 Freiburg (Breisgau), Germany
| |
Collapse
|
4
|
Datar A, Witman M, Lin L. Monte Carlo simulations for water adsorption in porous materials: Best practices and new insights. AIChE J 2021. [DOI: 10.1002/aic.17447] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Archit Datar
- William G. Lowrie Department of Chemical and Biomolecular Engineering The Ohio State University Columbus Ohio USA
| | | | - Li‐Chiang Lin
- William G. Lowrie Department of Chemical and Biomolecular Engineering The Ohio State University Columbus Ohio USA
- Department of Chemical Engineering National Taiwan University Taipei Taiwan
| |
Collapse
|
5
|
Mahynski NA, Hatch HW, Witman M, Sheen DA, Errington JR, Shen VK. Flat-histogram extrapolation as a useful tool in the age of big data. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1747617] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Nathan A. Mahynski
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, MD, USA
| | - Harold W. Hatch
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, MD, USA
| | | | - David A. Sheen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, MD, USA
| | - Jeffrey R. Errington
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, NY, USA
| | - Vincent K. Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, MD, USA
| |
Collapse
|
6
|
Witman M, Wright B, Smit B. Simulating Enhanced Methane Deliverable Capacity of Guest Responsive Pores in Intrinsically Flexible MOFs. J Phys Chem Lett 2019; 10:5929-5934. [PMID: 31532681 DOI: 10.1021/acs.jpclett.9b02449] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
A novel computational procedure, based on the principles of flat-histogram Monte Carlo, is developed for facile prediction of the adsorption thermodynamics of intrinsically flexible adsorbents. We then demonstrate how an accurate prediction of methane deliverable capacity in a metal-organic framework (MOF) with significant intrinsic flexibility requires use of such a method. Dynamic side chains in the framework respond to methane adsorbates and reorganize to exhibit a more conducive pore space at high adsorbate densities while simultaneously providing a less conducive pore space at low adsorbate densities. This "responsive pore" MOF achieves ∼20% higher deliverable capacity than if the framework were rigid and elucidates a strategy for designing high deliverable capacity MOFs in the future.
Collapse
Affiliation(s)
- Matthew Witman
- Department of Chemical and Biomolecular Engineering , University of California , Berkeley 94720 , United States
| | - Bradley Wright
- Department of Chemical and Biomolecular Engineering , University of California , Berkeley 94720 , United States
| | - Berend Smit
- Department of Chemical and Biomolecular Engineering , University of California , Berkeley 94720 , United States
- Laboratory of Molecular Simulation (LSMO), Institut des Sciences et Ingénierie Chimiques, Valais , Ecole Polytechnique Fédérale de Lausanne (EPFL) , Rue de l'Industrie 17 , CH-1951 Sion , Switzerland
| |
Collapse
|
7
|
Witman M, Mahynski NA, Smit B. Flat-Histogram Monte Carlo as an Efficient Tool To Evaluate Adsorption Processes Involving Rigid and Deformable Molecules. J Chem Theory Comput 2018; 14:6149-6158. [PMID: 30296088 DOI: 10.1021/acs.jctc.8b00534] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Monte Carlo simulations are the foundational technique for predicting thermodynamic properties of open systems where the process of interest involves the exchange of particles. Thus, they have been used extensively to computationally evaluate the adsorption properties of nanoporous materials and are critical for the in silico identification of promising materials for a variety of gas storage and chemical separation applications. In this work we demonstrate that a well-known biasing technique, known as "flat-histogram" sampling, can be combined with temperature extrapolation of the free energy landscape to efficiently provide significantly more useful thermodynamic information than standard open ensemble MC simulations. Namely, we can accurately compute the isosteric heat of adsorption and number of particles adsorbed for various adsorbates over an extremely wide range of temperatures and pressures from a set of simulations at just one temperature. We extend this derivation of the temperature extrapolation to adsorbates with intramolecular degrees of freedom when Rosenbluth sampling is employed. Consequently, the working capacity and isosteric heat can be computed for any given combined temperature/pressure swing adsorption process for a large range of operating conditions with both rigid and deformable adsorbates. Continuous thermodynamic properties can be computed with this technique at very moderate computational cost, thereby providing a strong case for its application to the in silico identification of promising nanoporous adsorbents.
Collapse
Affiliation(s)
- Matthew Witman
- Department of Chemical and Biomolecular Engineering , University of California , Berkeley 94720 , United States.,Laboratory of Molecular Simulation (LSMO) , Institut des Sciences et Ingénierie Chimiques, Valais, École Polytechnique Fédérale de Lausanne (EPFL) , Rue de l'Industrie 17 , CH-1951 Sion , Switzerland
| | - Nathan A Mahynski
- Chemical Sciences Division , National Institute of Standards and Technology , Gaithersburg , Maryland 20899-8320 , United States
| | - Berend Smit
- Department of Chemical and Biomolecular Engineering , University of California , Berkeley 94720 , United States.,Laboratory of Molecular Simulation (LSMO) , Institut des Sciences et Ingénierie Chimiques, Valais, École Polytechnique Fédérale de Lausanne (EPFL) , Rue de l'Industrie 17 , CH-1951 Sion , Switzerland
| |
Collapse
|
8
|
Mahynski NA, Errington JR, Shen VK. Temperature extrapolation of multicomponent grand canonical free energy landscapes. J Chem Phys 2018; 147:054105. [PMID: 28789543 DOI: 10.1063/1.4996759] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive a method for extrapolating the grand canonical free energy landscape of a multicomponent fluid system from one temperature to another. Previously, we introduced this statistical mechanical framework for the case where kinetic energy contributions to the classical partition function were neglected for simplicity [N. A. Mahynski et al., J. Chem. Phys. 146, 074101 (2017)]. Here, we generalize the derivation to admit these contributions in order to explicitly illustrate the differences that result. Specifically, we show how factoring out kinetic energy effects a priori, in order to consider only the configurational partition function, leads to simpler mathematical expressions that tend to produce more accurate extrapolations than when these effects are included. We demonstrate this by comparing and contrasting these two approaches for the simple cases of an ideal gas and a non-ideal, square-well fluid.
Collapse
Affiliation(s)
- Nathan A Mahynski
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| | - Jeffrey R Errington
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, USA
| | - Vincent K Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| |
Collapse
|
9
|
Mahynski NA, Errington JR, Shen VK. Multivariable extrapolation of grand canonical free energy landscapes. J Chem Phys 2018; 147:234111. [PMID: 29272947 DOI: 10.1063/1.5006906] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive an approach for extrapolating the free energy landscape of multicomponent systems in the grand canonical ensemble, obtained from flat-histogram Monte Carlo simulations, from one set of temperature and chemical potentials to another. This is accomplished by expanding the landscape in a Taylor series at each value of the order parameter which defines its macrostate phase space. The coefficients in each Taylor polynomial are known exactly from fluctuation formulas, which may be computed by measuring the appropriate moments of extensive variables that fluctuate in this ensemble. Here we derive the expressions necessary to define these coefficients up to arbitrary order. In principle, this enables a single flat-histogram simulation to provide complete thermodynamic information over a broad range of temperatures and chemical potentials. Using this, we also show how to combine a small number of simulations, each performed at different conditions, in a thermodynamically consistent fashion to accurately compute properties at arbitrary temperatures and chemical potentials. This method may significantly increase the computational efficiency of biased grand canonical Monte Carlo simulations, especially for multicomponent mixtures. Although approximate, this approach is amenable to high-throughput and data-intensive investigations where it is preferable to have a large quantity of reasonably accurate simulation data, rather than a smaller amount with a higher accuracy.
Collapse
Affiliation(s)
- Nathan A Mahynski
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| | - Jeffrey R Errington
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, USA
| | - Vincent K Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| |
Collapse
|
10
|
Mahynski NA, Jiao S, Hatch HW, Blanco MA, Shen VK. Predicting structural properties of fluids by thermodynamic extrapolation. J Chem Phys 2018; 148:194105. [PMID: 30307179 PMCID: PMC6183068 DOI: 10.1063/1.5026493] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We describe a methodology for extrapolating the structural properties of multicomponent fluids from one thermodynamic state to another. These properties generally include features of a system that may be computed from an individual configuration such as radial distribution functions, cluster size distributions, or a polymer's radius of gyration. This approach is based on the principle of using fluctuations in a system's extensive thermodynamic variables, such as energy, to construct an appropriate Taylor series expansion for these structural properties in terms of intensive conjugate variables, such as temperature. Thus, one may extrapolate these properties from one state to another when the series is truncated to some finite order. We demonstrate this extrapolation for simple and coarse-grained fluids in both the canonical and grand canonical ensembles, in terms of both temperatures and the chemical potentials of different components. The results show that this method is able to reasonably approximate structural properties of such fluids over a broad range of conditions. Consequently, this methodology may be employed to increase the computational efficiency of molecular simulations used to measure the structural properties of certain fluid systems, especially those used in high-throughput or data-driven investigations.
Collapse
Affiliation(s)
- Nathan A. Mahynski
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| | - Sally Jiao
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Harold W. Hatch
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| | - Marco A. Blanco
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
- Institute for Bioscience and Biotechnology Research, University of Maryland, Rockville, Maryland 20850, USA
| | - Vincent K. Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| |
Collapse
|
11
|
Shen VK, Siderius DW, Mahynski NA. Molecular simulation of capillary phase transitions in flexible porous materials. J Chem Phys 2018; 148:124115. [DOI: 10.1063/1.5022171] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Affiliation(s)
- Vincent K. Shen
- Chemical Informatics Research Group, Chemical Sciences Division, Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| | - Daniel W. Siderius
- Chemical Informatics Research Group, Chemical Sciences Division, Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| | - Nathan A. Mahynski
- Chemical Informatics Research Group, Chemical Sciences Division, Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| |
Collapse
|
12
|
Siderius DW, Krekelberg WP, Chiang WS, Shen VK, Liu Y. Quasi-Two-Dimensional Phase Transition of Methane Adsorbed in Cylindrical Silica Mesopores. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:14252-14262. [PMID: 29185779 PMCID: PMC5831196 DOI: 10.1021/acs.langmuir.7b03406] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Using Monte Carlo and molecular dynamics simulations, we examine the adsorption of methane in cylindrical silica mesopores in an effort to understand a possible phase transition of adsorbed methane in MCM-41 and SBA-15 silica that was previously identified by an unexpected increase in the adsorbed fluid density following capillary condensation, as measured by small-angle neutron scattering (SANS) [Chiang, W-S., et al., Langmuir 2016, 32, 8849]. Our initial simulation results identify a roughly 10 % increase in the density of the liquidlike adsorbed phase for either an isotherm with increasing pressure or an isobar with decreasing temperature and that this densification is associated with a local maximum in the isosteric enthalpy of adsorption. Subsequent analysis of the simulated fluid, via computation of bond-orientational order parameters of specific annular layers of the adsorbed fluid, showed that the layers undergo an ordering transition from a disordered, amorphous state to one with two-dimensional hexagonal structure. Furthermore, this two-dimensional restructuring of the fluid occurs at the same thermodynamic state points as the aforementioned densification and local maximum in the isosteric enthalpy of adsorption. We thus conclude that the densification of the fluid is the result of structural reorganization, which is signaled by the maximum in the isosteric enthalpy. Owing to the qualitative similarity of the structural transitions in the simulated and experimental methane fluids, we propose this hexagonal reorganization as a plausible explanation of the densification observed in SANS measurements. Lastly, we speculate how this structural transition may impact the transport properties of the adsorbed fluid.
Collapse
Affiliation(s)
- Daniel W. Siderius
- Chemical Sciences Division, National Institute of Standards and Technology, 100 Bureau Drive M.S. 8320, Gaithersburg, Maryland 20899, United States
| | - William P. Krekelberg
- Chemical Sciences Division, National Institute of Standards and Technology, 100 Bureau Drive M.S. 8320, Gaithersburg, Maryland 20899, United States
| | - Wei-Shan Chiang
- NIST Center for Neutron Research, National Institute of Standards and Technology, 100 Bureau Drive M.S. 6102, Gaithersburg, Maryland 20899, United States
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Colburn Laboratory, Newark, Delaware 19716, United States
| | - Vincent K. Shen
- Chemical Sciences Division, National Institute of Standards and Technology, 100 Bureau Drive M.S. 8320, Gaithersburg, Maryland 20899, United States
| | - Yun Liu
- NIST Center for Neutron Research, National Institute of Standards and Technology, 100 Bureau Drive M.S. 6102, Gaithersburg, Maryland 20899, United States
| |
Collapse
|
13
|
Mahynski NA, Blanco MA, Errington JR, Shen VK. Predicting low-temperature free energy landscapes with flat-histogram Monte Carlo methods. J Chem Phys 2017; 146:074101. [PMID: 28228029 DOI: 10.1063/1.4975331] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
We present a method for predicting the free energy landscape of fluids at low temperatures from flat-histogram grand canonical Monte Carlo simulations performed at higher ones. We illustrate our approach for both pure and multicomponent systems using two different sampling methods as a demonstration. This allows us to predict the thermodynamic behavior of systems which undergo both first order and continuous phase transitions upon cooling using simulations performed only at higher temperatures. After surveying a variety of different systems, we identify a range of temperature differences over which the extrapolation of high temperature simulations tends to quantitatively predict the thermodynamic properties of fluids at lower ones. Beyond this range, extrapolation still provides a reasonably well-informed estimate of the free energy landscape; this prediction then requires less computational effort to refine with an additional simulation at the desired temperature than reconstruction of the surface without any initial estimate. In either case, this method significantly increases the computational efficiency of these flat-histogram methods when investigating thermodynamic properties of fluids over a wide range of temperatures. For example, we demonstrate how a binary fluid phase diagram may be quantitatively predicted for many temperatures using only information obtained from a single supercritical state.
Collapse
Affiliation(s)
- Nathan A Mahynski
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| | - Marco A Blanco
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| | - Jeffrey R Errington
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, USA
| | - Vincent K Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| |
Collapse
|
14
|
Desgranges C, Delhommelle J. Evaluation of the grand-canonical partition function using expanded Wang-Landau simulations. V. Impact of an electric field on the thermodynamic properties and ideality contours of water. J Chem Phys 2016; 145:184504. [DOI: 10.1063/1.4967336] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, 151 Cornell Street Stop 9024, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, 151 Cornell Street Stop 9024, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
15
|
Lee YH, Yevick D. Renormalized multicanonical sampling in multiple dimensions. Phys Rev E 2016; 94:043323. [PMID: 27841523 DOI: 10.1103/physreve.94.043323] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2016] [Indexed: 11/07/2022]
Abstract
This paper extends to two and three dimensions the recently proposed renormalized multicanonical sampling procedure [D. Yevick, Int. J. Mod. Phys. C 27, 1650033 (2016)IJMPEO0129-183110.1142/S0129183116500339], which generates the density of states of a system from that of a smaller system. The accuracy and efficiency of the method is additionally enhanced through a transition matrix Monte Carlo method that incorporates multicanonical sampling. The performance of the resulting procedure is investigated through calculations of the specific heat of the standard Ising model for square lattices with up to 128^{2} spins and for cubic lattices with up to 24^{3} spins. A comparison of the scaling properties to those of other algorithms based on the transition matrix Monte Carlo method demonstrates that the relative advantage of the renormalized procedure increases rapidly with system size.
Collapse
Affiliation(s)
- Yong Hwan Lee
- Department of Physics, University of Waterloo, Waterloo, Ontario, Canada N2L 3G7
| | - David Yevick
- Department of Physics, University of Waterloo, Waterloo, Ontario, Canada N2L 3G7
| |
Collapse
|
16
|
Desgranges C, Delhommelle J. Evaluation of the grand-canonical partition function using expanded Wang-Landau simulations. IV. Performance of many-body force fields and tight-binding schemes for the fluid phases of silicon. J Chem Phys 2016; 144:124510. [PMID: 27036464 DOI: 10.1063/1.4944619] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We extend Expanded Wang-Landau (EWL) simulations beyond classical systems and develop the EWL method for systems modeled with a tight-binding Hamiltonian. We then apply the method to determine the partition function and thus all thermodynamic properties, including the Gibbs free energy and entropy, of the fluid phases of Si. We compare the results from quantum many-body (QMB) tight binding models, which explicitly calculate the overlap between the atomic orbitals of neighboring atoms, to those obtained with classical many-body (CMB) force fields, which allow to recover the tetrahedral organization in condensed phases of Si through, e.g., a repulsive 3-body term that favors the ideal tetrahedral angle. Along the vapor-liquid coexistence, between 3000 K and 6000 K, the densities for the two coexisting phases are found to vary significantly (by 5 orders of magnitude for the vapor and by up to 25% for the liquid) and to provide a stringent test of the models. Transitions from vapor to liquid are predicted to occur for chemical potentials that are 10%-15% higher for CMB models than for QMB models, and a ranking of the force fields is provided by comparing the predictions for the vapor pressure to the experimental data. QMB models also reveal the formation of a gap in the electronic density of states of the coexisting liquid at high temperatures. Subjecting Si to a nanoscopic confinement has a dramatic effect on the phase diagram with, e.g. at 6000 K, a decrease in liquid densities by about 50% for both CMB and QMB models and an increase in vapor densities between 90% (CMB) and 170% (QMB). The results presented here provide a full picture of the impact of the strategy (CMB or QMB) chosen to model many-body effects on the thermodynamic properties of the fluid phases of Si.
Collapse
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, 151 Cornell Street Stop 9024, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, 151 Cornell Street Stop 9024, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
17
|
Haber R, Hoffmann KH. Extending the parQ transition matrix method to grand canonical ensembles. Phys Rev E 2016; 93:063314. [PMID: 27415394 DOI: 10.1103/physreve.93.063314] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2015] [Indexed: 11/07/2022]
Abstract
Phase coexistence properties as well as other thermodynamic features of fluids can be effectively determined from the grand canonical density of states (DOS). We present an extension of the parQ transition matrix method in combination with the efasTM method as a very fast approach for determining the grand canonical DOS from the transition matrix. The efasTM method minimizes the deviation from detailed balance in the transition matrix using a fast Krylov-based equation solver. The method allows a very effective use of state space transition data obtained by different exploration schemes. An application to a Lennard-Jones system produces phase coexistence properties of the same quality as reference data.
Collapse
Affiliation(s)
- René Haber
- Institut für Physik, Technische Universität Chemnitz, 09107 Chemnitz, Germany
| | - Karl Heinz Hoffmann
- Institut für Physik, Technische Universität Chemnitz, 09107 Chemnitz, Germany
| |
Collapse
|
18
|
Hatch HW, Yang SY, Mittal J, Shen VK. Self-assembly of trimer colloids: effect of shape and interaction range. SOFT MATTER 2016; 12:4170-4179. [PMID: 27087490 PMCID: PMC4939708 DOI: 10.1039/c6sm00473c] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
Trimers with one attractive bead and two repulsive beads, similar to recently synthesized trimer patchy colloids, were simulated with flat-histogram Monte Carlo methods to obtain the stable self-assembled structures for different shapes and interaction potentials. Extended corresponding states principle was successfully applied to self-assembling systems in order to approximately collapse the results for models with the same shape, but different interaction range. This helps us directly compare simulation results with previous experiment, and good agreement was found between the two. In addition, a variety of self-assembled structures were observed by varying the trimer geometry, including spherical clusters, elongated clusters, monolayers, and spherical shells. In conclusion, our results help to compare simulations and experiments, via extended corresponding states, and we predict the formation of self-assembled structures for trimer shapes that have not been experimentally synthesized.
Collapse
Affiliation(s)
- Harold W. Hatch
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA;
| | - Seung-Yeob Yang
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA;
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, USA;
| | - Vincent K. Shen
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA;
| |
Collapse
|
19
|
Gor GY, Siderius DW, Rasmussen CJ, Krekelberg WP, Shen VK, Bernstein N. Relation between pore size and the compressibility of a confined fluid. J Chem Phys 2016; 143:194506. [PMID: 26590541 DOI: 10.1063/1.4935430] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
When a fluid is confined to a nanopore, its thermodynamic properties differ from the properties of a bulk fluid, so measuring such properties of the confined fluid can provide information about the pore sizes. Here, we report a simple relation between the pore size and isothermal compressibility of argon confined in such pores. Compressibility is calculated from the fluctuations of the number of particles in the grand canonical ensemble using two different simulation techniques: conventional grand-canonical Monte Carlo and grand-canonical ensemble transition-matrix Monte Carlo. Our results provide a theoretical framework for extracting the information on the pore sizes of fluid-saturated samples by measuring the compressibility from ultrasonic experiments.
Collapse
Affiliation(s)
- Gennady Y Gor
- NRC Research Associate, Resident at Center for Computational Materials Science, Naval Research Laboratory, Washington, DC 20375, USA
| | - Daniel W Siderius
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Christopher J Rasmussen
- DuPont Central Research and Development Experimental Station, Wilmington, Delaware 19803, USA
| | - William P Krekelberg
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Vincent K Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Noam Bernstein
- Center for Computational Materials Science, Naval Research Laboratory, Washington, DC 20375, USA
| |
Collapse
|
20
|
Hatch HW, Mittal J, Shen VK. Computational study of trimer self-assembly and fluid phase behavior. J Chem Phys 2016; 142:164901. [PMID: 25933785 DOI: 10.1063/1.4918557] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
The fluid phase diagram of trimer particles composed of one central attractive bead and two repulsive beads was determined as a function of simple geometric parameters using flat-histogram Monte Carlo methods. A variety of self-assembled structures were obtained including spherical micelle-like clusters, elongated clusters, and densely packed cylinders, depending on both the state conditions and shape of the trimer. Advanced simulation techniques were employed to determine transitions between self-assembled structures and macroscopic phases using thermodynamic and structural definitions. Simple changes in particle geometry yield dramatic changes in phase behavior, ranging from macroscopic fluid phase separation to molecular-scale self-assembly. In special cases, both self-assembled, elongated clusters and bulk fluid phase separation occur simultaneously. Our work suggests that tuning particle shape and interactions can yield superstructures with controlled architecture.
Collapse
Affiliation(s)
- Harold W Hatch
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, USA
| | - Vincent K Shen
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| |
Collapse
|
21
|
Whitmer JK, Fluitt AM, Antony L, Qin J, McGovern M, de Pablo JJ. Sculpting bespoke mountains: Determining free energies with basis expansions. J Chem Phys 2015; 143:044101. [DOI: 10.1063/1.4927147] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Jonathan K. Whitmer
- Department of Chemical and Biomolecular Engineering, University of Notre Dame du Lac, Notre Dame, Indiana 46556, USA
| | - Aaron M. Fluitt
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - Lucas Antony
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - Jian Qin
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
- Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Michael McGovern
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - Juan J. de Pablo
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
- Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
22
|
Hemmen A, Panagiotopoulos AZ, Gross J. Grand Canonical Monte Carlo Simulations Guided by an Analytic Equation of State—Transferable Anisotropic Mie Potentials for Ethers. J Phys Chem B 2015; 119:7087-99. [DOI: 10.1021/acs.jpcb.5b01806] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Andrea Hemmen
- Institute
of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | | | - Joachim Gross
- Institute
of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
23
|
Vink RLC. Crossover from a Kosterlitz-Thouless phase transition to a discontinuous phase transition in two-dimensional liquid crystals. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:062132. [PMID: 25615069 DOI: 10.1103/physreve.90.062132] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2014] [Indexed: 06/04/2023]
Abstract
Liquid crystals in two dimensions do not support long-range nematic order, but a quasinematic phase where the orientational correlations decay algebraically is possible. The transition from the isotropic to the quasinematic phase can be continuous and of the Kosterlitz-Thouless type, or it can be first order. We report here on a liquid-crystal model where the nature of the isotropic to quasinematic transition can be tuned via a single parameter p in the pair potential. For p<p(t), the transition is of the Kosterlitz-Thouless type, while for p>p(t), it is first order. Precisely at p=p(t), there is a tricritical point where, in addition to the orientational correlations, also the positional correlations decay algebraically. The tricritical behavior is analyzed in detail, including an accurate estimate of p(t). The results follow from extensive Monte Carlo simulations combined with a finite-size scaling analysis. Paramount in the analysis is a scheme to facilitate the extrapolation of simulation data in parameters that are not necessarily field variables (in this case, the parameter p), the details of which are also provided. This scheme provides a simple and powerful alternative for situations where standard histogram reweighting cannot be applied.
Collapse
Affiliation(s)
- Richard L C Vink
- Institute of Theoretical Physics, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, D-37077 Göttingen, Germany
| |
Collapse
|
24
|
Shen VK, Siderius DW. Elucidating the effects of adsorbent flexibility on fluid adsorption using simple models and flat-histogram sampling methods. J Chem Phys 2014; 140:244106. [PMID: 24985617 DOI: 10.1063/1.4884124] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Affiliation(s)
- Vincent K Shen
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| | - Daniel W Siderius
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| |
Collapse
|
25
|
Desgranges C, Delhommelle J. Evaluation of the grand-canonical partition function using expanded Wang-Landau simulations. III. Impact of combining rules on mixtures properties. J Chem Phys 2014; 140:104109. [DOI: 10.1063/1.4867498] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|
26
|
Vink RLC. A finite-temperature Monte Carlo algorithm for network forming materials. J Chem Phys 2014; 140:104509. [DOI: 10.1063/1.4867897] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Escobedo FA. Mapping coexistence lines via free-energy extrapolation: Application to order-disorder phase transitions of hard-core mixtures. J Chem Phys 2014; 140:094102. [DOI: 10.1063/1.4866764] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
|
28
|
Hicks E, Desgranges C, Delhommelle J. Adsorption and diffusion of the antiparkinsonian drug amantadine in carbon nanotubes. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2013.841908] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
29
|
Koenig A, Desgranges C, Delhommelle J. Adsorption of hydrogen in covalent organic frameworks using expanded Wang–Landau simulations. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.841907] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- A.R.V. Koenig
- Department of Chemistry, University of North Dakota, 151 Cornell Street Stop 9024, Grand Forks, ND58202, USA
| | - C. Desgranges
- Department of Chemistry, University of North Dakota, 151 Cornell Street Stop 9024, Grand Forks, ND58202, USA
| | - J. Delhommelle
- Department of Chemistry, University of North Dakota, 151 Cornell Street Stop 9024, Grand Forks, ND58202, USA
| |
Collapse
|
30
|
Chaimovich A, Shell MS. Length-scale crossover of the hydrophobic interaction in a coarse-grained water model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:052313. [PMID: 24329270 DOI: 10.1103/physreve.88.052313] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/28/2013] [Indexed: 06/03/2023]
Abstract
It has been difficult to establish a clear connection between the hydrophobic interaction among small molecules typically studied in molecular simulations (a weak, oscillatory force) and that found between large, macroscopic surfaces in experiments (a strong, monotonic force). Here, we show that both types of interaction can emerge with a simple, core-softened water model that captures water's unique pairwise structure. As in hydrophobic hydration, we find that the hydrophobic interaction manifests a length-scale dependence, exhibiting distinct driving forces in the molecular and macroscopic regimes. Moreover, the ability of this simple model to capture both regimes suggests that several features of the hydrophobic force can be understood merely through water's pair correlations.
Collapse
Affiliation(s)
- Aviel Chaimovich
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara, California 93106-5080, USA
| | - M Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara, California 93106-5080, USA
| |
Collapse
|
31
|
Rane KS, Murali S, Errington JR. Monte Carlo Simulation Methods for Computing Liquid–Vapor Saturation Properties of Model Systems. J Chem Theory Comput 2013; 9:2552-66. [DOI: 10.1021/ct400074p] [Citation(s) in RCA: 71] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Kaustubh S. Rane
- Department of Chemical
and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200,
United States
| | - Sabharish Murali
- Department of Chemical
and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200,
United States
| | - Jeffrey R. Errington
- Department of Chemical
and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200,
United States
| |
Collapse
|
32
|
Aleksandrov T, Desgranges C, Delhommelle J. Numerical estimate for boiling points via Wang–Landau simulations. MOLECULAR SIMULATION 2012. [DOI: 10.1080/08927022.2012.702906] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
33
|
Do H, Wheatley RJ. Density of States Partitioning Method for Calculating the Free Energy of Solids. J Chem Theory Comput 2012; 9:165-71. [DOI: 10.1021/ct3007056] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Hainam Do
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K
| | - Richard J. Wheatley
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K
| |
Collapse
|
34
|
Maerzke KA, Gai L, Cummings PT, McCabe C. Incorporating configurational-bias Monte Carlo into the Wang-Landau algorithm for continuous molecular systems. J Chem Phys 2012. [DOI: 10.1063/1.4766354] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|
35
|
Ngale KN, Desgranges C, Delhommelle J. Wang–Landau configurational bias Monte Carlo simulations: vapour–liquid equilibria of alkenes. MOLECULAR SIMULATION 2012. [DOI: 10.1080/08927022.2012.694432] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
36
|
Desgranges C, Delhommelle J. Evaluation of the grand-canonical partition function using expanded Wang-Landau simulations. I. Thermodynamic properties in the bulk and at the liquid-vapor phase boundary. J Chem Phys 2012; 136:184107. [DOI: 10.1063/1.4712023] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
|
37
|
Desgranges C, Delhommelle J. Evaluation of the grand-canonical partition function using expanded Wang-Landau simulations. II. Adsorption of atomic and molecular fluids in a porous material. J Chem Phys 2012; 136:184108. [DOI: 10.1063/1.4712025] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
|
38
|
Do H, Hirst JD, Wheatley RJ. Calculation of Partition Functions and Free Energies of a Binary Mixture Using the Energy Partitioning Method: Application to Carbon Dioxide and Methane. J Phys Chem B 2012; 116:4535-42. [DOI: 10.1021/jp212168f] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Hainam Do
- School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, United Kingdom
| | - Jonathan D. Hirst
- School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, United Kingdom
| | - Richard J. Wheatley
- School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, United Kingdom
| |
Collapse
|
39
|
Abstract
One of the central problems in statistical mechanics is that of finding the density of states of a system. Knowledge of the density of states of a system is equivalent to knowledge of its fundamental equation, from which all thermodynamic quantities can be obtained. Over the past several years molecular simulations have made considerable strides in their ability to determine the density of states of complex fluids and materials. In this review we discuss some of the more promising approaches proposed in the recent literature along with their advantages and limitations.
Collapse
Affiliation(s)
- Sadanand Singh
- Department of Chemical and Biological Engineering, University of Wisconsin, Madison, WI 53706, USA
| | | | | |
Collapse
|
40
|
Zhang P, Li B, Wang Q. Quantitative Study of Fluctuation Effects by Fast Lattice Monte Carlo Simulations. III. Homopolymer Brushes in an Explicit Solvent. Macromolecules 2012. [DOI: 10.1021/ma202454s] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pengfei Zhang
- School of Physics, Nankai University, Tianjin, P. R. China 300071
- Department of Chemical and Biological Engineering, Colorado State University, Fort Collins, Colorado 80523-1370,
United States
| | - Baohui Li
- School of Physics, Nankai University, Tianjin, P. R. China 300071
| | - Qiang Wang
- Department of Chemical and Biological Engineering, Colorado State University, Fort Collins, Colorado 80523-1370,
United States
| |
Collapse
|
41
|
Kumar V, Sridhar S, Errington JR. Monte Carlo simulation strategies for computing the wetting properties of fluids at geometrically rough surfaces. J Chem Phys 2011; 135:184702. [DOI: 10.1063/1.3655817] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
42
|
Do H, Hirst JD, Wheatley RJ. Rapid calculation of partition functions and free energies of fluids. J Chem Phys 2011; 135:174105. [DOI: 10.1063/1.3656296] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
Zhang P, Li B, Wang Q. Quantitative Study of Fluctuation Effects by Fast Lattice Monte Carlo Simulations. 2. Homopolymer Brushes in an Implicit, Good Solvent. Macromolecules 2011. [DOI: 10.1021/ma201310e] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pengfei Zhang
- School of Physics, Nankai University, Tianjin, P. R. China 300071
| | - Baohui Li
- School of Physics, Nankai University, Tianjin, P. R. China 300071
| | - Qiang Wang
- Department of Chemical and Biological Engineering, Colorado State University, Fort Collins, Colorado 80523-1370, United States
| |
Collapse
|
44
|
Paluch AS, Shah JK, Maginn EJ. Efficient Solvation Free Energy Calculations of Amino Acid Analogs by Expanded Ensemble Molecular Simulation. J Chem Theory Comput 2011; 7:1394-403. [DOI: 10.1021/ct1006746] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Andrew S. Paluch
- Department of Chemical and Biomolecular Engineering and ‡Center for Research Computing, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Jindal K. Shah
- Department of Chemical and Biomolecular Engineering and ‡Center for Research Computing, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Edward J. Maginn
- Department of Chemical and Biomolecular Engineering and ‡Center for Research Computing, University of Notre Dame, Notre Dame, Indiana 46556, United States
| |
Collapse
|
45
|
Wu D. Understanding free-energy perturbation calculations through a model of harmonic oscillators: theory and implications to improve the sampling efficiency by molecular simulation. J Chem Phys 2011; 133:244116. [PMID: 21197985 DOI: 10.1063/1.3511703] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Free-energy perturbation calculation is frequently used to calculate free-energy differences because it is easy to implement and the computation is fast. However, the calculation is subject to large inaccuracies in some circumstances due to the insufficient sampling of the relevant tails of the energy-difference distributions. Here we expand this knowledge of insufficient sampling into a two-dimensional (2D) energy space using a model of harmonic oscillators. We show analytically the relation between the energies of the sampling system and those of the desired target energy spaces, which provide the basis to understand the difficulties in free-energy perturbation calculations. We clarify the reasons of the inaccurate calculation in the different harmonic cases that stem from the spatial separations of the reference and the target energy pairs located in the two-dimensional energy space. The potential-energy space introduced into this 2D energy-space model provides additional clues to improve the sampling efficiency. Based on this understanding, we propose two ways to calculate the free-energy differences using the two schemes of the distribution method. We show that the distribution method implemented in the appropriate energy space--the energy-difference space and the potential-energy space, respectively--can improve the calculation of free energies in different circumstances. This analysis implies that the sampling can be improved if it is directed toward the appropriate region in the potential-energy space, which is easily implemented in various types of free-energy calculations. To test this, we calculate the free-energy surface of alanine dipeptide in gas phase and in aqueous phase, respectively. We demonstrate that the free-energy surface calculation is improved when the biased sampling of the potential energy is integrated into the sampling scheme.
Collapse
Affiliation(s)
- Di Wu
- Department of Physiology and Biophysics, School of Life Sciences, Fudan University, Shanghai 200433, China.
| |
Collapse
|
46
|
Schmidle H, Klapp SHL. Phase transitions of two-dimensional dipolar fluids in external fields. J Chem Phys 2011; 134:114903. [DOI: 10.1063/1.3564916] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
|
47
|
Paluch AS, Jayaraman S, Shah JK, Maginn EJ. A method for computing the solubility limit of solids: application to sodium chloride in water and alcohols. J Chem Phys 2011; 133:124504. [PMID: 20886947 DOI: 10.1063/1.3478539] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
We present an adaptable method to compute the solubility limit of solids by molecular simulation, which avoids the difficulty of reference state calculations. In this way, the method is highly adaptable to molecules of complex topology. Results are shown for solubility calculations of sodium chloride in water and light alcohols at atmospheric conditions. The pseudosupercritical path integration method is used to calculate the free energy of the solid and gives results that are in good agreement with previous studies that reference the Einstein crystal. For the solution phase calculations, the self-adaptive Wang-Landau transition-matrix Monte Carlo method is used within the context of an expanded isothermal-isobaric ensemble. The method shows rapid convergence properties and the uncertainty in the calculated chemical potential was 1% or less for all cases. The present study underpredicts the solubility limit of sodium chloride in water, suggesting a shortcoming of the molecular models. Importantly, the proper trend for the chemical potential in various solvents was captured, suggesting that relative solubilities can be computed by the method.
Collapse
Affiliation(s)
- Andrew S Paluch
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | | | | | | |
Collapse
|
48
|
Cao QL, Wang WL, Li YD, Liu CS. Correlations among residual multiparticle entropy, local atomic-level pressure, free volume and the phase-ordering rule in several liquids. J Chem Phys 2011; 134:044508. [DOI: 10.1063/1.3524206] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
49
|
Shimoyama H, Nakamura H, Yonezawa Y. Simple and effective application of the Wang–Landau method for multicanonical molecular dynamics simulation. J Chem Phys 2011; 134:024109. [DOI: 10.1063/1.3517105] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hiromitsu Shimoyama
- Laboratory of Protein Informatics, Research Center for Structural Biology, Institute for Protein Research, Osaka University 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | | | | |
Collapse
|
50
|
Do H, Wheatley RJ, Hirst JD. Molecular simulation of the binary mixture of 1–1–1–2–tetrafluoroethane and carbon dioxide. Phys Chem Chem Phys 2011; 13:15708-13. [DOI: 10.1039/c1cp21419e] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|