Microscopic structure and dynamics of a partial bilayer smectic liquid crystal

Comparative study of electrooptic, dielectric, and structural properties of two glassforming antiferroelectric mixtures with a high tilt angle

Vitrification of the antiferroelectric smectic-C_{A}^{*} phase is reported for the orthoconic mixture W-1000 and its new derivative W-356. The crystallization is not observed even upon slow cooling and the cold crystallization on subsequent heating is also absent. Molecular dynamics in the smectic phases of both mixtures is investigated down to 173 K and the fragility parameters are determined from the temperature behavior of the α-process. X-ray diffraction is applied to compare the structural parameters of W-356 and W-1000 as well as to study the structural changes during the glass transition of the Sm-C_{A}^{*} phase. The evolution of the smectic layer order within the Sm-C_{A}^{*} glass is reported, while the correlation length of the short-range order in the smectic layers is shown to be almost constant below the glass transition temperature. Electrooptic properties of W-356: spontaneous polarization, tilt angle and switching time are determined and compared with these of W-1000. The observed differences between the properties of W-356 and W-1000 might be explained by the dimer formation of components with the -C≡N terminal group, present only in the W-356 mixture, and by the different structure of the aromatic molecular core in one of the W-356 components.

Shape segregation in molecular organisation: a combined X-ray scattering and molecular dynamics study of smectic liquid crystals

Temperature-dependent X-ray scattering studies have been carried out on 4-undecyloxy-4'-cyanobiphenyl (11OCB) and 4-(12,12-dimethyltridecyloxy)-4'-cyanobiphenyl (t-Bu-11OCB) in the smectic A phase, from which their layer spacings and translational order parameters were obtained. Marked differences between the layer structures of the two compounds were demonstrated, showing that the addition of the t-Bu group results in a higher translational order parameter and wider layer spacing for t-Bu-11OCB than 11OCB. Fully atomistic MD simulations of both compounds run for >1000 ns demonstrated the spontaneous formation of smectic mesophases from isotropic starting geometries, and experimental trends in order parameters and absolute layer spacings were shown to be replicated well. Further analysis showed that both the aromatic head-groups and the alkyl tail-groups exhibit interdigitation in the simulated smectic phases of both compounds, and the difference in the layer structures between 11OCB and t-Bu-11OCB could be attributed mainly to a shape segregation effect arising from the addition of the bulky t-Bu end-group to the alkyl chain.

Simulation insights into the role of antiparallel molecular association in the formation of smectic A phases

A simple dissipative particle dynamics (DPD) model is introduced, which can be used to represent a broad range of calamitic mesogens. The model allows for antiparallel association that occurs naturally in a number of mesogens with terminal dipoles, including the 4-n-alkyl-4'-cyanobiphenyl (nCB) series. Favourable antiparallel interactions lead to the formation of SmA_d phases in which the layer spacing is intermediate between monolayer and bilayer. The model is easily tuned to vary the strength of antiparallel association and the SmA layer spacing, and to give either isotropic-smectic or isotropic-nematic-smectic phase sequences. The model allows for a range of other smectics: including SmA₁ phases exhibiting microphase separation within layers, and smectics A structures with more complicated repeat units. For large system sizes (≥50 000 molecules) in the nematic phase, we are able to demonstrate the formation of three distinct types of cybotactic domains depending on the local interactions. Cybotactic domains are found to grow in the nematic-smectic pretransitional region as the system moves closer to T_SN.

Evaluations of Mesogen Orientation in Thin Films of Polyacrylate with Cyanobiphenyl Side Chain

The orientation behavior of mesogens in a polyacrylate with cyanobiphenyl (CB) side chain in thin films was investigated in detail by UV-vis absorption spectroscopy and grazing incidence small-angle X-ray scattering (GI-SAXS) measurements using both high-energy X-rays of Cu Kα line (λ = 0.154 nm) and low-energy synchrotron X-rays (λ = 0.539 nm). By changing the film thickness ranging 7-200 nm, it is concluded that the planar orientation is predominant for thin films with thickness below 10-15 nm. This planar mesogen orientation near the substrate surface coexists with the homeotropically aligned CB mesogens in films thicker than 30 nm. For the thinnest 7 nm film, the planar orientation is unexpectedly lost, which is in consort with a disordering of smectic layer structure. Peculiar orienting characteristics of CB mesogen are suggested, which probably stem from the tendency to form an antiparallel arrangement of mesogens due to the strong dipole moment of the terminal cyano group.

