Liquid–liquid phase transition in an ionic model of silica

Crystallization and the liquid-liquid critical point in nonbonded modified-WAC models

For decades, it has been known that Liquid-Liquid Critical Points (LLCPs) can exist in one-component liquids, yet a comprehensive understanding of the conditions under which they arise remains elusive. To better comprehend the possible interplay between the LLCP and the crystalline phase, we conduct molecular dynamics simulations using the nonbonded family of modified-WAC (mWAC) models, which are known to exhibit a LLCP for certain parameter values. By comparing different versions of the mWAC model-those featuring a LLCP and those lacking one-we identify several key differences between the models relating to crystallization. Those models that do have a LLCP are found to have multiple stable crystalline phases, one of them being a solid-state ionic conductor similar to superionic ice. Moreover, we find that for models that do not have a LLCP, the liquid becomes a glass at a larger range of temperatures, possibly preventing the occurrence of a LLCP. Further studies are required to determine if these results are general or model-specific.

Ionic liquids as the effective technology for enhancing transdermal drug delivery: Design principles, roles, mechanisms, and future challenges

Ionic liquids (ILs) have been proven to be an effective technology for enhancing drug transdermal absorption. However, due to the unique structural components of ILs, the design of efficient ILs and elucidation of action mechanisms remain to be explored. In this review, basic design principles of ideal ILs for transdermal drug delivery system (TDDS) are discussed considering melting point, skin permeability, and toxicity, which depend on the molar ratios, types, functional groups of ions and inter-ionic interactions. Secondly, the contributions of ILs to the development of TDDS through different roles are described: as novel skin penetration enhancers for enhancing transdermal absorption of drugs; as novel solvents for improving the solubility of drugs in carriers; as novel active pharmaceutical ingredients (API-ILs) for regulating skin permeability, solubility, release, and pharmacokinetic behaviors of drugs; and as novel polymers for the development of smart medical materials. Moreover, diverse action mechanisms, mainly including the interactions among ILs, drugs, polymers, and skin components, are summarized. Finally, future challenges related to ILs are discussed, including underlying quantitative structure-activity relationships, complex interaction forces between anions, drugs, polymers and skin microenvironment, long-term stability, and in vivo safety issues. In summary, this article will promote the development of TDDS based on ILs.

Kinetic reconstruction of free energies as a function of multiple order parameters

A vast array of phenomena, ranging from chemical reactions to phase transformations, are analyzed in terms of a free energy surface defined with respect to a single or multiple order parameters. Enhanced sampling methods are typically used, especially in the presence of large free energy barriers, to estimate free energies using biasing protocols and sampling of transition paths. Kinetic reconstructions of free energy barriers of intermediate height have been performed, with respect to a single order parameter, employing the steady state properties of unconstrained simulation trajectories when barrier crossing is achievable with reasonable computational effort. Considering such cases, we describe a method to estimate free energy surfaces with respect to multiple order parameters from a steady state ensemble of trajectories. The approach applies to cases where the transition rates between pairs of order parameter values considered is not affected by the presence of an absorbing boundary, whereas the macroscopic fluxes and sampling probabilities are. We demonstrate the applicability of our prescription on different test cases of random walkers executing Brownian motion in order parameter space with an underlying (free) energy landscape and discuss strategies to improve numerical estimates of the fluxes and sampling. We next use this approach to reconstruct the free energy surface for supercooled liquid silicon with respect to the degree of crystallinity and density, from unconstrained molecular dynamics simulations, and obtain results quantitatively consistent with earlier results from umbrella sampling.

Anomalous Vapor and Ice Nucleation in Water at Negative Pressures: A Classical Density Functional Theory Study

In contrast to the abundance of work on the anomalous behavior of water, the relationship between the water's thermodynamic anomalies and kinetics of phase transition from metastable water is relatively unexplored. In this work, we have employed classical density functional theory to provide a unified and coherent picture of nucleation (both vapor and ice) from metastable water at negative pressure conditions. Our results suggest a peculiar nonmonotonic temperature dependence of vapor-liquid surface tension at temperatures where vapor-liquid coexistence is metastable with respect to the ice phase. The vapor nucleation barrier on isochoric cooling also shows a nonmonotonic temperature dependence. We further report that, for low density isochores, the temperature of the minimum vapor nucleation barrier (TΔΩv/min*) does not coincide with the temperature of maximum density (TMD) where metastability is maximum. The difference between the TΔΩv/min* and the TMD, however, decreases with increasing the density of the isochore. The vapor nucleation barrier along isobars shows an interesting crossover behavior in the vicinity of the Widom line on lowering the temperature. Our results on the ice nucleation suggest an anomalous retracing behavior of the nucleation barrier along isotherms at negative pressures and theoretically validate the recent findings that the reentrant ice(Ih)-liquid coexistence line can induce a drastic change in the kinetics of ice nucleation. Thus, this study establishes a direct connection between the metastable water's thermodynamic anomalies and the (vapor and ice) nucleation kinetics. In addition, this study provides deeper insights into the origin of the isothermal compressibility maximum on isochoric cooling.

