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

David R. Nelson

David R. Nelson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softPopulations and EvolutionBiological Physicscond-mat.mtrl-scicond-matnlin.CDphysics.flu-dynCell Behaviorcond-mat.dis-nncond-mat.mes-hallcond-mat.supr-conBiomoleculescond-mat.quant-gasmath-phmath.MPNeurons and CognitionSubcellular Processessupr-con

Catalog footprint

What is connected

47works
19topics
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

47 published item(s)

preprint2022arXiv

Anomalous thermal expansion in Ising-like puckered sheets

Motivated by efforts to create thin nanoscale metamaterials and understand atomically thin binary monolayers, we study the finite temperature statistical mechanics of arrays of bistable buckled dilations embedded in free-standing two-dimensional crystalline membranes that are allowed to fluctuate in three dimensions. The buckled nodes behave like discrete, but highly compressible, Ising spins, leading to a phase transition at $T_c$ with singularities in the staggered "magnetization," susceptibility, and specific heat, studied via molecular dynamics simulations. Unlike conventional Ising models, we observe a striking divergence and sign change of the coefficient of thermal expansion near $T_c$ caused by the coupling of flexural phonons to the buckled spin texture. We argue that a phenomenological model coupling Ising degrees of freedom to the flexural phonons in a thin elastic sheet can explain this unusual response.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Defect absorption and emission for $p$-atic liquid crystals on cones

We investigate the ground state configurations of $p$-atic liquid crystals on fixed curved surfaces. We focus on the intrinsic geometry and show that isothermal coordinates are particularly convenient as they explicitly encode a geometric contribution to the elastic potential. In the special case of a cone with half-angle $β$, the apex develops an effective topological charge of $-χ$, where $2πχ= 2π(1-\sinβ)$ is the deficit angle of the cone, and a topological defect of charge $σ$ behaves as if it had an effective topological charge $Q_\mathrm{eff} = (σ- σ^2/2)$ when interacting with the apex. The effective charge of the apex leads to defect absorption and emission at the cone apex as the deficit angle of the cone is varied. For total topological defect charge 1, e.g. imposed by tangential boundary conditions at the edge, we find that for a disk the ground state configuration consists of $p$ defects each of charge $+1/p$ lying equally spaced on a concentric ring of radius $d = (\frac{p-1}{3p-1})^{\frac{1}{2p}} R$, where $R$ is the radius of the disk. In the case of a cone with tangential boundary conditions at the base, we find three types of ground state configurations as a function of cone angle: (1) for sharp cones, all of the $+1/p$ defects are absorbed by the apex; (2) at intermediate cone angles, some of the $+1/p$ defects are absorbed by the apex and the rest lie equally spaced along a concentric ring on the flank; and (3) for nearly flat cones, all of the $+1/p$ defects lie equally spaced along a concentric ring on the flank. Here the defect positions and the absorption transitions depend intricately on $p$ and the deficit angle which we analytically compute. We check these results with numerical simulations for a set of commensurate cone angles and find excellent agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fractional defect charges in $p$-atic liquid crystals on cones

Conical surfaces, with a delta function of Gaussian curvature at the apex, are perhaps the simplest example of geometric frustration. We study two-dimensional liquid crystals with $p$-fold rotational symmetry ($p$-atics) on the surfaces of cones. For free boundary conditions at the base, we find both the ground state(s) and a discrete ladder of metastable states as a function of both the cone angle and the liquid crystal symmetry $p$. We find that these states are characterized by a set of fractional defect charges at the apex and that the ground states are in general frustrated due to effects of parallel transport along the azimuthal direction of the cone. We check our predictions for the ground state energies numerically for a set of commensurate cone angles (corresponding to a set of commensurate Gaussian curvatures concentrated at the cone apex), whose surfaces can be polygonized as a perfect triangular or square mesh, and find excellent agreement with our theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spatial population genetics with fluid flow

The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2d), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1d and 2d flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how these genetic drops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Phonon eigenfunctions of inhomogeneous lattices: Can you hear the shape of a cone?

We study the phonon modes of interacting particles on the surface of a truncated cone resting on a plane subject to gravity, inspired by recent colloidal experiments. We derive the ground state configuration of the particles under gravitational pressure in the small cone angle limit, and find an inhomogeneous triangular lattice with spatially varying density but robust local order. The inhomogeneity has striking effects on the normal modes such that an important feature of the cone geometry, namely its apex angle, can be extracted from the lattice excitations. The shape of the cone leads to energy crossings at long wavelengths and frequency-dependent quasi-localization at short wavelengths. We analytically derive the localization domain boundaries of the phonons in the limit of small cone angle and check our results with numerical results for eigenfunctions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Buckling and metastability in membranes with dilation arrays

