We introduce a novel numerical method for direct simulation of front propagation in the Fisher-KPP equation with a time-dependent parameter on an infinite domain. The method computes a time-dependent boundary condition that accurately captures the leading-edge dynamics by coupling the nonlinear simulation region to a linear approximation region in which the dynamics can be solved exactly via the Green's function of the linearized equation. This approach enables precise front velocity measurements on relatively small computational domains for a variety of nonautonomous regimes and initial conditions for which existing numerical methods break down. We apply the method to pulled and pushed fronts in the Fisher-KPP equation with quadratic and quadratic-cubic nonlinearities, finding that it improves the accuracy of the simulated front velocity even for constant parameters and a fixed domain size. For pulled fronts with a diffusion coefficient that increases algebraically in time, our results reveal a deviation from the natural asymptotic velocity predicted by linear theory, whose explanation requires nonlinear theory. For pushed fronts with constant parameters, the method reproduces the
Front propagation into unstable states is often determined by the linearization, that is, propagation speeds agree with predictions from the linearized equation at the unstable state. The leading edge behavior is then a Gaussian tail propagating with the linear spreading speed. Fronts following this leading edge are commonly referred to as pulled fronts, alluding to the idea that they are ``pulled'' by this leading-edge Gaussian tail. We describe here a class of examples that exhibits how these leading-order effects do not completely describe the dynamics in the wake of the front. In fact, leading edge behavior predicts at most two possible invasion scenarios, associated with positive and negative amplitudes of the Gaussian tail, but our examples exhibit three or more invasion fronts with different states in the wake. The resulting invasion process therefore leaves behind a state that is not solely determined by the leading edge, and thus not just pulled by the Gaussian tail.
We investigate polymers pulled away from an interacting surface, where the force is applied to the untethered endpoint and at an angle $θ$ to the surface. We use the canonical self-avoiding walk model of polymers and obtain the phase diagram of the model using Monte Carlo simulations for a range of angles, temperatures and force magnitudes. The phase diagram of the model displays re-entrance at low temperatures for three-dimensional walks when the pulling is more vertical than horizontal. Our results agree with various exactly solvable lattice models that have been previously studied.
It is vital to infer a signed distance function (SDF) in multi-view based surface reconstruction. 3D Gaussian splatting (3DGS) provides a novel perspective for volume rendering, and shows advantages in rendering efficiency and quality. Although 3DGS provides a promising neural rendering option, it is still hard to infer SDFs for surface reconstruction with 3DGS due to the discreteness, the sparseness, and the off-surface drift of 3D Gaussians. To resolve these issues, we propose a method that seamlessly merge 3DGS with the learning of neural SDFs. Our key idea is to more effectively constrain the SDF inference with the multi-view consistency. To this end, we dynamically align 3D Gaussians on the zero-level set of the neural SDF using neural pulling, and then render the aligned 3D Gaussians through the differentiable rasterization. Meanwhile, we update the neural SDF by pulling neighboring space to the pulled 3D Gaussians, which progressively refine the signed distance field near the surface. With both differentiable pulling and splatting, we jointly optimize 3D Gaussians and the neural SDF with both RGB and geometry constraints, which recovers more accurate, smooth, and complete su
Fronts that start from a local perturbation and propagate into a linearly unstable state come in two classes: pulled and pushed. ``Pulled'' fronts are ``pulled along'' by the spreading of linear perturbations about the unstable state, so their asymptotic speed $v^*$ equals the spreading speed of linear perturbations of the unstable state. The central result of this paper is that the velocity of pulled fronts converges universally for time $t\to\infty$ like $v(t)=v^*-3/(2λ^*t) + (3\sqrtπ/2) Dλ^*/(D{λ^*}^2t)^{3/2}+O(1/t^2)$. The parameters $v^*$, $λ^*$, and $D$ are determined through a saddle point analysis from the equation of motion linearized about the unstable invaded state. The interior of the front is essentially slaved to the leading edge, and we derive a simple, explicit and universal expression for its relaxation towards $φ(x,t)=Φ^*(x-v^*t)$. Our result, which can be viewed as a general center manifold result for pulled front propagation, is derived in detail for the well known nonlinear F-KPP diffusion equation, and extended to much more general (sets of) equations (p.d.e.'s, difference equations, integro-differential equations etc.). Our universal result for pulled fronts
The Bramson logarithmic shift of the position of pulled fronts is a universal feature common to a large class of monostable traveling wave equations. As one varies the non-linearities it so happens that one can observe, at some critical non linearity, a transition from pulled fronts to pushed fronts. At this transition the Bramson shift is modified. In the limit where time goes to infinity and the non-linearity becomes critical, the position of the front exhibits a cross-over. The goal of the present note is to give the expression of this cross-over function, for a particular model which is exactly soluble, with the hope that this expression would remain valid for more general traveling wave equations at the transition between pulled and pushed fronts. Other cross-over functions are also obtained, for this particular model, to describe the dependence on initial conditions or the effect of a cut-off.
We analyze the transition between pulled and pushed fronts both analytically and numerically from a model-independent perspective. Based on minimal conceptual assumptions, we show that pushed fronts bifurcate from a branch of pulled fronts with an effective speed correction that scales quadratically in the bifurcation parameter. Strikingly, we find that in this general context without assumptions on comparison principles, the pulled front loses stability and gives way to a pushed front when monotonicity in the leading edge is lost. Our methods rely on far-field core decompositions that identify explicitly asymptotics in the leading edge of the front. We show how the theoretical construction can be directly implemented to yield effective algorithms that determine spreading speeds and bifurcation points with exponentially small error in the domain size. Example applications considered here include an extended Fisher-KPP equation, a Fisher-Burgers equation, negative taxis in combination with logistic population growth, an autocatalytic reaction, and a Lotka-Volterra model.
We investigate the inside structure of one-dimensional reaction-diffusion traveling fronts. The reaction terms are of the monostable, bistable or ignition types. Assuming that the fronts are made of several components with identical diffusion and growth rates, we analyze the spreading properties of each component. In the monostable case, the fronts are classified as pulled or pushed ones, depending on the propagation speed. We prove that any localized component of a pulled front converges locally to 0 at large times in the moving frame of the front, while any component of a pushed front converges to a well determined positive proportion of the front in the moving frame. These results give a new and more complete interpretation of the pulled/pushed terminology which extends the previous definitions to the case of general transition waves. In particular, in the bistable and ignition cases, the fronts are proved to be pushed as they share the same inside structure as the pushed monostable critical fronts. Uniform convergence results and precise estimates of the left and right spreading speeds of the components of pulled and pushed fronts are also established.
Understanding the temporal spread of gene drive alleles -- alleles that bias their own transmission -- through modeling is essential before any field experiments. In this paper, we present a deterministic reaction-diffusion model describing the interplay between demographic and allelic dynamics, in a one-dimensional spatial context. We focused on the traveling wave solutions, and more specifically, on the speed of gene drive invasion (if successful). We considered various timings of gene conversion (in the zygote or in the germline) and different probabilities of gene conversion (instead of assuming 100$\%$ conversion as done in a previous work). We compared the types of propagation when the intrinsic growth rate of the population takes extreme values, either very large or very low. When it is infinitely large, the wave can be either successful or not, and, if successful, it can be either pulled or pushed, in agreement with previous studies (extended here to the case of partial conversion). In contrast, it cannot be pushed when the intrinsic growth rate is vanishing. In this case, analytical results are obtained through an insightful connection with an epidemiological SI model. We
We analyze spatial spreading in a population model with logistic growth and chemorepulsion. In a parameter range of short-range chemo-diffusion, we use geometric singular perturbation theory and functional-analytic farfield-core decompositions to identify spreading speeds with marginally stable front profiles. In particular, we identify a sharp boundary between between linearly determined, pulled propagation, and nonlinearly determined, pushed propagation, induced by the chemorepulsion. The results are motivated by recent work on singular limits in this regime using PDE methods.
The concept of pulled fronts with a cutoff $ε$ has been introduced to model the effects of discrete nature of the constituent particles on the asymptotic front speed in models with continuum variables (Pulled fronts are the fronts which propagate into an unstable state, and have an asymptotic front speed equal to the linear spreading speed $v^*$ of small linear perturbations around the unstable state). In this paper, we demonstrate that the introduction of a cutoff actually makes such pulled fronts weakly pushed. For the nonlinear diffusion equation with a cutoff, we show that the longest relaxation times $τ_m$ that govern the convergence to the asymptotic front speed and profile, are given by $τ_m^{-1} \simeq [(m+1)^2-1] π^2 / \ln^2 ε$, for $m=1,2,...$.
The derivation of a Moving Boundary Approximation or of the response of a coherent structure like a front, vortex or pulse to external forces and noise, is generally valid under two conditions: the existence of a separation of time scales of the dynamics on the inner and outer scale and the existence and convergence of solvability type integrals. We point out that these conditions are not satisfied for pulled fronts propagating into an unstable state: their relaxation on the inner scale is power law like and in conjunction with this, solvability integrals diverge. The physical origin of this is traced to the fact that the important dynamics of pulled fronts occurs in the leading edge of the front rather than in the nonlinear internal front region itself. As recent work on the relaxation and stochastic behavior of pulled fronts suggests, when such fronts are coupled to other fields or to noise, the dynamical behavior is often qualitatively different from the standard case in which fronts between two (meta)stable states or pushed fronts propagating into an unstable state are considered.
We argue that while fluctuating fronts propagating into an unstable state should be in the standard KPZ universality class when they are {\em pushed}, they should not when they are {\em pulled}: The universal $1/t$ velocity relaxation of deterministic pulled fronts makes it unlikely that the KPZ equation is the appropriate effective long-wavelength low-frequency theory in this regime. Simulations in 2$D$ confirm the proposed scenario, and yield exponents $β\approx 0.29\pm 0.01$, $ζ\approx 0.40\pm 0.02$ for fluctuating pulled fronts, instead of the KPZ values $β=1/3$, $ζ= 1/2$. Our value of $β$ is consistent with an earlier result of Riordan {\em et al.}
Traveling waves describe diverse natural phenomena from crystal growth in physics to range expansions in biology. Two classes of waves exist with very different properties: pulled and pushed. Pulled waves are driven by high growth rates at the expansion edge, where the number of organisms is small and fluctuations are large. In contrast, fluctuations are suppressed in pushed waves because the region of maximal growth is shifted towards the population bulk. Although it is commonly believed that expansions are either pulled or pushed, we found an intermediate class of waves with bulk-driven growth, but exceedingly large fluctuations. These waves are unusual because their properties are controlled by both the leading edge and the bulk of the front.
We establish selection of critical pulled fronts in invasion processes. Our result shows convergence to a pulled front with a logarithmic shift for open sets of steep initial data, including one-sided compactly supported initial conditions. We rely on robust, conceptual assumptions, namely existence and marginal spectral stability of a front traveling at the linear spreading speed. We demonstrate that the assumptions hold for open classes of spatially extended systems. Previous results relied on comparison principles or probabilistic tools with implied non-open conditions on initial data and structure of the equation. Technically, we describe the invasion process through the interaction of a Gaussian leading edge with the pulled front in the wake. Key ingredients are sharp linear decay estimates to control errors in the nonlinear matching and corrections from initial data.
We analyze the dynamics of pattern forming fronts which propagate into an unstable state, and whose dynamics is of the pulled type, so that their asymptotic speed is equal to the linear spreading speed v^*. We discuss a method that allows to derive bounds on the front velocity, and which hence can be used to prove for, among others, the Swift-Hohenberg equation, the Extended Fisher-Kolmogorov equation and the cubic Complex Ginzburg-Landau equation, that the dynamically relevant fronts are of the pulled type. In addition, we generalize the derivation of the universal power law convergence of the dynamics of uniformly translating pulled fronts to both coherent and incoherent pattern forming fronts. The analysis is based on a matching analysis of the dynamics in the leading edge of the front, to the behavior imposed by the nonlinear region behind it. Numerical simulations of fronts in the Swift-Hohenberg equation are in full accord with our analytical predictions.
It has recently been proposed that fluctuating ``pulled'' fronts propagating into an unstable state should not be in the standard KPZ universality class for rough interface growth. We introduce an effective field equation for this class of problems, and show on the basis of it that noisy pulled fronts in {\em d+1} bulk dimensions should be in the universality class of the {\em (d+1)+1}D KPZ equation rather than of the {\em d+1}D KPZ equation. Our scenario ties together a number of heretofore unexplained observations in the literature, and is supported by previous numerical results.
We study the propagation of a ``pulled'' front with multiplicative noise that is created by a local perturbation of an unstable state. Unlike a front propagating into a metastable state, where a separation of time scales for sufficiently large $t$ creates a diffusive wandering of the front position about its mean, we predict that for so-called pulled fronts, the fluctuations are subdiffusive with root mean square wandering $Δ(t) \sim t^{1/4}$, {\em not} $t^{1/2}$. The subdiffusive behavior is confirmed by numerical simulations: For $t \le 600$, these yield an effective exponent slightly larger than 1/4.
Streamer ionization fronts are pulled fronts propagating into a linearly unstable state; the spatial decay of the initial condition of a planar front selects dynamically one specific long time attractor out of a continuous family. A transverse stability analysis has to take these features into account. In this paper we introduce a framework for this transverse stability analysis, involving stable and unstable manifolds in a weighted space. Within this framework, a numerical dynamical systems method for the calculation of the dispersion relation as an eigenvalue problem is defined and dispersion curves for different values of the electron diffusion constant and of the electric field ahead of the front are derived. Numerical solutions of the initial value problem confirm the eigenvalue calculations. The numerical work is complemented with analytical expressions for the dispersion relation in the limit of small and large wave numbers and with a fit formula for intermediate wave numbers. This empirical fit supports the conjecture that the smallest unstable wave length of the Laplacian instability is proportional to the diffusion length that characterizes the leading edge of the pulled
We consider ponderomotive forces acting on small particles in propagating wave packets (pulses). Specifically, we analyze simple point particles as well as composite dipole and dumbbell particles in the fields of forward-propagating (parallel phase and group velocities) and backward-propagating (antiparallel phase and group velocities) wave packets. Depending on the characteristics of the wave packet, particles may be pushed away from the wave source or pulled toward it. We also examine particle dynamics in the field of a beat wave generated by two forward-propagating waves with slightly different frequencies. Such a beat wave can emulate a periodic sequence of either forward- or backward-propagating pulses. In particular, this provides a simple mechanism for realizing pulling forces as employed in optical and acoustic `tractor beams'.