Fickian crossover and length scales from two point functions in supercooled liquids

Particle motion of a Lennard-Jones supercooled liquid near the glass transition is studied by molecular dynamics simulations. We analyze the wave vector dependence of relaxation times in the incoherent self-scattering function and show that at least three different regimes can be identified and its scaling properties determined. The transition from one regime to another happens at characteristic length scales. The length scale associated with the onset of Fickian diffusion corresponds to the maximum size of heterogeneities in the system, and the characteristic time scale is several times larger than the alpha relaxation time. A second crossover length scale is observed, which corresponds to the typical time and length of heterogeneities, in agreement with results from four point functions. The different regimes can be traced back to the behavior of the van Hove distribution of displacements, which shows a characteristic exponential regime in the heterogeneous region before the crossover to Gaussian diffusion and should be observable in experiments. Our results show that it is possible to obtain characteristic length scales of heterogeneities through the computation of two point functions at different times.     

The short time regime of segmental dynamics of glass-forming polymers

The short time regime (10⁻¹³ s ≤ t ≤ 10⁻¹¹ s) of the segmental dynamics of several glass-forming polymers with different chemical microstructure and glass-transition temperature has been investigated by means of time of flight (TOF) neutron scattering techniques. In this paper we present a summary of the main results obtained as well as the theoretical framework proposed by us to interpret these results. The main conclusion is that the segmental dynamics of glass-forming polymers display in the very short time range (t = 2 ps) a crossover from a simple exponential behaviour towards a non-exponential regime. The first exponential decay, which is controlled by C-C rotational barriers, can be understood as a trace of the behaviour of the system in absence of the effects (correlations, cooperativity, memory effects…) which characterise the dense supercooled liquid like state against the normal liquid state. The non-exponential regime at t > 2 ps corresponds to what is usually understood as α and β relaxations. Some implications of these results are also discussed.     

Molecular view of the isothermal transformation of a stable glass to a liquid

We have used neutron reflectivity to measure translational motion on the nanometer length scale in exceptionally stable glasses of tris(naphthylbenzene). These glasses are prepared by vapor deposition onto a substrate held somewhat below the glass transition temperature (T(g) = 342 K). When the most stable samples are annealed at 345 K, no translational motion is observed on the 12 nm length scale for over 10,000 s and full mixing requires more than 60,000 s. For comparison, the equilibrium supercooled liquid mixes in 1000 s at this temperature and on this length scale. These measurements provide insight into the mechanism by which a stable glass transforms into a liquid. "Melting" of the stable glass appears to occur by the growth of liquid regions into the surrounding glassy matrix, perhaps by a surface-initiated growth process. At 345 K, translational motion in the stable glass is at least 100 times slower than motion in the supercooled liquid.     

Diffusion and Relaxation Dynamics of Supercooled Polymer Melts

The dynamic properties of polymer melts are investigated in the range of normal liquid regime to the supercooled liquid regime. The polymer is modeled as a coarse-grained bead-spring model with chain length ranging from 5 to 160. The mean squared displacement and non-Gaussian parameter are used to describe the self diffusion of polymer beads. We find slow dynamics with decreasing temperature and increasing chain length. The time evolution of non-Gaussian parameters shows two peaks(or one peak one shoulder) in the α-relaxation time, τα, regime and sub-diffusion time regime, respectively, where the first primary peak indicates the dynamic heterogeneity stemmed from the motion of beads, and the secondary peak is the result of correlated motion along a polymer chain. Moreover, the relaxation of polymer beads shows clear two-step decay in supercooled melts and the dynamics shows growing heterogeneity with decreasing temperature. As chain length is increased, a peak of the dynamic susceptibility occurs, and the peak height,χ*4, increases and then reaches a plateau. The curves of the height of the first peak of α_2, α _2~*, versus τ and the curves of χ_4~*α versus τα follow two master curves for different chain lengths. Our results indicate the similarity of dynamic heterogeneity dominated by the motion of single bead even the chain length is different. It is interesting to find that the Stokes-Einstein(SE) relation between τα and diffusion coefficient D, D~τ-1 q, is highly length-scale dependent. The SE relation breaks down in both normal melts regime and supercooled regime at large magnitude of wave vectors, attributed to the non-Brownian motion arising from the chain connectivity and growing heterogeneity due to supercooling. However, the SE relation is reconstructed when the probing length scale is large(at small magnitude of wave vectors). Our results show a hierarchical physical picture of the supercooled polymeric dynamics.     

