1
|
Baptista LA, Sevilla M, Wagner M, Kremer K, Cortes-Huerto R. Chilling alcohol on the computer: isothermal compressibility and the formation of hydrogen-bond clusters in liquid propan-1-ol. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2023; 46:117. [PMID: 38019330 PMCID: PMC10687148 DOI: 10.1140/epje/s10189-023-00380-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/22/2023] [Accepted: 11/15/2023] [Indexed: 11/30/2023]
Abstract
Molecular dynamics simulations have been performed to compute the isothermal compressibility [Formula: see text] of liquid propan-1-ol in the temperature range [Formula: see text] K. A change in behaviour, from normal (high T) to anomalous (low T), has been identified for [Formula: see text] at [Formula: see text] K. The average number of hydrogen bonds (H-bond) per molecule turns to saturation in the same temperature interval, suggesting the formation of a relatively rigid network. Indeed, simulation results show a strong tendency to form H-bond clusters with distinct boundaries, with the average largest size and width of the size distribution growing upon decreasing temperature, in agreement with previous theoretical and experimental studies. These results also emphasise a connection between the behaviour of [Formula: see text] and the formation of nanometric structures.
Collapse
Affiliation(s)
- Luis A Baptista
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - Mauricio Sevilla
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - Manfred Wagner
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | | |
Collapse
|
2
|
Simon JM, Krüger P, Schnell SK, Vlugt TJH, Kjelstrup S, Bedeaux D. Kirkwood-Buff integrals: From fluctuations in finite volumes to the thermodynamic limit. J Chem Phys 2022; 157:130901. [PMID: 36209013 DOI: 10.1063/5.0106162] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The Kirkwood-Buff theory is a cornerstone of the statistical mechanics of liquids and solutions. It relates volume integrals over the radial distribution function, so-called Kirkwood-Buff integrals (KBIs), to particle number fluctuations and thereby to various macroscopic thermodynamic quantities such as the isothermal compressibility and partial molar volumes. Recently, the field has seen a strong revival with breakthroughs in the numerical computation of KBIs and applications to complex systems such as bio-molecules. One of the main emergent results is the possibility to use the finite volume KBIs as a tool to access finite volume thermodynamic quantities. The purpose of this Perspective is to shed new light on the latest developments and discuss future avenues.
Collapse
Affiliation(s)
- J-M Simon
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR-6303 CNRS - Université de Bourgogne Franche-Comté, F-21078 Dijon, France
| | - P Krüger
- Graduate School of Science and Engineering, Molecular Chirality Research Center, Chiba University, Chiba 263-8522, Japan
| | - S K Schnell
- Department of Materials Science and Engineering, NTNU - Norwegian University of Science and Technology, N-7491 Trondheim, Norway
| | - T J H Vlugt
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - S Kjelstrup
- Center of Excellence PoreLab, Department of Chemistry, Faculty of Natural Sciences, NTNU - Norwegian University of Science and Technology, Trondheim, Norway
| | - D Bedeaux
- Center of Excellence PoreLab, Department of Chemistry, Faculty of Natural Sciences, NTNU - Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
3
|
Mozafar O, Denniston C. Effects of structural inhomogeneity on equilibration processes in Langevin dynamics. Phys Rev E 2022; 105:064109. [PMID: 35854545 DOI: 10.1103/physreve.105.064109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Accepted: 05/16/2022] [Indexed: 11/06/2022]
Abstract
In recent decades, computer experiments have led to an accurate and fundamental understanding of atomic and molecular mechanisms in fluids, such as different kinds of relaxation processes toward steady physical states. In this paper, we investigate how exactly the configuration of initial states in a molecular-dynamics simulation can affect the rates of decay toward equilibrium for the widely known Langevin canonical ensemble. For this purpose, we derive an original expression relating the system relaxation time τ_{sys} and the radial distribution function g(r) in the near-zero and high-density limit. We found that, for an initial state which is slightly marginally inhomogeneous in the number density of atoms, the system relaxation time τ_{sys} is much longer than that for the homogeneous case and an increasing function of the Langevin coupling constant, γ. We also found, during structural equilibration, g(r) at large distances approaches 1 from above for the inhomogeneous case and from below for the macroscopically homogeneous one.
Collapse
Affiliation(s)
- Omid Mozafar
- Applied Mathematics Department, The University of Western Ontario, London, Ontario, Canada N6A 5B7
| | - Colin Denniston
- Physics and Astronomy Department, The University of Western Ontario, London, Ontario, Canada N6A 3K7
| |
Collapse
|
4
|
Venetsanos F, Anogiannakis SD, Theodorou DN. Mixing Thermodynamics and Flory–Huggins Interaction Parameter of Polyethylene Oxide/Polyethylene Oligomeric Blends from Kirkwood–Buff Theory and Molecular Simulations. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00642] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fotis Venetsanos
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| | - Stefanos D. Anogiannakis
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| | - Doros N. Theodorou
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| |
Collapse
|
5
|
Tripathy M, Bharadwaj S, van der Vegt NFA. Solvation shell thermodynamics of extended hydrophobic solutes in mixed solvents. J Chem Phys 2022; 156:164901. [DOI: 10.1063/5.0090646] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The ability of various cosolutes and cosolvents to enhance or quench solvent density fluctuations at solute–water interfaces has crucial implications on the conformational equilibrium of macromolecules such as polymers and proteins. Herein, we use an extended hydrophobic solute as a model system to study the effect of urea and methanol on the density fluctuations in the solute’s solvation shell and the resulting thermodynamics. On strengthening the solute–water/cosolute repulsive interaction, we observe distinct trends in the mutual affinities between various species in, and the thermodynamic properties of, the solvation shell. These trends strongly follow the respective trends in the preferential adsorption of urea and methanol: solute–water/cosolute repulsion strengthens, urea accumulation decreases, and methanol accumulation increases. Preferential accumulation of urea is found to quench the density fluctuations around the extended solute, leading to a decrease in the compressibility of the solvation shell. In contrast, methanol accumulation enhances the density fluctuations, leading to an increase in the compressibility. The mode of action of urea and methanol seems to be strongly coupled to their hydration behavior. The observations from this simple model is discussed in relation to urea driven swelling and methanol induced collapse of some well-known thermo-responsive polymers.
Collapse
Affiliation(s)
- Madhusmita Tripathy
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Swaminath Bharadwaj
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Nico F. A. van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| |
Collapse
|
6
|
Banerjee A, Sevilla M, Rudzinski JF, Cortes-Huerto R. Finite-size scaling and thermodynamics of model supercooled liquids: long-range concentration fluctuations and the role of attractive interactions. SOFT MATTER 2022; 18:2373-2382. [PMID: 35258066 DOI: 10.1039/d2sm00089j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
We compute partial structure factors, Kirkwood-Buff integrals (KBIs) and chemical potentials of model supercooled liquids with and without attractive interactions. We aim at investigating whether relatively small differences in the tail of the radial distribution functions result in contrasting thermodynamic properties. Our results suggest that the attractive potential favours the nucleation of long-range structures. Indeed, upon decreasing temperature, Bathia-Thornton structure factors display anomalous behaviour in the k→0 limit. KBIs extrapolated to the thermodynamic limit confirm this picture, and excess coordination numbers identify the anomaly with long-range concentration fluctuations. By contrast, the purely repulsive system remains perfectly miscible for the same temperature interval and only reveals qualitatively similar concentration fluctuations in the crystalline state. Furthermore, differences in both isothermal compressibilities and chemical potentials show that thermodynamics is not entirely governed by the short-range repulsive part of the interaction potential, emphasising the nonperturbative role of attractive interactions. Finally, at higher density, where both systems display nearly identical dynamical properties and repulsive interactions become dominant, the anomaly disappears, and both systems also exhibit similar thermodynamic properties.
Collapse
Affiliation(s)
- Atreyee Banerjee
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| | - Mauricio Sevilla
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| | - Joseph F Rudzinski
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| | | |
Collapse
|
7
|
Sevilla M, Cortes-Huerto R. Connecting density fluctuations and Kirkwood–Buff integrals for finite-size systems. J Chem Phys 2022; 156:044502. [DOI: 10.1063/5.0076744] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Affiliation(s)
- Mauricio Sevilla
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | | |
Collapse
|
8
|
Bråten V, Bedeaux D, Wilhelmsen Ø, Schnell SK. Small size effects in open and closed systems: What can we learn from ideal gases about systems with interacting particles? J Chem Phys 2021; 155:244504. [PMID: 34972373 DOI: 10.1063/5.0076684] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Small systems have higher surface area-to-volume ratios than macroscopic systems. The thermodynamics of small systems therefore deviates from the description of classical thermodynamics. One consequence of this is that properties of small systems can be dependent on the system's ensemble. By comparing the properties in grand canonical (open) and canonical (closed) systems, we investigate how a small number of particles can induce an ensemble dependence. Emphasis is placed on the insight that can be gained by investigating ideal gases. The ensemble equivalence of small ideal gas systems is investigated by deriving the properties analytically, while the ensemble equivalence of small systems with particles interacting via the Lennard-Jones or the Weeks-Chandler-Andersen potential is investigated through Monte Carlo simulations. For all the investigated small systems, we find clear differences between the properties in open and closed systems. For systems with interacting particles, the difference between the pressure contribution to the internal energy, and the difference between the chemical potential contribution to the internal energy, are both increasing with the number density. The difference in chemical potential is, with the exception of the density dependence, qualitatively described by the analytic formula derived for an ideal gas system. The difference in pressure, however, is not captured by the ideal gas model. For the difference between the properties in the open and closed systems, the response of increasing the particles' excluded volume is similar to the response of increasing the repulsive forces on the system walls. This indicates that the magnitude of the difference between the properties in open and closed systems is related to the restricted movement of the particles in the system. The work presented in this paper gives insight into the mechanisms behind ensemble in-equivalence in small systems, and illustrates how a simple statistical mechanical model, such as the ideal gas, can be a useful tool in these investigations.
Collapse
Affiliation(s)
- Vilde Bråten
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| | - Dick Bedeaux
- PoreLab, Department of Chemistry, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| | - Øivind Wilhelmsen
- PoreLab, Department of Chemistry, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| | - Sondre Kvalvåg Schnell
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| |
Collapse
|
9
|
Luo S, Thachuk M. Conservative Potentials for a Lattice-Mapped Coarse-Grained Scheme. J Phys Chem A 2021; 125:6486-6497. [PMID: 34264666 DOI: 10.1021/acs.jpca.1c02000] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The conservative potential, arising from a coarse-grain (CG) mapping scheme for nonbonded atomistic particles, is studied. This is a bottom-up approach from first-principles that maps atomistic particles to fluid element-like subcells whose centers lie on a regular, cubic lattice. Unlike standard CG mapping schemes, the current one uses dynamic labeling which on-the-fly changes the CG labels of the particles. The subcells can also be different sizes and shapes, in principle. Equilibrium atomistic molecular dynamics trajectories for different Lennard-Jones fluids are calculated and converted to CG ones, from which CG probability distribution functions are calculated. Correlation studies show position and mass CG variables are uncoupled in a given subcell, as are different vector components of position. Furthermore, the strongest coupling occurs with neighboring cells in specific directions, and the resulting distribution is well described by a multivariate Gaussian. This implies the CG potential has a generalized quadratic form, whose derivative can be determined analytically. A microscopic rationalization is provided for the signs and relative magnitudes of different correlation coefficients, and in some cases, a connection is made with bulk properties of the fluid. We argue the generalized quadratic form should be robust to changes in the particulars of the CG scheme, as well as the nature of the atomistic intermolecular potential. Only a few potential parameters need to be calculated from the underlying atomistic system. This is significant because it indicates the transferability of this form to other, more complex systems. This transferability will be tested in future work, where mapping schemes with fuzzy boundaries will be considered.
Collapse
Affiliation(s)
- Siwei Luo
- Department of Chemistry, University of British Columbia,Vancouver V6T 1Z1, Canada
| | - Mark Thachuk
- Department of Chemistry, University of British Columbia,Vancouver V6T 1Z1, Canada
| |
Collapse
|
10
|
Baptista LA, Dutta RC, Sevilla M, Heidari M, Potestio R, Kremer K, Cortes-Huerto R. Density-functional-theory approach to the Hamiltonian adaptive resolution simulation method. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:184003. [PMID: 33690194 DOI: 10.1088/1361-648x/abed1d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2021] [Accepted: 03/09/2021] [Indexed: 06/12/2023]
Abstract
In the Hamiltonian adaptive resolution simulation method (H-AdResS) it is possible to simulate coexisting atomistic (AT) and ideal gas representations of a physical system that belong to different subdomains within the simulation box. The Hamiltonian includes a field that bridges both models by smoothly switching on (off) the intermolecular potential as particles enter (leave) the AT region. In practice, external one-body forces are calculated and applied to enforce a reference density throughout the simulation box, and the resulting external potential adds up to the Hamiltonian. This procedure suggests an apparent dependence of the final Hamiltonian on the system's thermodynamic state that challenges the method's statistical mechanics consistency. In this paper, we explicitly include an external potential that depends on the switching function. Hence, we build a grand canonical potential for this inhomogeneous system to find the equivalence between H-AdResS and density functional theory (DFT). We thus verify that the external potential inducing a constant density profile is equal to the system's excess chemical potential. Given DFT's one-to-one correspondence between external potential and equilibrium density, we find that a Hamiltonian description of the system is compatible with the numerical implementation based on enforcing the reference density across the simulation box. In the second part of the manuscript, we focus on assessing our approach's convergence and computing efficiency concerning various model parameters, including sample size and solute concentrations. To this aim, we compute the excess chemical potential of water, aqueous urea solutions and Lennard-Jones (LJ) mixtures. The results' convergence and accuracy are convincing in all cases, thus emphasising the method's robustness and capabilities.
Collapse
Affiliation(s)
- L A Baptista
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - R C Dutta
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - M Sevilla
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - M Heidari
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - R Potestio
- Physics Department, University of Trento, via Sommarive, 14 I-38123 Trento, Italy
- INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, I-38123 Trento, Italy
| | - K Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - R Cortes-Huerto
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| |
Collapse
|
11
|
Bråten V, Wilhelmsen Ø, Schnell SK. Chemical Potential Differences in the Macroscopic Limit from Fluctuations in Small Systems. J Chem Inf Model 2021; 61:840-855. [PMID: 33566592 PMCID: PMC8023585 DOI: 10.1021/acs.jcim.0c01367] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a new method for computing chemical potential differences of macroscopic systems by sampling fluctuations in small systems. The small system method, presented by Schnell et al. [Schnell et al., J. Phys. Chem. B, 2011, 115, 10911], is used to create small embedded systems from molecular dynamics simulations, in which fluctuations of the number of particles are sampled. The sampled fluctuations represent the Boltzmann distributed probability of the number of particles. The overlapping region of two such distributions, sampled from two different systems, is used to compute their chemical potential difference. Since the thermodynamics of small systems is known to deviate from the classical thermodynamic description, the particle distributions will deviate from the macroscopic behavior as well. We show how this can be utilized to calculate the size dependence of chemical potential differences and eventually extract the chemical potential difference in the thermodynamic limit. The macroscopic chemical potential difference is determined with a relative error of 3% in systems containing particles that interact through the truncated and shifted Lennard-Jones potential. In addition to computing chemical potential differences in the macroscopic limit directly from molecular dynamics simulation, the new method provides insights into the size dependency that is introduced to intensive properties in small systems.
Collapse
Affiliation(s)
- Vilde Bråten
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| | - Øivind Wilhelmsen
- SINTEF Energy Research, Trondheim NO-7465, Norway.,Department of Energy and Process Engineering, Norwegian University of Science and Technology, Trondheim NO-7491, Norway
| | - Sondre Kvalvåg Schnell
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, Trondheim N-7491, Norway
| |
Collapse
|
12
|
Tripathy M, Bharadwaj S, B. SJ, van der Vegt NFA. Characterizing Polymer Hydration Shell Compressibilities with the Small-System Method. NANOMATERIALS 2020; 10:nano10081460. [PMID: 32722500 PMCID: PMC7466400 DOI: 10.3390/nano10081460] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 07/21/2020] [Accepted: 07/22/2020] [Indexed: 12/03/2022]
Abstract
The small-system method (SSM) exploits the unique feature of finite-sized open systems, whose thermodynamic quantities scale with the inverse system size. This scaling enables the calculation of properties in the thermodynamic limit of macroscopic systems based on computer simulations of finite-sized systems. We herein extend the SSM to characterize the hydration shell compressibility of a generic hydrophobic polymer in water. By systematically increasing the strength of polymer-water repulsion, we find that the excess inverse thermodynamic correction factor (Δ1/Γs∞) and compressibility (Δχs) of the first hydration shell change sign from negative to positive. This occurs with a concurrent decrease in water hydrogen bonding and local tetrahedral order of the hydration shell water. The crossover lengthscale corresponds to an effective polymer bead diameter of 0.7 nm and is consistent with previous works on hydration of small and large hydrophobic solutes. The crossover lengthscale in polymer hydration shell compressibility, herein identified with the SSM approach, relates to hydrophobic interactions and macromolecular conformational equilibria in aqueous solution. The SSM approach may further be applied to study thermodynamic properties of polymer solvation shells in mixed solvents.
Collapse
|
13
|
Heidari M, Kremer K, Golestanian R, Potestio R, Cortes-Huerto R. Open-boundary Hamiltonian adaptive resolution. From grand canonical to non-equilibrium molecular dynamics simulations. J Chem Phys 2020; 152:194104. [DOI: 10.1063/1.5143268] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Maziar Heidari
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Ramin Golestanian
- Max Planck Institute for Dynamics and Self-Organization (MPIDS), 37077 Göttingen, Germany
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3PU, United Kingdom
| | - Raffaello Potestio
- Physics Department, University of Trento, Via Sommarive, 14, I-38123 Trento, Italy
- INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, I-38123 Trento, Italy
| | | |
Collapse
|
14
|
Petris PC, Anogiannakis SD, Tzounis PN, Theodorou DN. Thermodynamic Analysis of n-Hexane-Ethanol Binary Mixtures Using the Kirkwood-Buff Theory. J Phys Chem B 2019; 123:247-257. [PMID: 30516991 DOI: 10.1021/acs.jpcb.8b10425] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A complete thermodynamic analysis of mixtures consisting of molecules with complex chemical constitution can be rather demanding. The Kirkwood-Buff theory of solutions allows the estimation of thermodynamic properties, which cannot be directly extracted from atomistic simulations, such as the Gibbs energy of mixing (Δmix G). In this work, we perform molecular dynamics simulations of n-hexane-ethanol binary mixtures in the liquid state under two temperature-pressure conditions and at various mole fractions. On the basis of the recently published methodology of Galata [ Fluid Phase Equilib. 2018 , 470 , 25 - 37 ] , we first calculate the Kirkwood-Buff integrals in the isothermal-isobaric ( NpT) ensemble, identifying how system size affects their estimation. We then extract the activity coefficients, excess Gibbs energy, excess enthalpy, and excess entropy for the n-hexane-ethanol binary mixtures we simulate. We employ two approaches for quantifying composition fluctuations: one based on counting molecular centers of mass and a second one based on counting molecular segments. Results from the two approaches are practically indistinguishable. We compare our results against predictions of vapor-liquid equilibria obtained in a previous simulation work using the same force field, as well as with experimental data, and find very good agreement. In addition, we develop a simple methodology to identify the hydrogen bonds between ethanol molecules and analyze their effects on mixing properties.
Collapse
Affiliation(s)
- Panagiotis C Petris
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| | - Stefanos D Anogiannakis
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| | | | - Doros N Theodorou
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| |
Collapse
|