The Influence of Disorder on Thermotropic Nematic Liquid Crystals Phase Behavior

We review the theoretical research on the influence of disorder on structure and phase behavior of condensed matter system exhibiting continuous symmetry breaking focusing on liquid crystal phase transitions. We discuss the main properties of liquid crystals as adequate systems in which several open questions with respect to the impact of disorder on universal phase and structural behavior could be explored. Main advantages of liquid crystalline materials and different experimental realizations of random field-type disorder imposed on liquid crystal phases are described.


Introduction
Understanding of phase and structural behavior of various systems which are randomly perturbed by some static origin of disorder is of considerable interest for different branches of physics [1,2,3]. Presence of disorder is practically unavoidable in any system (e.g., it is enforced via "impurities", thermal fluctuations...) and could significantly modify its properties quantitatively or even qualitatively. Randomly perturbed configurations often display several universal features. Such systems are therefore related from the mathematical point of view, although they could be seemingly completely different physically. Of particular interest are phases and structures, which are in the absence of disorder reached via a continuous symmetry breaking (CSB) transition [4] and exhibit long range order (LRO). Such configurations are extremely susceptible to disorder due to the existence of a Goldstone mode in the broken symmetry phase. The mode's amplitude tend to diverge in the long wavelength limit because the energy costs of such excitations are negligible.
The pioneering studies have been mostly carried out in various randomly perturbed magnetic materials [1]. In order to understand and to control mean observed features minimal models have been developed in which magnetic spins interact through the Heisenberg or Ising form of the coupling interaction. In the absence of disorder, such systems exhibit continuous symmetry breaking when crossing from the paramagnetic to a magnetically ordered phase. The quenched disorder was either introduced via random field (RF) or random anisotropy magnetic (RAM) field-type coupling term exhibiting one fold and two fold axis symmetry [5,6], respectively. Let us list some results of universal validity (i.e., the systems of interest exibit CSB) which such models yield.
It was shown that the 1st order phase transitions become gradual as disorder exceeds the critical value [7]. With respect to the 2nd order phase transitions we refer to the Harris criterion [8]. It claims that the random bond disorder can change critical properties of a transition if α > 0. Here α stands for the critical coefficient describing thermal behavior of the specific heat. On the other hand, the renormalization study of the random anisotropy magnetic system [5] showed that the fixed point corresponding to the nonrandom critical behavior is unstable even with respect to infinitesimally weak disorder strength. This instability indicates either smoothing of the transition or a transformation into the 1st order transition. The latter effect could arise due to the static disorder affecting the Goldstone (also called gauge) type continuum field [9] in a reminiscent way as the thermal fluctuations of this field trigger the so called Halperin-Lubensky-Ma effect [10,11] (i.e., they transform the 2nd order phase transition into the 1st order). Note that the static disorder is much more efficient in comparison with the thermal one because of its persistent influence. It was shown that the static random field fluctuations could behave like thermal fluctuations with the lower marginal dimensionality increased by two [12]. Regarding the structure of the disordered magnetic phase, some studies predict that the symmetry broken phase exhibits a domain type structure characterized [12] by a single characteristic length scale ξ d . In that case the disordered phase exhibits a short range order (SRO). One of the cornerstone of statistical mechanics of disorder, the so called Imry-Ma theorem [12], claims that below the spatial dimension d = 4 the ordered state is unstable to an arbitrary weak random field that couples linearly with the order parameter. The latter exhibits continuous symmetry breaking transition on entering the lower symmetry phase. Furthermore, the scaling law ξ d ∝ w 2/(d−4) is predicted, where w measures the disorder strength. However, some studies using different approaches predict for random anisotropy models and a weak enough disorder strength the algebraic decay of correlations [13,14], the hallmark of quasi long range order (QLRO). But latter studies reveal that even in this case a characteristic length is present marking the distance above which the correlations decay rather weakly [15].
In recent years it has been shown that randomly perturbed liquid crystal (LC) phases and structures [16] present an adequate system in which several universal mechanism in the realm of statistical mechanics of disorder could be studied. Liquid crystal (LC) phases [17,18] are often referred to as a fourth state of matter. In the literature, however, one can find also plasma or quark-gluon plasma as candidates for the same ranking. But already such a candidature indicates the important role this state of matter plays in nature. Most LC phases represent an intermediate state between ordinary liquids and crystals. In most cases they flow like a liquid. On the other hand they possess some long-range order (orientational or translational, or both) usually characterizing crystalline solids. They are also typical representatives of soft-matter systems [17,18], in which a relatively small amount of locally supplied energy can cause a response on a macroscopic scale.
LCs have several extraordinary (mechanical, electrical and optical) properties which make them indispensable in biological systems and in our daily life [18]. But aside these facts LC phases offer a convenient tool to study or test several fundamental laws of physics [4,19,20]. For example, studies on the coarsening dynamics following a sudden phase transition from ordinary liquid to orientationally order nematic phase can shed light on phenomena in the early universe [21]. The reasons behind this are as follows. There exists a variety of LC phases, where every phase can exhibit various structures [17,22] depending on boundary conditions and/or direction and strength of external electric or magnetic fields. The resulting states display extremely rich diversity of physical phenomena. Furthermore, LCs are adequate for experimental observations [17,22]. The samples can be prepared relatively easily because LCs have liquid properties and can be shaped by boundaries. In addition their structure can be strongly influenced by conditions at the confining surface [23] or by external electric or magnetic fields. The characteristic lengths and time responses to locally induced perturbations are accessible to conventional experiments. In particular, several optical measurements can be used because LCs are transparent and show an optical anisotropy strongly linked to electric and magnetic anisotropy.
In order to study influence of random field-type disorder on phase and structural behavior one commonly enforces disorder to LCs either by mixing them with "impurities" or confining LCs to various porous matrices [16]. As impurities, aerosil nanoparticles [3,24,25,26] are typically used, which can form qualitatively different random networks as the concentration of aerosils is varied. In terms of porous confinements, aerogels [27,28], Russian glasses [29], Vycor glasses [30], and Controlled Pore Glasses (CPG) [31,32] are mostly used. Therefore, by choosing different LCs phases and perturbing porous matrices or "impurities" one can impose different types and strength of disorder.
The paper is organized as follows. In Section II, we present a brief overview on liquid crystals, containing the principal characteristics and properties of the three main classes and the definition of the nematic orientational order parameter. In Section III, the basic properties of principal models of nematicisotropic phase transition are shown. In Section IV, the influence of the random field on structural and phase transition properties is theoretically described in the framework of a mesoscopic phenomenological model. The experimental observations which support the theoretical model are also described. In the final Section V, we draw some conclusions.