A mode coupling theory description of the short- and long-time dynamics of nematogens in the isotropic phase

Optical heterodyne-detected optical Kerr effect (OHD-OKE) experimental data are pre-sented on nematogens 4-(trans-4-n-octylcyclohexyl)isothiocyanatobenzene (8-CHBT), and 4-(4-pentyl-cyclohexyl)-benzonitrile (5-PCH) in the isotropic phase. The 8-CHBT and 5-PCH data and previously published data on 4-pentyl-4-biphenylcarbonitrile (5-CB) are analyzed using a modification of a schematic mode coupling theory (MCT) that has been successful in describing the dynamics of supercooled liquids. At long time, the OHD-OKE data (orientational relaxation) are well described with the standard Landau-de Gennes (LdG) theory. The data decay as a single exponential. The decay time diverges as the isotropic to nematic phase transition is approached from above. Previously there has been no theory that can describe the complex dynamics that occur at times short compared to the LdG exponential decay. Earlier, it has been noted that the short-time nematogen dynamics, which consist of several power laws, have a functional form identical to that observed for the short time behavior of the orientational relaxation of supercooled liquids. The temperature-dependent orientational dynamics of supercooled liquids have recently been successfully described using a schematic mode coupling theory. The schematic MCT theory that fits the supercooled liquid data does not reproduce the nematogen data within experimental error. The similarities of the nematogen data to the supercooled liquid data are the motivation for applying a modification of the successful MCT theory to nematogen dynamics in the isotropic phase. The results presented below show that the new schematic MCT theory does an excellent job of reproducing the nematogen isotropic phase OHD-OKE data on all time scales and at all temperatures.     

A phenomenological approach to relaxation in disordered systems

We propose a phenomenological approach to relaxation in disordered systems which is modelled after the well-known percolation behaviour. The approach shows an interesting behaviour of the time correlation function of relevant variables which is characterized by three temporal regimes, an initial exponential decay, a power-law regime and a stretched exponential decay for long times. We give some relevant examples of this behaviour in a wide variety of systems such as gels, microemulsions and glasses.     

Liquid-liquid transition in supercooled water investigated by interaction with LiCl and Xe

The hypothesis that supercooled water consists of two distinct liquid phases has been explored on the basis of their ability to hydrate nonpolar (Xe) and electrolytic (LiCl) species. Xe incorporated in the bulk of amorphous solid water survives in the deeply supercooled regime above the glass-transition temperature of 136 K and is finally dehydrated at 165 K, whereas LiCl dissolves only in the liquid phase appearing above 165 K. The second liquid phase connects with normal water as inferred from high (poor) solubility of LiCl (Xe). This result also suggests that decoupling of translational diffusion and viscosity in the deeply supercooled regime is caused by domain structures of the two liquid phases formed during a possible liquid-liquid transition.     

Deformation of inherent structures to detect long-range correlations in supercooled liquids

We propose deformations of inherent structures as a suitable tool for detecting structural changes underlying the onset of cooperativity in supercooled liquids. The non-affine displacement (NAD) field resulting from the applied deformation shows characteristic differences between the high temperature liquid and supercooled state, which are typically observed in dynamic quantities. The average magnitude of the NAD is very sensitive to temperature changes in the supercooled regime and is found to be strongly correlated with the inherent structure energy. In addition, the NAD field is characterized by a correlation length that increases upon lowering the temperature towards the supercooled regime.     

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

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

Superdipole liquid scenario for the dielectric primary relaxation in supercooled polar liquids

We propose a dynamic structure of coupled dynamic molecular strings for supercooled small polar molecule liquids and accordingly we obtain the Hamiltonian of the rotational degrees of freedom of the system. From the Hamiltonian, the strongly correlated supercooled polar liquid state is renormalized to a normal superdipole liquid state. This scenario describes the following main features of the primary or alpha-relaxation dynamics in supercooled polar liquids: (1) the average relaxation time evolves from a high temperature Arrhenius to a low temperature non-Arrhenius or super-Arrhenius behavior; (2) the relaxation function crosses over from the high temperature exponential to low temperature nonexponential form; and (3) the temperature dependence of the relaxation strength shows non-Curie features. According to the present model, the crossover phenomena of the first two characteristics arise from the transition between the superdipole gas and the superdipole liquid. The model predictions are quantitatively compared with the experimental results of glycerol, a typical glass former.     

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.     

