In search of temporal power laws in the orientational relaxation near isotropic–nematic phase transition in model nematogens

Measuring order in disordered systems and disorder in ordered systems: Random matrix theory for isotropic and nematic liquid crystals and its perspective on pseudo-nematic domains

Surprisingly long-ranged intermolecular correlations begin to appear in isotropic (orientationally disordered) phases of liquid crystal forming molecules when the temperature or density starts to close in on the boundary with the nematic (ordered) phase. Indeed, the presence of slowly relaxing, strongly orientationally correlated, sets of molecules under putatively disordered conditions ("pseudo-nematic domains") has been apparent for some time from light-scattering and optical-Kerr experiments. Still, a fully microscopic characterization of these domains has been lacking. We illustrate in this paper how pseudo-nematic domains can be studied in even relatively small computer simulations by looking for order-parameter tensor fluctuations much larger than one would expect from random matrix theory. To develop this idea, we show that random matrix theory offers an exact description of how the probability distribution for liquid-crystal order parameter tensors converges to its macroscopic-system limit. We then illustrate how domain properties can be inferred from finite-size-induced deviations from these random matrix predictions. A straightforward generalization of time-independent random matrix theory also allows us to prove that the analogous random matrix predictions for the time dependence of the order-parameter tensor are similarly exact in the macroscopic limit, and that relaxation behavior of the domains can be seen in the breakdown of the finite-size scaling required by that random-matrix theory.

Glassiness of Thermotropic Liquid Crystals across the Isotropic−Nematic Transition

The orientational dynamics of thermotropic liquid crystals across the isotropic-nematic phase transition have traditionally been investigated at long times or low frequencies using frequency domain measurements. The situation has now changed significantly with the recent report of a series of interesting transient optical Kerr effect (OKE) experiments that probed orientational relaxation of a number of calamitic liquid crystals (which consist of rod-like molecules) directly in the time domain, over a wide time window ranging from subpicoseconds to tens of microseconds. The most intriguing revelation is that the decay of the OKE signal at short to intermediate times (from a few tens of picoseconds to several hundred nanoseconds) follows multiple temporal power laws. Another remarkable feature that has emerged from these OKE measurements is the similarity in the orientational relaxation behavior between the isotropic phase of calamitic liquid crystals near the isotropic-nematic transition and supercooled molecular liquids, notwithstanding their largely different macroscopic states. In this article, we present an overview of the understanding that has emerged from recent computational and theoretical studies of calamitic liquid crystals across the isotropic-nematic transition. Topics discussed include (a) single-particle as well as collective orientational dynamics at a short-to-intermediate time window, (b) heterogeneous dynamics in orientational degrees of freedom diagnosed by a non-Gaussian parameter, (c) fragility, and (d) temperature-dependent exploration of underlying energy landscapes as calamitic liquid crystals settle into increasingly ordered mesophases upon cooling from the high-temperature isotropic phase. A comparison of our results with those of supercooled molecular liquids reveals an array of analogous features in these two important classes of soft matter systems. We further find that the onset of growth of the orientational order in the parent nematic phase induces translational order, resulting in smectic-like layers in the potential energy minima of calamitic systems if the parent nematic phase is sandwiched between the high-temperature isotropic phase and the low-temperature smectic phase. We discuss implications of this startling observation. We also discuss recent results on the orientational dynamics of discotic liquid crystals that are found to be rather similar to those of calamitic liquid crystals.

Orientational relaxation in a discotic liquid crystal

We investigate orientational relaxation of a model discotic liquid crystal, consisting of disclike molecules, by molecular dynamics simulations along two isobars starting from the high temperature isotropic phase. The two isobars have been so chosen that (a) the phase sequence isotropic- (I-) nematic- (N-) columnar (C) appears upon cooling along one of them and (b) the sequence isotropic- (I-) columnar- (C) along the other. While the orientational relaxation in the isotropic phase near the I-N phase transition in system (a) 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 in system (b). In order to understand this difference (the existence or the absence of the power law decay), we calculated the growth of the orientational pair distribution functions (OPDFs) near the I-N phase boundary and also near the I-C phase boundary. We find that the OPDF shows a marked growth in long range correlation as the I-N phase boundary is approached in the I-N-C system (a), but such a growth is absent in the I-C system, which appears to be consistent with the result that I-N phase transition in the former is weakly first order while the I-C phase transition in the latter is not weak. As the system settles into the nematic phase, the decay of the single-particle second-rank orientational time correlation function follows a pattern that is similar to what is observed with calamitic liquid crystals and supercooled molecular liquids.