Liquid Crystals: A Brief Overview
Liquid crystals were discovered in 1861 by Planer, confirmed in 1888 by Reinitzer [33], named by Otto Lehmann in 1889 [34] and understood in terms of the local ordering of elongated particles and classified by Georges Friedel in 1922 [35].
Liquid crystals or mesophases constitute a state of matter whose physical properties are intermediate between those of an isotropic liquid and a crystalline solid [36,37,38,39,40,41,42,43,44,45]. The existing long-range order in a crystal consisting of anisotropic molecules is both positional (the molecules occupy specific sites in a lattice) and orientational (the molecular axes point in specific directions). When such a crystal is heated from the solid phase it is possible that both types of order (positional and orientational) disappear at the same temperature and the resulting phase is an isotropic liquid. Alternatively, it is possible that one type of ordering survives until a higher transition temperature. If the orientational order disappears leaving the positional order intact the corresponding phase is called plastic crystal. On the contrary, if the positional order either fully or partially disappears while some degree of orientational order is maintained, the resulted phase is called liquid crystal.
Liquid crystals possess many of the mechanical properties of a liquid, e.g., high fluidity and the inability to support shear, but on the other hand, they have some properties similar to crystals, e.g., they are birefringent and have anisotropic magnetic and electric susceptibilities. The relative large value of latent heat (∝ 250 J/g) at crystal-liquid crystal phase transition indicates that most of the order of a crystal is lost when it transform to a mesophase. On the contrary, the latent heat at liquid crystal-isotropic phase transition is much smaller , typically about 5 J/g.
The constituent organic molecules of liquid crystalline materials could be in general very different, but it is essential that they are anisotropic and rigid for at least some portion of the molecule length (since it must maintain an elongated shape in order to produce interactions that favor alignment). The liquid crystals composed from rod-like molecules are called "calamitics", while those formed by disk-shape molecules are known as "discotic" [41,42].
Two types of liquid crystals must be differentiated: (i) thermotropic and (ii) lyotropic. The transitions involving thermotropic liquid crystals are effected by changing temperature. Lyotropic liquid crystals are found in solutions and concentration is the important controllable parameter. Usually a lyotropic liquid crystal molecule combines a hydrophobic group at one end with a hydrophilic group at the other end. These amphiphilic molecules forms ordered structures (micelles, vesicle, lamellar phases or reversed phases) in both polar and non-polar solvents as in the case of soaps and various phospholipids [45].
Another important class of liquid crystals is derived from certain macromolecules, usually in solution but sometimes even in the pure state. These liquid crystals are known as "liquid crystal polymers" [43,46,47,48]. They can be divided corresponding either to molecular anisotropy or location of the mesogenic group. On the basis of molecular anisotropy two types of structures are possible. In the first type the mesomorphism is due to rigid macromolecule as a whole or from several units made up of rigid chain segments, whereas in the second class only the monomer unit possesses a mesogenic structure connected to the rest of molecule via a flexible spacer. On the basis of location of the mesogenic unit the liquid crystal polymers can be divided into: main-chain liquid crystal polymers and side-chain liquid crystal polymers [44].
Based upon their symmetry, Friedel distinguishes three major classes of liquid crystals-nematics, cholesterics, and smectics.

Nematic, Cholesteric and Smectic Liquid Crystals
An isotropic liquid possesses full translational and orientational symmetry T (3) × O (3). At the transition in the nematic phase the translational symmetry T (3) remains the same, but the rotational symmetry O(3) is broken. In the simplest uniaxial nematic structure the group O(3) is replaced by one of the uniaxial symmetry groups D ∞ or D ∞ h . The nematic molecular order is shown schematically in Figure 1. In this phase there exists no long-range positional order (between the centers of mass of the molecules). The molecules tend to be parallel to some common axis, labeled by a unit vector n, called "director". The direction of n is arbitrary in space; in practice it is imposed by external forces such as the guiding effect of the wall of the container or electric and magnetic fields. The states of director n and − n are indistinguishable even if the molecules may be polar (no ferroelectricity has been observed). Microscopically, an equal number of molecules points "up" and "down". The uniaxial nematic phase is characterized by rotational symmetry of the system around the director n. One important consequence of this is that there are two different principal components of the second rank tensorial properties as, e.g., magnetic susceptibility, dielectric constant, or refractive indices. In general two uniaxial nematic phases can be distinguished: (i) prolate nematic phase which occur for rodlike molecules and (ii) oblate nematic phase characteristic for disklike molecules.
The biaxial nematic phase results due to further breaking of the rotational symmetry and is characterized by three orthogonal directors, the Goldstone modes n, l, and m = n × l. Due to the lack of polarity of the known biaxial nematics, the directions n = − n, m = − m, and l = − l are equivalent. The biaxial nematic phase possess a D 2 h point group symmetry and the corresponding second rank tensor properties have three different principal components.
If the liquid crystal molecules are chiral (lacking inversion symmetry), the uniaxial nematic phase is replaced by the chiral nematic phase or twisted nematic phase or cholesteric phase. Locally, a cholesteric is very similar to a nematic. Again, the centers of mass have no long-range order and the molecular orientation shows a preferred axis labeled by a director n. However, n is not constant in space. The preferred conformation is helical (see Figure 2). The helical distortion can be found either in pure cholesterol esters (which are chiral) or by dissolving in a nematic a fluid composed by chiral molecules. Figure 2. The arrangement of molecules in the cholesteric phase; the plans have been drawn for convenience, but do not have any specific physical meaning.