We study periodic arrays of impurities that create localized regions of expansion, embedded in two-dimensional crystalline membranes. These arrays provide a simple elastic model of shape memory. As the size of each dilational impurity increases (or the relative cost of bending to stretching decreases), it becomes energetically favorable for each of the $M$ impurities to buckle up or down into the third dimension, thus allowing for of order $2^M$ metastable surface configurations corresponding to different impurity "spin" configurations. With both discrete simulations and the nonlinear continuum theory of elastic plates, we explore the buckling of both isolated dilations and dilation arrays at zero temperature, guided by analogies with Ising antiferromagnets. We conjecture ground states for systems with triangular and square impurity superlattices, and comment briefly on the possible behaviors at finite temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Statistical Mechanics of Dislocation Pileups in Two Dimensions

Dislocation pileups directly impact the material properties of crystalline solids through the arrangement and collective motion of interacting dislocations. We study the statistical mechanics of these ordered defect structures embedded in two dimensional crystals, where the dislocations themselves form one-dimensional lattices. In particular, pileups exemplify a new class of inhomogeneous crystals characterized by spatially varying lattice spacings. By analytically formulating key statistical quantities and comparing our theory to numerical experiments using an intriguing mapping of dislocation positions onto the eigenvalues of recently studied random matrix ensembles, we uncover two types of one-dimensional phase transitions in dislocation pileups: a thermal depinning transition out of long-range translational order from the pinned-defect phase, due to a periodic Peierls potential, to a floating-defect state, and finally the melting out of a quasi-long range ordered floating defect-solid phase to a defect-liquid. We also find the set of transition temperatures at which these transitions can be directly observed through the one-dimensional structure factor, where the delta function Bragg peaks, at the pinned-defect to floating-defect transition, broaden into algebraically diverging Bragg peaks, which then sequentially disappear as one approaches the two-dimensional melting transition of the host crystal. We calculate a set of temperature-dependent critical exponents for the structure factor and radial distribution function, and obtain their exact forms for both uniform and inhomogeneous pileups using random matrix theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evolutionary Dynamics with Fluctuating Population Sizes and Strong Mutualism

Game theory ideas provide a useful framework for studying evolutionary dynamics in a well-mixed environment. This approach, however, typically enforces a strictly fixed overall population size, deemphasizing natural growth processes. We study a competitive Lotka-Volterra model, with number fluctuations, that accounts for natural population growth and encompasses interaction scenarios typical of evolutionary games. We show that, in an appropriate limit, the model describes standard evolutionary games with both genetic drift and overall population size fluctuations. However, there are also regimes where a varying population size can strongly influence the evolutionary dynamics. We focus on the strong mutualism scenario and demonstrate that standard evolutionary game theory fails to describe our simulation results. We then analytically and numerically determine fixation probabilities as well as mean fixation times using matched asymptotic expansions, taking into account the population size degree of freedom. These results elucidate the interplay between population dynamics and evolutionary dynamics in well-mixed systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Non-Hermitian Localization in Biological Networks

We explore the spectra and localization properties of the N-site banded one-dimensional non-Hermitian random matrices that arise naturally in sparse neural networks. Approximately equal numbers of random excitatory and inhibitory connections lead to spatially localized eigenfunctions, and an intricate eigenvalue spectrum in the complex plane that controls the spontaneous activity and induced response. A finite fraction of the eigenvalues condense onto the real or imaginary axes. For large N, the spectrum has remarkable symmetries not only with respect to reflections across the real and imaginary axes, but also with respect to 90 degree rotations, with an unusual anisotropic divergence in the localization length near the origin. When chains with periodic boundary conditions become directed, with a systematic directional bias superimposed on the randomness, a hole centered on the origin opens up in the density-of-states in the complex plane. All states are extended on the rim of this hole, while the localized eigenvalues outside the hole are unchanged. The bias dependent shape of this hole tracks the bias independent contours of constant localization length. We treat the large-N limit by a combination of direct numerical diagonalization and using transfer matrices, an approach that allows us to exploit an electrostatic analogy connecting the "charges" embodied in the eigenvalue distribution with the contours of constant localization length. We show that similar results are obtained for more realistic neural networks that obey "Dale's Law" (each site is purely excitatory or inhibitory), and conclude with perturbation theory results that describe the limit of large bias g, when all states are extended. Related problems arise in random ecological networks and in chains of artificial cells with randomly coupled gene expression patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Striated Populations in Disordered Environments with Advection

Growth in static and controlled environments such as a Petri dish can be used to study the spatial population dynamics of microorganisms. However, natural populations such as marine microbes experience fluid advection and often grow up in heterogeneous environments. We investigate a generalized Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation describing single species population subject to a constant flow field and quenched random spatially inhomogeneous growth rates with a fertile overall growth condition. We analytically and numerically demonstrate that the non-equilibrium steady-state population density develops a flow-driven striation pattern. The striations are highly asymmetric with a longitudinal correlation length that diverges linearly with the flow speed and a transverse correlation length that approaches a finite velocity-independent value. Linear response theory is developed to study the statistics of the steady states. Theoretical predictions show excellent agreement with the numerical steady states of the generalized FKPP equation obtained from Lattice Boltzmann simulations. These findings suggest that, although the growth disorder can be spatially uncorrelated, correlated population structures with striations emerge naturally at sufficiently strong advection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Response of thermalized ribbons to pulling and bending

