Glass-forming liquids, anomalous liquids, and polyamorphism in liquids and biopolymers

Supercritical Gallium Trichloride in Oxidative Metal Recycling: Ga2Cl6 Dimers vs GaCl3 Monomers and Rheological Behavior

Oxidative recycling of metals is crucial for a circular economy, encompassing the preservation of natural resources, the reduction of energy consumption, and the mitigation of environmental impacts and greenhouse gas emissions associated with traditional mining and processing. Low-melting gallium trichloride appears to be a promising oxidative solvent for rare-earth metals, transuranium elements, platinum, pnictogens, and chalcogens. Typically, oxidative dissolution with GaCl3 occurs at relatively low temperatures over a few days, assuming the presence of tetrahedral Ga-Cl entities. While supercritical gallium trichloride holds the potential for advanced recycling, little is known about its structure and viscosity. Using high-energy X-ray diffraction and multiscale modeling, which includes first-principles simulations, we have revealed a dual molecular nature of supercritical gallium trichloride, consisting of tetrahedral dimers and flat trigonal monomers. The molecular geometry can be precisely tuned by adjusting the temperature and pressure, optimizing the recycling process for specific metals. The derived viscosity, consistent with the reported results in the vicinity of melting, decreases by a factor of 100 above the critical temperature, enabling fast molecular diffusion, and efficient recycling kinetics.

Universal correlations between the fragility and interparticle repulsion of glass-forming liquids

A recently published analytical model describing and predicting elasticity, viscosity, and fragility of metallic melts is applied for the analysis of about 30 nonmetallic glassy systems, ranging from oxide network glasses to alcohols, low-molecular-weight liquids, polymers, plastic crystals, and even ionic glass formers. The model is based on the power-law exponent λ representing the steepness parameter of the repulsive part of the inter-atomic or inter-molecular potential and the thermal-expansion parameter α_T determined by the attractive anharmonic part of the effective interaction. It allows fitting the typical super-Arrhenius temperature variation of the viscosity or dielectric relaxation time for various classes of glass-forming matter, over many decades. We discuss the relation of the model parameters found for all these different glass-forming systems to the fragility parameter m and detect a correlation of λ and m for the non-metallic glass formers, in accord with the model predictions. Within the framework of this model, the fragility of glass formers can be traced back to microscopic model parameters characterizing the intermolecular interactions.

Control of electron transfer by protein dynamics in photosynthetic reaction centers

Trehalose and glycerol are low molecular mass sugars/polyols that have found widespread use in the protection of native protein states, in both short- and long-term storage of biological materials, and as a means of understanding protein dynamics. These myriad uses are often attributed to their ability to form an amorphous glassy matrix. In glycerol, the glass is formed only at cryogenic temperatures, while in trehalose, the glass is formed at room temperature, but only upon dehydration of the sample. While much work has been carried out to elucidate a mechanistic view of how each of these matrices interact with proteins to provide stability, rarely have the effects of these two independent systems been directly compared to each other. This review aims to compile decades of research on how different glassy matrices affect two types of photosynthetic proteins: (i) the Type II bacterial reaction center from Rhodobacter sphaeroides and (ii) the Type I Photosystem I reaction center from cyanobacteria. By comparing aggregate data on electron transfer, protein structure, and protein dynamics, it appears that the effects of these two distinct matrices are remarkably similar. Both seem to cause a "tightening" of the solvation shell when in a glassy state, resulting in severely restricted conformational mobility of the protein and associated water molecules. Thus, trehalose appears to be able to mimic, at room temperature, nearly all of the effects on protein dynamics observed in low temperature glycerol glasses.

Effect of Chitin Whiskers on the Molecular Dynamics of Carrageenan-Based Nanocomposites

