Leslie Coefficients

Liquid crystals represent a fascinating state of matter, paradoxically combining the fluidity of a liquid with the long-range orientational order of a solid. This unique duality means that describing their flow, or hydrodynamics, is far more complex than for a simple fluid like water. Where a single viscosity number suffices for water, the internal friction of a liquid crystal depends intricately on the direction of flow relative to its aligned, rod-like molecules. This creates a rich and complex behavior that cannot be captured by classical fluid dynamics.

This article delves into the cornerstone of nematic liquid crystal hydrodynamics: the Leslie-Ericksen theory and its fundamental parameters, the Leslie viscosity coefficients. This framework provides the essential language to understand the coupled dance between fluid motion and molecular orientation. We will bridge the gap between abstract theory and tangible phenomena, showing how these coefficients are not just mathematical constructs but the architects of material behavior.

The article is structured to guide you from fundamental principles to real-world impact. The coming chapters, "Principles and Mechanisms" and "Applications and Interdisciplinary Connections," will explore this topic in depth. We will first define the Leslie coefficients, explaining how they contribute to viscous stress and torque, and reveal the critical distinction between flow-aligning and tumbling nematics. We will then explore the profound consequences of this theory, from the switching speed of displays to the formation of complex flow patterns and the microscopic origins of viscosity itself.

Imagine stirring honey. It's thick, it's viscous, and it resists your motion. Now, imagine stirring a pot full of uncooked spaghetti. It's a completely different story, isn't it? The resistance you feel depends dramatically on whether you stir along the length of the spaghetti strands or try to drag them sideways. This is the world of nematic liquid crystals. They are fluids, but they are fluids with a direction, a grain, much like our pot of spaghetti. To understand how they flow, we can't just use one number for viscosity like we do for water or honey. We need a richer, more descriptive language. This is the language of the Leslie-Ericksen theory, and its vocabulary is a set of six remarkable numbers known as the ​​Leslie viscosity coefficients​​.

At the heart of liquid crystal hydrodynamics lies a beautiful interplay between two main characters. The first is the familiar ​​velocity field​​, $\vec{v}$, which tells us how the fluid is moving at every point. The second is a new character: the ​​director field​​, $\vec{n}$, a unit vector that points along the average direction of the rod-like molecules at each point. The story of how a liquid crystal flows is the story of how $\vec{v}$ and $\vec{n}$ influence each other—a coupled dance where the flow tries to align the directors, and the directors, in turn, affect the flow's resistance.

To speak about this dance precisely, we first need to break down the fluid's motion. Any motion of a small fluid element can be thought of as a combination of stretching and rotating. We capture this with two mathematical tools. The ​​rate-of-strain tensor​​, $A_{ij} = \frac{1}{2}(\partial_j v_i + \partial_i v_j)$, describes how the fluid is being stretched or compressed. Imagine a tiny circle drawn in the fluid; $A_{ij}$ tells you how it deforms into an ellipse. The ​​vorticity tensor​​, $W_{ij} = \frac{1}{2}(\partial_j v_i - \partial_i v_j)$, describes how the fluid is locally swirling or rotating, like a tiny whirlpool.

Now, what about the director? It's carried along by the fluid, but it can also rotate relative to the local fluid swirl. This "rebellious" rotation, the motion of the director as seen by an observer tumbling along with the local fluid element, is fundamentally important. We give it a special name: the ​​co-rotational time derivative of the director​​, $\vec{N}$. In mathematical terms, $N_i = \dot{n}_i - W_{ik}n_k$, where $\dot{n}_i$ is the rate of change of the director as it's carried by the flow. It is this relative motion, this struggle between the director and the flow's vorticity, that generates much of the unique physics of liquid crystals.

In an ordinary (isotropic) fluid, the internal friction, or ​​viscous stress​​, is simply proportional to the rate of strain. But for our "fluid of rods," the situation is far more intricate. The stress depends not only on the strain but also on the director's orientation and its motion relative to the flow. The Leslie-Ericksen theory provides the complete "recipe" for this anisotropic friction. The viscous stress, $\sigma'_{ij}$, is a linear combination of all the physically allowed terms involving our key ingredients: the director $\vec{n}$, the strain rate $\mathbf{A}$, and the co-rotational rate $\vec{N}$.

The full recipe looks a bit intimidating at first:

The six coefficients, $\alpha_1$ through $\alpha_6$, are the Leslie viscosity coefficients. They are the personality traits of the liquid crystal, defining its unique rheological character. Let's not be scared by the equation; let's see what it tells us.

