The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed “fluctuation matching” is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.
more »
« less
Exploring the relationship between softness and excess entropy in glass-forming systems
We explore the relationship between a machine-learned structural quantity (softness) and excess entropy in simulations of supercooled liquids. Excess entropy is known to scale well the dynamical properties of liquids, but this quasi-universal scaling is known to breakdown in supercooled and glassy regimes. Using numerical simulations, we test whether a local form of the excess entropy can lead to predictions similar to those made by softness, such as the strong correlation with particles’ tendency to rearrange. In addition, we explore leveraging softness to compute excess entropy in the traditional fashion over softness groupings. Our results show that the excess entropy computed over softness-binned groupings is correlated with activation barriers to rearrangement.
more »
« less
- Award ID(s):
- 2309043
- PAR ID:
- 10500942
- Publisher / Repository:
- AIP
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 20
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Coarse-grained (CG) modeling has gained significant attention in recent years due to its wide applicability in enhancing the spatiotemporal scales of molecular simulations. While CG simulations, often performed with Hamiltonian mechanics, faithfully recapitulate structural correlations at equilibrium, they lead to ambiguously accelerated dynamics. In Paper I [J. Jin, K. S. Schweizer, and G. A. Voth, J. Chem. Phys. 158(3), 034103 (2023)], we proposed the excess entropy scaling relationship to understand the CG dynamics. Then, in Paper II [J. Jin, K. S. Schweizer, and G. A. Voth, J. Chem. Phys. 158(3), 034104 (2023)], we developed a theory to map the CG system into a dynamically consistent hard sphere system to analytically derive an expression for fast CG dynamics. However, many chemical and physical systems do not exhibit hard sphere-like behavior, limiting the extensibility of the developed theory. In this paper, we aim to generalize the theory to the non-hard sphere system based on the Weeks–Chandler–Andersen perturbation theory. Since non-hard sphere-like CG interactions affect the excess entropy term as it deviates from the hard sphere description, we explicitly account for the extra entropy to correct the non-hard sphere nature of the system. This approach is demonstrated for two different types of interactions seen in liquids, and we further provide a generalized description for any CG models using the generalized Gaussian CG models using Gaussian basis sets. Altogether, this work allows for extending the range and applicability of the hard sphere CG dynamics theory to a myriad of CG liquids.more » « less
-
Abstract Kauzmann paradox (KP) suggests that deeply supercooled liquids can have a lower entropy than the corresponding crystalline solids. While this entropy catastrophe has been thoroughly studied via equilibrium thermodynamics, the solidification process occurs far‐from‐equilibrium. By analyzing this process experimentally and theoretically, we show that surface chemical speciation (oxidation‐driven generation and self‐organization of different species of the alloy components) in core‐shell particles (CSPs) can perturb the entropy production to an extent that a continuum equilibrium phase transition is not possible. Speciation of the surface causes divergence of associated stress vectors that generate nonequilibrium fluxes and frustrates homogeneous nucleation hence deep undercooling. The asymmetry of the speciation‐derived surface tensor skews the minimum entropy production criterion. We analyze a set of nonequilibrium models, one showing and one averting the entropy catastrophe. Applying thermodynamic speed limits to these models, we show that the KP takes another form. Deviations from the speed limit diverge the configurational entropy of the glass, but adding an interfacial state avoids the entropy catastrophe with significantly large supercooling.more » « less
-
The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid–liquid critical point (LLCP) in the supercooled regime, at (Pc = 150 MPa, Tc = 190 K, and ρc = 1.04 g/cm3) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water and the corresponding location of the LLCP are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian, which allows us to calculate the configurational entropy of the system, Sconf. The Sconf of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam–Gibbs relation, which relates the diffusion coefficient D with Sconf, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high frequencies lead to unphysical (negative) entropy for q-TIP4P/F water when using classical statistical mechanics (yet, the PEL formalism can be applied successfully).more » « less
-
A fundamental and much-debated issue in glass science is the existence and nature of liquid–liquid transitions in glass-forming liquids. Here, we report the existence of a novel reentrant structural transition in a S-rich arsenic sulfide liquid of composition As 2.5 S 97.5 . The nature of this transition and its effect on viscosity are investigated in situ using a combination of differential scanning calorimetry and simultaneous Raman spectroscopic and rheometric measurements. The results indicate that, upon heating significantly above its glass transition temperature (261 K), the constituent [Formula: see text] sulfur chains in the structure of the supercooled liquid first undergo a [Formula: see text] chain-to-ring conversion near ∼383 K, which is exothermic in nature. Further heating above 393 K alters the equilibrium to shift in the opposite direction toward an endothermic ring-to-chain conversion characteristic of the well-known λ-transition in pure sulfur liquid. This behavior is attributed to the competing effects of enthalpy of mixing and conformational entropy of ring and chain elements in the liquid. The existence of reentrant structural transitions in glass-forming liquids could provide important insights into the thermodynamics of liquid–liquid transitions and may have important consequences for harnessing novel functionalities of derived glasses.more » « less
An official website of the United States government

