Leslie Effects: Demography, Ecology & Liquid Crystals

• Leslie effects are distinct phenomena in demography, ecology, and liquid-crystal hydrodynamics, each utilizing class- or orientation-specific coefficients to propagate local changes through iterative, transport, or stress coupling.
• In demography, Leslie matrices use fertility and survival coefficients to forecast population structures and growth rates, highlighting the impact of structural propagation in age or category dynamics.
• In ecology and liquid crystals, Leslie effects model prey-dependent regulation and anisotropic viscous flow coupling, respectively, producing stable long-term behaviors through nonlocal consequences.

“Leslie effects” is a context-dependent technical expression rather than a single universally standardized concept. In current arXiv usage represented here, it denotes at least three distinct but structurally related ideas: the propagation of demographic composition by Leslie matrices in population forecasting, prey-dependent predator regulation in Leslie–Gower ecological systems, and anisotropic viscous flow–director coupling in Ericksen–Leslie theories of nematic liquid crystals. A plausible unifying theme is that each usage describes how class- or orientation-specific coefficients generate nonlocal consequences under iteration, transport, or stress coupling (Malafeyev et al., 2024, Korobeinikov et al., 2011, Hieber et al., 17 May 2025).

1. Terminological scope and principal domains

The expression appears in three main technical settings. In demography, it refers to the cumulative population consequences of differences in reproduction, survival, and movement between classes. In ecology, it refers to Leslie–Gower regulation, where predator limitation is controlled by prey abundance through a prey-dependent carrying capacity. In liquid-crystal hydrodynamics, it refers to Leslie viscous effects: the full anisotropic coupling between strain, vorticity, director dynamics, and extra stress (Malafeyev et al., 2024, Korobeinikov et al., 2011, Hieber et al., 2016).

Domain	Leslie effect in the cited literature	Representative object
Demography	Future class structure generated from current structure by fertility and survival coefficients	Leslie matrix $A$ with fertility first row and survival subdiagonal
Ecology	Predator self-limitation through prey-dependent carrying capacity	Leslie–Gower term $qP\left(1-\frac{P}{K(H)}\right)$
Liquid crystals	Anisotropic viscous and orientational flow coupling	Leslie stress $\sigma_L$, corotational derivative $N$

This multiplicity is not accidental. The demographic and ecological usages are historically connected to Leslie’s structured-population modeling, whereas the liquid-crystal usage comes from the Ericksen–Leslie constitutive theory. The same surname therefore indexes different mathematical traditions.

2. Demographic Leslie effects and matrix population dynamics

In demographic analysis, the Leslie matrix is a nonnegative projection matrix whose first row contains class-specific fertility terms and whose subdiagonal contains survival probabilities. A continuous age variable $x$ is discretized into classes $i$, the population at time $t$ is recorded in a state vector $n(t)$, and one-step evolution is written as

$n(t+1)=A\,n(t).$

In the three-class example,

$n_2(t+1)=P_1 n_1(t),\qquad n_3(t+1)=P_2 n_2(t),$

and

$qP\left(1-\frac{P}{K(H)}\right)$0

so that

$qP\left(1-\frac{P}{K(H)}\right)$1

Within this framework, a change in any $qP\left(1-\frac{P}{K(H)}\right)$2 alters inflow into the youngest class, and a change in any $qP\left(1-\frac{P}{K(H)}\right)$3 alters cohort transmission into the next class; repeated iteration propagates those local changes through the full future age composition (Malafeyev et al., 2024).

If the matrix is constant, then

$qP\left(1-\frac{P}{K(H)}\right)$4

The same framework supports the standard long-run interpretation: for a constant nonnegative Leslie matrix, the dominant eigenvalue $qP\left(1-\frac{P}{K(H)}\right)$5 determines asymptotic growth, and if

$qP\left(1-\frac{P}{K(H)}\right)$6

then $qP\left(1-\frac{P}{K(H)}\right)$7 is the stable age distribution and asymptotically

$qP\left(1-\frac{P}{K(H)}\right)$8

The interpretation is standard: $qP\left(1-\frac{P}{K(H)}\right)$9 implies long-run growth, $\sigma_L$0 decline, and $\sigma_L$1 stationarity. The cited demographic paper does not derive a full eigenvalue theory, but it explicitly situates Leslie matrices as tools for calculating internal growth rates and forecasting demographic change over time (Malafeyev et al., 2024).

A central point in that literature is that the classification need not be biological age in the narrow sense. Populations may also be classified by gender, employment status, maturity, or other demographic or socioeconomic categories. In that generalized sense, Leslie effects include compositional shifts produced by social transitions as well as biological ageing. This suggests that “Leslie effect” in demography is best understood as a structural propagation effect: differences in class-specific coefficients do not remain local once the matrix is iterated (Malafeyev et al., 2024).