Slow dynamics of the high density Gaussian core model

We numerically study crystal nucleation and glassy slow dynamics of the one-component Gaussian core model (GCM) at high densities. The nucleation rate at a fixed supercooling is found to decrease as the density increases. At very high densities, the nucleation is not observed at all in the time window accessed by long molecular dynamics (MD) simulation. Concomitantly, the system exhibits typical slow dynamics of the supercooled fluids near the glass transition point. We compare the simulation results of the supercooled GCM with the predictions of mode-coupling theory (MCT) and find that the agreement between them is better than any other model glassformers studied numerically in the past. Furthermore, we find that a violation of the Stokes-Einstein relation is weaker and the non-Gaussian parameter is smaller than canonical glassformers. Analysis of the probability distribution of the particle displacement clearly reveals that the hopping effect is strongly suppressed in the high density GCM. We conclude from these observations that the GCM is more amenable to the mean-field picture of the glass transition than other models. This is attributed to the long-ranged nature of the interaction potential of the GCM in the high density regime. Finally, the intermediate scattering function at small wavevectors is found to decay much faster than its self part, indicating that dynamics of the large-scale density fluctuations decouples with the shorter-ranged caging motion.     

Dynamics in supercooled ionic organic liquids and mode coupling theory analysis

Optically heterodyne-detected optical Kerr effect experiments are applied to study the orientational dynamics of the supercooled ionic organic liquids N-propyl-3-methylpyridinium bis(trifluoromethylsulfonyl)imide (PMPIm) and 1-ethyl-3-methylimidazolium tosylate (EMImTOS). The orientational dynamics are complex with relaxation involving several power law decays followed by a final exponential decay. A mode coupling theory (MCT) schematic model, the Sj?gren model, was able to reproduce the PMPIm data very successfully over a wide range of times from 1 ps to hundreds of ns for all temperatures studied. Over the temperature range from room temperature down to the critical temperature Tc of 231 K, the OHD-OKE signal of PMPIm is characterized by the intermediate power law t(-1.00+/-0.04) at short times, a von Schweidler power law t(-0.51+/-0.03) at intermediate times, and a highly temperature-dependent exponential (alpha relaxation) at long times. This form of the decay is identical to the form observed previously for a large number of organic van der Waals liquids. MCT analysis indicates that the theory can explain the experimental data very well for a range of temperatures above Tc, but as might be expected, there are some deviations from the theoretical modeling at temperatures close to Tc. For EMImTOS, the orientational dynamics were studied on the ps time scale in the deeply supercooled region near its glass transition temperature. The orientational relaxation of EMImTOS clearly displays the feature associated with the boson peak at approximately 2 ps, which is the first time domain evidence of the boson peak in ionic organic liquids. Overall, all the dynamical features observed earlier for organic van der Waals liquids using the same experimental technique are also observed for organic ionic liquids.     

On melting dynamics and the glass transition. I. Glassy aspects of melting dynamics

The following properties are in the present literature associated with the behavior of supercooled glass-forming liquids: faster than exponential growth of the relaxation time, dynamical heterogeneities, growing point-to-set correlation length, crossover from mean-field behavior to activated dynamics. In this paper we argue that these properties are also present in a much simpler situation, namely the melting of the bulk of an ordered phase beyond a first order phase transition point. This is a promising path toward a better theoretical, numerical and experimental understanding of the above phenomena and of the physics of supercooled liquids. We discuss in detail the analogies and the differences between the glass and the bulk melting transitions.     

Polymer brushes in cylindrical pores: simulation versus scaling theory

