1
|
Marshall BD. A modified perturbed chain‐statistical associating fluid theory equation of state for water which includes an association dependent hard sphere diameter. AIChE J 2021. [DOI: 10.1002/aic.17342] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
2
|
van Westen T. Algebraic second virial coefficient of the Mie m - 6 intermolecular potential based on perturbation theory. J Chem Phys 2021; 154:234502. [PMID: 34241261 DOI: 10.1063/5.0050659] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We propose several simple algebraic approximations for the second virial coefficient of fluids whose molecules interact by a generic Mie m - 6 intermolecular pair potential. In line with a perturbation theory, the parametric equations are formulated as the sum of a contribution due to a reference part of the intermolecular potential and a perturbation. Thereby, the equations provide a convenient (low-density) starting point for developing equation-of-state models of fluids or for developing similar approximations for the virial coefficient of (polymeric-)chain fluids. The choice of Barker and Henderson [J. Chem. Phys. 47, 4714 (1967)] and Weeks, Chandler, and Andersen [Phys. Rev. Lett. 25, 149 (1970); J. Chem. Phys. 54, 5237 (1971); and Phys. Rev. A 4, 1597 (1971)] for the reference part of the potential is considered. Our analytic approximations correctly recover the virial coefficient of the inverse-power potential of exponent m in the high-temperature limit and provide accurate estimates of the temperatures for which the virial coefficient equals zero or takes on its maximum value. Our description of the reference contribution to the second virial coefficient follows from an exact mapping onto the second virial coefficient of hard spheres; we propose a simple algebraic equation for the corresponding effective diameter of the hard spheres, which correctly recovers the low- and high-temperature scaling and limits of the reference fluid's second virial coefficient.
Collapse
Affiliation(s)
- Thijs van Westen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, D-70569 Stuttgart, Germany
| |
Collapse
|
3
|
Karnes JJ, Benjamin I. Deconstructing the Local Intermolecular Ordering and Dynamics of Liquid Chloroform and Bromoform. J Phys Chem B 2021; 125:3629-3637. [PMID: 33792320 DOI: 10.1021/acs.jpcb.0c10407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Local intermolecular structure and dynamics of the polar molecular liquids chloroform and bromoform are studied by molecular dynamics simulation. Structural distribution functions, including 1- and 2-D pair correlations and dipole contour plots allow direct comparison and show agreement with recent analyses of diffraction experiments. Studies of the haloforms' reorientational dynamics and longevity of structural features resulting from intermolecular interaction extend previous work toward deeper understanding of the factors controlling these features. Analyses of ensemble average structures and dynamical properties isolate mass, electrostatics, and steric packing as driving forces or contributing factors for the observed ordering and dynamics.
Collapse
Affiliation(s)
- John J Karnes
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Ilan Benjamin
- Department of Chemistry and Biochemistry, University of California-Santa Cruz, Santa Cruz, California 95064, United States
| |
Collapse
|
4
|
van Westen T, Gross J. Accurate first-order perturbation theory for fluids: uf-theory. J Chem Phys 2021; 154:041102. [PMID: 33514104 DOI: 10.1063/5.0031545] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
We propose a new first-order perturbation theory that provides a near-quantitative description of the thermodynamics of simple fluids. The theory is based on the ansatz that the Helmholtz free energy is bounded below by a first-order Mayer-f expansion. Together with the rigorous upper bound provided by a first-order u-expansion, this brackets the actual free energy between an upper and (effective) lower bound that can both be calculated based on first-order perturbation theory. This is of great practical use. Here, the two bounds are combined into an interpolation scheme for the free energy. The scheme exploits the fact that a first-order Mayer-f perturbation theory is exact in the low-density limit, whereas the accuracy of a first-order u-expansion grows when density increases. This allows an interpolation between the lower "f"-bound at low densities and the upper "u" bound at higher liquid-like densities. The resulting theory is particularly well behaved. Using a density-dependent interpolating function of only two adjustable parameters, we obtain a very accurate representation of the full fluid-phase behavior of a Lennard-Jones fluid. The interpolating function is transferable to other intermolecular potential types, which is here shown for the Mie m-6 family of fluids. The extension to mixtures is simple and accurate without requiring any dependence of the interpolating function on the composition of the mixture.
Collapse
Affiliation(s)
- Thijs van Westen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, D-70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, D-70569 Stuttgart, Germany
| |
Collapse
|
5
|
Ghorai PK, Matyushov DV. Equilibrium Solvation, Electron-Transfer Reactions, and Stokes-Shift Dynamics in Ionic Liquids. J Phys Chem B 2020; 124:3754-3769. [DOI: 10.1021/acs.jpcb.0c01773] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Pradip Kr. Ghorai
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur 741246, India
| | - Dmitry V. Matyushov
- Department of Physics and School of Molecular Sciences, Arizona State University, PO Box 871504, Tempe, Arizona 85287, United States
| |
Collapse
|
6
|
D'Urso C, Celebre G, Cinacchi G. Phase behavior of hard C_{2h}-symmetric particle systems. Phys Rev E 2019; 100:012709. [PMID: 31499787 DOI: 10.1103/physreve.100.012709] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Indexed: 06/10/2023]
Abstract
Using Monte Carlo numerical simulation, this work sketches the phase diagram of systems of certain hard C_{2h}-symmetric particles, formed by gluing two aligned and displaced hard spherocylinders with a cylindrical-length-to-diameter ratio realistically, if viewed not only from the lyotropic colloidal liquid-crystal side but also from the thermotropic low-molecular-mass liquid-crystal side, equal to 5, as a function of the displacement. Several distinctive phases are observed, such as a nonperiodic smectic-B-like phase, a nonperiodic smectic-H-like phase, a smectic-C phase, and a short-layer-spacing uniaxial smectic-A phase but no biaxial nematic phase.
Collapse
Affiliation(s)
- Christian D'Urso
- Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, via P. Bucci, I-87036 Arcavacata di Rende (CS), Italy
| | - Giorgio Celebre
- Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, via P. Bucci, I-87036 Arcavacata di Rende (CS), Italy
| | - Giorgio Cinacchi
- Departamento de Física Teórica de la Materia Condensada, Instituto de Física de la Materia Condensada (IFIMAC), Instituto de Ciencias de Materiales "Nicolás Cabrera," Universidad Autónoma de Madrid, Ciudad Universitaria de Cantoblanco, E-28049 Madrid, Spain
| |
Collapse
|
7
|
Altabet YE, Fenley AL, Stillinger FH, Debenedetti PG. Cavitation transition in the energy landscape: Distinct tensile yielding behavior in strongly and weakly attractive systems. J Chem Phys 2018; 148:114501. [DOI: 10.1063/1.5019274] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Y. Elia Altabet
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Andreia L. Fenley
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Frank H. Stillinger
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Pablo G. Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
8
|
Gyawali G, Sternfield S, Kumar R, Rick SW. Coarse-Grained Models of Aqueous and Pure Liquid Alkanes. J Chem Theory Comput 2017. [DOI: 10.1021/acs.jctc.7b00389] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Gaurav Gyawali
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Samuel Sternfield
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Revati Kumar
- Department
of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70808, United States
| | - Steven W. Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| |
Collapse
|
9
|
de Souza RF, Zannoni C. Rotational diffusion of shape switching particles in nematic liquid crystals. Phys Rev E 2017; 94:062702. [PMID: 28085296 DOI: 10.1103/physreve.94.062702] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2016] [Indexed: 11/07/2022]
Abstract
The theory of rotational diffusion of particles of various symmetry embedded in a liquid crystal host, essential to interpret a variety of spectroscopic observables, has been available for some time, but only for the case of rigid molecules. Here we generalize the treatment and present a theory to describe the rotational diffusion of shape-changing particles dispersed in nematic liquid crystals. The interaction of the particles with the environment is modeled by an effective field potential, while the particles are allowed to assume an arbitrary discrete number of shapes. The transition between shapes is modeled by a Markovian process which is combined with rotational diffusion. Our model is applied to the simple case of a particle that can exchange between three shapes: a rod, a disk, and a sphere. We consider in detail the effect of shape transitions in some selected correlation functions which are relevant for experiments.
Collapse
Affiliation(s)
- R F de Souza
- Dipartimento di Chimica Industriale "Toso Montanari" and INSTM, Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy.,CAPES Foundation, Ministry of Education of Brazil, Brasília, DF 70040-020, Brazil
| | - C Zannoni
- Dipartimento di Chimica Industriale "Toso Montanari" and INSTM, Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| |
Collapse
|
10
|
Dunn NJH, Foley TT, Noid WG. Van der Waals Perspective on Coarse-Graining: Progress toward Solving Representability and Transferability Problems. Acc Chem Res 2016; 49:2832-2840. [PMID: 27993007 DOI: 10.1021/acs.accounts.6b00498] [Citation(s) in RCA: 71] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Low-resolution coarse-grained (CG) models provide the necessary efficiency for simulating phenomena that are inaccessible to more detailed models. However, in order to realize their considerable promise, CG models must accurately describe the relevant physical forces and provide useful predictions. By formally integrating out the unnecessary details from an all-atom (AA) model, "bottom-up" approaches can, at least in principle, quantitatively reproduce the structural and thermodynamic properties of the AA model that are observable at the CG resolution. In practice, though, bottom-up approaches only approximate this "exact coarse-graining" procedure. The resulting models typically reproduce the intermolecular structure of AA models at a single thermodynamic state point but often describe other state points less accurately and, moreover, tend to provide a poor description of thermodynamic properties. These two limitations have been coined the "transferability" and "representability" problems, respectively. Perhaps, the simplest and most commonly discussed manifestation of the representability problem regards the tendency of structure-based CG models to dramatically overestimate the pressure. Furthermore, when these models are adjusted to reproduce the pressure, they provide a poor description of the compressibility. More generally, it is sometimes suggested that CG models are fundamentally incapable of reproducing both structural and thermodynamic properties. After all, there is no such thing as a "free lunch"; any significant gain in computational efficiency should come at the cost of significant model limitations. At least in the case of structural and thermodynamic properties, though, we optimistically propose that this may be a false dichotomy. Accordingly, we have recently re-examined the "exact coarse-graining" procedure and investigated the intrinsic consequences of representing an AA model in reduced resolution. These studies clarify the origin and inter-relationship of representability and transferability problems. Both arise as consequences of transferring thermodynamic information from the high resolution configuration space and encoding this information into the many-body potential of mean force (PMF), that is, the potential that emerges from an exact coarse-graining procedure. At least in principle, both representability and transferability problems can be resolved by properly addressing this thermodynamic information. In particular, we have demonstrated that "pressure-matching" provides a practical and rigorous means for addressing the density dependence of the PMF. The resulting bottom-up models accurately reproduce the structure, equilibrium density, compressibility, and pressure equation of state for AA models of molecular liquids. Additionally, we have extended this approach to develop transferable potentials that provide similar accuracy for heptane-toluene mixtures. Moreover, these potentials provide predictive accuracy for modeling concentrations that were not considered in their parametrization. More generally, this work suggests a "van der Waals" perspective on coarse-graining, in which conventional structure-based methods accurately describe the configuration dependence of the PMF, while independent variational principles infer the thermodynamic information that is necessary to resolve representability and transferability problems.
Collapse
Affiliation(s)
- Nicholas J. H. Dunn
- Department
of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Thomas T. Foley
- Department
of Physics, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - William G. Noid
- Department
of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
11
|
Altabet YE, Stillinger FH, Debenedetti PG. A cavitation transition in the energy landscape of simple cohesive liquids and glasses. J Chem Phys 2016; 145:211905. [DOI: 10.1063/1.4959846] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Y. Elia Altabet
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Frank H. Stillinger
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Pablo G. Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
12
|
Glover WJ, Schwartz BJ. Short-Range Electron Correlation Stabilizes Noncavity Solvation of the Hydrated Electron. J Chem Theory Comput 2016; 12:5117-5131. [PMID: 27576177 DOI: 10.1021/acs.jctc.6b00472] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The hydrated electron, e-(aq), has often served as a model system to understand the influence of condensed-phase environments on electronic structure and dynamics. Despite over 50 years of study, however, the basic structure of e-(aq) is still the subject of controversy. In particular, the structure of e-(aq) was long assumed to be an electron localized within a solvent cavity, in a manner similar to halide solvation. Recently, however, we suggested that e-(aq) occupies a region of enhanced water density with little or no discernible cavity. The potential we developed was only subtly different from those that give rise to a cavity solvation motif, which suggests that the driving forces for noncavity solvation involve subtle electron-water attractive interactions at close distances. This leads to the question of how dispersion interactions are treated in simulations of the hydrated electron. Most dispersion potentials are ad hoc or are not designed to account for the type of close-contact electron-water overlap that might occur in the condensed phase, and where short-range dynamic electron correlation is important. To address this, in this paper we develop a procedure to calculate the potential energy surface between a single water molecule and an excess electron with high-level CCSD(T) electronic structure theory. By decomposing the electron-water potential into its constituent energetic contributions, we find that short-range electron correlation provides an attraction of comparable magnitude to the mean-field interactions between the electron and water. Furthermore, we find that by reoptimizing a popular cavity-forming one-electron model potential to better capture these attractive short-range interactions, the enhanced description of correlation predicts a noncavity e-(aq) with calculated properties in better agreement with experiment. Although much attention has been placed on the importance of long-range dispersion interactions in water cluster anions, our study reveals that largely unexplored short-range correlation effects are crucial in dictating the solvation structure of the condensed-phase hydrated electron.
Collapse
Affiliation(s)
- William J Glover
- NYU-ECNU Center for Computational Chemistry, New York University Shanghai , Shanghai, 200122, China.,Department of Chemistry, New York University , New York, New York 10003, United States
| | - Benjamin J Schwartz
- Department of Chemistry and Biochemistry, University of California, Los Angeles , Los Angeles, California 90095, United States
| |
Collapse
|
13
|
Kreis K, Potestio R. The relative entropy is fundamental to adaptive resolution simulations. J Chem Phys 2016; 145:044104. [DOI: 10.1063/1.4959169] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Karsten Kreis
- Max-Planck-Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany
- Graduate School Materials Science in Mainz, Staudingerweg 9, 55128 Mainz, Germany
| | - Raffaello Potestio
- Max-Planck-Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
14
|
Dunn NJH, Noid WG. Bottom-up coarse-grained models that accurately describe the structure, pressure, and compressibility of molecular liquids. J Chem Phys 2015; 143:243148. [DOI: 10.1063/1.4937383] [Citation(s) in RCA: 76] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Nicholas J. H. Dunn
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - W. G. Noid
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
15
|
Hlushak S. Perturbation theory for multicomponent fluids based on structural properties of hard-sphere chain mixtures. J Chem Phys 2015; 143:124906. [PMID: 26429040 DOI: 10.1063/1.4931816] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
An analytical expression for the Laplace transform of the radial distribution function of a mixture of hard-sphere chains of arbitrary segment size and chain length is used to rigorously formulate the first-order Barker-Henderson perturbation theory for the contribution of the segment-segment dispersive interactions into thermodynamics of the Lennard-Jones chain mixtures. Based on this approximation, a simple variant of the statistical associating fluid theory is proposed and used to predict properties of several mixtures of chains of different lengths and segment sizes. The theory treats the dispersive interactions more rigorously than the conventional theories and provides means for more accurate description of dispersive interactions in the mixtures of highly asymmetric components.
Collapse
Affiliation(s)
- Stepan Hlushak
- Institute for Condensed Matter Physics, Svientsitskoho 1, 79011 Lviv, Ukraine
| |
Collapse
|
16
|
Abstract
A theory based on the exponential approximation of the liquid-state theory is applied to study properties of several models of one-component Yukawa plasma characterized by different values of the screening parameter z. The results of the new theory are compared to the results of a conventional theory, which is based on the first-order mean spherical approximation, and to the results of a Monte Carlo simulation. The new theory shows improvements in the predictions for the thermodynamic and structural properties of Yukawa plasmas with high and intermediate values of the screening parameter, z, and coupling parameter, Γ. For low values of z and Γ, the new theory is comparable in accuracy to the conventional theory, which in turn agrees well with the results of the Monte Carlo simulation.
Collapse
Affiliation(s)
- Stepan Hlushak
- Institute for Condensed Matter Physics, Svientsitskoho 1, 79011 Lviv, Ukraine and Department of Chemical and Materials Engineering, University of Alberta, 9107-116 Street, Edmonton, Alberta T6G 2V4, Canada
| |
Collapse
|
17
|
Abstract
By focusing on essential features, while averaging over less important details, coarse-grained (CG) models provide significant computational and conceptual advantages with respect to more detailed models. Consequently, despite dramatic advances in computational methodologies and resources, CG models enjoy surging popularity and are becoming increasingly equal partners to atomically detailed models. This perspective surveys the rapidly developing landscape of CG models for biomolecular systems. In particular, this review seeks to provide a balanced, coherent, and unified presentation of several distinct approaches for developing CG models, including top-down, network-based, native-centric, knowledge-based, and bottom-up modeling strategies. The review summarizes their basic philosophies, theoretical foundations, typical applications, and recent developments. Additionally, the review identifies fundamental inter-relationships among the diverse approaches and discusses outstanding challenges in the field. When carefully applied and assessed, current CG models provide highly efficient means for investigating the biological consequences of basic physicochemical principles. Moreover, rigorous bottom-up approaches hold great promise for further improving the accuracy and scope of CG models for biomolecular systems.
Collapse
Affiliation(s)
- W G Noid
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
18
|
Urbic T, Dias CL. Properties of the Lennard-Jones dimeric fluid in two dimensions: an integral equation study. J Chem Phys 2014; 140:094703. [PMID: 24606372 PMCID: PMC3977749 DOI: 10.1063/1.4867289] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2014] [Accepted: 02/19/2014] [Indexed: 12/21/2022] Open
Abstract
The thermodynamic and structural properties of the planar soft-sites dumbbell fluid are examined by Monte Carlo simulations and integral equation theory. The dimers are built of two Lennard-Jones segments. Site-site integral equation theory in two dimensions is used to calculate the site-site radial distribution functions for a range of elongations and densities and the results are compared with Monte Carlo simulations. The critical parameters for selected types of dimers were also estimated. We analyze the influence of the bond length on critical point as well as tested correctness of site-site integral equation theory with different closures. The integral equations can be used to predict the phase diagram of dimers whose molecular parameters are known.
Collapse
Affiliation(s)
- Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Aškerčeva 5, SI-1000 Ljubljana, Slovenia
| | - Cristiano L Dias
- Physics Department, New Jersey Institute of Technology, Newark, New Jersey 07102-1982, USA
| |
Collapse
|
19
|
Hlushak SP, Hlushak PA, Trokhymchuk A. An improved first-order mean spherical approximation theory for the square-shoulder fluid. J Chem Phys 2013; 138:164107. [PMID: 23635111 DOI: 10.1063/1.4801659] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
The theory, which utilizes an exponential enhancement of the first-order mean spherical approximation (FMSA) for the radial distribution functions of the hard-core plus square-well fluid, is adopted to study the properties of the simplest model of the core-softened fluids, i.e., the hard spheres with a square-shoulder interaction. The results for structure and thermodynamic properties are reported and compared against both the Monte Carlo simulation data as well as with those obtained within the conventional FMSA theory. We found that in the region of low densities and low temperatures, where the conventional FMSA theory fails, the exponential-based FMSA theory besides being qualitatively correct also provides with a notable quantitative improvement of the theoretical description.
Collapse
Affiliation(s)
- S P Hlushak
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37235, USA.
| | | | | |
Collapse
|
20
|
Abstract
This chapter provides a primer on theories for coarse-grained (CG) modeling and, in particular, reviews several systematic methods for determining effective potentials for CG models. The chapter first reviews a statistical mechanics framework for relating atomistic and CG models. This framework naturally leads to a quantitative criterion for CG models that are "consistent" with a particular atomistic model for the same system. This consistency criterion is equivalent to minimizing the relative entropy between the two models. This criterion implies that a many-body PMF is the appropriate potential for a CG model that is consistent with a particular atomistic model. This chapter then presents a unified exposition of the theory and numerical methods for several approaches for approximating this many-body PMF. Finally, this chapter closes with a brief discussion of a few of the outstanding challenges facing the field of systematic coarse-graining.
Collapse
Affiliation(s)
- W G Noid
- Department of Chemistry, The Pennsylvania State University, University Park, PA, USA
| |
Collapse
|
21
|
Larsen RE, Glover WJ, Schwartz BJ. Response to Comments on “Does the Hydrated Electron Occupy a Cavity?”. Science 2011. [DOI: 10.1126/science.1197884] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Affiliation(s)
- Ross E. Larsen
- Center for Scientific Computing, National Renewable Energy Laboratory, Golden, CO 80401–3305, USA
| | - William J. Glover
- Department of Chemistry, Stanford University, Stanford, CA 94305–5080, USA
| | - Benjamin J. Schwartz
- Department of Chemistry and Biochemistry, University of California–Los Angeles, Los Angeles, CA 90095–1569, USA
| |
Collapse
|
22
|
Mullinax JW, Noid WG. Reference state for the generalized Yvon-Born-Green theory: application for coarse-grained model of hydrophobic hydration. J Chem Phys 2010; 133:124107. [PMID: 20886924 PMCID: PMC3188631 DOI: 10.1063/1.3481574] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2010] [Accepted: 08/02/2010] [Indexed: 01/07/2023] Open
Abstract
Coarse-grained (CG) models provide a computationally efficient means for investigating phenomena that remain beyond the scope of atomically detailed models. Although CG models are often parametrized to reproduce the results of atomistic simulations, it is highly desirable to determine accurate CG models from experimental data. Recently, we have introduced a generalized Yvon-Born-Green (g-YBG) theory for directly (i.e., noniteratively) determining variationally optimized CG potentials from structural correlation functions. In principle, these correlation functions can be determined from experiment. In the present work, we introduce a reference state potential into the g-YBG framework. The reference state defines a fixed contribution to the CG potential. The remaining terms in the potential are then determined, such that the combined potential provides an optimal approximation to the many-body potential of mean force. By specifying a fixed contribution to the potential, the reference state significantly reduces the computational complexity and structural information necessary for determining the remaining potentials. We also validate the quantitative accuracy of the proposed method and numerically demonstrate that the reference state provides a convenient framework for transferring CG potentials from neat liquids to more complex systems. The resulting CG model provides a surprisingly accurate description of the two- and three-particle solvation structures of a hydrophobic solute in methanol. This work represents a significant step in developing the g-YBG theory as a useful computational framework for determining accurate CG models from limited experimental data.
Collapse
Affiliation(s)
- J W Mullinax
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | | |
Collapse
|
23
|
|
24
|
Naumann KH, Chen YP, Leland TW. Conformal-Solution-Theorie für Mischungen konvexer Moleküle. ACTA ACUST UNITED AC 2010. [DOI: 10.1002/bbpc.19810851119] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
25
|
Re M, Laria D, Fernández-Prini R. The role of solvent structure in perturbation methods applied to the dissolution of an apolar solute in water. ACTA ACUST UNITED AC 2010. [DOI: 10.1002/bbpc.19961000810] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
26
|
Matyushov DV, Schmid R. Properties of liquids at the boiling point: Equation of state, internal pressure and vaporization entropy. ACTA ACUST UNITED AC 2010. [DOI: 10.1002/bbpc.19940981214] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
27
|
Yoshimura Y, Nakahara M. Additivity Rule of the Partial Molal Volume. 4. Its Proof Based on the Postulate of Independent Volume Changes for Unit Reactions. ACTA ACUST UNITED AC 2010. [DOI: 10.1002/bbpc.19860900109] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
28
|
Athawale MV, Jamadagni SN, Garde S. How hydrophobic hydration responds to solute size and attractions: Theory and simulations. J Chem Phys 2009; 131:115102. [DOI: 10.1063/1.3227031] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
29
|
Morozov KI. Soret effect in molecular mixtures. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 79:031204. [PMID: 19391928 DOI: 10.1103/physreve.79.031204] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2009] [Indexed: 05/27/2023]
Abstract
A theoretical approach to the description of the Soret effect in binary nonpolar liquids is proposed. The temperature gradient of the partial pressure is determined as the driving force of thermal diffusion. The hard-sphere fluid is chosen as a reference system and an explicit relation for the Soret coefficient is found. Two additive contributions owing to steric repulsions and attractive interactions form the so-called chemical contribution to S_{T} [C. Debuschewitz and W. Köhler, Phys. Rev. Lett. 87, 055901 (2001)]. The parameters of interparticle interactions are defined with the help of the solvation theory. In particular, the van der Waals constant of cross interactions is expressed via the excess volume of mixture. The proposed theory is applied to the benzene-cyclohexane system. A reasonable agreement of theoretical and experimental results is revealed for the Soret coefficient and its temperature dependence.
Collapse
Affiliation(s)
- Konstantin I Morozov
- Institute of Continuous Media Mechanics, Ural Branch of RAS, Perm 614013, Russia.
| |
Collapse
|
30
|
Tsai KH, Wu TM. Hard-sphere perturbation theory for a model of liquid Ga. J Chem Phys 2008; 129:024503. [DOI: 10.1063/1.2948950] [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
|
31
|
LeBard DN, Matyushov DV. Redox entropy of plastocyanin: Developing a microscopic view of mesoscopic polar solvation. J Chem Phys 2008; 128:155106. [DOI: 10.1063/1.2904879] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Matyushov DV. Model energy landscapes of low-temperature fluids: Dipolar hard spheres. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 76:011511. [PMID: 17677454 DOI: 10.1103/physreve.76.011511] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/02/2007] [Revised: 03/05/2007] [Indexed: 05/16/2023]
Abstract
An analytical model of non-Gaussian energy landscape of low-temperature fluids is developed based on the thermodynamics of the fluid of dipolar hard spheres. The entire excitation profile of the liquid, from the high-temperature liquid to the point of ideal-glass transition, has been obtained from Monte Carlo simulations. The fluid of dipolar hard spheres loses stability close to the point of ideal-glass transition transforming via a first-order transition into a columnar liquid phase of dipolar chains locally arranged in a body-centered-tetragonal order. Significant non-Gaussianity of the energy landscape is responsible for narrowing of the distribution of potential energies and energies of inherent structures with decreasing temperature. We suggest that the proposed functionality of the enumeration function is widely applicable to both polar and nonpolar low-temperature liquids.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- Center for Biological Physics, Arizona State University, P.O. Box 871604, Tempe, Arizona 85287-1604, USA.
| |
Collapse
|
33
|
Chakrabarti D, Jana B, Bagchi B. Orientational relaxation in a discotic liquid crystal. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 75:061703. [PMID: 17677279 DOI: 10.1103/physreve.75.061703] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2006] [Indexed: 05/16/2023]
Abstract
We investigate orientational relaxation of a model discotic liquid crystal, consisting of disclike molecules, by molecular dynamics simulations along two isobars starting from the high temperature isotropic phase. The two isobars have been so chosen that (a) the phase sequence isotropic- (I-) nematic- (N-) columnar (C) appears upon cooling along one of them and (b) the sequence isotropic- (I-) columnar- (C) along the other. While the orientational relaxation in the isotropic phase near the I-N phase transition in system (a) shows a power law decay at short to intermediate times, such power law relaxation is not observed in the isotropic phase near the I-C phase boundary in system (b). In order to understand this difference (the existence or the absence of the power law decay), we calculated the growth of the orientational pair distribution functions (OPDFs) near the I-N phase boundary and also near the I-C phase boundary. We find that the OPDF shows a marked growth in long range correlation as the I-N phase boundary is approached in the I-N-C system (a), but such a growth is absent in the I-C system, which appears to be consistent with the result that I-N phase transition in the former is weakly first order while the I-C phase transition in the latter is not weak. As the system settles into the nematic phase, the decay of the single-particle second-rank orientational time correlation function follows a pattern that is similar to what is observed with calamitic liquid crystals and supercooled molecular liquids.
Collapse
Affiliation(s)
- Dwaipayan Chakrabarti
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | | | | |
Collapse
|
34
|
Morozov KI. The dielectric virial expansion and the models of dipolar hard-sphere fluid. J Chem Phys 2007; 126:194506. [PMID: 17523821 DOI: 10.1063/1.2736370] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The virial expansion technique to determine the dielectric constant epsilon of dipolar hard-sphere fluid is developed. It is shown that the formalism allows to bring into agreement the results of Debye's, Onsager's, and Langevin's to the problem. The third virial coefficient of epsilon is considered as a series over dipolar parameter lambda=m(2)d(3)kT. The terms up to O(lambda(11)) are calculated analytically providing a correct description of the third virial coefficient for small and intermediate values of lambda (0<or=lambda<or=4). The results of the dielectric virial series are compared with the Monte Carlo data for epsilon found by Matyushov and Ladanyi [J. Chem. Phys. 110, 994 (1999)]. The theory is in agreement with simulations only at small values of lambda<or=2. At higher polarities, the virial series diverges. Realization of the renormalization procedure permits to enlarge the range of applicability of the virial series. In this way, the new expression for the dielectric constant as a function of two dipolar parameters, lambda and y=4 pi nm(2)9kT, has been found explicitly. The expression gives a perfect upper bound of the dielectric constant and is more reliable for determination of epsilon than the previously known ones.
Collapse
Affiliation(s)
- Konstantin I Morozov
- Institute of Continuous Media Mechanics, UB of Russian Academy of Sciences, 614013 Perm, Russia.
| |
Collapse
|
35
|
Parsafar GA, Shokouhi M. Modification of a new potential model used for calculation of the second virial coefficient and zero density transport properties. Mol Phys 2006. [DOI: 10.1080/00268970601082887] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
36
|
Affiliation(s)
- Bhalachandra L. Tembe
- a Department of Chemistry and Radiation Laboratory , University of Notre Dame , Notre Dame , Indiana , 46556 , U.S.A
| | - John J. Kozak
- a Department of Chemistry and Radiation Laboratory , University of Notre Dame , Notre Dame , Indiana , 46556 , U.S.A
| |
Collapse
|
37
|
Renzi DG, Stoico CO, Vericat F. Hydration properties and potentials of mean force of nonpolar amino acid residues in water: A pertubation theoretic approach. J Chem Phys 2005; 123:104502. [PMID: 16178605 DOI: 10.1063/1.2039089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We use a perturbation approach to calculate the solute-solvent and solute-solute pair correlations for infinitely dilute solutions of nonpolar amino acid side chains in water. In particular, from the solute-solute correlations we derive potentials of mean force for all the distances between different residues of amino acids. Comparison with molecular dynamics simulations performed using GROMOS package shows that the goodness of the approximation varies with the amino acid considered. We discuss in terms of hydrophobicity the different effect that the inclusion of the attractions causes on the solute-water and solute-solute correlations.
Collapse
Affiliation(s)
- Danilo G Renzi
- Facultad de Ciencias Veterinarias, Universidad Nacional de Rosario (UNR), Casilda, Argentina
| | | | | |
Collapse
|
38
|
Ripoll M, Mussawisade K, Winkler RG, Gompper G. Dynamic regimes of fluids simulated by multiparticle-collision dynamics. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 72:016701. [PMID: 16090128 DOI: 10.1103/physreve.72.016701] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2005] [Indexed: 05/03/2023]
Abstract
We investigate the hydrodynamic properties of a fluid simulated with a mesoscopic solvent model. Two distinct regimes are identified, the "particle regime" in which the dynamics is gaslike and the "collective regime" where the dynamics is fluidlike. This behavior can be characterized by the Schmidt number, which measures the ratio between viscous and diffusive transport. Analytical expressions for the tracer diffusion coefficient, which have been derived on the basis of a molecular-chaos assumption, are found to describe the simulation data very well in the particle regime, but important deviations are found in the collective regime. These deviations are due to hydrodynamic correlations. The model is then extended in order to investigate self-diffusion in colloidal dispersions. We study first the transport properties of heavy pointlike particles in the mesoscopic solvent, as a function of their mass and number density. Second, we introduce excluded-volume interactions among the colloidal particles and determine the dependence of the diffusion coefficient on the colloidal volume fraction for different solvent mean-free paths. In the collective regime, the results are found to be in good agreement with previous theoretical predictions based on Stokes hydrodynamics and the Smoluchowski equation.
Collapse
Affiliation(s)
- M Ripoll
- Institut für Festkörperforschung, Forschungszentrum Jülich, D-52425 Jülich, Germany
| | | | | | | |
Collapse
|
39
|
|
40
|
Ben-Amotz D, Stell G. Hard sphere perturbation theory for fluids with soft-repulsive-core potentials. J Chem Phys 2004; 120:4844-51. [PMID: 15267344 DOI: 10.1063/1.1647520] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The thermodynamic properties of fluids with very soft repulsive-core potentials, resembling those of some liquid metals, are predicted with unprecedented accuracy using a new first-order thermodynamic perturbation theory. This theory is an extension of Mansoori-Canfield/Rasaiah-Stell (MCRS) perturbation theory, obtained by including a configuration integral correction recently identified by Mon, who evaluated it by computer simulation. In this work we derive an analytic expression for Mon's correction in terms of the radial distribution function of the soft-core fluid, g(0)(r), approximated using Lado's self-consistent extension of Weeks-Chandler-Andersen (WCA) theory. Comparisons with WCA and MCRS predictions show that our new extended-MCRS theory outperforms other first-order theories when applied to fluids with very soft inverse-power potentials (n< or =6), and predicts free energies that are within 0.3 kT of simulation results up to the fluid freezing point.
Collapse
Affiliation(s)
- Dor Ben-Amotz
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | | |
Collapse
|
41
|
|
42
|
|
43
|
Evans R, Schirmacher W. The long wavelength limit of the liquid structure factor-a theory for the compressibility of liquid rare gases and metals. ACTA ACUST UNITED AC 2001. [DOI: 10.1088/0022-3719/11/12/007] [Citation(s) in RCA: 63] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
44
|
Bratkovsky AM, Vaks VG, Trefilov AV. On the accuracy of the liquid theory approximate methods for the description of liquid metal thermodynamics. ACTA ACUST UNITED AC 2000. [DOI: 10.1088/0305-4608/13/12/010] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
45
|
|
46
|
Evans R, Sluckin TJ. The role of attractive forces in the structure of simple liquids: a theory for small-angle scattering. ACTA ACUST UNITED AC 2000. [DOI: 10.1088/0022-3719/14/19/008] [Citation(s) in RCA: 51] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
47
|
Hoshino K. Effect of the softness of the pair potential on the liquid structure factor. ACTA ACUST UNITED AC 2000. [DOI: 10.1088/0022-3719/13/17/003] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
48
|
MacGowan D. A discontinuous mean spherical approximation for the strongly coupled one-component plasma. ACTA ACUST UNITED AC 2000. [DOI: 10.1088/0022-3719/16/1/010] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
49
|
|
50
|
Kahl G, Hafner J. A blip-function calculation of the structure of liquid binary alloys. ACTA ACUST UNITED AC 2000. [DOI: 10.1088/0305-4608/15/8/003] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|