t w i s t a x i s
When the crystalline order is lost in two dimensions, one obtains the two-dimensional liquids, called smectics. Smectics liquid crystals have layered structures with a well-defined interlayer spacing that can be measured by X-ray diffraction. In addition to the orientational order, the smectic molecules exhibit some correlations in their position. In most smectics the molecules are mobile in two directions and can rotate around one axis. The interlayer attractions are weak compared with the lateral forces between molecules and the layers are able to slide over one another relatively easily which gives rise to fluid property of the system with higher viscosity than nematics. When smectic and nematic phases are found in one compound, the nematic phase is almost always found at higher temperatures with the exception of reentrant nematic phase.
The observed smectic phases differ from each other in the way of layer formation and the existing order inside the layer [42,44,45]. The simplest is the smectic A phase (the phase recognized by G. Friedel (see Figure 3) with symmetry T (2) × D ∞ h . In the phase smectic A, the average molecular axis is normal to the smectics layer. The layer thickness can vary from a value close to the length of the molecule in thermotropic liquid crystals to the periodicity as large as several µm in lyotropics. Inside each layer, the centers of mass of molecules show no long-range order and the molecules have considerably freedom of rotation around their long axis; each layer is a two-dimensional fluid and the whole structure may be defined as an orientational ordered fluid on which is superimposed a one-dimensional density wave. The distortions due to flexibility of layers leads to optical patterns known as focal conic textures [45]. When temperature is decreased, the smectic A phase may transform into a phase possessing a lower symmetry. The breaking of D ∞ h may lead to the appearance of tilting of molecules relative to the smectic layers. This phase is called smectic C which possesses T (2) × C 2 h symmetry.
As an example of compound that can be in both nematic and smectic A phases we cite the very much used octyl-4-cyanobiphenyl (8CB) with chemical formula shown in Figure 4, which is smectic A between 24 • C and 34 • C, and nematic between 34 • C and 42.6 • C.

Order Parameters
The fundamental characteristics of a liquid crystal is the presence of long-range orientational order while the positional order is either absent (nematic phase) or limited (smectic phases). One phase differs from another by its symmetry. The transition between phases of different symmetry can be described in terms of an order parameter Q that must satisfy the following requirements: (i) Q = 0 in the more symmetric (less ordered) phase and (ii) Q = 0 in the less symmetric (more ordered) phase. These requirements do not define the order parameter in a unique way, but in spite of this arbitrariness, in many cases the order parameter can be defined in a natural way. In the case of liquid-vapor transition, the order parameter is a scalar and is defined as the difference in density between the two phases. In the case of ferromagnetic transitions, the order parameter is the magnetization, a vector with three components. In general the order parameter can be defined as follows. Identify the symmetry of a phase. Next take the basis functions of the irreducible representations of this group symmetry. Expand the density distribution in the basis functions. The coefficients of this expansion are the order parameters. The order parameter as a quantitative measure of order in the system is used in such a diverse phase transitions as solid-solid, liquid-solid, in liquid crystals, superconductivity, superfluid helium, ferromagnetic-paramagnetic and polymeric systems.

Nematic order parameter
The macroscopic orientational nematic order parameter is a traceless, symmetric, second rank tensor with components given by [42,49] where m = n × l. The eigenvalue of Q largest in absolute magnitude gives the degree of orientation in the preferred direction, the corresponding eigenvector identifying the preferred direction n (the director in the uniaxial nematic phase). The difference between the remaining two eigenvalues characterizes the degree of biaxiality with the biaxiality axis l specified by the eigenvector corresponding to the second largest eigenvalue. The general biaxial state is thus specified by five parameters, which may be considered either as (i) the five independent elements of Q, or as (ii) the strength of uniaxial ordering S, the biaxial scalar amplitude η, two angles to specify the orientation of the director n, and a third angle to specify the orientation of l about n. Because Q is symmetric, it is always possible to find a reference frame in which it is diagonal, with its eigenvalues as diagonal elements: A uniaxial nematic state is described by the condition that two eigenvalues of Q coincide, while a biaxial nematic state is characterized by the condition that all eigenvalues of Q be distinct.
The microscopic ordered parameters S and η are related to specific molecular models and are obtained as expansion coefficients of the singlet orientational distribution function in a suitable base set [42,44,50].
The uniaxial nematic state state is specified by only three parameters: the magnitude S and two angles to orient the director. To define the microscopic orientational order parameters of the nematic, we start with the singlet orientational distribution function, defined as which is normalized to unity U (X 1 , X 2 , ...X N ) is the potential energy of N particles and β = 1/k B T with k B the Boltzmann constant. Z N is the configurational partition function of the system Here the X i = ( r i , Ω i ) specify the position r i of the center of mass of the ith molecule and its orientation Ω i described by the two angles (θ i , φ i ).
If the mesogenic molecules possess cylindrical symmetry, f (Ω) depends only on the angle θ between the director and the molecular symmetry axis and the orientational order parameters are defined by where P l are the even Legendre polynomials. From knowledge of f (θ) all the orientational order parameters P 2 , P 4 , etc., can be calculated. The most important one is usually called the scalar order parameter S defined as S = P 2 = 1 2 (3 cos 2 θ − 1). In the isotropic liquid where all orientations are equally probable, S = 0. In a fully oriented nematic phase where all the symmetry axes of the molecules are parallel with the director S = 1, and in fully oriented nematic phase where the symmetry axes of the molecules are distributed in the plane orthogonal to the director, S = −1/2.

Nematic-Isotropic Phase Transition: Theoretical Models
The experimental observations [41,43] using various techniques [38,39] show that the scalar order parameter S decreases monotonically as the temperature is raised in the nematic range and drops abruptly to zero at the nematic-isotropic phase transition temperature (T N I ). In the case of uniaxial nematic the jump of S at nematic-isotropic phase transition has a value in the range of 0.25 − 0.5 depending on the mesogenic material and usually the latent heat is only 1 kJ/mol. Thus, the nematic-isotropic phase transition is weakly (small latent heat) first order (the jump of the order parameter). This leads to large pretransitional [39,41] abrupt increases of response functions as specific heat, thermal expansion and isothermal compressibility near T N I . The changes of entropy and density associated with this transition are typically only a few percent of the corresponding values for solid-nematic phase transition. The theories of nematic phase and of nematic-isotropic phase transition have been developed in two directions: (i) microscopic and (ii) mesoscopic.

Microscopic Models
In the microscopic models we deal with the system at the molecular level and we start the calculations on the basis of partition function or the exact density functionals of the free energy. Accordingly, for a complete and satisfactory microscopic theory a knowledge of the intermolecular interaction is necessary. Unfortunately such a knowledge is almost entirely lacking. Even if the essentials of the intermolecular interaction are known, the successful application of the theory is very difficult due to enormous calculation problems. Owing to these difficulties the model potentials have been introduced, i.e., the most relevant characteristics of the molecules and their interactions are represented in term of simple models. The purpose of the remaining of this section is not the detailed description of all the technical details of the microscopic models, but rather the presentation of their basic properties.

The Onsager Theory
The first extremely successful molecular model (of mean field type) of liquid crystal ordering was developed by Onsager [51] in order to describe the ordering in the Tobacco Mosaic Virus (TMV) solution. TMV are rigid cylindrically molecules of length L = 300 nm and diameter D = 18 nm, each composed of 2130 proteins arranged in a helical way and formed a hollow cylinder. The diameter of the cylinder hole in the cylinder is about 2 nm. TMV in water solution forms a birefringent liquid crystal phase when its concentration exceeds 2%. In the model the molecules are rigid rods of length L and diameter D, which interact only via the hard-core repulsive forces. The rigid rod is usually modeled by a spherocylinder, which is a straight circular cylinder capped on each end by a hemisphere of the same radius. Since in the system the attractive interactions are absent, the temperature sets only the energy scale. The calculation of Onsager is based on a cluster expansion for the free energy as a functional of the singlet orientational distribution function and truncating the series after the second virial coefficient [44,51,52].
The configurational free energy (per rod) of a system of hard rods can be expressed as where F 0 is an additive constant and c is the concentration of rods. The second term in Equation (7) describes the drop in entropy associated with molecular alignment, while the third term describes the excluded volume effect; is the mutual excluded volume of two rods with orientations Ω i and Ωj.
The singlet orientational distribution function can be determined by minimizing the free energy subject to the constraint f (Ω)dΩ = 1. The result is the self-consistent equation where λ is the Lagrange multiplier determined by the normalized condition.
Onsager obtained an approximate variational solution which is based on a trial function of the form f = (const) cosh(α cos θ), where α is a variational parameter and θ is the angle between the molecular axis and the nematic axis. It was found [51] that the system exhibits an abrupt first-order phase transition from the isotropic to the nematic phase characterized by the following values of volume fractions of the isotropic φ iso , nematic phase φ nem , and order parameter S at coexistence The Onsager model of nematic-isotropic phase transition is exact in the limit of the infinite length to width ratio L/D and is was argued that the truncation of the cluster expansion series after the second virial coefficient can be quantitatively justified only for very long rods L/D > 100, which are typical of polymeric system. For the shorter rods, the Onsager theory tends to become inaccurate at the high densities typical of nematic-isotropic transition, and more refined approaches have been developed [44].

The Maier-Saupe Model
The mean-field model of Maier and Saupe [42,44] assumed that the nematic orientational order is caused by the anisotropic part of the dispersion interaction between molecules. Now, it is well accepted that the exact nature of the interaction need not be specified for the development of the theory. To obtain the results of the Maier Saupe theory, only an anisotropic potential with a particular dependence on the molecular orientations is required [53].
The standard Maier-Saupe free energy difference between nematic and isotropic phases can be written as where S(S) is the decrease in entropy due to the alignment of the molecules and the last term represents the anisotropic contribution to the internal energy. We now recall briefly the mean field calculation of S(T ).
The trial singlet orientational distribution function is chosen to be given by with where m is a dimensionless measure of the mean field aligning a molecule and P 2 (cos θ) = (3 cos 2 θ − 1)/2 is the second order Legendre polynomial. The conjugated nematic order parameter S(m) is related to Z by The entropy change associated with the nematic orientational order is given by The order parameter is calculated by minimizing ∆F with respect to S which leads to where S(m) is calculated from Equation (13). The results lead to a first order nematic-isotropic (NI) phase transition at and The possible extensions of the Maier-Saupe theory have been presented in extenso in [44].

The van der Waals (vdW) and density functional theories(DFT)
A complete molecular theory of the mesophases must include both anisotropic short range repulsive and long-range attractive forces. The vdW types theories [44,54,55] recognize that the predominant factor in determining the liquid crystalline stability are the short range repulsive forces and that the role of the attractive long range forces is, to a first approximation, merely to provide a negative, spatially uniform mean field in which the molecules move. All the vdW type theories can be derived from a perturbation method, within the mean field approximation, considering only the first order approximation term. Despite its quantitative inadequacies, the vdW approach indicates that the anisotropy of the short range intermolecular repulsions plays a major role in determining nematic orientational order and stability and cannot be neglected, even to a first approximation.
The classical version and development in the present form have been discussed at length in other reviews [44,56,57,58]. The main idea of the theory is that in the presence of a potential, the density is in fact uniquely determined. The net result is that the grand thermodynamic potential can be written as a functional of the density. Minimizing the grand thermodynamic potential, the thermodynamic properties of the system are obtained.

The Phenomenological Landau-de Gennes Theory
The Landau-de Gennes theory of phase transitions in liquid crystals [36,37,42,59] is a phenomenological mesoscopic theory which apply the general rules of Landau theory [60] originally proposed by Landau for second order phase transitions in crystals. The basic idea is that sufficiently closed to the order-disorder phase transition the thermodynamic potential can be expanded in a series of powers of the order parameter. The motivation is derived from the continuity of the order parameter at the second order phase transition. Hence the Landau original procedure is, in principle, restricted to second order phase transition. In this section the application of Landau-de Gennes theory to nematic-isotropic phase transition will be summarized.
The most general form of the free energy F of the system is given by where V is the nematic volume and A is the area of the interface confining it. Near the nematic-isotropic transition point, in a mesoscopic approach the phenomenological Landaude Gennes form of the homogeneous free energy density f h term is given by where f (T, 0) is the free energy density of the isotropic phase (usually the free energy is normalized such as f (T, 0) = 0), T * is the undercooling limit temperature of the isotropic phase and the coefficients a, B, and C depend only on the substance (for 8CB they have the following values: a = 3.3 · 10 4 J/K cm 3 , B = 8.9 · 10 5 J/cm 3 , and C = 5.6 · 10 5 J/cm 3 [45]). They do not have a physical significance and their values are obtained by fitting the theoretical values of some quantities (i. e., the values of the order parameter at the transition, of T N I − T * , and of latent heat, respectively) with experimental data. The main properties of the expansion (19) can be summarized as follows: (i) There is no term linear in Q. This allows for the possibility of an isotropic phase (Q = 0). In the case of external fields, a linear term has to be included and the isotropic phase transforms into a paranematic phase (a phase with a very low degree of orientational order).
(ii) If B > 0, the transition is first order (the order parameter changes discontinuous at the transition), while if B = 0 the transition is second order (the order parameter is continuous at the transition, but its first derivative with respect to temperature is discontinuous). The other possible mechanism to mimic the first order phase transition is to consider B = 0 and C < 0. In this case a stabilizing six order term with E > 0 is required.
(iii) To ensure stability of the nematic phase, C must be positive.
Using the representation (2) for the order parameter, the corresponding homogeneous Landau-de Gennes free energy density (19) is given by: The uniaxial nematic state (with the choice of n as the director, i. e., the vector with the maximum eigenvalue of Q) corresponds to η = 0, value obtained minimizing f (T, S, η) with respect to η. The corresponding order parameter is defined as (see Equation (1)) while the corresponding homogeneous free energy density is given by which describes a weakly first order uniaxial nematic-isotropic phase transition. At nematic-isotropic phase transition temperature T = T N I = T * + B 2 /24aC the two phases, uniaxial nematic (S nem = B/6C) and isotropic (S iso = 0) coexist in equilibrium. The minimum corresponding to the nematic phase becomes metastable for T N I < T < T + = T * + 3B 2 /64aC and does not exist for T > T + , the superheated limit temperature of the nematic phase. The minimum corresponding to the isotropic phase is metastable for T * < T < T N I and does not occur when T < T * , the supercooling limit of the isotropic phase. The latent heat at the transition is given by The elastic free energy density (f e ) opposes the spatial variation in Q and is expressed as [42]: where L is a representative bare nematic elastic constant. In the Frank description of the nematic ordering, this form of f e corresponds to the approximation of equal Frank elastic constants. The most essential surface contribution (anchoring term) of the free energy is given by where W is the anchoring energy. The anchoring term tends to align liquid crystal molecules along the unit vector e s , referred to as the easy axis. For the uniaxial ordering, using the representation (21) for the order parameter, the homogeneous free energy density term is given by Equation (22), while the other two terms can be expressed as

Relation between Landau-de Gennes and Maier-Saupe Free Energies for Uniaxial Nematics
To obtain an insight into the structure and origin of the phenomenological parameters a, B, and C of the homogeneous Landau-de Gennes free energy, we start from the mean field Maier-Saupe approach on the molecular level. For uniaxial nematic composed of cylindrically symmetric molecules, the states with orientational order are distinguished by an axis of symmetry, the director n and an infinite set of scalar order parametersP where n > 0 is an integer, P 2n (cos θ) are the Legendre polynomials and f (Ω) is the singlet orientational distribution function defined in Equation (3). Among these order parameters,P 2 is by far the most important both from theoretical and experimental point of view, and is referred to as the nematic scalar order parameter S. In the isotropic phase all orientations are equivalent and f (Ω) = 1/4π; S = 0.
Expanding the singlet orientational distribution function in terms of Legendre polynomials and keeping only the leading term, we obtain Using Equation (14), Equation (10) can now be written as a power series in S [61,62] ∆f where T * = 125U/k B . In the spirit of Landau-type approach, the dependence on temperature in cubic and quadratic terms in Equation (29) can be neglected and the free energy density of the nematic uniaxial phase takes the well-known Landau-de Gennes expression (22). It is to be noted that the T * term in Landau-de Gennes expansion is the only term that originates from the interaction Maier-Saupe term, the other terms originate from the Maier-Saupe entropy.

Influence of Random Field on Structural and Phase Transition Properties: A Mesoscopic Phenomenological Model
Liquid crystals, by virtue of their fluidity, intrinsically soft elasticity and experimental accessibility, offer opportunities for the study of the structural and dynamical effects of quenched disorder, which can be introduced, for example, by confinement within appropriate random porous media [16]. Such studies are also of interest in connection with composite electro-optic materials composed of mixtures of liquid crystal and impurities, or colloids, or nanoparticles. In such mixtures, a degree of randomness is imposed on the liquid crystal by the second component. The changes of the physical behavior of liquid crystal in such systems are generated by three main causes: (i) surface effects, (ii) confinement, and (iii) randomness. In this section we focus only on the theoretical study of the influence of randomness on phase transitions properties in liquid crystals.
Numerous studies [3,9,63,64,65,66,67,68,69,70,71,72,73] have suggested that similar random field type models may provide a useful context within which to understand the influence of randomness on physical properties of liquid crystals. These models assume that the porous medium (or the second component in a mixture) imposes a local field which has a randomizing effect on the nematic director. Such a randomizing field depending neither on temperature nor on the phase structure of liquid crystal, is generally refereed to as quenched disorder. The new equilibrium structure of the liquid crystal is the consequence of the competition between the quenched disorder provided by the randomness and the elastic force which tends to minimize the distortions of the nematic director.
In the following we analyze structural and phase behavior of a nematic liquid crystal experiencing quenched random anisotropy field (RAN) disorder in the simplest possible mean-field type approach. We first consider a typical domain structure in the nematic director field. We show that a domain structure is practically always expected, at least temporarily, if a lower symmetry (i.e., the nematic) phase is reached via a continuous symmetry breaking phase transition. If disorder is present it could stabilize domains, the characteristic size of which is determined by the disorder strength. We show that such systems might display pronounced memory effects. Then we discuss phase behavior across the isotropic-nematic phase transition assuming the domain-type nematic structure. We consider case of (i) nematic LC, and (ii) nematic-non nematogenic mixture experiencing RAN. Finally we present experimental systems in which described behavior is observed.

Structural Behavior : Domain-Type Patterns
The main part of the paper is concerned with phase and structural behavior of randomly perturbed systems in which domain-type patterns are assumed. For this purpose we list evidences supporting this description in phases or structures exhibiting continuous symmetry breaking, where the lower symmetry phase possesses long range order (or at least quasi long range order).
In general, at the mesoscopic continuum level the degree of ordering in the broken symmetry configuration is determined by an order parameter field (OPF). On the other hand a selected symmetry breaking state is described by a symmetry breaking field (SBF), also referred to as a gauge field [17]. In absence of external ordering fields OPF and SBF in general respond to local perturbations on apparently different length scales (e.g., these sizes can be comparable in severely confined structures, where the geometrical size characterizing the confinement is comparable to the relevant order parameter correlation length). OPFs relax towards an equilibrium value on distances comparable to order parameter correlation length ξ, which reflects intrinsic material properties [17]. In a thermotropic LC phase values of ξ are typically comparable to several nm far bellow the phase transition temperature [17]. On the other hand SBF is susceptible to a large scale geometry of the sample [17,74]. For example, if a frustration is imposed on SBF via conflicting boundary conditions at surfaces separated by a distance R, then the field exhibits gradual changes over all the available space. This is due to the to the existence of Goldstone fluctuation modes in SBFs.
An example of a continuous symmetry breaking in orientational ordering represents the isotropic to nematic liquid crystal phase transition. In the nematic phase ordering can be well described by the nematic director field − → n (SBF) and the uniaxial nematic orientational order parameter S (OPF). An example of a continuous translational symmetry breaking represents the nematic-smectic phase transition. In the smectic A phase the ordering is typically described by the phase factor φ (SBF) and the translational order parameter η (OPF).
For illustration purposes we consider orientational ordering in the nematic phase in the remaining part of this Section. We first show that a domain-type pattern in the SBF is inevitably formed at least temporally in the lower symmetry configuration after a quench from the higher symmetry phase even in absence of disorder. Then we give evidences that such domain pattern could be stabilized by impurities, which impose a quenched random field-type of disorder. In both cases the domain patterns are well described by a single characteristic length. We determine regimes where such patterns are expected.

Domain Coarsening Following a Fast Enough Phase Transition
We consider the onset of orientational ordering in the thermotropic nematic phase which is reached via the temperature quench from the isotropic phase in absence of static disorder. Immediately after the quench in causally disconnected parts of the system, a random value of SBF is chosen depending on a local preference mediated by a fluctuation [75,76]. The resulting nematic structure is strongly orientationally frustrated and consequently spatially dependent variations in degree of ordering appear. The established configuration is extremely energetically costly and consequently it undergoes coarsening dynamics towards the uniform equilibrium nematic configuration. In this process the dominant role is played by topological defects [19]. These refer to regions in space where SBF is not uniquely defined. Defects can be treated as objects carrying a topological charge which is a conserved quantity [19,22,77]. In a bulk system one can get rid of defects only via mutual annihilation [78,79] of their topological charges. The region, in which apparent deviations from the bulk ordering exist due to the presence of a defect is referred to as the core of defect [22,74]. Its linear dimension is roughly given by the correlation length ξ. The core typically exhibits a different phase with higher energy density with respect to the surrounding phase.
The coarsening dynamics following the quench evolves via three qualitatively different stages. Immediately after the quench the early stage is entered in which the orientational order parameter OP F ≡ S exhibits an exponential growth [4,76]. The early stage terminates roughly at the so called Zurek time [4], at which S ∼ S eq , where S eq determines the equilibrium value of S. At this stage topological defects are formed and consequently domains become visible [4]. Then the domain regime is entered [76]. In it the domain-type pattern in SBF is clearly pronounced. The domain pattern is well characterized by a single domain length ξ d (t). The subsequent domain growth gradually enters the dynamic scaling regime, where the power law ξ d (t) ∝ t γ is obeyed [75,76]. The universal scaling coefficient γ strongly depends whether a conservation law for the order parameter exists or not [75]. The order parameter evolution across a domain wall obeys the geodesic rule, i.e., it follows the shortest possible path in the order parameter space [4]. The domain growth is enabled by annihilation of topological defects. After some time the late stage regime [76] is entered in which the pattern appearance is dominated by individual topological defects. Finally, in bulk samples the defectless homogeneous structure is reached.
The above described mechanism of domain formation is universal and is referred to as the Kibble-Zurek mechanism [4,80]. It was originally developed to explain coarsening dynamics in the early universe [80]. The only essential ingredients of this mechanism are (i) continuous symmetry breaking and (ii) causality. It has been shown that the size ξ (0) d of protodomains, which are first formed at the Zurek time, depends on the quench rate. For example, for the temperature driven 2nd order phase transitions at T = T c the scaling law is obeyed [4]. Here τ Q stands for the quench rate. The universal coefficients v and η describe critical behavior of the relevant order parameter correlation length ξ ∼ ξ 0 |r| −v and the order parameter relaxation time τ ∼ τ 0 |r| −η close to the transition as a function of the reduced temperature r = (T − T c )/T c . The quantities ξ 0 and τ 0 are bare relaxation length and relaxation time estimating size of ξ and τ at relatively low temperatures. The scaling law (30) is obeyed only for fast enough quenches. Theory [81] suggests that the condition should be obeyed, where T G stands for the Ginzburg temperature.

Imry-Ma Domain Pattern
Next, we estimate characteristic size ξ d of domains that are stabilized by random field type of disorder following original ideas of Imry and Ma [12]. The estimate is based on static ground, i.e., it is not influenced by the history of a sample. For illustration purposes we consider the nematic LC phase perturbed by a random anisotropy-type of disorder. To estimate ξ d we balance the elastic free energy costs and free energy gains if the system obeys the imposed random field tendencies. Roughly speaking, if the domains are too small, the system possesses a large number of boundaries, whose energy is unfavorable. But if they are too large, they cannot order locally to take advantage of the local random fields. The compromise is a universal domain pattern, which is characterized by short range ordering.
We express the elastic free energy and the random anisotropy field density [17,71] costs as f e ∼ K 2 |∇ − → n | 2 and f RA ∼ − w 2 P 2 ( − → n . − → e ), respectively. Here we assume that the unit vectors − → e ( − → r ) are spatially randomly orientationally distributed and w measures the random anisotropy field strength. We assume that the system breaks up into a domain pattern characterized by ξ d . Typical elastic distortion in an average domain is estimated by f e ∼ K 2ξ 2 d , where < ... > denotes spatial average within a domain.
In order to estimate f RA we set that inside the domain − → e ( − → r ) undergoes random spatial variations. We label an average distance, on which − → e significantly changes, by ξ r and estimate the averaging rate within this domain. On traversing a domain volume LC molecules experience [12,64] random alternations in − → e orientation. We refer to N d as the averaging number. For a large enough domain, where − → n is essentially homogeneously oriented along a symmetry breaking direction, the term P 2 ( − → n . − → e ) averages to zero. According to the central limit theorem the averaging effectiveness is proportional to 1/ √ N d . Therefore one expects Balancing free energy costs f RA and f e yields the Imry-Ma scaling law [12] ξ d ∝ w 2 d−4 (34) Simulation Results The Imry-Ma estimation of ξ d is based on statical grounds reflecting the balance between the elastic and random field free energy penalties. On the other hand the Kibble-Zurek mechanism suggest that the size of protodomains depends on the quench rate. It is known that impurities could pin [82] regions exhibiting relatively strong elastic distortions. Therefore, in a configuaration reached via a phase transition quench the (quasi) stable pattern migh depend on the relative size of protodomains and the average distance among neighboring impurities. The size of protodomains depends on the quench rate [4] and consequently the domain pattern could depend on history of samples. Note that there have been published contradicting results related to the fate of low symmetry phase in presence of disorder. Some authors predict presence of Imry-Ma clusters [12,83], and others algebraic decay of correlations [2].
For this purpose [62,84,85], a systematic analysis have been recently carried out focusing on the impact of system history on randomly perturbed nematic phase. Simple minimal models in d = 2 and d = 3 were used in order to test universality of the Imry-Ma scaling prediction. Due to time consuming simulations and relatively wide phase space to be explored most simulations were realized at T = 0. Some predictions of these simulations were then tested at finite temperatures to confirm validity of obtained results for more realistic conditions. In the following we present the main results of these studies.
The Imry-Ma scaling and memory effects at temperature zero : The minimal model used is the following. Consider an ensemble of N cylindrically symmetric particles within a cubic lattice of unit lattice size a 0 . Local orientational ordering of a particle at the i-th site is given by an unit vector − → S i , to which we refer as a spin. The ± − → S i orientations are equivalent in order to mimic mesoscopic symmetry of LC molecules. The spins experience a random anisotropy field enforced by impurities, the concentration of which is labeled by p. The impurities are randomly distributed within the lattice and locally force spins to orient along − → e i . The orientational distribution of − → e i is isotropic.
The interaction energy F of the system is expressed as The short range interaction J ij is equal to J > 0 for the first neighbors and is set to zero otherwise. It tends to orient spins parallel or antiparallel. At randomly chosen sites (their number is equal to N p) a finite anchoring strength is imposed w i = w, representing impurities. At the remaining sites it holds w i = 0. The last term takes into account presence of an ordering external electric or magnetic field − → e z represent unit vectors of the Cartesiam coordinate system {x, y, z} and the constraint − → S i = 1 is imposed. Periodic boundary conditions are implemented.
In simulations one originates either from spatially randomly distributed orientations of spins, or from a homogeneously aligned sample along a single direction. We refer to these system histories as the random and homogeneous initial condition, respectively. The random initial configuration simulates sudden temperature quench from the isotropic phase. The homogeneous initial configuration can be realized by applying first a strong homogeneous external electric or magnetic field in the nematic phase. After a well enough alignment is achieved the field is switched off.
A system configuration is obtained by minimizing F with respect to spin components [84]. Therefore, the role of thermal fluctuations is neglected. Such an assumption is sensible deep in the nematic phase. In order to check validity of the Imry-Ma theorem d = 2 (2D) and d = 3 (3D) systems are considered. In 2D the simulations are constrained to the {x, y} plane and one sets S From obtained configurations the orientational correlation function G(r) is calculated. In 2D and 3D it is defined as respectively. The brackets denote the average over all lattice sites that are separated for a distance r. If the spins are completely correlated (i.e., homogeneously aligned along a symmetry breaking direction) it follows G(r) = 1. On the other hand G(r) = 0 reflects completely uncorrelated spins. In order to obtain structural details from calculated dependencies G(r) the ansatz is used, where ξ d , m, and s are adjustable parameters. The quantity ξ d estimates the average domain length (the coherence length) of the system. Over this length the spins are relatively well correlated. The distribution width of ξ d values is measured by m. Dominance of a single coherence length in the system is signaled by m = 1. A magnitude and system size dependence of s reveals the degree of ordering within the system. The case s = 0 indicates the short range order (SRO). A finite value of s reveals either the long range order (LRO) or quasi long range order (QLRO). To distinguish between these two orderings a finite size analysis should be carried out. If s(N ) saturates at a finite value the system exhibits LRO. If s(N ) exhibits algebraic dependence on N the system possesses QLRO. Representative G(r) dependencies for the random and homogeneous initial conditions in 2D and 3D for B = 0 are shown Figure 5. For the random condition in all studied cases [85] SRO was obtained, i.e., s = 0. On the contrary configurations obtained from the homogeneous initial condition yield s > 0 at a fixed value of N for weak anchoring strengths w and low enough concentrations p. This suggests either LRO or QLRO.
To check the validity of the Imry-Ma scaling the obtained ξ d values were fitted with the ansatz Figure 5. G(r) for homogeneous and random initial configurations in 2D and 3D systems, p = 0.3, w = 2. Note the difference between saturated G(r) dependence for random initial conditions. In the random case the correlations vanish (G(r → ∞) = 0) and for homogeneous initial conditions G(r → ∞) reaches a finite plato. Results show that in the strong anchoring limit (w > 10) ξ d does not depend on the history of the system. In the weak anchoring regime the Imry-Ma scaling was perfectly obeyed only for random initial conditions. In 2D and 3D the simulations yield ξ d ∝ w −1±0. 15 and ξ d ∝ w −2±0.3 , respectively. On the contrary, for the homogeneous initial configurations significant departures from the Imry-Ma behavior are observed. In this case SRO is not obtained for weak enough anchoring conditions, although systems still reveal a characteristic length as already suggested by Giamarchi and Doussal [15]. The simulations in this case suggest [85] 25 for 3D, see Figure 6.
It was further analyzed [84] behavior of systems using standard zero field cooled and field cooled cycles which are used in order to probe glassy features of systems. For this purpose simulations were started from random initial conditions. The zero field cooled configurations were calculated for B = 0. To obtain the field cooled configuration one first calculates metastable state in the presence of B. Afterwards B is switched off and again the quasiequlibrium state is searched for. The corresponding values of s, reflecting range of ordering in these cycles, are presented in Figure 7 in 2D and 3D. In both cases qualitatively similar behavior is obtained. Note that the zero field cooled configurations yield s = 0. One sees that the field cooled configuration possess a finite value of s. With increasing value of the external field, s in the first stage gradually increases and then saturates at a fixed value. These simulations signal pronounced memory effects in systems of interest. Finite temperatures : In order to verify memory dependent size of domains at finite temperatures a Brownian molecular dynamic study [76,86] was performed [62]. The domain patterns formed by rodlike particles, representing a LC molecule, was studied as a function of the concentration p of impurities imposing a random anisotropy field-type disorder to the particles and the history of the samples. The interaction between the neighboring spins − → S i is given by the modified Lebwohl-Lasher potential [87] Here r = | − → r ij | = | − → r i − − → r j | denotes the separation between the i-th and j-th spin in the 3D cubic lattice, J is a positive interaction constant and the parameter η describes the degree of orientational anisotropy. The case η = 0 is equivalent to the Lebwohl-Lasher [88], also referred to as the Maier-Saupe lattice model, which corresponds in the continuum limit to the approximation of equal Frank nematic elastic constants [42]. For η = 1 one gets the induced-dipole-induced-dipole type potential. However, the interaction f ij yields the nematic-like properties only for η < 0.3 [87]. Despite its simplicity, the Lebwohl-Lasher model well mimics main static and dynamic properties of a typical isotropic-nematic (I-N) phase transition. By studying cases η > 0 one probes the impact of elastic anisotropies on the system properties. The system is enclosed within a cube of volume N a 3 0 , where N is the number of molecules in the system and a 0 is the characteristic size of the unit cubic cell of the lattice. At the systems' boundary the periodic boundary conditions are imposed. Impurities of concentration p are assigned to randomly chosen lattice sites − → r (0) i of the system, where the unit cell is the simple cubic one with the lattice constant a 0 . It is assumed, that the i-th impurity enforces the orientation − → e i to the neighboring j-th spin via the short range potential The orientations of unit vectors − → e i are randomly distributed in 3D, and w stands for the orientational anchoring constant. The interaction energy W int of the whole sample is given as a sum over all pair interactions. In calculations one limits to the interactions with neighbors within a sphere of radius The orientation of the i-th spin is parameterized in the laboratory frame as where ϑ i = ϑ( − → r i , t) and ϕ i = ϕ( − → r i , t) represent dynamic variables of the model. The rotational dynamics of the system is driven by the Brownian molecular dynamics [76,86]. At each time interval ∆t (one sweep) the molecular orientation at the i-th site is updated in the local frame − → r = (x , y , z ) obeying the equations [86] ϑ in which the orientational diffusion tensor is diagonal. Its eigenvalues are assumed to be degenerate, equal to D, and k B is the Boltzmann constant. The z -axis of the local frame is oriented along the long axis of the LC molecule. The angles ϑ    In Figure 8 the time evolution of the domain coarsening in the pure bulk sample (i.e., p = 0) is shown following a sudden isotropic-nematic phase transition. Soon after the quench the scaling regime is entered, where ξ d (t) ∝ t 0.49 . Simulations show negligible influence on η. The domain length time evolution ξ d (t) for different concentrations p for different histories of systems is shown in Figure 9. The system was either quenched (i) from the isotropic phase, and (ii) from perfectly homogeneously aligned sample. The latter case corresponds to a sudden switch-off of a strong external magnetic or electric ordering field. For these histories the ξ d (t) values monotonously (i) increase and (ii) decrease with time, respectively. The saturated values ξ d of patterns are shown in Figure 10. They exceed the average distance between impurities and show significant memory dependence in line with predictions obtained using much simpler approach described in the previous subsection. The patterns originated from isotropic configurations have shorter values of ξ d in comparison to the structures reached from homogeneously aligned samples. In all cases the scaling ξ d ∝ 1/p is well obeyed. This observation in line with experimental results obtained in mixtures of nematic LC molecules and aerosil particles [83], for which the model described is appropriate.

Phase Behavior
In the following we assume that the random anisotropy field gives rise to a domain type pattern. We analyze phase expected behavior for perturbed nematic LC and nematic-non nematogen mixture.

Phase behavior of RAN
The random anisotropy nematic type disordering field enters the expression for the free energy F (18) via the surface contribution f s (26). In the case of a liquid crystal confined in a porous medium or in a mixture composed of liquid crystal and a second component (impurities, or colloids, or carbon nanotubes), due to essentially randomly spatially varying liquid crystal-substrate interface orientation, the easy axis e s exhibits random spatial variations. If a relatively weak disorder is enforced to a liquid crystal phase, then according to the Imry-Ma prediction the competition between the ordering and random field disordering tendencies can result in a domain-like pattern [12]. Accepting this prediction, according with the central limit theorem, the surface contribution to the free energy becomes: where ξ d is the average linear length of the domain, ξ r is a typical scale on which the easy axis undergoes random spatial variations, and N ∼ (ξ d /ξ r ) 3 . Note that ξ d is also the length on which a typical director field distortions are expected to evolve if other external fields are absent. On the contrary locally induced perturbations in S are expected to take place at the liquid crystal-perturber interface on the nematic correlation length scale ξ n (T ), which is a microscopic quantity (typically in the range of nm). Therefore in the mesoscopic limit we can consider only the spatial variations of the director n, and for a distorted nematic one typically expects that the elastic free energy density is given by For numerical purposes we introduce instead the length ξ = ξ 2 d − ξ 2 r . Therefore the case ξ = 0 corresponds to ξ d = ξ r .
Taking into account these previous considerations, the free energy density can now be written as For convenience we introduce nondimensional quantities. The temperature is replaced by the reduced temperature τ = (T − T * )/(T N I − T * ) where T N I = T * + B 2 /24aC. In this new temperature scale τ N I = 1, τ * = 0, and τ + = 9/8. The nematic order parameter is normalized with respect to its value at the nematic-isotropic phase transition of the homogeneous systemS = 6CS/B and the nondimensional free energy density is defined asf = 24 2 C 3 f /B 4 . The nondimensional length, elastic constant and surface anchoring are given byξ = ξ/ξ r andL = 8CL/B 2 ξ 2 r and Λ = 48C 2 W/ξ d B 3 . Omitting the bar notation, the dimensionless free energy density is written in the following form: Taking the limit ξ → ∞ in Equation (48) the undistorted bulk free energy density of the nematic phase is reproduced. The finite value of ξ reflects the competition between the elastic term L 1+ξ 2 S 2 (favoring ξ → ∞) and the surface field term − Λ (1+ξ 2 ) 3/4 S (favoring ξ → 0). Minimizing Equation (48) with respect to ξ leads to where S c = 3Λ/4L. The corresponding free energy density is given by The subscripts P and N stand for the paranematic and speronematic ordering, respectively. The paranematic phase closely resembles the isotropic phase but exhibits a finite degree of nematic ordering. The speronematic phase represents a distorted nematic phase, i.e., the phase with a finite value of ξ. The (Λ, τ ) phase diagram for L = 1 is shown in Figure 11.  Figure 11) marks a continuous transition characterized by a discontinuity in the slop dS/dτ . Note that at this line the switching between the two solutions for ξ (49) takes place. In the literature we have not found any experimental evidence resembling the observed slope discontinuity. Therefore, we believe that this anomaly is probably an artifact of our method. In reality a gradual evolution of the nematic ordering is expected.
The behavior of the correlation length ξ and of the order parameter S as functions of reduced temperature τ are shown in Figures 12 and 13 respectively.    For low values of the order parameter S (i.e., S < S c ), the length ξ is zero and the director orientations are correlated only on ξ r . At larger values of S the cost of changing the orientation from point to point is increased. The system responds by increasing the length scale ξ over which the director orientations are correlated. In the low disorder limit Λ → 0, the true nematic phase with infinite ξ is recovered. We mention that ξ is a quantity of key importance for understanding the phase transition in the presence of disorder. It is in fact a direct measure of resistance or the stability of the nematic phase in the presence of quenched disorder. In principle, ξ could be larger than the typical length over which the disordered structure is correlated. As expected, ξ increases as τ decreases, and indeed as S increases.
The profiles of the order parameter for two values of Λ shown in Figure 13 indicate that the bulk first order nematic-isotropic phase transition is weakened and eventually suppressed by increasing randomness.

