Tricritical points in biaxial liquid crystal phases

Complex free-energy landscapes in biaxial nematic liquid crystals and the role of repulsive interactions: A Wang-Landau study

General quadratic Hamiltonian models, describing the interaction between liquid-crystal molecules (typically with D_{2h} symmetry), take into account couplings between their uniaxial and biaxial tensors. While the attractive contributions arising from interactions between similar tensors of the participating molecules provide for eventual condensation of the respective orders at suitably low temperatures, the role of cross coupling between unlike tensors is not fully appreciated. Our recent study with an advanced Monte Carlo technique (entropic sampling) showed clearly the increasing relevance of this cross term in determining the phase diagram (contravening in some regions of model parameter space), the predictions of mean-field theory, and standard Monte Carlo simulation results. In this context, we investigated the phase diagrams and the nature of the phases therein on two trajectories in the parameter space: one is a line in the interior region of biaxial stability believed to be representative of the real systems, and the second is the extensively investigated parabolic path resulting from the London dispersion approximation. In both cases, we find the destabilizing effect of increased cross-coupling interactions, which invariably result in the formation of local biaxial organizations inhomogeneously distributed. This manifests as a small, but unmistakable, contribution of biaxial order in the uniaxial phase. The free-energy profiles computed in the present study as a function of the two dominant order parameters indicate complex landscapes. On the one hand, these profiles account for the unusual thermal behavior of the biaxial order parameter under significant destabilizing influence from the cross terms. On the other, they also allude to the possibility that in real systems, these complexities might indeed be inhibiting the formation of a low-temperature biaxial order itself-perhaps reflecting the difficulties in their ready realization in the laboratory.

Hierarchy of orientational phases and axial anisotropies in the gauge theoretical description of generalized nematic liquid crystals

The paradigm of spontaneous symmetry breaking encompasses the breaking of the rotational symmetries O(3) of isotropic space to a discrete subgroup, i.e., a three-dimensional point group. The subgroups form a rich hierarchy and allow for many different phases of matter with orientational order. Such spontaneous symmetry breaking occurs in nematic liquid crystals, and a highlight of such anisotropic liquids is the uniaxial and biaxial nematics. Generalizing the familiar uniaxial and biaxial nematics to phases characterized by an arbitrary point-group symmetry, referred to as generalized nematics, leads to a large hierarchy of phases and possible orientational phase transitions. We discuss how a particular class of nematic phase transitions related to axial point groups can be efficiently captured within a recently proposed gauge theoretical formulation of generalized nematics [K. Liu, J. Nissinen, R.-J. Slager, K. Wu, and J. Zaanen, Phys. Rev. X 6, 041025 (2016)2160-330810.1103/PhysRevX.6.041025]. These transitions can be introduced in the model by considering anisotropic couplings that do not break any additional symmetries. By and large this generalizes the well-known uniaxial-biaxial nematic phase transition to any arbitrary axial point group in three dimensions. We find in particular that the generalized axial transitions are distinguished by two types of phase diagrams with intermediate vestigial orientational phases and that the window of the vestigial phase is intimately related to the amount of symmetry of the defining point group due to inherently growing fluctuations of the order parameter. This might explain the stability of the observed uniaxial-biaxial phases as compared to the yet to be observed other possible forms of generalized nematic order with higher point-group symmetries.

Reexamination of the mean-field phase diagram of biaxial nematic liquid crystals: Insights from Monte Carlo studies

Investigations of the phase diagram of biaxial liquid-crystal systems through analyses of general Hamiltonian models within the simplifications of mean-field theory (MFT), as well as by computer simulations based on microscopic models, are directed toward an appreciation of the role of the underlying molecular-level interactions to facilitate its spontaneous condensation into a nematic phase with biaxial symmetry. Continuing experimental challenges in realizing such a system unambiguously, despite encouraging predictions from MFT, for example, are requiring more versatile simulational methodologies capable of providing insights into possible hindering barriers within the system, typically gleaned through its free-energy dependences on relevant observables as the system is driven through the transitions. The recent paper from this group [Kamala Latha et al., Phys. Rev. E 89, 050501(R) (2014)], summarizing the outcome of detailed Monte Carlo simulations carried out employing an entropic sampling technique, suggested a qualitative modification of the MFT phase diagram as the Hamiltonian is asymptotically driven toward the so-called partly repulsive regions. It was argued that the degree of (cross) coupling between the uniaxial and biaxial tensor components of neighboring molecules plays a crucial role in facilitating a ready condensation of the biaxial phase, suggesting that this could be a plausible factor in explaining the experimental difficulties. In this paper, we elaborate this point further, providing additional evidence from curious variations of free-energy profiles with respect to the relevant orientational order parameters, at different temperatures bracketing the phase transitions.