Motivated by recent free-standing graphene experiments, we show how thermal fluctuations affect the mechanical properties of microscopically thin solid ribbons. A renormalization group analysis of flexural phonons reveals that elongated ribbons behave like highly anisotropic polymers, where the two dimensional nature of ribbons is reflected in non-trivial power law scalings of the persistence length and effective bending and twisting rigidities with the ribbon width. With a coarse-grained transfer matrix approach, we then show how thermalized ribbons respond to pulling and bending forces over a wide spectrum of temperatures, forces and ribbon lengths.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Survival Probabilities at Spherical Frontiers

Motivated by tumor growth and spatial population genetics, we study the interplay between evolutionary and spatial dynamics at the surfaces of three-dimensional, spherical range expansions. We consider range expansion radii that grow with an arbitrary power-law in time: $R(t)=R_0(1+t/t^*)^Θ$, where $Θ$ is a growth exponent, $R_0$ is the initial radius, and $t^*$ is a characteristic time for the growth, to be affected by the inflating geometry. We vary the parameters $t^*$ and $Θ$ to capture a variety of possible growth regimes. Guided by recent results for two-dimensional inflating range expansions, we identify key dimensionless parameters that describe the survival probability of a mutant cell with a small selective advantage arising at the population frontier. Using analytical techniques, we calculate this probability for arbitrary $Θ$. We compare our results to simulations of linearly inflating expansions ($Θ=1$ spherical Fisher-Kolmogorov-Petrovsky-Piscunov waves) and treadmilling populations ($Θ=0$, with cells in the interior removed by apoptosis or a similar process). We find that mutations at linearly inflating fronts have survival probabilities enhanced by factors of 100 or more relative to mutations at treadmilling population frontiers. We also discuss the special properties of "marginally inflating" $(Θ=1/2)$ expansions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Asymmetric Mutualism in Two- and Three-Dimensional Range Expansions

Genetic drift at the frontiers of two-dimensional range expansions of microorganisms can frustrate local cooperation between different genetic variants, demixing the population into distinct sectors. In a biological context, mutualistic or antagonistic interactions will typically be asymmetric between variants. By taking into account both the asymmetry and the interaction strength, we show that the much weaker demixing in three dimensions allows for a mutualistic phase over a much wider range of asymmetric cooperative benefits, with mutualism prevailing for any positive, symmetric benefit. We also demonstrate that expansions with undulating fronts roughen dramatically at the boundaries of the mutualistic phase, with severe consequences for the population genetics along the transition lines.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bending forces plastically deform growing bacterial cell walls

Cell walls define a cell shape in bacteria. They are rigid to resist large internal pressures, but remarkably plastic to adapt to a wide range of external forces and geometric constraints. Currently, it is unknown how bacteria maintain their shape. In this work, we develop experimental and theoretical approaches and show that mechanical stresses regulate bacterial cell-wall growth. By applying a precisely controllable hydrodynamic force to growing rod-shaped Escherichia coli and Bacillus subtilis cells, we demonstrate that the cells can exhibit two fundamentally different modes of deformation. The cells behave like elastic rods when subjected to transient forces, but deform plastically when significant cell wall synthesis occurs while the force is applied. The deformed cells always recover their shape. The experimental results are in quantitative agreement with the predictions of the theory of dislocation-mediated growth. In particular, we find that a single dimensionless parameter, which depends on a combination of independently measured physical properties of the cell, can describe the cell's responses under various experimental conditions. These findings provide insight into how living cells robustly maintain their shape under varying physical environments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum hexatic order in two-dimensional dipolar and charged fluids

Recent advances in cold atom experimentation suggest that studies of quantum two-dimensional melting of dipolar molecules, with dipoles aligned perpendicular to ordering plane, may be on the horizon. An intriguing aspect of this problem is that two-dimensional \emph{classical} aligned dipoles (already studied in great detail in soft matter experiments on magnetic colloids) are known to melt via a two-stage process, with an intermediate hexatic phase separating the usual crystal and isotropic fluid phases. We estimate here the effect of quantum fluctuations on this hexatic phase, for both dipolar systems and charged Wigner crystals. Our approximate phase diagrams rely on a pair of Lindemann criteria, suitably adapted to deal with effects of thermal fluctuations in two dimensions. As part of our analysis, we determine the phonon spectra of quantum particles on a triangular lattice interacting with repulsive $1/r^3$ and $1/r$ potentials. A large softening of the transverse and longitudinal phonon frequencies, due to both lattice effects and quantum fluctuations, plays a significant role in our analysis. The hexatic phase is predicted to survive down to very low temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Defects on cylinders: superfluid helium films and bacterial cell walls