• The ​​$\alpha_4$ term​​ is the most familiar. It's just like the viscosity of a normal fluid, where stress is directly proportional to strain. If the liquid crystal had no orientational order, this would be the only term that mattered.

• The ​​$\alpha_5$ and $\alpha_6$ terms​​ describe a coupling between the fluid's stretching ($A_{ij}$) and the director's orientation ($\vec{n}$). They represent the extra stress generated when the fluid tries to deform in a way that is misaligned with the director.

• The ​​$\alpha_2$ and $\alpha_3$ terms​​ are fascinating. They connect the stress directly to $\vec{N}$, the director's rotation relative to the fluid. This means that just forcing the director to spin against the local fluid vortex creates internal friction, even if there is no stretching!

• Finally, the ​​$\alpha_1$ term​​ is purely anisotropic. The term $(n_k n_p A_{kp})$ measures how much the fluid is being stretched along the direction of the director. The $\alpha_1$ coefficient then determines how much stress this specific type of stretching produces.

Together, these six coefficients orchestrate the complex viscous response of the fluid.

A shear flow doesn't just create stress; it also exerts a ​​viscous torque​​ on the director, trying to twist it into a new orientation. Think of a log in a river; the current will try to turn it. This torque has two competing parts. There's a part that comes from the fluid's stretching ($A_{ij}$), and a part that comes from its rotation ($W_{ij}$). A stable orientation is reached only when the total viscous torque is zero.

TheLeslie-Ericksen theory beautifully captures this balance. The viscous torque is governed by two special combinations of the Leslie coefficients:

• ​​$\gamma_1 = \alpha_3 - \alpha_2$​​, known as the ​​rotational viscosity​​. It acts like pure rotational friction. If you try to rotate the director, $\gamma_1$ determines how strongly it resists that rotation. It must always be positive for the liquid crystal to be stable.

• ​​$\gamma_2 = \alpha_6 - \alpha_5$​​, an equally important coefficient that doesn't have a simple name but has a profound job. It describes the torque exerted by the stretching part of the flow.

The fate of a director in a shear flow is decided by the competition between the torque from vorticity (related to $\gamma_1$) and the torque from strain (related to $\gamma_2$). When these torques balance, the director settles at a specific, stable angle to the flow, known as the ​​Leslie angle​​, $\theta_L$. The condition for this balance is remarkably simple: $\cos(2\theta_L) = -\gamma_1 / \gamma_2$.

This little equation holds a deep secret. If the ratio $|\gamma_1 / \gamma_2|$ is less than or equal to 1, a stable angle exists, and the liquid crystal is called ​​flow-aligning​​. The directors settle into a peaceful, streamlined orientation with the flow. If, however, $|\gamma_1 / \gamma_2|$ is greater than 1, there is no angle where the torques can balance. The rotational torque always wins. The directors can never find peace; they are doomed to continuously tumble and spin forever as the fluid flows. This tumbling behavior is just as real and important as flow alignment, and it's all decided by the values of those few $\alpha$ coefficients.

The $\alpha$ coefficients might seem abstract, but they have direct, measurable consequences. The classic experiment that makes this connection is the ​​Miesowicz experiment​​. The idea is simple: use a strong magnetic field to force the director to point in a specific direction, then apply a simple shear flow (like flow between two moving plates) and measure the effective viscosity $\eta$. By choosing three special orientations, we can isolate different combinations of the Leslie coefficients.

Let's say the flow is in the $x$-direction, with the velocity gradient in the $y$-direction.

1. ​​Director along the vorticity ($z$-axis, $\vec{n} = (0,0,1)$):​​ The director is perpendicular to both the flow and the gradient, like a log standing upright in a shallow river. In this orientation, the flow barely interacts with the director's length. The resulting viscosity, called ​​$\eta_c$​​, is the simplest of all: $\eta_c = \frac{1}{2}\alpha_4$. It reflects that "normal fluid" part of the friction.

2. ​​Director along the flow ($x$-axis, $\vec{n} = (1,0,0)$):​​ The directors are streamlined with the flow. The resistance should be low. The viscosity, ​​$\eta_a$​​, is found to be $\eta_a = \frac{1}{2}(\alpha_3 + \alpha_4 + \alpha_6)$.

3. ​​Director along the gradient ($y$-axis, $\vec{n} = (0,1,0)$):​​ The directors are broadside to the flow, "plowing" through the velocity layers. This should be the orientation of highest resistance. The viscosity, ​​$\eta_b$​​, is $\eta_b = \frac{1}{2}(-\alpha_2 + \alpha_4 + \alpha_5)$.

