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.

Study of liquid crystalline order by NMR

The liquid crystalline phase represents a unique state of matter where partial order exists on molecular and supra-molecular levels and is responsible for several interesting properties observed in this phase. Hence a detailed study of ordering in liquid crystals is of significant scientific and technological interest. NMR provides several parameters that can be used to obtain information about the liquid crystalline phase. Of these, the measurement of dipolar couplings between nuclei has proved to be a convenient way of obtaining liquid crystalline ordering since the coupling is dependent on the average orientation of the dipolar vector in the magnetic field which also aligns the liquid crystal.However, measurement of the dipolar coupling between a pair of selected nuclei is beset with problems that require special solutions. In this article the use of cross polarization for measuring dipolar couplings in liquid crystals is illustrated. Transient oscillations observed during cross polarization provide the dipolar couplings between essentially isolated nearest neighbor spins which can be extracted for several sites simultaneously by employing two-dimensional NMR techniques. The use of the method for obtaining heteronuclear dipolar couplings and hence the order parameters of liquid crystals is presented. Several modifications to the basic experiment are considered and their utility illustrated. A method for obtaining proton–proton dipolar couplings, by utilizing cross polarization from the dipolar reservoir, is presented. Some applications are also highlighted.

Computer simulation of liquid crystals

In this article we review the current status in the modelling of both thermotropic and lyotropic Liquid crystal. We discuss various coarse-graining schemes as well as simulation techniques such as Monte Carlo (MC) and Molecular dynamics (MD) simulations.In the area of MC simulations we discuss in detail the algorithm for simulating hard objects such as spherocylinders of various aspect ratios where excluded volume interaction enters in the simulation through overlap test. We use this technique to study the phase diagram, of a special class of thermotropic liquid crystals namely banana liquid crystals. Next we discuss a coarse-grain model of surfactant molecules and study the self-assembly of the surfactant oligomers using MD simulations. Finally we discuss an atomistically informed coarse-grained description of the lipid molecules used to study the gel to liquid crystalline phase transition in the lipid bilayer system.

Nucleation dynamics and pool boiling characteristics of high pressure refrigerant using thermochromic liquid crystals

This paper reports an experimental investigation of oscillating temperature field beneath a single isolated nucleation site using a non-invasive TLC (thermochromic liquid crystal) based thermography technique. Empirical correlations are presented to demonstrate the influence of system pressure and wall heat flux on different ebullition characteristics in the nucleate pool boiling regime of refrigerant R-134a. TLC transient response and two-phase flow structure are captured using synchronized, high resolution imaging. It is observed that the area of influence of nucleation site exhibits a two-part distinct transient behavior during the bubble growth period and broadens to a maximum of 1.57 times the bubble diameter at the instant of bubble departure. This is accompanied by a sharp fall of 2.5 degrees C in the local excess temperature at the nucleation site, which results in momentary augmentation (similar to 40%) in the local heat transfer coefficient at the nucleation origin. The enhanced heat transfer rate observed during the bubble peel-off event is primarily due to transient micro-convection in the wake of the retreating bubble. Further, the results indicate that a slight increase in system pressure from 813.6 to 882.5 kPa has no considerable effect on either the wall superheat or the overall heat transfer coefficient and ebullition frequency. In addition, correlations have been obtained for bubble Reynolds number, Jackob number and the dimensionless bubble generation frequency in terms of modified boiling number.

Molecular structure of the discotic liquid crystalline phase of hexa-peri-hexabenzocoronene/oligothiophene hybrid and their charge transport properties

Using atomistic molecular dynamics simulation, we study the discotic columnar liquid crystalline (LC) phases formed by a new organic compound having hexa-peri-Hexabenzocoronene (HBC) core with six pendant oligothiophene units recently synthesized by Nan Hu et al. Adv. Mater. 26, 2066 (2014)]. This HBC core based LC phase was shown to have electric field responsive behavior and has important applications in organic electronics. Our simulation results confirm the hexagonal arrangement of columnar LC phase with a lattice spacing consistent with that obtained from small angle X-ray diffraction data. We have also calculated various positional and orientational correlation functions to characterize the ordering of the molecules in the columnar arrangement. The molecules in a column are arranged with an average twist of 25 degrees having an average inter-molecular separation of similar to 5 angstrom. Interestingly, we find an overall tilt angle of 43 degrees between the columnar axis and HBC core. We also simulate the charge transport through this columnar phase and report the numerical value of charge carrier mobility for this liquid crystal phase. The charge carrier mobility is strongly influenced by the twist angle and average spacing of the molecules in the column. (C) 2015 AIP Publishing LLC.

