Hydrodynamics and Electrohydrodynamics of Liquid Crystals

Particle-based and continuum models for confined nematics in two dimensions

We use the particle-based stochastic multi-particle collision dynamics (N-MPCD) algorithm to simulate confined nematic liquid crystals in regular two-dimensional polygons such as squares, pentagons and hexagons. We consider a range of values of the nematicities, U, and simulation domain sizes, R, that canvass nano-sized polygons to micron-sized polygons. We use closure arguments to define mappings between the N-MPCD parameters and the parameters in the continuum deterministic Landau-de Gennes framework. The averaged N-MPCD configurations agree with those predicted by Landau-de Gennes theory, at least for large polygons. We study relaxation dynamics or the non-equilibrium dynamics of confined nematics in polygons, in the N-MPCD framework, and the kinetic traps bear strong resemblance to the unstable saddle points in the Landau-de Gennes framework. Finally, we study nematic defect dynamics inside the polygons in the N-MPCD framework and the finite-size effects slow down the defects and attract them to polygon vertices. Our work is a comprehensive comparison between particle-based stochastic N-MPCD methods and deterministic/continuum Landau-de Gennes methods, and this comparison is essential for new-age multiscale theories.

Active nematic liquid crystals simulated by particle-based mesoscopic methods

Two Multi-particle collision dynamics algorithms that simulate nematic liquid crystals are generalised to reproduce active behaviour. One of the algorithms is due to Shendruk and Yeomans and is based on particles that carry an orientation vector ordered by a mean-field energy [T. N. Shendruk and J. M. Yeomans, Soft Matter, 2015, 11, 5101]. In the other algorithm, due to Mandal and Mazza, particles possess an order parameter tensor which evolves according to the Qian-Sheng model of nematohydrodynamics [S. Mandal and M. G. Mazza, Phys. Rev. E, 2019, 99, 063319]. For both methods activity is incorporated through a force proportional to the divergence of the local average order parameter tensor. Both implementations produce disclination curves in the nematic fluid that undergo nucleation and self-annihilation dynamics. Topological defects are found to be consistent with those observed in recent experiments of three-dimensional active nematics. Results permit to compare the length-scales over which the different nematic Multi-particle collision dynamics methods operate. The structure and dynamics of the orientation and flow fields agree with those obtained recently in numerical studies of continuum three-dimensional active nematics. Overall, our results open the opportunity to use mesoscopic particle-based approaches to study active liquid crystals in situations such as nonequilibrium states driven by flow or colloidal particles in active anisotropic solvents.

Hidden traces of chirality in the fluctuations of a fully unwound cholesteric

Confinement and hydrodynamic interactions often play an important role in the fluctuation dynamics of soft matter systems, which can typically be studied using light scattering techniques. With experimental and theoretical methodologies, I demonstrate here that chirality is an additional critical parameter that leads to diverging decay times and correlation lengths in chiral liquid crystal cells with a fully unwound cholesteric helix. This study combines light scattering measurements made in a tailored microscope geometry and theoretical calculations of the decay dynamics of chiral orientational fluctuations-including hydrodynamics-to establish the existence of two soft chiral modes of fluctuations driving the destabilization of the unwound cholesteric. Despite the achirality of the equilibrium state of unwound cholesterics, this study indicates that chirality hides itself in the orientational fluctuation modes and plays a major role in their dynamics, which can be exploited to locally measure the strength of chirality in frustrated chiral liquid crystal cells.

When low-order expansions fail and all higher-order contributions matter-basic example of the mean squared displacement for Brownian motion

Hardly any theoretically formulated realistic problem can be solved exactly. Therefore, as a standard, we resort to approximations. In this context, expansions play a major role. We are used to relying on lowest-order expansions and confining our point of view accordingly. However, one should always bear in mind that such considerations may fail at some point. Here, we address a very common example situation, namely, the motion of a Brownian particle. We know that the associated mean-squared displacement in the long term increases linearly in time. Yet, when we take the Fokker-Planck approach in combination with a low-order expansion, the direct route towards this result fails. That is, in the expansion the term linear in time vanishes. Instead, the treatment requires consideration of all higher-order contributions. Together, they restore the linear increase in time. In this way, we stress that care is always mandatory when resorting to low-order expansions, and we present in a traceable way a route to solving the considered problem.

Symmetry aspects in the macroscopic dynamics of magnetorheological gels and general liquid crystalline magnetic elastomers