Induced smectic phases of stoichiometric liquid crystal mixtures

We revealed the detailed structures of induced smectic liquid crystal (LC) phases composed of a binary mixture of charge-transfer (CT) LC substances. Although neither of the constituents had highly ordered smectic phases, the mixture exhibited smectic-E (SmE) or smectic-B (SmB) phases when mixed at ratios of 1 : 1 and 2 : 3, respectively. The results of polarized optical microscopy, differential scanning calorimetry, X-ray diffraction, and infrared spectroscopy indicated that the induced smectic phases were stabilized by an exquisite balance between the CT interactions, dipolar interactions, and excluded volume effects. We proposed a possible model for the molecular arrangements in the SmE and SmB phases, which consistently explained the experimental results including the stoichiometric ratios.

Single-molecule diffusion in freely suspended smectic films

We present a study of the molecular diffusion in freely suspended smectic-A liquid crystal films with thicknesses ranging from 20 down to only two molecular layers. The molecular mobility is directly probed by determining the trajectories of single, fluorescent tracer molecules. We demonstrate, using several different smectic compounds, that a monotonic increase of the diffusion coefficient with decreasing film thickness is a general phenomenon. In two-layer films, the diffusion is enhanced by a factor of 3 to 5 compared to the corresponding bulk smectic phase. Molecular dynamics simulations of freely suspended smectic films are presented which support the experimental results.

Chiral heliconical ground state of nanoscale pitch in a nematic liquid crystal of achiral molecular dimers

Freeze-fracture transmission electron microscopy study of the nanoscale structure of the so-called "twist-bend" nematic phase of the cyanobiphenyl (CB) dimer molecule CB(CH2)7CB reveals stripe-textured fracture planes that indicate fluid layers periodically arrayed in the bulk with a spacing of d ~ 8.3 nm. Fluidity and a rigorously maintained spacing result in long-range-ordered 3D focal conic domains. Absence of a lamellar X-ray reflection at wavevector q ~ 2π/d or its harmonics in synchrotron-based scattering experiments indicates that this periodic structure is achieved with no detectable associated modulation of the electron density, and thus has nematic rather than smectic molecular ordering. A search for periodic ordering with d ~ in CB(CH2)7CB using atomistic molecular dynamic computer simulation yields an equilibrium heliconical ground state, exhibiting nematic twist and bend, of the sort first proposed by Meyer, and envisioned in systems of bent molecules by Dozov and Memmer. We measure the director cone angle to be θ(TB) ~ 25° and the full pitch of the director helix to be p(TB) ~ 8.3 nm, a very small value indicating the strong coupling of molecular bend to director bend.

Smectic order parameters via liquid crystal NMR spectroscopy: Application to a partial bilayer smectic A phase

Solute molecules were dissolved in the liquid crystal 4-cyano-4'-n-octyloxybiphenyl (8OCB), known to form a partial bilayer smectic-A phase. Through measurement of solutes' and solvent's orientational order parameters via nuclear magnetic resonance spectroscopy, and their analysis via a statistical thermodynamic density functional theory, values of the solvent's positional order parameters and solutes' positional-orientational distribution functions were obtained. Near to the transition to the nematic phase, the main positional order parameter of the smectic liquid crystal turned out to be comprised in the interval 0.4-0.6, though the quality of the fittings assuming the phase as nematic all across the temperature range investigated was only slightly worse. This may be ascribed to the looseness of the partial bilayer smectic structure. Solutes were found to preferentially lie in those regions where liquid crystal molecule terminal chains are located.

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.

Atomistic simulation of a model liquid crystal

We present atomistic molecular dynamics computer simulations of the bulk phases of a model liquid crystal system based on 8CB. The model differs from real 8CB because it employs a united-atom description to eliminate all hydrogen atoms, and neglects all long-range electrostatic interactions. Despite this simplification, the pressure-temperature phase diagram shows an order-disorder transition, in which isotropic, smectic, and nematiclike behaviors are observed. A detailed analysis of the inter- and intramolecular structures of the ordered phases is given, together with an examination of finite size effects and the equilibration times of the system. It is shown that, whereas a system may appear to be thermodynamically and mechanically equilibrated after a period of 10-15 ns, it is possible for an imprint of the starting configuration to persist for much longer time scales. In the present case, however, such an imprint does not appear to affect the observed phase behavior.

Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations

