# Numerical Results for the Blue Phases

###### Abstract

We review recent numerical work investigating the equilibrium phase diagram, and the dynamics, of the cholesteric blue phases. In equilibrium numerical results confirm the predictions of the classic analytical theories, and extend them to incorporate different values of the elastic constants, or the effects of an applied electric field. There is a striking increase in the stability of blue phase I in systems where the cholesteric undergoes helical sense inversion, and the anomalous electrostriction observed in this phase is reproduced. Solving the equations of motion allows us to present results for the phase transition kinetics of blue phase I under dielectric or flexoelectric coupling to an applied electric field. We also present simulations of the blue phases in a flow field, showing how the disclination network acts to oppose the flow. The results are based on the Landau-de Gennes exapnsion of the liquid crystal free energy: that such a simple and elegant theory can predict such complex and subtle physical behaviour is remarkable.

###### pacs:

61.30.Mp, 64.70.mf, 83.80.Xz, 61.30.Gd## I Introduction

One of the spectacular successes of de Gennes’ Landau theory of liquid crystals is that it captures, without need for extension or modification, many of the remarkable and subtle phase transitions and properties of liquid crystals. Blue phases provide an especially vivid example of this. The general features and properties of these unique mesophases are both ‘beautiful and mysterious’ degennes93 () and yet they can be understood qualitatively and quantitatively using de Gennes’ -tensor theory of liquid crystal phase transitions.

Blue phases are found in highly chiral liquid crystals between the high temperature isotropic fluid and the lower temperature cholesteric phases reinitzer88 (); degennes93 (); blumel85 (); yang87 (); crooker89 (); belyakov85 (); wright89 (). They are remarkable mesophases, exhibiting a brightly coloured texture of individual, micron sized platelets. The bright colour indicates selective reflection due to a periodic structure, much like in an ordinary crystal, but with a much larger characteristic length scale, and indeed the reflection spectra show Bragg peaks that can be indexed by cubic space groups with lattice constants of several hundred nanometers. Furthermore, individual platelets of monodomain crystals themselves show distinctive facetting corresponding to Miller planes of the lattice structure and with the faces growing in a sequence of steps pieranski85b (). Yet, blue phases are not crystals in the traditional sense: they have no long range positional order and are full three-dimensional fluids. The crystalline order is in the orientational degrees of freedom of the liquid crystal.

The key to understanding the properties of the blue phases was in realising that the locally preferred order in the chiral liquid crystal is one of double twist, with the molecules adopting helical ordering along two, perpendicular axes, as opposed to the usual single twist of an ordinary cholesteric helix meiboom81 (); meiboom83 (); hornreich81 (); hornreich82 (). However, it is not possible to construct a global state with helical ordering in two directions, without introducing disclination lines into the structure. Therefore the blue phases consist of local, cylindrical, regions of double twist separated by a regular disclination network.

A tutorial, two-dimensional example of a blue phase texture is shown in Fig. 1. This demonstrates how local regions of double twist can be pieced together only by introducing a square array of topological defects. In physical, three-dimensional, crystals the disclination structures are more complex; three distinct blue phase textures have been identified upon cooling from the isotropic liquid. Typical experimental phase diagrams are shown in Fig. 2. Blue phases I and II exhibit textures with cubic symmetry corresponding to the space groups and respectively wright89 (). These are illustrated in Fig. 3 which shows both the disclination networks and the topology of the lattice of double twist cylinders in the two phases. Blue phase III is less well understood; it has an amorphous structure with the same symmetry as the isotropic fluid lubensky96 (); anisimov98 ().

The unique combination of crystalline order with lattice constants comparable to the wavelength of visible light and full three dimensional fluidity make the blue phases ideal for technological uses such as fast light modulators, photonic crystals or tunable lasers cao02 (); hisakado05 (); yokoyama06 (); kitzerow06 (). The principal obstruction to their application was the very limited temperature range, , over which the phases were stable. However, recently, the stability range has been extended to as much as , including room temperature, by the addition of bimesogenic molecules or photo-crosslinking of polymers kikuchi02 (); coles05 (). This has opened the way to device applications and, in May 2008, Samsung Electronics unveiled a prototype ‘blue phase mode LCD’ wikiBPLCD08 (); kikuchi07b ().

The theoretical description of the general features and properties of the blue phases was developed in the 1980s in a series of important papers, summarised in the review wright89 (). Two approaches were followed, a low chirality theory, based on the Frank director field description meiboom81 (); meiboom83 () and a high chirality theory using the Landau-de Gennes -tensor grebel83 (); grebel84 (). The theoretical approaches were able to correctly predict the symmetries of the blue phases and to give a good account of the phase diagram, although with some discrepancies. For example, additional textures with space groups and were found to be stable within the parameter range corresponding to the experiments, but were not observed experimentally.

