Anomalous diffusion and cage effects in the isotropic phase of a liquid crystal

The translational motion of 4-n-hexyl-4'-cyanobiphenyl (6CB) in its isotropic phase has been studied by atomistic molecular dynamics simulation from 280 to 330 K. The mean square displacement shows evidence of a subdiffusive dynamics, with a plateau that becomes very apparent at the lowest temperatures. A three-time self-intermediate scattering function reveals that this plateau is connected with a homogeneous dynamics that, at longer times, becomes heterogeneous and finally exponential. These features are shared by, for example, a high-density system of hard spheres, which supports the universal character of the translational dynamics of liquids in their supercooled condition. As predicted by the idealized version of the mode-coupling theory (MCT), the diffusion coefficient dependence upon temperature is well described by a power law, with a critical temperature very close to that obtained by experimental measurements on orientational relaxation. This agreement might indicate a complete freezing of both rotational and translational intradomain dynamics. The time-temperature superposition principle also holds. The shape of the cage that surrounds a 6CB molecule has been reconstructed, and this analysis suggests a preferential side-by-side arrangement of molecules, which locally tend to align their long axes even in the isotropic phase.     

Dynamics in supercooled ionic organic liquids and mode coupling theory analysis

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

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

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

Orientational and interaction induced dynamics in the isotropic phase of a liquid crystal: polarization resolved ultrafast optical Kerr effect spectroscopy

The ultrafast dynamics of the isotropic phase of a liquid crystal 4'-pentyl-4-p-biphenylcarbonitrile (5CB) have been investigated using polarization resolved optical Kerr effect spectroscopy. Measurements were made as a function of both temperature and dilution in nonpolar solvents. To separate single molecule and interaction induced components to the relaxation of the induced birefringence, measurements of both the anisotropic and isotropic response were made. The isotropic response was found to be dominated by a damped low-frequency mode of intramolecular origin. There is a minor additional component assigned to an interaction induced contribution. There is at most an extremely weak isotropic signal beyond 1 ps, showing that the picosecond time scale dynamics of 5CB are dominated by orientational relaxation. The isotropic response is independent of temperature in the range studied (0.2-50 K above the nematic to isotropic phase-transition temperature). The anisotropic response exhibits relaxation dynamics on time scales spanning subpicosecond to several hundred picoseconds and beyond. The fastest components are dominated by a librational response, but there are smaller contributions from three low-frequency intramolecular modes, and a contribution from interaction induced effects. The low-frequency spectral density extracted from these data are independent of temperature in the range studied, 0.2-30 K above the phase-transition temperature, but shift to lower frequency on dilution in alkane solvents. In neat 5CB the picosecond time scale orientational dynamics are dominated by temperature-independent reorientation within the pseudonematic domains, while in solution these are disrupted, and the orientational response becomes faster and temperature dependent.     

Contribution to understanding of the molecular dynamics in liquids

