Role of molecular reorientation in the vitrification of molecular liquids: A Monte Carlo study of the liquid and glassy states of 2,2-dimethylbutane

The liquid and the glassy phases of 2,2-dimethylbutane have been investigated by isothermal isobaric ensemble Monte Carlo simulation. Thermodynamic Properties and radial distribution functions for both the liquid and the glass have been obtained. The radial distribution functions have been classified into three types based on the accessibility of the group. It has been shown that the structure of the Iiquid and the glass can be understood in terms of the above classification of the radial distribution functions. Molecular reorientation plays an important role in the structural rearrangement accompanying glass formation. As much as 35% of the contribution to the increase in the intermolecular interaction energy on vitrification is due to the reorientation of the neighbouring pairs of molecules. The observed changes in the dimerisation energy and the bonding energy distribution function are consistent with the observed structural changes.

Molecular interpretation of the linear relationship between the entropy and the enthalpy of activation of charge transfer reactions in polar liquids

The well-known linear relationship (T?S# =??H# +?, where 1 >? > 0,? > 0) between the entropy (?S#) and the enthalpy (?H#) of activation for reactions in polar liquids is investigated by using a molecular theory. An explicit derivation of this linear relation from first principles is presented for an outersphere charge transfer reaction. The derivation offers microscopic interpretation for the quantities? and?. It has also been possible to make connection with and justify the arguments of Bell put forward many years ago.

Collective Effects on Single Particle Orientational Relaxation in Slow Dipolar Liquids

Theoretical and computer simulation studies of orientational relaxation in dense molecular liquids are presented. The emphasis of the study is to understand the effects of collective orientational relaxation on the single-particle orientational dynamics. The theoretical analysis is based on a recently developed molecular hydrodynamic theory which allows a self-consistent description of both the collective and the single-particle orientational relaxation. The molecular hydrodynamic theory can be used to derive a relation between the memory function for the collective orientational correlation function and the frequency-dependent dielectric function. A novel feature of the present work is the demonstration that this collective memory function is significantly different from the single-particle rotational friction. However, a microscopic expression for the single-particle rotational friction can be derived from the molecular hydrodynamic theory where the collective memory function can be used to obtain the single-particle orientational friction. This procedure allows, us to calculate the single-particle orientational correlation function near the alpha-beta transition in the supercooled liquid. The calculated correlation function shows an interesting bimodal decay below the bifurcation temperature as the glass transition is approached from above. Brownian dynamics simulations have been carried out to check the validity of the above procedure of translating the memory function from the dielectric relaxation data. We have also investigated the following two issues important in understanding the orientational relaxation in slow liquids. First, we present an analysis of the ''orientational caging'' of translational motion. The value of the translational friction is found to be altered significantly by the orientational caging. Second, we address the question of the rank dependence of the dielectric friction using both simulation and the molecular hydrodynamic theory.

Solvation dynamics in slow, viscous liquids: Application to amides

As the viscosity of a liquid increases rapidly in the supercooled regime, the nature of molecular relaxation can exhibit dynamics rather different from the fast dynamics observed in the normal regime. In this article, we present theoretical studies of solvation dynamics and orientational relaxation in slow liquids. As the local short-range correlations are important in the slow liquids, we have extended our previous theory to take into account the shea-range pair correlations between the polar solute and the dipolar solvent molecules. Application of the generalized theory To the study of solvation dynamics of amide systems gives nice agreement with the experimental results of Maroncelli and co-workers (J. Phys. Chem. 1990, 94, 4929). The theory also provides valuable insight into the orientational relaxation precesses in the viscous liquids.

Structure-potential relationship and nature of disorder in liquids and glassy solids

A method based on the minimal-spanning tree is extended to a collection of points in three dimensions. Two parameters, the average edge length and its standard deviation characterize the disorder. The structural phase diagram for a monatomic system of particles and the characteristic values for the uniform random distribution of points have been obtained. The method is applied to hard spheres and Lennard-Jones systems. These systems occupy distinct regions in the structural phase diagram. The structure of the Lennard-Jones system approaches that of the defective close-packed arrangements at low temperatures whereas in the liquid regime, it deviates from the close-packed configuration.

Thermodynamics of flow and vaporization processes in long-chain liquids

The flow and vaporization behaviors of long-chain esters of varying molecular weights (300-900) ana branching (linear, Y-shaped, and +-shaped molecules) have been studied. The flow behavior is found to depend on the structure as well as the molecular weight. Below a molecular weight of 600, the molecules flow wholly but above this, segmental motion occurs, and the flow becomes independent of the molecular weight which is explained from the blob model. The blob concept demonstrates that the hole of a size of about 11 angstrom is needed for the flow to occur and it is much less than the size of the molecule. The blob size is observed to slightly decrease along the series linear and Y- and +-branched esters. The heat of vaporization is found to be independent of the molecular structure since the molecules acquire a coiled spherical shape during vaporization and hence depends only on the molecular weight. A significant structural effect is observed for the esters on their glass transition temperature (T(g)). The T(g) vs molecular weight plot displays contrasting trend for linear and +-branched esters, with Y esters showing an intermediate behavior. It is explained from their molecular packing and entanglement as visualized by the blob model.

Translational and rotational motion in molecular liquids: A computer simulation study of Lennard-Jones ellipsoids

In order to understand the translational and rotational motion in dense molecular liquids, detailed molecular dynamics simulations of Lennard-Jones ellipsoids have been carried out for three different values of the aspect ratio kappa. For ellipsoids with an aspect ratio equal to 2, the product of the translational diffusion coefficient (D-T) and the average orientational correlation time of the l-th rank harmonics (tau(lR)), converges to a nearly constant value at high density. Surprisingly, this density independent value of D-T tau(lR) is within 5% of the hydrodynamic prediction with the slip boundary condition. This is despite the fact that both D-T and tau(lR) themselves change nearly by an order of magnitude in the density range considered, and the rotational correlation function itself is strongly nonexponential. For small aspect ratios (kappa less than or equal to 1.5), the rotational correlation function remains largely Gaussian even at a very large density, while for a large aspect ratio (kappa greater than or equal to 3), the transition to the nematic liquid-crystalline phase precludes the hydrodynamic regime. Thus, the rotational dynamics of ellipsoids show great sensitivity to the aspect ratio. At low density, tau(lR) goes through a minimum value, indicating the role of interactions in enhancing the rate of orientational relaxation. (C) 1997 American Institute of Physics. [S0021-9606(97)50142-5].

Isomerization dynamics in viscous liquids: Microscopic investigation of the coupling and decoupling of the rate to and from solvent viscosity and dependence on the intermolecular potential

A detailed investigation of viscosity dependence of the isomerization rate is carried out for continuous potentials by using a fully microscopic, self-consistent mode-coupling theory calculation of both the friction on the reactant and the viscosity of the medium. In this calculation we avoid approximating the short time response by the Enskog limit, which overestimates the friction at high frequencies. The isomerization rate is obtained by using the Grote-Hynes formula. The viscosity dependence of the rate has been investigated for a large number of thermodynamic state points. Since the activated barrier crossing dynamics probes the high-frequency frictional response of the liquid, the barrier crossing rate is found to be sensitive to the nature of the reactant-solvent interaction potential. When the solute-solvent interaction is modeled by a 6-12 Lennard-Jones potential, we find that over a large variation of viscosity (eta), the rate (k) can indeed be fitted very well to a fractional viscosity dependence: (k similar to eta(-alpha)), with the exponent alpha in the range 1 greater than or equal to alpha >0. The calculated values of the exponent appear to be in very good agreement with many experimental results. In particular, the theory, for the first time, explains the experimentally observed high value of alpha even at the barrier frequency, omega(b). similar or equal to 9 X 10(12) s(-1) for the isomerization reaction of 2-(2'-propenyl)anthracene in liquid eta-alkanes. The present study can also explain the reason for the very low value of vb observed in another study for the isomerization reaction of trans-stilbene in liquid n-alkanes. For omega(b) greater than or equal to 2.0 X 10(13) s(-1), we obtain alpha similar or equal to 0, which implies that the barrier crossing rate becomes identical to the transition-state theory predictions. A careful analysis of isomerization reaction dynamics involving large amplitude motion suggests that the barrier crossing dynamics itself may become irrelevant in highly viscous liquids and the rate might again be coupled directly to the viscosity. This crossover is predicted to be strongly temperature dependent and could be studied by changing the solvent viscosity by the application of pressure. (C) 1999 American Institute of Physics. [S0021-9606(9950514-X].

Solvation dynamics of a quadrupolar solute in dipolar liquids

A microscopic calculation of solvation dynamics of dipolar and quadrupolar solutes in liquid water and acetonitrile is presented. The solvation is found to he biphasic. The calculated solvation time correlation function of ionic quadrupolar solute (K+) in water is in good agreement with re cent computer simulation results. Present study reveals some interesting aspects of quadrupolar solvation dynamics which differ significantly from that of ionic and dipolar solvation.

A chemical approach to understand fragilities of glass-forming liquids

Fragility is viewed as a measure of the loss of rigidity of a glass structure above its glass transition temperature. It is attributed to the weakness of directional bonding and to the presence of a high density of low-energy configurational states. An a priori fragility function of electronegativities and bond distances is proposed which quite remarkably reproduces the entire range of reported fragilities and demonstrates that the fragility of a melt is indeed encrypted in the chemistry of the parent material. It has also been shown that the use of fragility-modified activation barriers in the Arrhenius function account for the whole gamut of viscosity behavior of liquids. It is shown that fragility can be a universal scaling parameter to collapse all viscosity curves on to a master plot.

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.

Adam-Gibbs Relation for Glass-Forming Liquids in Two, Three, and Four Dimensions

The Adam-Gibbs relation between relaxation times and the configurational entropy has been tested extensively for glass formers using experimental data and computer simulation results. Although the form of the relation contains no dependence on the spatial dimensionality in the original formulation, subsequent derivations of the Adam-Gibbs relation allow for such a possibility. We test the Adam-Gibbs relation in two, three, and four spatial dimensions using computer simulations of model glass formers. We find that the relation is valid in three and four dimensions. But in two dimensions, the relation does not hold, and interestingly, no single alternate relation describes the results for the different model systems we study.