The same paper embeds Leslie-type projections in a two-country policy game. Each country’s population is divided into children, adults, and elderly, and the update rule is extended to

$\sigma_L$2

where $\sigma_L$3 is an immigration vector. Two strategies are available: investment in survival, which modifies the Leslie matrix itself, and investment in immigration, which modifies the additive immigration inflow. When country $\sigma_L$4 moves first, backward induction yields immigration by $\sigma_L$5 and survival by $\sigma_L$6 as the sequentially rational outcome, and the corresponding normal form contains two Nash equilibria,

$\sigma_L$7

The substantive point is that Leslie effects become policy-relevant once fertility, survival, and immigration are treated as strategic levers (Malafeyev et al., 2024).

3. Ecological Leslie–Gower effects

In ecological modeling, the Leslie–Gower effect is the replacement of direct predator growth by predator logistic growth with a prey-dependent carrying capacity. In the classical two-species form,

$\sigma_L$8

The predator equation can be written as

$\sigma_L$9

so predator carrying capacity is proportional to prey abundance. This is the core ecological Leslie–Gower mechanism: unlike Lotka–Volterra formulations, predator increase is bounded by a prey-controlled carrying capacity rather than only prey encounter (Korobeinikov et al., 2011).

That mechanism extends to trophic chains. For a chain with basal prey $N$0 and predator levels $N$1, each predator is limited by the abundance of the level below it: $N$2 For positive parameter values, the cited food-chain model has a unique positive coexistence equilibrium $N$3, and this equilibrium is globally asymptotically stable in the positive orthant. The proof uses a Lyapunov function with recursively chosen weights $N$4 so that mixed trophic terms cancel exactly. In this setting, Leslie–Gower effects regularize trophic dynamics: the long-run behavior is a unique stable coexistence state rather than oscillatory or multistable behavior (Korobeinikov et al., 2011).

A spatial extension appears in modified Leslie–Gower reaction–diffusion systems. There the predator equation takes the form

$N$5

with $N$6 representing alternative food or background carrying capacity. The modification is decisive: even when $N$7, predators may persist at the prey-free state $N$8. For traveling waves connecting the prey-present state $N$9 to either coexistence or prey-free states, the minimal wave speed is

$x$0

and this threshold is both necessary and sufficient for semi-traveling waves in the general framework studied there (Li et al., 2023).

The ecological literature represented here also combines Leslie–Gower regulation with strong Allee effects, cooperative hunting, and constant prey placement. In the nondimensional planar model,

$x$1

positive equilibria satisfy $x$2 and reduce to a cubic condition in $x$3. Under

$x$4

with specific parameter choices $x$5 and $x$6, the model has a triple positive equilibrium $x$7. The paper then analyzes node and higher-degeneracy behavior near that equilibrium and proves bifurcation with residual dimension $x$8 using cooperative hunting and placement rate as branching parameters (Zhao, 2024).

4. Leslie effects in nematic liquid crystals

In liquid-crystal hydrodynamics, Leslie effects are the full anisotropic viscous couplings between fluid deformation and director dynamics. In the classical incompressible-isothermal formulation, the stress is written as

$x$9

with Leslie stress

$i$0

where

$i$1

These terms encode rotation by vorticity, stretching by the symmetric strain, and the resulting anisotropic extra stress (Hieber et al., 17 May 2025).

A thermodynamically adapted parameterization used in the incompressible isotropic non-isothermal setting decomposes the total stress as

$i$2

with

$i$3

The director equation contains the explicit Leslie couplings

$i$4

and the general Leslie stress contains terms involving $i$5, $i$6, and $i$7. The key point is that this framework treats the full Leslie stress without imposing structural relations such as Parodi’s law for well-posedness (Hieber et al., 2016).

Not all papers labeled Ericksen–Leslie retain those full effects. A simplified reduced model on $i$8 keeps only

$i$9

with isotropic Newtonian viscosity in the fluid and no $t$0, no $t$1, no $t$2, no corotational derivative $t$3, and no Parodi relation. That distinction is fundamental: the simplified model is analytically useful for weak compactness and concentration-cancellation, but it does not represent full Leslie viscous effects in the classical constitutive sense (Kortum, 2019).

5. Thermodynamics, well-posedness, and limiting theories

The PDE theory of Leslie effects is shaped by energy and entropy structure. In the non-isothermal general Ericksen–Leslie system, entropy production satisfies

$t$4

so the classical viscous dissipation $t$5 is reinterpreted as part of the entropy production. The compressible derivation admits a general Leslie stress with coefficients $t$6, and the second law yields explicit algebraic inequalities on these coefficients together with

