Molecular dynamics simulation of the fragile glass-former orthoterphenyl: A flexible molecule model

Does the Naı̈ve Mode-Coupling Power Law Divergence Provide an Objective Determination of the Crossover Temperature in Glass Formation Behavior?

The glass formation temperature range is commonly divided into a weakly supercooled regime at higher temperatures and a deeply supercooled regime at lower temperatures, with a change in the physical mechanisms that govern dynamics often postulated to occur at the crossover between these regimes. This crossover temperature Tc is widely determined based on a fit of relaxation time vs temperature data to a power law divergence form predicted by the naı̈ve mode coupling theory (MCT). Here, we show, based on simulation data spanning polymeric, small molecule organic, metallic, and inorganic glass formers, that this approach does not yield an objective measure of a crossover temperature. Instead, the value of Tc is determined by the lowest temperature Tmin employed in the fit, and no regime of stationary or convergent Tc value is generally observed as Tmin is varied. Nor does the coefficient of determination R2 provide any robust means of selecting a fit range and thus a value of Tc. These results may require a re-evaluation of published results that have employed the fit MCT Tc value as a metric of temperature-dependent dynamics or a benchmark for depth of supercooling, and they highlight a need for the field to converge on a more objective determination of any posited crossover temperature between high and low temperature regimes of glass formation.

Construction of a quantitative relation between structural relaxation and dynamic heterogeneity by vibrational dynamics in glass-forming liquids and polymers

The structural relaxation slows down drastically upon approaching the glass transition, accompanied by the significant growth of dynamic heterogeneity. The fundamental question of elusiveness and interest is whether there exists an underlying quantitative relationship between structural relaxation and dynamic heterogeneity. Here, we reveal that b̃ which is related to the reduced mean square displacements to overcome the energy barriers of activated jumps, instead of the kinetic fragility m, is the genuine key parameter connecting dynamic heterogeneity with structural relaxation for varying types of glass formers. Furthermore, based on the dependence of dynamic heterogeneity on the Debye-Waller factor we obtained a direct quantitative relation between dynamic heterogeneity and structural relaxation is built for different glass-forming liquids. More importantly, a scaling collapse of structural relaxation and dynamic heterogeneity is achieved by the important parameter b̃. These results are of fundamental and critical importance for developing a unified theory of glassy dynamics.

Experimental evidence of mosaic structure in strongly supercooled molecular liquids

When a liquid is cooled to produce a glass its dynamics, dominated by the structural relaxation, become very slow, and at the glass-transition temperature T_g its characteristic relaxation time is about 100 s. At slightly elevated temperatures (~1.2 T_g) however, a second process known as the Johari-Goldstein relaxation, β_JG, decouples from the structural one and remains much faster than it down to T_g. While it is known that the β_JG-process is strongly coupled to the structural relaxation, its dedicated role in the glass-transition remains under debate. Here we use an experimental technique that permits us to investigate the spatial and temporal properties of the β_JG relaxation, and give evidence that the molecules participating in it are highly mobile and spatially connected in a system-spanning, percolating cluster. This correlation of structural and dynamical properties provides strong experimental support for a picture, drawn from theoretical studies, of an intermittent mosaic structure in the deeply supercooled liquid phase.

The Johari-Goldstein relaxation is a precursor of the glass transition in liquids. Caporaletti et al. use time-dependent interferometry data to substantiate its suggested structural appearance as a globally percolating, fluctuating mosaic.

Insight into the Structure and Dynamics of Polymers by Neutron Scattering Combined with Atomistic Molecular Dynamics Simulations

Combining neutron scattering and fully atomistic molecular dynamics simulations allows unraveling structural and dynamical features of polymer melts at different length scales, mainly in the intermolecular and monomeric range. Here we present the methodology developed by us and the results of its application during the last years in a variety of polymers. This methodology is based on two pillars: (i) both techniques cover approximately the same length and time scales and (ii) the classical van Hove formalism allows easily calculating the magnitudes measured by neutron scattering from the simulated atomic trajectories. By direct comparison with experimental results, the simulated cell is validated. Thereafter, the information of the simulations can be exploited, calculating magnitudes that are experimentally inaccessible or extending the parameters range beyond the experimental capabilities. We show how detailed microscopic insight on structural features and dynamical processes of various kinds has been gained in polymeric systems with different degrees of complexity, and how intriguing questions as the collective behavior at intermediate length scales have been faced.

