 Award ID(s):
 1905632
 NSFPAR ID:
 10385647
 Date Published:
 Journal Name:
 Soft Matter
 Volume:
 18
 Issue:
 41
 ISSN:
 1744683X
 Page Range / eLocation ID:
 7998 to 8007
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

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 LennardJones (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 unitedatom (UA) and allatom (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

The entropy change that occurs upon mixing two fluids has remained an intriguing topic since the dawn of statistical mechanics. In this work, we generalize the grandisobaric ensemble to mixtures and develop a Monte Carlo algorithm for the rapid determination of entropy in these systems. A key advantage of adiabatic ensembles is the direct connection they provide with entropy. Here, we show how the entropy of a binary mixture A–B can be readily obtained in the adiabatic grandisobaric ( μ_{A}, μ_{B}, P, R) ensemble, in which μ_{A}and μ_{B}denote the chemical potential of components A and B, respectively, P is the pressure, and R is the heat (Ray) function, that corresponds to the total energy of the system. This, in turn, allows for the evaluation of the entropy of mixing and the Gibbs free energy of mixing. We also demonstrate that our approach performs very well both on systems modeled with simple potentials and with complex manybody force fields. Finally, this approach provides a direct route to the determination of the thermodynamic properties of mixing and allows for the efficient detection of departures from ideal behavior in mixtures.

Coarsegrained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipidprotein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and WimleyWhite peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the allatom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membraneassociated peptides LS3 and WALP23. CgProt will find application in studies of lipidprotein interactions and the conformational properties of diverse membrane protein systems.more » « less

Abstract We investigate the rigidity of global minimizers u ≥ 0 u\ge 0 of the AltPhillips functional involving negative power potentials ∫ Ω ( ∣ ∇ u ∣ 2 + u − γ χ { u > 0 } ) d x , γ ∈ ( 0 , 2 ) , \mathop{\int }\limits_{\Omega }( \nabla u{ }^{2}+{u}^{\gamma }{\chi }_{\left\{u\gt 0\right\}}){\rm{d}}x,\hspace{1.0em}\gamma \in \left(0,2), when the exponent γ \gamma is close to the extremes of the admissible values. In particular, we show that global minimizers in R n {{\mathbb{R}}}^{n} are onedimensional if γ \gamma is close to 2 and n ≤ 7 n\le 7 , or if γ \gamma is close to 0 and n ≤ 4 n\le 4 .more » « less

The potential energy curves (PECs) for the homonuclear He–He, Ar–Ar, Cu–Cu, and Si–Si dimers, as well as heteronuclear Cu–He, Cu–Ar, Cu–Xe, Si–He, Si–Ar, and Si–Xe dimers, are obtained in quantum Monte Carlo (QMC) calculations. It is shown that the QMC method provides the PECs with an accuracy comparable with that of the stateoftheart coupled cluster singles and doubles with perturbative triples corrections [CCSD(T)] calculations. The QMC data are approximated by the Morse long range (MLR) and (126) LennardJones (LJ) potentials. The MLR and LJ potentials are used to calculate the deflection angles in binary collisions of corresponding atom pairs and transport coefficients of Cu and Si vapors and their mixtures with He, Ar, and Xe gases in the range of temperature from 100 K to 10 000 K. It is shown that the use of the LJ potentials introduces significant errors in the transport coefficients of hightemperature vapors and gas mixtures. The mixtures with heavy noble gases demonstrate anomalous behavior when the viscosity and thermal conductivity can be larger than that of the corresponding pure substances. In the mixtures with helium, the thermal diffusion factor is found to be unusually large. The calculated viscosity and diffusivity are used to determine parameters of the variable hard sphere and variable soft sphere molecular models as well as parameters of the powerlaw approximations for the transport coefficients. The results obtained in the present work include all information required for kinetic or continuum simulations of dilute Cu and Si vapors and their mixtures with He, Ar, and Xe gases.more » « less