Molecular orientational and dipolar correlation in the liquid crystal mixture E7: a molecular dynamics simulation study at a fully atomistic level

Molecular dynamics simulations are reported for the four component nematic liquid crystal mixture E7, which is used commercially. We are able to show the growth of a nematic phase directly from an isotropic liquid over a 100 ns period for an all-atom model, and study orientational and dipole order within the nematic phase. The simulations show that the cyanoterphenyl component of the mixture, 5CT, is more ordered than the three cyanobiphenyl components. The simulations show also that both parallel and anti-parallel dipole correlation take place in E7 but that the strong anti-parallel dipole correlation is localised to particular arrangements of molecules. It is possible to identify two key preferred configurations for molecular pairs in the fluid, which explain the form of the dipole correlation function, g(1)(r).     

Simulation studies of dipole correlation in the isotropic liquid phase

The Kirkwood correlation factor g₁ determines the preference for local parallel or antiparallel dipole association in the isotropic phase. Calamitic mesogens with longitudinal dipole moments and Kirkwood factors greater than 1 have an enhanced effective dipole moment along the molecular long axis. This leads to higher values of Δ? in the nematic phase. This paper describes state-of-the-art molecular dynamics simulations of two calamitic mesogens 4-(trans-4-n-pentylcyclohexyl)benzonitrile (PCH5) and 4-(trans-4-n-pentylcyclohexyl)chlorobenzene (PCH5-Cl) in the isotropic liquid phase using an all-atom force field and taking long range electrostatics into account using an Ewald summation. Using this methodology, PCH5 is seen to prefer antiparallel dipole alignment with a negative g₁ and PCH5-Cl is seen to prefer parallel dipole alignment with a positive g₁; this is in accordance with experimental dielectric measurements. Analysis of the molecular dynamics trajectories allows an assessment of why these molecules behave differently.     

Influence of electrostatic interactions on the properties of cyanobiphenyl liquid crystals predicted from atomistic molecular dynamics simulations

The influence of force field details in all-atom molecular dynamics (MD) simulations on the predicted thermodynamic, structural, and dynamic properties of bulk 4-cyano-4?-pentylbiphenyl (5CB) systems have been investigated in the 292–368 K temperature range. The effect of the molecular dipole moment and the details of dihedral potential for biphenyl unit were investigated using both polarisable (POL) and non-polarisable (NP) versions of the quantum chemistry-based force field. The predicted densities for the nematic and isotropic phases of bulk 5CB were found to be in excellent agreement with available experimental data. The nematic-isotropic transition temperature (T_NI) showed strong sensitivity to the force field details, MD simulations with partial atomic charge distributions and molecular dipole moment corresponding to high-level quantum chemistry calculations predicted an overestimation of the T_NI by about 30 K. Rescaling the charges to allow the molecular dipole to be closer to experimentally reported values of 5CB dipole in condensed phases, significantly improved the prediction of T_NI as well as other thermodynamic and dynamic properties of 5CB. We also discuss how the structural, thermodynamic, and dynamic properties of bulk 5CB are affected by the flexibility of the central biphenyl dihedral and the inclusion of induced polarisation effects.     

Molecular dynamics simulation study on the phase behavior of the Gay-Berne model with a terminal dipole and a flexible tail

To study the effect of the alkyl tail and the terminal dipole on the stability of the liquid crystalline phase of mesogens, we have carried out molecular dynamics simulations for 1CB(4-methyl-4'-cyanobiphenyl) and 5CB(4-n-pentyl-4'-cyanobiphenyl) by using a coarse-grained model. In the coarse-grained model, a 5CB molecule is divided into the rigid part of 1CB moiety, which is represented by an ellipsoid, and the remaining flexible part which is represented by a chain of united atoms. The nonbonded potential between coarse-grained segments is represented by the generalized Gay-Berne (GB) potential and the potential parameters are determined by directly comparing the GB potential with the atomistic potentials averaged over the rotation of the mesogen around its axis. In addition, a dipole moment is placed at one end of the ellipsoid opposite to the flexible tail. The ordered state obtained in the polar 5CB model was assigned as the nematic phase, and the experimental static and dynamical properties were reproduced well by using this coarse-grained model. Both the dipole-dipole interactions and the thermal fluctuation of the flexible tail increase the positional disorder in the director direction, and stabilize the nematic phase. Thus, the nematic phase in the polar 5CB is induced by a cooperative effect of the flexible tail and the terminal dipole. It is noted that a local bilayer structure with head-to-head association is formed in the nematic phase, as experimentally observed by x-ray diffraction measurements.     