The renewed experimental interest in the blue phases provided the motivation to revisit them theoretically. Since the original analytic investigations, computers have become a powerful tool for studying complex fluids, and it is now feasible not only to revisit and extend the calculations of the blue phase phase diagrams dupuis05a (); gareth06 (), but also to obtain results on blue phase kinetics and hydrodynamics dupuis05b (); gareth08 (). The aim of this article is to review recent research showing how numerical simulations of the blue phases are allowing us to gain further insights into their behaviour.

In Section II we describe the Landau-de Gennes equations that have proved so successful in modelling the blue phases and outline how they can be solved numerically. In particular we discuss how the size of the unit cell must be allowed to vary to correctly minimise the free energy. Then, in Section III, we revisit the classic calculations of the thermodynamic phase diagram showing the changes that result when the approximations necessary to make the analytic calculations feasible can be circumvented by a numerical solution. We also show that including additional terms in the Landau-de Gennes free energy can significantly change the phase diagram. Most dramatic among the results is the increase in the stability of blue phase I obtained in systems where the cholesteric undergoes helical sense inversion.

Many of the potential technological applications of the blue phases rely on the response of the material to an electric field. In the blue phases the effects of an electric field are particularly rich, including continuous distortions in the size and shape of the unit cell stegemeyer84 (); porsch86 (); porsch89 () and a series of field-induced transitions to new blue phase structures, not stable in zero field pieranski85 (); cladis86 (); pieranski87 (); chen87 (). In Section IIIC we show how the Landau-de Gennes expansion correctly predicts electrostriction and, in particular, the anomalous electrostriction of blue phase I.

We next, in Section IV, summarise the equations of motion which describe blue phase hydrodynamics. These can be used to investigate the kinetics of transitions between blue phases, and their response to an applied electric field or to an imposed flow. Section V is devoted to examples demonstrating the application of the dynamical equations. We describe an investigation into the viscoelastic properties of the blue phases as they are subject to a Poiseuille flow field. We then describe the kinetics of the phase transitions induced in blue phase I as it is coupled to an electric field, either via a dielectric, or a flexoelectric, term in the free energy. This allows us to propose a candidate structure for blue phase X.

## Ii Minimising the Landau-de Gennes free energy

### ii.1 Landau expansion in terms of the -tensor

One of de Gennes’ vital contributions to the theory of liquid crystals was to identify the tensor as a suitable order parameter degennes93 (). A Landau expansion of the free energy of liquid crystals in terms of has proved useful as a starting point for interpreting liquid crystal phase diagrams degennes93 (); wright89 (); gareth06 (); grebel83 (); brazovskii75 ()

(1) |

Here is a constant with the dimensions of an energy density, plays the role of an effective temperature for thermotropic liquid crystals, defines the helical pitch and the are elastic constants. The expansion, taken to second order in the derivatives of , can only account for two independent Frank elastic constants and imposes that the magnitude of splay and bend are equal wright89 (). In order to remove this constraint, and to allow for a temperature dependent helical pitch, it is necessary to consider, at least, terms cubic in and quadratic order in gradients berreman84 (); longa87 (); longa89 (). There are eight such terms and it is unrealistic to consider them all. We include one chiral and two achiral terms and . The term was chosen as it contributes equally to all three elastic constants and the term as it gives the largest distinction between splay and bend and because it has the largest contribution to the energetics of an isolated double twist cylinder. (However, the choice is arbitrary and other coefficients may be non-zero in any given material, and may have a considerable effect on the phase diagram.)

Rewriting the free energy in terms of dimensionless variables demonstrates that it depends only on the two dimensionless parameters gareth06 ()

(2) |

known as the reduced temperature and the chirality respectively, together with ratios of the elastic constants. In line with previous work, phase diagrams will be presented in terms of these parameters.

### ii.2 Numerical minimisation of the free energy

The free energy (1) is minimised by relaxing the -tensor according to a Ginzburg-Landau equation

(3) |

This equation can be solved using many different numerical approaches; the results presented here were obtained using a lattice Boltzmann algorithm denniston04 (); julia06 ().

To study the different blue phases it is necessary to implement appropriate initial conditions for the simulation. The -tensor is initialised using analytic expressions appropriate to the high chirality limit (), which act to define the symmetry of the chosen phase. For blue phase I grebel83 (); wright89 (); dupuis05a () we use

(4) |

where , with the lattice constant, and the other components are obtained by cyclic permutation. Similarly, for blue phase II the -tensor is initialised as grebel83 (); wright89 (); dupuis05a ()

(5) |

where and the other components are again obtained by cyclic permutation.

