BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

J. M. Schwarz

J. M. Schwarz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softBiological Physicscond-mat.dis-nncond-mat.stat-mechCell BehaviorMachine LearningNeural and Evolutionary ComputingSubcellular Processescond-mat.mtrl-scimath-phmath.MPNeurons and Cognitionphysics.chem-phQuantitative MethodsTissues and Organs

Catalog footprint

What is connected

30works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2023arXiv

A topologically-protected interior for three-dimensional confluent cellular collectives

Organoids are in vitro cellular collectives from which brain-like, or gut-like, or kidney-like structures emerge. To make quantitative predictions regarding the morphology and rheology of a cellular collective in its initial stages of development, we construct and study a three-dimensional vertex model. In such a model, the cells are represented as deformable polyhedrons with cells sharing faces such that there are no gaps between them, otherwise known as confluent. In a bulk model with periodic boundary conditions, we find a rigidity transition as a function of the target cell shape index $s_0$ with a critical value $s_0^*=5.39\pm0.01$. For a confluent cellular collective with a finite boundary, and in the presence of lateral extensile and in-plane, radial extensile deformations, we find a significant boundary-bulk effect that is one-cell layer thick. More specifically, for lateral extensile deformations, the cells in the bulk are much less aligned with the direction of the lateral deformation than the cells at the boundary. For in-plane, radial deformations, the cells in the bulk exhibit much less reorientation perpendicular to the radial direction than the cells at the boundary. In other words, for both deformations, the bulk, interior cells are topologically-protected from the deformations, at least over time scales much slower than the timescale for cellular rearrangements and up to reasonable amounts of strain. Our results provide an underlying mechanism for some observed cell shape patterning in organoids. Finally, we discuss the use of a cellular-based approach to designing organoids with new types of morphologies to study the intricate relationship between structure and function at the multi-cellular scale for example.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Effective medium theory of random regular networks

Disordered spring networks can exhibit rigidity transitions, due to either the removal of materials in over-constrained networks or the application of strain in under-constrained ones. While an effective medium theory (EMT) exists for the former, there is none for the latter. We, therefore, formulate an EMT for random regular networks, under-constrained spring networks with purely geometrical disorder, to predict their stiffness via the distribution of tensions. We find a linear dependence of stiffness on strain in the rigid phase and a nontrivial dependence on both the mean and standard deviation of the tension distribution. While EMT does not yield highly accurate predictions of shear modulus due to spatial heterogeneities, the noninvasiveness of this EMT makes it an ideal starting point for experimentalists quantifying the mechanics of such networks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Frequency propagation: Multi-mechanism learning in nonlinear physical networks

We introduce frequency propagation, a learning algorithm for nonlinear physical networks. In a resistive electrical circuit with variable resistors, an activation current is applied at a set of input nodes at one frequency, and an error current is applied at a set of output nodes at another frequency. The voltage response of the circuit to these boundary currents is the superposition of an `activation signal' and an `error signal' whose coefficients can be read in different frequencies of the frequency domain. Each conductance is updated proportionally to the product of the two coefficients. The learning rule is local and proved to perform gradient descent on a loss function. We argue that frequency propagation is an instance of a multi-mechanism learning strategy for physical networks, be it resistive, elastic, or flow networks. Multi-mechanism learning strategies incorporate at least two physical quantities, potentially governed by independent physical mechanisms, to act as activation and error signals in the training process. Locally available information about these two signals is then used to update the trainable parameters to perform gradient descent. We demonstrate how earlier work implementing learning via chemical signaling in flow networks also falls under the rubric of multi-mechanism learning.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning by non-interfering feedback chemical signaling in physical networks

Both non-neural and neural biological systems can learn. So rather than focusing on purely brain-like learning, efforts are underway to study learning in physical systems. Such efforts include equilibrium propagation (EP) and coupled learning (CL), which require storage of two different states-the free state and the perturbed state-during the learning process to retain information about gradients. Inspired by slime mold, we propose a new learning algorithm rooted in chemical signaling that does not require storage of two different states. Rather, the output error information is encoded in a chemical signal that diffuses into the network in a similar way as the activation/feedforward signal. The steady state feedback chemical concentration, along with the activation signal, stores the required gradient information locally. We apply our algorithm using a physical, linear flow network and test it using the Iris data set with 93% accuracy. We also prove that our algorithm performs gradient descent. Finally, in addition to comparing our algorithm directly with EP and CL, we address the biological plausibility of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermal Fluctuations of Singular Bar-Joint Mechanisms

