Meso-scale Algorithm for Flowing Nematics

Most of my projects investigate some microscopic object (or objects) such as polymers or microbial swimmers immersed in a simple background fluid. Since it’s the dynamics of the objects that we’re typically interested in, we treat the fluid as simply as possible. In order to do this, we often use the MPCD method (described here).

However, not all fluids are simple. Some biological fluids have intrinsic elasticity, while others are made of or include rod-like molecules. These fluids are called “nematic” — from the Greek word for “threadlike”. Since the molecules are rod-like, they can align; and when they align they are partially ordered. Ordered solids are called crystals and so these fluids are called liquid-crystals. Since many biological molecules are rod-like, such as filamentous-actin, their solutions can be liquid crystalline.


Dealing with flowing liquid crystals is notoriously difficult. When the liquid crystal is just the background fluid to the dynamics of the colloids, polymers, or microbes, significant computational effort must be spent to simulate the background. This is why we developed a computationally efficient and versatile multi-particle collision dynamics (MPCD)-based algorithm to simulate nematohydrodynamic baths.


Our key realization was that the MPCD algorithm isn’t limited to solving the Navier-Stokes equation for momentum transport in a fluid, but can really solve any transport equation so long as a suitable “collision” operation (now more generally referred to as an “exchange”) can be defined and the continuum field of interest can be discretized as a particle property. We gave each point-like MPCD particle an orientation that would interact by an order-collision operation based on Maier-Saupe theory for liquid crystals in equilibrium.

This simple algorithm works very well. It is accurate in 2D and 3D; it reproduces defect dynamics; it has an intrinsic elasticity; it can have tumbling or shear aligning regimes when it flows; and it can handle complex boundary conditions and mobile macromolecules, colloids and swimmers.

Lattice Boltzmann Simulations of Swimming of Microbial Organisms

Microorganisms live in a very different world than us — a world commonly submerged in water. But, at such microscopic sizes, swimming in water seems far more like swimming through tar. Understanding the dynamics of how microorganisms travel through liquids can help us control their spread.

Since the microbial world is a liquid one, flows play an important role in setting microbe motion. Bacteria, for instance, commonly have “run-tumble” dynamics: They swim forward for a time then briefly tumble in a new, random direction. Walls and obstacles can impact these dynamics, not only by simply getting in the way but also by changing the flows caused by the bacteria. In this way, even distant objects interact with the swimmers through flow.

Images by Joost de Graaf.

These hydrodynamic interactions are challenging to predict. We’ve worked with Joost de Graaf (University of Edinburgh), who has developed versatile models for computationally simulating the dynamics of self-propelled particles (such as living bacteria) within fluids. These simulations build swimmers out of beads, which appear raspberry-like in shape. The resulting models reproduce well the hydrodynamics of swimming cells, like E. coli or sperm, which means that in the future these models may be able to accurately predict swimming dynamics in complex geometries and flows.

Microbial Swimming in Films

Living near interfaces can be beneficial to microbes. Liquid-air interfaces may be the source of oxygen and sunlight. Solid surfaces accumulate sediments, which can include valuable nutrients. They can also offer anchoring points for the formation of protective biofilm. Films of liquid over a solid wall have both of these types of beneficial interfaces. However, an unavoidable consequence of being near surfaces is a shearing when the film flows. These flows can have a large impact on the swimming dynamics of microbes.

Image by Arnold Mathijssen.

In a series of papers with Arnold Mathijssen and Amin Doostmohammadi, we investigated the swimming dynamics of motile microbes in flowing films. We began by building an analytical model for the flows produced when microbes swim in a film. We were able to reproduce many of the different flows generated by swimmers. These include flagellated bacteria that swim by propelling fluid back and then pushing fluid forward (called “pusher-type swimmers”), and puller-type swimmers that pull fluid toward their head and drag their body forward. We also reproduced ciliated swimmers with their many hair-like flagella.

We then put our swimmer models in a flowing film and looked at the dynamics. We found there is a sweet spot for rinsing bacteria away. When the film isn’t flowing, the microbes are attracted most strongly to the solid surface and they stay near the floor until the flow is strong enough to peel microbial swimmers off the surface. That is to say, at moderate flows bacteria in a film are washed away but that’s not the end of it. If much stronger flows are applied, the flow turns the microbes back down towards the surface, where they are more likely to form biofilms.