Universal localization transition accompanying glass formation: insights from efficient molecular dynamics simulations of diverse supercooled liquids

The origin of the precipitous dynamic arrest known as the glass transition is a grand open question of soft condensed matter physics. It has long been suspected that this transition is driven by an onset of particle localization and associated emergence of a glassy modulus. However, progress towards an accepted understanding of glass formation has been impeded by an inability to obtain data sufficient in chemical diversity, relaxation timescales, and spatial and temporal resolution to validate or falsify proposed theories for its physics. Here we first describe a strategy enabling facile high-throughput simulation of glass-forming liquids to nearly unprecedented relaxation times. We then perform simulations of 51 glass-forming liquids, spanning polymers, small organic molecules, inorganics, and metallic glass-formers, with longest relaxation times exceeding one microsecond. Results identify a universal particle-localization transition accompanying glass formation across all classes of glass-forming liquid. The onset temperature of non-Arrhenius dynamics is found to serve as a normalizing condition leading to a master collapse of localization data. This transition exhibits a non-universal relationship with dynamic arrest, suggesting that the nonuniversality of supercooled liquid dynamics enters via the dependence of relaxation times on local cage scale. These results suggest that a universal particle-localization transition may underpin the glass transition, and they emphasize the potential for recent theoretical developments connecting relaxation to localization and emergent elasticity to finally explain the origin of this phenomenon. More broadly, the capacity for high-throughput prediction of glass formation behavior may open the door to computational inverse design of glass-forming materials.

Ultrafast Pump-Repump-Probe Photochemical Hole Burning as a Probe of Excited-State Reaction Pathway Branching

We demonstrate pump-repump-probe (PRP) transient hole burning as a spectroscopic tool for differentiating reactive from nonreactive deactivation of excited photochemical reactants observed by transient absorption spectroscopy (TAS). This method utilizes a time-delayed, wavelength-tunable ultrafast pulse to alter the excited reactant population, with the impact of "repumping" quantified through depletions in photoproduct absorption. We apply this approach to characterize dynamics affecting the nonadiabatic photocyclization efficiency to form S₀ dihydrotriphenylene (DHT) following 266 nm excitation of ortho-terphenyl (OTP). TAS studies revealed bimodal deactivation of OTP*, but neither relaxation time scale (700 fs and 3.0 ps) could be assigned unambiguously to DHT formation due to overlap of excited-state and product spectra. PRP studies reveal that S₁ OTP only cyclizes on the slower of these time scales, with the faster process attributable to nonreactive deactivation. We demonstrate that this method offers greater photochemical insights without assuming models to globally fit spectral transients collected by TAS.

Investigation of the dynamics of aqueous proline solutions using neutron scattering and molecular dynamics simulations

We applied quasielastic neutron scattering (QENS) techniques to samples with two different contrasts (deuterated solute/hydrogenated solvent and the opposite label) to selectively study the component dynamics of proline/water solutions. Results on diluted and concentrated solutions (31 and 6 water molecules/proline molecule, respectively) were analyzed in terms of the susceptibility and considering a recently proposed model for water dynamics [Arbe et al., Phys. Rev. Lett., 2016, 117, 185501] which includes vibrations and the convolution of localized motions and diffusion. We found that proline molecules not only reduce the average diffusion coefficient of water but also extend the time/frequency range of the crossover region ('cage') between the vibrations and purely diffusive behavior. For the high proline concentration we also found experimental evidence of water heterogeneous dynamics and a distribution of diffusion coefficients. Complementary molecular dynamics simulations show that water molecules start to perform rotational diffusion when they escape the cage regime but before the purely diffusive behavior is established. The rotational diffusion regime is also retarded by the presence of proline molecules. On the other hand, a strong coupling between proline and water diffusive dynamics which persists with decreasing temperature is directly observed using QENS. Not only are the temperature dependences of the diffusion coefficients of both components the same, but their absolute values also approach each other with increasing proline concentration. We compared our results with those reported using other techniques, in particular using dielectric spectroscopy (DS). A simple approach based on molecular hydrodynamics and a molecular treatment of DS allows rationalizing the a priori puzzling inconsistency between QENS and dielectric results regarding the dynamic coupling of the two components. The interpretation proposed is based on general grounds and therefore should be applicable to other biomolecular solutions.