A bar-joint mechanism is a deformable assembly of freely rotating joints connected by stiff bars. Here we develop a formalism to study the equilibration of common bar-joint mechanisms with a thermal bath. When the constraints in a mechanism cease to be linearly independent, singularities can appear in its shape space, which is the part of its configuration space after discarding rigid motions. We show that the free-energy landscape of a mechanism at low temperatures is dominated by the neighborhoods of points that correspond to these singularities. We consider two example mechanisms with shape-space singularities and find that they are more likely to be found in configurations near the singularities than others. These findings are expected to help improve the design of nanomechanisms for various applications.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How does the extracellular matrix affect the rigidity of an embedded spheroid?

Cellularized tissue and polymer networks can both transition from floppy to rigid as a function of their control parameters, and, yet, the two systems often mechanically interact, which may affect their respective rigidities. To study this interaction, we consider a vertex model with interfacial tension (a spheroid) embedded in a spring network in two dimensions. We identify two regimes with different global spheroid shapes, governed by the pressure resulting from competition between interfacial tension and tension in the network. In the first regime, the tissue remains compact, while in the second, a cavitation-like instability leads to the emergence of gaps at the tissue-network interface. Intriguingly, compression of the tissue promotes fluidization, while tension promotes cellular alignment and rigidification, with the mechanisms driving rigidification differing on either side of the instability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Loops versus lines and the compression stiffening of cells

Both animal and plant tissue exhibit a nonlinear rheological phenomenon known as compression stiffening, or an increase in moduli with increasing uniaxial compressive strain. Does such a phenomenon exist in single cells, which are the building blocks of tissues? One expects an individual cell to compression soften since the semiflexible biopolymer-based cytoskeletal network maintains the mechanical integrity of the cell and in vitro semiflexible biopolymer networks typically compression soften. To the contrary, we find that mouse embryonic fibroblasts (mEFs) compression stiffen under uniaxial compression via atomic force microscopy (AFM) studies. To understand this finding, we uncover several potential mechanisms for compression stiffening. First, we study a single semiflexible polymer loop modeling the actomyosin cortex enclosing a viscous medium modeled as an incompressible fluid. Second, we study a two-dimensional semiflexible polymer/fiber network interspersed with area-conserving loops, which are a proxy for vesicles and fluid-based organelles. Third, we study two-dimensional fiber networks with angular-constraining crosslinks, i.e. semiflexible loops on the mesh scale. In the latter two cases, the loops act as geometric constraints on the fiber network to help stiffen it via increased angular interactions. We find that the single semiflexible polymer loop model agrees well with our AFM experiments until approximately 35% compressive strain. We also find for the fiber network with area-conserving loops model that the stress-strain curves are sensitive to the packing fraction and size distribution of the area-conserving loops, thereby creating a mechanical fingerprint across different cell types. Finally, we make comparisons between this model and experiments on fibrin networks interlaced with beads as well as discuss the tissue-scale implications of cellular compression stiffening.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Small-scale demixing in confluent biological tissues

Surface tension governed by differential adhesion can drive fluid particle mixtures to sort into separate regions, i.e., demix. Does the same phenomenon occur in confluent biological tissues? We begin to answer this question for epithelial monolayers with a combination of theory via a vertex model and experiments on keratinocyte monolayers. Vertex models are distinct from particle models in that the interactions between the cells are shape-based, as opposed to distance-dependent. We investigate whether a disparity in cell shape or size alone is sufficient to drive demixing in bidisperse vertex model fluid mixtures. Surprisingly, we observe that both types of bidisperse systems robustly mix on large lengthscales. On the other hand, shape disparity generates slight demixing over a few cell diameters, a phenomenon we term micro-demixing. This result can be understood by examining the differential energy barriers for neighbor exchanges (T1 transitions). Experiments with mixtures of wild-type and E-cadherin-deficient keratinocytes on a substrate are consistent with the predicted phenomenon of micro-demixing, which biology may exploit to create subtle patterning. The robustness of mixing at large scales, however, suggests that despite some differences in cell shape and size, progenitor cells can readily mix throughout a developing tissue until acquiring means of recognizing cells of different types.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sponge-like rigid structures in frictional granular packings