We investigate theoretically the macroscopic dynamics of various types of ordered magnetic fluid, gel, and elastomeric phases. We take a symmetry point of view and emphasize its importance for a macroscopic description. The interactions and couplings among the relevant variables are based on their individual symmetry behavior, irrespective of the detailed nature of the microscopic interactions involved. Concerning the variables we discriminate between conserved variables related to a local conservation law, symmetry variables describing a (spontaneously) broken continuous symmetry (e.g., due to a preferred direction) and slowly relaxing ones that arise from special conditions of the system are considered. Among the relevant symmetries, we consider the behavior under spatial rotations (e.g., discriminating scalars, vectors or tensors), under spatial inversion (discriminating e.g., polar and axial vectors), and under time reversal symmetry (discriminating e.g., velocities from polarizations, or electric fields from magnetic ones). Those symmetries are crucial not only to find the possible cross-couplings correctly but also to get a description of the macroscopic dynamics that is compatible with thermodynamics. In particular, time reversal symmetry is decisive to get the second law of thermodynamics right. We discuss (conventional quadrupolar) nematic order, polar order, active polar order, as well as ferromagnetic order and tetrahedral (octupolar) order. In a second step, we show some of the consequences of the symmetry properties for the various systems that we have worked on within the SPP1681, including magnetic nematic (and cholesteric) elastomers, ferromagnetic nematics (also with tetrahedral order), ferromagnetic elastomers with tetrahedral order, gels and elastomers with polar or active polar order, and finally magnetorheological fluids and gels in a one- and two-fluid description.

Cubic and tetragonal liquid crystal phases composed of non-chiral molecules: Chirality and macroscopic properties

We discuss the symmetry properties as well as the macroscopic behavior of the cubic liquid crystal phases showing large chiral domains of either hand in some non-chiral compounds reported recently in the group of Tschierske. These phases are tricontinuous. While they have O or I432 symmetry in each chiral domain, the overall symmetry is [Formula: see text] as there is no net chirality for compounds composed of non-chiral molecules. It turns out that a rather similar type of phase has also been reported for triblock copolymers. Here we analyze in detail the macroscopic static and dynamic behavior of such phases and we predict, among other results, that they show the analog of static and dissipative Lehmann-type effects in their chiral domains. A description of a cubic liquid crystalline phase of T_h symmetry, which has not yet been found experimentally, is also included. Suggestions for experiments are outlined to identify such a phase. In addition, we discuss tetragonal liquid crystalline phases of D_4h and D₄ (I422) symmetry as they have been reported last year experimentally in connection to the Q phase.

Influence of tetrahedral order on ferromagnetic gel phases

We investigate the macroscopic dynamics of gels with tetrahedral/octupolar symmetry, which possess in addition a spontaneous permanent magnetization. We derive the corresponding static and dynamic macroscopic equations for a phase, where the magnetization is parallel to one of the improper fourfold tetrahedral symmetry axes. Apart from elastic strains, we take into account relative rotations between the magnetization and the elastic network. The influence of tetrahedral order on these degrees of freedom is investigated and some experiments are proposed that are specific for such a material and allow to indirectly detect tetrahedral order. We also consider the case of a transient network and predict that stationary elastic shear stresses arise when a temperature gradient is applied.

On the influence of a network on optically isotropic fluid phases with tetrahedral/octupolar order

We investigate the influence of transient or permanent elasticity on liquid phases with octupolar (tetrahedral) order, a question that has never been addressed before. The focus will be on optically isotropic liquid phases with tetrahedral order including the T _d phase and the chiral T phase introduced by Fel. It turns out that the presence of both, a network as well as tetrahedral order can lead to the formation of chiral domains of both hands in an optically isotropic fluid due to a completely novel mechanism, thus providing a possible macroscopic explanation for recent experimental observations. We study in detail how elasticity influences the macroscopic dynamics of both, the T _d and the T phase. The simultaneous presence of a transient network as well as of octupolar order is shown to lead to completely new cross-coupling phenomena for optically isotropic systems including transient dissipative elastic strains due to temperature gradients.

Thermophoretically induced large-scale deformations around microscopic heat centers

Selectively heating a microscopic colloidal particle embedded in a soft elastic matrix is a situation of high practical relevance. For instance, during hyperthermic cancer treatment, cell tissue surrounding heated magnetic colloidal particles is destroyed. Experiments on soft elastic polymeric matrices suggest a very long-ranged, non-decaying radial component of the thermophoretically induced displacement fields around the microscopic heat centers. We theoretically confirm this conjecture using a macroscopic hydrodynamic two-fluid description. Both thermophoretic and elastic effects are included in this theory. Indeed, we find that the elasticity of the environment can cause the experimentally observed large-scale radial displacements in the embedding matrix. Additional experiments confirm the central role of elasticity. Finally, a linearly decaying radial component of the displacement field in the experiments is attributed to the finite size of the experimental sample. Similar results are obtained from our theoretical analysis under modified boundary conditions.

Macroscopic behavior of polar nematic gels and elastomers