Are polymers standard glass-forming systems? The role of intramolecular barriers on the glass-transition phenomena of glass-forming polymers

Traditionally, polymer melts have been considered archetypal glass-formers. This has been mainly due to the fact that these systems can easily be obtained as glasses by cooling from the melt, even at low cooling rates. However, the macromolecules, i.e. the structural units of polymer systems in general, are rather different from the standard molecules. They are long objects ('chains') made by repetition of a given chemical motif (monomer) and have intra-macromolecular barriers that limit their flexibility. The influence of these properties on, for instance, the glass-transition temperature of polymers, is a topic that has been widely studied by the polymer community almost from the early times of polymer science. However, in the framework of the glass-community, the relevant influence of intra-macromolecular barriers and chain connectivity on glass-transition phenomena of polymers has started to be recognized only recently. The aim of this review is to give an overview and to critically revise the results reported on this topic over the last years. From these results, it seems to be evident that there are two different mechanisms involved in the dynamic arrest in glass-forming polymers: (i) the intermolecular packing effects, which dominate the dynamic arrest of low molecular weight glass-forming systems; and (ii) the effect of intra-macromolecular barriers combined with chain connectivity. It has also been shown that the mode coupling theory (MCT) is a suitable theoretical framework to discuss these questions. The values found for polymers for the central MCT parameter--the so-called λ-exponent--are of the order of 0.9, clearly higher than the standard values (λ ≈ 0.7) found in systems where the dynamic arrest is mainly driven by packing effects ('standard' glass-formers). Within the MCT, this is a signature of the presence of two competing mechanisms of dynamic arrest, as it has been observed in short-ranged attractive colloids or two component mixtures with dynamic asymmetry. Moreover, recent MD-simulations of a 'bead-spring' polymer model, but including intra-macromolecular potential of different strengths, confirm that the high λ-values found in polymers are due to the effect of intra-macromolecular barriers. Although there are still open questions, these results allow to conclude that there is a fundamental difference between the nature of the glass transition in polymers and in simple (standard) glass-formers.

Simulated glass-forming polymer melts: dynamic scattering functions, chain length effects, and mode-coupling theory analysis

We present molecular-dynamics simulations for a fully flexible model of polymer melts with different chain length N ranging from short oligomers (N = 4) to values near the entanglement length (N = 64). For these systems we explore the structural relaxation of the supercooled melt near the critical temperature T c of mode-coupling theory (MCT). Coherent and incoherent scattering functions are analyzed in terms of the idealized MCT. For temperatures T > T c we provide evidence for the space-time factorization property of the β relaxation and for the time-temperature superposition principle (TTSP) of the α relaxation, and we also discuss deviations from these predictions for T ≈ T c. For T larger than the smallest temperature where the TTSP holds we perform a quantitative analysis of the dynamics with the asymptotic MCT predictions for the late β regime. Within MCT a key quantity, in addition to T c, is the exponent parameter λ. For the fully flexible polymer models studied we find that λ is independent of N and has a value (λ = 0.735 ) typical of simple glass-forming liquids. On the other hand, the critical temperature increases with chain length toward an asymptotic value T c (∞) . This increase can be described by T c (∞) - T c(N) ∼ 1/N and may be interpreted in terms of the N dependence of the monomer density ρ, if we assume that the MCT glass transition is ruled by a soft-sphere-like constant coupling parameter Γ c = ρ c T c (-1/4), where ρ c is the monomer density at T c. In addition, we also estimate T c from a Hansen-Verlet-like criterion and MCT calculations based on structural input from the simulation. For our polymer model both the Hansen-Verlet criterion and the MCT calculations suggest T c to decrease with increasing chain length, in contrast to the direct analysis of the simulation data.

