Phase behavior of density-dependent pair potentials

Phase diagram is calculated by a recently proposed third-order thermodynamic perturbation theory (TPT) for fluid phase and a recently proposed first-order TPT for solid phases; the underlying interparticle potential consists of a hard sphere repulsion and a perturbation tail of an attractive inverse power law type or Yukawa type whose range varies with bulk densities. It is found that besides usual phase transitions associated with density-independent potentials, the density dependence of the perturbation tail evokes some additional novel phase transitions including isostructural solid-solid transition and liquid-liquid transition. Novel triple points are also exhibited which includes stable fluid (vapor or liquid)-face-centered cubic(fcc)-fcc and liquid-liquid-fcc, metastable liquid-body-centered cubic(bcc)-bcc. It also is found that the phase diagram sensitively depends on the density dependence and the concrete mathematical form of the underlying potentials. Some of the disclosed novel transitions has been observed experimentally in complex fluids and molecular liquids, while others still remain to be experimentally verified.     

From rolling ball to complete wetting: the dynamic tuning of liquids on nanostructured surfaces

In this work, for the first time, a dynamic electrical control of the wetting behavior of liquids on nanostructured surfaces, which spans the entire possible range from the superhydrophobic behavior to nearly complete wetting, has been demonstrated. Moreover, this kind of dynamic control was obtained at voltages as low as 22 V. We have demonstrated that the liquid droplet on a nanostructured surface exhibits sharp transitions between three possible wetting states as a function of applied voltage and liquid surface tension. We have examined experimentally and theoretically the nature of these transitions. The reported results provide novel methods of manipulating liquids at the microscale.     

Liquid-liquid phase transformations and the shape of the melting curve

The phase diagram of elemental liquids has been found to be surprisingly rich, including variations in the melting curve and transitions in the liquid phase. The effect of these transitions in the liquid state on the shape of the melting curve is analyzed. First-order phase transitions intersecting the melting curve imply piecewise continuous melting curves, with solid-solid transitions generating upward kinks or minima and liquid-liquid transitions generating downward kinks or maxima. For liquid-liquid phase transitions proposed for carbon, phosphorous selenium, and possibly nitrogen, we find that the melting curve exhibits a kink. Continuous transitions imply smooth extrema in the melting curve, the curvature of which is described by an exact thermodynamic relation. This expression indicates that a minimum in the melting curve requires the solid compressibility to be greater than that of the liquid, a very unusual situation. This relation is employed to predict the loci of smooth maxima at negative pressures for liquids with anomalous melting curves. The relation between the location of the melting curve maximum and the two-state model of continuous liquid-liquid transitions is discussed and illustrated by the case of tellurium.     

Water polyamorphism: reversibility and (dis)continuity

An understanding of water's anomalies is closely linked to an understanding of the phase diagram of water's metastable noncrystalline states. Despite the considerable effort, such an understanding has remained elusive and many puzzles regarding phase transitions in supercooled liquid water and their possible amorphous proxies at low temperatures remain. Here, decompression of very high density amorphous ice (VHDA) from 1.1 to 0.02 GPa at 140 K is studied by means of dilatometry and powder x-ray diffraction of quench-recovered states. It is shown that the three amorphous states of ice are reversibly connected to each other, i.e., LDA<-->e-HDA<-->VHDA. However, while the downstroke VHDA-->e-HDA transition takes place in the pressure range of 0.06 GPaLDA transition takes place quasi-discontinuously at p approximately 0.06 GPa. That is, two amorphous-amorphous transitions of a distinct nature are observed for the first time in a one-component system-a first-order-like transition (e-HDA-->LDA) and a transition which is not first-order like but possibly of higher order (VHDA-->e-HDA). VHDA and e-HDA are established as the most stable and limiting states in the course of the transition. We interpret this as evidence disfavoring the hypothesis of multiple first-order liquid-liquid transitions (and the option of a third critical point), but favoring a single first-order liquid-liquid transition (and the option of a second critical point).     

Structural origin of energetic heterogeneity in ionic liquids

How structural features observed computationally are connected to excitation-wavelength dependent kinetics observed experimentally remains an unanswered question in the field of ionic liquids. Using molecular dynamics simulation methods and approximated models that simplify the electrostatic interactions in ionic liquids, we discovered that on the timescale shorter than the relaxation time of the photochemical process, the energetic heterogeneity in terms of distribution of excitation energies is the consequence of the structure heterogeneity formed by local arrangement of ions around the solute probe.     

