Learning the Physics of Pattern Formation from Images

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries1 and electrocatalysts2. Experimental characterizations of such materials by operando microscopy produce rich image datasets3-6, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation7. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

Perspective: How to overcome dynamical density functional theory

We argue in favour of developing a comprehensive dynamical theory for rationalizing, predicting, designing, and machine learning nonequilibrium phenomena that occur in soft matter. To give guidance for navigating the theoretical and practical challenges that lie ahead, we discuss and exemplify the limitations of dynamical density functional theory (DDFT). Instead of the implied adiabatic sequence of equilibrium states that this approach provides as a makeshift for the true time evolution, we posit that the pending theoretical tasks lie in developing a systematic understanding of the dynamical functional relationships that govern the genuine nonequilibrium physics. While static density functional theory gives a comprehensive account of the equilibrium properties of many-body systems, we argue that power functional theory is the only present contender to shed similar insights into nonequilibrium dynamics, including the recognition and implementation of exact sum rules that result from the Noether theorem. As a demonstration of the power functional point of view, we consider an idealized steady sedimentation flow of the three-dimensional Lennard-Jones fluid and machine-learn the kinematic map from the mean motion to the internal force field. The trained model is capable of both predicting and designing the steady state dynamics universally for various target density modulations. This demonstrates the significant potential of using such techniques in nonequilibrium many-body physics and overcomes both the conceptual constraints of DDFT as well as the limited availability of its analytical functional approximations.

Perspective: New directions in dynamical density functional theory

Classical dynamical density functional theory (DDFT) has become one of the central modeling approaches in nonequilibrium soft matter physics. Recent years have seen the emergence of novel and interesting fields of application for DDFT. In particular, there has been a remarkable growth in the amount of work related to chemistry. Moreover, DDFT has stimulated research on other theories such as phase field crystal models and power functional theory. In this perspective, we summarize the latest developments in the field of DDFT and discuss a variety of possible directions for future research.

Bayesian modeling of pattern formation from one snapshot of pattern

Partial differential equations (PDEs) have been widely used to reproduce patterns in nature and to give insight into the mechanism underlying pattern formation. Although many PDE models have been proposed, they rely on the pre-request knowledge of physical laws and symmetries, and developing a model to reproduce a given desired pattern remains difficult. We propose a method, referred to as Bayesian modeling of PDEs (BM-PDEs), to estimate the best dynamical PDE for one snapshot of a objective pattern under the stationary state without ground truth. We apply BM-PDEs to nontrivial patterns, such as quasicrystals (QCs), a double gyroid, and Frank-Kasper structures. We also generate three-dimensional dodecagonal QCs from a PDE model. This is done by using the estimated parameters for the Frank-Kasper A15 structure, which closely approximates the local structures of QCs. Our method works for noisy patterns and the pattern synthesized without the ground-truth parameters, which are required for the application toward experimental data.

Interfacial Resistive Switching by Multiphase Polarization in Ion-Intercalation Nanofilms

Nonvolatile resistive-switching (RS) memories promise to revolutionize hardware architectures with in-memory computing. Recently, ion-interclation materials have attracted increasing attention as potential RS materials for their ion-modulated electronic conductivity. In this Letter, we propose RS by multiphase polarization (MP) of ion-intercalated thin films between ion-blocking electrodes, in which interfacial phase separation triggered by an applied voltage switches the electron-transfer resistance. We develop an electrochemical phase-field model for simulations of coupled ion-electron transport and ion-modulated electron-transfer rates and use it to analyze the MP switching current and time, resistance ratio, and current-voltage response. The model is able to reproduce the complex cyclic voltammograms of lithium titanate (LTO) memristors, which cannot be explained by existing models based on bulk dielectric breakdown. The theory predicts the achievable switching speeds for multiphase ion-intercalation materials and could be used to guide the design of high-performance MP-based RS memories.

From predicting to learning dissipation from pair correlations of active liquids