There is a deep analogy between the physics of crystalline solids and the behavior of superfluids, dating back to pioneering work of Phillip Anderson, Paul Martin and others. The stiffness to shear deformations in a periodic crystal resembles the superfluid density that controls the behavior of supercurrents in neutral superfluids such as He^4. Dislocations in solids have a close analogy with quantized vortices in superfluids. Remarkable recent experiments on the way rod-shaped bacteria elongate their cell walls have focused attention on the dynamics and interactions of point-like dislocation defects in partially ordered cylindrical crystalline monolayers. In these lectures, we review the physics of superfluid helium films on cylinders and discuss how confinement in one direction affects vortex interactions with supercurrents. Although there are similarities with the way dislocations respond to strains on cylinders, important differences emerge, due to the vector nature of the topological charges characterizing the dislocations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Elastic Platonic Shells

On microscopic scales, the crystallinity of flexible tethered or cross linked membranes determines their mechanical response. We show that by controlling the type, number and distribution of defects on a spherical elastic shell, it is possible to direct the morphology of these structures. Our numerical simulations show that by deflating a crystalline shell with defects, we can create elastic shell analogs of the classical Platonic solids. These morphologies arise via a sharp buckling transition from the sphere which is strongly hysteretic in loading-unloading. We construct a minimal Landau theory for the transition using quadratic and cubic invariants of the spherical harmonic modes. Our approach suggests methods to engineer shape into soft spherical shells using a frozen defect topology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mechanical Properties of Warped Membranes

We explore how a frozen background metric affects the mechanical properties of planar membranes with a shear modulus. We focus on a special class of "warped membranes" with a preferred random height profile characterized by random Gaussian variables $h({\bf q})$ in Fourier space with zero mean and variance $<|h({\bf q})|^2> \sim q^{-d_h}$ and show that in the linear response regime the mechanical properties depend dramatically on the system size $L$ for $d_h \ge 2$. Membranes with $d_h=4$ could be produced by flash polymerization of lyotropic smectic liquid crystals. Via a self consistent screening approximation we find that the renormalized bending rigidity increases as $κ_R \sim L^{(d_h-2)/2}$ for membranes of size $L$, while the Young and shear modulii decrease according to $Y_R,\ μ_R \sim L^{-(d_h-2)/2}$ resulting in a universal Poisson ratio. Numerical results show good agreement with analytically determined exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nutrient Shielding in Clusters of Cells

Cellular nutrient consumption is influenced by both the nutrient uptake kinetics of an individual cell and the cells' spatial arrangement. Large cell clusters or colonies have inhibited growth at the cluster's center due to the shielding of nutrients by the cells closer to the surface. We develop an effective medium theory that predicts a thickness $\ell$ of the outer shell of cells in the cluster that receives enough nutrient to grow. The cells are treated as partially absorbing identical spherical nutrient sinks, and we identify a dimensionless parameter $ν$ that characterizes the absorption strength of each cell. The parameter $ν$ can vary over many orders of magnitude between different cell types, ranging from bacteria and yeast to human tissue. The thickness $\ell$ decreases with increasing $ν$, increasing cell volume fraction $ϕ$, and decreasing ambient nutrient concentration $ψ_{\infty}$. The theoretical results are compared with numerical simulations and experiments. In the latter studies, colonies of budding yeast, \textit{Saccharomyces cerevisiae}, are grown on glucose media and imaged under a confocal microscope. We measure the growth inside the colonies via a fluorescent protein reporter and compare the experimental and theoretical results for the thickness $\ell$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radial Domany-Kinzel Models with Mutation and Selection

We study the effect of spatial structure, genetic drift, mutation, and selective pressure on the evolutionary dynamics in a simplified model of asexual organisms colonizing a new territory. Under an appropriate coarse-graining, the evolutionary dynamics is related to the directed percolation processes that arise in voter models, the Domany-Kinzel (DK) model, contact process, etc. We explore the differences between linear (flat front) expansions and the much less familiar radial (curved front) range expansions. For the radial expansion, we develop a generalized, off-lattice DK model that minimizes otherwise persistent lattice artifacts. With both simulations and analytical techniques, we study the survival probability of advantageous mutants, the spatial correlations between domains of neutral strains, and the dynamics of populations with deleterious mutations. "Inflation" at the frontier leads to striking differences between radial and linear expansions. For a colony with initial radius $R_0$ expanding at velocity $v$, significant genetic demixing, caused by local genetic drift, only occurs up to a finite time $t^* = R_0/v$, after which portions of the colony become causally disconnected due to the inflating perimeter of the expanding front. As a result, the effect of a selective advantage is amplified relative to genetic drift, increasing the survival probability of advantageous mutants. Inflation also modifies the underlying directed percolation transition, introducing novel scaling functions and modifications similar to a finite size effect. Finally, we consider radial range expansions with deflating perimeters, as might arise from colonization initiated along the shores of an island.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spinodal decomposition in homogeneous and isotropic turbulence