The putative liquid-liquid transition is a liquid-solid transition in atomistic models of water

We use numerical simulation to examine the possibility of a reversible liquid-liquid transition in supercooled water and related systems. In particular, for two atomistic models of water, we have computed free energies as functions of multiple order parameters, where one is density and another distinguishes crystal from liquid. For a range of temperatures and pressures, separate free energy basins for liquid and crystal are found, conditions of phase coexistence between these phases are demonstrated, and time scales for equilibration are determined. We find that at no range of temperatures and pressures is there more than a single liquid basin, even at conditions where amorphous behavior is unstable with respect to the crystal. We find a similar result for a related model of silicon. This result excludes the possibility of the proposed liquid-liquid critical point for the models we have studied. Further, we argue that behaviors others have attributed to a liquid-liquid transition in water and related systems are in fact reflections of transitions between liquid and crystal.     

On the Adam-Gibbs-Kirkpatrick-Thirumalai-Wolynes scenario for the viscosity increase in glasses

We reformulate the interpretation of the mean-field glass transition scenario for finite dimensional systems, proposed by Kirkpatrick, Thirumalai, and Wolynes (KTW). This allows us to establish clearly a temperature dependent length xi( *) above which the mean-field glass transition picture has to be modified. We argue in favor of the mosaic state introduced by KTW, which leads to the Adam-Gibbs relation between the viscosity and configurational entropy of glass forming liquids. Our argument is a mixture of thermodynamics and kinetics, partly inspired by the random energy model: small clusters of particles are thermodynamically frozen in low energy states, whereas large clusters are kinetically frozen by large activation energies. The relevant relaxation time is that of the smallest "liquid" clusters. Some physical consequences are discussed.     

Solid-liquid transitions of sodium chloride at high pressures

We investigate solid-liquid transitions in NaCl at high pressures using molecular dynamics simulations, including the melting curve and superheating/supercooling as well as solid-solid and liquid-liquid transitions. The first-order B1-B2 (NaCl-CsCl type) transition in solid is observed at high pressures besides continuous liquid structure transitions, which are largely analogous to the B1-B2 transition in solid. The equilibrium melting temperatures (T(m)) up to megabar pressure are obtained from the solid-liquid coexistence technique and the superheating-supercooling hysteresis method. Lindemann's vibrational and Born's mechanical instabilities are found upon melting. The Lindemann frequency is calculated from the vibrational density of states. The Lindemann parameter (fractional root-mean-squared displacement) increases with pressure and approaches a constant asymptotically, similar to the Lennard-Jones system. However, the Lindemann melting relation holds for both B1 and B2 phases to high accuracy as for the Lennard-Jonesium. The B1 and B2 NaCl solids can be superheated by 0.18T(m) and 0.24T(m), and the NaCl liquid, supercooled by 0.22T(m) and 0.32T(m), respectively, at heating or cooling rates of 1 K/s and 1 K/ps. The amount of maximum superheating or supercooling and its weak pressure dependence observed for NaCl are in accord with experiments on alkali halides and with simulations on the Lennard-Jones system and Al.     

Liquid-liquid phase transitions in supercooled water studied by computer simulations of various water models

Liquid-liquid and liquid-vapor coexistence regions of various water models were determined by Monte Carlo (MC) simulations of isotherms of density fluctuation-restricted systems and by Gibbs ensemble MC simulations. All studied water models show multiple liquid-liquid phase transitions in the supercooled region: we observe two transitions of the TIP4P, TIP5P, and SPCE models and three transitions of the ST2 model. The location of these phase transitions with respect to the liquid-vapor coexistence curve and the glass temperature is highly sensitive to the water model and its implementation. We suggest that the apparent thermodynamic singularity of real liquid water in the supercooled region at about 228 K is caused by an approach to the spinodal of the first (lowest density) liquid-liquid phase transition. The well-known density maximum of liquid water at 277 K is related to the second liquid-liquid phase transition, which is located at positive pressures with a critical point close to the maximum. A possible order parameter and the universality class of liquid-liquid phase transitions in one-component fluids are discussed.     

Pressure effects on the transitions between disordered phases in supercooled liquid silicon

