Instability of microbial droplets growing on viscous substrates

This paper presents a mathematical model of a microbial droplet growing on a viscous fluid, revealing that while growth forces stabilize its axisymmetric expansion, buoyancy forces induce morphological instability, a finding that aligns with experimental observations.

The Gist

Imagine a tiny, living island made of yeast cells floating on top of a thick, sticky pool of honey. This isn't just a static island; it's a growing city. The cells are eating the nutrients dissolved in the honey below them, getting bigger, and pushing outward.

This paper is a mathematical detective story about what happens when this living island grows. The authors, Vicente Gomez Herrera and Scott Weady, built a computer model to figure out why these microbial islands sometimes stay perfectly round and other times break apart into weird, finger-like shapes.

Here is the breakdown of their discovery, using some everyday analogies:

1. The Setup: The Sticky Pancake

Think of the nutrient-rich fluid as a giant, deep pool of thick syrup. The microbial colony is a flat pancake sitting on the surface.

• The Growth: The yeast cells are hungry. They eat the sugar (nutrients) in the syrup right under them. As they eat, they multiply and push outward, trying to expand the pancake.
• The Physics: Because the syrup is so thick (viscous), the yeast can't just slide easily. Their growth creates pressure, like someone trying to blow up a balloon inside a jar. This pressure pushes the syrup down and out.

2. The Two Forces at War

The paper identifies two main characters fighting for control over the shape of the colony:

• Character A: The "Growth Force" (The Stabilizer)
  Imagine the yeast cells are like a team of construction workers building a wall. They push outward evenly. The authors found that this growth pressure acts like a tightening belt. It wants to keep the colony perfectly round and smooth. If a little bump tries to form, the growth pressure squishes it back down. It's a stabilizing force.

• Character B: The "Buoyancy Force" (The Destabilizer)
  This is the tricky part. As the yeast eats the sugar, the syrup directly under them becomes "lighter" (less dense) because the sugar is gone. Meanwhile, the syrup further away is still heavy with sugar.
  Think of this like a hot air balloon or oil floating on water. The lighter syrup wants to rise, and the heavier syrup wants to sink. This creates a swirling current (a vortex) underneath the colony.
  The authors found that these swirling currents act like a bully. They grab the edges of the colony and pull them into jagged fingers. If this force gets too strong, it overpowers the "tightening belt" of the growth, and the perfect circle breaks apart.

3. The Mathematical Magic Trick

Solving how a 3D fluid moves under a 2D growing pancake is usually a nightmare for mathematicians. It involves solving complex equations for the entire volume of the syrup.

The authors used a clever trick. Instead of tracking every drop of syrup in the deep pool, they realized they could describe the whole system using only the surface of the pancake.

• The Analogy: Imagine you are trying to predict how a trampoline bounces. Instead of calculating the tension of every single spring inside the trampoline, you just look at the fabric on top. The authors turned the problem into a set of equations that only "live" on the surface of the colony. This made the math much cleaner and allowed them to find the exact tipping point where the colony goes from round to ragged.

4. The Big Discovery: The Tipping Point

The most exciting part of the paper is the "Stability Threshold."

The authors calculated that there is a specific "recipe" for when the colony stays round and when it breaks.

• If the syrup is very thick (high viscosity), the buoyancy currents are weak. The "Growth Force" wins, and the colony stays a smooth, expanding circle.
• If the syrup is thinner (lower viscosity) or the yeast grows very fast, the "Buoyancy Force" gets stronger. Eventually, it crosses a critical line (they calculated a specific number, 96, in their units).
• The Result: Once you cross that line, the currents become too strong. The smooth circle becomes unstable, and the colony starts sprouting fingers or breaking into smaller droplets.

5. Why Does This Matter?

You might wonder, "Who cares about yeast on honey?"

• Real World: This happens in nature with algae blooms in the ocean, in industrial fermentation (like making kombucha), and even in cleaning up oil spills with bacteria.
• The Connection: The authors compared their math to real experiments with yeast. The experiments showed that when the fluid gets less viscous, the yeast colonies do exactly what the math predicted: they stop being round and start forming messy, finger-like patterns.

Summary

In short, this paper explains that a growing microbial colony is a tug-of-war between growth (which wants to keep things round and tidy) and buoyancy currents (which want to stir things up and make them messy).

The authors created a mathematical "rulebook" that predicts exactly when the messy currents will win. It's a beautiful example of how math can explain why living things change shape, turning a complex fluid dynamics problem into a clear story about balance and instability.

Technical Summary

1. Problem Statement

The paper investigates the morphological stability of a microbial colony (modeled as a 2D "droplet") growing on the surface of a semi-infinite, viscous fluid substrate. This system is relevant to various biological and industrial contexts, such as yeast colonies on viscous media, algal blooms, and biofilms on oil droplets.

The core physical challenge involves a coupled feedback loop between:

1. Growth: Microbial consumption of nutrients leads to local expansion and the generation of an isotropic growth pressure.
2. Fluid Dynamics: The growth pressure drives viscous flows in the underlying fluid.
3. Buoyancy: Nutrient consumption creates density gradients in the fluid (denser where nutrients are depleted), inducing buoyancy-driven flows (Rayleigh-Bénard type convection).

The authors aim to determine how these competing forces—growth-induced expansion and buoyancy-driven convection—affect the stability of the colony's shape, specifically whether the colony remains a radially expanding disk or develops interfacial instabilities (fingering/roughness).

2. Methodology

Mathematical Formulation

The authors develop a free-boundary problem governed by:

• Bulk Fluid: The 3D Stokes equations (low Reynolds number) with a Boussinesq approximation for buoyancy forces arising from nutrient concentration gradients.
• Nutrient Transport: Advection-diffusion equation for nutrient concentration.
• Surface Kinematics: A constitutive law relating surface divergence to the growth rate $\gamma(c)$, and a kinematic condition for the moving boundary.