Under numerical evolution using Eq. (3) the system relaxes to the structure of the same symmetry that locally minimises the free energy. We are therefore able to obtain, for any value of the parameters, local minima of the free energy corresponding to each of the cholesteric and blue phases. The global free energy minimum was taken to be the smallest of these calculated local minima.

To achieve a full minimisation of the free energy it is necessary to set the correct unit cell size in the simulation. This is not known a priori, but rather depends on the magnitude of the order parameter, a quantity that is only determined by the numerical minimisation itself. Therefore we must introduce a means of determining, and setting, the unit cell size as the -tensor evolves during the simulation. It is possible to account for a change in unit cell size by rescaling the gradient contributions to the free energy and molecular field. This is accomplished in practice by changing the elastic constants to

(6) |

where a superscript ‘init’ denotes the initial value of a simulation parameter and is the appropriate rescaling factor, which is identical to the ‘redshift’ described in grebel83 (); grebel84 ().

To calculate the optimal size of the unit cell we note that, since the free energy is quadratic in gradients, it may be written formally in -space as

(7) |

where the coefficients and depend on the -tensor, but not on . The optimum wavevector is given by , and since the coefficients and are determined by the simulation it is straightforward to use these values to determine the exact value for the size of the unit cell at every timestep, thereby obtaining a full minimisation of the free energy.

## Iii The equilibrium phase diagram

### iii.1 Revisiting the analytic calculations

The phase diagram for chiral liquid crystals obtained for a selection of parameter values using the Landau-de Gennes free energy, Eq. (1), is shown in Fig. 4. Fig. 4(a), adapted from grebel84 (), shows the phase diagram calculated analytically in the high chirality limit, in the one elastic constant approximation, and Fig. 4(b) compares numerical results for the same parameters in the free energy. The differences between the analytic and numerical results show that including higher order harmonics in the minimisation of the free energy is significant, as might be expected for phases with small free energy differences. In the numerical (exact) minimisation is stable over a larger range of parameters and at lower chirality values, in better agreement with the experimental phase diagram, Fig 2. Moreover, the regions of stability of the and textures found in the early studies grebel83 (); grebel84 () are shifted to unphysically high values of the chirality dupuis05a (). This result is again consistent with experiment, where blue phase structures of these symmetries have not been observed.

### iii.2 Varying the elastic constants

We now compare the phase diagrams obtained as the ratios of the elastic constants are varied. To investigate the effect of the bend elastic constant we chose parameter values which corresponds to a ratio of splay to bend of about 0.5, while splay and twist remain degenerate. The resulting phase diagram is shown in Fig. 4(c). Comparing to the case of equal elastic constants the stability of blue phase I is seen to decrease quite significantly relative to the cholesteric phase while at the same time there is a small increase in stability over blue phase II. There is only a minor shift in the cholesteric-blue phase I phase boundary at the transition temperature, however, as the temperature decreases the shift becomes larger.

The value of the twist elastic constant is controlled by the Landau-de Gennes parameter . In most liquid crystals the twist elastic constant is smaller than either splay or bend. In order to match this, we constructed the phase diagram for parameter values , which is shown in Fig. 4(d). This choice of parameters resulted in a ratio of splay to twist of about 1.5, while splay and bend remained degenerate. Again we observe that the stability of blue phase I is reduced relative to the cholesteric phase by an amount similar to that seen by varying the bend elastic constant.

Finally, we consider the effect of the chiral cubic invariant on the blue phases. We chose parameter values of , which gives a ratio of bend to splay of about . Moreover, for these coefficients the cholesteric undergoes helical sense inversion huff00 (); slaney92 () at a reduced temperature of about . The phase diagram is shown in Fig. 4(e). What is remarkable is the dramatic increase in stability of blue phase I relative to the cholesteric phase. The region of stability has been increased down to chiralities as low as and at such low chiralities the phase boundary is essentially independent of for all . In addition, there is a very small region of stability for blue phase II located close to the isotropic transition (Fig. 4(f)). Since blue phase I is now stable over a much larger temperature range it displays a significant variation in unit cell size as the temperature is lowered, with the lattice parameter more than doubling between and .

### iii.3 Electrostriction

The response of the cubic blue phases to the application of an external field has been a topic of interest since the mid 1980s and continues to be so because of the importance to potential blue phase based devices kitzerow91 (); kitzerow06 (); kikuchi07b (). The principal features include electrostriction, a continuous distortion of the shape and size of the unit cell with increasing field, and a series of field induced textural transitions. The electrostriction involves a shift of the back-scattered Bragg peak of 5–10% and is quadratic in the field strength stegemeyer84 (); porsch86 (); porsch89 (). The direction of the shift changes sign with the sign of the dielectric anisotropy, but blue phase I also displays an unusual response referred to as anomalous electrostriction, where an expansion along the field direction is seen when the field is applied parallel to the direction, but a contraction for fields parallel to . At larger field strengths new blue phases appear. Three distinct field induced textures have been identified, possessing tetragonal, screw hexagonal and two-dimensional hexagonal symmetry with increasing field strength pieranski85 (); cladis86 (); pieranski87 (); chen87 ().

