We develop a discrete differential geometry for surfaces of non-constant negative curvature, which can be used to model various phenomena from the growth of flower petals to marine invertebrate swimming. Specifically, we derive and numerically integrate a version of the classical Lelieuvre formulas that apply to immersions of C1,1 hyperbolic surfaces of non-constant curvature. In contrast to the constant curvature case, these formulas do not provide an explicit method for constructing an immersion but rather describe an immersion via an implicit set of equations. We propose an iterative method for resolving these equations. Because we are interested in scenarios where the curvature is a function of the intrinsic material coordinates, in particular, on the geodesic distance from an origin or from an edge, we suggest a fast marching method for computing geodesic distance on manifolds. We apply our methods to generate surfaces of non-constant curvature and demonstrate how one can introduce branch points to account for the multi-generational buckling and subwrinkling observed in many applications.

We develop a discrete differential geometry for surfaces of non-constant negative curvature, which can be used to model various phenomena from the growth of flower petals to marine invertebrate swimming. Specifically, we derive and numerically integrate a version of the classical Lelieuvre formulas that apply to immersions of $C^{1,1}$ hyperbolic surfaces of non-constant curvature. In contrast to the constant curvature case, these formulas do not provide an explicit method for constructing an immersion but rather describe an immersion via an implicit set of equations. We propose an iterative method for resolving these equations. Because we are interested in scenarios where the curvature is a function of the intrinsic material coordinates, in particular, on the geodesic distance from an origin, we suggest a fast marching method for computing geodesic distance on manifolds. We apply our methods to generate surfaces of non-constant curvature and demonstrate how one can introduce branch points to account for the sort of multi-generational buckling and subwrinkling observed in many applications.

There are three aims of this paper. The first is to explain the reasons for behavior we had long suspected to be true but the real reasons for which we could never quite nail down. Modulated striped patterns arising from a wide class of gradient microscopic pattern forming systems display universal behavior. Their order parameter evolves according to a phase diffusion equation that derives from an averaged energy that consists of coordinate invariant combinations of the coefficients in the metric and curvature two forms of a well-defined phase surface. The evolution towards universality is asymptotic in that the pattern evolves in such a way that, over longer and longer time scales, many terms from the microscopic energy simply become negligible leaving behind canonical forms common to a wide class of microscopic pattern-forming systems. The second aim is to emphasize with some new results the key role that the Jacobian matrix of the map from physical to order parameter space plays in both two and three dimensions. In two dimensions, it is closely related to the Gaussian curvature of the phase surface. It is a conserved density whose integral over space in two dimensions or on cross-sections in three become boundary terms that measure the topological indices of point defects in two dimensions and loop defects in three. In all dimensions, the Jacobian matrix acting on the order parameter vector, the gradient of the phase, is zero when the local pattern wavenumber is close to its preferred value and this leads to the effective linearization of the phase diffusion equation. The third aim is to honor Hermann Flaschka, a close friend and scientific colleague for over fifty years, an outstanding mathematician, a true gentleman and scholar with an uncanny knack of explaining the most complicated of ideas in the simplest of ways, who passed away last year. Hermann was one of the founding editors of Physica D and served as the coordinating editor for almost twenty years.

We investigate stripe patterns formation far from threshold using a combination of topological, analytic, and numerical methods. We first give a definition of the mathematical structure of `multi-valued' phase functions that are needed for describing layered structures or stripe patterns containing defects. This definition yields insight into the appropriate `gauge symmetries' of patterns, and leads to the formulation of variational problems, in the class of special functions with bounded variation, to model patterns with defects. We then discuss approaches to discretize and numerically solve these variational problems. These energy minimizing solutions support defects having the same character as seen in experiments.