We show how rigidity emerges in experiments of sheared frictional granular materials by using generalizations of two methods for identifying rigid structures. Both approaches, the force-based dynamical matrix and the topology-based rigidity percolation, agree with each other and identify similar rigid structures. As the system becomes jammed, at a contact number of $z=2.4\pm 0.1$, a rigid backbone interspersed with floppy, particle-filled holes of a broad range of sizes emerges, creating a sponge-like morphology. We also find that the pressure within rigid structures always exceeds the pressure outside the rigid structures, i.e. that the backbone is load-bearing. These findings shows that it is necessary to go beyond mean-field theory to capture the physics of frictional jamming and also suggests that mechanical stability arises through arch structures and hinges at the mesoscale.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A density-independent glass transition in biological tissues

Cell migration is important in many biological processes, including embryonic development, cancer metastasis, and wound healing. In these tissues, a cell's motion is often strongly constrained by its neighbors, leading to glassy dynamics. While self-propelled particle models exhibit a density-driven glass transition, this does not explain liquid-to-solid transitions in confluent tissues, where there are no gaps between cells and therefore the density is constant. Here we demonstrate the existence of a new type of rigidity transition that occurs in the well-studied vertex model for confluent tissue monolayers at constant density. We find the onset of rigidity is governed by a model parameter that encodes single-cell properties such as cell-cell adhesion and cortical tension, providing an explanation for a liquid-to-solid transitions in confluent tissues and making testable predictions about how these transitions differ from those in particulate matter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rigid cluster decomposition reveals criticality in frictional jamming

We study the nature of the frictional jamming transition within the framework of rigidity percolation theory. Slowly sheared frictional packings are decomposed into rigid clusters and floppy regions with a generalization of the pebble game including frictional contacts. We discover a second-order transition controlled by the emergence of a system-spanning rigid cluster accompanied by a critical cluster size distribution. Rigid clusters also correlate with common measures of rigidity. We contrast this result with frictionless jamming, where the rigid cluster size distribution is noncritical.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shape-shifting droplet networks

We consider a three-dimensional network of aqueous droplets joined by single lipid bilayers to form a cohesive, tissue-like material. The droplets in these networks can be programmed to have distinct osmolarities so that osmotic gradients generate internal stresses via local fluid flows to cause the network to change shape. We discover, using molecular dynamics simulations, a reversible folding-unfolding process by adding an osmotic interaction with the surrounding environment which necessarily evolves dynamically as the shape of the network changes. This discovery is the next important step towards osmotic robotics in this system. We also explore analytically and numerically how the networks become faceted via buckling and how quasi-one-dimensional networks become three-dimensional.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Active elastic dimers: Cells moving on rigid tracks

Experiments suggest that the migration of some cells in the three-dimensional extra cellular matrix bears strong resemblance to one-dimensional cell migration. Motivated by this observation, we construct and study a minimal one-dimensional model cell made of two beads and an active spring moving along a rigid track. The active spring models the stress fibers with their myosin-driven contractility and alpha-actinin-driven extendability, while the friction coefficients of the two beads describe the catch/slip bond behavior of the integrins in focal adhesions. In the absence of active noise, net motion arises from an interplay between active contractility (and passive extendability) of the stress fibers and an asymmetry between the front and back of the cell due to catch bond behavior of integrins at the front of the cell and slip bond behavior of integrins at the back. We obtain reasonable cell speeds with independently estimated parameters. We also study the effects of hysteresis in the active spring, due to catch bond behavior and the dynamics of cross-linking, and the addition of active noise on the motion of the cell. Our model highlights the role of alpha-actinin in three-dimensional cell motility and does not require Arp2/3 actin filament nucleation for net motion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Elastic instabilities in a layered cerebral cortex: A revised axonal tension model for cortex folding

We model the elasticity of the cerebral cortex as a layered material with bending energy along the layers and elastic energy between them in both planar and polar geometries. The cortex is also subjected to axons pulling from the underlying white matter. Above a critical threshold force, a "flat" cortex configuration becomes unstable and periodic unduluations emerge, i.e. a buckling instability occurs. These undulations may indeed initiate folds in the cortex. We identify analytically the critical force and the critical wavelength of the undulations. Both quantities are physiologically relevant values. Our model is a revised version of the axonal tension model for cortex folding, with our version taking into account the layered structure of the cortex. Moreover, our model draws a connection with another competing model for cortex folding, namely the differential growth-induced buckling model. For the polar geometry, we study the relationship between brain size and the critical force and wavelength to understand why small mice brains exhibit no folds, while larger human brains do, for example. Finally, an estimate of the bending rigidity constant for the cortex can be made based on the critical wavelength.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy barriers govern glassy dynamics in tissues