Films of carrageenan (KC) and glycerol (g) with different contents of chitin nanowhiskers (CHW) were prepared by a solution casting process. The molecular dynamics of pure carrageenan (KC), carrageenan/glycerol (KCg) and KCg with different quantities of CHWs as a filler was studied using dielectric relaxation spectroscopy. The analysis of the CHW effect on the molecular mobility at the glass transition, T_g, indicates that non-attractive intermolecular interactions between KCg and CHW occur. The fragility index increased upon CHW incorporation, due to a reduction in the polymer chains mobility produced by the CHW confinement of the KCg network. The apparent activation energy associated with the relaxation dynamics of the chains at T_g slightly increased with the CHW content. The filler nature effect, CHW or montmorillonite (MMT), on the dynamic mobility of the composites was analyzed by comparing the dynamic behavior of both carrageenan-based composites (KCg/xCHW, KCg/xMMT).

Real-time observation of the isothermal crystallization kinetics in a deeply supercooled liquid

Below the melting temperature T_m, crystals are the stable phase of typical elemental or molecular systems. However, cooling down a liquid below T_m, crystallization is anything but inevitable. The liquid can be supercooled, eventually forming a glass below the glass transition temperature T_g. Despite their long lifetimes and the presence of strong barriers that produces an apparent stability, supercooled liquids and glasses remain intrinsically a metastable state and thermodynamically unstable towards the crystal. Here we investigated the isothermal crystallization kinetics of the prototypical strong glassformer GeO₂ in the deep supercooled liquid at 1100 K, about half-way between T_m and T_g. The crystallization process has been observed through time-resolved neutron diffraction for about three days. Data show a continuous reorganization of the amorphous structure towards the alpha-quartz phase with the final material composed by crystalline domains plunged into a low-density, residual amorphous matrix. A quantitative analysis of the diffraction patterns allows determining the time evolution of the relative fractions of crystal and amorphous, that was interpreted through an empirical model for the crystallization kinetics. This approach provides a very good description of the experimental data and identifies a predator-prey-like mechanism between crystal and amorphous, where the density variation acts as a blocking barrier.

Connectivity, dynamics, and structure in a tetrahedral network liquid

We report a detailed computational study by Brownian dynamics simulations of the structure and dynamics of a liquid of patchy particles which forms an amorphous tetrahedral network upon decreasing the temperature. The highly directional particle interactions allow us to investigate the system connectivity by discriminating the total set of particles into different populations according to a penta-modal distribution of bonds per particle. With this methodology we show how the particle bonding process is not randomly independent but it manifests clear bond correlations at low temperatures. We further explore the dynamics of the system in real space and establish a clear relation between particle mobility and particle connectivity. In particular, we provide evidence of anomalous diffusion at low temperatures and reveal how the dynamics is affected by the short-time hopping motion of the weakly bounded particles. Finally we widely investigate the dynamics and structure of the system in Fourier space and identify two quantitatively similar length scales, one dynamic and the other static, which increase upon cooling the system and reach distances of the order of few particle diameters. We summarize our findings in a qualitative picture where the low temperature regime of the viscoelastic liquid is understood in terms of an evolving network of long time metastable cooperative domains of particles.

Relaxation and physical aging in network glasses: a review

Recent progress in the description of glassy relaxation and aging are reviewed for the wide class of network-forming materials such as GeO2, Ge x Se1-x , silicates (SiO2-Na2O) or borates (B2O3-Li2O), all of which have an important usefulness in domestic, geological or optoelectronic applications. A brief introduction of the glass transition phenomenology is given, together with the salient features that are revealed both from theory and experiments. Standard experimental methods used for the characterization of the slowing down of the dynamics are reviewed. We then discuss the important role played by aspects of network topology and rigidity for the understanding of the relaxation of the glass transition, while also permitting analytical predictions of glass properties from simple and insightful models based on the network structure. We also emphasize the great utility of computer simulations which probe the dynamics at the molecular level, and permit the calculation of various structure-related functions in connection with glassy relaxation and the physics of aging which reveal the non-equilibrium nature of glasses. We discuss the notion of spatial variations of structure which leads to the concept of 'dynamic heterogeneities', and recent results in relation to this important topic for network glasses are also reviewed.

Characterization of the dynamics of glass-forming liquids from the properties of the potential energy landscape