We present the derivation of the macroscopic equations for uniaxial polar nematic gels and elastomers. We include the strain field as well as relative rotations as independent dynamic macroscopic degrees of freedom. As a consequence, special emphasis is laid on possible static and dynamic cross-couplings between these macroscopic degrees of freedom associated with the network, and the other macroscopic degrees of freedom including reorientations of the macroscopic polarization. In particular, we find static and dissipative dynamic cross-couplings between strain fields and relative rotations on one hand and the macroscopic polarization on the other that allow for new possibilities to manipulate polar nematics. To give one example each for the effects of a static and a dissipative cross-coupling: we find that a static electric field applied perpendicularly to the polar preferred direction leads to relative rotations while dynamically relative rotations can lead to transverse electric currents. In addition to a permanent network, we also consider the effect of a transient network, which is particularly important for the case of gels, melts and concentrated polymer solutions. A section on the influence of macroscopic chirality is included as well.

Macroscopic behavior of ferrocholesteric liquid crystals and ferrocholesteric gels and elastomers

We study the influence of macroscopic chirality on the macroscopic properties of superparamagnetic liquid crystals and gels. Specifically we derive macroscopic dynamic equations for ferrocholesteric low molecular weight (LMW) liquid crystals and for ferrocholesteric gels and elastomers in the local description using the director field as macroscopic variable. The magnetization is treated as a macroscopic dynamic degree of freedom and its coupling to all other macroscopic variables is examined in detail. We incorporate into our dynamic analysis terms that are linear in a magnetic field giving rise to a number of cross-coupling terms not possible otherwise. A number of properties that are unique to the class of systems studied arise. As an example for a static property we find a term in the generalized energy which is linear in the electric field and quadratic in the magnetic field. We find that applying a magnetic field to a ferrocholesteric can lead to reversible electric currents, heat currents and concentration currents, which change their sign with a sign change of macroscopic chirality. As an example of a rather intriguing dissipative dynamic contribution we point out that for ferrocholesterics and for ferrocholesteric gels and elastomers in a magnetic field extensional flow leads to electric and heat currents.

Macroscopic behavior of ferronematic gels and elastomers

We present the derivation of the macroscopic equations for uniaxial ferronematic gels and elastomers. We deal with the superparamagnetic case, where no permanent magnetization is present and the anisotropy is provided by the nematic director. We include the magnetization as an independent dynamic degree of freedom. As a consequence special emphasis is laid on possible static and dynamic cross-couplings between magnetization and the non-magnetic degrees of freedom, as director reorientations, flow, elastic strains and relative rotations between director and the elastic network. In particular, we find reversible dynamic cross-couplings among rotations of the magnetization, the director, relative rotations, and deformational flow that allow for new possibilities to manipulate such materials. Application of simple (oscillatory) shear induces, in general, a finite magnetization normal to the shear plane and a relative rotation in the shear plane, whose amplitudes are linear in the shear rate. Induced magnetization, induced relative rotation and the director are mutually orthogonal, with the director aligned obliquely to the flow direction. This orientation is independent of the shear rate and is a material property.

Reversible and dissipative macroscopic contributions to the stress tensor: active or passive?

The issue of dynamic contributions to the macroscopic stress tensor has been of high interest in the field of bio-inspired active systems over the last few years. Of particular interest is a direct coupling ("active term") of the stress tensor with the order parameter, the latter describing orientational order induced by active processes. Here we analyze more generally possible reversible and irreversible dynamic contributions to the stress tensor for various passive and active macroscopic systems. This includes systems with tetrahedral/octupolar order, polar and non-polar (chiral) nematic and smectic liquid crystals, as well as active fluids with a dynamic preferred (polar or non-polar) direction. We show that it cannot a priori be seen, neither from the symmetry properties of the macroscopic variables involved, nor from the structure of the cross-coupling contributions to the stress tensor, whether the system studied is active or passive. Rather, that depends on whether the variables that give rise to those cross-couplings in the stress tensor are driven or not. We demonstrate that several simplified descriptions of active systems in the literature that neglect the necessary counter term to the active term violate linear irreversible thermodynamics and lead to an unphysical contribution to the entropy production.

Low symmetry tetrahedral nematic liquid crystal phases: Ambidextrous chirality and ambidextrous helicity

We discuss the symmetry properties as well as the dynamic behavior of various non-polar nematic liquid crystal phases with tetrahedral order. We concentrate on systems that show biaxial nematic order coexisting with octupolar (tetrahedral) order. Non-polar examples are phases with D2 and S4 symmetries, which can be characterized as biaxial nematics lacking inversion symmetry. It is this combination that allows for new features in the statics and dynamics of these phases. The D2-symmetric phase is chiral, even for achiral molecules, and shows ambidextrous chirality in all three preferred directions. The achiral S4-symmetric phase allows for ambidextrous helicity, similar to the higher-symmetric D2d-symmetric phase. Such phases are candidates for nematic phases made from banana-shaped molecules.