Liquid-liquid phase transition in deeply supercooled Stillinger-Weber silicon

The existence of a phase transition between two distinct liquid phases in single-component network-forming liquids (e.g. water, silica, silicon) has elicited considerable scientific interest. The challenge, both for experiments and simulations, is that the liquid-liquid phase transition (LLPT) occurs under deeply supercooled conditions, where crystallization occurs very rapidly. Thus, early evidence from numerical equation of state studies was challenged with the argument that slow spontaneous crystallization had been misinterpreted as evidence of a second liquid state. Rigorous free-energy calculations have subsequently confirmed the existence of a LLPT in some models of water, and exciting new experimental evidence has since supported these computational results. Similar results have so far not been found for silicon. Here, we present results from free-energy calculations performed for silicon modeled with the classical, empirical Stillinger-Weber-potential. Through a careful study employing state-of-the-art constrained simulation protocols and numerous checks for thermodynamic consistency, we find that there are two distinct metastable liquid states and a phase transition. Our results resolve a long-standing debate concerning the existence of a liquid-liquid transition in supercooled liquid silicon and address key questions regarding the nature of the phase transition and the associated critical point.

Modeling fluid polyamorphism through a maximum-valence approach

We suggest a simple model to describe polyamorphism in single-component fluids using a maximum-valence approach. The model contains three types of interactions: (i) Atoms attract each other by van der Waals forces that generate a liquid-gas transition at low pressures, (ii) atoms may form covalent bonds that induce association, and (iii) atoms with maximal valence attract or repel each other stronger than other atoms, thus generating liquid-liquid separation. As an example, we qualitatively compare this model with the behavior of liquid sulfur and show that condition (iii) generates a liquid-liquid phase transition in addition to the liquid-gas phase transition.

Signatures of sluggish dynamics and local structural ordering during ice nucleation

We investigate the microscopic pathway of spontaneous crystallization in the ST2 model of water under deeply supercooled conditions via unbiased classical molecular dynamics simulations. After quenching below the liquid–liquid critical point, the ST2 model spontaneously separates into low-density liquid (LDL) and high-density liquid phases, respectively. The LDL phase, which is characterized by lower molecular mobility and enhanced structural order, fosters the formation of a sub-critical ice nucleus that, after a stabilization time, develops into the critical nucleus and grows. Polymorphic selection coincides with the development of the sub-critical nucleus and favors the formation of cubic (Ic) over hexagonal (Ih) ice. We rationalize polymorphic selection in terms of geometric arguments based on differences in the symmetry of second neighbor shells of ice Ic and Ih, which are posited to favor formation of the former. The rapidly growing critical nucleus absorbs both Ic and Ih crystallites dispersed in the liquid phase, a crystal with stacking faults. Our results are consistent with, and expand upon, recent observations of non-classical nucleation pathways in several systems.

Avalanches, Clusters, and Structural Change in Cyclically Sheared Silica Glass

We investigate avalanches and clusters associated with plastic rearrangements and the nature of structural change in the prototypical strong glass, silica, computationally. We perform a detailed analysis of avalanches, and of spatially disconnected clusters that constitute them, for a wide range of system sizes. Although qualitative aspects of yielding in silica are similar to other glasses, the statistics of clusters exhibits significant differences, which we associate with differences in local structure. Across the yielding transition, anomalous structural change and densification, associated with a suppression of tetrahedral order, is observed to accompany strain localization.

The role of annealing in determining the yielding behavior of glasses under cyclic shear deformation

Yielding behavior in amorphous solids has been investigated in computer simulations using uniform and cyclic shear deformation. Recent results characterize yielding as a discontinuous transition, with the degree of annealing of glasses being a significant parameter. Under uniform shear, discontinuous changes in stresses at yielding occur in the high annealing regime, separated from the poor annealing regime in which yielding is gradual. In cyclic shear simulations, relatively poorly annealed glasses become progressively better annealed as the yielding point is approached, with a relatively modest but clear discontinuous change at yielding. To understand better the role of annealing on yielding characteristics, we perform athermal quasistatic cyclic shear simulations of glasses prepared with a wide range of annealing in two qualitatively different systems-a model of silica (a network glass) and an atomic binary mixture glass. Two strikingly different regimes of behavior emerge. Energies of poorly annealed samples evolve toward a unique threshold energy as the strain amplitude increases, before yielding takes place. Well-annealed samples, in contrast, show no significant energy change with strain amplitude until they yield, accompanied by discontinuous energy changes that increase with the degree of annealing. Significantly, the threshold energy for both systems corresponds to dynamical cross-over temperatures associated with changes in the character of the energy landscape sampled by glass-forming liquids.