We develop a framework for understanding the difference between strong and fragile behavior in the dynamics of glass-forming liquids from the properties of the potential energy landscape. Our approach is based on a master equation description of the activated jump dynamics among the local minima of the potential energy (the so-called inherent structures) that characterize the potential energy landscape of the system. We study the dynamics of a small atomic cluster using this description as well as molecular dynamics simulations and demonstrate the usefulness of our approach for this system. Many of the remarkable features of the complex dynamics of glassy systems emerge from the activated dynamics in the potential energy landscape of the atomic cluster. The dynamics of the system exhibits typical characteristics of a strong supercooled liquid when the system is allowed to explore the full configuration space. This behavior arises because the dynamics is dominated by a few lowest-lying minima of the potential energy and the potential energy barriers between these minima. When the system is constrained to explore only a limited region of the potential energy landscape that excludes the basins of attraction of a few lowest-lying minima, the dynamics is found to exhibit the characteristics of a fragile liquid.

Polymorphism, thermodynamic anomalies, and network formation in an atomistic model with two internal states

Using molecular dynamics simulations we study, in two dimensions, the temperature-density phase diagram of a simple model with two internal states labeled 1 and -1. The particles interact with a modified Lennard-Jones potential, which depends on relative positions of the particles as well as on their state. Working in an ensemble where the system reduces to a 50:50 binary mixture of two species of particles, we obtain a rich phase diagram showing many features of common network-forming liquids such as water and silica, including polymorphism and thermodynamic anomalies. We believe our model may be useful for studies concerning generic features of such complex liquids.

Thermophysical properties of undercooled alloys: an overview of the molecular simulation approaches

We review the studies on the thermophysical properties of undercooled metals and alloys by molecular simulations in recent years. The simulation methods of melting temperature, enthalpy, specific heat, surface tension, diffusion coefficient and viscosity are introduced and the simulated results are summarized. By comparing the experimental results and various theoretical models, the temperature and the composition dependences of the thermophysical properties in undercooled regime are discussed.

Temperature development of glassy α-relaxation dynamics determined by broadband dielectric spectroscopy

We present the temperature dependence of α -relaxation times of 13 glass formers determined from broadband dielectric spectroscopy, also including data from aging measurements. The data sets partly cover relaxation-time ranges of up to 16 decades enabling a critical test of the validity of model predictions. For this purpose, the data are provided for electronic download. Here we employ these results to test the applicability of the Vogel-Fulcher-Tammann equation and a recently proposed new approach that was demonstrated to provide superior fits of a vast collection of viscosity data.

Glass transition and fragility in the simple molecular glassformer CS(2) from CS(2)-S(2)Cl(2) solution studies

With an interest in finding the fragility for a simple, single component, molecular glassformer, we have determined the dielectric relaxation and glass transition behavior for a series of glasses in the CS(2)-S(2)Cl(2) and CS(2)-toluene systems. Crystallization of CS(2) can be completely avoided down to the composition 20 mol% second component, and the fragility proves almost independent of CS(2) content in each system. Since the glass temperature T(g) obtained from both thermal studies and from dielectric relaxation (using T(g,diel)=T(tau=100 s)) is quite linear over the whole composition range in each system, and since relaxation time data for pure CS(2) fall on the same master plot when scaled by the linearly extrapolated T(g) value, we deduce that pure CS(2) has the same high fragility as the binary solutions. The value is m=86, as for ortho-terphenyl (OTP). Based on observations of independent studies for the vibrational density of states (VDoS) (of inherent structures for OTP and instantaneous, at-temperature structures for CS(2)), we attribute the high fragility to an excess vibrational heat capacity (defined by C(p) (vib, excess)=dS(vib, excess)/d ln T) originating in the behavior of the low frequency modes of the VDoS (the boson peak modes). Both low frequency DoS and anharmonicity increase with increasing temperature, augmenting the configurational entropy drive to the top of the system energy landscape. The surprising implication is that fragility is determined in the vibrational, not configurational, manifold of microstates.

Solidification of Ionic Liquids: Theory and Techniques

