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]

Effect of orientational motion of mobile chromophores on the dynamics of Forster energy transfer in polymers

Brownian dynamics (BD) simulations have been carried out to explore the effects of the orientational motion of the donor-acceptor (D-A) chromophore pair on the Forster energy transfer between the D-A pair embedded in a polymer chain in solution. It is found that the usually employed orientational averaging (that is, replacing the orientational factor, kappa, by kappa (2) = 2/3) may lead to an error in the estimation of the rate of the reaction by about 20%. In the limit of slow orientational relaxation, the preaveraging of the orientational factor leads to an overestimation of the rate, while in the opposite limit of very fast orientational relaxation, the usual scheme underestimates the rate. The latter results from an interesting interplay between reaction and diffusion. On the other hand, when one of the chromophores is fixed, the preaveraged rate is found to be fairly reliable if the rotational relaxation of the chromophore is sufficiently fast. The present study also reveals a power law dependence of the FRET rate on the chain length (rate proportional to N- alpha, with alpha approximate to 2.6).

Universality in the fast orientational relaxation near isotropic-nematic transition

Detailed molecular dynamics simulations of Lennard-Jones ellipsoids have been carried out to investigate the emergence of criticality in the single-particle orientational relaxation near the isotropic-nematic (IN) phase transition. The simulations show a sudden appearance of a power-law behavior in the decay of the second-rank orientational relaxation as the IN transition is approached. The simulated value of the power-law exponent is 0.56, which is larger than the mean-field value (0.5) but less than the observed value (0.63) and may be due to the finite size of the simulated system. The decay of the first-rank orientational time correlation function, on the other hand, is nearly exponential but its decay becomes very slow near the isotropic-nematic transition, The zero-frequency rotational friction, calculated from the simulated angular Velocity correlation function, shows a marked increase near the IN transition.

1H and 19F NMR relaxation time studies in (NH4)2ZrF6 superionic conductor

1H and 19F spin-lattice relaxation times in polycrystalline diammonium hexafluorozirconate have been measured in the temperature range of 10–400 K to elucidate the molecular motion of both cation and anion. Interesting features such as translational diffusion at higher temperatures, molecular reorientational motion of both cation and anion groups at intermediate temperatures and quantum rotational tunneling of the ammonium group at lower temperatures have been observed. Nuclear magnetic resonance (NMR) relaxation time results correlate well with the NMR second moment and conductivity studies reported earlier.

H-1 and F-19 NMR relaxation time studies in (NH4)(2)ZrF6 superionic conductor

H-1 and F-19 spin-lattice relaxation times in polycrystalline diammonium hexafluorozirconate have been measured in the temperature range of 10-400 K to elucidate the molecular motion of both cation and anion. Interesting features such as translational diffusion at higher temperatures, molecular reorientational motion of both cation and anion groups at intermediate temperatures and quantum rotational tunneling of the ammonium group at lower temperatures have been observed. Nuclear magnetic resonance (NMR) relaxation time results correlate well with the NMR second moment and conductivity studies reported earlier.

1H NMR study of internal motions and quantum rotational tunneling in (CH3)4NGeCl3

(CH3)4NGeCl3 is prepared, characterized and studied using 1H NMR spin lattice relaxation time and second moment to understand the internal motions and quantum rotational tunneling. Proton second moment is measured at 7 MHz as function of temperature in the range 300-77 K and spin lattice relaxation time (T1) is measured at two Larmor frequencies, as a function of temperature in the range 270-17 K employing a homemade wide-line/pulsed NMR spectrometers. T1 data are analyzed in two temperature regions using relevant theoretical models. The relaxation in the higher temperatures (270-115 K) is attributed to the hindered reorientations of symmetric groups (CH3 and (CH3)4N). Broad asymmetric T1 minima observed below 115 K down to 17 K are attributed to quantum rotational tunneling of the inequivalent methyl groups.

Study of molecular reorientation and quantum rotational tunneling in tetramethylammonium selenate by H-1 NMR

H-1 NMR spin-lattice relaxation time measurements have been carried out in [(CH3)(4)N](2)SeO4 in the temperature range 389-6.6K to understand the possible phase transitions, internal motions and quantum rotational tunneling. A broad T, minimum observed around 280K is attributed to the simultaneous motions of CH3 and (CH3)(4)N groups. Magnetization recovery is found to be stretched exponential below 72 K with varying stretched exponent. Low-temperature T-1 behavior is interpreted in terms of methyl groups undergoing quantum rotational tunneling. (c) 2007 Elsevier Inc. All rights reserved.

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.

Molecular theory of dielectric relaxation in a dense binary dipolar liquid

A molecular theory of dielectric relaxation in a dense binary dipolar liquid is presented. The theory takes into account the effects of intra- and interspecies intermolecular interactions. It is shown that the relaxation is, in general, nonexponential. In certain limits, we recover the biexponential form traditionally used to analyze the experimental data of dielectric relaxation in a binary mixture. However, the relaxation times are widely different from the prediction of the noninteracting rotational diffusion model of Debye for a binary system. Detailed numerical evaluation of the frequency-dependent dielectric function epsilon-(omega) is carried out by using the known analytic solution of the mean spherical approximation (MSA) model for the two-particle direct correlation function for a polar mixture. A microscopic expression for both wave vector (k) and frequency (omega) dependent dielectric function, epsilon-(k,omega), of a binary mixture is also presented. The theoretical predictions on epsilon-(omega) (= epsilon-(k = 0, omega)) have been compared with the available experimental results. In particular, the present theory offers a molecular explanation of the phenomenon of fusing of the two relaxation channels of the neat liquids, observed by Schallamach many years ago.