Recent experiments have shown that liquid crystals can be used to image mammalian cell membranes and to amplify structural reorganization in phospholipid-laden liquid crystal-aqueous interfaces. In this work, molecular dynamics simulations were employed to explore the interactions between commonly used liquid crystal-forming molecules and phospholipid bilayers. In particular, umbrella sampling was used to obtain the potential of mean force of 4-cyano-4'-pentylbiphenyl (5CB) and 4'-(3,4-difluor-phenyl)-4-pentyl-bicylohexyl (5CF) molecules partitioning into a dipalmitoylphosphatidylcholine bilayer. In addition, results of simulations are presented for systems consisting of a fully hydrated bilayer with 5CB or 5CF molecules at the lowest (4.5 mol %) and highest (20 mol %) concentrations used in recent laboratory experiments. It is found that mesogens preferentially partition from the aqueous phase into the membrane; the potential of mean force exhibits highly favorable free energy differences for partitioning (-18 k(B)T for 5CB and -26 k(B)T for 5CF). The location and orientation of mesogens associated with the most stable free energies in umbrella sampling simulations of dilute systems were found to be consistent with those observed in liquid-crystal-rich bilayers. It is found that the presence of mesogens in the bilayer enhances the order of lipid acyl tails, and changes the spatial and orientational arrangement of lipid headgroup atoms. These effects are more pronounced at higher liquid-crystal concentrations. In comparing the behavior of 5CB and 5CF, a stronger spatial correlation (i.e., possibly leading to aggregation) is observed between 5CB molecules within a bilayer than between 5CF molecules. Also, the range of molecular orientations and positions along the bilayer normal is larger for 5CB molecules. At the same time, 5CF molecules were found to bind more strongly to lipid headgroups, thereby slowing the lateral motion of lipid molecules.

Atomistic Simulation of a Nematogen Using a Force Field Derived from Quantum Chemical Calculations

Bulk phase atomistic computer simulations of 4-n-pentyl-4'-cyanobiphenyl (5CB) were performed with a specific force field obtained from ab initio and DFT calculations. The intermolecular potential was previously derived through the fragmentation reconstruction method (FRM), developed in our group. The description of some intramolecular interactions, like the torsional potential between the phenyl rings and at the aryl-alkyl linkage, is achieved through accurate DFT studies. Lengthy ( approximately 40 ns) molecular dynamics (MD) simulations were then carried out at constant atmospheric pressure and different temperatures. The system was stable in the experimental crystalline structure up to 285 K, where the early stage of the melting process appears with the loss of positional order. At higher temperatures (between 290 and 305 K) a kinetically stable, orientationally ordered phase is obtained. This nematic phase was reached starting with three initial configurations, differing in their orientational order parameter. The calculated values of thermodynamic and structural properties of each phase were in fairly good agreement with the relevant experimental data.

Anisotropic self-diffusion in thermotropic liquid crystals studied by 1H and 2H pulse-field-gradient spin-echo NMR

The molecular self-diffusion coefficients in nematic and smectic-A thermotropic liquid crystals are measured using stimulated-echo-type 2H and 1H pulse-field-gradient spin-echo nuclear magnetic resonance (PGSE NMR) combined with multiple-pulse dipolar decoupling and slice selection. The temperature dependence of the principal components of the diffusion tensor in the nematic phase follows a simple Arrhenius relationship except in the region of nematic-isotropic phase transition where it reflects, merely, the decrease of the molecular orientational order. The average of the principal diffusion coefficients in the isotropic-nematic phase transition region is close to the diffusion coefficient in the isotropic phase. At the nematic-smectic-A phase transition the diffusion coefficients change continuously. The results in nematic phase are best described in terms of the affine transformation model for diffusion in nematics formed by hard ellipsoids. In the smectic-A phase the data are interpreted using a modified model for diffusion in presence of a periodic potential along the director.