Recent observations demonstrate that densely packed tissues exhibit features of glassy dynamics, such as caging behavior and dynamical heterogeneities, although it has remained unclear how single-cell properties control this behavior. Here we develop numerical and theoretical models to calculate energy barriers to cell rearrangements, which help govern cell migration in cell monolayers. In contrast to work on sheared foams, we find that energy barrier heights are exponentially distributed and depend systematically on the cell's number of neighbors. Based on these results, we predict glassy two-time correlation functions for cell motion, with a timescale that increases rapidly as cell activity decreases. These correlation functions are used to construct simple random walks that reproduce the caging behavior observed for cell trajectories in experiments. This work provides a theoretical framework for predicting collective motion of cells in wound-healing, embryogenesis and cancer tumorigenesis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the modeling of endocytosis in yeast

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The mechanics of anisotropic spring networks

We construct and analyze a model for a disordered linear spring network with anisotropy. The modeling is motivated by, for example, granular systems, nematic elastomers, and ultimately cytoskeletal networks exhibiting some underlying anisotropy. The model consists of a triangular lattice with two different bond occupation probabilities, $p_x$ and $p_y$, for the linear springs. We develop an effective medium theory (EMT) to describe the network elasticity as a function of $p_x$ and $p_y$. We find that the onset of rigidity in the EMT agrees with Maxwell constraint counting. We also find beyond linear behavior in the shear and bulk modulus as a function of occupation probability in the rigid phase for small strains, which differs from the isotropic case. We compare our EMT with numerical simulations to find rather good agreement. Finally, we discuss the implications of extending the reach of effective medium theory as well as draw connections with prior work on both anisotropic and isotropic spring networks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Contact processes in crowded environments

Periodically sheared colloids at low densities demonstrate a dynamical phase transition from an inactive to active phase as the strain amplitude is increased. The inactive phase consists of no collisions/contacts between particles in the steady state limit, while in the active phase collisions persist. To investigate this system at higher densities, we construct and study a conserved-particle-number contact process with novel three-body interactions, which are potentially more likely than two-body interactions at higher densities. For example, consider one active (diffusing) particle colliding with two inactive (non-diffusing) particles such that they become active, in addition to spontaneous inactivation. In mean-field, this system exhibits a continuous dynamical phase transition belonging to the conserved directed percolation universality class. Simulations on square lattices support the mean field result. In contrast, the three-body interaction requiring two active particles to activate one inactive particle exhibits a discontinuous transition. Finally, inspired by kinetically-constrained models of the glass transition, we investigate the "caging effect" at even higher particle densities to look for a second dynamical phase transition back to an inactive phase. Square lattice simulations suggest a continuous transition with a new set of exponents differing from conserved directed percolation, i.e. a new universality class for contact processes with conserved particle number.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Jamming graphs: A local approach to global mechanical rigidity

We revisit the concept of minimal rigidity as applied to soft repulsive, frictionless sphere packings in two-dimensions with the introduction of the jamming graph. Minimal rigidity is a purely combinatorial property encoded via Laman's theorem in two-dimensions. It constrains the global, average coordination number of the graph, for example. However, minimal rigidity does not address the geometry of local mechanical stability. The jamming graph contains both properties of global mechanical stability at the onset of jamming and local mechanical stability. We demonstrate how jamming graphs can be constructed using local moves via the Henneberg construction such that these graphs fall under the jurisdiction of correlated percolation. We then probe how jamming graphs destabilize, or become unjammed, by deleting a bond and computing the resulting rigid cluster distribution. We also study how the system restabilizes with the addition of new contacts and how a jamming graph with extra/redundant contacts destabilizes. The latter endeavor allows us to probe a disc packing in the rigid phase and uncover a potentially new diverging lengthscale associated with the random deletion of contacts as compared to the study of cut-out (or frozen in) subsystems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Vicious accelerating walkers

A vicious walker system consists of N random walkers on a line with any two walkers annihilating each other upon meeting. We study a system of N vicious accelerating walkers with the velocity undergoing Gaussian fluctuations, as opposed to the position. We numerically compute the survival probability exponent, α, for this system, which characterizes the probability for any two walkers not to meet. For example, for N = 3, α = 0.71 \pm 0.01. Based on our numerical data, we conjecture that 1/8N(N - 1) is an upper bound on α. We also numerically study N vicious Levy flights and find, for instance, for N = 3 and a Levy index μ = 1 that α = 1.31 \pm 0.03. Vicious accelerating walkers relate to no-crossing configurations of semiflexible polymer brushes and may prove relevant for a non-Markovian extension of Dyson's Brownian motion model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Comment on: "On the consistency of solutions of the space fractional Schrödinger equation"