We investigate the pressure effects on the transitions between the disordered phases in supercooled liquid silicon through Monte Carlo simulations and efficient methods to compute free energies. Our calculations, using an environment dependent interatomic potential for Si, indicate that at zero pressure the liquid-liquid phase transition, between the high density liquid and the low density liquid, occurs at a temperature 325K below melting. We found that the liquid-liquid transition temperature decreases with increasing pressure, following the liquid-solid coexistence curve. As pressure increases, the liquid-liquid coexistence curve approaches the region where the glass transition between the low density liquid and the low density amorphous takes place. Above 5 GPa, our calculations show that the liquid-liquid transition is suppressed by the glassy dynamics of the system. We also found that above 5 GPa, the glass transition temperature is lower than that at lower pressures, suggesting that under these conditions the glass transition occurs between the high density liquid and the high density amorphous.     

On the role of frustration on the glass transition and polyamorphism of mesoscopically heterophase liquids

The model of heterophase fluctuations is developed accounting frustration of the mesoscopic solidlike fluctuons. Within the framework of this model, the glass transition and polyamorphous transformations are considered. It is shown that the frustration increases the temperature range in which the heterophase liquid state exists. the upper and lower boundaries of this temperature range are determined. These boundaries separate different phase states-amorphous solid, heterophase liquid, and fluid phases. Polyamorphous liquid-liquid transitions in the liquid are investigated. Frustration can call forth continuous fluid-solid phase transformation avoiding the first- or second-order phase transition. Conditions under which the first-order phase transition fraction takes place are formulated. Two scenarios of the first-order liquid-liquid polyamorphous transformation are described. As an example the glacial phase formation and the first-order liquid-liquid phase transition in triphenyl phosphate are considered and discussed. Impact of frustration on the liquid crystallization and crystallinity of the glassy state is studied.     

Theoretical study of the absorption and nonradiative deactivation of 1-nitronaphthalene in the low-lying singlet and triplet excited states including methanol and ethanol solvent effects

The photophysics of the 1-nitronaphthalene molecular system, after the absorption transition to the first singlet excited state, is theoretically studied for investigating the ultrafast multiplicity change to the triplet manifold. The consecutive transient absorption spectra experimentally observed in this molecular system are also studied. To identify the electronic states involved in the nonradiative decay, the minimum energy path of the first singlet excited state is obtained using the complete active space self-consistent field∕∕configurational second-order perturbation approach. A near degeneracy region was found between the first singlet and the second triplet excited states with large spin-orbit coupling between them. The intersystem crossing rate was also evaluated. To support the proposed deactivation model the transient absorption spectra observed in the experiments were also considered. For this, computer simulations using sequential quantum mechanic-molecular mechanic methodology was used to consider the solvent effect in the ground and excited states for proper comparison with the experimental results. The absorption transitions from the second triplet excited state in the relaxed geometry permit to describe the transient absorption band experimentally observed around 200 fs after the absorption transition. This indicates that the T(2) electronic state is populated through the intersystem crossing presented here. The two transient absorption bands experimentally observed between 2 and 45 ps after the absorption transition are described here as the T(1)→T(3) and T(1)→T(5) transitions, supporting that the intermediate triplet state (T(2)) decays by internal conversion to T(1).     

Two-Gaussian excitations model for the glass transition

We develop a modified "two-state" model with Gaussian widths for the site energies of both ground and excited states, consistent with expectations for a disordered system. The thermodynamic properties of the system are analyzed in configuration space and found to bridge the gap between simple two-state models ("logarithmic" model in configuration space) and the random energy model ("Gaussian" model in configuration space). The Kauzmann singularity given by the random energy model remains for very fragile liquids but is suppressed or eliminated for stronger liquids. The sharp form of constant-volume heat capacity found by recent simulations for binary mixed Lennard-Jones and soft-sphere systems is reproduced by the model, as is the excess entropy and heat capacity of a variety of laboratory systems, strong and fragile. The ideal glass in all cases has a narrow Gaussian, almost invariant among molecular and atomic glassformers, while the excited-state Gaussian depends on the system and its width plays a role in the thermodynamic fragility. The model predicts the possibility of first-order phase transitions for fragile liquids. The analysis of laboratory data for toluene and o-terphenyl indicates that fragile liquids resolve the Kauzmann paradox by a first-order transition from supercooled liquid to ideal-glass state at a temperature between T(g) and Kauzmann temperature extrapolated from experimental data. We stress the importance of the temperature dependence of the energy landscape, predicted by the fluctuation-dissipation theorem, in analyzing the liquid thermodynamics.     

