Subquadratic wavenumber dependence of the structural relaxation of supercooled liquid in the crossover regime

As a liquid is progressively supercooled toward its glass transition temperature, an intriguing weakening of the wavenumber (q) dependence of the structural relaxation time tau(q) in the intermediate-to-large q limit is observed both in experiments and simulation studies. Neither continuous Brownian diffusive dynamics nor discontinuous activated events can alone explain the anomalous wavenumber dependence. Here we use our recently developed theory that unifies the mode coupling theory for continuous dynamics, with the random first order transition theory treatment of activated discontinuous motion as a nucleationlike instanton process to understand the wavenumber dependence of density relaxation. The predicted smooth change in mechanism of relaxation from diffusive to activated, in the crossover regime, is wavevector dependent and appears to be responsible for the observed subquadratic,nalmost linear, q dependence of the relaxation time.

NMR relaxation study of disorder in condensed matter: solid solutions of betaine phosphate and phosphite

Proton NMR relaxation measurements have been carried out in anti-ferroelectric Betaine phosphate (BP), ferroelectric Betaine phosphite (BPI) and the mixed system BPI(1-x)BPx, at 11.4MHz and 23.3MHz from 300K to 80K for x=0.0, 0.25, 0.45, 0.85, and 1.0. The temperature dependence of spin lattice relaxation time T, exhibits two minima as expected from the BPP model in BP and BPI. The Larmor frequency dependence of T, in the mixed system is rather unusual and exhibits different slopes for the low temperature wings at the two frequencies, which is a clear experimental evidence of the presence of different methyl groups with different activation energies (E-a) indicating disorder.

A molecular theory of collective orientational relaxation in pure and binary dipolar liquids

A molecular theory of collective orientational relaxation of dipolar molecules in a dense liquid is presented. Our work is based on a generalized, nonlinear, Smoluchowski equation (GSE) that includes the effects of intermolecular interactions through a mean‐field force term. The effects of translational motion of the liquid molecules on the orientational relaxation is also included self‐consistently in the GSE. Analytic expressions for the wave‐vector‐dependent orientational correlation functions are obtained for one component, pure liquid and also for binary mixtures. We find that for a dipolar liquid of spherical molecules, the correlation function ϕ(k,t) for l=1, where l is the rank of the spherical harmonics, is biexponential. At zero wave‐vector, one time constant becomes identical with the dielectric relaxation time of the polar liquid. The second time constant is the longitudinal relaxation time, but the contribution of this second component is small. We find that polar forces do not affect the higher order correlation functions (l>1) of spherical dipolar molecules in a linearized theory. The expression of ϕ(k,t) for a binary liquid is a sum of four exponential terms. We also find that the wave‐vector‐dependent relaxation times depend strongly on the microscopic structure of the dense liquid. At intermediate wave vectors, the translational diffusion greatly accelerates the rate of orientational relaxation. The present study indicates that one must pay proper attention to the microscopic structure of the liquid while treating the translational effects. An analysis of the nonlinear terms of the GSE is also presented. An interesting coupling between the number density fluctuation and the orientational fluctuation is uncovered.

Polarization relaxation, dielectric dispersion, and solvation dynamics in dense dipolar liquid

A unified treatment of polarization relaxation, dielectric dispersion and solvation dynamics in a dense, dipolar liquid is presented. It is shown that the information of solvent polarization relaxation that is obtained by macroscopic dielectric dispersion experiments is not sufficient to understand dynamics of solvation of a newly created ion or dipole. In solvation, a significant contribution comes from intermediate wave vector processes which depend critically on the short range (nearest‐neighbor) spatial and orientational order that are present in a dense, dipolar liquid. An analytic expression is obtained for the time dependent solvation energy that depends, in addition to the translational and rotational diffusion coefficients of the liquid, on the ratio of solute–solvent molecular sizes and on the microscopic structure of the polar liquid. Mean spherical approximation (MSA) theory is used to obtain numerical results for polarization relaxation, for wave vector and frequency dependent dielectric function and for time dependent solvation energy. We find that in the absence of translational contribution, the solvation of an ion is, in general, nonexponential. In this case, the short time decay is dominated by the longitudinal relaxation time but the long time decay is dominated by much slower large wave vector processes involving nearest‐neighbor molecules. The presence of a significant translational contribution drastically alters the decay behavior. Now, the long‐time behavior is given by the longitudinal relaxation time constant and the short time dynamics is controlled by the large wave vector processes. Thus, although the continuum model itself is conceptually wrong, a continuum model like result is recovered in the presence of a sizeable translational contribution. The continuum model result is also recovered in the limit of large solute to solvent size ratio. In the opposite limit of small solute size, the decay is markedly nonexponential (if the translational contribution is not very large) and a complete breakdown of the continuum model takes place. The significance of these results is discussed.

Dynamics of polar solvation: Route to single exponential relaxation via translational diffusion

A microscopic theoretical calculation of time-dependent solvation energy shows that the solvation of an ion or a dipole is dominated by a single relaxation time if the translational contribution to relaxation is significant.