Detection of an intermediate biaxial phase in the phase diagram of biaxial liquid crystals: entropic sampling study

We investigate the phase sequence of biaxial liquid crystals, based on a general quadratic model Hamiltonian over the relevant parameter space, with a Monte Carlo simulation which constructs equilibrium ensembles of microstates, overcoming possible (free) energy barriers (combining entropic and frontier sampling techniques). The resulting phase diagram qualitatively differs from the universal phase diagram predicted earlier from mean-field theory (MFT), as well as the Monte Carlo simulations with the Metropolis algorithm. The direct isotropic-to-biaxial transition predicted by the MFT is replaced in certain regions of the space by the onset of an additional intermediate biaxial phase of very low order, leading to the sequence N(B)-N(B1)-I. This is due to inherent barriers to fluctuations of the components comprising the total energy, and may explain the difficulties in the experimental realization of these phases.

Molecular field theory for biaxial smectic A liquid crystals

Thermotropic biaxial nematic phases seem to be rare, but biaxial smectic A phases less so. Here we use molecular field theory to study a simple two-parameter model, with one parameter promoting a biaxial phase and the second promoting smecticity. The theory combines the biaxial Maier-Saupe and McMillan models. We use alternatively the Sonnet-Virga-Durand (SVD) and geometric mean approximations (GMA) to characterize molecular biaxiality by a single parameter. For non-zero smecticity and biaxiality, the model always predicts a ground state biaxial smectic A phase. For a low degree of smectic order, the phase diagram is very rich, predicting uniaxial and biaxial nematic and smectic phases, with the addition of a variety of tricritical and tetracritical points. For higher degrees of smecticity, the region of stability of the biaxial nematic phase is restricted and eventually disappears, yielding to the biaxial smectic phase. Phase diagrams from the two alternative approximations for molecular biaxiality are similar, except inasmuch that SVD allows for a first-order isotropic-biaxial nematic transition, whereas GMA predicts a Landau point separating isotropic and biaxial nematic phases. We speculate that the rarity of thermotropic biaxial nematic phases is partly a consequence of the presence of stabler analogous smectic phases.

Calamitic and antinematic orientational order produced by the generalized Straley lattice model

We consider here a classical model, consisting of D_{2h}-symmetric particles in a three-dimensional simple-cubic lattice; the pair potential is isotropic in orientation space, and restricted to nearest neighbors. The simplest potential model is written in terms of the squares of the scalar products between unit vectors describing the three interacting arms of the molecules, as proposed in previous literature. Two predominant antinematic couplings of equal strength (+1) are perturbed by a comparatively weaker calamitic one, parameterized by a coupling constant -z ranging in [-1,0]. This choice rules out thermodynamically stable phases endowed with macroscopic biaxiality. The antinematic terms favor states with the corresponding molecular axes mutually orthogonal. Although the low-temperature phase of the special case with null calamitic term (PP0) is uniaxial and antinematically ordered, in the general case presented here both Monte Carlo and molecular-field approaches show that, for z close to zero, the models exhibit a low-temperature uniaxial nematic phase, followed by an antinematic one, and finally by the orientationally disordered one. On the other hand, for sufficiently large values of z, we only find evidence of uniaxial calamitic behavior, as expected by following the limiting cases.

Molecular-field-theory approach to the Landau theory of liquid crystals: uniaxial and biaxial nematics

Nematic liquid crystal phase diagrams in temperature-biaxiality space are usually complex. We construct a Landau theory based on the analogous molecular-field theory for orthorhombic biaxial nematic fluids. A formal procedure yields coefficients (some of which, unusually, can be tensorial) in this Landau expansion, correctly predicts the complete set of invariants formed from the ordering tensors, and avoids ad hoc parametrization of the molecular biaxiality. By regularizing the Landau expansion to avoid unwanted order parameter divergences at low temperatures, we predict phase behavior over the whole range of biaxiality. The resulting phase diagrams have the same topology as those of molecular-field theory.

Effect of Molecular Flexibility on the Nematic-to-Isotropic Phase Transition for Highly Biaxial Molecular Non-Symmetric Liquid Crystal Dimers

