-
Multiphase Field Model of Cells on a Substrate
Authors:
Michael Chiang,
Austin Hopkins,
Benjamin Loewe,
Davide Marenduzzo,
M. Cristina Marchetti
Abstract:
Multiphase field models have emerged as an important computational tool for understanding biological tissue while resolving single-cell properties. While they have successfully reproduced many experimentally observed behaviors of living tissue, the theoretical underpinnings have not been fully explored. We show that a two-dimensional version of the model, which is commonly employed to study tissue…
▽ More
Multiphase field models have emerged as an important computational tool for understanding biological tissue while resolving single-cell properties. While they have successfully reproduced many experimentally observed behaviors of living tissue, the theoretical underpinnings have not been fully explored. We show that a two-dimensional version of the model, which is commonly employed to study tissue monolayers, can be derived from a three-dimensional version in the presence of a substrate. We also show how viscous forces, which arise from friction between different cells, can be included in the model. Finally, we numerically simulate a tissue monolayer, and find that intercellular friction tends to solidify the tissue.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Mitigating Density Fluctuations in Particle-based Active Nematic Simulations
Authors:
Timofey Kozhukhov,
Benjamin Loewe,
Tyler N. Shendruk
Abstract:
Understanding active matter has led to new perspectives on biophysics and non-equilibrium dynamics. However, the development of numerical tools for simulating active fluids capable of incorporating non-trivial boundaries or inclusions has lagged behind. Active particle-based methods, which typically excel at this, suffer from large density fluctuations that affect the dynamics of inclusions. To th…
▽ More
Understanding active matter has led to new perspectives on biophysics and non-equilibrium dynamics. However, the development of numerical tools for simulating active fluids capable of incorporating non-trivial boundaries or inclusions has lagged behind. Active particle-based methods, which typically excel at this, suffer from large density fluctuations that affect the dynamics of inclusions. To this end, we advance the Active-Nematic Multi-Particle Collision Dynamics algorithm, a particle-based method for simulating active nematics, by addressing the large density fluctuations that arise from activity. This paper introduces three novel activity formulations that mitigate the coupling between activity and local density. Local density fluctuations are decreased to a level comparable to the passive limit while retaining the phenomenology of active nematics and increasing the active turbulence regime four-fold. These developments extend the technique into a flexible tool for modeling active systems, including solutes and inclusions, with broad applications for the study of biophysical systems.
△ Less
Submitted 31 January, 2024;
originally announced January 2024.
-
Intercellular Friction and Motility Drive Orientational Order in Cell Monolayers
Authors:
Michael Chiang,
Austin Hopkins,
Benjamin Loewe,
M. Cristina Marchetti,
Davide Marenduzzo
Abstract:
Spatiotemporal patterns in multicellular systems are important to understanding tissue dynamics, for instance, during embryonic development and disease. Here, we use a multiphase field model to study numerically the behavior of a near-confluent monolayer of deformable cells with intercellular friction. Varying friction and cell motility drives a solid-liquid transition, and near the transition bou…
▽ More
Spatiotemporal patterns in multicellular systems are important to understanding tissue dynamics, for instance, during embryonic development and disease. Here, we use a multiphase field model to study numerically the behavior of a near-confluent monolayer of deformable cells with intercellular friction. Varying friction and cell motility drives a solid-liquid transition, and near the transition boundary, we find the emergence of nematic order of cell deformation driven by shear-aligning cellular flows. Intercellular friction endows the monolayer with a finite viscosity, which significantly increases the spatial correlation in the flow and, concomitantly, the extent of nematic order. We also show that hexatic and nematic order are tightly coupled and propose a mechanical-geometric model for the colocalization of +1/2 nematic defects and 5-7 disclination pairs, which are the structural defects in the hexatic phase. Such topological defects coincide with regions of high cell-cell overlap, suggesting that they may mediate cellular extrusion from the monolayer, as found experimentally. Our results delineate a mechanical basis for the recent observation of nematic and hexatic order in multicellular collectives in experiments and simulations and pinpoint a generic pathway to couple topological and physical effects in these systems.
△ Less
Submitted 31 October, 2023;
originally announced October 2023.
-
Motility induced phase separation of deformable cells
Authors:
Austin Hopkins,
Benjamin Loewe,
Michael Chiang,
Davide Marenduzzo,
M. Cristina Marchetti
Abstract:
Using a multi-phase field model, we examine how particle deformability, which is a proxy for cell stiffness, affects motility induced phase separation (MIPS). We show that purely repulsive deformable, i.e., squishy, cells phase separate more effectively than their rigid counterparts. This can be understood as due to the fact that deformability increases the effective duration of collisions. In add…
▽ More
Using a multi-phase field model, we examine how particle deformability, which is a proxy for cell stiffness, affects motility induced phase separation (MIPS). We show that purely repulsive deformable, i.e., squishy, cells phase separate more effectively than their rigid counterparts. This can be understood as due to the fact that deformability increases the effective duration of collisions. In addition, the dense regions become increasingly disordered as deformability increases. Our results contextualize the applicability of MIPS to biological systems and have implications for how cells in biological systems may self-organize
△ Less
Submitted 9 August, 2023;
originally announced August 2023.
-
Selective Activation of Aromatic C-H Bonds Catalyzed by Single Gold Atoms at Room Temperature
Authors:
Benjamin Lowe,
Jack Hellerstedt,
Adam Matěj,
Pingo Mutombo,
Dhaneesh Kumar,
Martin Ondráček,
Pavel Jelinek,
Agustin Schiffrin
Abstract:
Selective activation and controlled functionalization of C-H bonds in organic molecules is one of the most desirable processes in synthetic chemistry. Despite progress in heterogeneous catalysis using metal surfaces, this goal remains challenging due to the stability of C-H bonds and their ubiquity in precursor molecules, hampering regioselectivity. Here, we examine the interaction between 9,10-di…
▽ More
Selective activation and controlled functionalization of C-H bonds in organic molecules is one of the most desirable processes in synthetic chemistry. Despite progress in heterogeneous catalysis using metal surfaces, this goal remains challenging due to the stability of C-H bonds and their ubiquity in precursor molecules, hampering regioselectivity. Here, we examine the interaction between 9,10-dicyanoanthracene (DCA) molecules and Au adatoms on a Ag(111) surface at room temperature (RT). Characterization via low-temperature scanning tunneling microscopy, spectroscopy, and noncontact atomic force microscopy, supported by theoretical calculations, revealed the formation of organometallic DCA-Au-DCA dimers, where C atoms at the ends of the anthracene moieties are bonded covalently to single Au atoms. The formation of this organometallic compound is initiated by a regioselective cleaving of C-H bonds at RT. Hybrid quantum mechanics/molecular mechanics calculations show that this regioselective C-H bond cleaving is enabled by an intermediate metal-organic complex which significantly reduces the dissociation barrier of a specific C-H bond. Harnessing the catalytic activity of single metal atoms, this regioselective on-surface C-H activation reaction at RT offers promising routes for future synthesis of functional organic and organometallic materials.
△ Less
Submitted 8 August, 2023;
originally announced August 2023.
-
Local gate control of Mott metal-insulator transition in a 2D metal-organic framework
Authors:
Benjamin Lowe,
Bernard Field,
Jack Hellerstedt,
Julian Ceddia,
Henry L. Nourse,
Ben J. Powell,
Nikhil V. Medhekar,
Agustin Schiffrin
Abstract:
Electron-electron interactions in materials lead to exotic many-body quantum phenomena including Mott metal-insulator transitions (MITs), magnetism, quantum spin liquids, and superconductivity. These phases depend on electronic band occupation and can be controlled via the chemical potential. Flat bands in two-dimensional (2D) and layered materials with a kagome lattice enhance electronic correlat…
▽ More
Electron-electron interactions in materials lead to exotic many-body quantum phenomena including Mott metal-insulator transitions (MITs), magnetism, quantum spin liquids, and superconductivity. These phases depend on electronic band occupation and can be controlled via the chemical potential. Flat bands in two-dimensional (2D) and layered materials with a kagome lattice enhance electronic correlations. Although theoretically predicted, correlated-electron Mott insulating phases in monolayer 2D metal-organic frameworks (MOFs) with a kagome structure have not yet been realised experimentally. Here, we synthesise a 2D kagome MOF on a 2D insulator. Scanning tunnelling microscopy (STM) and spectroscopy reveal a MOF electronic energy gap of ~200 meV, consistent with dynamical mean field theory predictions of a Mott insulator. Combining template-induced (via work function variations of the substrate) and STM probe-induced gating, we locally tune the electron population of the MOF kagome bands and induce Mott MITs. These findings enable technologies based on electrostatic control of many-body quantum phases in 2D MOFs.
△ Less
Submitted 1 May, 2024; v1 submitted 24 May, 2023;
originally announced May 2023.
-
Anisotropic run-and-tumble-turn dynamics
Authors:
Benjamin Loewe,
Tyler N. Shendruk
Abstract:
Run-and-tumble processes successfully model several living systems. While studies have typically focused on particles with isotropic tumbles, recent examples exhibit "tumble-turns", in which particles undergo 90° tumbles and so possess explicitly anisotropic dynamics. We study the consequences of such tumble-turn anisotropicity at both short and long-time scales. We model run-and-tumble-turn parti…
▽ More
Run-and-tumble processes successfully model several living systems. While studies have typically focused on particles with isotropic tumbles, recent examples exhibit "tumble-turns", in which particles undergo 90° tumbles and so possess explicitly anisotropic dynamics. We study the consequences of such tumble-turn anisotropicity at both short and long-time scales. We model run-and-tumble-turn particles as self-propelled particles subjected to an angular potential. Using agent-based simulations, we study the interplay of noise and potential on the particles' trajectories, demonstrating that the long-time effect is to alter the tumble-turn time, which governs the long-time dynamics. In particular, when normalized by this timescale, trajectories become independent of the underlying details of the potential. As such, we develop a simplified continuum theory, which quantitatively agrees with agent simulations. The hydrodynamic limit reveals that the transition to diffusive dynamics precedes the transition to isotropic dynamics as the hydrodynamic limit, while purely diffusive, can drive anisotropicity at intermediate times.
△ Less
Submitted 5 May, 2023;
originally announced May 2023.
-
Imaging the breakdown and restoration of topological protection in magnetic topological insulator MnBi$_2$Te$_4$
Authors:
Qile Li,
Iolanda Di Bernardo,
Johnathon Maniatis,
Daniel McEwen,
Liam Watson,
Benjamin Lowe,
Thi-Hai-Yen Vu,
Chi Xuan Trang,
Jinwoong Hwang,
Sung-Kwan Mo,
Michael S. Fuhrer,
Mark T. Edmonds
Abstract:
Quantum anomalous Hall (QAH) insulators transport charge without resistance along topologically protected chiral one-dimensional edge states. Yet, in magnetic topological insulators (MTI) to date, topological protection is far from robust, with the zero-magnetic field QAH effect only realised at temperatures an order of magnitude below the Néel temperature TN, though small magnetic fields can stab…
▽ More
Quantum anomalous Hall (QAH) insulators transport charge without resistance along topologically protected chiral one-dimensional edge states. Yet, in magnetic topological insulators (MTI) to date, topological protection is far from robust, with the zero-magnetic field QAH effect only realised at temperatures an order of magnitude below the Néel temperature TN, though small magnetic fields can stabilize QAH effect. Understanding why topological protection breaks down is therefore essential to realising QAH effect at higher temperatures. Here we use a scanning tunnelling microscope to directly map the size of the exchange gap (Eg,ex) and its spatial fluctuation in the QAH insulator 5-layer MnBi$_2$Te$_4$. We observe long-range fluctuations of Eg,ex with values ranging between 0 (gapless) and 70 meV, uncorrelated to individual point defects. We directly image the breakdown of topological protection, showing that the chiral edge state, the hallmark signature of a QAH insulator, hybridizes with extended gapless metallic regions in the bulk. Finally, we unambiguously demonstrate that the gapless regions originate in magnetic disorder, by demonstrating that a small magnetic field restores Eg,ex in these regions, explaining the recovery of topological protection in magnetic fields. Our results indicate that overcoming magnetic disorder is key to exploiting the unique properties of QAH insulators.
△ Less
Submitted 16 January, 2023;
originally announced January 2023.
-
Yield Stress and Compliance in Active Cell Monolayers
Authors:
Austin Hopkins,
Michael Chiang,
Benjamin Loewe,
Davide Marenduzzo,
M. Cristina Marchetti
Abstract:
The rheology of biological tissue plays an important role in many processes, from organ formation to cancer invasion. Here, we use a multi-phase field model of motile cells to simulate active microrheology within a tissue monolayer. When unperturbed, the tissue exhibits a transition between a solid-like state and a fluid-like state tuned by cell motility and deformability - the ratio of the energe…
▽ More
The rheology of biological tissue plays an important role in many processes, from organ formation to cancer invasion. Here, we use a multi-phase field model of motile cells to simulate active microrheology within a tissue monolayer. When unperturbed, the tissue exhibits a transition between a solid-like state and a fluid-like state tuned by cell motility and deformability - the ratio of the energetic costs of steric cell-cell repulsion and cell surface tension. When perturbed, solid tissues exhibit yield-stress behavior, with a threshold force for the onset of motion of a probe particle that vanishes upon approaching the solid-to-liquid transition. This onset of motion is qualitatively different in the low and high deformability regimes. At high deformability, the tissue is amorphous when solid, it responds compliantly to deformations, and the probe transition to motion is smooth. At low deformability, the monolayer is more ordered translationally and stiffer, and the onset of motion appears discontinuous. Our results suggest that cellular or nanoparticle transport in different types of tissues can be fundamentally different, and point to ways in which it can be controlled.
△ Less
Submitted 3 March, 2022;
originally announced March 2022.
-
Helical Flow States in Active Nematics
Authors:
Ryan Keogh,
Santhan Chandragiri,
Benjamin Loewe,
Tapio Ala-Nissila,
Sumesh Thampi,
Tyler N. Shendruk
Abstract:
We show that confining extensile nematics in 3D channels leads to the emergence of two self-organized flow states with nonzero helicity. The first is a pair of braided anti-parallel streams - this double helix occurs when the activity is moderate, anchoring negligible and reduced temperature high. The second consists of axially aligned counter-rotating vortices - this grinder train arises between…
▽ More
We show that confining extensile nematics in 3D channels leads to the emergence of two self-organized flow states with nonzero helicity. The first is a pair of braided anti-parallel streams - this double helix occurs when the activity is moderate, anchoring negligible and reduced temperature high. The second consists of axially aligned counter-rotating vortices - this grinder train arises between spontaneous axial streaming and the vortex lattice. These two unanticipated helical flow states illustrate the potential of active fluids to break symmetries and form complex but organized spatio-temporal structures in 3D fluidic devices.
△ Less
Submitted 17 June, 2022; v1 submitted 2 December, 2021;
originally announced December 2021.
-
Passive Janus Particles Are Self-propelled in Active Nematics
Authors:
Benjamin Loewe,
Tyler N. Shendruk
Abstract:
While active systems possess notable potential to form the foundation of new classes of autonomous materials, designing systems that can extract functional work from active surroundings has proven challenging. In this work, we extend these efforts to the realm of designed active liquid crystal/colloidal composites. We propose suspending colloidal particles with Janus anchoring conditions in an act…
▽ More
While active systems possess notable potential to form the foundation of new classes of autonomous materials, designing systems that can extract functional work from active surroundings has proven challenging. In this work, we extend these efforts to the realm of designed active liquid crystal/colloidal composites. We propose suspending colloidal particles with Janus anchoring conditions in an active nematic medium. These passive Janus particles become effectively self-propelled once immersed into an active nematic bath. The self-propulsion of passive Janus particles arises from the effective $+1/2$ topological charge their surface enforces on the surrounding active fluid. We analytically study their dynamics and the orientational dependence on the position of a companion $-1/2$ defect. We predict that at sufficiently small activity, the colloid and companion defect remain bound to each other, with the defect strongly orienting the colloid to propel either parallel or perpendicular to the nematic. At sufficiently high activity, we predict an unbinding of the colloid/defect pair. This work demonstrates how suspending engineered colloids in active liquid crystals may present a path to extracting activity to drive functionality.
△ Less
Submitted 21 January, 2022; v1 submitted 14 November, 2021;
originally announced November 2021.
-
Manifestation of strongly correlated electrons in a 2D kagome metal-organic framework
Authors:
Dhaneesh Kumar,
Jack Hellerstedt,
Bernard Field,
Benjamin Lowe,
Yuefeng Yin,
Nikhil V. Medhekar,
Agustin Schiffrin
Abstract:
The kagome lattice, whose electronic valence band (VB) structure includes two Dirac bands and one flat band, offers a rich space to realise tuneable topological and strongly correlated electronic phases in two-dimensional (2D) and layered materials. While strong electron correlations have been observed in inorganic kagome crystals, they remain elusive in organic systems. The latter offer versatile…
▽ More
The kagome lattice, whose electronic valence band (VB) structure includes two Dirac bands and one flat band, offers a rich space to realise tuneable topological and strongly correlated electronic phases in two-dimensional (2D) and layered materials. While strong electron correlations have been observed in inorganic kagome crystals, they remain elusive in organic systems. The latter offer versatile synthesis protocols via molecular self-assembly and metal-ligand coordination. Here, we report the synthesis of a 2D metal-organic framework (MOF) where di-cyano-anthracene (DCA) molecules form a kagome arrangement via coordination with copper (Cu) atoms on a silver surface [Ag(111)]. Low-temperature scanning tunnelling spectroscopy revealed Kondo screening--by the Ag(111) conduction electrons--of magnetic moments localised at Cu and DCA sites of the MOF. Density functional theory and mean-field Hubbard modelling show that these local moments emerge from strong interactions between kagome MOF electrons, enabling organic 2D materials as viable platforms for strongly correlated nanoelectronics and spintronics.
△ Less
Submitted 23 April, 2021;
originally announced April 2021.
-
Shape and size changes of adherent elastic epithelia
Authors:
Benjamin Loewe,
Francesco Serafin,
Suraj Shankar,
Mark J. Bowick,
M. Cristina Marchetti
Abstract:
Epithelial tissues play a fundamental role in various morphogenetic events during development and early embryogenesis. Although epithelial monolayers are often modeled as two-dimensional (2D) elastic surfaces, they distinguish themselves from conventional thin elastic plates in three important ways - the presence of an apical-basal polarity, spatial variability of cellular thickness, and their non…
▽ More
Epithelial tissues play a fundamental role in various morphogenetic events during development and early embryogenesis. Although epithelial monolayers are often modeled as two-dimensional (2D) elastic surfaces, they distinguish themselves from conventional thin elastic plates in three important ways - the presence of an apical-basal polarity, spatial variability of cellular thickness, and their nonequilibrium active nature. Here, we develop a minimal continuum model of a planar epithelial tissue as an active elastic material that incorporates all these features. We start from a full three-dimensional (3D)description of the tissue and derive an effective 2D model that captures, through the curvature of the apical surface, both the apical-basal asymmetry and the spatial geometry of the tissue. Crucially, variations of active stresses across the apical-basal axis lead to active torques that can drive curvature transitions. By identifying four distinct sources of activity, we find that bulk active stresses arising from actomyosin contractility and growth compete with boundary active tensions due to localized actomyosin cables and lamellipodial activity to generate the various states spanning the morphospace of a planar epithelium. Our treatment hence unifies 3D shape deformations through the coupled mechanics of apical curvature change and in-plane expansion/contraction of substrate-adhered tissues. Finally, we discuss the implications of our results for some biologically relevant processes such as tissue folding at the onset of lumen formation.
△ Less
Submitted 11 June, 2020; v1 submitted 10 February, 2020;
originally announced February 2020.
-
Solid-Liquid Transition of Deformable and Overlapping Active Particles
Authors:
Benjamin Loewe,
Michael Chiang,
Davide Marenduzzo,
M. Cristina Marchetti
Abstract:
Experiments and theory have shown that cell monolayers and epithelial tissues exhibit solid-liquid and glass-liquid transitions. These transitions are biologically relevant to our understanding of embryonic development, wound healing, and cancer. Current models typically consider purely two-dimensional monolayers with no overlaps between neighboring cells. In reality, overlaps are important, and t…
▽ More
Experiments and theory have shown that cell monolayers and epithelial tissues exhibit solid-liquid and glass-liquid transitions. These transitions are biologically relevant to our understanding of embryonic development, wound healing, and cancer. Current models typically consider purely two-dimensional monolayers with no overlaps between neighboring cells. In reality, overlaps are important, and they may be precursors of cell extrusion -- a key biophysical process to maintain homeostasis in epithelial tissues. Here, we use a multi-phase field model to study the solid-liquid transition in a confluent monolayer of deformable cells which can overlap. When cells overlap rather than deform, we find that the melting transition changes from continuous to discontinuous, and that there is an intermittent regime close to the transition, where solid and liquid states alternate over time. By studying the dynamics of $5$- and $7$-fold disclinations in the hexagonal lattice formed by the cell centers, we observe that these correlate with spatial fluctuations in the cellular overlap, and that cell extrusion tends to initiate near $5$-fold disclinations.
△ Less
Submitted 19 May, 2020; v1 submitted 22 December, 2019;
originally announced December 2019.
-
Flocking from a quantum analogy: Spin-orbit coupling in an active fluid
Authors:
Benjamin Loewe,
Anton Souslov,
Paul M. Goldbart
Abstract:
Systems composed of strongly interacting self-propelled particles can form a spontaneously flowing polar active fluid. The study of the connection between the microscopic dynamics of a single such particle and the macroscopic dynamics of the fluid can yield insights into experimentally realizable active flows, but this connection is well understood in only a few select cases. We introduce a model…
▽ More
Systems composed of strongly interacting self-propelled particles can form a spontaneously flowing polar active fluid. The study of the connection between the microscopic dynamics of a single such particle and the macroscopic dynamics of the fluid can yield insights into experimentally realizable active flows, but this connection is well understood in only a few select cases. We introduce a model of self-propelled particles based on an analogy with the motion of electrons that have strong spin-orbit coupling. We find that, within our model, self-propelled particles are subject to an analog of the Heisenberg uncertainty principle that relates translational and rotational noise. Furthermore, by coarse-graining this microscopic model, we establish expressions for the coefficients of the Toner-Tu equations---the hydrodynamic equations that describe an active fluid composed of these "active spins." The connection between self-propelled particles and quantum spins may help realize exotic phases of matter using active fluids via analogies with systems composed of strongly correlated electrons.
△ Less
Submitted 2 August, 2017;
originally announced August 2017.
-
Emergent tilt order in Dirac polymer liquids
Authors:
Anton Souslov,
Benjamin Loewe,
Paul M. Goldbart
Abstract:
We study a liquid of zigzagging two-dimensional directed polymers with bending rigidity, i.e., polymers whose conformations follow checkerboard paths. In the continuum limit the statistics of such polymers obey the Dirac equation for particles of imaginary mass. We exploit this observation to investigate a liquid of these polymers via a quantum many-fermion analogy. A self-consistent approximation…
▽ More
We study a liquid of zigzagging two-dimensional directed polymers with bending rigidity, i.e., polymers whose conformations follow checkerboard paths. In the continuum limit the statistics of such polymers obey the Dirac equation for particles of imaginary mass. We exploit this observation to investigate a liquid of these polymers via a quantum many-fermion analogy. A self-consistent approximation predicts a phase of tilted order, in which the polymers may develop a preference to zig rather than zag. We compute the phase diagram and key response functions for the polymer liquid, and comment on the role played by fluctuations.
△ Less
Submitted 3 October, 2015; v1 submitted 23 December, 2014;
originally announced December 2014.
-
Identification of vacancy defects in a thin film perovskite oxide
Authors:
D. J. Keeble,
R. A. Mackie,
W. Egger,
B. Lowe,
P. Pikart,
C. Hugenschmidt,
T. J. Jackson
Abstract:
Vacancy defects in thin film laser ablated SrTiO3 on SrTiO3 were identified using variable energy positron annihilation lifetime measurements. Strontium vacancy related defects were the dominant positron traps and, apart from in the top ~ 50 nm, were found to be uniformly distributed. The surface layer showed an increase in annihilation from larger open-volume defects, large vacancy clusters or…
▽ More
Vacancy defects in thin film laser ablated SrTiO3 on SrTiO3 were identified using variable energy positron annihilation lifetime measurements. Strontium vacancy related defects were the dominant positron traps and, apart from in the top ~ 50 nm, were found to be uniformly distributed. The surface layer showed an increase in annihilation from larger open-volume defects, large vacancy clusters or nanovoids.
△ Less
Submitted 30 June, 2009;
originally announced June 2009.