A kinetic Ising model study of dynamical correlations in confined fluids: Emergence of both fast and slow time scales

Experiments and computer simulation studies have revealed existence of rich dynamics in the orientational relaxation of molecules in confined systems such as water in reverse micelles, cyclodextrin cavities, and nanotubes. Here we introduce a novel finite length one dimensional Ising model to investigate the propagation and the annihilation of dynamical correlations in finite systems and to understand the intriguing shortening of the orientational relaxation time that has been reported for small sized reverse micelles. In our finite sized model, the two spins at the two end cells are oriented in the opposite directions to mimic the effects of surface that in real system fixes water orientation in the opposite directions. This produces opposite polarizations to propagate inside from the surface and to produce bulklike condition at the center. This model can be solved analytically for short chains. For long chains, we solve the model numerically with Glauber spin flip dynamics (and also with Metropolis single-spin flip Monte Carlo algorithm). We show that model nicely reproduces many of the features observed in experiments. Due to the destructive interference among correlations that propagate from the surface to the core, one of the rotational relaxation time components decays faster than the bulk. In general, the relaxation of spins is nonexponential due to the interplay between various interactions. In the limit of strong coupling between the spins or in the limit of low temperature, the nature of relaxation of the spins undergoes a qualitative change with the emergence of a homogeneous dynamics where decay is predominantly exponential, again in agreement with experiments. (C) 2010 American Institute of Physics. doi: 10.1063/1.3474948]

Relationship between microscopic and macroscopic orientational relaxation times in polar liquids

Microscopic relations between single-particle orientational relaxation time (T, ) , dielectric relaxation time ( T ~ )a,n d many-body orientational relaxation time ( T ~o)f a dipolar liquid are derived. We show that both T~ and T~ are influenced significantly by many-body effects. In the present theory, these many-body effects enter through the anisotropic part of the two-particle direct correlation function of the polar liquid. We use mean-spherical approximation (MSA) for dipolar hard spheres for explicit numerical evaluation of the relaxation times. We find that, although the dipolar correlation function is biexponential, the frequency-dependent dielectric constant is of simple Debye form, with T~ equal to the transverse polarization relaxation time. The microscopic T~ falls in between Debye and Onsager-Glarum expressions at large values of the static dielectric constant.

Carbon-13 and proton relaxation study of backbone dynamics of poly(vinyl acetate) in solution

The dynamics of poly(vinyl acetate) in toluene solution has been examined by C-13 and proton relaxation. C-13 spin-lattice relaxation time and nuclear Overhauser enhancement measurements were carried out as a function of temperature at 50.3 and 100.6 MHz. The spin-lattice relaxation times for backbone protons were measured at different temperatures at 200 MHz. The relaxation data have been analyzed using the Hall-Weber-Helfand (HWH) model, which describes backbone dynamics in terms of conformational transitions and the Dejean-Laupretre-Monnerie (DLM) model, which includes bond librations in addition to conformational transitions. The parameters obtained from the analysis of C-13 relaxation data were utilized to predict the proton relaxation data. The DLM model was found to be more successful in reproducing the experimental results. To study the influence of libration further, proton relaxation data for poly(vinyl acetate) over a wider range of temperature reported in the literature were analyzed by these two models. The DLM model could reproduce the experimental data at all temperatures whereas the HWH model was found to be successful only in accounting for the experimental data at high temperatures. The results demonstrate the importance of including the librational mode in the description of the backbone dynamics in polymers.

Carbon-13 relaxation study of chain segmental motion and side-chain internal rotations of poly(isobutyl methacrylate) in solution

The dynamics of poly(isobutyl methacrylate) in toluene solution has been examined by C-13 spin-lattice relaxation time and NOE measurements as a function of temperature. The experiments were performed at 50.3 and 100.6 MHz. The backbone carbon relaxation data have been analyzed using the Dejean-Laupretre-Monnerie (DLM) model, which describes the dynamical processes in the backbone in terms of conformational transitions and bond librations. The relaxation data of the side chain nuclei have been analyzed by assuming different motional models, namely, unrestricted rotational diffusion, three site jumps, and restricted rotational diffusion. The different models have been compared for their ability to reproduce the experimental spin-lattice relaxation times and also to predict the behavior of NOE as a function of temperature. Conformational energy calculations have been carried out on a model compound by using the semiempirical quantum chemical method, AM1, and the results confirm the validity of the motional models used to describe the side-chain motion.

Dynamic Behavior of Poly(2-vinylpyridine) in Solution Studied by 13C Nuclear Magnetic Relaxation

Abstract: The dynamics of poly(2-vinylpyridine) in chloroform solution has been examined by C-13 spin-lattice relaxation time and NOE measurements as a function of temperature. The experiments were performed at 50.3 and 100.6 MHz. The backbone carbon relaxation data have been analyzed in terms of six motional models. Among these models, the models which consider conformational transitions and bond librations for the backbone were found to be more successful. Pyridyl ring motion has been modeled as a restricted rotation with the rotational amplitude varying with temperature. The activation energy parameters obtained from the relaxation data of the pyridyl ring carbon have been compared with the energy barrier for ring rotation estimated from conformational energy calculations using the AM1 semiempirical quantum chemical method. The results of the conformational energy calculations support the description of pyridyl ring motion as a restricted rotation.