Dynamics of a rigid rod in a disordered medium with long-range spatial correlation

We investigate the diffusion of a rigid rodlike object in a two-dimensional disordered host medium, which consists of static pointlike sources of force. The points are distributed with long-range spatial correlation and interact with the rod via a repulsive potential. The time dependence of the rod's center-of-mass mean-squared displacement and its rotational mean-squared displacement are obtained for various degrees of long-range spatial correlation and rod's lengths. These transport characteristics are compared to those obtained in previous studies for the case of homogeneous distribution of force points. It is shown that existence of long-range correlation among force points makes the center of mass diffusion anomalous.

Applicability of mode-coupling theory to polyisobutylene: a molecular dynamics simulation study

The applicability of Mode Coupling Theory (MCT) to the glass-forming polymer polyisobutylene (PIB) has been explored by using fully atomistic molecular dynamics simulations. MCT predictions for the so-called asymptotic regime have been successfully tested on the dynamic structure factor and the self-correlation function of PIB main-chain carbons calculated from the simulated cell. The factorization theorem and the time-temperature superposition principle are satisfied. A consistent fitting procedure of the simulation data to the MCT asymptotic power-laws predicted for the α-relaxation regime has delivered the dynamic exponents of the theory-in particular, the exponent parameter λ-the critical non-ergodicity parameters, and the critical temperature T(c). The obtained values of λ and T(c) agree, within the uncertainties involved in both studies, with those deduced from depolarized light scattering experiments [A. Kisliuk et al., J. Polym. Sci. Part B: Polym. Phys. 38, 2785 (2000)]. Both, λ and T(c)/T(g) values found for PIB are unusually large with respect to those commonly obtained in low molecular weight systems. Moreover, the high T(c)/T(g) value is compatible with a certain correlation of this parameter with the fragility in Angell's classification. Conversely, the value of λ is close to that reported for real polymers, simulated "realistic" polymers and simple polymer models with intramolecular barriers. In the framework of the MCT, such finding should be the signature of two different mechanisms for the glass-transition in real polymers: intermolecular packing and intramolecular barriers combined with chain connectivity.

Slow processes in supercooled o-terphenyl: relaxation and decoupling

We mapped the relaxation times of inter- and intramolecular correlations in o-terphenyl by a quasielastic scattering method using nuclear resonant scattering of synchrotron radiation. From the obtained map, we found that the slow β process is decoupled from the α process at 278 K, and this temperature is clearly below the previous decoupling temperature of 290 K, at which the α-relaxation dynamics changes. Then, it was also concluded that sufficient solidlike condition achieved by further cooling from 290 K is required to decouple the slow β process from the α process and, due to the difference of the length scales between the α and the slow β processes, these two averaged relaxation times <τ> are concluded not to cross as an extrapolation assumed so far. Furthermore, evidence of the restricted dynamics of the slow β process could be obtained as an anomalous momentum transfer (q) dependence of <τ>(<τ> ∝q(-2.9)) at 265 K, observed at q values of 18-48  nm(-1).

Relaxation properties in a diffusive model of k-mers with constrained movements on a triangular lattice