We study the competition between domain coarsening in a symmetric binary mixtures below the critical temperature and turbulent fluctuations. We find that the coarsening process is arrested in presence of turbulence. The physics of the process shares remarkable similarities with the behaviour of diluted turbulent emulsions and the arrest length scale can be estimated with an argument similar to the one proposed by Kolmogorov and Hinze for the maximal stability diameter of droplets in turbulence. Although in the absence of flow the microscopic diffusion constant is negative, turbulence does effectively arrest the inverse cascade of concentration fluctuations by making the low wavelength diffusion constant positive for scales above the Hinze length.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Theory of interacting dislocations on cylinders

We study the mechanics and statistical physics of dislocations interacting on cylinders, motivated by the elongation of rod-shaped bacterial cell walls and cylindrical assemblies of colloidal particles subject to external stresses. The interaction energy and forces between dislocations are solved analytically, and analyzed asymptotically. The results of continuum elastic theory agree well with numerical simulations on finite lattices even for relatively small systems. Isolated dislocations on a cylinder act like grain boundaries. With colloidal crystals in mind, we show that saddle points are created by a Peach-Koehler force on the dislocations in the circumferential direction, causing dislocation pairs to unbind. The thermal nucleation rate of dislocation unbinding is calculated, for an arbitrary mobility tensor and external stress, including the case of a twist-induced Peach-Koehler force along the cylinder axis. Surprisingly rich phenomena arise for dislocations on cylinders, despite their vanishing Gaussian curvature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermal Excitations of Warped Membranes

We explore thermal fluctuations of thin planar membranes with a frozen spatially-varying background metric and a shear modulus. We focus on a special class of $D$-dimensional ``warped membranes'' embedded in a $d-$dimensional space with $d\ge D+1$ and a preferred height profile characterized by quenched random Gaussian variables $\{h_α({\bf q})\}$, $α=D+1,\ldots, d$, in Fourier space with zero mean and a power law variance $\overline{ h_α({\bf q}_1) h_β({\bf q}_2) } \sim δ_{α, β} \, δ_{{\bf q}_1, -{\bf q}_2} \, q_1^{-d_h}$. The case $D=2$, $d=3$ with $d_h = 4$ could be realized by flash polymerizing lyotropic smectic liquid crystals. For $D < \max\{4, d_h\}$ the elastic constants are non-trivially renormalized and become scale dependent. Via a self consistent screening approximation we find that the renormalized bending rigidity increases for small wavevectors ${\bf q}$ as $κ_R \sim q^{-η_f}$, while the in-hyperplane elastic constants decrease according to $λ_R,\ μ_R \sim q^{+η_u}$. The quenched background metric is relevant (irelevant) for warped membranes characterized by exponent $d_h > 4 - η_f^{(F)}$ ($d_h < 4 - η_f^{(F)}$), where $η_f^{(F)}$ is the scaling exponent for tethered surfaces with a flat background metric, and the scaling exponents are related through $η_u + η_f = d_h - D$ ($η_u + 2 η_f=4-D$).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dislocation-mediated growth of bacterial cell walls