The dielectric relaxation spectroscopy is used for studying the orientational molecular dynamics in the isotropic (I) and nematic (N) phases of two mesogenic liquids composed of the molecules of similar structure and length, but of an essentially different polarity: n-heptylcyanobiphenyl, C(7)H(15)PhPhCN, 7CB (molecular dipole moment mu approximately 5D) and 4-(trans-4'-n-hexylcyclohexyl)isothiocyanatobenzene, C(6)H(13)CyHxPhNCS, 6CHBT (mu approximately 2.5D); advantageously, the temperatures of the I-N phase transition for the two compounds are very close to each other (T(NI) = 316.6 +/- 0.2 K). It is shown that regardless of the differences in polarity of 7CB and 6CHBT molecules and their abilities in dipolar aggregation, the values and temperature dependences of the relaxation time (corresponding to the rotational diffusion of the molecules around their short axis) are very close to each other, in both the isotropic and nematic phases of the liquids studied. Therefore, the data show that the dielectric relaxation processes occurring in dipolar liquids in the isotropic and nematic states lead through the rotational diffusion of individual molecules and the diffusion seems to be not influenced by the intermolecular interactions.     

Conformational and Orientational Order of Molecules in Nematic Liquid Crystals

The statistical properties of the conformational and orientational distribution of molecules with internal rotation in the isotropic and nematic phases of liquid crystal (LC) are investigated in terms of molecular statistical theory. The paper discusses the effect of mutual correlation between the molecular conformational and orientational degrees of freedom on the conformational, orientational, and mixed order parameters of molecules in LC and on the nematic–isotropic liquid transition temperature. For these order parameters, the recurrent relation method is suggested and used to derive the partial functions of the conformational and orientational distributions of molecules in LC.     

Coupling of orientational and translational motions in viscoelastic fluids. Low frequency depolarized (VH) light scattering

Physical processes in supercooled liquids are discussed and it is explained how elasticity might prevail. Consequently a theory is presented for the VH spectrum of low frequency light scattering from an isotropic solid-like visco-elastic medium containing molecules with orientational degrees of freedom. In contrast with the low frequency VH spectrum from non-elastic fluids the present one contains under specified conditions two central lorentzians instead of only one in addition to two pairs of shifted lorentzians and dispersion curves. The information that can be extracted from comparing this theory to experimental data contains a constant which designates the degree of elasticity and also orientational relaxation times, shear viscosities and coupling constants between the translational and the orientational motions.     

Pretransitional behavior of a water in liquid crystal microemulsion close to the demixing transition: evidence for intermicellar attraction mediated by paranematic fluctuations

We present a study of a water-in-oil microemulsion in which surfactant coated water nanodroplets are dispersed in the isotropic phase of the thermotropic liquid-crystal penthyl-cyanobiphenyl (5CB). As the temperature is lowered below the isotropic to nematic phase transition of pure 5CB, the system displays a demixing transition leading to a coexistence of a droplet-rich isotropic phase with a droplet-poor nematic. The transition is anticipated, in the high T side, by increasing pretransitional fluctuations in 5CB molecular orientation and in the nanodroplet concentration. The observed phase behavior supports the notion that the nanosized droplets, while large enough for their statistical behavior to be probed via light scattering, are also small enough to act as impurities, disturbing the local orientational ordering of the liquid crystal and thus experiencing pretransitional attractive interaction mediated by paranematic fluctuations. The pretransitional behavior, together with the topology of the phase diagram, can be understood on the basis of a diluted Lebwohl-Lasher model which describes the nanodroplets simply as holes in the liquid crystal.     

Short time dynamics in the isotropic phase of liquid crystals: the aspect ratio and the power law decay

Optical heterodyne detected optical Kerr effect (OHD-OKE) experiments are used to study the orientational dynamics of the liquid crystal 4^′-octyl-4-biphenylcarbonitrile (8CB) in the isotropic phase near the isotropic to nematic phase transition. The results are compared to those for three other liquid crystals. The 8CB data display a short time scale temperature independent power law decay and a long time scale exponential decay with a temperature dependence described by Landau–de Gennes theory. The power law exponent is −0.56. Combining this result with previous results for three other liquid crystals [J. Chem. Phys. 116 (2002) 6339; J. Chem. Phys. 116 (2002) 360], it is found that the power law exponent depends linearly on the aspect ratio of the liquid crystal.     

Collective orientation dynamics in semi-rigid polymers

The orientational relaxation function in non-dilute isotropic solutions of semi-rigid polymers with different ladder-type oligophenylene monomer structures is recorded by time domain depolarized light scattering. The shape and the peculiar dependence of this function on the scattering angle arise from the coupling between orientational and shear translational motions. At sufficiently high concentrations, this coupling facilitates the opening of a faster relaxation channel better resolved in the forward scattering. The overall shape of the relaxation function is concentration-dependent, and the strength of the coupling appears to be system-specific.     

Multiple short time power laws in the orientational relaxation of nematic liquid crystals

Relaxation in the nematic liquid crystalline phase is known to be sensitive to its proximity to both isotropic and smectic phases. Recent transient optical Kerr effect (OKE) studies have revealed, rather surprisingly, two temporal power laws at short to intermediate times and also an apparent absence of the expected exponential decay at longer times. In order to understand this unusual dynamics, we have carried out extensive molecular dynamics simulations of transient OKE and related orientational time correlation functions in a system of prolate ellipsoids (with aspect ratio equal to 3). The simulations find two distinct power laws, with a crossover region, in the decay of the orientational time correlation function at short to intermediate times (in the range of a few picoseconds to a few nanoseconds). In addition, the simulation results fail to recover any long time exponential decay component. The system size dependence of the exponents suggests that the first power law may originate from the local orientational density fluctuations (like in a glassy liquid). The origin of the second power law is less clear and may be related to the long range fluctuations (such as smecticlike density fluctuations)--these fluctuations are expected to involve small free energy barriers. In support of the latter, the evidence of pronounced coupling between orientational and spatial densities at intermediate wave numbers is presented. This coupling is usually small in normal isotropic liquids, but it is large in the present case. In addition to slow collective orientational relaxation, the single particle orientational relaxation is also found to exhibit slow dynamics in the nematic phase in the long time.     

Liquid crystal properties of the n-alkyl-cyanobiphenyl series from atomistic simulations with ab initio derived force fields

Lengthy molecular dynamics (MD) simulations were performed at constant atmospheric pressure and different temperatures for the series of the 4-n-alkyl-4'-cyanobiphenyls (nCB) with n = 6, 7, and 8. The accurate atomistic force field (Bizzarri, M.; Cacelli, I.; Prampolini, G; Tani, A. J. Phys. Chem. A 2004, 108, 10336), successfully employed to reproduce thermodynamic and transport properties of the 5CB molecule, has here been extended to higher homologues. Nematic and isotropic phases were found for all members of the series, and also, a smectic phase was (tentatively) identified for 8CB at 1 atm and 300 K. Transition temperatures reproduce the experimental values within +/-10 K. Also, structural properties as second and fourth rank orientational order parameters are in good agreement with the corresponding experimental quantities. This means that the well-known odd-even effect, observed for many properties along the nCB series, is well reproduced, despite the narrow range of oscillations, e.g., in clearing temperatures. A detailed analysis of the correlation between molecular properties and odd-even effects is presented.     

Magnetic field induced dipolar couplings in the pretransitional region of a nematic liquid crystal

13C and 1H NMR spectroscopy have been used to study the orientational order which develops when a nematogenic compound, 4-pentyl-4'-cyanobiphenyl (5CB), approaches the transition from the isotropic to the nematic phase at T(NI). The experiments yield values of field induced dipolar couplings, (1), between all of the directly bonded carbon and hydrogen nuclei in the molecule, and 2D(HH)B, the geminal dipolar coupling between protons in the first methylene group in the alkyl chain. The temperature dependence of these couplings shows that in each case they follow a divergence behavior governed by (T - T*)(-1),where T* is a temperature determined from the experimental data and which is close to but less than T(NI). Experiments performed at spectrometer field strengths of 11.75 and 18.79 T confirm the prediction that the induced couplings should depend on the square of the applied field strength. It was found that, within experimental error, T* is the same for each field-induced coupling, and that T(NI) - T* is the same at 11.75 and 18.79 T. It is shown that the set of field-induced couplings 1D(CHi)B obtained at a temperature close to T(NI) can be used to derive a conformer distribution for 5CB in the isotropic phase.     

α-Scale decoupling of the mechanical relaxation and diverging shear wave propagation length scale in triphenylphosphite

We have performed depolarized impulsive stimulated scattering experiments to observe shear acoustic phonons in supercooled triphenylphosphite (TPP) from ~10-500 MHz. These measurements, in tandem with previously performed longitudinal and shear measurements, permit further analyses of the relaxation dynamics of TPP within the framework of the mode coupling theory. Our results provide evidence of α coupling between the shear and longitudinal degrees of freedom up to a decoupling temperature T(c) = 231 K. A lower bound length scale of shear wave propagation in liquids verified the exponent predicted by theory in the vicinity of the decoupling temperature.     

Decoupling charge transport from the structural dynamics in room temperature ionic liquids

Light scattering and dielectric spectroscopy measurements were performed on the room temperature ionic liquid (RTIL) [C4mim][NTf2] in a broad temperature and frequency range. Ionic conductivity was used to estimate self-diffusion of ions, while light scattering was used to study structural relaxation. We demonstrate that the ionic diffusion decouples from the structural relaxation process as the temperature of the sample decreases toward T(g). The strength of the decoupling appears to be significantly lower than that expected for a supercooled liquid of similar fragility. The structural relaxation process in the RTIL follows well the high-temperature mode coupling theory (MCT) scenario. Using the MCT analysis we estimated the dynamic crossover temperature in [C4mim][NTf2] to be T(c) ~ 225 ± 5 K. However, our analysis reveals no sign of the dynamic crossover in the ionic diffusion process.     

Relaxation processes of the liquid crystal ME6N in the isotropic phase studied by Raman scattering experiments

We have investigated the Raman profiles of the nu(C[Triple Bond]N) and nu(C=O) vibrational modes of the nematic liquid crystal ME6N (4-cyanophenyl-4(')-hexylbenzoate) in the isotropic phase at different temperatures and used them as probes of the dynamics and structural organization of this liquid. The vibrational time correlation functions of the nu(C[Triple Bond]N) mode, rather adequately interpreted within the assumption of exponential modulation function (the Kubo-Rothschild theory), indicate that the system experiences an intermediate dynamical regime that gets only slightly faster with increasing temperature. However, this theory fails in predicting the non-exponential behavior that the time correlation functions manifest in the long time range (t>3 ps). For this reason we have additionally approached the interpretation of vibrational correlation functions in terms of the theory formulated by Rothschild and co-workers for locally structured liquids. The application of this theory reveals that the molecular dynamics in this liquid crystal in the isotropic phase is that deriving from a distribution of differently sized clusters, which narrows as the temperature increases. Even at the highest temperature reached in this study (87 degrees C above the nematic-isotropic transition), the liquid has not yet achieved the structure of the simple liquid and the dynamics has not reached the limit of the single channel process. The vibrational and orientational relaxations occur in very different time scales. The temperature independence of the orientational dynamics in the whole range from 55 degrees C to 135 degrees C has been referred to the nonhydrodynamic behavior of the system, arising when local pseudonematic structures persist for times longer than the orientational relaxation. The occurrence of the process of resonant vibrational energy transfer between the C=O groups of adjacent molecules has been revealed in the isotropic phase by a slightly positive Raman noncoincidence effect in the band associated with the nu(C=O) mode. A qualitative interpretation is tentatively given in terms of partial cancellation of contributions deriving from structures having opposite orientations of their C=O groups.     

Onset of slow dynamics in difluorotetrachloroethane glassy crystal

Complementary neutron spin-echo and x-ray experiments and molecular-dynamics simulations have been performed on difluorotetrachloroethane (CFCl2-CFCl2) glassy crystal. Static, single-molecule reorientational dynamics and collective dynamics properties are investigated. Our results confirm the strong analogy between molecular liquids and plastic crystals. The orientational disorder is characterized at different temperatures and a change in the nature of rotational dynamics is observed. A careful check of the rotational diffusion model is performed using self-angular correlation functions Cl with high l values and compared to results obtained on molecular liquids composed of A-B dumbbells. Below the crossover temperature at which slow dynamics emerge, we show that some scaling predictions of the mode coupling theory hold and that alpha-relaxation times and nonergodicity parameters are controlled by the nontrivial static correlations.     

Structural relaxation and rheological response of a driven amorphous system

The interplay between the structural relaxation and the rheological response of a simple amorphous system {a 80:20 binary Lennard-Jones mixture [W. Kob and H. C. Andersen, Phys. Rev. Lett. 73, 1376 (1994)]} is studied via molecular dynamics simulations. In the quiescent state, the model is well known for its sluggish dynamics and a two step relaxation of correlation functions at low temperatures. An ideal glass transition temperature of Tc=0.435 has been identified in the previous studies via the analysis of the system's dynamics in the framework of the mode coupling theory of the glass transition [W. Kob and H. C. Andersen, Phys. Rev. E 51, 4626 (1995)]. In the present work, we focus on the question whether a signature of this ideal glass transition can also be found in the case where the system's dynamics is driven by a shear motion. Indeed, the following distinction in the structural relaxation is found: In the supercooled state, the structural relaxation is dominated by the shear at relatively high shear rates gamma, whereas at sufficiently low gamma the (shear-independent) equilibrium relaxation is recovered. In contrast to this, the structural relaxation of a glass is always driven by shear. This distinct behavior of the correlation functions is also reflected in the rheological response. In the supercooled state, the shear viscosity eta decreases with increasing shear rate (shear thinning) at high shear rates, but then converges toward a constant as the gamma is decreased below a (temperature-dependent) threshold value. Below Tc, on the other hand, the shear viscosity grows as eta proportional, etax 1/gamma, suggesting a divergence at gamma=0. Thus, within the accessible observation time window, a transition toward a nonergodic state seems to occur in the driven glass as the driving force approaches zero. As to the flow curves (stress versus shear rate), a plateau forms at low shear rates in the glassy phase. A consequence of this stress plateau for Poiseuille-type flows is demonstrated.     

Computer simulation study of bulk atactic polystyrene in the vicinity of the glass transition

Molecular dynamics (MD) simulations of bulk atactic polystyrene have been performed in a temperature range from 100 K to 650 K at atmospheric pressure. Local translational mobility has been investigated by measuring the mean square translational displacements of monomers. The long-time asymptotic slope of these dependencies is 0.54 at T>T_g, showing Rouse behavior. Cross-over from motion in the cage to Rouse like dynamics has been studied at T>T_g with a characteristic crossover time follows a power law behavior as a function of T, as predicted by mode-coupling theory (MCT). Local orientational mobility has been studied via the orientational autocorrelation functions, ACFs, (Legendre polynomials of the first and second, order) of both the main-chain and side-group bonds. The relaxation times of the orientational α-relaxation follow the same power law (γ∼2.9) as the characteristic translational diffusion time. Below T>T_g both types of dynamics are described by the same activated law. The ACFs time-distribution functions reveal the existence of activated local rearrangements already above T>T_g.     

Molecular dynamics investigation of a density-driven glass transition in a liquid crystal system

Molecular dynamics simulations are carried out to address the density-driven glass transition in a system of rodlike particles that interact with the Gay-Berne potential. Since crystallization occurs in this system on the time scale of the simulations, direct simulation of the glass transition is not possible. Instead, glasses with isotropic orientational order are heated to a temperature T, and the relaxation times by which nematic orientational order develops are determined. These relaxation times appear to diverge at a critical density rho(c); i.e., the system can equilibrate at rhorho(c) (at the temperature T). The relaxation times follow a power-law scaling as the critical density is approached, suggesting that this density-driven glass transition concurs with mode coupling theory.