The rα structure, chemical shift anisotropy, and the abnormal orientation of methyldichlorophosphine from the NMR spectrum

NMR studies of methyldichlorophosphine have been undertaken in the nematic phase of mixed liquid crystals of opposite diamagnetic anisotropies. The rα structure is derived. The proton chemical-shift anisotropy has been determined from the studies without the use of a reference compound and without a change of experimental conditions. It is shown that the molecule orients in the liquid crystal with positive diamagnetic anisotropy in such a way that the C3 symmetry axis of the CH3P moiety is preferentially aligned perpendicular to the direction of the magnetic field, unlike other similar systems. This is interpreted in terms of the formation of a weak solvent-solute molecular complex. The heteronuclear indirect spin-spin coupling constants are determined. The sign of the two-bond JPH is found to be positive.

Visualization of enantiomers using natural abundant 13C-filtered single and double quantum selective refocusing experiments: Application to small chiral molecules

The routine use of proton NMR for the visualization of enantiomers, aligned in the chiral liquid crystal solvent poly-γ-benzyl-l-glutamate (PBLG), is restricted due to severe loss of resolution arising from large number of pair wise interaction of nuclear spins. In the present study, we have designed two experimental techniques for their visualization utilizing the natural abundance 13C edited selective refocusing of single quantum (CH-SERF) and double quantum (CH-DQSERF) coherences. The methods achieve chiral discrimination and aid in the simultaneous determination of homonuclear couplings between active and passive spins and heteronuclear couplings between the excited protons and the participating 13C spin. The CH-SERF also overcomes the problem of overlap of central transitions of the methyl selective refocusing (SERF) experiment resulting in better chiral discrimination. Theoretical description of the evolution of magnetization in both the sequences has been discussed using polarization operator formalism.

Deuterium-carbon correlation in static oriented systems using DAPT - A windowless multiple pulse sequence

We present a new method for establishing correlation between deuterium and its attached carbon in a deuterated liquid crystal. The method is based on transfer of polarization using the DAPT pulse sequence proposed originally for two spin half nuclei, now extended to a spin-1 and a spin-1/2 nuclei. DAPT utilizes the evolution of magnetization of the spin pair under two blocks of phase shifted BLEW-12 pulses on one of the spins separated by a 90 degree pulse on the other spin. The method is easy to implement and does not need to satisfy matching conditions unlike the Hartmann-Hahn cross-polarization. Experimental results presented demonstrate the efficacy of the method.

2(n)-SEMA-a robust solid state nuclear magnetic resonance experiment for measuring heteronuclear dipolar couplings in static oriented systems using effective transverse spin-lock

Separated local field (SLF) spectroscopy is a powerful technique to measure heteronuclear dipolar couplings. The method provides site-specific dipolar couplings for oriented samples such as membrane proteins oriented in lipid bilayers and liquid crystals. A majority of the SLF techniques utilize the well-known Polarization Inversion Spin Exchange at Magic Angle (PISEMA) pulse scheme which employs spin exchange at the magic angle under Hartmann-Hahn match. Though PISEMA provides a relatively large scaling factor for the heteronuclear dipolar coupling and a better resolution along the dipolar dimension, it has a few shortcomings. One of the major problems with PISEMA is that the sequence is very much sensitive to proton carrier offset and the measured dipolar coupling changes dramatically with the change in the carrier frequency. The study presented here focuses on modified PISEMA sequences which are relatively insensitive to proton offsets over a large range. In the proposed sequences, the proton magnetization is cycled through two quadrants while the effective field is cycled through either two or four quadrants. The modified sequences have been named as 2(n)-SEMA where n represents the number of quadrants the effective field is cycled through. Experiments carried out on a liquid crystal and a single crystal of a model peptide demonstrate the usefulness of the modified sequences. A systematic study under various offsets and Hartmann-Hahn mismatch conditions has been carried out and the performance is compared with PISEMA under similar conditions.

Geometric quantum computation using fictitious spin-1/2 subspaces of strongly dipolar coupled nuclear spins