The structure of flexible polymers endgrafted in cylindrical pores of diameter D is studied as a function of chain length N and grafting density sigma, assuming good solvent conditions. A phenomenological scaling theory, describing the variation of the linear dimensions of the chains with sigma, is developed and tested by molecular dynamics simulations of a bead-spring model. Different regimes are identified, depending on the ratio of D to the size of a free polymer N(3/5). For D>N(3/5) a crossover occurs for sigma=sigma*=N(-6/5) from the "mushroom" behavior (R(gx)=R(gy)=R(gz)=N(35)) to the behavior of a flat brush (R(gz)=sigma(1/3)N,R(gx)=R(gy)=sigma(-1/12)N(1/2)), until at sigma**=(D/N)3 a crossover to a compressed state of the brush, [R(gz)=D,R(gx)=R(gy)=(N(3)D/4sigma)(1/8)相似文献   

Thermodynamic, structural, and dynamic properties of supercooled water confined in mesoporous MCM-41 studied with calorimetric, neutron diffraction, and neutron spin echo measurements

Thermodynamic, structural, and dynamic properties of heavy water (D(2)O) confined in mesoporous silica glass MCM-41 C10, C12, and C14 were investigated by differential scanning calorimetry, neutron diffraction, and neutron spin echo (NSE) measurements, respectively. The DSC data showed that no crystallization of D(2)O confined in C10 occurs in a temperature range between 298 and 180 K, and that crystalline ice is formed at 204 and 221 K for C12 and C14, respectively. For C10, the neutron radial distribution functions of confined D(2)O suggested a structural change in the supercooled state between 223 and 173 K. For C10 sample, it has been found that the tetrahedral-like water structure is partially enhanced in the central part of pores at 173 K. For all the samples, the intermediate scattering functions from the NSE measurements are fitted by the Kohlrausch-Williams-Watts stretched exponential function which implies that confined supercooled D(2)O exhibits a wide distribution of relaxation times. For C10, C12, and C14 samples, between 298 and 240 K, the relaxation times of supercooled D(2)O follow remarkably well the Vogel-Fulcher-Tamman equation; for C10 sample, below 240 K, the relaxation times of nonfreezing D(2)O show an Arrhenius type behavior. From the present experimental results on calorimetric, structural, and dynamic properties, it has been concluded that supercooled D(2)O confined in MCM-41 C10 experiences a transition from high-density to low-density hydrogen-bonded structure at around 229 K.     

Decay of correlation functions in hard-sphere mixtures: structural crossover

We investigate the decay of pair correlation functions in a homogeneous (bulk) binary mixture of hard spheres. At a given state point the asymptotic decay r-->infinity of all three correlation functions is governed by a common exponential decay length and a common wavelength of oscillations. Provided the mixture is sufficiently asymmetric, size ratios q less than or approximately 0.7, we find that the common wavelength reflects either the size of the small or that of the big spheres. By analyzing the (complex) poles of the partial structure factors we find a sharp structural crossover line in the phase diagram. On one side of this line the common wavelength is approximately the diameter of the smaller sized spheres whereas on the other side it is approximately the diameter of the bigger ones; the wavelength of the longest ranged oscillations changes discontinuously at the structural crossover line. Using density functional theory and Monte Carlo simulations we show that structural crossover also manifests itself in the intermediate range behavior of the pair correlation functions and we comment on the relevance of this observation for real (colloidal) mixtures. In highly asymmetric mixtures, q< or =0.1, where there is metastable fluid-fluid transition, we find a Fisher-Widom line with two branches. This line separates a region of the phase diagram where the decay of pair correlations is oscillatory from one in which it is monotonic.     

Collective dynamics of supercooled water close to the liquid-liquid coexistence lines

We report molecular dynamics simulation results for the collective dynamical properties of supercooled bulk water at 180 K at three different densities, corresponding to different phases whose coexistence has recently been discovered in the supercooled regime. In this study, we focus on the behaviour of the longitudinal and transverse current correlation functions and their relative spectra, which we analyze in detail to understand the dynamical processes responsible for the main features observed. Despite the considerable differences in the structure and densities of the three thermodynamic states considered, the obtained current correlation functions show rather similar behaviour in every case. We show that the longitudinal spectra can only be described in terms of three Lorentzian lines, while the accurate reproduction of the transverse spectra requires at least four separate spectral lines. In fact, the behaviour of the peak frequency of the modes necessary to reproduce the spectra as a function of the wavevector indicates in a clear way the nature of the dynamical process. We demonstrate the presence of collective modes associated with the propagation of both longitudinal and transverse sound along with the important contribution of "optical-like" modes, which point out the relevance of localized motions for a right interpretation of the spectral line shape. The wavevector dependence of the relative contributions of the various modes to the total spectral area is also discussed.