1
|
Grimm RT, Eaves JD. Direct Numerical Solutions to Stochastic Differential Equations with Multiplicative Noise. PHYSICAL REVIEW LETTERS 2024; 132:267101. [PMID: 38996310 DOI: 10.1103/physrevlett.132.267101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Accepted: 03/25/2024] [Indexed: 07/14/2024]
Abstract
Inspired by path integral solutions to the quantum relaxation problem, we develop a numerical method to solve classical stochastic differential equations with multiplicative noise that avoids averaging over trajectories. To test the method, we simulate the dynamics of a classical oscillator multiplicatively coupled to non-Markovian noise. When accelerated using tensor factorization techniques, it accurately estimates the transition into the bifurcation regime of the oscillator and outperforms trajectory-averaging simulations with a computational cost that is orders of magnitude lower.
Collapse
Affiliation(s)
- Ryan T Grimm
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Joel D Eaves
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| |
Collapse
|
2
|
Omar AK, Klymko K, GrandPre T, Geissler PL, Brady JF. Tuning nonequilibrium phase transitions with inertia. J Chem Phys 2023; 158:074904. [PMID: 36813709 DOI: 10.1063/5.0138256] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Abstract
In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation-dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. This density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.
Collapse
Affiliation(s)
- Ahmad K Omar
- Department of Materials Science and Engineering, University of California, Berkeley, California 94720, USA
| | - Katherine Klymko
- NERSC, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Trevor GrandPre
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Phillip L Geissler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John F Brady
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
3
|
Abstract
The human genome is arranged in the cell nucleus nonrandomly, and phase separation has been proposed as an important driving force for genome organization. However, the cell nucleus is an active system, and the contribution of nonequilibrium activities to phase separation and genome structure and dynamics remains to be explored. We simulated the genome using an energy function parametrized with chromosome conformation capture (Hi-C) data with the presence of active, nondirectional forces that break the detailed balance. We found that active forces that may arise from transcription and chromatin remodeling can dramatically impact the spatial localization of heterochromatin. When applied to euchromatin, active forces can drive heterochromatin to the nuclear envelope and compete with passive interactions among heterochromatin that tend to pull them in opposite directions. Furthermore, active forces induce long-range spatial correlations among genomic loci beyond single chromosome territories. We further showed that the impact of active forces could be understood from the effective temperature defined as the fluctuation-dissipation ratio. Our study suggests that nonequilibrium activities can significantly impact genome structure and dynamics, producing unexpected collective phenomena.
Collapse
Affiliation(s)
- Zhongling Jiang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| | - Yifeng Qi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| | - Kartik Kamat
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, United States
| |
Collapse
|
4
|
Ishimoto K, Moreau C, Yasuda K. Self-organized swimming with odd elasticity. Phys Rev E 2022; 105:064603. [PMID: 35854482 DOI: 10.1103/physreve.105.064603] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Accepted: 03/23/2022] [Indexed: 06/15/2023]
Abstract
We theoretically investigate self-oscillating waves of an active material, which were recently introduced as a nonsymmetric part of the elastic moduli, termed odd elasticity. Using Purcell's three-link swimmer model, we reveal that an odd-elastic filament at low Reynolds number can swim in a self-organized manner and that the time-periodic dynamics are characterized by a stable limit cycle generated by elastohydrodynamic interactions. Also, we consider a noisy shape gait and derive a swimming formula for a general elastic material in the Stokes regime with its elasticity modulus being represented by a nonsymmetric matrix, demonstrating that the odd elasticity produces biased net locomotion from random noise.
Collapse
Affiliation(s)
- Kenta Ishimoto
- Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606-8502, Japan
| | - Clément Moreau
- Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606-8502, Japan
| | - Kento Yasuda
- Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
5
|
O'Byrne J, Tailleur J. Lamellar to Micellar Phases and Beyond: When Tactic Active Systems Admit Free Energy Functionals. PHYSICAL REVIEW LETTERS 2020; 125:208003. [PMID: 33258650 DOI: 10.1103/physrevlett.125.208003] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Accepted: 10/14/2020] [Indexed: 06/12/2023]
Abstract
We consider microscopic models of active particles whose velocities, rotational diffusivities, and tumbling rates depend on the gradient of a local field that is either externally imposed or depends on all particle positions. Despite the fundamental differences between active and passive dynamics at the microscopic scale, we show that a large class of such tactic active systems admit fluctuating hydrodynamics equivalent to those of interacting Brownian colloids in equilibrium. We exploit this mapping to show how taxis may lead to the lamellar and micellar phases observed for soft repulsive colloids. In the context of chemotaxis, we show how the competition between chemoattractant and chemorepellent may lead to a bona fide equilibrium liquid-gas phase separation in which a loss of thermodynamic stability of the fluid signals the onset of a chemotactic collapse.
Collapse
Affiliation(s)
- J O'Byrne
- Université de Paris, Laboratoire Matière et Systèmes Complexes (MSC), UMR 7057 CNRS, F-75205 Paris, France
| | - J Tailleur
- Université de Paris, Laboratoire Matière et Systèmes Complexes (MSC), UMR 7057 CNRS, F-75205 Paris, France
| |
Collapse
|
6
|
Alvarado J, Sheinman M, Sharma A, MacKintosh FC, Koenderink GH. Force percolation of contractile active gels. SOFT MATTER 2017; 13:5624-5644. [PMID: 28812094 DOI: 10.1039/c7sm00834a] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Living systems provide a paradigmatic example of active soft matter. Cells and tissues comprise viscoelastic materials that exert forces and can actively change shape. This strikingly autonomous behavior is powered by the cytoskeleton, an active gel of semiflexible filaments, crosslinks, and molecular motors inside cells. Although individual motors are only a few nm in size and exert minute forces of a few pN, cells spatially integrate the activity of an ensemble of motors to produce larger contractile forces (∼nN and greater) on cellular, tissue, and organismal length scales. Here we review experimental and theoretical studies on contractile active gels composed of actin filaments and myosin motors. Unlike other active soft matter systems, which tend to form ordered patterns, actin-myosin systems exhibit a generic tendency to contract. Experimental studies of reconstituted actin-myosin model systems have long suggested that a mechanical interplay between motor activity and the network's connectivity governs this contractile behavior. Recent theoretical models indicate that this interplay can be understood in terms of percolation models, extended to include effects of motor activity on the network connectivity. Based on concepts from percolation theory, we propose a state diagram that unites a large body of experimental observations. This framework provides valuable insights into the mechanisms that drive cellular shape changes and also provides design principles for synthetic active materials.
Collapse
Affiliation(s)
- José Alvarado
- Systems Biophysics Department, AMOLF, 1098 XG Amsterdam, The Netherlands.
| | | | | | | | | |
Collapse
|
7
|
Dasbiswas K, Majkut S, Discher DE, Safran SA. Substrate stiffness-modulated registry phase correlations in cardiomyocytes map structural order to coherent beating. Nat Commun 2015; 6:6085. [PMID: 25597833 DOI: 10.1038/ncomms7085] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2014] [Accepted: 12/11/2014] [Indexed: 11/09/2022] Open
Abstract
Recent experiments show that both striation, an indication of the structural registry in muscle fibres, as well as the contractile strains produced by beating cardiac muscle cells can be optimized by substrate stiffness. Here we show theoretically how the substrate rigidity dependence of the registry data can be mapped onto that of the strain measurements. We express the elasticity-mediated structural registry as a phase-order parameter using a statistical physics approach that takes the noise and disorder inherent in biological systems into account. By assuming that structurally registered myofibrils also tend to beat in phase, we explain the observed dependence of both striation and strain measurements of cardiomyocytes on substrate stiffness in a unified manner. The agreement of our ideas with experiment suggests that the correlated beating of heart cells may be limited by the structural order of the myofibrils, which in turn is regulated by their elastic environment.
Collapse
Affiliation(s)
- K Dasbiswas
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovot 76100, Israel
| | - S Majkut
- 1] Department of Molecular and Biophysical Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA [2] Physics and Astronomy Graduate Group, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - D E Discher
- 1] Department of Molecular and Biophysical Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA [2] Physics and Astronomy Graduate Group, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA [3] Cell and Molecular Biology Graduate Group, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Samuel A Safran
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovot 76100, Israel
| |
Collapse
|
8
|
Muñoz JJ, Albo S. Physiology-based model of cell viscoelasticity. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:012708. [PMID: 23944493 DOI: 10.1103/physreve.88.012708] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2012] [Revised: 06/08/2013] [Indexed: 06/02/2023]
Abstract
The measured viscoelastic properties of biological tissues is the result of the passive and active response of the cells. We propose an evolution law of the remodeling process in the cytoskeleton which is able to mimic the viscous properties of biological cellular tissues. Our model is based on dynamical changes of the resting length. We show that under the small strain regime, the linear rheology models are recovered, with the relaxation time being replaced by the cell resistance to remodel. We implement the one-dimensional model into network systems of two and three dimensions, and show that the same conclusions may be drawn for those systems.
Collapse
Affiliation(s)
- José J Muñoz
- Department of Applied Mathematics III, Laboratori de Càlcul Numèric (LaCàN) and Universitat Politècnica de Catalunya (UPC), 08036 Barcelona, Spain.
| | | |
Collapse
|
9
|
Ben-Isaac E, Park Y, Popescu G, Brown FLH, Gov NS, Shokef Y. Effective temperature of red-blood-cell membrane fluctuations. PHYSICAL REVIEW LETTERS 2011; 106:238103. [PMID: 21770546 DOI: 10.1103/physrevlett.106.238103] [Citation(s) in RCA: 78] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2011] [Indexed: 05/31/2023]
Abstract
Biologically driven nonequilibrium fluctuations are often characterized by their non-Gaussianity or by an "effective temperature", which is frequency dependent and higher than the ambient temperature. We address these two measures theoretically by examining a randomly kicked particle, with a variable number of kicking motors, and show how these two indicators of nonequilibrium behavior can contradict. Our results are compared with new experiments on shape fluctuations of red-blood cell membranes, and demonstrate how the physical nature of the motors in this system can be revealed using these global measures of nonequilibrium.
Collapse
Affiliation(s)
- Eyal Ben-Isaac
- Department of Chemical Physics, The Weizmann Institute of Science, Rehovot, Israel
| | | | | | | | | | | |
Collapse
|
10
|
Morozov KI, Pismen LM. Cytoskeleton fluidization versus resolidification: prestress effect. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:051920. [PMID: 21728584 DOI: 10.1103/physreve.83.051920] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2010] [Revised: 02/15/2011] [Indexed: 05/31/2023]
Abstract
The differential elastic modulus of an active actomyosin network is computed as a function of applied stress, taking into account both thermal and motor contributions to filament compliance in the low-frequency domain. It is shown that, due to a dual nature of motor activity, increasing motor concentration may either stiffen the network due to stronger prestress or soften it due to motor agitation, in accordance with experimental data. Prestress anisotropy, which may be induced by redistribution of motors triggered by external force, causes anisotropy of the elastic moduli. This helps to explain the contradictory phenomena of cell fluidization and resolidification in response to transient stretch observed in recent experiments.
Collapse
Affiliation(s)
- Konstantin I Morozov
- Department of Chemical Engineering, Technion-Israel Institute of Technology, Haifa, Israel
| | | |
Collapse
|