Computer simulation study of thermodynamic scaling of dynamics of 2Ca(NO3)2·3KNO3

Molecular dynamics (MD) simulations of the glass-former 2Ca(NO(3))(2·3KNO(3), CKN, were performed as a function of temperature at pressures 0.1 MPa, 0.5 GPa, 1.0 GPa, and 2.0 GPa. Diffusion coefficient, relaxation time of the intermediate scattering function, and anion reorientational time were obtained as a function of temperature and densitiy ρ. These dynamical properties of CKN scale as ρ(γ)∕T with a common value γ = 1.8 ± 0.1. The scaling parameter γ is consistent with the exponent of the repulsive part of an effective intermolecular potential for the repulsion between the atoms at shortest distance in the equilibrium structure of liquid CKN, Ca(2+), and oxygen atoms of NO(3)(-). Correlation between potential energy and virial is obeyed for the short-range terms of the potential function, but not for the whole potential including coulombic interactions. Decoupling of diffusion coefficient and reorientational relaxation time from relaxation time take place at a given ρ(γ)∕T value, i.e., breakdown of Stokes-Einstein and Debye-Stokes-Einstein equations result from combined thermal and volume effects. The MD results agree with correlations proposed between long-time relaxation and short-time dynamics, lnτ ∝ 1∕, where the mean square displacement concerns a time window of 10.0 ps. It has been found that scales as ρ(γ)∕T above and below the glass transition temperature, so that thermodynamic scaling of liquid dynamics can be thought as a consequence of theories relating short- and long-time dynamics, and the more fundamental scaling concerns short-time dynamical properties.     

Activation energy-activation volume master plots for ion transport behavior in polymer electrolytes and supercooled molten salts

We demonstrate the use of activation energy versus activation volume "master plots" to explore ion transport in typical fragile glass forming systems exhibiting non-Arrhenius behavior. These systems include solvent-free salt complexes in poly(ethylene oxide) (PEO) and low molecular weight poly(propylene oxide) (PPO) and molten 2Ca(NO3)2.3KNO3 (CKN). Plots showing variations in apparent activation energy EA versus apparent activation volume VA are straight lines with slopes given by M = DeltaEA/DeltaVA. A simple ion transport mechanism is described where the rate determining step involves a dilatation (expressed as VA) around microscopic cavities and a corresponding work of expansion (EA). The slopes of the master plots M are equated to internal elastic moduli, which vary from 1.1 GPa for liquid PPO to 5.0 GPa for molten CKN on account of differing intermolecular forces in these materials.     

Orientational dynamics and energy landscape features of thermotropic liquid crystals: An analogy with supercooled liquids

Recent optical kerr effect (OKE) studies have revealed that orientational relaxation of rodlike nematogens near the isotropic-nematic (I-N) phase boundary and also in the nematic phase exhibit temporal power law decay at intermediate times. Such behaviour has drawn an intriguing analogy with supercooled liquids. Here, we have investigated the single-particle and collective orientational dynamics of a family of model system of thermotropic liquid crystals using extensive computer simulations. Several remarkable features of glassy dynamics are on display including non-exponential relaxation, dynamical heterogeneity, and non-Arrhenius temperature dependence of the orientational relaxation time. Over a temperature range near the I-N phase boundary, the system behaves like a fragile glass-forming liquid. Using proper scaling, we construct the usual relaxation time versus inverse temperature plot and explicitly demonstrate that one can successfully define a density dependent fragility of liquid crystals. The fragility of liquid crystals shows a temperature and density dependence which is remarkably similar to the fragility of glass forming supercooled liquids. Energy landscape analysis of inherent structures shows that the breakdown of the Arrhenius temperature dependence of relaxation rate occurs at a temperature that marks the onset of the growth of the depth of the potential energy minima explored by the system.     

Conformational equilibrium of cytochrome P450 BM-3 complexed with N-palmitoylglycine: a replica exchange molecular dynamics study

UV-vis absorbance measurements and associated studies of cytochrome P450 BM-3 in complex with N-palmitoylglycine (NPG) indicate that a conformational change occurs in the active site of the complex where the terminal atoms of the ligand move from a site distant from the heme iron, as seen in the low temperature crystal structure to a site proximal to the heme iron at biological temperatures. We employ replica exchange molecular dynamics simulations to study this conformational change. The population of the proximal state is found to increase with temperature in agreement with UV-vis absorbance and NMR measurements. In addition to the conformations characterized by X-ray crystallography and computer modeling, this study shows that a new conformational state is significantly populated at room temperature. The observed increase in the population of conformations where the terminal atoms of NPG are proximal to the heme iron with increasing temperature indicates that the proximal state is stabilized by conformational entropy. A proposal for the origin of this entropic stabilization is provided on the basis of the structure of the newly identified state. We use the temperature weighted histogram (T-WHAM) method to characterize the transition state regions of the conformational ensemble and propose a mechanism of interconversion between these low free energy conformational states.     

Decoupling charge transport from the structural dynamics in room temperature ionic liquids

Light scattering and dielectric spectroscopy measurements were performed on the room temperature ionic liquid (RTIL) [C4mim][NTf2] in a broad temperature and frequency range. Ionic conductivity was used to estimate self-diffusion of ions, while light scattering was used to study structural relaxation. We demonstrate that the ionic diffusion decouples from the structural relaxation process as the temperature of the sample decreases toward T(g). The strength of the decoupling appears to be significantly lower than that expected for a supercooled liquid of similar fragility. The structural relaxation process in the RTIL follows well the high-temperature mode coupling theory (MCT) scenario. Using the MCT analysis we estimated the dynamic crossover temperature in [C4mim][NTf2] to be T(c) ~ 225 ± 5 K. However, our analysis reveals no sign of the dynamic crossover in the ionic diffusion process.     

Gaussian excitations model for glass-former dynamics and thermodynamics

We describe a model for the thermodynamics and dynamics of glass-forming liquids in terms of excitations from an ideal glass state to a Gaussian manifold of configurationally excited states. The quantitative fit of this three parameter model to the experimental data on excess entropy and heat capacity shows that "fragile" behavior, indicated by a sharply rising excess heat capacity as the glass transition is approached from above, occurs in anticipation of a first-order transition--usually hidden below the glass transition--to a "strong" liquid state of low excess entropy. The distinction between fragile and strong behavior of glass formers is traced back to an order of magnitude difference in the Gaussian width of their excitation energies. Simple relations connect the excess heat capacity to the Gaussian width parameter, and the liquid-liquid transition temperature, and strong, testable, predictions concerning the distinct properties of energy landscape for fragile liquids are made. The dynamic model relates relaxation to a hierarchical sequence of excitation events each involving the probability of accumulating sufficient kinetic energy on a separate excitable unit. Super-Arrhenius behavior of the relaxation rates, and the known correlation of kinetic with thermodynamic fragility, both follow from the way the rugged landscape induces fluctuations in the partitioning of energy between vibrational and configurational manifolds. A relation is derived in which the configurational heat capacity, rather than the configurational entropy of the Adam-Gibbs equation, controls the temperature dependence of the relaxation times, and this gives a comparable account of the experimental observations without postulating a divergent length scale. The familiar coincidence of zero mobility and Kauzmann temperatures is obtained as an approximate extrapolation of the theoretical equations. The comparison of the fits to excess thermodynamic properties of laboratory glass formers, and to configurational thermodynamics from simulations, reveals that the major portion of the excitation entropy responsible for fragile behavior resides in the low-frequency vibrational density of states. The thermodynamic transition predicted for fragile liquids emerges from beneath the glass transition in case of laboratory water and the unusual heat capacity behavior observed for this much studied liquid can be closely reproduced by the model.     

An upper limit to kinetic fragility in glass-forming liquids

The kinetic fragility of a liquid is correlated to the magnitude of enthalpy hysteresis in various glass-forming materials during thermal cycling across the glass transition. While the lower bound of liquid fragility is well known, there has been little research into the possibility of an inherent upper limit to fragility. In this paper, we present a theoretical argument for the existence of a maximum fragility and show that the correlation between fragility and enthalpy hysteresis allows for an empirical evaluation of the upper limit of fragility. This upper limit occurs as the enthalpy hysteresis involved in thermal cycling about the glass transition approaches zero, leading to m(max)≈175. This result agrees remarkably well with our previous estimate. The dynamics of maximum fragility liquids are discussed, and a critical temperature of ～1.5 T(g) (where T(g) is the glass transition temperature) is revealed where a transition from nonexponential to exponential structural relaxation occurs.     

Thermodynamics,viscous flow and relaxation dynamics of bulk glass-forming Pd alloys

Experimental results concerning the thermodynamic excess properties and the equilibrium shear viscosity of the undercooled liquid state are presented, which were consistently obtained on Pd-rich bulk metallic glasses of different stability by applying the concept of the limiting fictive temperature. These results are evaluated with respect to the stability criteria that are controlling vitrification. In addition, modulated-temperature calorimetry has been performed to obtain the average relaxation time of the cooperative liquid modes near the glass transition directly in dependence of temperature. The results are compared to data obtained from static non-equilibrium measurements and discussed with respect to the relation between fragility and thermal stability of Pd-rich bulk glass-forming alloys.     

Light-emission and excited-state dynamics in Tm2+ doped CsCaCl3, CsCaBr3, and CsCaI3

The light-emission and photophysical properties of CsCaCl3:Tm2+ (1.04%), CsCaBr3:Tm2+ (0.48%), and CsCaI3:Tm2+ (0.76%) are presented. We find that Tm2+ is a multiple emitter under 21,834 cm-1 laser excitation at low temperatures in all three compounds. Several distinct types of emission are observed and characterized: sharp and long-lived 4f-4f emission in the infrared (IR) and up to four broad and fast decaying emission bands in the near-IR and visible, originating from the 4f-5d states of Tm2+. The optical spectroscopic properties of the samples are compared, and we find that the measured differences in the relative intensities and the shifts in the position of the emissions can be related to the chemical influence on the absorption and emission properties of Tm2+. Thus, it nicely illustrates the principle of chemical variation on the optical spectroscopic properties. An investigation of the temperature dependence of the luminescence yields important information about the dynamics of the excited states. The interplay and competition between radiative and nonradiative pathways is explained and modeled using a single configurational coordinate approach.     

Thermodynamic scaling of α-relaxation time and viscosity stems from the Johari-Goldstein β-relaxation or the primitive relaxation of the coupling model

By now it is well established that the structural α-relaxation time, τ(α), of non-associated small molecular and polymeric glass-formers obey thermodynamic scaling. In other words, τ(α) is a function Φ of the product variable, ρ(γ)/T, where ρ is the density and T the temperature. The constant γ as well as the function, τ(α) = Φ(ρ(γ)/T), is material dependent. Actually this dependence of τ(α) on ρ(γ)/T originates from the dependence on the same product variable of the Johari-Goldstein β-relaxation time, τ(β), or the primitive relaxation time, τ(0), of the coupling model. To support this assertion, we give evidences from various sources itemized as follows. (1) The invariance of the relation between τ(α) and τ(β) or τ(0) to widely different combinations of pressure and temperature. (2) Experimental dielectric and viscosity data of glass-forming van der Waals liquids and polymer. (3) Molecular dynamics simulations of binary Lennard-Jones (LJ) models, the Lewis-Wahnstr?m model of ortho-terphenyl, 1,4 polybutadiene, a room temperature ionic liquid, 1-ethyl-3-methylimidazolium nitrate, and a molten salt 2Ca(NO(3))(2)·3KNO(3) (CKN). (4) Both diffusivity and structural relaxation time, as well as the breakdown of Stokes-Einstein relation in CKN obey thermodynamic scaling by ρ(γ)/T with the same γ. (5) In polymers, the chain normal mode relaxation time, τ(N), is another function of ρ(γ)/T with the same γ as segmental relaxation time τ(α). (6) While the data of τ(α) from simulations for the full LJ binary mixture obey very well the thermodynamic scaling, it is strongly violated when the LJ interaction potential is truncated beyond typical inter-particle distance, although in both cases the repulsive pair potentials coincide for some distances.     

Novel method to detect the volumetric glass --> liquid transition at high pressures: glycerol as a test case

We report a novel method of detecting the glass --> liquid transition at high pressures, which comprises measuring the relative volume change incurred upon heating glassy samples into the liquid state. We show data on glycerol in the pressure range 0.050-1.00 GPa to demonstrate the viability of the method. The reversible glass --> liquid transition is observed by means of a kink in the relative volume change on heating the sample isobarically, which is attributed to the glass --> liquid transition temperature Tg. This kink can only be observed in the second and subsequent heating cycles since it is superposed by a compaction in the first heating cycle. The isobaric thermal expansivity beta, which is closely related to the first derivative of this curve, shows the features expected for a glass --> liquid transition, including a sharp rise of beta(glass) in a narrow temperature interval to beta(viscous liquid) and an accompanying overshoot effect. Both Tg and the size of the overshoot effect vary in accordance with theory upon changing the ratio of cooling to heating rates. From the shape of this curve the onset, inflection, overshoot peak, and endpoint of the glass --> liquid transition can be extracted, which can be employed to calculate the reduced glass transition width as a measure for the fragility of the liquid. Comparison with literature data allows quantifying the accuracy of the liquid's thermal expansivity beta to be at least +/-10%, while the error in beta is significantly larger for the expansivity of the glassy state. The reproducibility of the glass --> liquid transition temperature Tg is better than +/-2 K. Our glycerol data confirms literature studies showing a nonlinear increase of Tg with increasing pressure (approximately 35 K/GPa on average), which is accompanied by an increase in fragility.     

Seven-dimensional quantum dynamics study of the O(3P)+CH4 reaction

The initial state selected time-dependent wave packet calculations have been carried out to study the title reaction with seven degrees of freedom included by restricting the nonreacting CH(3) group under C(3V) symmetry and the CH bond length in the group. Total reaction probabilities as well as integral cross sections were calculated for the ground and four vibrationally excited reagent states. Our calculation shows that the reactivity is very small for the reaction for collision energy up to 1.0 eV for all the initial states. Initial vibration excitation of CH(4), in particular, the CH stretch excitation, enhances the reactivity, but only part of the excitation energy deposited can be used to reduce the reaction threshold. The rate constant for the ground initial state agrees rather well with that from a recent quasiclassical trajectory study and is larger than that from the semirigid vibrating rotor target calculations, in particular, in the low temperature region. On the other hand, the thermal rate constant calculated from the integral cross sections for these five vibrational states is about a factor of 20 smaller than that from the multiconfiguration time-dependent Hartree study.     

Seven-dimensional quantum dynamics study of the H+NH3-->H2+NH2 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H+NH3-->H2+NH2 reaction using a seven-dimensional model and an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values. The total reaction probabilities are calculated for the initial ground and seven excited states of NH3 with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants are calculated for the temperature range 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the total reaction probabilities are overall very small, (b) the symmetric and asymmetric NH stretch excitations enhance the reaction significantly and almost all of the excited energy deposited was used to reduce the reaction threshold, (c) the excitation of the umbrella and bending motion have a smaller contribution to the enhancement of reactivity, (d) the main contribution to the thermal rate constants is thought to come from the ground state at low temperatures and from the stretch excited states at high temperatures, and (e) the calculated thermal rate constants are three to ten times smaller than the experimental data and transition state theory results.     

The effects of surface temperature on the gas-liquid interfacial reaction dynamics of O(3P)+squalane

OH/OD product state distributions arising from the reaction of gas-phase O(3P) atoms at the surface of the liquid hydrocarbon squalane C30H62/C30D62 have been measured. The O(3P) atoms were generated by 355 nm laser photolysis of NO2 at a low pressure above the continually refreshed liquid. It has been shown unambiguously that the hydroxyl radicals detected by laser-induced fluorescence originate from the squalane surface. The gas-phase OH/OD rotational populations are found to be partially sensitive to the liquid temperature, but do not adapt to it completely. In addition, rotational temperatures for OH/OD(v'=1) are consistently colder (by 34+/-5 K) than those for OH/OD(v'=0). This is reminiscent of, but less pronounced than, a similar effect in the well-studied homogeneous gas-phase reaction of O(3P) with smaller hydrocarbons. We conclude that the rotational distributions are composed of two different components. One originates from a direct abstraction mechanism with product characteristics similar to those in the gas phase. The other is a trapping-desorption process yielding a thermal, Boltzmann-like distribution close to the surface temperature. This conclusion is consistent with that reached previously from independent measurements of OH product velocity distributions in complementary molecular-beam scattering experiments. It is further supported by the temporal profiles of OH/OD laser-induced fluorescence signals as a function of distance from the surface observed in the current experiments. The vibrational branching ratios for (v'=1)/(v'=0) for OH and OD have been found to be (0.07+/-0.02) and (0.30+/-0.10), respectively. The detection of vibrationally excited hydroxyl radicals suggests that secondary and/or tertiary hydrogen atoms may be accessible to the attacking oxygen atoms.     

Direct gas-liquid interfacial dynamics: the reaction between O(3P) and a liquid hydrocarbon

We report the first measurements of internal energy distributions of the OH produced via a direct mechanism, isolated from other components on the basis of time-of-flight, in the interfacial reaction between gas-phase O((3)P) atoms and the liquid hydrocarbon squalane, C(30)H(62). O((3)P) atoms were generated by laser photolysis of NO(2) above the liquid. Resulting hydroxyl radicals that escape from the surface were detected by laser-induced fluorescence. Time-of-flight profiles demonstrate that the kinetic energy of the fastest OH (nu' = 1) is lower than that of (nu' = 0). Rotational distributions were measured at the rising edge of their appearance for both OH (nu' = 0) and (nu' = 1). They were found to differ substantially more than at the peak of their profiles. They were also less dependent on the bulk liquid temperature. We conclude that the new data confirm strongly that at least two mechanisms contribute to the production of OH. The higher-velocity component has translational and rotational energy distributions, observed cleanly for the first time, consistent with a direct mechanism. The close correspondence of these rotational distributions to those from the corresponding homogeneous gas-phase reaction of O((3)P) with smaller hydrocarbons suggests a very similar, near collinear direct abstraction. This is accompanied by a slower component with kinetic energy and rotational (but not vibrational) distributions reflecting the temperature of the liquid, consistent with a distinct trapping-desorption mechanism.     

Local dielectric spectroscopy of near-surface glassy polymer dynamics

A noncontact scanning probe microscopy method was used to probe local near-surface dielectric susceptibility and dielectric relaxation in polyvinyl acetate near the glass transition. Dielectric spectra were measured from 10(-4) to 10(2) Hz as a function of temperature. The measurements probed a 20 nm thick layer below the free surface of a bulk film. A significant change in the fragility index and moderate narrowing of the distribution of relaxation times were found in the near-surface layer. In contrast to results for ultrathin films confined on or between metallic electrodes, no reduction in the dielectric strength was found, inconsistent with the immobilization of slower modes.     

Crystal growth kinetics exhibit a fragility-dependent decoupling from viscosity

In this paper we establish the temperature dependence of the kinetic coefficient associated with crystal growth into the supercooled liquid for a wide range of organic and inorganic materials. We show that the kinetic coefficient for crystal growth scales with the shear viscosity eta as eta(-xi) and that the exponent depends systematically on the fragility of the liquid. The greater the fragility (i.e., deviation away from an Arrhenius temperature dependence for eta), the larger the difference 1-xi. We argue that this breakdown in scaling between the crystal growth kinetics and the viscosity is a manifestation of heterogeneous dynamics in supercooled liquids. In addition, we show that the absolute growth rate at intermediate viscosities is correlated with the entropy difference between the liquid and the crystal.     

Energy relaxation of intermolecular motions in supercooled water and ice: a molecular dynamics study

We investigate the energy relaxation of intermolecular motions in liquid water at temperatures ranging from 220 K to 300 K and in ice at 220 K using molecular dynamics simulations. We employ the recently developed frequency resolved transient kinetic energy analysis, which provides detailed information on energy relaxation in condensed phases like two-color pump-probe spectroscopy. It is shown that the energy cascading in liquid water is characterized by four processes. The temperature dependences of the earlier three processes, the rotational-rotational, rotational-translational, and translational-translational energy transfers, are explained in terms of the density of states of the intermolecular motions. The last process is the slow energy transfer arising from the transitions between potential energy basins caused by the excitation of the low frequency translational motion. This process is absent in ice because the hydrogen bond network rearrangement, which accompanies the interbasin transitions in liquid water, cannot take place in the solid phase. We find that the last process in supercooled water is well approximated by a stretched exponential function. The stretching parameter, β, decreases from 1 to 0.72 with decreasing temperature. This result indicates that the dynamics of liquid water becomes heterogeneous at lower temperatures.     

Dependence of the fragility of a glass former on the softness of interparticle interactions

We study the influence of the softness of the interparticle interactions on the fragility of a glass former by considering three model binary mixture glass formers. The interaction potential between particles is a modified Lennard-Jones type potential, with the repulsive part of the potential varying with an inverse power q of the interparticle distance, and the attractive part varying with an inverse power p. We consider the combinations (12,11) (model I), (12,6) (model II), and (8,5) (model III) for (q,p) such that the interaction potential becomes softer from model I to III. We evaluate the kinetic fragilities from the temperature variation of diffusion coefficients and relaxation times, and a thermodynamic fragility from the temperature variation of the configurational entropy. We find that the kinetic fragility increases with increasing softness of the potential, consistent with previous results for these model systems, but at variance with the thermodynamic fragility, which decreases with increasing softness of the interactions, as well as expectations from earlier results. We rationalize our results by considering the full form of the Adam-Gibbs relation, which requires, in addition to the temperature dependence of the configurational entropy, knowledge of the high temperature activation energies in order to determine fragility. We show that consideration of the scaling of the high temperature activation energy with the liquid density, analyzed in recent studies, provides a partial rationalization of the observed behavior.     

A Phosphorescent Molecular “Butterfly” that undergoes a Photoinduced Structural Change allowing Temperature Sensing and White Emission

A butterfly‐like phosphorescent platinum(II) binuclear complex can undergo a molecular structure change in which the Pt–Pt distance shortens upon photoexcitation, which leads to the formation of two distinct excited states and dual emission in the steady state, that is, greenish‐blue emission from the high‐energy excited state at the long Pt–Pt distance and red emission from the low‐energy excited state at the short Pt–Pt distance. This photoinduced molecular structure change has a strong dependence on the molecule’s surrounding environment, allowing its application as self‐referenced luminescent sensor for solid–liquid phase change, viscosity, and temperature, with greenish‐blue emission in solid matrix and rising red emission in molten liquid phase. With proper control of the surrounding media to manipulate the structural change and photophysical properties, a broad white emission can be achieved from this molecular butterfly.