Relaxation behavior of twin nonlinear optical chromophores: Effect of the spacer length

A new series of twin nonlinear optical (NLO) molecules, having two 4-nitrophenol chromophores that are linked via a flexible polymethylene spacer of varying length [(CH2)(n), n = 1-12], were synthesized. Powder second harmonic generation measurements of these twin samples indicated a pronounced odd-even oscillation, with the odd twins exhibiting a high SHG value while the even ones gave no measurable SH signal. This behavior reflects the crystal packing preferences in such twin NLO systems that have odd and even numbers of atoms linking them - the even ones appear to prefer a centrosymmetric packing arrangement. The orientational/disordering dynamics of these twin NLO molecules, doped in a polymer (poly(methyl methacrylate)) matrix, has also been studied using SHG in electric field poled samples. Interestingly, the maximum attainable SH signal, chi((2)), in, the poled samples also showed an odd-even oscillation; the odd ones again having a higher value of chi((2)) This unprecedented odd-even oscillation in such molecularly doped systems is rationalized as being due to the intrinsically greater ease of a parallel alignment of the two chromophores in the twins with an odd spacer than in those with an even one. Further, the temporal stability of the SHG intensity at 70 degrees C, after the removal of the applied corona, was also studied. The relaxation of all the twin chromophores followed a biexponential decay; the characteristic relaxation time (tau(2)) for the slow decay component suggests that while the twin with a single methylene unit relaxes relatively slowly, the relaxation is significantly faster in cases where n = 2 and 3. In the twins with even longer spacer segments, the relaxation again becomes slower and reaches a saturation value. The observed minimum appears to reflect the interplay of two competing factors that affect the chromophore alignment in such twin systems, namely, the electrostatic repulsion between neighboring oriented dipoles and the intrinsic flexibility of the spacer.

Molecular Dynamics of Thermotropic Liquid Crystals: Anomalous relaxation dynamics of calamitic and discotic liquid crystals

Recent optical kerr effect (OKE) studies have demonstrated that orientational relaxation of rod-like nematogens exhibits temporal power law decay at intermediate times not only near the isotropic–nematic (I–N) phase boundary but also in the nematic phase. Such behaviour has drawn an intriguing analogy with supercooled liquids. We have investigated both collective and single-particle orientational dynamics of a family of model system of thermotropic liquid crystals using extensive computer simulations. Several remarkable features of glassy dynamics are on display including non-exponential relaxation, dynamical heterogeneity, and non-Arrhenius temperature dependence of the orientational relaxation time. Over a temperature range near the I–N phase boundary, the system behaves remarkably like a fragile glass-forming liquid. Using proper scaling, we construct the usual relaxation time versus inverse temperature plot and explicitly demonstrate that one can successfully define a density dependent fragility of liquid crystals. The fragility of liquid crystals shows a temperature and density dependence which is remarkably similar to the fragility of glass forming supercooled liquids. Energy landscape analysis of inherent structures shows that the breakdown of the Arrhenius temperature dependence of relaxation rate occurs at a temperature that marks the onset of the growth of the depth of the potential energy minima explored by the system. A model liquid crystal, consisting of disk-like molecules, has also been investigated in molecular dynamics simulations for orientational relaxation along two isobars starting from the high temperature isotropic phase. The isobars have been so chosen that the phase sequence isotropic (I)–nematic (N)–columnar (C) appears upon cooling along one of them and the sequence isotropic (I)–columnar(C) along the other. While the orientational relaxation in the isotropic phase near the I–N phase transition shows a power law decay at short to intermediate times, such power law relaxation is not observed in the isotropic phase near the I–C phase boundary. The origin of the power law decay in the single-particle second-rank orientational time correlation function (OTCF) is traced to the growth of the orientational pair distribution functions near the I–N phase boundary. As the system settles into the nematic phase, the decay of the single-particle second-rank orientational OTCF follows a pattern that is similar to what is observed with calamitic liquid crystals and supercooled molecular liquids.

Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K ion are comparable to that observed in bulk salt solutions. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4717710]

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R-G,R-bulk can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time tau(r) of a Gaussian chain of polymerization index N and persistence length l(0). The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D similar to dN(-1) ln(N/d(2)) and tau(r) similar to N(2)d(-1) ln(N/d(2))(-1) in the regime of moderate confinement, where l(0) << d < R-G,R-bulk. The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed.

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R-G,R-bulk can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time tau(r) of a Gaussian chain of polymerization index N and persistence length l(0). The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D similar to dN(-1) ln(N/d(2)) and tau(r) similar to N(2)d(-1) ln(N/d(2))(-1) in the regime of moderate confinement, where l(0) << d < R-G,R-bulk. The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed. (C) 2013 AIP Publishing LLC.