Ionic liquids (ILs) have become an important class of solvents and soft materials over the past decades. Despite being salts built by discrete cations and anions, many of them are liquid at room temperature and below. They have been used in a wide variety of applications such as electrochemistry, separation science, chemical synthesis and catalysis, for breaking azeotropes, as thermal fluids, lubricants and additives, for gas storage, for cellulose processing, and photovoltaics. It has been realized that the true advantage of ILs is their modular character. Each specific cation–anion combination is characterized by a unique, characteristic set of chemical and physical properties. Although ILs have been known for roughly a century, they are still a novel class of compounds to exploit due to the vast number of possible ion combinations and one fundamental question remains still inadequately answered: why do certain salts like ILs have such a low melting point and do not crystallize readily? This Review aims to give an insight into the liquid–solid phase transition of ILs from the viewpoint of a solid-state chemist and hopes to contribute to a better understanding of this intriguing class of compounds. It will introduce the fundamental theories of liquid–solid-phase transition and crystallization from melt and solution. Aside form the formation of ideal crystals the development of solid phases with disorder and of lower order like plastic crystals and liquid crystals by ionic liquid compounds are addressed. The formation of ionic liquid glasses is discussed and finally practical techniques, strategies and methods for crystallization of ionic liquids are given.

Viscosity of glass-forming liquids

The low-temperature dynamics of ultraviscous liquids hold the key to understanding the nature of glass transition and relaxation phenomena, including the potential existence of an ideal thermodynamic glass transition. Unfortunately, existing viscosity models, such as the Vogel-Fulcher-Tammann (VFT) and Avramov-Milchev (AM) equations, exhibit systematic error when extrapolating to low temperatures. We present a model offering an improved description of the viscosity-temperature relationship for both inorganic and organic liquids using the same number of parameters as VFT and AM. The model has a clear physical foundation based on the temperature dependence of configurational entropy, and it offers an accurate prediction of low-temperature isokoms without any singularity at finite temperature. Our results cast doubt on the existence of a Kauzmann entropy catastrophe and associated ideal glass transition.

Study of solvent-protein coupling effects by neutron scattering

The present work aims to characterize the dynamical behavior of proteins immersed in bio-preserving liquids and glasses. For this purpose, the protein dUTPase was chosen, while the selected solvents were glycerol, a triol, and some homologous disaccharides, i.e., trehalose, maltose, and sucrose, which are known to be very effective bio-preserving agents. The results highlight that the disaccharides show a slowing down effect on the water dynamics, which is stronger for trehalose than in the case of the other disaccharides. Furthermore, a characterization of the medium which hosts the protein is performed by using an operative definition of fragility based on the mean square displacement extracted by elastic incoherent neutron scattering, which is directly connected to Angell's kinetic fragility based on the viscosity. Finally, a study of the dynamics of the protein sequestered within the solvents is performed. The result shows that the protein dynamics is coupled with that of the surrounding matrix.

Glass transitions in one-, two-, three-, and four-dimensional binary Lennard-Jones systems

We investigate the calorimetric liquid-glass transition by performing simulations of a binary Lennard-Jones mixture in one through four dimensions. Starting at a high temperature, the systems are cooled to T = 0 and heated back to the ergodic liquid state at constant rates. Glass transitions are observed in two, three and four dimensions as a hysteresis between the cooling and heating curves. This hysteresis appears in the energy and pressure diagrams, and the scanning rate dependence of the area and height of the hysteresis can be described using power laws. The one-dimensional system does not experience a glass transition but its specific heat curve resembles the shape of the D≥2 results in the supercooled liquid regime above the glass transition. As D increases, the radial distribution functions reflect reduced geometric constraints. Nearest neighbor distances become smaller with increasing D due to interactions between nearest and next-nearest neighbors. Simulation data for the glasses are compared with crystal and melting data obtained with a Lennard-Jones system with only one type of particle and we find that with increasing D crystallization becomes increasingly more difficult.

Molecular Dynamics Simulation of the Ionic Liquid N-Ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium Bis(trifluoromethanesulfonyl)imide