Excitation and relaxation energies of trans-stilbene: Confined singlet, triplet, and charged bipolarons

The π-electronic excitations and excited-state geometries of trans-stilbene (tS) are found by combining exact solutions of the Pariser-Parr-Pople (PPP) model and semiempirical Parametric Method 3 (PM3) calculations. Comprehensive comparisons with tS spectra are obtained and related to the fluorescence and topological alternation of poly(paraphenylenevinylene) (PPV). The one-photon absorption and triplet of tS correspond, respectively, to singlet and triplet bipolarons confined to two phenyls, while the tS2- ground state is a confined charged bipolaron. Independent estimates of the relaxation energy between vertical and adiabatic excitation show the bipolaron binding energy to depend on both charge and spin, as expected for interacting π electrons in correlated or molecular states. Complete configuration interaction within the PPP model of tS accounts for the singlet-triplet gap, for the fine-structure constants and triplet-triplet spectra, for two-photon transitions and intensities, and for one-photon spectra and the radiative lifetime, although the relative position of nearly degenerate covalent and ionic singlets is not resolved. The planar PM3 geometry and low rotational barrier of tS agree with resolved rotational and vibrational spectra in molecular beams. PM3 excitation and relaxation energies for tS bipolarons are consistent with experiment and with PPP results. Instead of the exciton model, we interpret tS excitations in terms of states that are localized on each ring or extended over an alternating chain, as found exactly in Hückel theory, and find nearly degenerate transitions between extended and localized states in the singlet, triplet, and dianion manifolds. The large topological alternation of the extended system increases the ionicity and interchanges the order of the lowest one- and two-photon absorption of PPV relative to polyenes.

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.

Collective Effects on Single Particle Orientational Relaxation in Slow Dipolar Liquids

Theoretical and computer simulation studies of orientational relaxation in dense molecular liquids are presented. The emphasis of the study is to understand the effects of collective orientational relaxation on the single-particle orientational dynamics. The theoretical analysis is based on a recently developed molecular hydrodynamic theory which allows a self-consistent description of both the collective and the single-particle orientational relaxation. The molecular hydrodynamic theory can be used to derive a relation between the memory function for the collective orientational correlation function and the frequency-dependent dielectric function. A novel feature of the present work is the demonstration that this collective memory function is significantly different from the single-particle rotational friction. However, a microscopic expression for the single-particle rotational friction can be derived from the molecular hydrodynamic theory where the collective memory function can be used to obtain the single-particle orientational friction. This procedure allows, us to calculate the single-particle orientational correlation function near the alpha-beta transition in the supercooled liquid. The calculated correlation function shows an interesting bimodal decay below the bifurcation temperature as the glass transition is approached from above. Brownian dynamics simulations have been carried out to check the validity of the above procedure of translating the memory function from the dielectric relaxation data. We have also investigated the following two issues important in understanding the orientational relaxation in slow liquids. First, we present an analysis of the ''orientational caging'' of translational motion. The value of the translational friction is found to be altered significantly by the orientational caging. Second, we address the question of the rank dependence of the dielectric friction using both simulation and the molecular hydrodynamic theory.

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.

Translational and rotational motion in molecular liquids: A computer simulation study of Lennard-Jones ellipsoids

In order to understand the translational and rotational motion in dense molecular liquids, detailed molecular dynamics simulations of Lennard-Jones ellipsoids have been carried out for three different values of the aspect ratio kappa. For ellipsoids with an aspect ratio equal to 2, the product of the translational diffusion coefficient (D-T) and the average orientational correlation time of the l-th rank harmonics (tau(lR)), converges to a nearly constant value at high density. Surprisingly, this density independent value of D-T tau(lR) is within 5% of the hydrodynamic prediction with the slip boundary condition. This is despite the fact that both D-T and tau(lR) themselves change nearly by an order of magnitude in the density range considered, and the rotational correlation function itself is strongly nonexponential. For small aspect ratios (kappa less than or equal to 1.5), the rotational correlation function remains largely Gaussian even at a very large density, while for a large aspect ratio (kappa greater than or equal to 3), the transition to the nematic liquid-crystalline phase precludes the hydrodynamic regime. Thus, the rotational dynamics of ellipsoids show great sensitivity to the aspect ratio. At low density, tau(lR) goes through a minimum value, indicating the role of interactions in enhancing the rate of orientational relaxation. (C) 1997 American Institute of Physics. [S0021-9606(97)50142-5].

Non-monotonic, distance-dependent relaxation of water in reverse micelles: Propagation of surface induced frustration along hydrogen bond networks

Layer-wise, distance-dependent orientational relaxation of water confined in reverse micelles (RM) is studied using theoretical and computational tools. We use both a newly constructed ``spins on a ring'' (SOR) Ising-type model (with Shore-Zwanzig rotational dynamics) and atomistic simulations with explicit water. Our study explores the effect of reverse micelle size and role of intermolecular correlations, compromised by the presence of a highly polar surface, on the distance (from the interface) dependence of water relaxation. The ``spins on a ring'' model can capture some aspects of distance dependence of relaxation, such as acceleration of orientational relaxation at intermediate layers. In atomistic simulations, layer-wise decomposition of hydrogen bond formation pattern clearly reveals that hydrogen bond arrangement of water at a certain distance away from the surface can remain frustrated due to the interaction with the polar surface head groups. This layer-wise analysis also reveals the presence of a non-monotonic slow relaxation component which can be attributed to this frustration effect and which is accentuated in small to intermediate size RMs. For large size RMs, the long time component decreases monotonically from the interface to the interior of the RMs with slowest relaxation observed at the interface. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4732095]