Much of the electric field behaviour has been understood theoretically via extensions of the Landau-de Gennes theory, including the qualitative features of the electrostriction and the field induced textural transitions hornreich85 (); lubin87 (); stark91 (); hornreich90 (); longa96 (); zelazna98 (). However, the approximations inherent in the analytic calculations limited the quantitative comparison that was possible and a number of features, including the anomalous electrostriction of blue phase I and the tetragonal field induced texture blue phase X, could not be accounted for zelazna98 (); hornreich90 ().

The main additional difficulty in minimising the free energy of the blue phases in the presence of an electric field gareth08 () is in accounting for electrostriction, as there is a change not only in the size, but also in the shape, of the unit cell as the field is applied. This distortion can be accounted for in two steps: first the shape of the simulation unit cell is fixed and the corresponding size, which minimises the free energy, is determined. This is then repeated for a set of varying shapes of the unit cell, e.g. from cubic to tetragonal. Hence the free energy is determined for several values of the ratio parameterising the cubic-tetragonal distortion and fitted to a quadratic. The actual distortion is then given by the miminum of this fit.

The field dependence of the blue phase lattice parameters is shown in Fig. 5 for an applied electric field along the direction. Note, in particular, that the unit cell expands along the field direction in blue phase II, but contracts in blue phase I, in agreement with experiment kitzerow91 (). When the field is instead applied along to the direction both blue phases undergo an expansion parallel to the field, a precursor to the transition to blue phase X that is observed at larger field strengths. It is very pleasing that a numerical approach can predict anomalous electrostriction in blue phase I as the effect is lost in the truncations needed in analytic calculations zelazna98 (). Mapping from physical to simulation units gives a magnitude of the elctrostriction in the range for both blue phases, again in good agreement with experiments kitzerow91 ().

## Iv Hydrodynamic equations

The Ginzburg-Landau equation, Eq. (3), describes the relaxation of a liquid crystal to the minimum free energy, but does not describe physical dynamics in situations where flow is important. Several authors have recently used numerical approaches to solve the full hydrodynamic equations of motion for liquid crystals in the nematic and cholesteric phases. For example, it has been possible to simulate defect hydrodynamics toth02 (); toth03 (); svensek02 (); svensek03 (), phase ordering denniston00 (); denniston01c (); nidhal06 (), the kinetics of transitions between different liquid crystal phases fukuda98 (); gareth08 (), cholesteric rheology davide04a (); davide04b (); orlandini05 (); davide05b (); davide06 (); orlandini07 () and the effect of flow on device switching qian97 (); qian01 (); denniston01a (); toth02b (); denniston02 (); davide03 (); davide05a (); stromer06 (); james08 (). Because of their disclination structure the kinetics and the rheology of the blue phases is an exciting, but demanding, numerical problem which requires intensive numerical resources. However, such simulations are rapidly becoming feasible: we summarise the equations of motion and some of the results obtained so far, and then discuss possible directions for future work.

The hydrodynamic equations of motion of liquid crystals are complex, both because of the anisotropy of the molecules and because of the coupling between the order parameter field and the flow field. In general, flow leads to a rotation of the local orientation, which in turn influences the flow. Similarly, if a disturbance is initiated in the director, its reorientation is generally accompanied by fluid motion, an effect sometimes referred to as backflow. If a system is close to a phase transition, or contains disclinations, variations in the magnitude of the order parameter can be significant and therefore a hydrodynamic description based on the -tensor is needed beris94 (); grmela97 (); ottinger97 (); ottinger05 (); qian98 ().

The order parameter evolves towards the minimum of the free energy, but with a convective time derivative to account for the advection with the fluid

(8) |

The term in brackets on the right hand side is called the molecular field, , and is a collective rotational diffusion constant. The material derivative for rod-like molecules is given by beris94 ()

(9) |

where and are the symmetric and antisymmetric parts, respectively, of the velocity gradient tensor . The constant depends on the molecular details of a given liquid crystal.

The fluid velocity field is taken to obey the continuity equation and a Navier-Stokes equation with a stress tensor generalised to describe liquid crystal hydrodynamics beris94 ()

(10) | |||

(11) | |||

(12) |

Eqs. (8), (10) and (11) can be solved either using a lattice Boltzmann approach denniston01d (); denniston04 (), or finite difference schemes svensek02 (); james06 (); qian01 () or, giving somewhat improved stability while retaining the advantages of lattice Boltzmann, a hybrid scheme where a lattice Boltzmann solution of the flow equation is coupled to a finite difference solution of the order parameter field davide07 ().