In this work, a study of the nematic (N)-isotropic (I) phase transition has been made in a series of odd non-symmetric liquid crystal dimers, the α-(4-cyanobiphenyl-4'-yloxy)-ω-(1-pyrenimine-benzylidene-4'-oxy) alkanes, by means of accurate calorimetric and dielectric measurements. These materials are potential candidates to present the elusive biaxial nematic (NB) phase, as they exhibit both molecular biaxiality and flexibility. According to the theory, the uniaxial nematic (NU)-isotropic (I) phase transition is first-order in nature, whereas the NB-I phase transition is second-order. Thus, a fine analysis of the critical behavior of the N-I phase transition would allow us to determine the presence or not of the biaxial nematic phase and understand how the molecular biaxiality and flexibility of these compounds influences the critical behavior of the N-I phase transition.

129Xe and 2H nuclear magnetic resonance (NMR) of xenon and deuterated-chloroform solutes in a thermotropic biaxial nematic liquid crystal

Properties of a thermotropic nematic liquid crystal (LC) TM35 (organosiloxane tetrapode) were investigated by means of nuclear magnetic resonance (NMR) spectroscopy of probe atoms and molecules. 129Xe isotope enriched xenon gas and 2H isotope enriched chloroform were used as probes. Both 129Xe and 2H NMR measurements were carried out over a wide temperature range, covering isotropic, nematic, and glass phases. A clear change in the behavior of the LC is observed within the nematic phase at 289 K by changing the temperature from low to high and inspecting the spectral line width. The derivative curve of the 129Xe chemical shift also reveals a change at 289 K. These changes may be an indication of hindered molecular rotation around the long axis, and the appearance of a biaxial nematic phase.

Molecular field theory for biaxial nematic liquid crystals composed of molecules with C(2h) point group symmetry

The biaxial nematic phase is generally taken, either explicitly or implicitly, to have D(2h) point group symmetry. However, it is possible for the biaxial phase to have a lower symmetry depending on that of its constituent molecules. Here we develop a molecular field theory for a nematogen composed of C(2h) molecules in terms of the nine independent second rank orientational order parameters defining the C(2h) biaxial nematic. In addition, there is a rank one order parameter constructed from two pseudovectors which is only nonzero in the C(2h) phase. The theory is simplified by removing all but the three dominant order parameters. The predicted phase behavior is found to be rich with three possible biaxial nematic phases and with the transitions involving a biaxial nematic phase exhibiting tricritical points.

Orientationally ordered phase produced by fully antinematic interactions: a simulation study

We consider here a classical model, consisting of D(2h) symmetric particles, whose centers of mass are associated with a three-dimensional simple-cubic lattice; the pair potential is isotropic in orientation space, and restricted to nearest neighbors. Two orthonormal triads define orientations of a pair of interacting particles; the simplest potential models proposed in the literature can be written as a linear combination involving the squares of the scalar products between corresponding unit vectors only, thus depending on three parameters, and making the interaction model rather versatile. A coupling constant with negative sign tends to keep the two interacting unit vectors parallel to each other, whereas a positive sign tends to keep them mutually orthogonal (antinematic coupling). We address here a special, extreme case of the above family, involving only antinematic couplings: more precisely, three antinematic terms whose coefficients are set to a common positive value (hence the name PPP model). The model under investigation produces a doubly degenerate pair ground state; the nearest-neighbor range of the interaction and the bipartite character of the lattice can propagate the pair ground state and increase the overall degeneracy, but without producing frustration. The model was investigated by a simplified molecular field treatment as well as by Monte Carlo simulation, whose results suggested a second-order transition to a low-temperature biaxially ordered phase; ground-state configurations producing orientational order have been selected by thermal fluctuations. The molecular field treatment also predicted a continuous transition, and was found to overestimate the transition temperature by a factor 2.

Ferroelectricity in low-symmetry biaxial nematic liquid crystals

Order parameters and phenomenological theory for both high- and low-symmetry biaxial nematic phases are presented and it is predicted that the chiral low-symmetry biaxial phase must be ferroelectric. This conclusion is based on general symmetry arguments and on the results of the Landau-de Gennes theory. The microscopic mechanism of the ferroelectric ordering in this chiral biaxial phase is illustrated using a simple molecular model based on dispersion interactions between biaxial molecules of low symmetry. Similar to the chiral smectic C* phase, the ferroelectricity in the chiral biaxial nematic phase is improper, i.e., polarization is not a primary order parameter and is not determined by dipolar interactions. Ferroelectric ordering in biaxial nematics may be found, in principle, in materials composed of chiral analogues of the tetrapod molecules which are known to exhibit biaxial phases.