Comparative study of temperature dependent orientational relaxation in a model thermotropic liquid crystal and in a model supercooled liquid

Recent optical Kerr effect experiments have revealed a power law decay of the measured signal with a temperature independent exponent at short-to-intermediate times for a number of liquid crystals in the isotropic phase near the isotropic-nematic transition and supercooled molecular liquids above the mode coupling theory critical temperature. In this work, the authors investigate the temperature dependence of short-to-intermediate time orientational relaxation in a model thermotropic liquid crystal across the isotropic-nematic transition and in a binary mixture across the supercooled liquid regime in molecular dynamics simulations. The measure of the experimentally observable optical Kerr effect signal is found to follow a power law decay at short-to-intermediate times for both systems in agreement with recent experiments. In addition, the temperature dependence of the power law exponent is found to be rather weak. As the model liquid crystalline system settles into the nematic phase upon cooling, the decay of the single-particle second-rank orientational time correlation function exhibits a pattern that is similar to what has been observed for supercooled liquids.

Structural relaxation in complex liquids: Non-Markovian dynamics in a bistable potential

The time correlation function C(t) identical with <x(0)x(t)> of the distance fluctuations of a particle moving in a bistable potential under the action of fractional Gaussian noise (fGn) is calculated from a Smoluchowski-type equation derived from a generalized Langevin equation (GLE). The time derivative of this function, dC(t)dt, is compared with data from optical Kerr effect measurements of liquid crystal dynamics in the vicinity of the isotropic-to-nematic transition, which are related to the time derivative of an orientational correlation function. A number of characteristic features of the experimental decay curves, including short and intermediate time power law behavior and long time exponential relaxation, are qualitatively reproduced by the analytical calculations, even though the latter do not explicitly treat orientational degrees of freedom. The GLE formalism with fGn was, in fact, originally proposed as a model of protein conformational fluctuations, so the present results suggest that it may also serve more generally as a model of structural relaxation in complex condensed phase media.

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.

Power law relaxation and glassy dynamics in Lebwohl-Lasher model near the isotropic-nematic phase transition

Orientational dynamics in a liquid crystalline system near the isotropic-nematic (IN) phase transition is studied using molecular dynamics simulations of the well-known Lebwohl-Lasher model. As the IN transition temperature is approached from the isotropic side, we find that the decay of the orientational time correlation functions (OTCF) slows down noticeably, giving rise to a power law decay at intermediate time scales. The angular velocity time correlation function also exhibits a rather pronounced power law decay near the IN boundary. In the mean squared angular displacement at comparable time scales, we observe the emergence of a subdiffusive regime which is followed by a superdiffusive regime before the onset of the long-time diffusive behavior. We observe signature of dynamical heterogeneity through pronounced non-Gaussian behavior in orientational motion particularly at lower temperatures. This behavior closely resembles what is usually observed in supercooled liquids. We obtain the free energy as a function of orientational order parameter by the use of the transition matrix Monte Carlo method. The free energy surface is flat for the system considered here and the barrier between isotropic and nematic phases is vanishingly small for this weakly first-order phase transition, hence allowing large scale, collective, and correlated orientational density fluctuations. This might be responsible for the observed power law decay of the OTCFs.

Complete breakdown of the Debye model of rotational relaxation near the isotropic-nematic phase boundary: effects of intermolecular correlations in orientational dynamics