## V Blue phase rheology

### v.1 Response of the disclination lattice to a Poisueille flow

To investigate the rheological response of the blue phases we placed a unit cell between fixed plates and imposed constant force on the fluid, together with no-slip boundary conditions on the velocity field at the plates. In a Newtonian fluid this geometry leads to a quadratic, Poiseuille flow profile. The choice of boundary conditions for the director field was to assume that the disclinations are fixed at the boundaries. Fig. 6 compares the apparent viscosity (obtained comparing to the Poiseuille velocity) in blue phases I and II, in a phase comprising a square array of disclinations with two-dimensional cross section corresponding to Fig. 1, and in the isotropic phase. The corresponding disclination configurations, comparing zero and a finite velocity field, are shown in Fig. 7.

For small forcing the blue phase viscosities increase by a factor over that of the isotropic fluid. This is because the disclination network acts to oppose the flow and dissipate energy. The blue phases reach a stationary state in which the disclination network is bent and twisted by the flow. The viscosities of blue phases I and II are approximately constant over a range of forcing, but the square lattice structure shows shear thickening. This is because each defect line, of topological strength 1, opens to a disclination ring, comprising disclinations of strength 1/2. The ring then twists and bends on itself as the flow increases, as shown in Fig. 7.

As the forcing increases there is significant shear thinning in all three blue phase structures. This occurs because the disclination network is destroyed by the flow, and the viscosity drops to that of an isotropic liquid crystal.

These results indicate that blue phase rheology is extremely rich, and worthy of further study, both experimentally and numerically. There is a need for careful experiments, with good control over boundary conditions. On the numerical side, work is in progress to move towards larger numbers of unit cells in the simulations and to assess the effects of different flow geometries and boundary conditions henrich09 ().

### v.2 Blue phase X

The ability to solve the dynamical equations of motion means that it is possible to follow the way in which disclinations rearrange during continuous textural transitions. Moreover, the simulations provide a way to identify possible candidate structures for blue phases which have not yet been identified.

An example of this approach is simulations in which blue phase I is placed under an intermediate electric field applied parallel to the direction. This is the set-up corresponding to the blue phase I-blue phase X transition observed experimentally in the 1980s cladis86 (); pieranski87 (). The electrostriction distorts the shape of the unit cell until it becomes tetragonal, at which point there is a transition from the blue phase I texture to a new texture known as blue phase X. Starting from this geometry, and an unperturbed blue phase I, the evolution of the disclination network is simulated numerically. Initially, the disclinations in the network twist, they then merge to form a transiently branched structure, which finally reorganises into a new defect network, not stable at zero field as shown in Fig. 8. This is a candidate structure for blue phase X, as it is a new network, found via a continuous reorganisation starting from blue phase I, and only stable in a field. The results predict that in the candidate blue phase X (i) the disclinations perpendicular to the field are largely unaffected, (ii) the network conforms to the space group identified in Reference pieranski87 (), and (iii) the double twist cylinders deform but do not break during the transition; this observation may be verified by experiments along the lines of those in References kikuchi07 (); higashiguchi08 ().

### v.3 Flexoelectric blue phases

In 1969 Meyer introduced the concept of flexoelectric coupling to an external electric field and showed that this could lead to a one-dimensional splay-bend distortion of the nematic director field meyer69 (). These results have recently been extended to show that, near the isotropic-nematic transition and with sufficiently strong coupling, two-dimensional splay-bend structures with hexagonal symmetry can be stable gareth07 ().

Fig. 9 shows the evolution of blue phase I when flexoelectric coupling to an applied field is increased quasistatically. As in the dielectric case, the electric field induces a twist in the disclination lines allowing them to transiently merge, thus fascilitating textural transitions. For the largest applied field strengths a transition is indeed observed to the two-dimensional hexagonal flexoelectric blue phase, as expected. However, at intermediate values, two further transitions precede it, each yielding distinct textures that are stable over a small range of field strengths. The first is to a centred tetragonal texture with space group . This has the same disclination network as the blue phase X structure found under dielectric switching and described above. At slightly higher field strengths there is a second transition to a distinct tetragonal texture with space group . In this texture, the disclinations are all parallel to the field, occuring in pairs that wrap around each other to form a double helix, and with the axes of the double helices themselves then arranged on a square lattice. Finally, as the flexoelectric coupling is inceased, the two members of a given double helix transiently merge and re-separate, allowing them to both straighten out and to adopt an hexagonal configuration in the plane perpendicular to the field.

## Vi Discussion