Asymptotic Reduction

Through dimensional analysis, the authors identify that for typical microbial experiments (slow growth, high viscosity), the Reynolds number ($Re$) and Péclet number ($Pe$) are small ($Re, Pe \ll 1$). This allows the system to be reduced to linear, elliptic PDEs in the bulk fluid volume, with all nonlinearity and time dependence confined to the boundary conditions and geometry.

Boundary Integral Reformulation

To avoid the computational complexity of solving 3D PDEs coupled with a moving 2D boundary, the authors reformulate the problem as a system of integro-differential equations defined solely on the microbial domain ($\Omega$).

• They utilize Fourier transforms in the horizontal plane to solve the bulk equations analytically.
• The solution is expressed in terms of surface quantities: growth pressure ($p$), surface velocity ($\mathbf{u}$), and a nutrient density source term ($\sigma$).
• The resulting system involves singular integral operators (Single Layer Potential $S_\Omega$, Hypersingular Operator $N_\Omega$, and Volume Potentials $B_\Omega, V_\Omega$) acting on the domain.

Linear Stability Analysis

The authors analyze the stability of the axisymmetric solution (a radially expanding disk) by:

1. Deriving exact axisymmetric solutions for two limiting cases: Growth-dominated ($Ra = 0$) and Buoyancy-dominated ($Ra \to \infty$).
2. Introducing small perturbations to the boundary shape using a conformal mapping approach on a unit disk reference domain.
3. Computing the shape derivatives of the integral operators to determine the evolution of perturbation modes ($m$).
4. Deriving a stability criterion based on the Metabolic Rayleigh Number ($Ra$), which compares buoyancy forces to growth forces.

3. Key Contributions

1. Integro-Differential Framework: The paper provides a rigorous mathematical reformulation of the microbial droplet problem, reducing a complex 3D free-boundary problem to a 2D system of integro-differential equations. This significantly reduces degrees of freedom and isolates the spectral properties of the integral operators.
2. Spectral Analysis: The authors utilize projected spherical harmonics to characterize the eigenfunctions and eigenvalues of the integral operators ($S_D, N_D, B_D, V_D$) on a unit disk. This allows for the analytical derivation of axisymmetric solutions and the calculation of stability coefficients.
3. Mechanism of Instability: The study identifies a distinct competition between two physical mechanisms:
   • Growth forces act to stabilize the axisymmetric shape.
   • Buoyancy forces (driven by nutrient depletion) act to destabilize the shape.
4. Critical Threshold: The authors derive a critical threshold for the metabolic Rayleigh number ($Ra^*$) above which the axisymmetric solution becomes unstable. They show that the critical value depends on the perturbation mode $m$ and asymptotically approaches a constant ($\beta Ra^* \approx 96$) for high-frequency modes.

4. Key Results

• Axisymmetric Solutions:

  • In the growth-dominated regime, the pressure is positive, and the velocity field resembles a stagnation point flow (fluid pulled from below, pushed outward).
  • In the buoyancy-dominated regime, the pressure becomes negative (compressive), and the flow forms a vortex ring beneath the droplet, with the core near the droplet edge.
  • The general solution is a linear superposition of these two regimes.

• Stability Analysis:

  • Growth Stabilization: For $Ra=0$, all non-trivial perturbation modes ($m \ge 2$) have negative growth rates, meaning the circular shape is linearly stable. The stabilizing effect scales as $O(m^{1/2})$.
  • Buoyancy Destabilization: For $Ra \to \infty$, high-frequency modes ($m \ge 2$) have positive growth rates, leading to instability. The destabilizing effect also scales as $O(m^{1/2})$.
  • Competition: The net growth rate of a perturbation mode $m$ is given by $\dot{\varepsilon}/\varepsilon \approx -1/2 + \sigma_m^g + \beta Ra \sigma_m^b$. Instability occurs when the buoyancy term overcomes the growth stabilization.

• Transition to Instability:

  • The system becomes unstable when the metabolic Rayleigh number exceeds a critical value: $Ra > Ra_m^*$.
  • As $m \to \infty$, the critical threshold converges to $\beta Ra^* \approx 96$.
  • At this threshold, the pressure inside the droplet becomes negative everywhere, and a local maximum in radial velocity appears outside the droplet edge.

• Experimental Validation:

  • The predicted critical threshold ($\beta Ra \approx 96$) is of the correct order of magnitude for experimental observations of yeast colonies on viscous substrates (e.g., Atis et al., 2019), where instability occurs at viscosities around 450 Pa·s ($Ra \approx 50$).
  • The model explains the observed "fingering" instability as a result of high-frequency modes growing fastest due to buoyancy.

5. Significance

• Theoretical Insight: The paper provides the first rigorous mechanical explanation for the onset of pattern formation in microbial colonies on viscous fluids. It resolves the ambiguity of how metabolic flows connect to colony morphology.
• Predictive Power: The identification of a single dimensionless parameter ($\beta Ra$) governing stability allows for predictive testing in laboratory settings by varying viscosity, droplet size, or growth rates.
• Methodological Advancement: The use of boundary integral methods and spectral analysis of singular operators offers a robust framework for studying other active matter systems on fluid interfaces, such as lipid membranes or active particle layers.
• Future Directions: The authors note that while the linear theory explains the onset of instability, future work is needed to address large deformations (non-linear regime), finite domain effects, and the potential breakup of droplets into multiple colonies.

In summary, this work bridges the gap between microscopic nutrient consumption and macroscopic morphological instability, demonstrating that while growth tends to smooth the interface, buoyancy-driven flows induced by nutrient depletion can drive the formation of complex patterns.