Nematic-non-Nematic Mixture
In this section we present some results concerning the influence of a random anisotropy type disorder on the phase separation of the nematogen-non-nematogen mixture. We assume that the impurities (i.e., the non-nematic component) via the impurity-liquid crystal interface orientational anchoring interaction enforce to the liquid crystal phase a kind of random anisotropy field [72]. A combination of the random anisotropy nematic model and Flory-Huggins [89] model is used for the theoretical study of the phase stability of the mixture as a function of the absolute temperature T , impurity volume fraction Φ, and random anisotropy field strength Λ. We show that the random anisotropy field qualitatively changes the topology of the (Φ, T ) phase diagram.
The mixture is characterized by the volume fractions of the two components: where N i is the number of molecules of component i (i = 1 defines the liquid crystal and i = 2 the non-nematogenic fluid) and v i is the volume of a particle of component i. In what follows, we consider that the volumes of a nematic and impurity molecule are equal (v 1 = v 2 = v) and Φ 2 = Φ (due to the conservation of the number of particles Φ 1 = 1 − Φ). The orientational order of the mixture is characterized by the nematic order parameter Q αβ which, in the case of uniaxial nematic state, is given by Equation (21). The free energy F of the system is expressed as where V is the volume of the system and A is the surface area of the liquid crystal-impurity interface. The first term is the free energy density of the isotropic mixing for the two components [89] where k B is the Boltzmann constant and χ = (U 0 /k B T ) is the Flory-Huggins interaction parameter related to isotropic interaction between unlike molecular species. The second term in Equation (52) is the Landau-de Gennes free energy density of the isotropicnematic phase transition generalized to include the interaction between nematic and non-mesogenic molecules: The coupling between Φ and Q αβ in Equation (54) results from microscopic considerations [62,90]). The elastic free energy density f e and the anchoring term f s have been discussed in sections (3.2.) and (4.2.) and are given by Equations (23) and (24) respectively.
Using the non-dimensional variables defined in section (4.2.) the dimensionless free energy density is given by Note that the local ordering tendency of impurities induces a finite value of S even in the isotropic phase. One commonly refers to such the case as the paranematic ordering. The equilibrium condition at the nematic-paranematic phase transition is given by the following set of Equations [91]: where is the difference in grand potential density between the two phases and is the chemical potential. The ordering in the nematic and paranematic phase is determined by pairs (Φ, S) and (Φ p , S p ), respectively.
Phase behavior in the absence of randomness : We first consider the phase behavior of the system in the absence of the random filed. In this case the elastic and surface terms are absent (i.e., f s = f e = 0 and Equation (55) becomes: Figure 14. Phase diagram for nematic -non-nematic binary mixture in the absence of disorder for Γ = λ = χ = 1. In this case the paranematic phase is replaced by the isotropic one and S p = S i ≡ 0, where S i stands for the degree of ordering in the isotropic phase. Taking into account Equations (56) we calculate a representative phase diagram of the nematic -impurity (i.e., nematic-non-nematic) mixture, that is plotted in Figure 14 for Γ = λ = χ = 1. The solid curve refers to the binodal, that constitutes the actual phase boundary. For temperatures below T N I = 317.5K (value corresponding to pure 5CB [45]), there exists a two-phase coexistence region (N+I) between the isotropic (I) and nematic phase (N). On decreasing the temperature the biphasic region broadens.
To show how the onset of nematic ordering triggers phase separation in conventional liquid crystal, we rewrite f Equation (59)) in the following form where τ ef f = τ and χ ef f = χ + λS 2 /Γ stands for the effective Flory-Huggins interaction parameter. We see that the nematic ordering effectively increases the Flory-Huggins interaction parameter, that can potentially lead to order-induced phase separation.
Phase behavior in the presence of random field: In this section some theoretical results concerning the influence of disorder Λ > 0 on phase behavior of thermotropic nematic-non-nematic mixture are presented. The phase separation in the presence of random field takes place between the paranematic and speronematic phases. The paranematic phase closely resembles the isotropic phase but exhibits a finite degree of nematic ordering and ξ = 0. The speronematic phase represents distorted nematic phase that is characterized by a finite value of ξ (for an ordinary nematic phase ξ → ∞).
Minimizing Equation (55) where S c = 3ΛΦ/4L. The corresponding free energy density is given by, The subscripts P and N stand for the paranematic and speronematic ordering, respectively. To calculate the phase diagram we use the equilibrium conditions given by the Equations (56). The solution of these equilibrium conditions is determined by four quantities: (i) the paranematic phase is determined by Φ p and S p and (ii) the speronematic phase by Φ s and S s .
In the vicinity of T N I the phase diagram is very similar to that in the absence of random field. But at some value of temperature T − T N I = −1.6K (which depends on the value of Λ), the paranematic order parameter does not fulfill the condition S p < S c = 3ΛΦ/4L. As a consequence, the system is in only one state, the speronematic phase and correspondingly the phase separation cancels. This fact can have important consequences in the thermodynamic properties of the mixture, for example can explain the double pick in the specific heat experimentally observed in similar systems. This work is now in progress.