Molecular model of biaxial ordering in nematic liquid crystals composed of flat molecules with four mesogenic groups

Relative stability of uniaxial and biaxial nematic phases is analyzed in a model nematic liquid crystal composed of flat molecules of C2h symmetry with four mesogenic groups rigidly linked to the same center. The generalized effective quadrupole mean-field potential is proposed and its constants are evaluated numerically for the pair intermolecular potential based on Gay-Berne interaction between mesogenic groups. The dependencies of the constants on molecular shape parameters are systematically analyzed. Order parameters of the uniaxial and biaxial nematic phases are evaluated by direct minimization of the free energy at different temperatures. The corresponding phase diagrams are obtained enabling one to study the effects of molecular model parameters on the stability regions of uniaxial and biaxial phases. The results are used to clarify the nature of experimentally observed biaxial ordering in nematic liquid crystals composed of tetrapode molecules with the same symmetry.

Mesogenic lattice models with partly antinematic interactions producing uniaxial nematic phases

The present paper considers nematogenic lattice models, involving particles of D_{2h} symmetry, whose centers of mass are associated with a three-dimensional simple cubic lattice; the pair potential is isotropic in orientation space and restricted to nearest neighbors. Let two orthonormal triads define orientations of a pair of interacting particles; the simplest potential models proposed in the literature can be reduced to a linear combination involving the squares of the scalar products between corresponding unit vectors only and depending on three parameters. By now, various sets of potential parameters have been proposed and studied in the literature, some of which capable of producing biaxial orientational order at sufficiently low temperature. On the other hand, in experimental terms, mesogenic biaxial molecules mostly produce uniaxial mesophases; thus we address here two very simple cases, involving a nematic (calamitic) term as well as one (model P0M) or two (model PPM) antinematic ones, whose coefficients are set equal in magnitude; when only one antinematic coefficient is used, the third one is set to zero. The calamitic term favors the alignment of two corresponding molecular axes, whereas antinematic terms or geometric constraints tend to keep two other pairs of axes mutually orthogonal. The models were investigated by molecular-field treatments and Monte Carlo simulation and found to predict a first- or second-order transitions between uniaxial nematic and isotropic phases; the molecular-field treatments yielded results in reasonable agreement with simulation.

Computer simulations of biaxial nematics

Biaxial nematic (N(b)) liquid crystals are a fascinating condensed matter phase that has baffled, for more than thirty years, scientists engaged in the challenge of demonstrating its actual existence, and which has only recently been experimentally found. During this period computer simulations of model N(b) have played an important role, both in providing the basic physical properties to be expected from these systems, and in giving clues about the molecular features essential for the thermodynamic stability of N(b) phases. However, simulation studies are expected to be even more crucial in the future for unravelling the structural features of biaxial mesogens at the molecular level, and for helping in the design and optimization of devices towards the technological deployment of N(b) materials. This review article gives an overview of the simulation work performed so far, and relying on the recent experimental findings, focuses on the still unanswered questions which will determine the future challenges in the field.

Biaxial and uniaxial phases produced by partly repulsive mesogenic models involving D2h molecular symmetries

The present paper considers biaxial nematogenic lattice models, involving particles of D2h symmetry, whose centers of mass are associated with a three-dimensional simple-cubic lattice. The pair potential is isotropic in orientation space and restricted to nearest neighbors. Let two orthonormal triads define orientations of a pair of interacting particles. The investigated potential models are quadratic with respect to the nine scalar products between the two sets of unit vectors. Actually, based on available geometric identities, these expressions can be reduced to diagonal form containing only the scalar products between corresponding unit vectors and depending on three parameters. Over the years, this comparatively simple functional form has also proven to be rather versatile. By now, various sets of potential parameters capable of producing mesogenic behavior of some kind have been proposed and studied in the literature. A new and simplified form was recently proposed and investigated by Sonnet, Virga, Durand, and De Matteis [A. M. Sonnet, E. G. Virga, and G. E. Durand, Phys. Rev. E 67, 061701 (2003); G. De Matteis and E. G. Virga, Phys. Rev. E 71, 061703 (2005)] and is known to support a biaxial phase at sufficiently low temperature. Following the idea of the above authors, we have studied a more extended range of parameters, including cases where biaxiality cannot be sustained in the pair ground state. In cases where a biaxial phase survives, an appropriate mean-field analysis may predict the existence of a direct second-order transition to the isotropic phase as well as a second-order sequence isotropic-to-uniaxial-to-biaxial. A second-order phase transition is also predicted, which involves isotropic and uniaxial phases only. Monte Carlo simulations have been carried out as well, for a few points in the parameter space, and found to produce results which partly confirm mean-field predictions.

