1
|
Vo P, Forsman J, Woodward CE. Boundary-Monte Carlo Method for Neutral and Charged Confined Fluids. J Chem Theory Comput 2022; 18:3766-3780. [PMID: 35575645 DOI: 10.1021/acs.jctc.1c01146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, we describe a new Monte Carlo (MC) simulation method to investigate highly coupled fluids in confined geometries at a constant chemical potential. This method is based on so-called multi-scale Hamiltonian methods, wherein the chemical potential is determined using a more amenable Hamiltonian for a fluid in an "outer" region, which facilitates standard methods, such as grand canonical MC simulations or Widom's particle insertion method. The (inner region) fluid of interest is placed in diffusive contact with the simpler outer fluid via a boundary zone wherein the Hamiltonian is transformed. The current method utilizes an ideal fluid for the outer regions, which allows for implicit rather than explicit simulations. Only the boundary and inner region need explicit consideration; hence, the nomenclature used is boundary-Monte Carlo. We illustrate the utility of the method for simple neutral and charged fluids in cylindrical and planar pores. In the latter case, we use a dense room-temperature ionic liquid model and illustrate how the boundary zone establishes a proper Donnan equilibrium between inner and outer fluids in the presence of charged planar electrodes. Thus, the method allows direct calculation of properties such as the differential capacitance, without the need for additional difficult calculations of the requisite Donnan potential.
Collapse
Affiliation(s)
- Phuong Vo
- School of Science, University of New South Wales, Canberra, Canberra ACT 2600, Australia
| | - Jan Forsman
- Department of Theoretical Chemistry, Chemical Centre, Lund University, Lund S-22100, Sweden
| | - Clifford E Woodward
- School of Science, University of New South Wales, Canberra, Canberra ACT 2600, Australia
| |
Collapse
|
2
|
Vo P, Lu H, Ma K, Forsman J, Woodward CE. Local Grand Canonical Monte Carlo Simulation Method for Confined Fluids. J Chem Theory Comput 2019; 15:6944-6957. [PMID: 31665596 DOI: 10.1021/acs.jctc.9b00804] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
We describe a new local grand canonical Monte Carlo method to treat fluids in pores in chemical equilibrium with a reference bulk. The method is applied to Lennard-Jones particles in pores of different geometry and is shown to be much more accurate and efficient than other techniques such as traditional grand canonical simulations or Widom's particle insertion method. It utilizes a penalty potential to create a gas phase, which is in equilibrium with a more dense liquid component in the pore. Grand canonical Monte Carlo moves are employed in the gas phase, and the system then maintains chemical equilibrium by "diffusion" of particles. This creates an interface, which means that the confined fluid needs to occupy a large enough volume so that this is not an issue. We also applied the method to confined charged fluids and show how it can be used to determine local electrostatic potentials in the confined fluid, which are properly referenced to the bulk. This precludes the need to determine the Donnan potential (which controls electrochemical equilibrium) explicitly. Prior approaches have used explicit bulk simulations to measure this potential difference, which are significantly costly from a computational point of view. One outcome of our analysis is that pores of finite cross-section create a potential difference with the bulk via a small but nonzero linear charge density, which diminishes as ∼1/ln(L), where L is the pore length.
Collapse
Affiliation(s)
- Phuong Vo
- School of Science , University of New South Wales, Canberra , Canberra ACT 2600 , Australia
| | - Hongduo Lu
- Department of Theoretical Chemistry, Chemical Centre , Lund University P. O. Box 124, S-22100 Lund , Sweden
| | - Ke Ma
- School of Materials Science and Engineering , Tianjin University of Technology Tianjin 300384 , People's Republic of China
| | - Jan Forsman
- Department of Theoretical Chemistry, Chemical Centre , Lund University P. O. Box 124, S-22100 Lund , Sweden
| | - Clifford E Woodward
- School of Science , University of New South Wales, Canberra , Canberra ACT 2600 , Australia
| |
Collapse
|
3
|
Hatch HW, Hall SW, Errington JR, Shen VK. Improving the efficiency of Monte Carlo simulations of ions using expanded grand canonical ensembles. J Chem Phys 2019; 151:144109. [PMID: 31615250 PMCID: PMC7254863 DOI: 10.1063/1.5123683] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
While ionic liquids have promising applications as industrial solvents, predicting their fluid phase properties and coexistence remains a challenge. Grand canonical Monte Carlo simulation is an effective method for such predictions, but equilibration is hampered by the apparent requirement to insert and delete neutral sets of ions simultaneously in order to maintain charge neutrality. For relatively high densities and low temperatures, previously developed methods have been shown to be essential in improving equilibration by gradual insertion and deletion of these neutral sets of ions. We introduce an expanded ensemble approach which may be used in conjunction with these existing methods to further improve efficiency. Individual ions are inserted or deleted in one Monte Carlo trial rather than simultaneous insertion/deletion of neutral sets. We show how charge neutrality is maintained and show rigorous quantitative agreement between the conventional and the proposed expanded ensemble approaches, but with up to an order of magnitude increase in efficiency at high densities. The expanded ensemble approach is also more straightforward to implement than simultaneous insertion/deletion of neutral sets, and its implementation is demonstrated within open source software.
Collapse
Affiliation(s)
- Harold W. Hatch
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| | - Steven W. Hall
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC 29634, USA
| | - Jeffrey R. Errington
- Department of Chemical and Biological Engineering, University at Buffalo, Buffalo, NY 14260, USA
| | - Vincent K. Shen
- Chemical Informatics Research Group, Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8380, USA
| |
Collapse
|
4
|
Jiménez-Serratos G, Gil-Villegas A, Vega C, Blas FJ. Monte Carlo simulation of flexible trimers: From square well chains to amphiphilic primitive models. J Chem Phys 2013; 139:114901. [PMID: 24070305 DOI: 10.1063/1.4820530] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Guadalupe Jiménez-Serratos
- Departamento de Ingeniería Física, División de Ciencias e Ingenierías Campus León, Universidad de Guanajuato, Colonia Lomas del Campestre, León 37150, Mexico
| | | | | | | |
Collapse
|
5
|
Zwanikken JW, Olvera de la Cruz M. Tunable soft structure in charged fluids confined by dielectric interfaces. Proc Natl Acad Sci U S A 2013; 110:5301-8. [PMID: 23487798 PMCID: PMC3619344 DOI: 10.1073/pnas.1302406110] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Fluids of charged particles act as the supporting medium for chemical reactions and physical, dynamical, and biological processes. The local structure in an electrolytic background is deformed by micro- and nanoscopic polarizable objects. Vice versa, the forces between the objects are regulated by the cohesive properties of the background. We study here the range and strength of these forces and the microscopic origin from which they emerge. We find the forces to be sensitively dependent on the material properties of the charged fluid and the immersed solutes. The induced interactions can be varied over decades, offering high tunability and aided by accurate theory, control in experiments and applications. To distinguish correlational effects from simple ionic screening, we describe electrolyte-induced forces between neutral objects. The interplay of thermal motion, short-range repulsions, and electrostatic forces is responsible for a soft structure in the fluid. This structure changes near polarizable interfaces and causes diverse attractions between confining walls that seem well-exploited by microbiological systems. For parameters that correspond to monovalent electrolytes in biologically and technologically relevant aqueous environments, we find induced forces between nanoscopic areas of the order of piconewtons over a few nanometers.
Collapse
Affiliation(s)
| | - Monica Olvera de la Cruz
- Departments of Materials Science and Engineering and
- Chemistry, Northwestern University, Evanston, IL 60208
| |
Collapse
|
6
|
de Graaf J, Boon N, Dijkstra M, van Roij R. Electrostatic interactions between Janus particles. J Chem Phys 2012; 137:104910. [DOI: 10.1063/1.4751482] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
7
|
Bastea S. Thermodynamics and diffusion in size-symmetric and asymmetric dense electrolytes. J Chem Phys 2011; 135:084515. [PMID: 21895207 DOI: 10.1063/1.3629782] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
MD simulation results for model size-symmetric and asymmetric electrolytes at high densities and temperatures (well outside the liquid-gas coexistence region) are generated and analyzed focusing on thermodynamic and diffusion properties. An extension of the mean spherical approximation for electrolytes originally derived for charged hard sphere fluids is adapted to these systems by exploiting the separation of short range and Coulomb interaction contributions intrinsic to these theoretical models and is found to perform well for predicting equation of state quantities. The diffusion coefficients of these electrolytes can also be reasonably well predicted using entropy scaling ideas suitably adapted to charged systems and mixtures. Thus, this approach may provide an avenue for studying dense electrolytes or complex molecular systems containing charged species at high pressures and temperatures.
Collapse
Affiliation(s)
- Sorin Bastea
- Lawrence Livermore National Laboratory, 7000 East Ave., Livermore, California 94550, USA.
| |
Collapse
|
8
|
Zwanikken JW, Jha PK, de la Cruz MO. A practical integral equation for the structure and thermodynamics of hard sphere Coulomb fluids. J Chem Phys 2011; 135:064106. [DOI: 10.1063/1.3624809] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Coslovich D, Hansen JP, Kahl G. Ultrasoft primitive model of polyionic solutions: Structure, aggregation, and dynamics. J Chem Phys 2011; 134:244514. [DOI: 10.1063/1.3602469] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Zwanikken JW, de la Cruz MO. Correlated electrolyte solutions and ion-induced attractions between nanoparticles. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:050401. [PMID: 21230424 DOI: 10.1103/physreve.82.050401] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/10/2010] [Indexed: 05/30/2023]
Abstract
Simple expressions are presented for the equations of state of a correlated electrolyte solution, calculated straightforwardly within a full nonlinear Debye-Hückel approach in terms of the mean potential at contact, that predict quantitatively different phase behavior from the popular Debye-Hückel limiting law. The theory includes pair correlations accurately and may provide a basis for a quantitative theoretical study of organic or multivalent ionic solutions. As an example, cohesive effects are addressed of strong couplings between ions on the effective interactions between nanoparticles.
Collapse
Affiliation(s)
- Jos W Zwanikken
- Department of Materials Science and Engineering, Northwestern University, 2220 Campus Drive, Evanston, Illinois 60208, USA
| | | |
Collapse
|