Recent experiments have illuminated a remarkable growth mechanism of rod-shaped bacteria: proteins associated with cell wall extension move at constant velocity in circles oriented approximately along the cell circumference (Garner et al., Science (2011), Dominguez-Escobar et al. Science (2011), van Teeffelen et al. PNAS (2011). We view these as dislocations in the partially ordered peptidoglycan structure, activated by glycan strand extension machinery, and study theoretically the dynamics of these interacting defects on the surface of a cylinder. Generation and motion of these interacting defects lead to surprising effects arising from the cylindrical geometry, with important implications for growth. We also discuss how long range elastic interactions and turgor pressure affect the dynamics of the fraction of actively moving dislocations in the bacterial cell wall.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fluctuating shells under pressure

Thermal fluctuations strongly modify the large length-scale elastic behavior of crosslinked membranes, giving rise to scale-dependent elastic moduli. While thermal effects in flat membranes are well understood, many natural and artificial microstructures are modeled as thin elastic {\it shells}. Shells are distinguished from flat membranes by their nonzero curvature, which provides a size-dependent coupling between the in-plane stretching modes and the out-of-plane undulations. In addition, a shell can support a pressure difference between its interior and exterior. Little is known about the effect of thermal fluctuations on the elastic properties of shells. Here, we study the statistical mechanics of shape fluctuations in a pressurized spherical shell using perturbation theory and Monte Carlo computer simulations, explicitly including the effects of curvature and an inward pressure. We predict novel properties of fluctuating thin shells under point indentations and pressure-induced deformations. The contribution due to thermal fluctuations increases with increasing ratio of shell radius to thickness, and dominates the response when the product of this ratio and the thermal energy becomes large compared to the bending rigidity of the shell. Thermal effects are enhanced when a large uniform inward pressure acts on the shell, and diverge as this pressure approaches the classical buckling transition of the shell. Our results are relevant for the elasticity and osmotic collapse of microcapsules.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Growth, competition and cooperation in spatial population genetics

We study an individual based model describing competition in space between two different alleles. Although the model is similar in spirit to classic models of spatial population genetics such as the stepping stone model, here however space is continuous and the total density of competing individuals fluctuates due to demographic stochasticity. By means of analytics and numerical simulations, we study the behavior of fixation probabilities, fixation times, and heterozygosity, in a neutral setting and in cases where the two species can compete or cooperate. By concluding with examples in which individuals are transported by fluid flows, we argue that this model is a natural choice to describe competition in marine environments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How Does A Porous Shell Collapse? Delayed Buckling And Guided Folding Of Inhomogeneous Capsules

Colloidal capsules can sustain an external osmotic pressure; however, for a sufficiently large pressure, they will ultimately buckle. This process can be strongly influenced by structural inhomogeneities in the capsule shells. We explore how the time delay before the onset of buckling decreases as the shells are made more inhomogeneous; this behavior can be quantitatively understood by coupling shell theory with Darcy's law. In addition, we show that the shell inhomogeneity can dramatically change the folding pathway taken by a capsule after it buckles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Population dynamics in compressible flows

Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible turbulence is relevant for population dynamics in the ocean and we consider cases both where the velocity field is turbulent and when it is static. Furthermore, we investigate populations in terms of a continuos density field and when the populations are treated via discrete particles. In the last case we focus on the competition and fixation of one species compared to another

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Quantitative Test of Population Genetics Using Spatio-Genetic Patterns in Bacterial Colonies

It is widely accepted that population genetics theory is the cornerstone of evolutionary analyses. Empirical tests of the theory, however, are challenging because of the complex relationships between space, dispersal, and evolution. Critically, we lack quantitative validation of the spatial models of population genetics. Here we combine analytics, on and off-lattice simulations, and experiments with bacteria to perform quantitative tests of the theory. We study two bacterial species, the gut microbe Escherichia coli and the opportunistic pathogen Pseudomonas aeruginosa, and show that spatio-genetic patterns in colony biofilms of both species are accurately described by an extension of the one-dimensional stepping-stone model. We use one empirical measure, genetic diversity at the colony periphery, to parameterize our models and show that we can then accurately predict another key variable: the degree of short-range cell migration along an edge. Moreover, the model allows us to estimate other key parameters including effective population size (density) at the expansion frontier. While our experimental system is a simplification of natural microbial community, we argue it is a proof of principle that the spatial models of population genetics can quantitatively capture organismal evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Competition and cooperation in one-dimensional stepping stone models

Cooperative mutualism is a major force driving evolution and sustaining ecosystems. Although the importance of spatial degrees of freedom and number fluctuations is well-known, their effects on mutualism are not fully understood. With range expansions of microbes in mind, we show that, even when mutualism confers a distinct selective advantage, it persists only in populations with high density and frequent migrations. When these parameters are reduced, mutualism is generically lost via a directed percolation process, with a phase diagram strongly influenced by an exceptional DP2 transition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cumulative compressibility effects on population dynamics in turbulent flows

Bacteria and plankton populations living in oceans and lakes reproduce and die under the in- fluence of turbulent currents. Turbulent transport interacts in a complex way with the dynamics of populations because the typical reproduction time of microorganism is within the inertial range of turbulent time scales. In the present manuscript we quantitatively investigate the effect of flow compressibility on the dynamics of populations. While a small compressibility can be induced by several physical mechanisms, like density mismatch or the finite size of microorganisms with respect to the fluid turbulence, its effect on the the carrying capacity of the ecosystem can be dramatic. We report, for the first time, how a small compressibility can produce a sizeable reduction in the carrying capacity, due to an integrated effect made possible by the long replication times of the organisms with respect to turbulent time scales. A full statistical quantification of the fluctuations of population concentration field leads to data collapse over a broad range in parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Genetic Demixing and Evolutionary Forces in the One-Dimensional Stepping Stone Model

We review and extend results for mutation, selection, genetic drift, and migration in a one-dimensional continuous population. The population is described by a continuous limit of the stepping stone model, which leads to the stochastic Fisher-Kolmogorov-Petrovsky-Piscounov equation with additional terms describing mutations. Although the stepping stone model was first proposed for population genetics, it is closely related to "voter models" of interest in nonequilibrium statistical mechanics. The stepping stone model can also be regarded as an approximation to the dynamics of a thin layer of actively growing pioneers at the frontier of a colony of microorganisms undergoing a range expansion on a Petri dish. We find that the population tends to segregate into monoallelic domains. This segregation slows down genetic drift and selection because these two evolutionary forces can only act at the boundaries between the domains; the effects of mutation, however, are not significantly affected by the segregation. Although fixation in the neutral well-mixed (or "zero dimensional") model occurs exponentially in time, it occurs only algebraically fast in the one-dimensional model. If selection is weak, selective sweeps occur exponentially fast in both well-mixed and one-dimensional populations, but the time constants are different. We also find an unusual sublinear increase in the variance of the spatially averaged allele frequency with time. Although we focus on two alleles or variants, q-allele Potts-like models of gene segregation are considered as well. We also investigate the effects of geometry at the frontier by considering growth of circular colonies. Our analytical results are checked with simulations, and could be tested against recent spatial experiments on range expansions off linear inoculations of Escherichia coli and Saccharomyces cerevisiae.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Population genetics in compressible flows

We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We then study three examples of compressible advecting fields: a shell model of turbulence, a sinusoidal velocity field and a linear velocity sink. In all cases, advection leads to a striking drop in the fixation time, as well as a large reduction in the global carrying capacity. Despite localization on convergence zones, one species goes extinct much more rapidly than in well-mixed populations. For a weak harmonic potential, one finds a bimodal distribution of fixation times. The long-lived states in this case are demixed configurations with a single boundary, whose location depends on the fitness advantage.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Population dynamics at high Reynolds number

We study the statistical properties of population dynamics evolving in a realistic two-dimensional compressible turbulent velocity field. We show that the interplay between turbulent dynamics and population growth and saturation leads to quasi-localization and a remarkable reduction in the carrying capacity. The statistical properties of the population density are investigated and quantified via multifractal scaling analysis. We also investigate numerically the singular limit of negligibly small growth rates and delocalization of population ridges triggered by uniform advection.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structure and dynamics of topological defects in a glassy liquid on a negatively curved manifold

We study the low-temperature regime of an atomic liquid on the hyperbolic plane by means of molecular dynamics simulation and we compare the results to a continuum theory of defects in a negatively curved hexagonal background. In agreement with the theory and previous results on positively curved (spherical) surfaces, we find that the atomic configurations consist of isolated defect structures, dubbed "grain boundary scars", that form around an irreducible density of curvature-induced disclinations in an otherwise hexagonal background. We investigate the structure and the dynamics of these grain boundary scars.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fisher equation with turbulence in one dimension

We investigate the dynamics of the Fisher equation for the spreading of micro-organisms in one dimension subject to both turbulent convection and diffusion. We show that for strong enough turbulence, bacteria, for example, track in a quasilocalized fashion (with remakably long persistance times) sinks in the turbulent field. An important consequence is a large reduction in the carrying capacity of the fluid medium. We determine analytically the regimes where this quasi-localized behavior occurs and test our predictions by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Collapse and folding of pressurized rings in two dimensions

Hydrostatically pressurized circular rings confined to two dimensions (or cylinders constrained to have only z-independent deformations) undergo Euler type buckling when the outside pressure exceeds a critical value. We perform a stability analysis of rings with arc-length dependent bending moduli and determine how weakened bending modulus segments affect the buckling critical pressure. Rings with a 4-fold symmetric modulation are particularly susceptible to collapse. In addition we study the initial post-buckling stages of the pressurized rings to determine possible ring folding patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Vortices in Superfluid Films on Curved Surfaces

We present a systematic study of how vortices in superfluid films interact with the spatially varying Gaussian curvature of the underlying substrate. The Gaussian curvature acts as a source for a geometric potential that attracts (repels) vortices towards regions of negative (positive) Gaussian curvature independently of the sign of their topological charge. Various experimental tests involving rotating superfluid films and vortex pinning are first discussed for films coating gently curved substrates that can be treated in perturbation theory from flatness. An estimate of the experimental regimes of interest is obtained by comparing the strength of the geometrical forces to the vortex pinning induced by the varying thickness of the film which is in turn caused by capillary effects and gravity. We then present a non-perturbative technique based on conformal mappings that leads an exact solution for the geometric potential as well as the geometric correction to the interaction between vortices. The conformal mapping approach is illustrated by means of explicit calculations of the geometric effects encountered in the study of some strongly curved surfaces and by deriving universal bounds on their strength.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Geometric Theory of Columnar Phases on Curved Substrates

We study thin self-assembled columns constrained to lie on a curved, rigid substrate. The curvature presents no local obstruction to equally spaced columns in contrast to curved crystals for which the crystalline bonds are frustrated. Instead, the vanishing compressional strain of the columns implies that their normals lie on geodesics which converge (diverge) in regions of positive (negative) Gaussian curvature, in analogy to the focussing of light rays by a lens. We show that the out of plane bending of the cylinders acts as an effective ordering field.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Unzipping Vortices in Type-II Superconductors

The unzipping of vortex lines using magnetic-force microscopy from extended defects is studied theoretically. We study both the unzipping isolated vortex from common defects, such as columnar pins and twin-planes, and the unzipping of a vortex from a plane in the presence of other vortices. We show, using analytic and numerical methods, that the universal properties of the unzipping transition of a single vortex depend only on the dimensionality of the defect in the presence and absence of disorder. For the unzipping of a vortex from a plane populated with many vortices is shown to be very sensitive to the properties of the vortices in the two-dimensional plane. In particular such unzipping experiments can be used to measure the ``Luttinger liquid parameter'' of the vortices in the plane. In addition we suggest a method for measuring the line tension of the vortex directly using the experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Sequence heterogeneity and the dynamics of molecular motors

The effect of sequence heterogeneity on the dynamics of molecular motors is reviewed and analyzed using a set of recently introduced lattice models. First, we review results for the influence of heterogenous tracks such as a single-strand of DNA or RNA on the dynamics of the motors. We stress how the predicted behavior might be observed experimentally in anomalous drift and diffusion of motors over a wide range of parameters near the stall force and discuss the extreme limit of strongly biased motors with one-way hopping. We then consider the dynamics in an environment containing a variety of different fuels which supply chemical energy for the motor motion, either on a heterogeneous or on a periodic track. The results for motion along a periodic track are relevant to kinesin motors in a solution with a mixture of different nucleotide triphosphate fuel sources.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Non-Hermitian Luttinger Liquids and Vortex Physics

As a model of two thermally excited flux liquids connected by a weak link, we study the effect of a single line defect on vortex filaments oriented parallel to the surface of a thin planar superconductor. When the applied field is tilted relative to the line defect, the physics is described by a nonhermitian Luttinger liquid of interacting quantum bosons in one spatial dimension with a point defect. We analyze this problem using a combination of analytic and numerical density matrix renormalization group methods, uncovering a delicate interplay between enhancement of pinning due to Luttinger liquid effects and depinning due to the tilted magnetic field. Interactions dramatically improve the ability of a single columnar pin to suppress vortex tilt when the Luttinger liquid parameter g is less than or equal to one.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Towards a Tetravalent Chemistry of Colloids

We propose coating spherical particles or droplets with anisotropic nano-sized objects to allow micron-scale colloids to link or functionalize with a four-fold valence, similar to the sp3 hybridized chemical bonds associated with, e.g., carbon, silicon and germanium. Candidates for such coatings include triblock copolymers, gemini lipids, metallic or semiconducting nanorods and conventional liquid crystal compounds. We estimate the size of the relevant nematic Frank constants, discuss how to obtain other valences and analyze the thermal distortions of ground state configurations of defects on the sphere.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Polymer Winding Numbers and Quantum Mechanics

The winding of a single polymer in thermal equilibrium around a repulsive cylindrical obstacle is perhaps the simplest example of statistical mechanics in a multiply connected geometry. As shown by S.F. Edwards, this problem is closely related to the quantum mechanics of a charged particle interacting with a Aharonov-Bohm flux. In another development, Pollock and Ceperley have shown that boson world lines in 2+1 dimensions with periodic boundary conditions, regarded as ring polymers on a torus, have a mean square winding number given by $<W^2> = 2n_s\hbar^2/mk_BT$, where $m$ is the boson mass and $n_s$ is the superfluid number density. Here, we review the mapping of the statistical mechanics of polymers with constraints onto quantum mechanics, and show that there is an interesting generalization of the Pollock-Ceperley result to directed polymer melts interacting with a repulsive rod of radius $a$. When translated into boson language, the mean square winding number around the rod for a system of size $R$ perpendicular to the rod reads $<W^2> = {n_s\hbar^2\over 2πmk_BT}\ln(R/a)$. This result is directly applicable to vortices in Type II superconductors in the presence of columnar defects. An external current passing through the rod couples directly to the winding number in this case.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Irreversibility, Mechanical Entanglement and Thermal Melting in Superconducting Vortex Crystals with Point Impurities

We discuss the onset of irreversibility and entanglement of vortex lines in high Tc superconductors due to point disorder and thermal fluctuations using a simplified cage model. A combination of Flory arguments, known results from directed polymers in random media, and a Lindemann criterion are used to estimate the field and temperature dependence of irreversibility, mechanical entanglement and thermal melting. The qualitative features of this dependence, including its nonmonotonicity when disorder is sufficiently strong, are in good agreement with recent experiments.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Localization transitions in non-Hermitian quantum mechanics

We study the localization transitions which arise in both one and two dimensions when quantum mechanical particles described by a random Schrödinger equation are subjected to a constant imaginary vector potential. A path-integral formulation relates the transition to flux lines depinned from columnar defects by a transverse magnetic field in superconductors. The theory predicts that the transverse Meissner effect is accompanied by stretched exponential relaxation of the field into the bulk and a diverging penetration depth at the transition.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Quantum Smectic and Supersolid Order in Helium Films and Vortex Arrays

A flux liquid can condense into a smectic crystal in a pure layered superconductor with the magnetic field oriented nearly parallel to the layers. Similar order can arise in low temperature $^4$He films with a highly anisotropic periodic substrate potential. If the smectic order is commensurate with the layering, this periodic array is {\sl stable} to quenched random point--like disorder. By tilting and adjusting the magnitude of the applied field, both incommensurate and tilted smectic and crystalline phases are found for vortex arrays. Related variations are possible by changing the chemical potential in the helium system. We discuss transport near the second order smectic freezing transition, and show that permeation modes in superconductors lead to a small non--zero resistivity and a large but finite tilt modulus in the smectic crystal. In helium films, the theory predicts a nonzero superfluid density and propagating third sound modes, showing that the quantum smectic is always simultaneously crystalline and superfluid.

0 citations0 reviewsNo code linkedOpen access