Pulsed electric fields induce modulation of protein liquid-liquid phase separation

The time-resolved dynamic assembly and the structures of protein liquid dense clusters (LDCs) were analyzed under pulsed electric fields (EFs) applying complementary polarized and depolarized dynamic light scattering (DLS/DDLS), optical microscopy, and transmission electron microscopy (TEM). We discovered that pulsed EFs substantially affected overall morphologies and spatial distributions of protein LDCs and microcrystals, and affected the phase diagrams of LDC formation, including enabling protein solutions to overcome the diffusive flux energy barrier to phase separate. Data obtained from DLS/DDLS and TEM showed that LDCs appeared as precursors of protein crystal nuclei, followed by the formation of ordered structures within LDCs applying a pulsed EF. Experimental results of circular dichroism spectroscopy provided evidence that the protein secondary structure content is changing under EFs, which may consequently modulate protein-protein interactions, and the morphology, dimensions, and internal structure of LDCs. Data and results obtained unveil options to modulate the phase diagram of crystallization, and physical morphologies of protein LDCs and microcrystals by irradiating sample suspensions with pulsed EFs.

Effects of the bond polarity on the structural and dynamical properties of silica-like liquids

Silica is a network-forming liquid that shares many properties with water due to its tetrahedral structure. It undergoes a transition from a fragile to a strong liquid as the temperature is decreased, which is accompanied by a structural change to lower density and higher tetrahedral order. In order to disentangle the effects of Coulomb and van der Waals interactions on the structure and dynamics of liquid silica, we modify the bond polarity by changing the partial charges assigned to each atom. Using molecular dynamics simulations, we show that density, tetrahedral order, and structural relaxation times decrease when reducing bond polarity. Moreover, we find that the density maximum and the fragile-to-strong transition move to lower temperatures until they eventually vanish when the partial charges are decreased below approximately 75% of their regular value. Irrespective of whether strong or fragile behavior exists, structural relaxation is governed by hopping motion at sufficiently low temperatures. As long as there is a strong regime, the energy barrier associated with strong dynamics decreases with decreasing partial charges, but the dependence on the bond polarity differs from that of the activation energy in the Arrhenius regime at high temperatures. We show that the fragile-to-strong transition is associated with structural changes occurring between the first and second coordination shells that lead to a decrease in density and an increase in tetrahedral order. In particular, independent of the value of the partial charges, the distribution of the local structures is the same at this dynamic crossover, but we find no evidence that the effect occurs upon crossing the Widom line. In the fragile regime at intermediate temperatures, the relaxation times are well described by a previously proposed model which decomposes the apparent activation energy into a constant single-particle contribution and a temperature-dependent collective contribution. However, our results for silica-like melts do not obey several common relations of the model parameters reported for molecular glass formers.

Anomalous phase behavior of first-order fluid-liquid phase transition in phosphorus

Although the existence of liquid-liquid phase transition has become more and more convincing, whether it will terminate at a critical point and what is the order parameter are still open. To explore these questions, we revisit the fluid-liquid phase transition (FLPT) in phosphorus (P) and study its phase behavior by performing extensive first-principles molecular dynamics simulations. The FLPT observed in experiments is well reproduced, and a fluid-liquid critical point (FLCP) at T = 3000 ∼ 3500 K, P = 1.5-2.0 Kbar is found. With decreasing temperature from the FLCP along the transition line, the density difference (Δρ) between two coexisting phases first increases from zero and then anomalously decreases; however, the entropy difference (ΔS) continuously increases from zero. These features suggest that an order parameter containing contributions from both the density and the entropy is needed to describe the FLPT in P, and at least at low temperatures, the entropy, instead of the density, governs the FLPT.

Fluctuations near the liquid–liquid transition in a model of silica

Molecular dynamics simulations reveal anomalous small-angle scattering and liquid–liquid phase separation in an ionic model of silica.

Abnormal separation of the silicon–oxygen bond in the liquid layering transition of silicon dioxide in a nanoslit

Layering transition and separation of silicon and oxygen in liquid SiO₂ become obvious due to the strengthening of the nanoconfined effect.