In [J. Math. Phys. 53, 042105 (2012)], Bayın claims to prove the consistency of the purported piece-wise solutions to the fractional Schrödinger equation for an infinite square well. However, his calculation uses standard contour integral techniques despite the absence of an analytic integrand. The correct calculation is presented and supports our earlier work proving that the purported piece-wise solutions do not solve the fractional Schrödinger equation for an infinite square well [M. Jeng, S.-L.-Y. Xu, E. Hawkins, and J. M. Schwarz, J. Math. Phys. 51, 062102 (2010)].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Correlated percolation and tricriticality

The recent proliferation of correlated percolation models---models where the addition of edges/vertices is no longer independent of other edges/vertices---has been motivated by the quest to find discontinuous percolation transitions. The leader in this proliferation is what is known as explosive percolation. A recent proof demonstrates that a large class of explosive percolation-type models does not, in fact, exhibit a discontinuous transition[O. Riordan and L. Warnke, Science, {\bf 333}, 322 (2011)]. We, on the other hand, discuss several correlated percolation models, the $k$-core model on random graphs, and the spiral and counter-balance models in two-dimensions, all exhibiting discontinuous transitions in an effort to identify the needed ingredients for such a transition. We then construct mixtures of these models to interpolate between a continuous transition and a discontinuous transition to search for a tricritical point. Using a powerful rate equation approach, we demonstrate that a mixture of $k=2$-core and $k=3$-core vertices on the random graph exhibits a tricritical point. However, for a mixture of $k$-core and counter-balance vertices, heuristic arguments and numerics suggest that there is a line of continuous transitions as the fraction of counter-balance vertices is increased from zero with the line ending at a discontinuous transition only when all vertices are counter-balance. Our results may have potential implications for glassy systems and a recent experiment on shearing a system of frictional particles to induce what is known as jamming.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Level statistics for quantum $k$-core percolation

Quantum $k$-core percolation is the study of quantum transport on $k$-core percolation clusters where each occupied bond must have at least $k$ occupied neighboring bonds. As the bond occupation probability, $p$, is increased from zero to unity, the system undergoes a transition from an insulating phase to a metallic phase. When the lengthscale for the disorder, $l_d$, is much greater than the coherence length, $l_c$, earlier analytical calculations of quantum conduction on the Bethe lattice demonstrate that for $k=3$ the metal-insulator transition (MIT) is discontinuous, suggesting a new universality class of disorder-driven quantum MITs. Here, we numerically compute the level spacing distribution as a function of bond occupation probability $p$ and system size on a Bethe-like lattice. The level spacing analysis suggests that for $k=0$, $p_q$, the quantum percolation critical probability, is greater than $p_c$, the geometrical percolation critical probability, and the transition is continuous. In contrast, for $k=3$, $p_q=p_c$ and the transition is discontinuous such that these numerical findings are consistent with our previous work to reiterate a new universality class of disorder-driven quantum MITs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimal orientation in branched cytoskeletal networks

Actin cytoskeletal protrusions in crawling cells, or lamellipodia, exhibit various morphological properties such as two characteristic peaks in the distribution of filament orientation with respect to the leading edge. To understand these properties, using the dendritic nucleation model as a basis for cytoskeletal restructuring, a kinetic-population model with orientational-dependent branching (birth) and capping (death) is constructed and analyzed. Optimizing for growth yields a relation between the branch angle and filament orientation that explains the two characteristic peaks. The model also exhibits a subdominant population that allows for more accurate modeling of recent measurements of filamentous actin density along the leading edge of lamellipodia in keratocytes. Finally, we explore the relationship between orientational and spatial organization of filamentous actin in lamellipodia and address recent observations of a prevalence of overlapping filaments to branched filaments---a finding that is claimed to be in contradiction with the dendritic nucleation model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Redundancy and cooperativity in the mechanics of compositely crosslinked filamentous networks