Active systems, which are driven out of equilibrium by local non-conservative forces, can adopt unique behaviors and configurations. Towards designing such materials, an important challenge is to precisely connect the static structure of active systems to the dissipation of energy induced by the local driving. Here, we use tools from liquid-state theories and machine learning to take on these challenges. We first demonstrate analytically for an isotropic active matter system that dissipation and pair correlations are closely related when driving forces behave like an active temperature. We then extend a nonequilibrium mean-field framework for predicting these pair correlations which, unlike most existing approaches, is applicable even for strongly interacting particles and far from equilibrium, to predict dissipation in these systems. Based on this theory, we reveal analytically a robust relation between dissipation and structure which holds even as the system approaches a nonequilibrium phase transition. Finally, we construct a neural network which maps static configurations of particles to their dissipation rate without any prior knowledge of the underlying dynamics. Our results open novel perspectives on the interplay between dissipation and organization out-of-equilibrium.

Correlative image learning of chemo-mechanics in phase-transforming solids

Constitutive laws underlie most physical processes in nature. However, learning such equations in heterogeneous solids (for example, due to phase separation) is challenging. One such relationship is between composition and eigenstrain, which governs the chemo-mechanical expansion in solids. Here we developed a generalizable, physically constrained image-learning framework to algorithmically learn the chemo-mechanical constitutive law at the nanoscale from correlative four-dimensional scanning transmission electron microscopy and X-ray spectro-ptychography images. We demonstrated this approach on Li_XFePO₄, a technologically relevant battery positive electrode material. We uncovered the functional form of the composition-eigenstrain relation in this two-phase binary solid across the entire composition range (0 ≤ X ≤ 1), including inside the thermodynamically unstable miscibility gap. The learned relation directly validates Vegard's law of linear response at the nanoscale. Our physics-constrained data-driven approach directly visualizes the residual strain field (by removing the compositional and coherency strain), which is otherwise impossible to quantify. Heterogeneities in the residual strain arise from misfit dislocations and were independently verified by X-ray diffraction line profile analysis. Our work provides the means to simultaneously quantify chemical expansion, coherency strain and dislocations in battery electrodes, which has implications on rate capabilities and lifetime. Broadly, this work also highlights the potential of integrating correlative microscopy and image learning for extracting material properties and physics.

Identification of approximate symmetries in biological development

Virtually all forms of life, from single-cell eukaryotes to complex, highly differentiated multicellular organisms, exhibit a property referred to as symmetry. However, precise measures of symmetry are often difficult to formulate and apply in a meaningful way to biological systems, where symmetries and asymmetries can be dynamic and transient, or be visually apparent but not reliably quantifiable using standard measures from mathematics and physics. Here, we present and illustrate a novel measure that draws on concepts from information theory to quantify the degree of symmetry, enabling the identification of approximate symmetries that may be present in a pattern or a biological image. We apply the measure to rotation, reflection and translation symmetries in patterns produced by a Turing model, as well as natural objects (algae, flowers and leaves). This method of symmetry quantification is unbiased and rigorous, and requires minimal manual processing compared to alternative measures. The proposed method is therefore a useful tool for comparison and identification of symmetries in biological systems, with potential future applications to symmetries that arise during development, as observed in vivo or as produced by mathematical models. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Artificial Intelligence Applied to Battery Research: Hype or Reality?

This is a critical review of artificial intelligence/machine learning (AI/ML) methods applied to battery research. It aims at providing a comprehensive, authoritative, and critical, yet easily understandable, review of general interest to the battery community. It addresses the concepts, approaches, tools, outcomes, and challenges of using AI/ML as an accelerator for the design and optimization of the next generation of batteries—a current hot topic. It intends to create both accessibility of these tools to the chemistry and electrochemical energy sciences communities and completeness in terms of the different battery R&D aspects covered.

Roadmap on emerging concepts in the physical biology of bacterial biofilms: from surface sensing to community formation