TSNet: predicting transition state structures with tensor field networks and transfer learning

Transition states are among the most important molecular structures in chemistry, critical to a variety of fields such as reaction kinetics, catalyst design, and the study of protein function. However, transition states are very unstable, typically only existing on the order of femtoseconds. The transient nature of these structures makes them incredibly difficult to study, thus chemists often turn to simulation. Unfortunately, computer simulation of transition states is also challenging, as they are first-order saddle points on highly dimensional mathematical surfaces. Locating these points is resource intensive and unreliable, resulting in methods which can take very long to converge. Machine learning, a relatively novel class of algorithm, has led to radical changes in several fields of computation, including computer vision and natural language processing due to its aptitude for highly accurate function approximation. While machine learning has been widely adopted throughout computational chemistry as a lightweight alternative to costly quantum mechanical calculations, little research has been pursued which utilizes machine learning for transition state structure optimization. In this paper TSNet is presented, a new end-to-end Siamese message-passing neural network based on tensor field networks shown to be capable of predicting transition state geometries. Also presented is a small dataset of S_N2 reactions which includes transition state structures – the first of its kind built specifically for machine learning. Finally, transfer learning, a low data remedial technique, is explored to understand the viability of pretraining TSNet on widely available chemical data may provide better starting points during training, faster convergence, and lower loss values. Aspects of the new dataset and model shall be discussed in detail, along with motivations and general outlook on the future of machine learning-based transition state prediction.

Transition states are among the most important molecular structures in chemistry, critical to a variety of fields such as reaction kinetics, catalyst design, and the study of protein function.     

Thermal properties of ionic systems near the liquid-liquid critical point

Isobaric heat capacity per unit volume, C(p), and excess molar enthalpy, h(E), were determined in the vicinity of the critical point for a set of binary systems formed by an ionic liquid and a molecular solvent. Moreover, and, since critical composition had to be accurately determined, liquid-liquid equilibrium curves were also obtained using a calorimetric method. The systems were selected with a view on representing, near room temperature, examples from clearly solvophobic to clearly coulombic behavior, which traditionally was related with the electric permittivity of the solvent. The chosen molecular compounds are: ethanol, 1-butanol, 1-hexanol, 1,3-dichloropropane, and diethylcarbonate, whereas ionic liquids are formed by imidazolium-based cations and tetrafluoroborate or bis-(trifluromethylsulfonyl)amide anions. The results reveal that solvophobic critical behavior-systems with molecular solvents of high dielectric permittivity-is very similar to that found for molecular binary systems. However, coulombic systems-those with low permittivity molecular solvents-show strong deviations from the results usually found for these magnitudes near the liquid-liquid phase transition. They present an extremely small critical anomaly in C(p)-several orders of magnitude lower than those typically obtained for binary mixtures-and extremely low h(E)-for one system even negative, fact not observed, up to date, for any liquid-liquid transition in the nearness of an upper critical solution temperature.     

Comparative study of temperature dependent orientational relaxation in a model thermotropic liquid crystal and in a model supercooled liquid

Recent optical Kerr effect experiments have revealed a power law decay of the measured signal with a temperature independent exponent at short-to-intermediate times for a number of liquid crystals in the isotropic phase near the isotropic-nematic transition and supercooled molecular liquids above the mode coupling theory critical temperature. In this work, the authors investigate the temperature dependence of short-to-intermediate time orientational relaxation in a model thermotropic liquid crystal across the isotropic-nematic transition and in a binary mixture across the supercooled liquid regime in molecular dynamics simulations. The measure of the experimentally observable optical Kerr effect signal is found to follow a power law decay at short-to-intermediate times for both systems in agreement with recent experiments. In addition, the temperature dependence of the power law exponent is found to be rather weak. As the model liquid crystalline system settles into the nematic phase upon cooling, the decay of the single-particle second-rank orientational time correlation function exhibits a pattern that is similar to what has been observed for supercooled liquids.     

Structural studies and polymorphism in amorphous solids and liquids at high pressure

