-
Cell-level modelling of homeostasis in confined epithelial monolayers
Authors:
KVS Chaithanya,
Jan Rozman,
Andrej Košmrlj,
Rastko Sknepnek
Abstract:
Tissue homeostasis, the biological process of maintaining a steady state in tissue via control of cell proliferation, death, and metabolic function, is essential for the development, growth, maintenance, and proper function of living organisms. Disruptions to this process can lead to serious diseases and even death. In this study, we use the vertex model for the cell-level description of tissue me…
▽ More
Tissue homeostasis, the biological process of maintaining a steady state in tissue via control of cell proliferation, death, and metabolic function, is essential for the development, growth, maintenance, and proper function of living organisms. Disruptions to this process can lead to serious diseases and even death. In this study, we use the vertex model for the cell-level description of tissue mechanics to investigate the impact of the tissue microenvironment and local mechanical properties of cells on homeostasis in confined epithelial tissues. We find a dynamic steady state, where the balance between cell divisions and removals sustains homeostasis. By characterising homeostasis in terms of cell count, tissue area, and the cells' neighbour count distribution, we identify the factors that govern regulated and ordered tissue growth. This work, therefore, sheds light on the mechanisms underlying tissue homeostasis and highlights the importance of mechanics in the control of biological processes such as tissue development and disease pathology.
△ Less
Submitted 23 March, 2024;
originally announced March 2024.
-
From Dry to Wet Vertex Model Dynamics: Generating Sustained Flows
Authors:
Jan Rozman,
Chaithanya K. V. S.,
Julia M. Yeomans,
Rastko Sknepnek
Abstract:
Complex tissue flows in epithelia are driven by intra- and inter-cellular processes that generate, maintain, and coordinate mechanical forces. There has been growing evidence that cell shape anisotropy, manifested as nematic order, plays an important role in this process. Here we extend a nematic vertex model by replacing substrate friction with internal viscous dissipation, of relevance to epithe…
▽ More
Complex tissue flows in epithelia are driven by intra- and inter-cellular processes that generate, maintain, and coordinate mechanical forces. There has been growing evidence that cell shape anisotropy, manifested as nematic order, plays an important role in this process. Here we extend a nematic vertex model by replacing substrate friction with internal viscous dissipation, of relevance to epithelia not supported by a substrate or the extracellular matrix, such as many early-stage embryos. When coupled to cell shape anisotropy, the internal viscous dissipation allows for long-range velocity correlations and thus enables the spontaneous emergence of flows with a large degree of spatiotemporal organisation. We demonstrate sustained flow in epithelial sheets confined to a channel, thus providing a link between the dynamical behaviour of continuum active nematics and the cell-level vertex model of tissue dynamics.
△ Less
Submitted 18 December, 2023;
originally announced December 2023.
-
PyMembrane: A flexible framework for efficient simulations of elastic and liquid membranes
Authors:
D. A. Matoz-Fernandez,
Siyu Li,
Monica Olvera de la Cruz,
Rastko Sknepnek
Abstract:
PyMembrane is a software package for simulating liquid and elastic membranes using a discretisation of the continuum description based on unstructured triangulated two-dimensional meshes embedded in three-dimensional space. The package is written in C++, with a flexible and intuitive Python interface, allowing for a quick setup, execution and analysis of complex simulations. PyMembrane follows mod…
▽ More
PyMembrane is a software package for simulating liquid and elastic membranes using a discretisation of the continuum description based on unstructured triangulated two-dimensional meshes embedded in three-dimensional space. The package is written in C++, with a flexible and intuitive Python interface, allowing for a quick setup, execution and analysis of complex simulations. PyMembrane follows modern software engineering principles and features a modular design that allows for straightforward implementation of custom extensions while ensuring consistency and enabling inexpensive maintenance. A hallmark feature of this design is the use of a standardized C++ interface which streamlines adding new functionalities. Furthermore, PyMembrane uses data structures optimised for unstructured meshes, ensuring efficient mesh operations and force calculations. By providing several templates for typical simulations supplemented by extensive documentation, the users can seamlessly set up and run research-level simulations and extend the package to integrate additional features, underscoring PyMembrane's commitment to user-centric design.
△ Less
Submitted 24 August, 2023;
originally announced August 2023.
-
Shape-tension coupling produces nematic order in an epithelium vertex model
Authors:
Jan Rozman,
Rastko Sknepnek,
Julia M. Yeomans
Abstract:
We study the vertex model for epithelial tissue mechanics extended to include coupling between the cell shapes and tensions in cell-cell junctions. This coupling represents an active force which drives the system out of equilibrium and leads to the formation of nematic order interspersed with prominent, long-lived $+1$ defects. The defects in the nematic ordering are coupled to the shape of the ce…
▽ More
We study the vertex model for epithelial tissue mechanics extended to include coupling between the cell shapes and tensions in cell-cell junctions. This coupling represents an active force which drives the system out of equilibrium and leads to the formation of nematic order interspersed with prominent, long-lived $+1$ defects. The defects in the nematic ordering are coupled to the shape of the cell tiling, affecting cell areas and coordinations. This intricate interplay between cell shape, size, and coordination provides a possible mechanism by which tissues could spontaneously develop long-range polarity through local mechanical forces without resorting to long-range chemical patterning.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
Linear Viscoelastic Response of the Vertex Model with Internal and External Dissipation: Normal Modes Analysis
Authors:
Sijie Tong,
Rastko Sknepnek,
Andrej Kosmrlj
Abstract:
We use the normal mode formalism to study the shear rheology of the vertex model for epithelial tissue mechanics in the overdamped linear response regime. We consider systems with external (e.g., cell-substrate) and internal (e.g., cell-cell) dissipation mechanisms, and derive expressions for stresses on cells due to mechanical and dissipative forces. The semi-analytical method developed here is,…
▽ More
We use the normal mode formalism to study the shear rheology of the vertex model for epithelial tissue mechanics in the overdamped linear response regime. We consider systems with external (e.g., cell-substrate) and internal (e.g., cell-cell) dissipation mechanisms, and derive expressions for stresses on cells due to mechanical and dissipative forces. The semi-analytical method developed here is, however, general and can be directly applied to study the linear response of a broad class of soft matter systems with internal and external dissipation. It involves normal mode decomposition to calculate linear loss and storage moduli of the system. Specifically, displacements along each normal mode produce stresses due to elastic deformation and internal dissipation, which are in force balance with loads due to external dissipation. Each normal mode responds with a characteristic relaxation timescale, and its rheological behavior can be described as a combination of a standard linear solid element due to elastic stresses and a Jeffreys model element due to the internal dissipative stresses. The total response of the system is then fully determined by connecting in parallel all the viscoelastic elements corresponding to individual normal modes. This allows full characterization of the potentially complex linear rheological response of the system at all driving frequencies and identification of collective excitations. We show that internal and external dissipation mechanisms lead to qualitatively different rheological behaviors due to the presence or absence of Jeffreys elements, which is particularly pronounced at high driving frequencies. Our findings, therefore, underscore the importance of microscopic dissipation mechanisms in understanding the rheological behavior of soft materials and tissues, in particular.
△ Less
Submitted 2 October, 2022; v1 submitted 7 February, 2022;
originally announced February 2022.
-
Generating active T1 transitions through mechanochemical feedback
Authors:
Rastko Sknepnek,
Ilyas Djafer-Cherif,
Manli Chuai,
Cornelis J. Weijer,
Silke Henkes
Abstract:
Convergence-extension in embryos is controlled by chemical and mechanical signalling. A key cellular process is the exchange of neighbours via T1 transitions. We propose and analyse a model with positive feedback between recruitment of myosin motors and mechanical tension in cell junctions. The model produces active T1 events, which act to elongate the tissue perpendicular to external pulling. Usi…
▽ More
Convergence-extension in embryos is controlled by chemical and mechanical signalling. A key cellular process is the exchange of neighbours via T1 transitions. We propose and analyse a model with positive feedback between recruitment of myosin motors and mechanical tension in cell junctions. The model produces active T1 events, which act to elongate the tissue perpendicular to external pulling. Using an idealized tissue patch comprising several active cells embedded in a matrix of passive hexagonal cells we identified an optimal range of pulling forces to trigger an active T1 event. We show that pulling also generates tension chains in a realistic patch made entirely of active cells of random shapes, and leads to convergence-extension over a range of parameters. Our findings show that active intercalations can generate stress that activates T1 events in neighbouring cells resulting in tension dependent tissue re-organisation, in qualitative agreement with experiments on gastrulation in chick embryos.
△ Less
Submitted 3 April, 2022; v1 submitted 23 June, 2021;
originally announced June 2021.
-
Linear Viscoelastic Properties of the Vertex Model for Epithelial Tissues
Authors:
Sijie Tong,
Navreeta K. Singh,
Rastko Sknepnek,
Andrej Kosmrlj
Abstract:
Epithelial tissues act as barriers and, therefore, must repair themselves, respond to environmental changes and grow without compromising their integrity. Consequently, they exhibit complex viscoelastic rheological behavior where constituent cells actively tune their mechanical properties to change the overall response of the tissue, e.g., from solid-like to fluid-like. Mesoscopic mechanical prope…
▽ More
Epithelial tissues act as barriers and, therefore, must repair themselves, respond to environmental changes and grow without compromising their integrity. Consequently, they exhibit complex viscoelastic rheological behavior where constituent cells actively tune their mechanical properties to change the overall response of the tissue, e.g., from solid-like to fluid-like. Mesoscopic mechanical properties of epithelia are commonly modeled with the vertex model. While previous studies have predominantly focused on the rheological properties of the vertex model at long time scales, we systematically studied the full dynamic range by applying small oscillatory shear and bulk deformations in both solid-like and fluid-like phases for regular hexagonal and disordered cell configurations. We found that the shear and bulk responses in the fluid and solid phases can be described by standard spring-dashpot viscoelastic models. Furthermore, the solid-fluid transition can be tuned by applying pre-deformation to the system. Our study provides insights into the mechanisms by which epithelia can regulate their rich rheological behavior.
△ Less
Submitted 26 January, 2022; v1 submitted 22 February, 2021;
originally announced February 2021.
-
Wrinkle patterns in active viscoelastic thin sheets
Authors:
D. A. Matoz-Fernandez,
Fordyce A. Davidson,
Nicola R. Stanley-Wall,
Rastko Sknepnek
Abstract:
We show that a viscoelastic thin sheet driven out of equilibrium by active structural remodelling develops a rich variety of shapes as a result of a competition between viscous relaxation and activity. In the regime where active processes are faster than viscoelastic relaxation, wrinkles that are formed due to remodelling are unable to relax to a configuration that minimises the elastic energy and…
▽ More
We show that a viscoelastic thin sheet driven out of equilibrium by active structural remodelling develops a rich variety of shapes as a result of a competition between viscous relaxation and activity. In the regime where active processes are faster than viscoelastic relaxation, wrinkles that are formed due to remodelling are unable to relax to a configuration that minimises the elastic energy and the sheet is inherently out of equilibrium. We argue that this non-equilibrium regime is of particular interest in biology as it allows the system to access morphologies that are unavailable if restricted to the adiabatic evolution between configurations that minimise the elastic energy alone. Here, we introduce activity using the formalism of evolving target metric and showcase the diversity of wrinkling morphologies arising from out of equilibrium dynamics.
△ Less
Submitted 18 April, 2019;
originally announced April 2019.
-
Dense active matter model of motion patterns in confluent cell monolayers
Authors:
Silke Henkes,
Kaja Kostanjevec,
J. Martin Collinson,
Rastko Sknepnek,
Eric Bertin
Abstract:
Epithelial cell monolayers show remarkable displacement and velocity correlations over distances of ten or more cell sizes that are reminiscent of supercooled liquids and active nematics. We show that many observed features can be described within the framework of dense active matter, and argue that persistent uncoordinated cell motility coupled to the collective elastic modes of the cell sheet is…
▽ More
Epithelial cell monolayers show remarkable displacement and velocity correlations over distances of ten or more cell sizes that are reminiscent of supercooled liquids and active nematics. We show that many observed features can be described within the framework of dense active matter, and argue that persistent uncoordinated cell motility coupled to the collective elastic modes of the cell sheet is sufficient to produce swirl-like correlations. We obtain this result using both continuum active linear elasticity and a normal modes formalism, and validate analytical predictions with numerical simulations of two agent-based cell models, soft elastic particles and the self-propelled Voronoi model together with in-vitro experiments of confluent corneal epithelial cell sheets. Simulations and normal mode analysis perfectly match when tissue-level reorganisation occurs on times longer than the persistence time of cell motility. Our analytical model quantitatively matches measured velocity correlation functions over more than a decade with a single fitting parameter.
△ Less
Submitted 11 February, 2020; v1 submitted 15 January, 2019;
originally announced January 2019.
-
Dynamical patterns in active nematics on a sphere
Authors:
Silke Henkes,
M. Cristina Marchetti,
Rastko Sknepnek
Abstract:
Using agent-based simulations of self-propelled particles subject to short-range repulsion and nematic alignment we explore the dynamical phases of a dense active material confined to the surface of a sphere. We map the dynamical phase diagram as a function of curvature, alignment strength and activity and reproduce phases seen in recent experiments on active microtubules moving on the surfaces of…
▽ More
Using agent-based simulations of self-propelled particles subject to short-range repulsion and nematic alignment we explore the dynamical phases of a dense active material confined to the surface of a sphere. We map the dynamical phase diagram as a function of curvature, alignment strength and activity and reproduce phases seen in recent experiments on active microtubules moving on the surfaces of vesicles. At low driving, we recover the equilibrium nematic ground state with four +1/2 defects. As the driving is increased, geodesic forces drive the transition to a band of polar matter wrapping around an equator, with large bald spots corresponding to two +1 defects at the poles. Finally, bands fold onto themselves, followed by the system moving into a turbulent state marked by active proliferation of pairs of topological defects. We highlight the role of nematic persistence length and time for pattern formation in these confined systems with finite curvature.
△ Less
Submitted 15 May, 2017;
originally announced May 2017.
-
Active Vertex Model for Cell-Resolution Description of Epithelial Tissue Mechanics
Authors:
Daniel L. Barton,
Silke Henkes,
Cornelis J. Weijer,
Rastko Sknepnek
Abstract:
We introduce an Active Vertex Model (AVM) for cell-resolution studies of the mechanics of confluent epithelial tissues consisting of tens of thousands of cells, with a level of detail inaccessible to similar methods. The AVM combines the Vertex Model for confluent epithelial tissues with active matter dynamics. This introduces a natural description of the cell motion and accounts for motion patter…
▽ More
We introduce an Active Vertex Model (AVM) for cell-resolution studies of the mechanics of confluent epithelial tissues consisting of tens of thousands of cells, with a level of detail inaccessible to similar methods. The AVM combines the Vertex Model for confluent epithelial tissues with active matter dynamics. This introduces a natural description of the cell motion and accounts for motion patterns observed on multiple scales. Furthermore, cell contacts are generated dynamically from positions of cell centres. This not only enables efficient numerical implementation, but provides a natural description of the T1 transition events responsible for local tissue rearrangements. The AVM also includes cell alignment, cell-specific mechanical properties, cell growth, division and apoptosis. In addition, the AVM introduces a flexible, dynamically changing boundary of the epithelial sheet allowing for studies of phenomena such as the fingering instability or wound healing. We illustrate these capabilities with a number of case studies.
△ Less
Submitted 18 December, 2016;
originally announced December 2016.
-
Cell division and death inhibit glassy behaviour of confluent tissues
Authors:
D. A. Matoz-Fernandez,
K. Martens,
Rastko Sknepnek,
J. L. Barrat,
S. Henkes
Abstract:
We investigate the effects of cell division and apopotosis on collective dynamics in two-dimensional epithelial tissues. Our model includes three key ingredients observed across many epithelia, namely cell-cell adhesion, cell death and a cell division process that depends on the surrounding environment. We show a rich non-equilibrium phase diagram depending on the ratio of cell death to cell divis…
▽ More
We investigate the effects of cell division and apopotosis on collective dynamics in two-dimensional epithelial tissues. Our model includes three key ingredients observed across many epithelia, namely cell-cell adhesion, cell death and a cell division process that depends on the surrounding environment. We show a rich non-equilibrium phase diagram depending on the ratio of cell death to cell division and on the adhesion strength. For large apopotosis rates, cells die out and the tissue disintegrates. As the death rate decreases, however, we show, consecutively, the existence of a gas-like phase, a gel-like phase, and a dense confluent (tissue) phase. Most striking is the observation that the tissue is self-melting through its own internal activity, ruling out the existence of any glassy phase.
△ Less
Submitted 17 April, 2017; v1 submitted 28 October, 2016;
originally announced October 2016.
-
Non-Hookean statistical mechanics of clamped graphene ribbons
Authors:
Mark J. Bowick,
Andrej Kosmrlj,
David R. Nelson,
Rastko Sknepnek
Abstract:
Thermally fluctuating sheets and ribbons provide an intriguing forum in which to investigate strong violations of Hooke's Law: large distance elastic parameters are in fact not constant, but instead depend on the macroscopic dimensions. Inspired by recent experiments on free-standing graphene cantilevers, we combine the statistical mechanics of thin elastic plates and large-scale numerical simulat…
▽ More
Thermally fluctuating sheets and ribbons provide an intriguing forum in which to investigate strong violations of Hooke's Law: large distance elastic parameters are in fact not constant, but instead depend on the macroscopic dimensions. Inspired by recent experiments on free-standing graphene cantilevers, we combine the statistical mechanics of thin elastic plates and large-scale numerical simulations to investigate the thermal renormalization of the bending rigidity of graphene ribbons clamped at one end. For ribbons of dimensions $W\times L$ (with $L\geq W$), the macroscopic bending rigidity $κ_R$ determined from cantilever deformations is independent of the width when $W<\ell_\textrm{th}$, where $\ell_\textrm{th}$ is a thermal length scale, as expected. When $W>\ell_\textrm{th}$, however, this thermally renormalized bending rigidity begins to systematically increase, in agreement with the scaling theory, although in our simulations we were not quite able to reach the system sizes necessary to determine the fully developed power law dependence on $W$. When the ribbon length $L > \ell_p$, where $\ell_p$ is the $W$-dependent thermally renormalized ribbon persistence length, we observe a scaling collapse and the beginnings of large scale random walk behavior.
△ Less
Submitted 19 February, 2017; v1 submitted 15 August, 2016;
originally announced August 2016.
-
Dynamically Generated Patterns in Dense Suspensions of Active Filaments
Authors:
K. R. Prathyusha,
Silke Henkes,
Rastko Sknepnek
Abstract:
We use Langevin dynamics simulations to study dynamical behaviour of a dense planar layer of active semi-flexible filaments. Using the strength of active force and the thermal persistence length as parameters, we map a detailed phase diagram and identify several non-equilibrium phases in this system. In addition to a slowly flowing melt phase, we observe that for sufficiently high activity, collec…
▽ More
We use Langevin dynamics simulations to study dynamical behaviour of a dense planar layer of active semi-flexible filaments. Using the strength of active force and the thermal persistence length as parameters, we map a detailed phase diagram and identify several non-equilibrium phases in this system. In addition to a slowly flowing melt phase, we observe that for sufficiently high activity, collective flow accompanied by signatures of local polar and nematic order appears in the system. This state is also characterised by strong density fluctuations. Furthermore, we identify an activity-driven cross-over from this state of coherently flowing bundles of filaments to a phase with no global flow, formed by individual filaments coiled into rotating spirals. This suggests a mechanism where the system responds to activity by changing the shape of active agents, an effect with no analogue in systems of active particles without internal degrees of freedom.
△ Less
Submitted 20 November, 2017; v1 submitted 10 August, 2016;
originally announced August 2016.
-
Stiffening Thermal Membranes by Cutting
Authors:
Emily Russell,
Rastko Sknepnek,
Mark J. Bowick
Abstract:
Two-dimensional crystalline membranes have recently been realized experimentally in such systems as graphene and molybdenum disulfide, sparking a resurgence in interest in their statistical properties. Thermal fluctuations can significantly affect the effective mechanical properties of properly thermalized membranes, renormalizing both bending rigidity and elastic moduli so that in particular they…
▽ More
Two-dimensional crystalline membranes have recently been realized experimentally in such systems as graphene and molybdenum disulfide, sparking a resurgence in interest in their statistical properties. Thermal fluctuations can significantly affect the effective mechanical properties of properly thermalized membranes, renormalizing both bending rigidity and elastic moduli so that in particular they become stiffer to bending than their bare bending rigidity would suggest. We use molecular dynamics simulations to examine how the mechanical behavior of thermalized two-dimensional clamped ribbons (cantilevers) depends on their precise topology and geometry. We find that a simple slit smooths roughness as measured by the variance of height fluctuations. This counterintuitive effect may be due to the counter-posed coupling of the lips of the slit to twist in the intact regions of the ribbon.
△ Less
Submitted 2 June, 2017; v1 submitted 15 December, 2015;
originally announced December 2015.
-
Effects of scars on crystalline shell stability under external pressure
Authors:
Duanduan Wan,
Mark J. Bowick,
Rastko Sknepnek
Abstract:
We study how the stability of spherical crystalline shells under external pressure is influenced by the defect structure. In particular, we compare stability for shells with a minimal set of topologically-required defects to shells with extended defect arrays (grain boundary "scars" with non-vanishing net disclination charge). We perform Monte Carlo simulations to compare how shells with and witho…
▽ More
We study how the stability of spherical crystalline shells under external pressure is influenced by the defect structure. In particular, we compare stability for shells with a minimal set of topologically-required defects to shells with extended defect arrays (grain boundary "scars" with non-vanishing net disclination charge). We perform Monte Carlo simulations to compare how shells with and without scars deform quasi-statically under external hydrostatic pressure. We find that the critical pressure at which shells collapse is lowered for scarred configurations that break icosahedral symmetry and raised for scars that preserve icosahedral symmetry. The particular shapes which arise from breaking of an initial icosahedrally-symmetric shell depend on the Föppl-von Kármán number.
△ Less
Submitted 22 October, 2014;
originally announced October 2014.
-
Active swarms on a sphere
Authors:
Rastko Sknepnek,
Silke Henkes
Abstract:
Here we show that coupling to curvature has profound effects on collective motion in active systems, leading to patterns not observed in flat space. Biological examples of such active motion in curved environments are numerous: curvature and tissue folding are crucial during gastrulation, epithelial and endothelial cells move on constantly growing, curved crypts and vili in the gut, and the mammal…
▽ More
Here we show that coupling to curvature has profound effects on collective motion in active systems, leading to patterns not observed in flat space. Biological examples of such active motion in curved environments are numerous: curvature and tissue folding are crucial during gastrulation, epithelial and endothelial cells move on constantly growing, curved crypts and vili in the gut, and the mammalian corneal epithelium grows in a steady-state vortex pattern. On the physics side, droplets coated with actively driven microtubule bundles show active nematic patterns. We study a model of self-propelled particles with polar alignment on a sphere. Hallmarks of these motion patterns are a polar vortex and a circulating band arising due to the incompatibility between spherical topology and uniform motion - a consequence of the hairy ball theorem. We present analytical results showing that frustration due to curvature leads to stable elastic distortions storing energy in the band.
△ Less
Submitted 31 July, 2014;
originally announced July 2014.
-
Defect dynamics in active nematics
Authors:
Luca Giomi,
Mark J. Bowick,
Prashant Mishra,
Rastko Sknepnek,
M. Cristina Marchetti
Abstract:
Topological defects are distinctive signatures of liquid crystals. They profoundly affect the viscoelastic behavior of the fluid by constraining the orientational structure in a way that inevitably requires global changes not achievable with any set of local deformations. In active nematic liquid crystals topological defects not only dictate the global structure of the director, but also act as lo…
▽ More
Topological defects are distinctive signatures of liquid crystals. They profoundly affect the viscoelastic behavior of the fluid by constraining the orientational structure in a way that inevitably requires global changes not achievable with any set of local deformations. In active nematic liquid crystals topological defects not only dictate the global structure of the director, but also act as local sources of motion, behaving as self-propelled particles. In this article we present a detailed analytical and numerical study of the mechanics of topological defects in active nematic liquid crystals.
△ Less
Submitted 26 March, 2014; v1 submitted 20 March, 2014;
originally announced March 2014.
-
Optimal shapes and stresses of adherent cells on patterned substrates
Authors:
Shiladitya Banerjee,
Rastko Sknepnek,
M. Cristina Marchetti
Abstract:
We investigate a continuum mechanical model for an adherent cell on two dimensional adhesive micropatterned substrates. The cell is modeled as an isotropic and homogeneous elastic material subject to uniform internal contractile stresses. The build-up of tension from cortical actin bundles at the cell periphery is incorporated by introducing an energy cost for bending of the cell boundary, resulti…
▽ More
We investigate a continuum mechanical model for an adherent cell on two dimensional adhesive micropatterned substrates. The cell is modeled as an isotropic and homogeneous elastic material subject to uniform internal contractile stresses. The build-up of tension from cortical actin bundles at the cell periphery is incorporated by introducing an energy cost for bending of the cell boundary, resulting to a resistance to changes in local curvature. Integrin-based adhesions are modeled as harmonic springs, that pin the cell to adhesive patches of a predefined geometry. Using Monte Carlo simulations and analytical techniques we investigate the competing effects of bulk contractility and cortical bending rigidity in regulating cell shapes on non-adherent regions. We show that the crossover from convex to concave cell edges is controlled by the interplay between contractile stresses and boundary bending rigidity. In particular, the cell boundary becomes concave beyond a critical value of the contractile stress that is proportional to the cortical bending rigidity. Furthermore, the intracellular stresses are found largely concentrated at the concave edge of the cell. The model can be used to generate a cell-shape phase diagram for each specific adhesion geometry.
△ Less
Submitted 24 December, 2013;
originally announced December 2013.
-
On the modeling of endocytosis in yeast
Authors:
T. Zhang,
R. Sknepnek,
M. J. Bowick,
J. M. Schwarz
Abstract:
The cell membrane deforms during endocytosis to surround extracellular material and draw it into the cell. Experiments on endocytosis in yeast all agree that (i) actin polymerizes into a network of filaments exerting active forces on the membrane to deform it and (ii) the large scale membrane deformation is tubular in shape. There are three competing proposals, in contrast, for precisely how the a…
▽ More
The cell membrane deforms during endocytosis to surround extracellular material and draw it into the cell. Experiments on endocytosis in yeast all agree that (i) actin polymerizes into a network of filaments exerting active forces on the membrane to deform it and (ii) the large scale membrane deformation is tubular in shape. There are three competing proposals, in contrast, for precisely how the actin filament network organizes itself to drive the deformation. We use variational approaches and numerical simulations to address this competition by analyzing a meso-scale model of actin-mediated endocytosis in yeast. The meso-scale model breaks up the invagination process into three stages: (i) initiation, where clathrin interacts with the membrane via adaptor proteins, (ii) elongation, where the membrane is then further deformed by polymerizing actin filaments, followed by (iii) pinch-off. Our results suggest that the pinch-off mechanism may be assisted by a pearling-like instability. We rule out two of the three competing proposals for the organization of the actin filament network during the elongation stage. These two proposals could possibly be important in the pinch-off stage, however, where additional actin polymerization helps break off the vesicle. Implications and comparisons with earlier modeling of endocytosis in yeast are discussed.
△ Less
Submitted 4 December, 2014; v1 submitted 31 October, 2013;
originally announced October 2013.
-
Shapes of pored membranes
Authors:
Zhenwei Yao,
Rastko Sknepnek,
Creighton K. Thomas,
Monica Olvera de la Cruz
Abstract:
We study the shapes of pored membranes within the framework of the Helfrich theory under the constraints of fixed area and pore size. We show that the mean curvature term leads to a budding- like structure, while the Gaussian curvature term tends to flatten the membrane near the pore; this is corroborated by simulation. We propose a scheme to deduce the ratio of the Gaussian rigidity to the bendin…
▽ More
We study the shapes of pored membranes within the framework of the Helfrich theory under the constraints of fixed area and pore size. We show that the mean curvature term leads to a budding- like structure, while the Gaussian curvature term tends to flatten the membrane near the pore; this is corroborated by simulation. We propose a scheme to deduce the ratio of the Gaussian rigidity to the bending rigidity simply by observing the shape of the pored membrane. This ratio is usually difficult to measure experimentally. In addition, we briefly discuss the stability of a pore by relaxing the constraint of a fixed pore size and adding the line tension. Finally, the flattening effect due to the Gaussian curvature as found in studying pored membranes is extended to two-component membranes. We find that sufficiently high contrast between the components' Gaussian rigidities leads to budding which is distinct from that due to the line tension.
△ Less
Submitted 7 May, 2013;
originally announced May 2013.
-
Pathways to faceting of vesicles
Authors:
Mark Bowick,
Rastko Sknepnek
Abstract:
The interplay between geometry, topology and order can lead to geometric frustration that profoundly affects the shape and structure of a curved surface. In this commentary we show how frustration in this context can result in the faceting of elastic vesicles. We show that, under the right conditions, an assortment of regular and irregular polyhedral structures may be the low energy states of elas…
▽ More
The interplay between geometry, topology and order can lead to geometric frustration that profoundly affects the shape and structure of a curved surface. In this commentary we show how frustration in this context can result in the faceting of elastic vesicles. We show that, under the right conditions, an assortment of regular and irregular polyhedral structures may be the low energy states of elastic membranes with spherical topology. In particular, we show how topological defects, necessarily present in any crystalline lattice confined to spherical topology, naturally lead to the formation of icosahedra in a homogeneous elastic vesicle. Furthermore, we show that introducing heterogeneities in the elastic properties, or allowing for non-linear bending response of a homogeneous system, opens non-trivial pathways to the formation of faceted, yet non-icosahedral, structures.
△ Less
Submitted 9 July, 2013; v1 submitted 7 March, 2013;
originally announced March 2013.
-
Nonlinear elastic model for faceting of vesicles with soft grain boundaries
Authors:
Rastko Sknepnek,
Monica Olvera de la Cruz
Abstract:
We use an elastic model to explore faceting of solid-wall vesicles with elastic heterogeneities. We show that faceting occurs in regions where the vesicle wall is softer, such as areas of reduced wall thicknesses or concentrated in crystalline defects. The elastic heterogeneities are modeled as a second component with reduced elastic parameters. Using simulated annealing Monte Carlo simulations we…
▽ More
We use an elastic model to explore faceting of solid-wall vesicles with elastic heterogeneities. We show that faceting occurs in regions where the vesicle wall is softer, such as areas of reduced wall thicknesses or concentrated in crystalline defects. The elastic heterogeneities are modeled as a second component with reduced elastic parameters. Using simulated annealing Monte Carlo simulations we obtain the vesicle shape by optimizing the distributions of facets and boundaries. Our model allows us to reduce the effects of the residual stress generated by crystalline defects, and reveals a robust faceting mechanism into polyhedra other than the icosahedron.
△ Less
Submitted 29 May, 2012;
originally announced May 2012.
-
A curvature-driven effective attraction in multicomponent membranes
Authors:
Matthew F. Demers,
Rastko Sknepnek,
Monica Olvera de la Cruz
Abstract:
We study closed liquid membranes that segregate into three phases due to differences in the chemical and physical properties of its components. The shape and in-plane membrane arrangement of the phases are coupled through phase-specific bending energies and line tensions. We use simulated annealing Monte Carlo simulations to find low-energy structures, allowing both phase arrangement and membrane…
▽ More
We study closed liquid membranes that segregate into three phases due to differences in the chemical and physical properties of its components. The shape and in-plane membrane arrangement of the phases are coupled through phase-specific bending energies and line tensions. We use simulated annealing Monte Carlo simulations to find low-energy structures, allowing both phase arrangement and membrane shape to relax. The three-phase system is the simplest one in which there are multiple interface pairs, allowing us to analyze interfacial preferences and pairwise distinct line tensions. We observe the system's preference for interface pairs that maximize differences in spontaneous curvature. From a pattern selection perspective, this acts as an effective attraction between phases of most disparate spontaneous curvature. We show that this effective attraction is robust enough to persist even when the interface between these phases is the most penalized by line tension. This effect is driven by geometry and not by any explicit component-component interaction.
△ Less
Submitted 31 July, 2012; v1 submitted 28 May, 2012;
originally announced May 2012.
-
Planar sheets meet negative curvature liquid interfaces
Authors:
Zhenwei Yao,
Mark Bowick,
Xu Ma,
Rastko Sknepnek
Abstract:
If an inextensible thin sheet is adhered to a substrate with a negative Gaussian curvature it will experience stress due to geometric frustration. We analyze the consequences of such geometric frustration using analytic arguments and numerical simulations. Both concentric wrinkles and eye-like folds are shown to be compatible with negative curvatures. Which pattern will be realized depends on the…
▽ More
If an inextensible thin sheet is adhered to a substrate with a negative Gaussian curvature it will experience stress due to geometric frustration. We analyze the consequences of such geometric frustration using analytic arguments and numerical simulations. Both concentric wrinkles and eye-like folds are shown to be compatible with negative curvatures. Which pattern will be realized depends on the curvature of the substrate. We discuss both types of folding patterns and determine the phase diagram governing their appearance.
△ Less
Submitted 20 November, 2012; v1 submitted 23 April, 2012;
originally announced April 2012.
-
Spectral analysis for the iron-based superconductors: Anisotropic spin fluctuations and fully gapped s^{\pm}-wave superconductivity
Authors:
Junhua Zhang,
Rastko Sknepnek,
Jörg Schmalian
Abstract:
Spin fluctuations are considered to be one of the candidates that drive a sign-reversed s^{\pm} superconducting state in the iron pnictides. In the magnetic scenario, whether the spin fluctuation spectrum exhibits certain unique fine structures is an interesting aspect for theoretical study in order to understand experimental observations. We investigate the detailed momentum dependence of the sho…
▽ More
Spin fluctuations are considered to be one of the candidates that drive a sign-reversed s^{\pm} superconducting state in the iron pnictides. In the magnetic scenario, whether the spin fluctuation spectrum exhibits certain unique fine structures is an interesting aspect for theoretical study in order to understand experimental observations. We investigate the detailed momentum dependence of the short-range spin fluctuations using a 2-orbital model in the self-consistent fluctuation exchange approximation and find that a common feature of those fluctuations that are capable of inducing a fully gapped s^{\pm} state is the momentum anisotropy with lengthened span along the direction transverse to the antiferromagnetic momentum transfer. Performing a qualitative analysis based on the orbital character and the deviation from perfect nesting of the electronic structure for the 2-orbital and a more complete 5-orbital model, we gain the insight that this type of anisotropic spin fluctuations favor superconductivity due to their enhancement of intra-orbital, but inter-band, pair scattering processes. The momentum anisotropy leads to elliptically shaped magnetic responses which have been observed in inelastic neutron scattering measurements. Meanwhile, our detailed study on the magnetic and the electronic spectrum shows that the dispersion of the magnetic resonance mode in the nearly isotropic s^{\pm} superconducting state exhibits anisotropic propagating behavior in an upward pattern and the coupling of the resonance mode to fermions leads to a dip feature in the spectral function.
△ Less
Submitted 21 October, 2010; v1 submitted 29 July, 2010;
originally announced July 2010.
-
Orbital coupling and superconductivity in the iron pnictides
Authors:
Junhua Zhang,
Rastko Sknepnek,
Rafael M. Fernandes,
Joerg Schmalian
Abstract:
We demonstrate that strong inter-orbital interaction is very efficient to achieve superconductivity due to magnetic fluctuations in the iron pnictides. Fermi surface states that are coupled by the antiferromagnetic wave vector are often of different orbital nature, causing pair-hopping interactions between distinct Fe-3d orbitals to become important. Performing a self-consistent FLEX calculation…
▽ More
We demonstrate that strong inter-orbital interaction is very efficient to achieve superconductivity due to magnetic fluctuations in the iron pnictides. Fermi surface states that are coupled by the antiferromagnetic wave vector are often of different orbital nature, causing pair-hopping interactions between distinct Fe-3d orbitals to become important. Performing a self-consistent FLEX calculation below $T_{c}$ we determine the superconducting order parameter as function of intra- and inter-orbital couplings. We find an $s^{\pm}$-pairing state with $T_{c}\simeq80\mathrm{K}$ for realistic parameters.
△ Less
Submitted 20 May, 2009; v1 submitted 25 March, 2009;
originally announced March 2009.
-
Anisotropic pairing in the iron pnictides
Authors:
Rastko Sknepnek,
German Samolyuk,
Yong-bin Lee,
Joerg Schmalian
Abstract:
We determine the anisotropy of the spin fluctuation induced pairing gap on the Fermi surface of the FeAs based superconductors as function of the exchange and Hund's coupling $J_{H}$. We find that for sufficiently large $J_{H}$, nearly commensurate magnetic fluctuations yield a fully gapped $s^{\pm}$-pairing state with small anisotropy of the gap amplitude on each Fermi surface sheet, but signif…
▽ More
We determine the anisotropy of the spin fluctuation induced pairing gap on the Fermi surface of the FeAs based superconductors as function of the exchange and Hund's coupling $J_{H}$. We find that for sufficiently large $J_{H}$, nearly commensurate magnetic fluctuations yield a fully gapped $s^{\pm}$-pairing state with small anisotropy of the gap amplitude on each Fermi surface sheet, but significant variations of the gap amplitude for different sheets of the Fermi surface. In particular, we obtain the large variation of the gap amplitude on different Fermi surface sheets, as seen in ARPES experiments. For smaller values of Hund's coupling incommensurate magnetic fluctuations yield an $s^{\pm}$-pairing state with line nodes. Such a state is also possible once the anisotropy of the material is reduced and three dimensional effects come into play.
△ Less
Submitted 25 January, 2009; v1 submitted 28 July, 2008;
originally announced July 2008.
-
Universality of liquid-gas Mott transitions at finite temperatures
Authors:
Stefanos Papanikolaou,
Rafael Monteiro Fernandes,
Eduardo Fradkin,
Philip W. Phillips,
Joerg Schmalian,
Rastko Sknepnek
Abstract:
We explain in a consistent manner the set of seemingly conflicting experiments on the finite temperature Mott critical point, and demonstrate that the Mott transition is in the Ising universality class. We show that, even though the thermodynamic behavior of the system near such critical point is described by an Ising order parameter, the global conductivity can depend on other singular observab…
▽ More
We explain in a consistent manner the set of seemingly conflicting experiments on the finite temperature Mott critical point, and demonstrate that the Mott transition is in the Ising universality class. We show that, even though the thermodynamic behavior of the system near such critical point is described by an Ising order parameter, the global conductivity can depend on other singular observables and, in particular, on the energy density. Finally, we show that in the presence of weak disorder the dimensionality of the system has crucial effects on the size of the critical region that is probed experimentally.
△ Less
Submitted 8 October, 2007;
originally announced October 2007.
-
Quantum phase transitions of the diluted O(3) rotor model
Authors:
Thomas Vojta,
Rastko Sknepnek
Abstract:
We study the phase diagram and the quantum phase transitions of a site-diluted two-dimensional O(3) quantum rotor model by means of large-scale Monte-Carlo simulations. This system has two quantum phase transitions, a generic one for small dilutions, and a percolation transition across the lattice percolation threshold. We determine the critical behavior for both transitions and for the multicri…
▽ More
We study the phase diagram and the quantum phase transitions of a site-diluted two-dimensional O(3) quantum rotor model by means of large-scale Monte-Carlo simulations. This system has two quantum phase transitions, a generic one for small dilutions, and a percolation transition across the lattice percolation threshold. We determine the critical behavior for both transitions and for the multicritical point that separates them. In contrast to the exotic scaling scenarios found in other random quantum systems, all these transitions are characterized by finite-disorder fixed points with power-law scaling. We relate our findings to a recent classification of phase transitions with quenched disorder according to the rare region dimensionality, and we discuss experiments in disordered quantum magnets.
△ Less
Submitted 13 September, 2006; v1 submitted 6 June, 2006;
originally announced June 2006.
-
Spin Structure Factor of the Frustrated Quantum Magnet Cs_2CuCl_4
Authors:
Denis Dalidovich,
Rastko Sknepnek,
A. John Berlinsky,
Junhua Zhang,
Catherine Kallin
Abstract:
The ground state properties and neutron structure factor for the two-dimensional antiferromagnet on the triangular lattice, with uni-directional anisotropy in the nearest-neighbor exchange couplings and a weak Dzyaloshinskii-Moriya (DM) interaction, are studied. This Hamiltonian has been used to interpret neutron scattering measurements on the spin 1/2 spiral spin-density-wave system, Cs_2CuCl_4…
▽ More
The ground state properties and neutron structure factor for the two-dimensional antiferromagnet on the triangular lattice, with uni-directional anisotropy in the nearest-neighbor exchange couplings and a weak Dzyaloshinskii-Moriya (DM) interaction, are studied. This Hamiltonian has been used to interpret neutron scattering measurements on the spin 1/2 spiral spin-density-wave system, Cs_2CuCl_4, [R. Coldea, et al., Phys. Rev. B 68, 134424 (2003)]. Calculations are performed using a 1/S expansion, taking into account interactions between spin-waves. The ground state energy, the shift of the ordering wave-vector, Q, and the local magnetization are all calculated to order 1/S^2. The neutron structure factor, obtained using anharmonic spin-wave Green's functions to order 1/S, is shown to be in reasonable agreement with published neutron data, provided that slightly different parameters are used for the exchange and DM interactions than those inferred from measurements in high magnetic field.
△ Less
Submitted 20 January, 2006;
originally announced January 2006.
-
Dynamics at a smeared phase transition
Authors:
Bernard Fendler,
Rastko Sknepnek,
Thomas Vojta
Abstract:
We investigate the effects of rare regions on the dynamics of Ising magnets with planar defects, i.e., disorder perfectly correlated in two dimensions. In these systems, the magnetic phase transition is smeared because static long-range order can develop on isolated rare regions. We first study an infinite-range model by numerically solving local dynamic mean-field equations. Then we use extrema…
▽ More
We investigate the effects of rare regions on the dynamics of Ising magnets with planar defects, i.e., disorder perfectly correlated in two dimensions. In these systems, the magnetic phase transition is smeared because static long-range order can develop on isolated rare regions. We first study an infinite-range model by numerically solving local dynamic mean-field equations. Then we use extremal statistics and scaling arguments to discuss the dynamics beyond mean-field theory. In the tail region of the smeared transition the dynamics is even slower than in a conventional Griffiths phase: the spin autocorrelation function decays like a stretched exponential at intermediate times before approaching the exponentially small equilibrium value following a power law at late times.
△ Less
Submitted 2 March, 2005; v1 submitted 27 September, 2004;
originally announced September 2004.
-
Critical points and quenched disorder: From Harris criterion to rare regions and smearing
Authors:
Thomas Vojta,
Rastko Sknepnek
Abstract:
We consider the influence of quenched spatial disorder on phase transitions in classical and quantum systems. We show that rare strong disorder fluctuations can have dramatic effects on critical points. In classical systems with sufficiently correlated disorder or in quantum systems with overdamped dynamics they can completely destroy the sharp phase transition by smearing. This is caused by eff…
▽ More
We consider the influence of quenched spatial disorder on phase transitions in classical and quantum systems. We show that rare strong disorder fluctuations can have dramatic effects on critical points. In classical systems with sufficiently correlated disorder or in quantum systems with overdamped dynamics they can completely destroy the sharp phase transition by smearing. This is caused by effects similar to but stronger than Griffiths phenomena: True static order can develop on a rare region while the bulk system is still in the disordered phase. We discuss the thermodynamic behavior in the vicinity of such a smeared transition using optimal fluctuation theory, and we present numerical results for a two-dimensional model system.
△ Less
Submitted 1 September, 2004; v1 submitted 4 May, 2004;
originally announced May 2004.
-
Exotic vs. conventional scaling and universality in a disordered bilayer quantum Heisenberg antiferromagnet
Authors:
Rastko Sknepnek,
Thomas Vojta,
Matthias Vojta
Abstract:
We present large-scale Monte-Carlo simulations of a two-dimensional (2d) bilayer quantum Heisenberg antiferromagnet with random dimer dilution. In contrast to the exotic scaling scenarios found in many other random quantum systems, the quantum phase transition in this system is characterized by a finite-disorder fixed point with power-law scaling. After accounting for strong corrections to scali…
▽ More
We present large-scale Monte-Carlo simulations of a two-dimensional (2d) bilayer quantum Heisenberg antiferromagnet with random dimer dilution. In contrast to the exotic scaling scenarios found in many other random quantum systems, the quantum phase transition in this system is characterized by a finite-disorder fixed point with power-law scaling. After accounting for strong corrections to scaling, characterized by a leading irrelevant exponent of ω= 0.48, we find universal, i.e., disorder-independent, critical exponents z=1.310(6) and ν=1.16(3). We discuss the consequences of these findings and suggest new experiments.
△ Less
Submitted 1 September, 2004; v1 submitted 12 February, 2004;
originally announced February 2004.
-
Smeared phase transition in a three-dimensional Ising model with planar defects: Monte-Carlo simulations
Authors:
Rastko Sknepnek,
Thomas Vojta
Abstract:
We present results of large-scale Monte Carlo simulations for a three-dimensional Ising model with short range interactions and planar defects, i.e., disorder perfectly correlated in two dimensions. We show that the phase transition in this system is smeared, i.e., there is no single critical temperature, but different parts of the system order at different temperatures. This is caused by effect…
▽ More
We present results of large-scale Monte Carlo simulations for a three-dimensional Ising model with short range interactions and planar defects, i.e., disorder perfectly correlated in two dimensions. We show that the phase transition in this system is smeared, i.e., there is no single critical temperature, but different parts of the system order at different temperatures. This is caused by effects similar to but stronger than Griffiths phenomena. In an infinite-size sample there is an exponentially small but finite probability to find an arbitrary large region devoid of impurities. Such a rare region can develop true long-range order while the bulk system is still in the disordered phase. We compute the thermodynamic magnetization and its finite-size effects, the local magnetization, and the probability distribution of the ordering temperatures for different samples. Our Monte-Carlo results are in good agreement with a recent theory based on extremal statistics.
△ Less
Submitted 13 May, 2004; v1 submitted 17 November, 2003;
originally announced November 2003.
-
Order parameter symmetry and mode coupling effects at dirty superconducting quantum phase transitions
Authors:
Rastko Sknepnek,
Thomas Vojta,
Rajesh Narayanan
Abstract:
We derive an order-parameter field theory for a quantum phase transition between a disordered metal and an exotic (non-s-wave) superconductor. Mode coupling effects between the order parameter and other fermionic soft modes lead to an effective long-range interaction between the anomalous density fluctuations which is reflected in singularities in the free energy functional. However, this long-r…
▽ More
We derive an order-parameter field theory for a quantum phase transition between a disordered metal and an exotic (non-s-wave) superconductor. Mode coupling effects between the order parameter and other fermionic soft modes lead to an effective long-range interaction between the anomalous density fluctuations which is reflected in singularities in the free energy functional. However, this long-range interaction is not strong enough to suppress disorder fluctuations. The asymptotic critical region is characterized by run-away flow to large disorder. For weak coupling, this asymptotic region is very narrow. It is preempted by a wide crossover regime with mean-field critical behavior and, in the p-wave case, logarithmic corrections to scaling in all dimensions.
△ Less
Submitted 27 September, 2004; v1 submitted 22 November, 2002;
originally announced November 2002.
-
Infinite-randomness critical point in the itinerant quantum antiferromagnet
Authors:
Rastko Sknepnek,
Thomas Vojta
Abstract:
We study the quantum phase transition in the three-dimensional disordered itinerant antiferromagnet by Monte-Carlo simulations of the order-parameter field theory. We find strong evidence for the transition being controlled by an infinite-randomness fixed point: The dynamical scaling is activated, i.e., the logarithm of the energy scales like a power of the length, implying a dynamical exponent…
▽ More
We study the quantum phase transition in the three-dimensional disordered itinerant antiferromagnet by Monte-Carlo simulations of the order-parameter field theory. We find strong evidence for the transition being controlled by an infinite-randomness fixed point: The dynamical scaling is activated, i.e., the logarithm of the energy scales like a power of the length, implying a dynamical exponent of infinity. The probability distribution of the energy gaps is very broad and becomes broader with increasing system size, even on a logarithmic scale.
△ Less
Submitted 11 June, 2002;
originally announced June 2002.
-
The quantum phase transition of itinerant helimagnets
Authors:
Thomas Vojta,
Rastko Sknepnek
Abstract:
We investigate the quantum phase transition of itinerant electrons from a paramagnet to a state which displays long-period helical structures due to a Dzyaloshinskii instability of the ferromagnetic state. In particular, we study how the self-generated effective long-range interaction recently identified in itinerant quantum ferromagnets is cut-off by the helical ordering. We find that for a suf…
▽ More
We investigate the quantum phase transition of itinerant electrons from a paramagnet to a state which displays long-period helical structures due to a Dzyaloshinskii instability of the ferromagnetic state. In particular, we study how the self-generated effective long-range interaction recently identified in itinerant quantum ferromagnets is cut-off by the helical ordering. We find that for a sufficiently strong Dzyaloshinskii instability the helimagnetic quantum phase transition is of second order with mean-field exponents. In contrast, for a weak Dzyaloshinskii instability the transition is analogous to that in itinerant quantum ferromagnets, i.e. it is of first order, as has been observed in MnSi.
△ Less
Submitted 1 October, 2000; v1 submitted 28 September, 2000;
originally announced September 2000.