Solvation dynamics in nonassociated polar solvents

Ultrafast solvation dynamics in three nonassociated polar solvents, namely, acetonitrile, dimethyl sulfoxide, and acetone, have been studied by using the molecular hydrodynamic theory. For solvation in acetonitrile, the solvent memory function required for this study has been obtained from recent dielectric relaxation measurements of Venabales and Schuttenmaer; earlier theoretical studies used only the Kerr relaxation data. As the latter provides only an indirect information regarding the polar dynamical response of the dipolar liquid, it fails to provide a fully quantitative description of the solvation time correlation function, S(t). The present study with full dielectric data, on the other hand, gives excellent agreement with the experimental results. The theory shows that the ultrafast part of the solvation dynamics originates almost entirely from the high-frequency component of dielectric relaxation (with time constant 0.177 ps), although the latter represents only a small part of the latter. For DMSO and acetone, however, the present theory predicts a decay slower than the experimental observation. It is proposed that for these two solvents specific chromophore-solvent interactions might be responsible for the-large discrepancy. On the basis of the theory, two experimental studies have also been proposed.

Solvation dynamics of a quadrupolar solute in dipolar liquids

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

Slow orientational dynamics of water molecules at a micellar surface

The dynamics of water molecules near an aqueous micellar interface is studied in an atomistic molecular dynamics simulation of cesium pentadecafluorooctanoate (CsPFO) in water. The dipolar orientational time correlation function (tcf) and the translational diffusion of the water molecules are investigated. Results show that both the reorientational and the translational motion of water molecules near the micelle are restricted. In particular, the orientational tcf exhibits a very slow component in the long time which is slower than its bulk value by 2 orders of magnitude. This slow decay seems to be related to the slow decay often observed in experiments. The origin of the slow decay is analyzed.

String-like propagation of the 5-coordinated defect state in supercooled water: molecular origin of dynamic and thermodynamic anomalies

We find that at low temperature water, large amplitude (similar to 60 degrees) rotational jumps propagate like a string, with the length of propagation increasing with lowering temperature. The strings are formed by mobile 5-coordinated water molecules which move like a Glarum defect (J. Chem. Phys., 1960, 33, 1371), causing water molecules on the path to change from 4-coordinated to 5-coordinated and again back to 4-coordinated water, and in the process cause the tagged water molecule to jump, by following essentially the Laage-Hynes mechanism (Science, 2006, 311, 832-835). The effects on relaxation of the propagating defect causing large amplitude jumps are manifested most dramatically in the mean square displacement (MSD) and also in the rotational time correlation function of the O-H bond of the molecule that is visited by the defect (transient transition to the 5-coordinated state). The MSD and the decay of rotational time correlation function, both remain quenched in the absence of any visit by the defect, as postulated by Glarum long time ago. We establish a direct connection between these propagating events and the known thermodynamic and dynamic anomalies in supercooled water. These strings are found largely in the regions that surround the relatively rigid domains of 4-coordinated water molecules. The propagating strings give rise to a noticeable dynamical heterogeneity, quantified here by a sharp rise in the peak of the four-point density response function, chi(4)(t). This dynamics heterogeneity is also responsible for the breakdown of the Stokes-Einstein relation.

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.

Distribution of transverse chain fluctuations in harmonically confined semiflexible polymers

Two different experimental studies of polymer dynamics based on single-molecule fluorescence imaging have recently found evidence of heterogeneities in the widths of the putative tubes that surround filaments of F-actin during their motion in concentrated solution. In one J. Glaser, D. Chakraborty, K. Kroy, I. Lauter, M. Degawa, N. Kirchesner, B. Hoffmann, R. Merkel, and M. Giesen, Phys. Rev. Lett. 105, 037801 (2010)], the observations were explained in terms of the statistics of a worm-like chain confined to a potential determined self-consistently by a binary collision approximation, and in the other B. Wang, J. Guan, S. M. Anthony, S. C. Bae, K. S. Schweizer, and S. Granick, Phys. Rev. Lett. 104, 118301 (2010)], they were explained in terms of the scaling properties of a random fluid of thin rods. In this paper, we show, using an exact path integral calculation, that the distribution of the length-averaged transverse fluctuations of a harmonically confined weakly bendable rod (one possible realization of a semiflexible chain in a tube), is in good qualitative agreement with the experimental data, although it is qualitatively different in analytic structure from the earlier theoretical predictions. We also show that similar path integral techniques can be used to obtain an exact expression for the time correlation function of fluctuations in the tube cross section. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4712306]