The Landau-de Gennes expansion of the free energy of liquid crystals has proved a vital tool in understanding their thermodynamics. The theory has proven successful even when very subtle energy-entropy balances, such as those that stabilise the blue phases, come into play. Beautiful early calculations, using approximate theories based on the Landau-de Gennes formalism that are analytically tractable, helped understand many of the features of the blue phases. Here we have shown that it is possible to make further progress by exploiting modern computational resources to minimise the free energy exactly.

We have shown that significant quantitative differences in the phase diagram arise from retaining cubic order terms in the free energy expansion. In particular, choosing expansion coefficients appropriate to systems where the cholesteric undergoes helical sense inversion gives rise to a very significant increase in the stability of blue phase I. It would be interesting to see if this link can be established experimentally. Moreover, by minimising the free energy with respect to both the order parameter field and the size and shape of the unit cell it is possible to predict the electrostriction of the blue phases, obtaining good qualitative and quantitative agreement with experiments.

The major drawback of the Landau-de Gennes expansion is in the difficulty of knowing values of the expansion coefficients for any particular compound. Inevitably this programme becomes increasingly difficult as additional terms are added. Therefore the value of the theory is primarily in suggesting trends, and in identifying regions of parameter space where novel behaviour might be observed.

Although the major features of the equilibrium behaviour of the blue phases are understood, much less is known about their hydrodynamics. The current simulations are on tiny samples and it is necessary to guess suitable boundary conditions on the director field. However, the reponse of the blue phases and their disclination networks to an imposed flow is fascinating, and a programme of developing the numerical approach in tandem with experiment is likely to uncover novel physics.

Solving the equations of motion of the blue phases also makes it possible to investigate their dynamics under changes in an applied electric field. Understanding this is an important step to assisting device design. Moreover it is possible to follow the kinetics of transitions between the blue phases, as an external parameter such as the electric field is varied. This allows prediction of possible sequences of phases that will, in general, depend on the values of the free energy, whether there is a convenient path through phase space allowing any given phase to be accessed, and the speed with which the perturbation is applied. For the examples described here the structural evolution depended primarily on relaxation kinetics, with hydrodynamics playing a minor role. However, this may not always be the case.

Finally we return to blue phases that have been stabilised over a larger temperature range. One approach has been to use bimesogenic molecules, and there remain interesting questions as to whether the resulting blue phase is stable or metastable and in identifying the physical mechanisms behind the increased temperature range. The polymer stabilised blue phases can be interpreted in terms of the polymers pinning the disclinations. As numerical work on ever more complex fluids develops, there are many exciting questions about blue phases in mixtures of liquid crystals and polymers, colloids and nanoparticles, that will become increasingly accessible.

###### Acknowledgements.

The authors would like to thank Davide Marenduzzo, Alexandre Dupuis and Enzo Orlandini for their help, advice and contributions to the research presented here. GPA would also like to thank Randall Kamien for discussions and acknowledges partial support from NSF Grant DMR05-47320.## References