Typically, for rod-like nematics, we find $\eta_b > \eta_c > \eta_a$, which matches our intuition perfectly. These three measurements provide us with three equations that help us determine the values of the underlying Leslie coefficients. If we let the director sit at any other angle $\theta$ in the shear plane, we find that the effective viscosity $\eta_{eff}(\theta)$ varies continuously between these extremes, a direct manifestation of the fluid's anisotropy. When the director is free to choose its own orientation, it settles at the Leslie angle $\theta_L$, resulting in an effective viscosity that is a specific, complex combination of the Miesowicz viscosities and $\alpha_1$. The system naturally finds a state of motion, and the energy it dissipates as heat is a direct consequence of this self-organized state.

This leaves us with a question. We started with six coefficients. Are they all truly independent? It seems like a lot. Physics often seeks simplicity and underlying unity. It turns out there is a hidden relationship lurking among them, one that comes not from mechanics but from a deeper level of physics: thermodynamics.

The theory of irreversible processes, pioneered by Lars Onsager, states that for any system near thermal equilibrium, there is a fundamental symmetry that relates the coupling between different dissipative processes. When we apply this principle to the entropy production in a liquid crystal, considering the coupling between fluid strain ($A_{ij}$) and director rotation ($N_i$), a beautiful constraint emerges. This is the famous ​​Parodi relation​​:

This is remarkable! It's the same combination of coefficients we saw in our torque balance: $\gamma_2 = \alpha_6 - \alpha_5$. The Parodi relation tells us that $\gamma_2 = \alpha_2 + \alpha_3$. This isn't just a coincidence; it's a profound statement about the thermodynamic consistency of the theory. It reduces the number of independent Leslie coefficients from six to five, revealing a hidden simplicity.

But can we go deeper? Where do these coefficients come from in the first place? They depend on temperature and pressure, but fundamentally they must arise from the collective behavior of the billions of molecules that make up the fluid. More microscopic theories, like the Beris-Edwards model, describe the liquid crystal using an order parameter tensor that captures not only the direction but also the degree of molecular alignment. By comparing the stress tensor from such a model with the Leslie-Ericksen form, we can indeed derive expressions for the $\alpha$ coefficients in terms of more fundamental quantities like the scalar order parameter $S$ and microscopic interaction parameters. This completes our journey, connecting the macroscopic flow we can see and measure all the way down to the world of the molecules themselves, all described by this elegant set of five (not six!) fundamental coefficients.

We have seen that the liquid crystal is a strange and beautiful state of matter, a fluid that flows yet possesses a direction. We have learned that its response to being pushed and stirred is not as simple as that of water or honey. This complexity is captured by a handful of numbers, the Leslie viscosity coefficients. At first glance, these coefficients, the six $\alpha_i$, might seem like a mere set of phenomenological parameters, a dry catalog of material properties. But to think that would be to miss the music for the notes. These coefficients are the choreographers of a magnificent and intricate dance that plays out across technology, biology, and the fundamental physics of matter. They decide whether a material will gracefully align with a flow or stubbornly tumble within it; they dictate the switching speed of the screen on which you might be reading this; and they even emerge from the deepest, most random jiggles of the molecules themselves. Let us now explore this world of applications, to see the profound and often surprising unity that the Leslie-Ericksen theory reveals.

Imagine you throw a long, thin stick into the middle of a river. Near the bank, the water flows slower than in the middle. This difference in speed, this shear, will both drag the stick downstream and try to turn it. The stick feels two competing urges: the swirling part of the flow (the vorticity) wants to spin it around like a log-roller, while the stretching part of the flow (the strain) wants to grab its ends and align it with the flow. Which urge wins?

For a simple stick in a river, the answer is complex. But for the rod-like molecules in a nematic liquid crystal, the answer is beautifully simple and is decided entirely by the Leslie coefficients. A nematic material subjected to a shear flow has a choice: it can be ​​flow-aligning​​, where its director finds a stable, quiet angle with respect to the flow, or it can be a ​​tumbler​​, where the director is doomed to endlessly rotate, forever caught between the conflicting demands of vorticity and strain.