We study the relaxation process in a two-dimensional lattice gas model, based on the concept of geometrical frustration. In this model the particles are k-mers that can both randomly translate and rotate on the planar triangular lattice. In the absence of rotation, the diffusion of hard-core particles in crossed single-file systems is investigated. We monitor, for different densities, several quantities: mean-square displacement, the self-part of the van Hove correlation function, and the self-intermediate scattering function. We observe a considerable slowing of diffusion on a long-time scale when suppressing the rotational motion of k-mers; our system is subdiffusive at intermediate times between the initial transient and the long-time diffusive regime. We show that the self-part of the van Hove correlation function exhibits, as a function of particle displacement, a stretched exponential decay at intermediate times. The self-intermediate scattering function (SISF), displaying slower than exponential relaxation, suggests the existence of heterogeneous dynamics. For each value of density, the SISF is well described by the Kohlrausch-Williams-Watts law; the characteristic timescale τ(q(n)) is found to decrease with the wave vector q(n) according to a simple power law. Furthermore, the slowing of the dynamics with density ρ(0) is consistent with the scaling law 1/τ(q(n);ρ(0))∝(ρ(c)-ρ(0))(ϰ), with the same exponent ϰ=3.34±0.12 for all wave vectors q(n). The density ρ(c) is approximately equal to the closest packing limit, θ(CPL)≲1, for dimers on the two-dimensional triangular lattice. The self-diffusion coefficient D(s) scales with the same power-law exponent and critical density.

Atomic motions in the alphabeta-merging region of 1,4-polybutadiene: a molecular dynamics simulation study

We present fully atomistic molecular dynamics simulations on 1,4-polybutadiene in a wide temperature range from 200 to 280 K, i.e., in the region where the alpha- and beta-relaxations merge and above. A big computational effort has been performed-especially for the lowest temperatures investigated-to extend the simulation runs to very long times (up to 1 mus for 200 K). The simulated sample has been carefully validated by using previous neutron scattering data on the real sample with similar microstructure. Inspecting the trajectories of the different hydrogens in real space, we have observed a heterogeneous dynamical behavior (each kind of hydrogen moves in a different way) with signatures of combined hopping and diffusive motions in the whole range investigated. The application of a previously proposed model [Colmenero et al., Europhys. Lett. 71, 262 (2005)] is successful and a characterization of the local motions and diffusion is possible. The comparison of our results to those reported in the literature provides a consistent scenario for polybutadiene dynamics and puts into a context the different experimental observations. We also discuss the impact of the hopping processes on the observation and interpretation of experimentally accessible magnitudes and the origin of the deviations from Gaussian behavior in this system.

Computational probes of molecular motion in the Lewis-Wahnström model for ortho-terphenyl

We use molecular dynamics simulations to investigate translational and rotational diffusion in a rigid three-site model of the fragile glass former ortho-terphenyl, at 260 K< or =T< or =346 K and ambient pressure. An Einstein formulation of rotational motion is presented, which supplements the commonly used Debye model. The latter is shown to break down at supercooled temperatures as the mechanism of molecular reorientation changes from small random steps to large infrequent orientational jumps. We find that the model system exhibits non-Gaussian behavior in translational and rotational motion, which strengthens upon supercooling. Examination of particle mobility reveals spatially heterogeneous dynamics in translation and rotation, with a strong spatial correlation between translationally and rotationally mobile particles. Application of the Einstein formalism to the analysis of translation-rotation decoupling results in a trend opposite to that seen in conventional approaches based on the Debye formalism, namely, an enhancement in the effective rate of rotational motion relative to translation upon supercooling.

Relaxation scenarios in a mixture of large and small spheres: Dependence on the size disparity

We present a computational investigation on the slow dynamics of a mixture of large and small soft spheres. By varying the size disparity at a moderate fixed composition different relaxation scenarios are observed for the small particles. For small disparity density-density correlators exhibit moderate stretching. Only small quantitative differences are observed between dynamic features for large and small particles. On the contrary, large disparity induces a clear time scale separation between the large and small particles. Density-density correlators for the small particles become extremely stretched and display logarithmic relaxation by properly tuning the temperature or the wave vector. Self-correlators decay much faster than density-density correlators. For very large size disparity, a complete separation between self- and collective dynamics is observed for the small particles. Self-correlators decay to zero at temperatures where density-density correlations are frozen. The dynamic picture obtained by varying the size disparity resembles features associated with mode coupling transition lines of the types B and A at, respectively, small and very large size disparities. Both lines might merge, at some intermediate disparity, at a higher-order point, to which logarithmic relaxation would be associated. This picture resembles predictions of a recent mode coupling theory for fluids confined in matrices with interconnected voids [V. Krakoviack, Phys. Rev. Lett. 94, 065703 (2005)].