Active polar two-fluid macroscopic dynamics

We study the dynamics of systems with a polar dynamic preferred direction. Examples include the pattern-forming growth of bacteria as well as shoals of fish, flocks of birds and migrating insects. Due to the fact that the preferred direction only exists dynamically, but not statically, the macroscopic variable of choice is the macroscopic velocity associated with the motion of the active units, which are typically biological in nature. We derive the macroscopic equations for such a system and discuss novel static, reversible and irreversible cross-couplings connected to a second velocity as a variable. We analyze in detail how the macroscopic behavior of an active system with a polar dynamic preferred direction compares to other systems with two velocities including immiscible liquids and electrically neutral quantum liquids such as superfluid (4)He and (3)He . We critically discuss changes in the normal mode spectrum when comparing uncharged superfluids, immiscible liquids and active system with a polar dynamic preferred direction. We investigate the influence of a macroscopic hand (collective effects of chirality) on the macroscopic behavior of such active media.

Macroscopic behavior of non-polar tetrahedratic nematic liquid crystals

We discuss the symmetry properties and the macroscopic behavior of a nematic liquid crystal phase with D(2d) symmetry. Such a phase is a prime candidate for nematic phases made from banana-shaped molecules where the usual quadrupolar order coexists with octupolar (tetrahedratic) order. The resulting nematic phase is nonpolar. While this phase could resemble the classic D (infinityh) nematic in the polarizing microscope, it has many static as well as reversible and irreversible properties unknown to nonpolar nematics without octupolar order. In particular, there is a linear gradient term in the free energy that selects parity leading to ambidextrously helical ground states when the molecules are achiral. In addition, there are static and irreversible coupling terms of a type only met otherwise in macroscopically chiral liquid crystals, e.g. the ambidextrous analogues of Lehmann-type effects known from cholesteric liquid crystals. We also discuss the role of hydrodynamic rotations about the nematic director. For example, we show how strong external fields could alter the D(2d) symmetry, and describe the non-hydrodynamic aspects of the dynamics, if the two order structures, the nematic and the tetrahedratic one, rotate relative to each other. Finally, we discuss certain nonlinear aspects of the dynamics related to the non-commutativity of three-dimensional finite rotations as well as other structural nonlinear hydrodynamic effects.

Nonisothermal model for the direct isotropic/smectic-A liquid-crystalline transition

An extension to a high-order model for the direct isotropic/smectic-A liquid-crystalline phase transition was derived to take into account thermal effects including anisotropic thermal diffusion and latent heat of phase ordering. Multiscale multitransport simulations of the nonisothermal model were compared to isothermal simulation, showing that the presented model extension corrects the standard Landau-de Gennes prediction from constant growth to diffusion-limited growth under shallow quench/undercooling conditions. Nonisothermal simulations, where metastable nematic preordering precedes smectic-A growth, were also conducted, and novel nonmonotonic phase-transformation kinetics were observed.

Reversible macroscopic dynamics of polar nematic liquid crystals: reversible currents and their experimental consequences

Polar liquid crystalline phases are relevant for fluid liquid crystal phases observed in banana liquid crystals as well as for a class of polymeric liquid crystalline materials investigated recently. In this Brief Report we present the reversible dynamics for polar nematic liquid crystals with C_{infinityv} symmetry, where the macroscopic polarization representing polar order acts as an independent macroscopic variable. We find reversible coupling terms, for example, between flow and temperature and concentration gradients specific for the existence of a polar preferred direction. We suggest concrete experiments to check the importance of the reversible dynamic cross-coupling terms presented here.

Shear-induced undulation of smectic-A: Molecular Dynamics simulations vs. analytical theory

Experiments on a variety of systems have shown that layered liquids are unstable under shear even if the liquid layers are planes of constant velocity. We investigate the stability of smectic- A like liquids under shear using Molecular Dynamics simulations and a macroscopic hydrodynamic theory (including the layer normal and the director as independent variables). Both methods show an instability of the layers, which sets in above a critical shear rate. We find a remarkable qualitative and reasonable quantitative agreement between both methods for the spatial homogeneous state and the onset of the instability.

Hydrodynamic instabilities in ferronematics

In the hydrodynamic description of ferronematics there are various dynamic magnetic-field effects, linear in the field strength, that are negligible in usual nematics, but can play a role in ferronematics. Here, we investigate theoretically the influence of these new terms on the thermal convection (Rayleigh-Bénard) and the viscous fingering (Saffman-Taylor) instability in ferronematics in the presence of a strong magnetic field. We find that the instabilities are qualitatively changed due to the occurrence of a finite vorticity component--a feature that is known from simple liquids in the case of a superimposed mechanical rotation. We suggest to use the additional effects (cross-flow within convection rolls, oblique rolls, rotating fingers) for measuring the phenomenological coefficients involved.