Bacterial biofilms are communities of bacteria that exist as aggregates that can adhere to surfaces or be free-standing. This complex, social mode of cellular organization is fundamental to the physiology of microbes and often exhibits surprising behavior. Bacterial biofilms are more than the sum of their parts: single-cell behavior has a complex relation to collective community behavior, in a manner perhaps cognate to the complex relation between atomic physics and condensed matter physics. Biofilm microbiology is a relatively young field by biology standards, but it has already attracted intense attention from physicists. Sometimes, this attention takes the form of seeing biofilms as inspiration for new physics. In this roadmap, we highlight the work of those who have taken the opposite strategy: we highlight the work of physicists and physical scientists who use physics to engage fundamental concepts in bacterial biofilm microbiology, including adhesion, sensing, motility, signaling, memory, energy flow, community formation and cooperativity. These contributions are juxtaposed with microbiologists who have made recent important discoveries on bacterial biofilms using state-of-the-art physical methods. The contributions to this roadmap exemplify how well physics and biology can be combined to achieve a new synthesis, rather than just a division of labor.

Continuum-scale modelling of polymer blends using the Cahn-Hilliard equation: transport and thermodynamics

The Cahn-Hilliard equation is commonly used to study multi-component soft systems such as polymer blends at continuum scales. We first systematically explore various features of the equation system, which give rise to a deep connection between transport and thermodynamics-specifically that the Gibbs free energy of mixing function is central to formulating a well-posed model. Accordingly, we explore how thermodynamic models from three broad classes of approach (lattice-based, activity-based and perturbation methods) can be incorporated within the Cahn-Hilliard equation and examine how they impact the numerical solution for two model polymer blends, noting that although the analysis presented here is focused on binary mixtures, it is readily extensible to multi-component mixtures. It is observed that, although the predicted liquid-liquid interfacial tension is quite strongly affected, the choice of thermodynamic model has little influence on the development of the morphology.

Non-invasive Inference of Thrombus Material Properties with Physics-Informed Neural Networks

We employ physics-informed neural networks (PINNs) to infer properties of biological materials using synthetic data. In particular, we successfully apply PINNs on inferring permeability and viscoelastic modulus from thrombus deformation data, which can be described by the fourth-order Cahn-Hilliard and Navier-Stokes Equations. In PINNs, the partial differential equations are encoded into a loss function, where partial derivatives can be obtained through automatic differentiation (AD). In addition to tackling the challenge of calculating the fourth-order derivative in the Cahn-Hilliard equation with AD, we introduce an auxiliary network along with the main neural network to approximate the second-derivative of the energy potential term. Our model can simultaneously predict unknown material parameters and velocity, pressure, and deformation gradient fields by merely training with partial information among all data, i.e., phase field and pressure measurements, while remaining highly flexible in sampling within the spatio-temporal domain for data acquisition. We validate our model by numerical solutions from the spectral/hp element method (SEM) and demonstrate its robustness by training it with noisy measurements. Our results show that PINNs can infer the material properties from noisy synthetic data, and thus they have great potential for inferring these properties from experimental multi-modality and multi-fidelity data.

Landscape-induced spatial oscillations in population dynamics

We study the effect that disturbances in the ecological landscape exert on the spatial distribution of a population that evolves according to the nonlocal FKPP equation. Using both numerical and analytical techniques, we characterize, as a function of the interaction kernel, the three types of stationary profiles that can develop near abrupt spatial variations in the environmental conditions vital for population growth: sustained oscillations, decaying oscillations and exponential relaxation towards a flat profile. Through the mapping between the features of the induced wrinkles and the shape of the interaction kernel, we discuss how heterogeneities can reveal information that would be hidden in a flat landscape.

A data-driven computational scheme for the nonlinear mechanical properties of cellular mechanical metamaterials under large deformation

Cellular mechanical metamaterials are a special class of materials whose mechanical properties are primarily determined by their geometry. However, capturing the nonlinear mechanical behavior of these materials, especially those with complex geometries and under large deformation, can be challenging due to inherent computational complexity. In this work, we propose a data-driven multiscale computational scheme as a possible route to resolve this challenge. We use a neural network to approximate the effective strain energy density as a function of cellular geometry and overall deformation. The network is constructed by "learning" from the data generated by finite element calculation of a set of representative volume elements at cellular scales. This effective strain energy density is then used to predict the mechanical responses of cellular materials at larger scales. Compared with direct finite element simulation, the proposed scheme can reduce the computational time up to two orders of magnitude. Potentially, this scheme can facilitate new optimization algorithms for designing cellular materials of highly specific mechanical properties.