This entire drama is governed by a single dimensionless quantity, the flow-alignment parameter, typically denoted by $\lambda$. This parameter is a specific ratio of Leslie coefficients, most commonly $\lambda = (\alpha_6 - \alpha_5)/(\alpha_3 - \alpha_2)$. The director finds peace and aligns with the flow only if the aligning influence of strain is strong enough to overcome the rotation. This happens when the magnitude of the alignment parameter is greater than or equal to one: $|\lambda| \ge 1$. If this condition is met, the director settles at a specific "Leslie angle" relative to the flow, an angle whose value is simply given by $\cos(2\theta) = -1/\lambda$. If $|\lambda| \lt 1$, no such steady angle is possible, and the director must tumble.

For many common rod-like nematics, the coefficient $\alpha_2$ is negative, and the rotational viscosity $\gamma_1 = \alpha_3 - \alpha_2$ is positive. In this common scenario, the great divide between aligning and tumbling often rests on the sign of a single coefficient, $\alpha_3$. If $\alpha_3$ is negative, the material typically flow-aligns. If $\alpha_3$ is positive, it tumbles. The boundary case, $\alpha_3 = 0$, marks the critical transition between these two fundamental behaviors. This simple criterion is not just an academic curiosity; it is a primary design principle for formulating liquid crystal mixtures for applications ranging from displays to rheological modifiers.

Why do some materials tumble while others align? What is the microscopic origin of the magical number $\lambda$? The beauty of physics is that we can connect these macroscopic behaviors to the world of molecules. The Leslie coefficients are not arbitrary; they are determined by the very nature of the particles and their collective behavior.

The most important factors are particle ​​shape​​ and the degree of ​​orientational order​​. It turns out that prolate, or rod-like, particles (think microscopic cigars) naturally lead to a positive flow-alignment parameter ($\lambda \gt 0$), while oblate, or disc-like, particles (microscopic frisbees) give a negative one ($\lambda \lt 0$). Furthermore, as the liquid crystal becomes more ordered—that is, as the molecules become more parallel to each other, a quantity measured by the order parameter $S$—the anisotropic effects become stronger, and the magnitude of $\lambda$ increases. A more ordered nematic is more "aware" of its directionality and thus couples more strongly to the flow's anisotropy.

This provides a wonderful contrast between different types of liquid crystals. Consider a thermotropic nematic, made of small organic molecules that align upon cooling, which often has a high order parameter ($S \approx 0.65$). Compare it to a lyotropic nematic, formed by dissolving long, rod-like polymers in a solvent, which might be only weakly ordered just above the concentration needed to form the nematic phase ($S \approx 0.5$). The small-molecule thermotropic, with its higher order and simpler microscopic dynamics, often has $|\lambda| \gt 1$ and readily flow-aligns. In contrast, the lyotropic polymer system is a far more complex beast. The long, floppy chains create "viscoelastic" effects; the stress in the fluid is coupled to the orientation of the polymer network. This coupling effectively resists alignment, reducing the value of $|\lambda|$ and making tumbling a much more common behavior. This distinction is crucial for processing polymer solutions, from spinning high-strength fibers to fabricating optical films.

So far, we have discussed how flow orients the director. But the coupling is a two-way street: the director's motion also creates flow. This phenomenon, known as ​​backflow​​, is a subtle but profoundly important consequence of the Leslie-Ericksen theory.

Imagine trying to quickly turn a large paddle in a tub of honey. As you turn the paddle, the honey is dragged along with it, creating a flow that resists the paddle's motion. The same thing happens in a liquid crystal display (LCD). Each pixel works by reorienting the director field with an electric field. As the directors turn, they drag the fluid along. This induced backflow creates a viscous torque that opposes the reorientation, effectively increasing the rotational viscosity. This "effective rotational viscosity," $\gamma_{1, \text{eff}}$, is always larger than the intrinsic rotational viscosity $\gamma_1$. The size of this backflow correction depends on a specific combination of Leslie coefficients. It is this effective viscosity that ultimately determines how fast a pixel can switch from dark to light, and thus governs the response time of our screens and projectors. A materials scientist designing a fast-switching LCD must therefore select materials with Leslie coefficients that minimize this backflow effect.

This "unseen hand" of backflow also choreographs the dance of topological defects. Liquid crystals are rarely perfectly ordered; they are often threaded by "disclination lines," which are like seams or scars in the orientational fabric. When a pair of disclinations with opposite topological "charge" (say, a $+1/2$ and a $-1/2$ defect) find each other, they are pulled together by the elastic energy stored in the surrounding distortion. As they move to annihilate, they must push the fluid out of their way. Their motion is a delicate balance between the elastic driving force and the hydrodynamic drag from the fluid they displace. This drag force is not simple; it is mediated by the complex anisotropic viscosity described by the Leslie coefficients. Understanding this dynamic is key to controlling the texture and annealing of liquid crystal materials during device fabrication.