Dominant biaxial quadrupolar contribution to the nematic potential of mean torque

Within the general quadrupolar model for biaxial nematic liquid crystals, whose potential of mean torque extends that in the Maier-Saupe theory with two extra interaction terms, we propose a quantitative criterion to identify the dominant biaxial interaction. We show that the ratio of the biaxial-to-uniaxial and uniaxial-to-isotropic transition temperatures is almost independent of one interaction parameter, thus indicating the other as dominant. We also show that there is a significant mismatch between the principal orientational order parameters predicted by the theory and those measured for the biaxial phase of a tetrapode.

Landau-de Gennes theory of biaxial nematics reexamined

Recent experiments report that the long-looked-for thermotropic biaxial nematic phase has been finally detected in some thermotropic liquid crystalline systems. Inspired by these experimental observations, we concentrate on some elementary theoretical issues concerned with the classical sixth-order Landau-de Gennes free energy expansion in terms of the symmetric and traceless tensor order parameter Q alpha beta. In particular, we fully explore the stability of the biaxial nematic phase giving analytical solutions for all distinct classes of the phase diagrams that theory allows. This includes diagrams with triple, critical, and tricritical points and with multiple (reentrant) biaxial and uniaxial phase transitions. A brief comparison with predictions of existing molecular theories is also given.

High-resolution calorimetric study of a liquid crystalline organo-siloxane tetrapode with a biaxial nematic phase

High-resolution adiabatic scanning calorimetry and differential scanning calorimetry have been employed to study the thermal behavior of an organo-siloxane tetrapode reported to exhibit a biaxial nematic phase. No signature of the uniaxial to biaxial nematic phase transition could be retraced in sequential heating and cooling runs under different scanning rates, within the experimental resolution. The results obtained reveal that an extremely small heat should be involved in the uniaxial to biaxial nematic phase transition. The isotropic to uniaxial nematic transition at 318+/-0.01 K is very stable, and it is weakly first order with a rather small latent heat of 0.20+/-0.02 J/g .

Uniaxial rebound at the nematic biaxial transition

Over the last few years, renewed interest has been raised by the simplified general interaction models proposed by Straley for mesogenic molecules possessing the D{2h} symmetry and capable of producing biaxial nematic order. It has already been shown that, in the presence of certain special symmetries, just two out of the four order parameters that are in general necessary, suffice for the description of a biaxial phase. For some other range of parameters, these reducing symmetries do not hold, and, moreover, a mean-field treatment has to be suitably changed into a minimax strategy, still producing a transition to a low-temperature biaxial phase. Upon studying the general parameter range, we identify as a common feature the behavior of a uniaxial order parameter, attaining a local minimum at the biaxial-to-uniaxial transition temperature, and recognizably increasing away from it. This finding is confirmed by a Monte Carlo simulation.

Quadrupolar projection of excluded-volume interactions in biaxial nematic liquid crystals

We compute the quadrupolar approximation to the excluded-volume interaction between hard spherocuboids, which applies to both platelets and spheroplatelets as special cases. We show that this approximation can be written as the superposition of two London interactions: one attractive and the other repulsive. This conclusion also proves why the phase diagram for the excluded-volume interaction of spherocuboids is expected to feature a direct isotropic-to-biaxial transition at a single Landau point.

Universal mean-field phase diagram for biaxial nematics obtained from a minimax principle

We study a class of quadratic Hamiltonians which describe both fully attractive and partly repulsive molecular interactions, characteristic of biaxial liquid crystal molecules. To treat the partly repulsive interactions we establish a minimax principle for the associated mean-field free energy. We show that the phase diagram described by Sonnet [Phys. Rev. E 67, 061701 (2003)] is universal. Our predictions are in good agreement with the recent observations on both V-shaped and tetrapodal molecules.

Bifurcation analysis and computer simulation of biaxial liquid crystals

We extend the analysis of a mean-field model for biaxial liquid crystals recently proposed by Sonnet et al. [Phys. Rev. E 67, 061701 (2003)]. In particular, we perform a bifurcation analysis of the equilibrium equations and derive the complete phase diagram. We show that two order parameters suffice to label all equilibrium phases, though they exhibit different bifurcation patterns. A Monte Carlo simulation study is performed as well, confirming qualitatively the predictions of this analysis.