An official website of the United States government
Predicting the mixing free energy of mixing for binary mixtures using simulations is challenging. We present a novel molecular dynamics (MD) simulation method to extract the chemical potential μ ( X ) for mixtures of species A and B. Each molecule of species A and B is placed in equal and opposite harmonic potentials ±(1/2) U ex ( x ) centered at the middle of the simulation box, resulting in a nonuniform mole fraction profile X ( z ) in which A is concentrated at the center, and B at the periphery. Combining these, we obtain U ex ( X ), the exchange chemical potential required to induce a given deviation of the mole fraction from its average. Simulation results for U ex ( X ) can be fitted to simple free energy models to extract the interaction parameter χ for binary mixtures. To illustrate our method, we investigate benzene–pyridine mixtures, which provide a good example of regular solution behavior, using both TraPPE united-atom and OPLS all-atom potentials, both of which have been validated for pure fluid properties. χ values obtained with the new method are consistent with values from other recent simulation methods. However, the TraPPE-UA results differ substantially from the χ obtained from VLE experimental data, while the OPLS-AA results are in reasonable agreement with experiment, highlighting the importance of accurate potentials in correctly representing mixture behavior.
more »
« less
Predicting mixing free energy using mutual ghosting
The excess free energy of mixing Δ G ex governs the phase behavior of mixtures and controls material properties. It is challenging, however, to measure Δ G ex in simulations. Previously, we developed a method that combines molecular dynamics (MD) simulations with thermodynamic integration along the path of transformation of chains to predict the Flory Huggins interaction parameter χ for polymer mixtures and block copolymers. However, this method is best applied when the constituent molecules of the blends are structurally related. To overcome this limitation, we have developed a new method to predict Δ G ex for mixtures. We perform simulations to induce phase separation within a mixture by gradually weakening the interaction between different species. To compute Δ G ex we measure the thermodynamic work required to modify the interactions and the interfacial energy between the separated phases. We validate our method by applying it first to equimolar mixtures of labeled and unlabeled Lennard-Jones (LJ) beads, and labeled and unlabeled benzene, which results in good agreement with ideal solution theory. Then we compute the excess free energy of mixing for equimolar mixtures of benzene and pyridine, using both united-atom (UA) and all-atom (AA) potentials. Our results using UA potentials predict a value for Δ G ex about four times the experimental value, whereas using AA potentials gives results consistent with experiment, highlighting the need for good potentials to faithfully represent mixture behavior.
more »
« less
- PAR ID:
- 10380436
- Date Published:
- Journal Name:
- Molecular Systems Design & Engineering
- Volume:
- 7
- Issue:
- 11
- ISSN:
- 2058-9689
- Page Range / eLocation ID:
- 1529 to 1537
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract By eliminating unnecessary details, coarse-grained (CG) models provide the necessary efficiency for simulating scales that are inaccessible to higher resolution models. However, because they average over atomic details, the effective potentials governing CG degrees of freedom necessarily incorporate significant entropic contributions, which limit their transferability and complicate the treatment of thermodynamic properties. This work employs a dual-potential approach to consider the energetic and entropic contributions to effective interaction potentials for CG models. Specifically, we consider one- and three-site CG models for ortho-terphenyl (OTP) both above and below its glass transition. We employ the multiscale coarse-graining (MS-CG) variational principle to determine interaction potentials that accurately reproduce the structural properties of an all-atom (AA) model for OTP at each state point. We employ an energy-matching variational principle to determine an energy operator that accurately reproduces the intra- and inter-molecular energy of the AA model. While the MS-CG pair potentials are almost purely repulsive, the corresponding pair energy functions feature a pronounced minima that corresponds to contacting benzene rings. These energetic functions then determine an estimate for the entropic component of the MS-CG interaction potentials. These entropic functions accurately predict the MS-CG pair potentials across a wide range of liquid state points at constant density. Moreover, the entropic functions also predict pair potentials that quite accurately model the AA pair structure below the glass transition. Thus, the dual-potential approach appears a promising approach for modeling AA energetics, as well as for predicting the temperature-dependence of CG effective potentials.more » « less
-
To achieve high power density and thermodynamic cycle efficiency, the working pressures of liquid-propellant rocket engines, diesel engines, and gas turbines (based on deflagration or detonation) are continuously increasing, which could reach or go beyond the thermodynamic critical pressure of the liquid propellant. For this reason, the studies of trans- and super-critical injection are getting more and more attention. However, the simulation of transcritical phase change is still a challenging topic. The phase boundary, especially near the mixture critical point, needs to be accurately determined to investigate the multicomponent effects on transcritical injection and atomization. This work used our previously developed thermodynamic model based on the vapor-liquid equilibrium (VLE) theory, which can predict the phase separation near the mixture critical point. An \textit{in situ} adaptive tabulation (ISAT) method was developed to accelerate the computationally expensive multicomponent VLE computation such that it can be cheap enough for CFD. The new thermodynamic model was integrated into OpenFOAM to develop a VLE-based CFD solver. In this work, shock-droplet interaction and two-phase mixing simulations are conducted using our new VLE-based CFD solver. The shock-droplet interaction simulation results capture the thermodynamic condition of the surface entering the supercritical state after shock passes through. The atomization of droplets could be triggered by vorticity formed at the droplets' surface. 2D temporal mixing layer simulations show the evolution of the transcritical mixing layer and capture the phase split effect at the mixing layer.more » « less
-
Abstract. Secondary organic aerosol (SOA) constitutes a largefraction of atmospheric aerosol. To assess its impacts on climate and airpollution, knowledge of the number of phases in internal mixtures ofdifferent SOA types is required. Atmospheric models often assume thatdifferent SOA types form a single phase when mixed. Here, we present visualobservations of the number of phases formed after mixing differentanthropogenic and biogenic SOA types. Mixing SOA types generated inenvironmental chambers with oxygen-to-carbon (O/C) ratios between 0.34 and 1.05, we found 6 out of 15 mixtures of two SOA types to result in two phase particles. We demonstrate that the number of phases depends on thedifference in the average O/C ratio between the two SOA types (Δ(O/C)). Using a threshold Δ(O/C) of 0.47, we can predict the phasebehavior of over 90 % of our mixtures, with one- and two-phase particlespredicted for Δ(O/C)<0.47 and Δ(O/C)≥0.47,respectively. This threshold ΔO/C value provides a simple parameterto predict whether mixtures of fresh and aged SOA form one- or two-phase particles in the atmosphere. In addition, we show that phase-separated SOAparticles form when mixtures of volatile organic compounds emitted from realtrees are oxidized.more » « less
-
While the excess chemical potential is the key quantity in determining phase diagrams, its direct computation for high-density liquids of long polymer chains has posed a significant challenge. Computationally, the excess chemical potential is calculated using the Widom insertion method, which involves monitoring the change in internal energy as one incrementally introduces individual molecules into the liquid. However, when dealing with dense polymer liquids, inserting long chains requires generating trial configurations with a bias that favors those at low energy on a unit-by-unit basis: a procedure that becomes more challenging as the number of units increases. Thus, calculating the excess chemical potential of dense polymer liquids using this method becomes computationally intractable as the chain length exceeds N ≥ 30. Here, we adopt a coarsegrained model derived from the integral equation theory for which inserting long polymer chains becomes feasible. The integral equation theory of coarse graining (IECG) represents a polymer as a sphere or a collection of blobs interacting through a soft potential. We employ the IECG approach to compute the excess chemical potential using Widom’s method for polymer chains of increasing lengths, extending up to N = 720 monomers, and at densities reaching up to ρ = 0.767 g/cm3. From a fundamental perspective, we demonstrate that the excess chemical potentials remain nearly constant across various levels of coarse graining, offering valuable insights into the consistency of this type of procedure. Ultimately, we argue that current Monte Carlo algorithms, originally designed for atomistic simulations, such as configurational bias Monte Carlo (CBMC) methods, can significantly benefit from the integration of the IECG approach, thereby enhancing their performance in the study of phase diagrams of polymer liquids.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D2ME00109H
