1
|
Álvarez Domenech I, Rodríguez-Laguna J, Cuerno R, Córdoba-Torres P, Santalla SN. Shape effects in the fluctuations of random isochrones on a square lattice. Phys Rev E 2024; 109:034104. [PMID: 38632797 DOI: 10.1103/physreve.109.034104] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 01/19/2024] [Indexed: 04/19/2024]
Abstract
We consider the isochrone curves in first-passage percolation on a 2D square lattice, i.e., the boundary of the set of points which can be reached in less than a given time from a certain origin. The occurrence of an instantaneous average shape is described in terms of its Fourier components, highlighting a crossover between a diamond and a circular geometry as the noise level is increased. Generally, these isochrones can be understood as fluctuating interfaces with an inhomogeneous local width which reveals the underlying lattice structure. We show that once these inhomogeneities have been taken into account, the fluctuations fall into the Kardar-Parisi-Zhang universality class with very good accuracy, where they reproduce the Family-Vicsek Ansatz with the expected exponents and the Tracy-Widom histogram for the local radial fluctuations.
Collapse
Affiliation(s)
- Iván Álvarez Domenech
- Dto. Física Matemática y de Fluidos, Universidad Nacional de Educación a Distancia, Las Rozas 28232, Spain
| | | | - Rodolfo Cuerno
- Dto. Matemáticas & GISC, Universidad Carlos III de Madrid, Leganés 28911, Spain
| | - Pedro Córdoba-Torres
- Dto. Física Matemática y de Fluidos, Universidad Nacional de Educación a Distancia, Las Rozas 28232, Spain
| | - Silvia N Santalla
- Dto. Física & GISC, Universidad Carlos III de Madrid, Leganés 28911, Spain
| |
Collapse
|
2
|
Santalla SN, Domenech IÁ, Villarrubia D, Cuerno R, Rodríguez-Laguna J. Universal fluctuations of global geometrical measurements in planar clusters. Phys Rev E 2024; 109:034127. [PMID: 38632763 DOI: 10.1103/physreve.109.034127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2023] [Accepted: 02/26/2024] [Indexed: 04/19/2024]
Abstract
We characterize universal features of the sample-to-sample fluctuations of global geometrical observables, such as the area, width, length, and center-of-mass position, in random growing planar clusters. Our examples are taken from simulations of both continuous and discrete models of kinetically rough interfaces, including several universality classes, such as Kardar-Parisi-Zhang. We mostly focus on the scaling behavior with time of the sample-to-sample deviation for those global magnitudes, but we have also characterized their histograms and correlations.
Collapse
Affiliation(s)
- Silvia N Santalla
- Dto. Física & GISC, Universidad Carlos III de Madrid, Leganés 28911, Spain
| | - Ivan Álvarez Domenech
- Dto. Física Matemática y de Fluidos, Universidad Nacional de Educación a Distancia, Las Rozas 28232, Spain
| | - Daniel Villarrubia
- Dto. Matemáticas & GISC, Universidad Carlos III de Madrid, Leganés 28911, Spain
| | - Rodolfo Cuerno
- Dto. Matemáticas & GISC, Universidad Carlos III de Madrid, Leganés 28911, Spain
| | | |
Collapse
|
3
|
Chan ST, Fried E. Marangoni spreading on liquid substrates in new media art. PNAS NEXUS 2024; 3:pgae059. [PMID: 38725527 PMCID: PMC11079615 DOI: 10.1093/pnasnexus/pgae059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/12/2023] [Accepted: 01/31/2024] [Indexed: 05/12/2024]
Abstract
With the advent of new media art, artists have harnessed fluid dynamics to create captivating visual narratives. A striking technique known as dendritic painting employs mixtures of ink and isopropanol atop paint, yielding intricate tree-like patterns. To unravel the intricacies of that technique, we examine the spread of ink/alcohol droplets over liquid substrates with diverse rheological properties. On Newtonian substrates, the droplet size evolution exhibits two power laws, suggesting an underlying interplay between viscous and Marangoni forces. The leading edge of the droplet spreads as a precursor film with an exponent of 3/8, while its main body spreads with an exponent of 1/4. For a weakly shear-thinning acrylic resin substrate, the same power laws persist, but dendritic structures emerge, and the texture of the precursor film roughens. The observed roughness and growth exponents (3/4 and 3/5) suggest a connection to the quenched Kardar-Parisi-Zhang universality class, hinting at the existence of quenched disorder in the liquid substrate. Mixing the resin with acrylic paint renders it more viscous and shear-thinning, refining the dendrite edges and further roughening the precursor film. At larger paint concentrations, the substrate becomes a power-law fluid. The roughness and growth exponents then approach 1/2 and 3/4, respectively, deviating from known universality classes. The ensuing structures have a fractal dimension of 1.68, characteristic of diffusion-limited aggregation. These findings underscore how the nonlinear rheological properties of the liquid substrate, coupled with the Laplacian nature of Marangoni spreading, can overshadow the local kinetic roughening of the droplet interface.
Collapse
Affiliation(s)
- San To Chan
- Mechanics and Materials Unit, Okinawa Institute of Science and Technology Graduate University, Onna, Okinawa 904-0495, Japan
| | - Eliot Fried
- Mechanics and Materials Unit, Okinawa Institute of Science and Technology Graduate University, Onna, Okinawa 904-0495, Japan
| |
Collapse
|
4
|
Mazarei M, Åström J, Westerholm J, Karttunen M. Effect of substrate heterogeneity and topology on epithelial tissue growth dynamics. Phys Rev E 2023; 108:054405. [PMID: 38115499 DOI: 10.1103/physreve.108.054405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Accepted: 09/20/2023] [Indexed: 12/21/2023]
Abstract
Tissue growth kinetics and interface dynamics depend on the properties of the tissue environment and cell-cell interactions. In cellular environments, substrate heterogeneity and geometry arise from a variety factors, such as the structure of the extracellular matrix and nutrient concentration. We used the CellSim3D model, a kinetic cell division simulator, to investigate the growth kinetics and interface roughness dynamics of epithelial tissue growth on heterogeneous substrates with varying topologies. The results show that the presence of quenched disorder has a clear effect on the colony morphology and the roughness scaling of the interface in the moving interface regime. In a medium with quenched disorder, the tissue interface has a smaller interface roughness exponent, α, and a larger growth exponent, β. The scaling exponents also depend on the topology of the substrate and cannot be categorized by well-known universality classes.
Collapse
Affiliation(s)
- Mahmood Mazarei
- Department of Physics and Astronomy, Western University, 1151 Richmond Street, London, Ontario, Canada N6A 3K7
| | - Jan Åström
- CSC Scientific Computing Ltd, Kägelstranden 14, 02150 Esbo, Finland
| | - Jan Westerholm
- Faculty of Science and Engineering, Åbo Akademi University, Vattenborgsvägen 3, FI-20500 Åbo, Finland
| | - Mikko Karttunen
- Department of Physics and Astronomy, Western University, 1151 Richmond Street, London, Ontario, Canada N6A 3K7
- Department of Chemistry, Western University, 1151 Richmond Street, London, Ontario, Canada N6A 5B7
| |
Collapse
|
5
|
Mazarei M, Åström J, Westerholm J, Karttunen M. In silico testing of the universality of epithelial tissue growth. Phys Rev E 2022; 106:L062402. [PMID: 36671099 DOI: 10.1103/physreve.106.l062402] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Accepted: 11/22/2022] [Indexed: 12/14/2022]
Abstract
The universality of interfacial roughness in growing epithelial tissue has remained a controversial issue. Kardar-Parisi-Zhang (KPZ) and molecular beam epitaxy (MBE) universality classes have been reported among other behaviors including a total lack of universality. Here, we simulate tissues using the cellsim3d kinetic division model for deformable cells to investigate cell-colony scaling. With seemingly minor model changes, it can reproduce both KPZ- and MBE-like scaling in configurations that mimic the respective experiments. Tissue growth with strong cell-cell adhesion in a linear geometry is KPZ like, while weakly adhesive tissues in a radial geometry are MBE like. This result neutralizes the apparent scaling controversy.
Collapse
Affiliation(s)
- Mahmood Mazarei
- Department of Physics and Astronomy, Western University, 1151 Richmond Street, London, Ontario, Canada N6A 3K7
| | - Jan Åström
- CSC Scientific Computing Limited, Kägelstranden 14, FI-02150 Esbo, Finland
| | - Jan Westerholm
- Faculty of Science and Engineering, Åbo Akademi University, Vattenborgsvägen 3, FI-20500 Åbo, Finland
| | - Mikko Karttunen
- Department of Physics and Astronomy, Western University, 1151 Richmond Street, London, Ontario, Canada N6A 3K7.,Department of Chemistry, Western University, 1151 Richmond Street, London, Ontario, Canada N6A 5B7
| |
Collapse
|
6
|
Abstract
The morphogenesis of two-dimensional bacterial colonies has been well studied. However, little is known about the colony morphologies of bacteria growing in three dimensions, despite the prevalence of three-dimensional environments (e.g., soil, inside hosts) as natural bacterial habitats. Using experiments on bacteria in granular hydrogel matrices, we find that dense multicellular colonies growing in three dimensions undergo a common morphological instability and roughen, adopting a characteristic broccoli-like morphology when they exceed a critical size. Analysis of a continuum “active fluid” model of the expanding colony reveals that this behavior originates from an interplay of competition for nutrients with growth-driven colony expansion, both of which vary spatially. These results shed light on the fundamental biophysical principles underlying growth in three dimensions. How do growing bacterial colonies get their shapes? While colony morphogenesis is well studied in two dimensions, many bacteria grow as large colonies in three-dimensional (3D) environments, such as gels and tissues in the body or subsurface soils and sediments. Here, we describe the morphodynamics of large colonies of bacteria growing in three dimensions. Using experiments in transparent 3D granular hydrogel matrices, we show that dense colonies of four different species of bacteria generically become morphologically unstable and roughen as they consume nutrients and grow beyond a critical size—eventually adopting a characteristic branched, broccoli-like morphology independent of variations in the cell type and environmental conditions. This behavior reflects a key difference between two-dimensional (2D) and 3D colonies; while a 2D colony may access the nutrients needed for growth from the third dimension, a 3D colony inevitably becomes nutrient limited in its interior, driving a transition to unstable growth at its surface. We elucidate the onset of the instability using linear stability analysis and numerical simulations of a continuum model that treats the colony as an “active fluid” whose dynamics are driven by nutrient-dependent cellular growth. We find that when all dimensions of the colony substantially exceed the nutrient penetration length, nutrient-limited growth drives a 3D morphological instability that recapitulates essential features of the experimental observations. Our work thus provides a framework to predict and control the organization of growing colonies—as well as other forms of growing active matter, such as tumors and engineered living materials—in 3D environments.
Collapse
|
7
|
Barreales BG, Meléndez JJ, Cuerno R, Ruiz-Lorenzo JJ. Large-scale kinetic roughening behavior of coffee-ring fronts. Phys Rev E 2022; 106:044801. [PMID: 36397471 DOI: 10.1103/physreve.106.044801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Accepted: 09/08/2022] [Indexed: 06/16/2023]
Abstract
We have studied the kinetic roughening behavior of the fronts of coffee-ring aggregates via extensive numerical simulations of the off-lattice model considered for this context [Dias et al., Soft Matter 14, 1903 (2018)1744-683X10.1039/C7SM02136D]. This model describes ballistic aggregation of patchy colloids and depends on a parameter r_{AB} which controls the affinity of the two patches, A and B. Suitable boundary conditions allow us to elucidate a discontinuous pinning-depinning transition at r_{AB}=0, with the front displaying intrinsic anomalous scaling, but with unusual exponent values α≃1.2, α_{loc}≃0.5, β≃1, and z≃1.2. For 0 0.01 and the system suffers a strong crossover dominated by the r_{AB}=0 behavior for r_{AB}≤0.01. A detailed analysis of correlation functions shows that the aggregate fronts are always in the moving phase for 0
Collapse
Affiliation(s)
- B G Barreales
- Departamento de Física, Universidad de Extremadura, 06006 Badajoz, Spain
- Instituto de Computación Científica Avanzada de Extremadura (ICCAEx), Universidad de Extremadura, 06006 Badajoz, Spain
| | - J J Meléndez
- Departamento de Física, Universidad de Extremadura, 06006 Badajoz, Spain
- Instituto de Computación Científica Avanzada de Extremadura (ICCAEx), Universidad de Extremadura, 06006 Badajoz, Spain
| | - R Cuerno
- Departamento de Matemáticas and Grupo Interdisciplinar de Sistemas Complejos (GISC), Universidad Carlos III de Madrid, 28911 Leganés, Spain
| | - J J Ruiz-Lorenzo
- Departamento de Física, Universidad de Extremadura, 06006 Badajoz, Spain
- Instituto de Computación Científica Avanzada de Extremadura (ICCAEx), Universidad de Extremadura, 06006 Badajoz, Spain
- Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), 50018 Zaragoza, Spain
| |
Collapse
|
8
|
Oguma T, Takigawa-Imamura H, Miura T. Mechanism underlying dynamic scaling properties observed in the contour of spreading epithelial monolayer. Phys Rev E 2020; 102:062408. [PMID: 33466041 DOI: 10.1103/physreve.102.062408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Accepted: 10/14/2020] [Indexed: 11/07/2022]
Abstract
We found evidence of dynamic scaling in the spreading of Madin-Darby canine kidney (MDCK) cell monolayer, which can be characterized by the Hurst exponent α=0.86 and the growth exponent β=0.73, and theoretically and experimentally clarified the mechanism that governs the contour shape dynamics. Dynamic scaling refers to the roughness of the surface scales, both spatially and temporally. During the spreading of the monolayer, it is known that so-called leader cells generate the driving force and lead the other cells. Our time-lapse observations of cell behavior showed that these leader cells appeared at the early stage of the spreading and formed the monolayer protrusion. Informed by these observations, we developed a simple mathematical model that included differences in cell motility, cell-cell adhesion, and random cell movement. The model reproduced the quantitative characteristics obtained from the experiment, such as the spreading speed, the distribution of the increment, and the dynamic scaling law. Analysis of the model equation shows that the model can reproduce different scaling laws from (α=0.5,β=0.25) to (α=0.9,β=0.75), where the exponents α and β are determined by two dimensionless quantities determined by the microscopic cell behavior. From the analytical result, parameter estimation from the experimental results was achieved. The monolayer on the collagen-coated dishes showed a different scaling law, α=0.74,β=0.68, suggesting that cell motility increased ninefold. This result was consistent with the assay of the single-cell motility. Our study demonstrated that the dynamics of the contour of the monolayer were explained by the simple model, and we propose a mechanism that exhibits the dynamic scaling property.
Collapse
Affiliation(s)
- Toshiki Oguma
- Department of Anatomy and Cell Biology, Graduate School of Medical Sciences, Kyushu University, Japan
| | - Hisako Takigawa-Imamura
- Department of Anatomy and Cell Biology, Graduate School of Medical Sciences, Kyushu University, Japan
| | - Takashi Miura
- Department of Anatomy and Cell Biology, Graduate School of Medical Sciences, Kyushu University, Japan
| |
Collapse
|
9
|
Santalla SN, Koroutchev K, Korutcheva E, Rodríguez-Laguna J. Power accretion in social systems. Phys Rev E 2019; 100:012143. [PMID: 31499905 DOI: 10.1103/physreve.100.012143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2019] [Indexed: 11/07/2022]
Abstract
We consider a model of power distribution in a social system where a set of agents plays a simple game on a graph: The probability of winning each round is proportional to the agent's current power, and the winner gets more power as a result. We show that when the agents are distributed on simple one-dimensional and two-dimensional networks, inequality grows naturally up to a certain stationary value characterized by a clear division between a higher and a lower class of agents. High class agents are separated by one or several lower class agents which serve as a geometrical barrier preventing further flow of power between them. Moreover, we consider the effect of redistributive mechanisms, such as proportional (nonprogressive) taxation. Sufficient taxation will induce a sharp transition towards a more equal society, and we argue that the critical taxation level is uniquely determined by the system geometry. Interestingly, we find that the roughness and Shannon entropy of the power distributions are a very useful complement to the standard measures of inequality, such as the Gini index and the Lorenz curve.
Collapse
Affiliation(s)
- Silvia N Santalla
- Departamento de Física and Grupo Interdisciplinar de Sistemas Complejos, Universidad Carlos III de Madrid, E-28911 Madrid, Spain
| | - Kostadin Koroutchev
- Escuela Politécnica Superior, Universidad Autónoma de Madrid, E-28049 Cantoblanco, Spain
| | - Elka Korutcheva
- Departamento de Física Fundamental, Universidad Nacional de Educación a Distancia (UNED), E-28040 Spain.,Department of Theoretical Physics, G. Nadjakov Institute of Solid State Physics, Bulgarian Academy of Sciences, 72 Tzarigradsko Shaussee Boulevard, B-1784 Sofia, Bulgaria
| | - Javier Rodríguez-Laguna
- Departamento de Física Fundamental, Universidad Nacional de Educación a Distancia (UNED), E-28040 Spain
| |
Collapse
|
10
|
Santalla SN, Ferreira SC. Eden model with nonlocal growth rules and kinetic roughening in biological systems. Phys Rev E 2018; 98:022405. [PMID: 30253509 DOI: 10.1103/physreve.98.022405] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2018] [Indexed: 11/07/2022]
Abstract
We investigate an off-lattice Eden model where the growth of new cells is performed with a probability dependent on the availability of resources coming externally towards the growing aggregate. The concentration of nutrients necessary for replication is assumed to be proportional to the voids connecting the replicating cells to the outer region, introducing therefore a nonlocal dependence on the replication rule. Our simulations point out that the Kadar-Parisi-Zhang (KPZ) universality class is a transient that can last for long periods in plentiful environments. For conditions of nutrient scarcity, we observe a crossover from regular KPZ to unstable growth, passing by a transient consistent with the quenched KPZ class at the pinning transition. Our analysis sheds light on results reporting on the universality class of kinetic roughening in akin experiments of biological growth.
Collapse
Affiliation(s)
- Silvia N Santalla
- Departamento de Física and Grupo Interdisciplinar de Sistemas Complejos, Universidad Carlos III de Madrid, Leganés, Spain
| | - Silvio C Ferreira
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, Minas Gerais, Brazil.,National Institute of Science and Technology for Complex Systems, Rio de Janeiro, Brazil
| |
Collapse
|