We consider the scaling behavior of the range and $p$-multiple range, that is the number of points visited and the number of points visited exactly $p\geq 1$ times, of simple random walk on ${\mathbb Z}^d$, for dimensions $d\geq 2$, up to time of exit from a domain $D_N$ of the form $D_N = ND$ where $D\subset {\mathbb R}^d$, as $N\uparrow\infty$. Recent papers have discussed connections of the range and related statistics with the Gaussian free field, identifying in particular that the distributional scaling limit for the range, in the case $D$ is a cube in $d\geq 3$, is proportional to the exit time of Brownian motion. The purpose of this note is to give a concise, different argument that the scaled range and multiple range, in a general setting in $d\geq 2$, both weakly converge to proportional exit times of Brownian motion from $D$, and that the corresponding limit moments are `polyharmonic', solving a hierarchy of Poisson equations.

We correct the statement of lemma 3.5 in [Jaramillo et al (2018 Nonlinearity 31 4162)]. The lemma establishes Fredholm properties for a specific radial operator, which in the original paper is incorrectly defined. The error does not alter our main theorem, but does affect the proof of a subsequent result, lemma 4.11, which we fix in this corrigendum.

We develop a theory for distributed branch points and investigate their role in determining the shape and influencing the mechanics of thin hyperbolic objects. We show that branch points are the natural topological defects in hyperbolic sheets, they carry a topological index which gives them a degree of robustness, and they can influence the overall morphology of a hyperbolic surface without concentrating energy. We develop a discrete differential geometric approach to study the deformations of hyperbolic objects with distributed branch points. We present evidence that the maximum curvature of surfaces with geodesic radius R containing branch points grow sub-exponentially, $$O(e^{c\sqrt{R}})$$ in contrast to the exponential growth $$O(e^{c' R})$$ for surfaces without branch points. We argue that, to optimize norms of the curvature, i.e., the bending energy, distributed branch points are energetically preferred in sufficiently large pseudospherical surfaces. Further, they are distributed so that they lead to fractal-like recursive buckling patterns.

The goal of this paper is to describe the oscillatory microstructure that can emerge from minimizing sequences for nonconvex energies. We consider integral functionals that are defined on real valued (scalar) functions $u(x)$ which are nonconvex in the gradient $\nabla u$ and possibly also in $u$. To characterize the microstructures for these nonconvex energies, we minimize the associated relaxed energy using two novel approaches: i) a semi-analytical method based on control systems theory, ii) and a numerical scheme that combines convex splitting together with a modified version of the split Bregman algorithm. These solutions are then used to gain information about minimizing sequences of the original problem and the spatial distribution of microstructure.

Defects are a ubiquitous feature of ordered media. They have certain universal features, independent of the underlying physical system, reflecting their topological origins. While the topological properties of defects are robust, they appear as ‘unphysical’ singularities, with non-integrable energy densities in coarse-grained macroscopic models. We develop a principled approach for enriching coarse-grained theories with enough of the ‘micro-physics’ to obtain thermodynamically consistent, well-set models that allow for the investigations of dynamics and interactions of defects in extended systems. We also develop associated numerical methods that are applicable to computing energy driven behaviors of defects across the amorphous-soft-crystalline materials spectrum. Our methods can handle order parameters that have a head-tail symmetry, i.e. director fields, in systems with a continuous translation symmetry, as in nematic liquid crystals, and in systems where the translation symmetry is broken, as in smectics and convection patterns. We illustrate our methods with explicit computations.

A ubiquitous motif in nature is the self-similar hierarchical buckling of a thin lamina near its margins. This is seen in leaves, flowers, fungi, corals, and marine invertebrates. We investigate this morphology from the perspective of non-Euclidean plate theory. We identify a novel type of defect, a branch-point of the normal map, that allows for the generation of such complex wrinkling patterns in thin elastic hyperbolic surfaces, even in the absence of stretching. We argue that branch points are the natural defects in hyperbolic sheets, they carry a topological charge which gives them a degree of robustness, and they can influence the overall morphology of a hyperbolic surface without concentrating elastic energy. We develop a theory for branch points and investigate their role in determining the mechanical response of hyperbolic sheets to weak external forces.

We develop reduced, stochastic models for high dimensional, dissipative dynamical systems that relax very slowly to equilibrium and can encode long term memory. We present a variety of empirical and first principles approaches for model reduction, and build a mathematical framework for analyzing the reduced models. We introduce the notions of universal and asymptotic filters to characterize ‘optimal’ model reductions for sloppy linear models. We illustrate our methods by applying them to the practically important problem of modeling evaporation in oil spills.

Galaxies are built by complex physical processes with significant inherent stochasticity. It is therefore surprising that the inferred dark matter distributions in galaxies are correlated with the observed baryon distributions leading to various ‘Baryon-Halo conspiracies’. The fact that no dark matter candidate has been definitively identified invites a search for alternative explanations for such correlations and we present an approach motivated by the behaviors of self organized patterns. We propose a nonlocal relativistic Lagrangian theory for a ‘pattern field’ which acts as an ‘effective dark matter’, built on the idea that defects in this pattern field couple to the baryonic matter distribution.The model applies to rotation supported systems and, for them, we compute galactic rotation curves, obtain a radial acceleration relation with two branches, and deduce the Freeman limit for central surface brightness.

We argue that a natural explanation for a variety of robust galaxy scaling relations comes from the perspective of pattern formation and self-organization as a result of symmetry breaking. We propose a simple Lagrangian model that combines a conventional model for normal matter in a galaxy with a conventional model for stripe pattern formation in systems that break continuous translation invariance. We show that the energy stored in the pattern field acts as an effective dark matter. Our theory reproduces the gross features of elliptic galaxies as well as disk galaxies (HSB and LSB) including their detailed rotation curves, the radial acceleration relation, and the Freeman law. We investigate the stability of disk galaxies in the context of our model and obtain scaling relations for the central dispersion for elliptical galaxies. A natural interpretation of our results is that (1) ‘dark matter’ is potentially a collective, emergent phenomenon and not necessarily an as yet undiscovered particle, and (2) MOND is an effective theory for the description of a self-organized complex system rather than a fundamental description of nature that modifies Newton’s second law.

Growth and other dynamical processes in soft materials can create novel types of mesoscopic defects including discontinuities for the second and higher derivatives of the deformation, and terminating defects for these discontinuities. These higher-order defects move “easily", and can thus confer a great degree of flexibility to the material. We develop a general continuum mechanical framework from which we can derive the dynamics of higher order defects in a thermodynamically consistent manner. We illustrate our framework by obtaining the explicit dynamical equations for the next higher order defects in an elastic body beyond dislocations, phase boundaries, and disclinations, namely, surfaces of inflection and branch lines.

In this essay we explore analogies between macroscopic patterns, which result from a sequence of phase transitions/instabilities starting from a homogeneous state, and similar phenomena in cosmology, where a sequence of phase transitions in the early universe is believed to have separated the fundamental forces from each other, and also shaped the structure and distribution of matter in the universe. We discuss three distinct aspects of this analogy: (i) Defects and topological charges in macroscopic patterns are analogous to spins and charges of quarks and leptons; (ii) Generic (3+1) stripe patterns carry an (energy) density that accounts for phenomena that are currently attributed to dark matter; (iii) Space-time patterns of interacting nonlinear waves display behaviors reminiscent of quantum phenomena including inflation, entanglement and dark energy.

This is an invited commentary on "Geometrically incompatible confinement of solids", B. Davidovitch, Y. Sun and G. M. Grason (PNAS, doi:10.1073/pnas.1815507116 , arxiv:1809.06919).

Growth processes in many living organisms create thin, soft materials with an intrinsically hyperbolic geometry. These objects support novel types of mesoscopic defects - discontinuity lines for the second derivative and branch points - terminating defects for these line discontinuities. These higher-order defects move easily, and thus confer a great degree of flexibility to thin hyperbolic elastic sheets. We develop a general, higher-order, continuum mechanical framework from which we can derive the dynamics of higher order defects in a thermodynamically consistent manner. We illustrate our framework by obtaining the explicit equations for the dynamics of branch points in an elastic body.

We consider an evolving preferential attachment random graph model where at discrete times a new node is attached to an old node, selected with probability proportional to a superlinear function of its degree. For such schemes, it is known that the graph evolution condenses, that is a.s. in the limit graph there will be a single random node with infinite degree, while all others have finite degree. In this note, we establish a.s. law of large numbers type limits and fluctuation results, as $n\uparrow\infty$, for the counts of the number of nodes with degree $k\geq 1$ at time $n\geq 1$. These limits rigorously verify and extend a physical picture of Krapivisky, Redner and Leyvraz (2000) on how the condensation arises with respect to the degree distribution.