- (1) de Gennes, P.G.; Prost, J. The Physics of Liquid Crystals; Clarendon: Oxford, 1993.
- (2) Reinitzer, F. Monatsh. Chem. 1888, 9, 421.
- (3) Blümel, T.; Collings, P.J.; Onusseit, H.; Stegemeyer, H. Chem. Phys. Lett. 1985, 116, 529.
- (4) Yang, D.K.; Crooker, P.P. Phys. Rev. A 1987, 35, 4419.
- (5) Crooker, P.P. Liq. Cryst. 1989, 5, 751.
- (6) Belyakov, V.A.; Dmitrienko, V.E. Usp. Fiz. Nauk 1985, 146, 369 [Sov. Phys. Usp. 1985, 28, 535].
- (7) Wright, D.C.; Mermin, N.D. Rev. Mod. Phys 1989, 61, 385.
- (8) Pieranski, P.; Barbet-Massin, R.; Cladis, P.E. Phys. Rev. A 1985, 31, 3912.
- (9) Meiboom, S.; Sethna, J.P.; Anderson, P.W.; Brinkman, W.F. Phys. Rev. Lett. 1981, 46, 1216.
- (10) Meiboom, S.; Sammon, M.; Brinkman, W.F. Phys. Rev. A 1983, 27, 438.
- (11) Hornreich, R.M.; Shtrikman, S. Phys. Rev. A 1981, 24, 635.
- (12) Hornreich, R.M.; Kugler, M.; Shtrikman, S. Phys. Rev. Lett. 1982, 48, 1404.
- (13) Lubensky, T.C.; Stark, H. Phys. Rev. E 1996, 53, 714.
- (14) Anisimov, M.A.; Agayan, V.A.; Collings, P.J. Phys. Rev. E 1998, 57, 582.
- (15) Cao, W.; Muñoz, A.; Palffy-Muhoray, P.; Taheri, B. Nat. Mater. 2002, 1, 111.
- (16) Hisakado, Y.; Kikuchi, H.; Nagamura, T.; Kajiyama, T. Adv. Mater. 2005, 17, 96.
- (17) Yokoyama, S.; Mashiko, S.; Kikuchi, H.; Uchida, K.; Nagamura, T. Adv. Mater. 2006, 18, 48.
- (18) Kitzerow, H.-S. ChemPhysChem 2006, 7, 63.
- (19) Coles, H.J.; Pivnenko, M.N. Nature 2005, 436, 997.
- (20) Kikuchi, H.; Yokota, M.; Hisakado, Y.; Yang, H.; Kajiyama, T. Nat. Mater. 2002, 1, 64.
- (21) http://en.wikipedia.org/wiki/Blue_Phase_Mode_LCD, retrieved 11 June 2008.
- (22) Kikuchi, H.; Haseba, Y.; Iwata, T.; Higuchi, H. SID07 Digest 2007, 38, 1737.
- (23) Grebel, H.; Hornreich, R.M.; Shtrikman, S. Phys. Rev. A 1983, 28, 1114.
- (24) Grebel, H.; Hornreich, R.M.; Shtrikman, S. Phys. Rev. A 1984, 30, 3264.
- (25) Dupuis, A.; Marenduzzo, D.; Yeomans, J.M. Phys. Rev. E 2005, 71, 011703.
- (26) Alexander, G.P.; Yeomans, J.M. Phys. Rev. E 2006, 74, 061706.
- (27) Dupuis, A.; Marenduzzo, D.; Orlandini, E.; Yeomans, J.M. Phys. Rev. Lett. 2005, 95, 097801.
- (28) Alexander, G.P.; Marenduzzo, D. Europhys. Lett. 2008, 81, 66004.
- (29) Stegemeyer, H.; Porsch, F. Phys. Rev. A 1984, 30, 3369.
- (30) Porsch, F.; Stegemeyer, F. Chem. Phys. Lett. 1986, 125, 319.
- (31) Porsch, F.; Stegemeyer, F. Chem. Phys. Lett. 1989, 155, 620.
- (32) Pieranski, P.; Cladis, P.E.; Barbet-Massin, R. J. Physique Lett. 1985, 46, L-973.
- (33) Cladis, P.E.; Garel, T.; Pieranski, P. Phys. Rev. Lett. 1986, 57, 2841.
- (34) Pieranski, P.; Cladis, P.E. Phys. Rev. A 1987, 35, 355.
- (35) Chen, N.; Ho, J. Phys. Rev. A 1987, 35, 4886.
- (36) Brazovskii, S.A.; Dmitriev, S.G. Zh. Eksp. Teor. Fiz. 1975, 69, 979 [Sov. Phys. JETP 1975, 42, 497].
- (37) Berreman D.W.; Meiboom, S. Phys. Rev. A 1984, 30, 1955.
- (38) Longa, L.; Monselesan, D.; Trebin, H.-R. Liq. Cryst. 1987, 2, 769.
- (39) Longa, L.; Trebin, H.-R. Phys. Rev. A 1989, 39, 2160.
- (40) Denniston, C.; Marenduzzo, D.; Orlandini, E.; Yeomans, J.M. Philos. Trans. R. Soc. London, Ser. A 2004, 362, 1745.
- (41) Yeomans, J.M. Physica A 2006, 369, 159.
- (42) Huff, B.P.; Krich, J.J.; Collings, P.J. Phys. Rev. E 2000, 61, 5372.
- (43) Slaney, A.J.; Nishiyama, I.; Styring, P.; Goodby, J.W. J. Mater. Chem 1992, 2, 805.
- (44) Kitzerow, H.-S. Mol. Cryst. Liq. Cryst. 1991, 202, 51.
- (45) Hornreich, R.M.; Kugler, M.; Shtrikman, S. Phys. Rev. Lett. 1985, 54, 2099.
- (46) Lubin, D.; Hornreich, R.M.; Phys. Rev. A 1987, 36, 849.
- (47) Stark, H.; Trebin, H.-R. Phys. Rev. A 1991, 44, 2752.
- (48) Hornreich, R.M.; Shtrikman, S. Phys. Rev. A 1990, 41, 1978.
- (49) Longa, L.; Zelazna, M.; Trebin, H.-R.; Mościcki, J. Phys. Rev. E 1996, 53, 6067.
- (50) Zelazna, M.; Longa, L.; Trebin, H.-R.; Stark, H. Phys. Rev. E 1998, 57, 6711.
- (51) Tóth, G.; Denniston, C.; Yeomans, J.M. Phys. Rev. Lett. 2002, 88, 105504.
- (52) Tóth, G.; Denniston, C.; Yeomans, J.M. Phys. Rev. E, 2003, 67, 051705.
- (53) Svensek, D.; Zumer, S. Phys. Rev. E 2002, 66, 021712.
- (54) Svensek, D.; Zumer, S. Phys. Rev. Lett. 2003, 90, 155501.
- (55) Denniston, C.; Orlandini, E.; Yeomans, J.M. Europhys. Lett. 2000, 52, 481.
- (56) Denniston, C.; Orlandini, E.; Yeomans, J.M. Phys. Rev. E 2001, 64, 021701.
- (57) Sulaiman, N.; Marenduzzo, D.; Yeomans, J.M. Phys. Rev. E 2006, 74, 041708.
- (58) Fukuda, J.I. Eur. Phys. J. B 1998, 1, 173.
- (59) Marenduzzo, D.; Orlandini, E.; Yeomans, J.M. Phys. Rev. Lett. 2004, 92, 188301.
- (60) Marenduzzo, D.; Orlandini, E.; Yeomans, J.M. J. Chem. Phys. 2004, 121, 582.
- (61) Orlandini, E.; Marenduzzo, D.; Yeomans, J.M. Comp. Phys. Commun. 2005, 169, 122.
- (62) Marenduzzo, D.; Dupuis, A.; Yeomans, J.M.; Orlandini, E. Mol. Cryst. Liq. Cryst. 2005, 435, 845.
- (63) Marenduzzo, D.; Orlandini, E.; Yeomans, J.M. J. Chem. Phys. 2006, 124, 204906.
- (64) Orlandini, E.; Marenduzzo, D.; Yeomans, J.M. Mol. Cryst. Liq. Cryst. 2007, 465, 1.
- (65) Qian, T.Z.; Xie, Z.L.; Kwok, H.S.; Sheng, P. Appl. Phys. Lett. 1997, 71, 596.
- (66) Qian, T.Z.; Xie, Z.L.; Kwok, H.S.; Sheng, P. J. Appl. Phys. 2001, 90, 3121.
- (67) Denniston, C.; Yeomans, J.M. Phys. Rev. Lett. 2001, 87, 275505.
- (68) Tóth, G.; Denniston, C.; Yeomans, J.M. Comp. Phys. Commun. 2002, 147, 7.
- (69) Denniston, C.; Tóth, G.; Yeomans, J.M. J. Stat. Phys. 2002, 107, 187.
- (70) Marenduzzo, D.; Orlandini, E.; Yeomans, J.M. Europhys. Lett. 2003, 64, 406.
- (71) Marenduzzo, D.; Orlandini, E.; Yeomans, J.M. Europhys. Lett. 2005, 71, 604.
- (72) Stromer, J.F.; Marenduzzo, D.; Brown, C.V.; Yeomans, J.M.; Raynes, E.P. J. Appl. Phys. 2006, 99, 064911.
- (73) James, R.; Willman, E.; Fernandez, A.; Day, S.E. IEEE Transactions on Magnetics 2008, 44, 814.
- (74) Beris, A.N.; Edwards, B.J. Thermodynamics of Flowing Systems; Oxford University Press: Oxford, 1994.
- (75) Grmela, M.; Öttinger, H.C. Phys. Rev. E 1997, 56, 6620.
- (76) Öttinger, H.C.; Grmela, M. Phys. Rev. E 1997, 56, 6633.
- (77) Öttinger, H.C. Beyond Equilibrium Thermodynamics; WileyBlackwell: 2005.
- (78) Qian, T.Z.; Sheng, P. Phys. Rev. E 1998, 58, 7475.
- (79) Denniston, C.; Orlandini, E.; Yeomans, J.M. Phys. Rev. E 2001, 63, 056702.
- (80) James, R.; Willman, E.; Fernandez, A.; Day, S.E. IEEE Transactions on Electron Devices 2006, 53, 1575.
- (81) Marenduzzo, D.; Orlandini, E.; Cates, M.E.; Yeomans, J.M. Phys. Rev. E 2007, 76, 031921.
- (82) Henrich, O.; Marenduzzo, D.; Stratford, K.; Cates, M.E. arXiv:0901.3293v1 [cond-mat.soft] 2009.
- (83) Kikuchi, H.; Hirata, S.; Uchida, K. Mol. Cryst. Liq. Cryst., 2007, 465, 283.
- (84) Higashiguchi, K.; Yasui, K.; Kikuchi, H. J. Am. Chem. Soc. 2008, 130, 6326.
- (85) Meyer, R.B. Phys. Rev. Lett. 1969, 22, 918.
- (86) Alexander, G.P.; Yeomans, J.M. Phys. Rev. Lett. 2007, 99, 067801.