A comparative molecular simulation study of the glass former ortho-terphenyl in bulk and freestanding films

We performed molecular dynamics simulations of the low-molecular weight organic glass former ortho-terphenyl in bulk and freestanding films. The main motivation is to provide molecular insight into the confinement effect without explicit interfaces. Based on earlier models of ortho-terphenyl we developed an atomistic model for bulk simulations. The model reproduces literature data both from simulations and experiments starting from specific volume and diffusivity to mean square displacement and radial distribution functions. After characterizing the bulk model we form freestanding films by the elongation and expansion method. These films give us the opportunity to study the dynamical heterogeneity near the glass transition through in-plane mobility and reorientation dynamics. We finally compare the model in bulk and under confinement. We found qualitatively a lower glass transition temperature for the freestanding film compared to the bulk.

Molecular Dynamics Study of Translational and Rotational Diffusion in Liquid Ortho-terphenyl

NVT molecular dynamics simulations were performed on liquid o-terphenyl as a function of temperature in the range 320-480 K. Computed translational diffusion coefficients displayed the non-Arrhenius behavior expected of a fragile glass-forming liquid and were in good, semiquantitative agreement with experimental results. Rotational correlation functions calculated for various vectors within the molecule exhibited a very short time (0-1 ps) initial decay, followed by a reversal, which corresponds to free reorientation within the "solvent" cage prior to collision with a wall. Rotational correlation times of three orthogonal vectors fixed on the central benzene were close to equal at all temperatures, indicating nearly isotropic overall molecular reorientation. The average correlation times exhibited a non-Arrhenius temperature dependence and were in very good agreement with experimental values derived from 2D and 1H NMR relaxation times. Correlation times of vectors located on the lateral phenyl rings were used to calculate the "spinning" internal rotation diffusion coefficients, which were approximately twice as great as the overall rotational diffusion constants, indicating rapid internal rotation of the phenyl side groups over wide ranges of angle in the liquid.

Hydrogen motions in the α-relaxation regime of poly(vinyl ethylene): A molecular dynamics simulation and neutron scattering study

The hydrogen motion in poly(vinyl ethylene) (1,2-polybutadiene) in the alpha-relaxation regime has been studied by combining neutron spin echo (NSE) measurements on a fully protonated sample and fully atomistic molecular dynamics simulations. The almost perfect agreement between experiment and simulation results validates the simulated cell. A crossover from Gaussian to non-Gaussian behavior is observed for the intermediate scattering function obtained from both NSE measurements and simulations. This crossover takes place at unusually low Q values, well below the first maximum of the static structure factor. Such anomalous deviation from Gaussian behavior can be explained by the intrinsic dynamic heterogeneity arising from the differences in the dynamics of the different protons in this system. Side group hydrogens show a markedly higher mobility than main chain protons. Taking advantage of the simulations we have investigated the dynamic features of all different types of hydrogens in the sample. Considering each kind of proton in an isolated way, deviations from Gaussian behavior are also found. These can be rationalized in the framework of a simple picture based on the existence of a distribution of discrete jumps underlying the atomic motions in the alpha process.

Structural relaxation in supercooled orthoterphenyl

We report molecular-dynamics simulation results performed for a model of molecular liquid orthoterphenyl in supercooled states, which we then compare with both experimental data and mode-coupling-theory (MCT) predictions, aiming at a better understanding of structural relaxation in orthoterphenyl. We pay special attention to the wave number dependence of the collective dynamics. It is shown that the simulation results for the model share many features with experimental data for real system, and that MCT captures the simulation results at the semiquantitative level except for intermediate wave numbers connected to the overall size of the molecule. Theoretical results at the intermediate wave number region are found to be improved by taking into account the spatial correlation of the molecule's geometrical center. This supports the idea that unusual dynamical properties at the intermediate wave numbers, reported previously in simulation studies for the model and discernible in coherent neutron-scattering experimental data, are basically due to the coupling of the rotational motion to the geometrical-center dynamics. However, there still remain qualitative as well as quantitative discrepancies between theoretical prediction and corresponding simulation results at the intermediate wave numbers, which call for further theoretical investigation.

