The new insight into dynamic crossover in glass forming liquids from the apparent enthalpy analysis

One of the most intriguing phenomena in glass forming systems is the dynamic crossover (T(B)), occurring well above the glass temperature (T(g)). So far, it was estimated mainly from the linearized derivative analysis of the primary relaxation time τ(T) or viscosity η(T) experimental data, originally proposed by Stickel et al. [J. Chem. Phys. 104, 2043 (1996); J. Chem. Phys. 107, 1086 (1997)]. However, this formal procedure is based on the general validity of the Vogel-Fulcher-Tammann equation, which has been strongly questioned recently [T. Hecksher et al. Nature Phys. 4, 737 (2008); P. Lunkenheimer et al. Phys. Rev. E 81, 051504 (2010); J. C. Martinez-Garcia et al. J. Chem. Phys. 134, 024512 (2011)]. We present a qualitatively new way to identify the dynamic crossover based on the apparent enthalpy space (H(a) (')=dlnτ/d(1/T)) analysis via a new plot lnH(a) (') vs. 1∕T supported by the Savitzky-Golay filtering procedure for getting an insight into the noise-distorted high order derivatives. It is shown that depending on the ratio between the "virtual" fragility in the high temperature dynamic domain (m(high)) and the "real" fragility at T(g) (the low temperature dynamic domain, m = m(low)) glass formers can be splitted into two groups related to f < 1 and f > 1, (f = m(high)∕m(low)). The link of this phenomenon to the ratio between the apparent enthalpy and activation energy as well as the behavior of the configurational entropy is indicated.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Is a "homogeneous" description of dynamic heterogeneities possible?

We study the simplest model of dynamic heterogeneities in glass forming liquids: one-spin facilitated kinetic Ising model introduced by Fredrickson and Andersen [G. H. Fredrickson and H. C. Andersen, Phys. Rev. Lett. 53, 1244 (1984); J. Chem. Phys. 83, 5822 (1985)]. We show that the low-temperature, long-time behavior of the density autocorrelation function predicted by a scaling approach can be obtained from a self-consistent mode-couplinglike approximation.     

Thermodynamic interpretation of the scaling of the dynamics of supercooled liquids

The recently discovered scaling law for the relaxation times, tau(T,upsilon) = I(Tupsilon(gamma)), where T is temperature and upsilon the specific volume, is derived by a revision of the entropy model of the glass transition dynamics originally proposed by Avramov [J. Non-Cryst. Solids 262, 258 (2000)]. In this modification the entropy is calculated by an alternative route. The resulting expression for the variation of the relaxation time with T and upsilon is shown to accurately fit experimental data for several glass-forming liquids and polymers over an extended range encompassing the dynamic crossover. From this analysis, which is valid for any model in which the relaxation time is a function of the entropy, we find that the scaling exponent gamma can be identified with the Gruneisen constant.     

On the effects of ion-wall chemical association on the electric double layer: a density functional approach for the restricted primitive model at a charged wall

We describe a density functional theory for the restricted primitive model of ionic fluid at a charged wall with active sites to which ions can bond. The theory is an extension of our recent approach [Pizio et al., J. Chem. Phys. 121, 11957 (2004)] and is focused in the effects of specific adsorption of ions on the wall, besides the electrostatic phenomena. In order to solve the problem, we use the first-order thermodynamic perturbation theory of chemical association developed by Wertheim [J. Chem. Phys. 87, 7323 (1987)]. The microscopic structure of the electric double layer and the amount of adsorbed charge are investigated. Also, the temperature dependence of capacitance is analyzed. The capacitance depends on the kind of ions that form associative bonds with the surface sites and is determined by a net charge acting on the diffuse layer. The shape of the temperature dependence of capacitance essentially depends on the association energy and the density of bonding sites.     

Recovering the Crooks equation for dynamical systems in the isothermal-isobaric ensemble: a strategy based on the equations of motion

The Crooks equation [Eq. (10) in J. Stat. Phys. 90, 1481 (1998)] relates the work done on a system during a nonequilibrium transformation to the free energy difference between the final and the initial state of the transformation. Recently, the authors have derived the Crooks equation for systems in the canonical ensemble thermostatted by the Nose-Hoover or Nose-Hoover chain method [P. Procacci et al., J. Chem. Phys. 125, 164101 (2006)]. That proof is essentially based on the fluctuation theorem by Evans and Searles [Adv. Phys. 51, 1529 (2002)] and on the equations of motion. Following an analogous approach, the authors derive here the Crooks equation in the context of molecular dynamics simulations of systems in the isothermal-isobaric (NPT) ensemble, whose dynamics is regulated by the Martyna-Tobias-Klein algorithm [J. Chem. Phys. 101, 4177 (1994)]. Their present derivation of the Crooks equation correlates to the demonstration of the Jarzynski identity for NPT systems recently proposed by Cuendet [J. Chem. Phys. 125, 144109 (2006)].     

Ultrafast far-infrared dynamics probed by terahertz pulses: a frequency domain approach. I. Model systems

Time-resolved terahertz spectroscopy has become a widely used experimental tool for the investigation of ultrafast dynamics of polar systems in the far infrared. We have recently proposed an analytical method for the extraction of a transient two-dimensional susceptibility from the experimental data [N?mec, Kadlec, and Kuzel, J. Chem. Phys. 117, 8454 (2002)]. In the present paper the methodology of optical pump-terahertz probe experiments is further developed for direct application in realistic experimental situations. The expected two-dimensional transient response function is calculated for a number of model cases (including Drude dynamics of free carriers, harmonic and anharmonic oscillator modes); these results serve as a basis for the interpretation of experimental results. We discuss also the cases where only partial (one-dimensional) information about the system dynamics can be experimentally obtained.     

Affinity and its derivatives in the glass transition process

The thermodynamic treatment of the glass transition remains an issue of intense debate. When associated with the formalism of non-equilibrium thermodynamics, the lattice-hole theory of liquids can provide new insight in this direction, as has been shown by Schmelzer and Gutzow [J. Chem. Phys. 125, 184511 (2006)], by Mo?ller et al. [J. Chem. Phys. 125, 094505 (2006)], and more recently by Tropin et al. [J. Non-Cryst. Solids 357, 1291 (2011); ibid. 357, 1303 (2011)]. Here, we employ a similar approach. We include pressure as an additional variable, in order to account for the freezing-in of structural degrees of freedom upon pressure increase. Second, we demonstrate that important terms concerning first order derivatives of the affinity-driving-force with respect to temperature and pressure have been previously neglected. We show that these are of crucial importance in the approach. Macroscopic non-equilibrium thermodynamics is used to enlighten these contributions in the derivation of C(p),κ(T), and α(p). The coefficients are calculated as a function of pressure and temperature following different theoretical protocols, revealing classical aspects of vitrification and structural recovery processes. Finally, we demonstrate that a simple minimalist model such as the lattice-hole theory of liquids, when being associated with rigorous use of macroscopic non-equilibrium thermodynamics, is able to account for the primary features of the glass transition phenomenology. Notwithstanding its simplicity and its limits, this approach can be used as a very pedagogical tool to provide a physical understanding on the underlying thermodynamics which governs the glass transition process.     

Complex wave patterns in an effective reaction-diffusion model for chemical reactions in microemulsions

An effective medium theory is employed to derive a simple qualitative model of a pattern forming chemical reaction in a microemulsion. This spatially heterogeneous system is composed of water nanodroplets randomly distributed in oil. While some steps of the reaction are performed only inside the droplets, the transport through the extended medium occurs by diffusion of intermediate chemical reactants as well as by collisions of the droplets. We start to model the system with heterogeneous reaction-diffusion equations and then derive an equivalent effective spatially homogeneous reaction-diffusion model by using earlier results on homogenization in heterogeneous reaction-diffusion systems [S.Alonso, M.Ba?r, and R.Kapral, J. Chem. Phys. 134, 214102 (2009)]. We study the linear stability of the spatially homogeneous state in the resulting effective model and obtain a phase diagram of pattern formation, that is qualitatively similar to earlier experimental results for the Belousov-Zhabotinsky reaction in an aerosol OT (AOT)-water-in-oil microemulsion [V.K.Vanag and I.R.Epstein, Phys. Rev. Lett. 87, 228301 (2001)]. Moreover, we reproduce many patterns that have been observed in experiments with the Belousov-Zhabotinsky reaction in an AOT oil-in-water microemulsion by direct numerical simulations.     

One-dimensional model for water and aqueous solutions. I. Pure liquid water

Two simplified one-dimensional models for waterlike particles are studied. One is referred to as the primitive model which is a simplified version of a model introduced by Ben-Naim in 1992 [Statistical Thermodynamics for Chemists and Biochemists (Plenum, New York, 1992)]. The second, referred to as the primitive cluster model, is a simplified version of the model used by Lovett and Ben-Naim in 1969 [J. Chem. Phys. 51, 3108 (1969)]. The two models are shown to be nearly equivalent and both exhibit some of the most characteristic behavior of liquid water. It is argued that a key feature of the molecular interactions--the correlation between the strong binding energy and low local density--is essential for the manifestation of the anomalous behavior of liquid water. It is also essential for the understanding of the outstanding behavior of liquid water.     

Phase behavior of the restricted primitive model of ionic fluids with association in slitlike pores. Density-functional approach

We present results of investigations of condensation of restricted primitive model of electrolyte solutions with association between oppositely charged ions confined to slitlike pores. The associative interaction leads to the formation of ionic pairs. It is accounted for by incorporating the first-order thermodynamic perturbation theory into the free energy functional. In order to elucidate the role of association, the phase diagrams are compared with those obtained by us recently [O. Pizio et al., J. Chem. Phys. 121, 11957 (2004)] for the restricted primitive model. The inclusion of the association into the theory leads to lowering the critical temperature for the fluid confined to pores with uncharged and with charged walls. We have observed that the average fraction of bonded ions is high along the coexistence envelope.     

Dynamic coupling between coordinates in a model for biomolecular isomerization

To understand a complex reaction, it is necessary to project the dynamics of the system onto a low-dimensional subspace of physically meaningful coordinates. We recently introduced an automatic method for identifying coordinates that relate closely to stable-state commitment probabilities and successfully applied it to a model for biomolecular isomerization, the C(7eq)-->alpha(R) transition of the alanine dipeptide [A. Ma and A. R. Dinner, J. Phys. Chem. B 109, 6769 (2005)]. Here, we explore approximate means for estimating diffusion tensors for systems subject to restraints in one and two dimensions and then use the results together with an extension of Kramers theory for unimolecular reaction rates [A. Berezhkovskii and A. Szabo, J. Chem. Phys. 122, 014503 (2005)] to show explicitly that both the potential of mean force and the diffusion tensor are essential for describing the dynamics of the alanine dipeptide quantitatively. In particular, the signficance of off-diagonal elements of the diffusion tensor suggests that the coordinates of interest are coupled by the hydrodynamic-like response of the bath of remaining degrees of freedom.     

Temperature of systems out of thermodynamic equilibrium

Two phenomenological approaches are currently used in the study of the vitreous state. One is based on the concept of fictive temperature introduced by Tool [J. Res. Natl. Bur. Stand. 34, 199 (1945)] and recently revisited by Nieuwenhuizen [Phys. Rev. Lett. 80, 5580 (1998)]. The other is based on the thermodynamics of irreversible processes initiated by De Donder at the beginning of the last century [L'Affinite (Gauthier-Villars, Paris, 1927)] and recently used by Moller et al. for a thorough study of the glass transition [J. Chem. Phys. 125, 094505 (2006)]. This latter approach leads to the possibility of describing the glass transition by means of the freezing-in of one or more order parameters connected to the internal structural degrees of freedom involved in the vitrification process. In this paper, the equivalence of the two preceding approaches is demonstrated, not only for glasses but in a very general way for any system undergoing an irreversible transformation. This equivalence allows the definition of an effective temperature for all systems departed from equilibrium generating a positive amount of entropy. In fact, the initial fictive temperature concept of Tool leads to the generalization of the notion of temperature for systems out of thermodynamic equilibrium, for which glasses are just particular cases.     

Using the thermal Gaussian approximation for the Boltzmann operator in semiclassical initial value time correlation functions

The thermal Gaussian approximation (TGA) recently developed by Frantsuzov et al. [Chem. Phys. Lett. 381, 117 (2003)] has been demonstrated to be a practical way for approximating the Boltzmann operator exp(-betaH) for multidimensional systems. In this paper the TGA is combined with semiclassical (SC) initial value representations (IVRs) for thermal time correlation functions. Specifically, it is used with the linearized SC-IVR (LSC-IVR, equivalent to the classical Wigner model), and the "forward-backward semiclassical dynamics" approximation developed by Shao and Makri [J. Phys. Chem. A 103, 7753 (1999); 103, 9749 (1999)]. Use of the TGA with both of these approximate SC-IVRs allows the oscillatory part of the IVR to be integrated out explicitly, providing an extremely simple result that is readily applicable to large molecular systems. Calculation of the force-force autocorrelation for a strongly anharmonic oscillator demonstrates its accuracy, and calculation of the velocity autocorrelation function (and thus the diffusion coefficient) of liquid neon demonstrates its applicability.     

Brownian dynamics algorithm for bead-rod semiflexible chain with anisotropic friction

A model of semiflexible bead-rod chain with anisotropic friction can mimic closely the hydrodynamics of a slender filament. We present an efficient algorithm for Brownian dynamics simulations of this model with configuration dependent anisotropic bead friction coefficients. The algorithm is an extension of that given previously for the case of configuration independent isotropic friction coefficients by Grassia and Hinch [J. Fluid Mech. 308, 255 (1996)]. We confirm that the algorithm yields predicted values for various equilibrium properties. We also present a stochastic algorithm for evaluation of the stress tensor, and we show that in the limit of stiff chains the algorithm recovers the results of Kirkwood and Plock [J. Chem. Phys. 24, 665 (1956)] for rigid rods with hydrodynamic interactions.     

Molecular dynamics simulation study of nonconcatenated ring polymers in a melt. I. Statics

Molecular dynamics simulations were conducted to investigate the structural properties of melts of nonconcatenated ring polymers and compared to melts of linear polymers. The longest rings were composed of N = 1600 monomers per chain which corresponds to roughly 57 entanglement lengths for comparable linear polymers. For the rings, the radius of gyration squared, [linear span]R(g)(2)[linear span], was found to scale as N(4/5) for an intermediate regime and N(2/3) for the larger rings indicating an overall conformation of a crumpled globule. However, almost all beads of the rings are "surface beads" interacting with beads of other rings, a result also in agreement with a primitive path analysis performed in the next paper [J. D. Halverson, W. Lee, G. S. Grest, A. Y. Grosberg, and K. Kremer, J. Chem. Phys. 134, 204905 (2011)]. Details of the internal conformational properties of the ring and linear polymers as well as their packing are analyzed and compared to current theoretical models.     

Configurational entropy of binary hard-disk glasses: nonexistence of an ideal glass transition

We study the thermodynamics of a binary hard-disk mixture in which the ratio of disk diameters is kappa=1.4. We use a recently developed molecular dynamics algorithm to calculate the free-volume entropy of glassy configurations and obtain the configurational entropy (degeneracy) of the supercompressed liquid as a function of density. We find that the configurational entropy of the glasses near the kinetic glass transition is very close to the mixing entropy, suggesting that the degeneracy is zero only for the phase-separated crystal. We explicitly construct an exponential number of jammed packings with densities spanning the spectrum from the accepted "amorphous" glassy state to the phase-separated crystal, thus showing that there is no ideal glass transition in binary hard-disk mixtures. This construction also demonstrates that the ideal glass, defined as having zero configurational entropy, is not amorphous, but instead is nothing more than a phase-separated crystal. This critique of the presumed existence of an ideal glass parallels our previous critique of the idea that there is a most-dense random (close) packing for hard spheres [Torquato et al., Phys. Rev. Lett. 84, 2064 (2000)]. We also perform free-energy calculations to determine the equilibrium phase behavior of the system. The calculations predict a first-order freezing transition at a density below the kinetic glass transition. However, this transition appears to be strongly kinetically suppressed and is not observed directly. New simulation techniques are needed in order to gain a more complete understanding of the thermodynamic and kinetic behavior of the binary disk mixture and, in particular, of the demixing process during crystallization.