The Debye-Stokes-Einstein (DSE) model of rotational diffusion predicts that the orientational correlation times tau l vary as [l(l+1)]-1, where l is the rank of the orientational time correlation function (given in terms of the Legendre polynomial of rank l). One often finds significant deviation from this prediction, in either direction. In supercooled molecular liquids where the ratio tau 1/tau 2 falls considerably below 3 (the Debye limit), one usually invokes a jump diffusion model to explain the approach of the ratio tau 1/tau 2 to unity. Here we show in a computer simulation study of a standard model system for thermotropic liquid crystals that this ratio becomes much less than unity as the isotropic-nematic phase boundary is approached from the isotropic side. Simultaneously, the ratio tau 2/eta, eta, being the shear viscosity of the liquid, becomes much larger than the hydrodynamic value near the I-N transition. We also analyze the breakdown of the Debye model of rotational diffusion in ratios of higher order orientational correlation times. We show that the breakdown of the DSE model is due to the growth of orientational pair correlation and provide a mode coupling theory analysis to explain the results.

Three homeotropically aligned nematic liquid crystals: Comparison of ultrafast to slow time-scale dynamics

The dynamics of two nematic liquid crystals, 4-(trans-4(')-n-octylcyclohexyl)isothiocyanatobenzene and 4-(4-pentyl-cyclohexyl)-benzonitrile, are investigated as a function of temperature both in the homeotropically aligned nematic phase and in the isotropic phase using optical heterodyne-detected optical Kerr effect experiments, which measures the time derivative of the polarizability-polarizability-correlation function (orientational relaxation). Data are presented over a time range of 500 fs-70 micros for the nematic phase and 500 fs to a few hundred nanoseconds for the isotropic phase. The nematic dynamics are compared with a previously studied liquid crystal in the nematic phase. All three liquid crystals have very similar dynamics in the nematic phase that are very different from the isotropic phase. On the slowest time scale (20 ns-70 micros), a temperature-independent power law, the final power law, t(-f) with f approximately 0.5, is observed. On short time scales (approximately 3 ps to approximately 1 ns), a temperature-dependent intermediate power law is observed with an exponent that displays a linear dependence on the nematic order parameter. Between the intermediate power law and the final power law, there is a crossover region that has an inflection point. For times that are short compared to the intermediate power law (approximately <2 ps), the data decay much faster, and can be described as a third power law, although this functional form is not definitive. The isotopic phase data have the same features as found in previous studies of nematogens in the isotropic phase, i.e., the temperature-independent intermediate power law and von Schweidler power law at short to intermediate times, and a highly temperature-dependent long time exponential decay that is well described by the Landau-de Gennes theory. The results show that liquid-crystal dynamics in the nematic phase exhibit universal behavior.

Orientational Dynamics in the Isotropic Phase of a Calamitic Liquid-Crystal Model

We report a molecular dynamics simulation study on the isotropic phase of an idealized calamitic liquid crystal model with a length-to-width ratio of approximately 5-6. The study focuses on the characterization of single-particle and collective orientational dynamics on approaching the phase transition to the nematic phase. Recent experimental and simulation works have suggested that a power law behavior exists at relatively short times in the decay of the time derivative of the orientational correlation functions. Qualitatively, our simulation data are consistent with these findings. Both single-particle and collective time correlation function derivatives possess, in their respective log-log plots, a linear region at very short times, whose slope is essentially independent from the thermodynamic state. Nevertheless, the single-particle orientational correlation functions are better described by a function which is the sum of a fast exponential, an intermediate stretched-exponential and a slow exponential, while the collective orientational correlation functions are satisfactorily described by a sum of two exponentials, at higher density, or by just one exponential, at lower density.

Universal power law in the orientational relaxation in thermotropic liquid crystals

We observe a surprisingly general power law decay at short to intermediate times in orientational relaxation in a variety of model systems (both calamitic and discotic, and also lattice) for thermotropic liquid crystals. As all these systems transit across the isotropic-nematic phase boundary, two power law relaxation regimes, separated by a plateau, emerge, giving rise to a steplike feature (well known in glassy liquids) in the single-particle second-rank orientational time correlation function. In contrast to its probable dynamical origin in supercooled liquids, we show that the power law here can originate from the thermodynamic fluctuations of the orientational order parameter, driven by the rapid growth in the second-rank orientational correlation length.

Dynamics of the guest-host orientational interaction in dye-doped liquid-crystalline materials

We present a comprehensive study on the dynamics of laser-induced molecular reorientation in a dye-doped liquid crystalline (LC) medium that exhibits significant enhancement of the optical Kerr nonlinearity due to guest-host interaction. Using various techniques, we separately characterized the dynamical responses of the relevant molecular species present in the medium following photoexcitation and, thus, were able to follow the transient process in which photoexcitation of the dye molecules exert through guest-host interaction a net torque on the host LC material, leading to the observed enhanced optical Kerr nonlinearity. Experimental results agree quantitatively with the time-dependent theory based on a mean-field model of the guest-host interaction.

Anomalous glassy relaxation near the isotropic-nematic phase transition

Dynamical heterogeneity in a system of Gay-Berne ellipsoids near its isotropic-nematic (I-N) transition, and also in an equimolar mixture of Lennard-Jones spheres and Gay-Berne ellipsoids in deeply supercooled regime, is probed by the time evolution of non-Gaussian parameters (NGP). The appearance of a dominant second peak in the rotational NGP near the I-N transition signals the growth of pseudonematic domains. Surprisingly, such a second peak is instead observed in the translational NGP for the glassy binary mixture. Localization of orientational motion near the I-N transition is found to be responsible for the observed anomalous orientational relaxation.

Ultrafast to Slow Orientational Dynamics of a Homeotropically Aligned Nematic Liquid Crystal

The orientational dynamics of a homeotropically aligned nematic liquid crystal, 4'-pentyl-4-biphenylcarbonitrile (5-CB), is studied over more than six decades of time (500 fs to 2 mus) using optical heterodyne detected optical Kerr effect experiments. In contrast to the dynamics of nematogens in the isotropic phase, the data do not decay as a highly temperature-dependent exponential on the longest time scale, but rather, a temperature-independent power law spanning more than two decades of time, the final power law, is observed. On short time scales (approximately 3 ps to approximately 1 ns) another power law, the intermediate power law, is observed that is temperature dependent. The power law exponent of the correlation function associated with the intermediate power law displays a linear dependence on the change in the nematic order parameter with temperature. Between the intermediate power law and the final power law, there is a crossover region that displays an inflection point. The temperature-dependent orientational dynamics in the nematic phase are shown to be very different than those observed in the isotropic phase.

Anomalous viscoelasticity near the isotropic-nematic phase transition in liquid crystals

Recent optical Kerr effect experiments have shown that orientational relaxation of nematogens shows a pronounced slow down of the response function at intermediate times and also a power law decay near the isotropic-nematic (I-N) transition. In many aspects, this behavior appears to be rather similar to the ones observed in the supercooled liquid near-glass transition. We have performed molecular dynamics simulations of model nematogens (Gay-Berne with aspect ratio 3) to explore the viscoelasticity near the I-N transition and also investigated the correlation of viscoelasticity (if any) with orientational relaxation. It is found that although the viscosity indeed undergoes a somewhat sharper than normal change near the I-N transition, it is not characterized by any divergence-like behavior (like the ones observed in the supercooled liquid). The rotational friction, on the other hand, shows a much sharper rise as the I-N transition is approached. Interestingly, the probability distribution of the amplitude of the three components of the stress tensor shows anisotropy near the I-N transition-similar anisotropy has also been seen in the deeply supercooled liquid. Frequency dependence of viscosity shows several unusual behaviors: (a) There is a weak, power law dependence on frequency [eta(')(omega) approximately omega(-alpha)] at low frequencies and (b) there is a rapid increase in the sharp peak observed in eta(')(omega) in the intermediate frequency on approach to the I-N transition density. These features can be explained from the stress-stress time correlation function. The angular velocity correlation function also exhibits a power law decay in time. The reason for this is discussed.