When amorphous materials are compressed their structures are expected to change in response to densification. In some cases, the changes in amorphous structure can be discontinuous and they can even have the character of first-order phase transitions. This is a phenomenon referred to as polyamorphism. Most evidence for polyamorphic transitions between low and high density liquids or analogous transformations between amorphous forms of the same substance to date has been indirect and based on the changes in thermodynamic and other structure-related properties with pressure. Recent studies using advanced X-ray and neutron scattering methods combined with molecular dynamics simulations are now revealing the details of structural changes in polyamorphic systems as a function of pressure. Various "two state" or "two species" models are used to understand the anomalous densification behaviour of liquids with melting curve maxima or regions of negative melting slope. Thermodynamic analysis of the two state model leads to the possibility of low- to high-density liquid transitions caused by differences in bulk thermodynamic properties between different amorphous forms and on the degree of cooperativity between low- and high-density structural configurations. The potential occurrence of first-order transitions between supercooled liquids is identified as a critical-like phenomenon. In this tutorial review we discuss the background to polyamorphism, incorporating the experimental observations, simulation studies and the two-state models. We also describe work carried on several systems that are considered to be polyamorphic.     

Comparison between protein-polyethylene glycol (PEG) interactions and the effect of PEG on protein-protein interactions using the liquid-liquid phase transition

Phase transitions of protein aqueous solutions are important for protein crystallization and biomaterials science in general. One source of thermodynamic complexity in protein solutions and their phase transitions is the required presence of additives such as polyethylene glycol (PEG). To investigate the effects of PEG on the thermodynamic behavior of protein solutions, we report measurements on the liquid-liquid phase separation (LLPS) of aqueous bovine serum albumin (BSA) in the presence of relatively small amounts of PEG with an average molecular weight of 1450 g/mol (PEG1450) as a model system. We experimentally characterize two thermodynamically independent properties of the phase boundary: (1) the effect of PEG1450 concentration on the LLPS temperature, (2) BSA/PEG1450 partitioning in the two liquid coexisting phases. We then use a thermodynamic perturbation theory to relate the first property to the effect of PEG concentration on protein-protein interactions and the second property to protein-PEG interactions. As criteria to determine the accuracy of a microscopic model, we examine the model's ability to describe both experimental thermodynamic properties. We believe that the parallel examination of these two properties is a valuable tool for verifying the validity of existing models and for developing more accurate ones. For our system, we have found that a depletion-interaction model satisfactorily explains both protein-PEG interactions and the effect of PEG concentration on protein-protein interactions. Finally, due to the general importance of LLPS, we will experimentally show that protein-PEG-buffer mixtures can exhibit two distinct types of liquid-liquid phase transitions.     

Density functional theory study of the nematic-isotropic transition in an hybrid cell

We have employed the density functional theory formalism to investigate the nematic-isotropic capillary transitions of a nematogen confined by walls that favor antagonist orientations to the liquid crystal molecules (hybrid cell). We analyze the behavior of the capillary transition as a function of the fluid-substrate interactions and the pore width. In addition to the usual capillary transition between isotropiclike to nematiclike states, we find that this transition can be suppressed when one substrate is wet by the isotropic phase and the other by the nematic phase. Under this condition the system presents interfacelike states which allow us to continuously transform the nematiclike phase to the isotropiclike phase without undergoing a sharp phase transition. Two different mechanisms for the disappearance of the capillary transition are identified. When the director of the nematiclike state is homogeneously planar-anchored with respect to the substrates, the capillary transition ends up in a critical point. This scenario is analogous to the observed in Ising models when confined in slit pores with opposing surface fields which have critical wetting transitions. When the nematiclike state has a linearly distorted director field, the capillary transition continuously transforms in a transition between two nematiclike states.     

Ab Initio simulations of nonstoichiometric Cd(x)Te(1-x) liquids

We present ab initio molecular-dynamics simulations for Cd(x)Te(1-x) liquids where the composition is nonstoichiometric. The simulations are performed following Born-Oppenheimer molecular dynamics. The required forces are obtained from a solution of the Kohn-Sham equation using ab initio pseudopotentials. We consider stoichiometries of the form: Cd(x)Te(1-x), where x=0.2, 0.4, 0.6, and 0.8. For each composition of the melt, we consider a range of temperatures near the experimentally determined liquid temperatures. We examine the microstructural properties of the melt, the viscosity, and self-diffusion properties of the liquid as a function of the stoichiometry and temperature. We also perform an analysis of the distribution of the electronic density of states in these liquids. We find that structural changes in the local order, experimentally predicted to occur when the concentration of Cd is increased, are closely related to changes in the electronic properties of the melt.