Thermodynamics, structure, and dynamics of an ionic liquid based on a quaternary ammonium salt with ether side chain, namely, N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide, MOENM2E TFSI, are investigated by molecular dynamics (MD) simulations. Average density and configurational energy of simulated MOENM2E TFSI are interpreted with models that take into account empirical ionic volumes. A throughout comparison of the equilibrium structure of MOENM2E TFSI with previous results for the more common ionic liquids based on imidazolium cations is provided. Several time correlation functions are used to reveal the microscopic dynamics of MOENM2E TFSI. Structural relaxation is discussed by the calculation of simultaneous space-time correlation functions. Temperature effects on transport coefficients (diffusion, conductivity, and viscosity) are investigated. The ratio between the actual conductivity and the estimate from ionic diffusion by the Nernst-Einstein equation indicates that correlated motion of neighboring ions in MOENM2E TFSI is similar to imidazolium ionic liquids. In line with experiment, Walden plot of conductivity and viscosity indicates that simulated MOENM2E TFSI should be classified as a poor ionic liquid.

Effect of chain length on fragility and thermodynamic scaling of the local segmental dynamics in poly(methylmethacrylate)

Local segmental relaxation properties of poly(methylmethacrylate) (PMMA) of varying molecular weight are measured by dielectric spectroscopy and analyzed in combination with the equation of state obtained from PVT measurements. Significant variations of glass transition temperature and fragility with molecular weight are observed. In accord with the general properties of glass-forming materials, single molecular weight dependent scaling exponent gamma is sufficient to define the mean segmental relaxation time taualpha and its distribution. This exponent can be connected to the Gruneisen parameter and related thermodynamic quantities, thus demonstrating the interrelationship between dynamics and thermodynamics in PMMA. Changes in the relaxation properties ("dynamic crossover") are observed as a function of both temperature and pressure, with taualpha serving as the control parameter for the crossover. At longer taualpha another change in the dynamics is apparent, associated with a decoupling of the local segmental process from ionic conductivity.

The Glass Transition Temperature of Polymer Melts

We develop an analytic theory to estimate the glass transition temperature T(g) of polymer melts as a function of the relative rigidities of the chain backbone and side groups, the monomer structure, pressure, and polymer mass. Our computations are based on an extension of the semiempirical Lindemann criterion of melting to locate T(g) and on the use of the advanced mean field lattice cluster theory (LCT) for treating the thermodynamics of systems containing structured monomer, semiflexible polymer chains. The Lindemann criterion is translated into a condition for T(g) by expressing this relation in terms of the specific volume, and this free volume condition is used to calculate T(g) from our thermodynamic theory. The mass dependence of T(g) is compared to that of other characteristic temperatures of glass-formation. These additional characteristic temperatures are determined from the temperature variation of the LCT configurational entropy, in conjunction with the Adam-Gibbs model for long wavelength structural relaxation. Our theory explains generally observed trends in the variation of T(g) with polymer microstructure, and we find that T(g) can be tuned either upward or downward by increasing the length of the side chains, depending on the relative rigidities of the side groups and the chain backbone. The elucidation of the molecular origins of T(g) in polymer liquids should be useful in designing and processing new synthetic materials and for understanding the dynamics and controlling the preservation of biological substances.

Mean-square displacement relationship in bioprotectant systems by elastic neutron scattering

Neutron intensity elastic scans on trehalose, maltose, and sucrose/H(2)O mixtures as a function of concentration, temperature, and exchanged wave vector are presented. The experimental findings show a crossover in molecular fluctuations between harmonic and anharmonic dynamical regimes. A new operative definition for the degree of fragility of glass-forming systems is furnished by using explicitly the connection between viscosity and mean-square displacement. The procedure is tested for the investigated mixtures and for a set of glass-forming systems. In this frame, the stronger character of trehalose/H(2)O mixture indicates a better attitude in respect to maltose and sucrose/H(2)O mixtures to encapsulate biostructures in a more rigid matrix.

Fragility by elastic incoherent neutron scattering