The coupling between flow and orientation, mediated by the Leslie coefficients, can do more than just modify existing flows—it can create entirely new and unexpected phenomena.

A stunning example is ​​electrohydrodynamic convection (EHD)​​. When you apply a voltage across a thin layer of a nematic liquid crystal with the right properties (negative dielectric anisotropy and positive conductivity anisotropy), a remarkable thing happens. Above a critical voltage, the initially placid fluid erupts into a regular pattern of rolling convective cells, visible under a microscope as a beautiful array of stripes known as "Williams domains." This pattern arises from a delicate feedback loop: the electric field drives charge separation, which causes ions to flow. This ion flow drags the fluid, creating a shear flow. The shear flow then exerts a viscous torque on the directors, causing them to tilt. This tilt, in turn alters the local conductivity and dielectric permittivity, reinforcing the initial charge separation. The crucial link in this chain is the viscous torque that tilts the director. This destabilizing torque only exists if the Leslie coefficient $\alpha_2$ is negative. If you use a material where $\alpha_2$ is positive, the torque acts in the opposite, stabilizing sense, the feedback loop is broken, and the beautiful patterns never form. The existence of this entire class of instabilities is thus predicated on the signs of the Leslie coefficients.

The influence of Leslie coefficients also extends to flows in more complex geometries. Consider a fluid flowing through a curved pipe. Even for a simple Newtonian fluid, the centrifugal force induces a secondary flow—a pair of counter-rotating vortices in the pipe's cross-section known as Dean vortices. In a nematic liquid crystal, this effect becomes far richer. If the nematic is flow-aligned, its directors point primarily along the pipe. The primary flow along the pipe feels one viscosity, while the weak secondary flow in the cross-section, being perpendicular to the directors, feels a different viscosity. These viscosities, known as Miesowicz viscosities, are just different linear combinations of the $\alpha_i$. For instance, the viscosity for flow parallel to the director is $\eta_a = \frac{1}{2}(\alpha_3 + \alpha_4 + \alpha_6)$, while a key viscosity for flow perpendicular to the director is $\eta_c = \frac{1}{2}\alpha_4$. Since the strength of the secondary flow depends on the balance between the driving centrifugal force (proportional to the primary flow) and the viscous damping in the cross-plane, the final intensity of the Dean vortices is exquisitely sensitive to the ratio of these anisotropic viscosities, and thus to the underlying Leslie coefficients.

We have seen the Leslie coefficients at work in a dozen different contexts, acting as the master parameters of nematic hydrodynamics. But a persistent question remains: where do these coefficients, these $\alpha$'s, fundamentally come from? The answer takes us to the very heart of statistical mechanics and reveals a profound connection between the deterministic world of fluid dynamics and the chaotic, random world of molecular motion.

Any transport coefficient, like viscosity, is ultimately a measure of how the system dissipates energy and transports momentum. At the microscopic level, these processes are the result of countless random molecular collisions and interactions—the thermal "jiggles" of the system's constituents. The ​​Fluctuation-Dissipation Theorem​​, one of the deepest results in modern physics, provides the bridge. It states that the way a system responds to a small external push (dissipation) is completely determined by the way its internal random fluctuations behave at equilibrium.

For a simple liquid, the Green-Kubo relation (a specific form of the theorem) tells us that the shear viscosity is proportional to the time integral of the equilibrium autocorrelation function of the shear stress. In plain English: imagine watching the random, momentary fluctuations in pressure on the walls of a container. Viscosity can be calculated by seeing how long a random spike in stress at one moment stays correlated with the stress a short time later.

For a nematic liquid crystal, the story is richer. The system has more ways to fluctuate. Not only can the stress fluctuate, but the local director orientation can also jiggle randomly around its average direction. The Leslie coefficients arise not only from the autocorrelation of stress fluctuations, but also from the ​​cross-correlation​​ between stress fluctuations and director fluctuations. For example, a coefficient like $\alpha_2$ is related to how a random jiggle in the director's orientation is correlated with a momentary shear stress in the fluid. In this beautiful picture, the macroscopic Leslie coefficients, which describe the deterministic dance of the director and flow, are revealed to be the echoes of a frantic, random microscopic ballet of fluctuating molecules. This connection marries the phenomenological brilliance of Leslie and Ericksen with the statistical foundations laid by Einstein, Onsager, and Kubo, providing a truly unified picture of these remarkable materials.