Influence of the anisometry of magnetic particles on the isotropic–nematic phase transition

The influence of the shape anisotropy of magnetic particles on the isotropic–nematic phase transition was studied in ferronematics based on the nematic liquid crystal (LC) 4-(trans-4-n-hexylcyclohexyl)-isothiocyanato-benzene (6CHBT). The LC was doped with spherical or rod-like magnetic particles of different size and volume concentrations. The phase transition from isotropic to nematic phase was observed by polarising microscope as well as by capacitance measurements. The influence of the concentration and the shape anisotropy of the magnetic particles on the isotropic–nematic phase transition in LC are demonstrated here. The results are in a good agreement with recent theoretical predictions.     

Dielectric and viscous properties of 6CHBT in the isotropic and nematic phases

This paper presents the results of measurements of the principal electric permittivities ε ∥_*(T, ω and ε⊥*(T, ω) and the viscosity for 4-(trans-4-n-hexylcyclohexyl)isothiocyanatobenzene (C₆H₁₃-CyHx-Bz-N=C=S, 6CHBT). In the nematic phase, the Miesowicz η₂ viscosity coefficient was measured in a sample oriented due to the flow. On the basis of the temperature dependence of the static permittivities, using the Maier-Meier equations, the angle β between the dipole moment vector and the long axis of the 6CHBT molecule, the square of the molecular apparent dipole moment μ² _app ² and the nematic order parameter S(T), were determined. From the temperature dependence of the viscosity and the relaxation time corresponding to the molecular rotation around the short axis, the strength of the nematic potential and the effective length of the 6CHBT molecule (in the isotropic phase) were estimated.     

Creation of a nematic phase in mixtures of smectic A1 phases

Phase diagrams were determined for binary mixtures consisting of two 5-n-alkyl-2-(4'-isothiocyanatophenyl)-1,3-dioxane compounds (n-DBT) or 4'-isothio-cyanatophenyl 4-(trans-4'-n-decylcyclohexyl)benzoate and n-DBT. All compounds investigated have monolayer smectic A phases. A nematic phase in the upper temperature range and a nematic gap between two smectic regions also were observed, with the smectic layer spacing ratio, d/d', of 1.23 and 1.87 respectively. The variation of the enthalpy of transition with mixture composition in relation to changes of layer spacing ratio are also discussed for these systems.     

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.     

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.     

Molecular theory of dielectric relaxation in nematic dimers

This paper reports a theory for the dielectric relaxation of dimeric mesogenic molecules in a nematic liquid crystal phase. Liquid crystal dimers consist of two mesogenic groups linked by a flexible chain. Recent experimental studies [D. A. Dunmur, G. R. Luckhurst, M. R. de la Fuente, S. Diez, and M. A. Perez Jubindo, J. Chem. Phys. 115, 8681 (2001)] of the dielectric properties of polar liquid crystal dimers have found unexpected results for both the static (low frequency) and variable frequency dielectric response of these materials. The theory developed in this paper provides a quantitative model with which to understand the observed experimental results. The mean-square dipole moments of alpha,omega-bis[(4-cyanobiphenyl-4'-yl]alkanes in a nematic phase have been calculated using both the rotational isomeric state model and a full torsional potential for the carbon-carbon bonds of the flexible chain. The orienting effect of the nematic phase is taken into account by a parametrized potential of mean torque acting on the mesogenic groups and the segments in the flexible chain. Results of calculations using the full torsional potential are in excellent agreement with experimental results for comparable systems. The probability density p(eq)(beta(A),beta(B)) for the orientation of the mesogenic groups (A,B) along the nematic director is also calculated. The resultant potential of mean torque is a surface characterized by four deep energy wells or sites equivalent to alignment of the terminal groups A and B approximately parallel and antiparallel to the director; of course, the reversal of the director leads to equivalent sites. This potential energy surface provides the basis for a kinetic model of dielectric relaxation in nematic dimers. Solution of the Fokker-Planck equation corresponding to this four-site model gives the time dependence of the site populations, and hence the time-correlation functions for the total dipole moment along the director. In this model the end-over-end rotation of the molecule, corresponding to simultaneous reversal of both mesogenic groups, is excluded because the activation energy is too large. Results are presented for a number of cases, in which a dipole is located on one or both of the mesogenic groups, and additionally where the groups differ in size. For the latter, under particular conditions, the correlation function exhibits a biexponential decay, which corresponds to two low frequency absorptions in the dielectric spectrum. This is exactly what has been observed for nonsymmetric nematic dimers having different groups terminating a flexible chain. Experimental results over a range of temperature for the nonsymmetric dimer alpha-[(4-cyanobiphenyl)-4'-yloxy]-omega-(4-decylanilinebenzylidene-4'-oxy)nonane can be fitted precisely to the theory, which provides new insight into the orientational and conformational dynamics of molecules in ordered liquid crystalline phases.     

Localization or delocalization in the electronic structure of Creutz-Taube-type complexes in aqueous solution

Creutz-Taube complex, [(NH3)5Ru-pyrazine-Ru(NH3)5]5+ (1), and its analogues, [(NH3)5Os-pyrazine-Os(NH3)5]5+ (2), [(NH3)5Ru(4,4'-bipyridine)Ru(NH3)5]5+ (3), and [(NH3)5Os(4,4'-bipyridine)Os(NH3)5]5+ (4), were theoretically investigated by the combination of a two-state model and the dielectric continuum model. Their electronic structures are very sensitive to the metal, ligand, and solvent. In the gas phase, the electronic structures of 1-4 would be completely delocalized. In aqueous solution, that of 3 becomes localized because the polar solvent stabilizes the localized electronic structure with the large dipole moment. However, 1 and 2 are still delocalized in aqueous solution. In 4, the electronic structure would be localized when the dihedral angle between two pyridyl rings is 80 degrees , while it would become delocalized when the angle is small. The origins of the difference are the smaller overlap integral and larger energy difference between two diabatic states, of which electronic structure is almost localized on each metal center.     

Influence of trans-cis photoisomerizations of solubilized compounds on transition temperatures in the system cetyltrimethylammonium bromide-water: Light-induced phase transitions

Transition temperatures (TN1) from the nematic lyotropic liquid-crystalline phase to the isotropic phase were measured for the system cetyltrimethylammonium bromide (CTAB) water in the presence of small amounts of 3-stilbene carboxylic acid (3SC), 4-stilbene carboxylic acid (4SC) and Δ2/2'-bi-(2H-l,4-benzothiazine) (BT). TNI, increases as a function of trans-3SC or trans-4SC concentration, ranging from 01 to 08 wt %, by up to 12°C. A further increase in TNI between 2 and 5°C can be achieved by photochemically converting the solubilized trans stilbene derivatives to the cis isomers. Irradiation of a trans-3SC containing sample at a temperature just above TNl leads to a light-induced phase transition to the lyotropic liquid-crystalline phase. Solubilization of trans-BT causes a slight decrease of TNI while photoisomerization to cis-BT increases T_NI by 1°C.     

Condensation of free volume in structures of nematic and hexatic liquid crystals

Eight novel liquid crystalline materials were prepared containing highly branched terminal chains, either 2,4,4-trimethylpentyl or 3,5,5-trimethylhexyl. All materials exhibit nematic mesophases, with additional smectic (Sm) C, hexatic B and SmI phases for certain homologues. Analysis by small- and wide-angle X-ray scattering reveals continual build-up of the correlation length within the nematic phases, where we also observe splitting of the small angle peak into four lobes, indicating pretransitional Sm fluctuations. Connoscopy confirms the nematic phase to be uniaxial and optically positive. We observe that in the solid state, the molecules exist as staggered antiparallel pairs as a consequence of the sterically demanding bulky terminal group, and this would also appear to manifest in the hexatic B phase, where the layer spacing was found to be greater than the molecular length. If true, this is an example of pair formation driven by sterics rather than dipole–dipole interactions and suggests that reentrant systems driven purely by steric frustration may be found.     

Broadband dielectric studies of weakly polar and non-polar liquid crystals

Complex dielectric permittivities, for two orientations of the director n, parallel (E || n) and perpendicular (E ⊥ n) to the probing electric field E, of the weakly polar liquid crystals (LCs) 4,4'-dihexylazoxybenzene (D6AOB) and 4,4'-diheptylazoxybenzene (D7AOB) as well as the non-polar LC diheptylazobenzene (D7AB) have been measured in the frequency range 75 kHz to 1 GHz. The measurements were performed in the nematic, smectic and isotropic phases of the LCs. The dielectric anisotropies Δε (=ε_||-ε_⊥) obtained from the values of dielectric permittivities at 100 kHz in the nematic phase were found to increase with decreasing temperature. However, for the DnAOBs, the Δε values are somewhat smaller than that for D7AB which does not have a permanent dipole moment. In the nematic phase two molecular relaxation processes were observed for both DnAOBs in each of the orientations—parallel and perpendicular. The four processes merge into two separate processes in the isotropic phase. For D7AB no orientational relaxations were observed in the experimental frequency range.     

Molecular field theory for polar,biaxial bent-core nematics

We have studied a polar, biaxial nematic liquid crystal formed from bent-core molecules using molecular field theory. The model includes a simple Heisenberg-form dipolar intermolecular interaction in addition to the usual quadrupolar nematic interaction, and mimics a system consisting of nematogenic bent-core molecules with a large transverse dipole along the bisector of the two molecular arms. Such systems are regarded as good candidates for biaxial nematic liquid crystals. In principle, the molecular dipoles can align, thus stabilising the ordering of the minor axes. Our calculations predict that, for suitable values of the bent-core interarm angle, the biaxial nematic phase can be stabilised at higher temperatures than in the absence of the transverse dipole. In general, the transverse macroscopic polar order stabilises the biaxial nematic phase. In particular, for a large enough dipolar interaction, the Landau point in the pure biaxial nematic develops into a line of first-order polar biaxial nematic-to-isotropic phase transitions.     

Effect of properties of the smectic Ad phase on the induction of a nematic phase in binary mixtures of smectics A1 and Ad

The phase diagrams of two series of binary mixtures composed of 5-trans-n-butyl-2-(4-isothiocyanatophenyl)-1,3-dioxane (smectic A₁) and 4-(trans-4'-n-alkyl-cyclohexyl)-1-(2-cyano-ethenyl)benzene or 4-cyanophenyl-4'-n-alkoxybenzoates (both smectic A_d) are determined. It is shown that the smectic layer spacing ratio and smectic phase transition enthalpy are the factors influencing the width and the position of the nematic gap created between the two smectic regions. It is found that their influence on this ability is in the opposite direction.     

Dipole moments and dielectric properties of fluorine substituted nematic liquid crystals

The effect of lateral fluorine substituents on the dielectric properties of a range of nematic liquid crystals is reported. Measurements of dipole moments and electric permittivities have been made and used to calculate Kirkwood correlation factors over a range of temperature. The results show that the extent and nature of dipole correlations are important considerations in determining the magnitude of the dielectric anisotropy in these materials. Correlation factors up to 1.4 are found, providing good evidence for the presence of parallel dipole associations.     

Creation of a nematic phase in mixtures of smectic A1 phases

Abstract

Phase diagrams were determined for binary mixtures consisting of two 5-n-alkyl-2-(4'-isothiocyanatophenyl)-1,3-dioxane compounds (n-DBT) or 4'-isothio-cyanatophenyl 4-(trans-4'-n-decylcyclohexyl)benzoate and n-DBT. All compounds investigated have monolayer smectic A phases. A nematic phase in the upper temperature range and a nematic gap between two smectic regions also were observed, with the smectic layer spacing ratio, d/d', of 1.23 and 1.87 respectively. The variation of the enthalpy of transition with mixture composition in relation to changes of layer spacing ratio are also discussed for these systems.     

Monte Carlo simulations on mesophase formation using dipolar Gay-Berne model

We report the results of a Monte Carlo simulation of polar particles interacting via the Gay-Berne potential combining dipole-dipole interactions. Simulations were carried out on a system of 256 particles with either a zero dipole moment or longitudinal dipole moment located at the centre of the molecule. The system was found to spontaneously form nematic, smectic and crystal phases from an isotropic phase with a random configuration as temperature was decreased, irrespective of values of the dipole moment. The results do not give any indication of a net polarization even in the system with a strong dipole moment (μ* = 2.00). The transition temperature from the isotropic to nematic phase is not sensitive to the value of the dipole moment within the limits of statistical error, while the transition from the nematic to smectic phase depends on the strength of dipole moment. At lower temperatures forming the smectic or the crystal phase, the translational order along the director increases with increasing dipole moment. The dipolar interactions contribute to the long range ordering.     

The mesomorphic anomaly of 4-[2-(4-alkoxy-2,3,5,6-tetrafluorophenyl)ethynyl] phenyl trans-4-pentylcyclohexyl-1-carboxylates

Four homologous series of 4-[2-[4-alkoxy-2,3,5,6-tetrafluorophenyl)ethynyl] phenyl trans-4-alkylcyclohexyl-1-carboxylates have been prepared. Their liquid crystalline behaviour was investigated by optical polarizing microscopy and DSC. Series A, B and C exhibit the nematic phase. The phase transition behaviour of series D is anomalous. The compounds of series D with a short alkoxy chain exhibit enantiotropic smectic A and nematic phases, while those with a long alkoxy chain exhibit only the nematic phase.