Experimental Observations
The described phenomenological RAN model roughly describes behavior of nematic LC experiencing a quenched random anisotropy field disorder. In the following we show some experimental evidences supporting this approach and discuss in which studied physical properties the influence of disorder might be significant for a temperature interval covering isotropic and nematic ordering in bulk LC samples.
In experimental studies the disorder is typically introduced geometrically by different perturbers via a varying LC-perturber interface [16]. For perturbers one commonly chooses inert porous matrices hosting LC phases or networks formed by aerosil [24,25] particles in the aerosil-LC mixtures. As matrices aerogels and Controlled-Pore Glasses (CPG) [31] are conventionally used. Aerosil mixtures [26] are particularly useful because one can obtain qualitatively different random field-type regimes by changing the concentration of aerosil particles.
All these perturbers introduce into systems a new characteristic length scale R. In the case of aerosil nanoparticles this is the mean aerosil void [24] size R ∼ 2/(aρ s ), where ρ s stands for the mass density of the aerosil particles and a is their active surface. In case of porous materials R represents a characteristic void radius. On commonly assumes that with decreasing value of R the disorder strength monotonously increases. However, a more detail analysis suggests that the dependence between R and the disorder strength is more complicated, at least in LC-aerosil mixtures [92]. Note that values of R are in most samples well below 1µm. In such systems, confinement, surface wetting and finite size effects could also play a significant role.
Most studies have been performed in mixtures of LC phases and aerosil nanoparticles. The spherular aerosil particles have diameter 2R s ∼ 7 nm. The specific surface area of aerosils [24,26] is a ∼ 300 m 2 /g for hydrophilic (type 300), and a ∼ 280 m 2 /g for hydrophobic (type R812) aerosils. The neighboring aerosils could be linked via covalent bonds forming networks exhibiting moderate branching. By varying the concentration of aerosils ρ s three qualitatively different networks [24,26] can be obtained. For ρ s < 0.01 g/cm 3 the aerosil particles are more or less randomly distributed in the system. For ρ s > 0.01 g/cm 3 the aerosils form a gel-like thixotropic network structure. In the so called annealed regime, ranging roughly between ρ s ∼ 0.01 g/cm 3 and ρ s ∼ 0.1 g/cm 3 , the aerosil network is relatively responsive. It can rearrange in order to partially anneal the elastic stress imposed by the surrounding LC phase. In the stiff regime, corresponding to ρ s > 0.1 g/cm 3 , the aerosil network becomes a rigid-like, imposing a quenched type of disorder to surrounding LC molecules.
Several investigations have been also carried out in LC confined to aerogel matrices. Aerogels [27] form random networks of silica backbones in an open connected void space. The void pore size distribution is rather broad and the average characteristic pore size R can be far below 1 µm, depending on the density ρ a of the gel. The porosity of aerogels can be extremely high therefore providing host lattices of relatively high surface to volume ratio. The porosity can be varied between 0% and 98% by partial desification [93].
X-rays studies [24] reveal that the aerogel matrices and aerosil networks exhibit similar structures for high enough densities. Roughly equivalent structures obey equality [94] ρ s ∼ ρ a φ p where φ p stands for the aerogel pore volume fraction. Furthermore, aerogels-LC and aerosil-LC mixtures are transparent and can be therefore studied also by means of optical spectroscopy.
Much work focusing on the influence of disorder on LC behavior has also been carried out in Controlledpore glasses (CPG) [32,95,96]. A CPG matrix [32] consists of strongly curved and interconnected voids, that introduce a certain degree of randomness into the systems. Each void resembles a curved cylinder of radius R with a rather narrow distribution of R-values. Samples with radia from several nm up to µm can be obtained.
In all these samples the orientational anchoring condition at the perturber-LC interface depends on interface treatment and also on the type of LC used. However as a rule, LC-perturber interfaces tend to increase degree of LC ordering locally [96]. Consequently, e.g., above the bulk I-N phase transition temperature T IN the isotropic phase is replaced by the paranematic phase, exhibiting a finite degree of nematic orientational ordering. On the other hand disorder gives rise to frustration in orientational ordering and therefore suppresses the degree of ordering. We will henceforth refer to the randomly frustrated nematic phase as the speronematic phase.
Regarding the range of ordering of the disordered phase, most studies confirm the prediction that the broken phase exhibits a domain-type short range order (SRO) in LC ordering. SRO was reported in samples using aerogels [3], CPGs [95] and aerosil particles [24] as LC perturbers. Most elaborate experimental investigations confirming the Imry-Ma hypothesis was performed by Bellini et al. [83] in LC-aerosil mixtures. Note that several analytic studies predict that quasi long range ordering is established instead or even LRO for weak enough degree of disorder [2]. However, these predictions have a relatively weak experimental support [97].
In randomly perturbed LCs on commonly observes a suppressed value of the paranematic-speronematic phase transition temperature T  [32,71,95]. It has been shown that in addition to randomness also surface ordering and in some cases even finite size effect should be taken into account [95]. In particular surface local ordering gives rise to enhanced LC ordering at the LC-perturber interface the thickness of which equals few nm. Such ordering could be in some cases well reproduced using a relatively simple bicomponent model [71], as demonstrated in the T (p) IN (R) analysis in CPG and aerosil samples. In this model the average LC ordering is described by two order parameters, describing the average degree of ordering within the LC-perturber interfaces and in the remaining LC body.
Furthermore, for small enough values of R, i.e., for a strong enough disorder, the 1st order speronematicparanematic is replaced by gradual evolution of nematic ordering on varying temperature. The corresponding critical value of R ≡ R c is in general comparable to the nematic correlation lengh at T IN . This effect has been observed in LC-aerosil [24], LC-aerogel [27] and LC-CPG samples [96]. More detail analyses of RAN-type models yield approximate reationship between disorder anchoring strength and material-geometrical characteristics of samples. Furthermore, numerical estimates give reasonable quantitative estimates on critical values of R c [71,96].
Note that the transition from the critical to the noncritical LC temperature behavior on decreasing R could be also triggered by an ordered surface due to noncritical character of the LC-interface coupling [98]. However, in such case T