Experimental evidence by neutron scattering of a crossover from Gaussian to non-Gaussian behavior in the alpha relaxation of polyisoprene

We report incoherent quasielastic neutron scattering experiments exploring the alpha-relaxation range in polyisoprene over an unprecedented range in momentum transfer Q. The data corroborate and validate earlier molecular dynamics simulations and reveal the existence of a crossover from Gaussian to non-Gaussian character of the main chain protons self-correlation function in the alpha-relaxation regime. The real challenge of the experiment was to push the neutron techniques to cover a Q range as wide as possible. By combining two neutron spin echo spectrometers and a backscattering instrument, we have been able to study the dynamics in a Q range of 0.1 < or = Q < or = 4.7 A(-1). In the low-Q regime the shape of the relaxation function was found to be related to the dispersion of the relaxation times as predicted by the Gaussian assumption. At short distances or large Q, this relationship is strongly violated indicating a non-Gaussian regime. We have performed a detailed comparison between the experiments and simulations at different temperatures and found, apart from a temperature shift, complete agreement. Combining experiments and simulations led to a consistent interpretation in terms of a distribution of jumps underlying the diffusive motion of protons in the alpha process. This model leads to a time-dependent non-Gaussianity parameter that agrees nearly quantitatively with the simulations and exhibits all features resolved so far from various simulations.

Non-Gaussian nature of the alpha relaxation of glass-forming polyisoprene

We report quasielastic neutron scattering experiments exploring the alpha relaxation in polyisoprene over an unprecedented range in momentum transfer. Corroborating and validating earlier molecular dynamics simulations, the measurements reveal a crossover from a Gaussian regime of sublinear diffusion to a strongly non-Gaussian regime at short distances. We show that a consistent interpretation in terms of a distribution of finite jumps underlying the alpha process is possible. This model leads to a time-dependent non-Gaussian parameter exhibiting all features revealed so far from various simulations.

Dynamics and configurational entropy in the Lewis-Wahnström model for supercooled orthoterphenyl

We study thermodynamic and dynamic properties of a rigid model of the fragile glass-forming liquid orthoterphenyl. This model, introduced by Lewis and Wahnström in 1993, collapses each phenyl ring to a single interaction site; the intermolecular site-site interactions are described by the Lennard-Jones potential whose parameters have been selected to reproduce some bulk properties of the orthoterphenyl molecule. A system of N=343 molecules is considered in a wide range of densities and temperatures, reaching simulation times up to 1 micros. Such long trajectories allow us to equilibrate the system at temperatures below the mode coupling temperature T(c) at which the diffusion constant reaches values of order 10(-10) cm(2)/s and thereby to sample in a significant way the potential energy landscape in the entire temperature range. Working within the inherent structures thermodynamic formalism, we present results for the temperature and density dependence of the number, depth and shape of the basins of the potential energy surface. We evaluate the total entropy of the system by thermodynamic integration from the ideal-noninteracting-gas state and the vibrational entropy approximating the basin free energy with the free energy of 6N-3 harmonic oscillators. We evaluate the configurational part of the entropy as a difference between these two contributions. We study the connection between thermodynamical and dynamical properties of the system. We confirm that the temperature dependence of the configurational entropy and of the diffusion constant, as well as the inverse of the characteristic structural relaxation time, are strongly connected in supercooled states; we demonstrate that this connection is well represented by the Adam-Gibbs relation, stating a linear relation between logD and the quantity 1/TS(c). This relation is found to hold both above and below the critical temperature T(c)-as previously found in the case of silica-supporting the hypothesis that a connection exists between the number of basins and the connectivity properties of the potential energy surface.