Solvation dynamics of tryptophan in water-dimethyl sulfoxide binary mixture: In search of molecular origin of composition dependent multiple anomalies

Experimental and simulation studies have uncovered at least two anomalous concentration regimes in water-dimethyl sulfoxide (DMSO) binary mixture whose precise origin has remained a subject of debate. In order to facilitate time domain experimental investigation of the dynamics of such binary mixtures, we explore strength or extent of influence of these anomalies in dipolar solvation dynamics by carrying out long molecular dynamics simulations over a wide range of DMSO concentration. The solvation time correlation function so calculated indeed displays strong composition dependent anomalies, reflected in pronounced non-exponential kinetics and non-monotonous composition dependence of the average solvation time constant. In particular, we find remarkable slow-down in the solvation dynamics around 10%-20% and 35%-50% mole percentage. We investigate microscopic origin of these two anomalies. The population distribution analyses of different structural morphology elucidate that these two slowing down are reflections of intriguing structural transformations in water-DMSO mixture. The structural transformations themselves can be explained in terms of a change in the relative coordination number of DMSO and water molecules, from 1DMSO:2H(2)O to 1H(2)O:1DMSO and 1H(2)O:2DMSO complex formation. Thus, while the emergence of first slow down (at 15% DMSO mole percentage) is due to the percolation among DMSO molecules supported by the water molecules (whose percolating network remains largely unaffected), the 2nd anomaly (centered on 40%-50%) is due to the formation of the network structure where the unit of 1DMSO:1H(2)O and 2DMSO:1H(2)O dominates to give rise to rich dynamical features. Through an analysis of partial solvation dynamics an interesting negative cross-correlation between water and DMSO is observed that makes an important contribution to relaxation at intermediate to longer times.

Anomalous power law decay in solvation dynamics of DNA: a mode coupling theory analysis of ion contribution

Several time dependent fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long-range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect, yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Finally, we propose several experiments to check the prediction of the present theoretical work.

Dynamics of the reaction between the free end of a tethered self-avoiding polymer and a flat penetrable surface: A renormalization group study

The average time tau(r) for one end of a long, self-avoiding polymer to interact for the first time with a flat penetrable surface to which it is attached at the other end is shown here to scale essentially as the square of the chain's contour length N. This result is obtained within the framework of the Wilemski-Fixman approximation to diffusion-limited reactions, in which the reaction time is expressed as a time correlation function of a ``sink'' term. In the present work, this sink-sink correlation function is calculated using perturbation expansions in the excluded volume and the polymer-surface interactions, with renormalization group methods being used to resum the expansion into a power law form. The quadratic dependence of tau(r) on N mirrors the behavior of the average time tau(c) of a free random walk to cyclize, but contrasts with the cyclization time of a free self-avoiding walk (SAW), for which tau(r) similar to N-2.2. A simulation study by Cheng and Makarov J. Phys. Chem. B 114, 3321 (2010)] of the chain-end reaction time of an SAW on a flat impenetrable surface leads to the same N-2.2 behavior, which is surprising given the reduced conformational space a tethered polymer has to explore in order to react. (C) 2014 AIP Publishing LLC.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

Laser light scattering of the molecular weight distribution of unfractionated phenolphthalein poly(aryl ether sulfone)