The actin cytoskeleton in living cells has many types of crosslinkers. The mechanical interplay between these different crosslinker types is an open issue in cytoskeletal mechanics. We develop a framework to study the cooperativity and redundancy in the mechanics of filamentous networks with two types of crosslinkers: crosslinkers that allow free rotations of filaments and crosslinkers that do not. The framework consists of numerical simulations and an effective medium theory on a percolating triangular lattice. We find that the introduction of angle-constraining crosslinkers significantly lowers the filament concentrations required for these networks to attain mechanical integrity. This cooperative effect also enhances the stiffness of the network and suppresses non-affine deformations at a fixed filament concentration. We further find that semiflexible networks with only freely-rotating crosslinks are mechanically very similar to compositely crosslinked flexible networks with both networks exhibiting the same scaling behavior. We show that the network mechanics can either be redundant or cooperative depending on the relative energy scale of filament bending to the energy stored in the angle-constraining crosslinkers, and the relative concentration of crosslinkers. Our results may have implications for understanding the role of multiple crosslinkers even in a system without bundle formation or other structural motifs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Force network analysis of jammed solids

Using a system of repulsive, soft particles as a model for a jammed solid, we analyze its force network as characterized by the magnitude of the contact force between two particles, the local contact angle subtended between three particles, and the local coordination number. In particular, we measure the local contact angle distribution as a function of the magnitude of the local contact force. We find the suppression of small contact angles for locally larger contact forces, suggesting the existence of chain-like correlations in the locally larger contact forces. We couple this information with a coordination number-spin state mapping to arrive at a Potts spin model with frustration and correlated disorder to draw a potential connection between jammed solids (no quenched disorder) and spin glasses (quenched disorder). We use this connection to measure chaos due to marginality in the jammed system. In addition, we present the replica solution of the one-dimensional, long-range Potts glass as a potential toy building block for a jammed solid, where a sea of weakly interacting spins provide for long-range interactions along a chain-like backbone of more strongly interacting spins.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling the formation of in vitro filopodia

Filopodia are bundles of actin filaments that extend out ahead of the leading edge of a crawling cell to probe its upcoming environment. {\it In vitro} experiments [D. Vignjevic {\it et al.}, J. Cell Biol. {\bf 160}, 951 (2003)] have determined the minimal ingredients required for the formation of filopodia from the dendritic-like morphology of the leading edge. We model these experiments using kinetic aggregation equations for the density of growing bundle tips. In mean field, we determine the bundle size distribution to be broad for bundle sizes smaller than a characteristic bundle size above which the distribution decays exponentially. Two-dimensional simulations incorporating both bundling and cross-linking measure a bundle size distribution that agrees qualitatively with mean field. The simulations also demonstrate a nonmonotonicity in the radial extent of the dendritic region as a function of capping protein concentration, as was observed in experiments, due to the interplay between percolation and the ratcheting of growing filaments off a spherical obstacle.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum $k$-core conduction on the Bethe lattice

Classical and quantum conduction on a bond-diluted Bethe lattice is considered. The bond dilution is subject to the constraint that every occupied bond must have at least $k-1$ neighboring occupied bonds, i.e. $k$-core diluted. In the classical case, we find the onset of conduction for $k=2$ is continuous, while for $k=3$, the onset of conduction is discontinuous with the geometric random first-order phase transition driving the conduction transition. In the quantum case, treating each occupied bond as a random scatterer, we find for $k=3$ that the random first-order phase transition in the geometry also drives the onset of quantum conduction giving rise to a new universality class of Anderson localization transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

On the study of force-balance percolation

We study models of correlated percolation where there are constraints on the occupation of sites that mimic force-balance, i.e. for a site to be stable requires occupied neighboring sites in all four compass directions in two dimensions. We prove rigorously that $p_c<1$ for the two-dimensional models studied. Numerical data indicate that the force-balance percolation transition is discontinuous with a growing crossover length, with perhaps the same form as the jamming percolation models, suggesting the same underlying mechanism driving the transition in both cases. In other words, force-balance percolation and jamming percolation may indeed belong to the same universality class. We find a lower bound for the correlation length in the connected phase and that the correlation function does not appear to be a power law at the transition. Finally, we study the dynamics of the culling procedure invoked to obtain the force-balance configurations and find a dynamical exponent similar to that found in sandpile models.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Comment on ``Jamming Percolation and Glass Transitions in Lattice Models.''

Toninelli, Biroli, and Fisher recently introduced the knights model, a correlated percolation system (Phys. Rev. Lett. 96, 035702 (2006)). They claimed to prove that the critical point of this model was the same as that for directed percolation, and then used this to show (assuming a conjecture about directed percolation, described later in a more detailed proof) that this model has a discontinuous phase transition with a diverging correlation length. However, there is an error in their work, so that these results are unproven for their model. Their proofs can, however, be modified to work for a similar model.

0 citations0 reviewsNo code linkedOpen access