Geometric phases have been used in NMR to implement controlled phase shift gates for quantum-information processing, only in weakly coupled systems in which the individual spins can be identified as qubits. In this work, we implement controlled phase shift gates in strongly coupled systems by using nonadiabatic geometric phases, obtained by evolving the magnetization of fictitious spin-1/2 subspaces, over a closed loop on the Bloch sphere. The dynamical phase accumulated during the evolution of the subspaces is refocused by a spin echo pulse sequence and by setting the delay of transition selective pulses such that the evolution under the homonuclear coupling makes a complete 2 pi rotation. A detailed theoretical explanation of nonadiabatic geometric phases in NMR is given by using single transition operators. Controlled phase shift gates, two qubit Deutsch-Jozsa algorithm, and parity algorithm in a qubit-qutrit system have been implemented in various strongly dipolar coupled systems obtained by orienting the molecules in liquid crystal media.

Enantiodiscrimination and extraction of short and long range homo- and hetero-nuclear residual dipolar couplings by a spin selective correlation experiment

A two dimensional correlation experiment for the measurement of short and long range homo- and hetero- nuclear residual dipolar couplings (RDCs) from the broad and featureless proton NMR spectra including C-13 satellites is proposed. The method employs a single natural abundant C-13 spin as a spy nucleus to probe all the coupled protons and permits the determination of RDCs of negligible strengths. The technique has been demonstrated for the study of organic chiral molecules aligned in chiral liquid crystal, where additional challenge is to unravel the overlapped spectrum of enantiomers. The significant advantage of the method is demonstrated in better chiral discrimination using homonuclear RDCs as additional parameters. (C) 2010 Elsevier B.V. All rights reserved.

Visualization of enantiomers and determination of homo- and hetero-nuclear residual dipolar and scalar couplings: The natural abundant C-13 edited J/D-resolved NMR techniques

The simple two dimensional C-13-satellite J/D-resolved experiments have been proposed for the visualization of enantiomers, extraction of homo- and hetero-nuclear residual dipolar couplings and also H-1 chemical shift differences between the enantiomers in the anisotropic medium. The significant advantages of the techniques are in the determination of scalar couplings of bigger organic molecules. The scalar couplings specific to a second abundant spin such as F-19 can be selectively extracted from the severely overlapped spectrum. The methodologies are demonstrated on a chiral molecule aligned in the chiral liquid crystal medium and two different organic molecules in the isotropic solutions. (C) 2010 Elsevier B.V. All rights reserved.

Solvent-solute and solute-solute interactions from NMR in nematic phases

The use of NMR spectroscopy of molecules oriented in liquid-crystalline media to study solvent-solute and solute-solute interactions in π-systems such as benzene-chloroform and in charge transfer complexes, for example pyridineiodine, is illustrated. Changes in molecular order and chemical shifts as a result of complexation are employed in such studies. The extraordinary symmetry of C60 has also been investigated by using a mixture of liquid crystals of opposite diamagnetic anisotropies indicating, thereby, negligible solvent-solute/solute-solute interactions.

Quantification of enantiomeric excess by H-1-detected heteronuclear refocusing and homonuclear multiple quantum NMR experiments

The omega(1)-heterodecoupled-C-13-filtered proton detected NMR experiments are reported for the accurate quantification of enantiomeric excess in chiral molecules embedded in chiral liquid crystal. The differential values of both H-1-H-1 and C-13-H-1 dipolar couplings in the direct dimension and only H-1-H-1 dipolar couplings in the indirect dimension enable unraveling of overlapped enantiomeric peaks. The creation of unequal C-13-bound proton signal for each enantiomer in the INEPT block and non-uniform excitation of coherences in homonuclear multiple quantum experiments do not yield accurate quantification of enantiomeric excess. In circumventing these difficulties, a coupling dependent intensity correction factor has been invoked. (C) 2010 Elsevier B.V. All rights reserved.

Application of a modified z-cosy experiment for the analysis of spectra of oriented molecules - the spectrum of cis,cis-mucanonitrile

NMR spectra of cis,cis-mucanonitrile oriented in a liquid crystal have been analysed using the connectivity information obtained from a modified Z-COSY experiment which provided crucial clues for the starting parameters for the iterative analysis. The proton spectra with and without C-13 satellites and the C-13 spectra have thus been interpreted. The indirect spin-spin couplings required for the analyses have been obtained from the corresponding isotropic spectra. The H-1-H-1 and C-13-H-1 dipolar couplings so obtained have been utilized to determine the relative internuclear distances. The results indicate that the molecule is planar. (C) 1994 Academic Press, Inc.