Glassy behavior of a homopolymer from molecular dynamics simulations

We study at- and out-of-equilibrium dynamics of a single homopolymer chain at low temperature using molecular dynamics. The main quantities of interest are the average root mean square displacement of the monomers below the theta point, and the structure factor, as a function of time. The observation of these quantities show a close resemblance to those measured in structural glasses and suggest that the polymer chain in its low temperature phase is in a glassy phase, with its dynamics dominated by traps. In equilibrium, at low temperature, we observe the trapping of the monomers and a slowing down of the overall motion of the polymer as well as nonexponential relaxation of the structure factor. Out of equilibrium, at low temperatures, we compute the two-time quantities and observe breaking of ergodicity in a range of waiting times, with the onset of aging.

Molecular dynamics simulation of the fragile glass former orthoterphenyl: a flexible molecule model. II. Collective dynamics

We present a molecular dynamics study of the collective dynamics of a model for the fragile glass former orthoterphenyl. In this model, introduced by Mossa, Di Leonardo, Ruocco, and Sampoli [Phys. Rev. E 62, 612 (2000)], the intramolecular interaction among the three rigid phenyl rings is described by a set of force constants whose value has been fixed in order to obtain a realistic isolated molecule spectrum. The interaction between different molecules is described by a Lennard Jones site-site potential. We study the behavior of the coherent scattering functions F(t)(q,t), considering the density fluctuations of both molecular and phenyl-ring centers of mass; moreover we directly simulate the neutron scattering spectra taking into account both the contributions due to carbon and hydrogens atoms. We compare our results with the main predictions of the mode-coupling theory and with the available coherent neutron scattering experimental data.

Mean-squared displacement of a molecule moving in a glassy system

The mean-squared displacement (MSD) of a hard sphere and of a dumbbell molecule consisting of two fused hard spheres immersed in a dense hard-sphere system is calculated within the mode-coupling theory for ideal liquid-glass transitions. It is proven that the velocity correlator, which is the second time derivative of the MSD, is the negative of a completely monotone function for times within the structural-relaxation regime. The MSD is found to exhibit a large time interval for structural relaxation prior to the onset of the alpha process, which cannot be described by the asymptotic formulas for the mode-coupling-theory-bifurcation dynamics. The alpha process for molecules with a large elongation is shown to exhibit an anomalously wide crossover interval between the end of the von Schweidler decay and the beginning of normal diffusion. The diffusivity of the molecule is predicted to vary nonmonotonically as a function of its elongation.

Dynamics in a supercooled molecular liquid: theory and simulations

We report extensive simulations of liquid supercooled states for a simple three-site molecular model, introduced by Lewis and Wahnström [Phys. Rev. E 50, 3865 (1994)] to mimic the behavior of orthoterphenyl. The large system size and the long simulation length allow us to calculate very precisely (in a large q-vector range) self-correlation and collective correlation functions, providing a clean and simple reference model for theoretical descriptions of molecular liquids in supercooled states. The time and wave-vector dependence of the site-site correlation functions are compared (neglecting the molecular constraints) with detailed ideal mode-coupling theory predictions. Except for the wave-vector region where the dynamics are controlled by the center of mass (around 9 nm(-1)), the theoretical predictions compare very well with the simulation data.

Dynamic heterogeneity of relaxations in glasses and liquids

We report an investigation of the heterogeneity in supercooled liquids and glasses using the non-Gaussianity parameter. We simulate selenium and a binary Lennard-Jones system by molecular dynamics. In the non-Gaussianity three time domains can be distinguished: an increase on the ps scale due to the vibrational (ballistic) motion of the atoms, followed by a growth, due to local relaxations ( beta relaxation) at not too high temperatures, and finally a slow drop at long times. The non-Gaussianity follows in the intermediate time domain a sqrt[t] law. This is explained by collective hopping and dynamic heterogeneity. We support this finding by a model calculation.