Concluding Remarks
We have presented typical phase and structural behavior of randomly perturbed systems exhibiting continuous symmetry breaking using nematic liquid crystal phase as a testing ground. In liquid crystalline materials several features can be experimentally probed due to relatively good experimental accessibility, which can be exposed to different strengths and types of disorder. The impact of disorder on liquid crystal ordering is pronounced due to the softness of liquid crystal phases.
We have reviewed behavior of a nematic liquid crystal experiencing quenched random anisotropy field disorder assuming domain-type orientational ordering. Typical ordering of such systems has been demonstrated using relatively simple numerical simulations based on Lebwohl-Lasher lattice type model. It has been shown that a domain-type pattern characterized by a single characteristic length in the symmetry breaking field is inevitably formed at least temporally in the lower symmetry state (nematic) after a quench from the higher symmetry phase (isotropic) even in the absence of disorder. Such domain pattern can be stabilized by impurities. If impurities impose a random field-type disorder a domain pattern again characterized by a single characteristic length obeying the Imry-Ma scaling can be formed. We have presented cases in which the Imry-Ma scaling is expected and have shown that memory effects could be pronounced.
Assuming the domain-type nematic structure, we have used the Random Anisotropy Nematic (RAN) phenomenological model to study the phase behavior of nematic (speronematic)-isotropic (paranematic) phase transition. It has been shown that the main features which such a simple approach predicts are indeed realized in several experimental systems, i.e., liquid crystal immersed in Controlled Pore Glasses, liquid crystal-aerosil dispersions, and other liquid crystal-impurities mixtures. The review indicates that such relatively simple RAN-type approaches are useful to predict qualitative behavior and make quantitative estimates on several physical properties of randomly perturbed systems, the phases of which are reached via a continuous symmetry breaking phase transition.