$t$7

for anisotropic heat conduction (Anna et al., 2017).

For the full Ericksen–Leslie system with general Oseen–Frank elasticity, global weak existence in $t$8 has been proved under physical constraints on the Leslie coefficients, including Parodi’s relation

$t$9

and positivity conditions expressed through the combinations $n(t)$0. In that setting, the Leslie stress contributes a nonnegative quadratic form in the strain tensor $n(t)$1, and this dissipative structure underlies both local strong theory and global weak compactness (Wang et al., 2013).

A related $n(t)$2 well-posedness theory constructs an approximation that preserves the dissipated energy of the full Leslie system rather than using a Ginzburg–Landau penalty. The exact energy law contains

$n(t)$3

and the corresponding algebraic condition on $n(t)$4 is necessary and sufficient for nonnegative Leslie dissipation (Wang et al., 2012).

Weak-solution uniqueness is more delicate. For the full $n(t)$5 Ericksen–Leslie system with general Leslie stress and general Oseen–Frank energy, uniqueness of weak solutions in $n(t)$6 has been proved provided the Frank coefficients are close to a positive constant. The proof works one order below the natural energy level and relies on rewriting the Leslie stress into principal and perturbative parts so that cancellation with the director equation becomes explicit (Li et al., 2014).

A separate line of work derives Leslie effects from tensorial theories. In a $n(t)$7 singular limit from Beris–Edwards/Landau–de Gennes dynamics, weak solutions converge to a director-based Ericksen–Leslie system with explicitly identified coefficients

$n(t)$8

and specific formulas for $n(t)$9. This shows that Leslie-type stress and director-flow coupling can emerge rigorously from a $n(t+1)=A\,n(t).$0-tensor model in a weak singular regime, although the limiting weak solution may contain singular points (Xin et al., 2021).

6. Computation, parametrization, and recurrent misconceptions

Computational work on Leslie-type structures is comparably diverse. In human demography, direct estimation of an age/sex Leslie-type matrix from population counts can be formulated as a constrained quadratic program. “Wood’s Method” fits a $n(t+1)=A\,n(t).$1 transition matrix from age-sex counts alone, using demographic constraints on fertility, survival, and migration. Applied to 3,120 U.S. counties, the reported mean absolute percentage error is $n(t+1)=A\,n(t).$2 overall, $n(t+1)=A\,n(t).$3 for counties above 50,000 population, and $n(t+1)=A\,n(t).$4 for counties below 50,000 (Sprague, 2012).

A related algebraic result shows that general Lefkovitch or Usher population matrices are similar to pseudo-Leslie matrices. If

$n(t+1)=A\,n(t).$5

then $n(t+1)=A\,n(t).$6 has Leslie form but its first-row entries may be negative, so it is a computational rather than directly biological object. Spectral quantities are preserved under this similarity transformation, but Kolmogorov–Sinai entropies of the associated Markov chains need not coincide even when the dominant growth ratio is the same (Alves et al., 2015).

In liquid crystals, recent numerical analysis has focused on preserving Leslie structure at the discrete level. A geometric reformulation of the full Ericksen–Leslie model rewrites the director equation in rotational form so that the evolution is intrinsically tangent to the unit sphere. On that basis, rotational discrete-gradient schemes preserve $n(t+1)=A\,n(t).$7 and satisfy an unconditional discrete energy law for the full model with general anisotropic Oseen–Frank elasticity and full Leslie stress (Wang et al., 15 Jun 2026).

Several misconceptions recur across these literatures. One is that Leslie effects form a single concept; the evidence here instead shows a surname shared by distinct mathematical traditions. Another is that demographic Leslie effects are restricted to biological age classes; the generalized matrix framework includes demographic and socioeconomic categories as well (Malafeyev et al., 2024). A third is that every Ericksen–Leslie paper studies the full Leslie constitutive law; simplified systems may omit all six Leslie coefficients and retain only transport plus Ericksen elastic stress (Kortum, 2019). A fourth is that predator persistence in Leslie–Gower models requires prey persistence; the modified reaction–diffusion theory with $n(t+1)=A\,n(t).$8 admits prey-free predator states and traveling waves connecting $n(t+1)=A\,n(t).$9 to $n_2(t+1)=P_1 n_1(t),\qquad n_3(t+1)=P_2 n_2(t),$0 (Li et al., 2023).

Taken together, these usages show that “Leslie effects” names a family of structurally rich coefficient-driven mechanisms rather than a single disciplinary object. In demography the effect is compositional propagation under matrix iteration; in ecology it is prey-dependent predator regulation and its spatial consequences; in liquid-crystal hydrodynamics it is anisotropic viscous and orientational coupling constrained by mechanics, thermodynamics, and geometry.