Two unfractionated samples of phenolphthalein poly( aryl ether sulfone) (PES-C) were characterized in CHCl3 at 25 degrees C by applying a recently developed laser light-scattering (LLS) procedure. The Laplace inversion of precisely measured intensity-intensity time correlation function lead us first to an estimate of the characteristic line-width distribution G(Gamma) and then to the translational diffusion coefficient distribution G(D). A combination of static and dynamic LLS results enabled us to determine D = (2.69 x 10(-4))M(-0.553), which agrees with the calibration of D = (2.45 x 10(-4))M(-0.55) previously established by a set of narrowly distributed PES-C samples. Using this newly obtained scaling between D and M, we were able to convert G(D) into a differential weight distribution f(w)(M) for the two PES-C samples. The weight-average molecular weights calculated from f(w)(M) are comparable to that obtained directly from static LLS. Our results showed that using two broadly distributed samples instead of a set of narrowly distributed samples have provided not only similar final results, but also a more practical method for the PES-C characterization. (C) 1997 John Wiley & Sons, Inc.

Dynamic light-scattering characterization of the molecular weight distribution of a broadly distributed phenolphthalein poly(aryl ether ketone)

Using a recently developed laser light-scattering (LLS) procedure, we accomplished the characterization of a broadly distributed unfractionated phenolphthalein poly(aryl ether ketone) (PEK-C) in CHCl3 at 25 degrees C. The laplace inversion of precisely measured intensity-intensity time correlation function from dynamic LLS leads us first to an estimate of the characteristic line-width distribution G(Gamma) and then to the translational diffusion coefficient distribution G(D). By using a previously established calibration of D (cm(2)/s) = 2.37 X 10(-4)M(-0.57), were able to convert G(D) into a differential weight distribution f(w)(M). The weight-average molecular weight M(w) calculated from f(w)(M) agrees well with that directly measured in static LLS. Our results indicate that both the calibration and LLS procedure used in this study are ready to be applied as a routine method for the characterization of the molecular weight distribution of PEK-C. (C) 1996 John Wiley & Sons, Inc.

A DYNAMIC LASER LIGHT-SCATTERING STUDY OF CHITOSAN IN AQUEOUS-SOLUTION

The solution behavior of four chitosans (91% deacetylated chitin) with different molecular weights in 0.2M CH3COOH/0.1M CH3COONa aqueous solution was investigated at 25 degrees C by dynamic laser light scattering (LLS). The Laplace inversion of the precisely measured intensity-intensity time correlation function leads us to an estimate of the line-width distribution G(Gamma), which could be further reduced to a translational diffusion coefficient distribution G(D). By using a combination of static and dynamic LLS results, i.e. Mw and G(D), we were able to establish a calibration of D = k(D)M(-alpha D) with k(D) = (3.14 +/- 0.20) X 10(-4) and alpha(D) = 0.655 +/- 0.015. By using this calibration, we successfully converted G(D) into a molecular weight distribution f(w)(M). The larger alpha(D) value confirms that the chitosan chain is slightly extended in aqueous solution even in the presence of salts. This is mainly due to its backbone and polyelectrolytes nature. As a very sensitive technique, our dynamic LLS results also revealed that even in dilute solution chitosan still forms a small amount of larger sized aggregates that have ben overlooked in previous studies. The calibration obtained in this study will provide another way to characterize the molecular weight distribution of chitosan in aqueous solution at room temperature. (C) 1995 John Wiley & Sons, Inc.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.

The stochastic nature of predator-prey cycles

We study by numerical simulations the time correlation function of a stochastic lattice model describing the dynamics of coexistence of two interacting biological species that present time cycles in the number of species individuals. Its asymptotic behavior is shown to decrease in time as a sinusoidal exponential function from which we extract the dominant eigenvalue of the evolution operator related to the stochastic dynamics showing that it is complex with the imaginary part being the frequency of the population cycles. The transition from the oscillatory to the nonoscillatory behavior occurs when the asymptotic behavior of the time correlation function becomes a pure exponential, that is, when the real part of the complex eigenvalue equals a real eigenvalue. We also show that the amplitude of the undamped oscillations increases with the square root of the area of the habitat as ordinary random fluctuations. (C) 2009 Elsevier B.V. All rights reserved.