The present work furnishes an operative definition for the fragility degree by using elastic incoherent neutron scattering (EINS). Such a definition is based on the relation between viscosity, a macroscopic quantity, and the atomic mean-square displacement, which refers to a nanoscopic property. This procedure has been used to analyze a set of glass-forming systems and it allows to obtain a linear dependence of the fragility parameter M, obtained by EINS, on the fragility parameter m, obtained by viscosity measurements.

AMORPHOUS WATER

After providing some background material to establish the interest content of this subject, we summarize the many different ways in which water can be prepared in the amorphous state, making clear that there seems to be more than one distinct amorphous state to be considered. We then give some space to structural and spectroscopic characterization of the distinct states, recognizing that whereas there seems to be unambiguously two distinct states, there may be in fact be more, the additional states mimicking the structures of the higher-density crystalline polymorphs. The low-frequency vibrational properties of the amorphous solid states are then examined in some detail because of the gathering evidence that glassy water, while difficult to form directly from the liquid like other glasses, may have some unusual and almost ideal glassy features, manifested by unusually low states of disorder. This notion is pursued in the following section dealing with thermodynamic and relaxational properties, where the uniquely low excess entropy of the vitreous state of water is confirmed by three different estimates. The fact that the most nearly ideal glass known has no properly established glass transition temperature is highlighted, using known dielectric loss data for amorphous solid water (ASW) and relevant molecular glasses. Finally, the polyamorphism of glassy water, and the kinetic aspects of transformation from one form to the other, are reviewed.

Correlation between nonexponential relaxation and non-Arrhenius behavior under conditions of high compression

Photon correlation spectroscopy was used to investigate the behavior of the dynamical properties of 1,1'-di(4-methoxy-5-methyl-phenyl)cyclohexane (BMMPC) at elevated pressures. The fragility of BMMPC measured by the steepness index m(T) is decreasing and the nonexponentiality parameter beta(KWW) is increasing with increasing pressure. This result strongly suggests that the phenomenological correlation between the steepness index and nonexponentionality is also preserved under high compression. The pressure dependence of the structural relaxation times is well characterized by a simple activation volume form. The activation volume continuously increases with decreasing temperature, which is probably due to the increase of cooperativity of the structural relaxation process. Moreover, we found that the glass-transition temperature exhibits a significant dependence on pressure.

Glassiness and constrained dynamics of a short-range nondisordered spin model

We study the low temperature dynamics of a two-dimensional short-range spin system with uniform ferromagnetic interactions, which displays glassiness at low temperatures despite the absence of disorder or frustration. The model has a dual description in terms of free defects subject to dynamical constraints, and is an explicit realization of the "hierarchically constrained dynamics" scenario for glassy systems. We give a number of exact results for the statics of the model, and study in detail the dynamical behavior of one-time and two-time quantities. We also consider the role played by the configurational entropy, which can be computed exactly, in the relation between fluctuations and response.

Glasslike dynamical behavior of the plastocyanin hydration water

The dynamical behavior of water around plastocyanin has been investigated in a wide temperature range by molecular dynamics simulation. The mean square displacements of water oxygen atoms show, at long times, a t(alpha) trend for all temperatures. Below 150 K, alpha is constant and equal to 1; at higher temperatures it drops to a value significantly smaller than 1, and thereafter decreases with increasing temperature. The occurrence of such an anomalous diffusion matches the onset of the dynamical transition observed in the protein. The intermediate scattering function of water is characterized, at high temperature, by a stretched exponential decay evolving, at low temperature, toward a two step relaxation behavior, which becomes more evident on increasing the exchanged wave vector q. Both the mean square displacements and the intermediate scattering functions show, beyond the ballistic regime, a plateau, which progressively extends for longer times as long as the temperature is lowered, such behavior reflecting trapping of water molecules within a cage formed by the nearest neighbors. At low temperature, a low frequency broad inelastic peak is observed in the dynamical structure factor of hydration water; such an excess of vibrational modes being reminiscent of the boson peak, characteristic of disordered, amorphous systems. All these features, which are typical of complex systems, can be traced back to the glassy character of the hydration water and suggest a dynamical coupling occurring at the macromolecule-solvent interface.