skip to main content

Title: Learning hydrodynamic equations for active matter from particle simulations and experiments

Recent advances in high-resolution imaging techniques and particle-based simulation methods have enabled the precise microscopic characterization of collective dynamics in various biological and engineered active matter systems. In parallel, data-driven algorithms for learning interpretable continuum models have shown promising potential for the recovery of underlying partial differential equations (PDEs) from continuum simulation data. By contrast, learning macroscopic hydrodynamic equations for active matter directly from experiments or particle simulations remains a major challenge, especially when continuum models are not known a priori or analytic coarse graining fails, as often is the case for nondilute and heterogeneous systems. Here, we present a framework that leverages spectral basis representations and sparse regression algorithms to discover PDE models from microscopic simulation and experimental data, while incorporating the relevant physical symmetries. We illustrate the practical potential through a range of applications, from a chiral active particle model mimicking nonidentical swimming cells to recent microroller experiments and schooling fish. In all these cases, our scheme learns hydrodynamic equations that reproduce the self-organized collective dynamics observed in the simulations and experiments. This inference framework makes it possible to measure a large number of hydrodynamic parameters in parallel and directly from video data.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Hydrodynamic theories effectively describe many-body systems out of equilibrium in terms of a few macroscopic parameters. However, such parameters are difficult to determine from microscopic information. Seldom is this challenge more apparent than in active matter, where the hydrodynamic parameters are in fact fields that encode the distribution of energy-injecting microscopic components. Here, we use active nematics to demonstrate that neural networks can map out the spatiotemporal variation of multiple hydrodynamic parameters and forecast the chaotic dynamics of these systems. We analyze biofilament/molecular-motor experiments with microtubule/kinesin and actin/myosin complexes as computer vision problems. Our algorithms can determine how activity and elastic moduli change as a function of space and time, as well as adenosine triphosphate (ATP) or motor concentration. The only input needed is the orientation of the biofilaments and not the coupled velocity field which is harder to access in experiments. We can also forecast the evolution of these chaotic many-body systems solely from image sequences of their past using a combination of autoencoders and recurrent neural networks with residual architecture. In realistic experimental setups for which the initial conditions are not perfectly known, our physics-inspired machine-learning algorithms can surpass deterministic simulations. Our study paves the way for artificial-intelligence characterization and control of coupled chaotic fields in diverse physical and biological systems, even in the absence of knowledge of the underlying dynamics.

    more » « less
  2. Abstract

    Whereas self-propelled hard discs undergo motility-induced phase separation, self-propelled rods exhibit a variety of nonequilibrium phenomena, including clustering, collective motion, and spatio-temporal chaos. In this work, we present a theoretical framework representing active particles by continuum fields. This concept combines the simplicity of alignment-based models, enabling analytical studies, and realistic models that incorporate the shape of self-propelled objects explicitly. By varying particle shape from circular to ellipsoidal, we show how nonequilibrium stresses acting among self-propelled rods destabilize motility-induced phase separation and facilitate orientational ordering, thereby connecting the realms of scalar and vectorial active matter. Though the interaction potential is strictly apolar, both, polar and nematic order may emerge and even coexist. Accordingly, the symmetry of ordered states is a dynamical property in active matter. The presented framework may represent various systems including bacterial colonies, cytoskeletal extracts, or shaken granular media.

    more » « less
  3. Agent-based models provide a flexible framework that is frequently used for modelling many biological systems, including cell migration, molecular dynamics, ecology and epidemiology. Analysis of the model dynamics can be challenging due to their inherent stochasticity and heavy computational requirements. Common approaches to the analysis of agent-based models include extensive Monte Carlo simulation of the model or the derivation of coarse-grained differential equation models to predict the expected or averaged output from the agent-based model. Both of these approaches have limitations, however, as extensive computation of complex agent-based models may be infeasible, and coarse-grained differential equation models can fail to accurately describe model dynamics in certain parameter regimes. We propose that methods from the equation learning field provide a promising, novel and unifying approach for agent-based model analysis. Equation learning is a recent field of research from data science that aims to infer differential equation models directly from data. We use this tutorial to review how methods from equation learning can be used to learn differential equation models from agent-based model simulations. We demonstrate that this framework is easy to use, requires few model simulations, and accurately predicts model dynamics in parameter regions where coarse-grained differential equation models fail to do so. We highlight these advantages through several case studies involving two agent-based models that are broadly applicable to biological phenomena: a birth–death–migration model commonly used to explore cell biology experiments and a susceptible–infected–recovered model of infectious disease spread. 
    more » « less
  4. Learning multi-agent system dynamics has been extensively studied for various real-world applications, such as molecular dynamics in biology, multi-body system in physics, and particle dynamics in material science. Most of the existing models are built to learn single system dynamics, which learn the dynamics from observed historical data and predict the future trajectory. In practice, however, we might observe multiple systems that are generated across different environments, which differ in latent exogenous factors such as temperature and gravity. One simple solution is to learn multiple environment-specific models, but it fails to exploit the potential commonalities among the dynamics across environments and offers poor prediction results where per-environment data is sparse or limited. Here, we present GG-ODE (Generalized Graph Ordinary Differential Equations), a machine learning framework for learning continuous multi-agent system dynamics across environments. Our model learns system dynamics using neural ordinary differential equations (ODE) parameterized by Graph Neural Networks (GNNs) to capture the continuous interaction among agents. We achieve the model generalization by assuming the dynamics across different environments are governed by common physics laws that can be captured via learning a shared ODE function. The distinct latent exogenous factors learned for each environment are incorporated into the ODE function to account for their differences. To improve model performance, we additionally design two regularization losses to (1) enforce the orthogonality between the learned initial states and exogenous factors via mutual information minimization; and (2) reduce the temporal variance of learned exogenous factors within the same system via contrastive learning. Experiments over various physical simulations show that our model can accurately predict system dynamics, especially in the long range, and can generalize well to new systems with few observations. 
    more » « less
  5. null (Ed.)
    In active matter systems, self-propelled particles can self-organize to undergo collective motion, leading to persistent dynamical behavior out of equilibrium. In cells, cytoskeletal filaments and motor proteins form complex structures important for cell mechanics, motility, and division. Collective dynamics of cytoskeletal systems can be reconstituted using filament gliding experiments, in which cytoskeletal filaments are propelled by surface-bound motor proteins. These experiments have observed diverse dynamical states, including flocks, polar streams, swirling vortices, and single-filament spirals. Recent experiments with microtubules and kinesin motor proteins found that the collective behavior of gliding filaments can be tuned by altering the concentration of the crowding macromolecule methylcellulose in solution. Increasing the methylcellulose concentration reduced filament crossing, promoted alignment, and led to a transition from active, isotropically oriented filaments to locally aligned polar streams. This emergence of collective motion is typically explained as an increase in alignment interactions by Vicsek-type models of active polar particles. However, it is not yet understood how steric interactions and bending stiffness modify the collective behavior of active semiflexible filaments. Here we use simulations of driven filaments with tunable soft repulsion and rigidity in order to better understand how the interplay between filament flexibility and steric effects can lead to different active dynamic states. We find that increasing filament stiffness decreases the probability of filament alignment, yet increases collective motion and long-range order, in contrast to the assumptions of a Vicsek-type model. We identify swirling flocks, polar streams, buckling bands, and spirals, and describe the physics that govern transitions between these states. In addition to repulsion and driving, tuning filament stiffness can promote collective behavior, and controls the transition between active isotropic filaments, locally aligned flocks, and polar streams. 
    more » « less