Image by Arnold Mathijssen.

Run-tumble Dynamics vs. Hydrodynamic Interactions in Films

Swimmers like E. coli have “run-tumble” dynamics, meaning they swim forward for a time (run stage) then rapidly and randomly change to new direction (tumble stage). While they swim in the run stage, they perturb the surrounding fluid and interact with their surroundings through the flow they generate. This causes swimmers to be attracted to walls.

We asked “When does the hydrodynamic attraction win over the randomizing run-tumble dynamics in a film?” We found if we modelled the dynamics as just an effective constant noise then the swimming strategy of the microbes makes quite a difference. However, if the swimmers have a proper run stage followed by a brief and random tumble stage then the swimming strategy is not nearly as important.

We believe this will be of interest to other scientists in our field because it means that when they make simulations of swimmers they should prioritize getting the run-tumble dynamics correct before attempting to reproduce the more complicated but less important hydrodynamics of microbial swimming.

DNA in Nanoslits

DNA doesn’t just float around. It it tugged on and wrapped up by the protein machinary within cells. All of this pushing and pulling occurs within the confined space of the nucleus. In artificial systems, experimentalists pull on DNA and other semi-flexible biomolecules, stretching them within channels or slits.

Simulation snapshots by Hendrick de Haan.

However, the model relating the extension of DNA to the amount of force applied is for fully 3D or perfectly 2D, like in a film. In strongly confining slits, the walls make it easier to stretch the DNA but not as easy as in 2D. The confining walls make it easier because they obstruct the DNA and so limit the number of conformational states that the DNA can have. Since the number of states sets the entropy and the entropy works against the applied stretching force, confinement makes stretching easier.

With Hendrick de Haan (University of Ontario Institute of Technology), we developed a theory for describing the extension of DNA in a channel of any height as a function of force. This theory is built on the classical theory of Marco-Siggia for DNA in 2D and 3D and uses an effective dimensionality between 2 or 3 when the forces are not too strong. When the forces are very strong a mathematical approximation must be applied for strong and weak confinements.

Entropic Collapse of Coarse-Grained Chromosomes

Figure 1
Illustration by David Goodsell from Chromosome, Cell Cycle, and Entropy.

Entropy can be tricky to understand, but one rule often quoted in pop-culture is that entropy increases as time goes on, which leads to more and more disorder. Depletion forces can therefore be surprising since they seem to break this rule-of-thumb; the depletion force is an entropically induced attraction between large colloids due many smaller surrounding “depletant” particles. By bringing the large colloids together, the attraction orders the larger colloids but this is at the benefit of the many small depletants, which gain increased entropy and are more disordered.


Such entropic forces exist in biological systems. For example, depletion forces lead to fibre bundling, aggregation of red blood cells, and change the reaction kinetics of polymerase chain reactions. They are also ubiquitous in the complex and crowded interior of cells.

Prior to our work, it was hypothesized that depletion forces acting on chromosomes in simple bacteria by surrounding cytoplasmic proteins could explain the collapse from a swollen conformation to a compact state. We performed simulations to explicitly test whether depletion-induced attraction is actually enough to collapse model chromosomes. The resulting simulations gave us strong evidence that depletion effects are enough to collapse chromosomes; however, this collapse is a continuous phase transition rather than a discontinuous one (as is seen in experiments). This suggests that some proteins must act in specific biological ways (rather than as generic depletants) to fully account for the compaction of nucleoids.


To read more, see our paper on the entropic collapse of coarse-grained chromosomes, or a News and Views by Suckjoon Jun discussing our work.


FFF in Microfluidic Channels

The ideal, unified Field-Flow Fractionation (FFF) analysis, suggested that a novel operational-mode exists at the largest particle sizes. Therefore, our collaborators in Dr. Michel Godin’s laboratory constructed a 18micron microfluidic channel to to test this prediction for polystyrene colloids subject to a gravitational field and to determine how well the retention theory would match experimental data across multiple operational-modes.

Our work has recently been submitted for publication. The manuscript presents experimental data, MPCD simulations and a modified, numerical retention theory to build up a consistent picture of FFF of large colloids eluting through microfluidic devices.