Freidlin-Gärtner formula and asymptotic profile in reaction-diffusion equations

Abstract: We address the question of the large-time behavior of solutions to reaction-diffusion equations in periodic media. We start with the description of the asymptotic shape of the invasion set, which is characterized by the Freidlin-Gärtner formula. We outline a proof of the formula that holds true for general types of reaction terms. We then present some recent results, obtained in collaboration with H. Guo and F. Hamel, for (weakly) bistable equations. They include a regular version of the Freidlin-Gärtner formula and the convergence in profile towards pulsating traveling fronts for solutions with either bounded or unbounded initial support.

• The paper establishes that the Freidlin-Gärtner formula determines the asymptotic spreading set and invasion speeds in periodic reaction-diffusion equations.
• It employs rigorous upper and lower bound methods using pulsating traveling fronts to capture direction-dependent speeds and local convergence properties.
• The analysis shows how the global spreading geometry and asymptotic profiles are linked to the medium’s periodic structure and specific reaction nonlinearity.

Freidlin-Gärtner Formula and Asymptotic Profile in Reaction-Diffusion Equations

Introduction

This paper investigates the long-time dynamics of reaction-diffusion equations posed in periodic heterogeneous media, focusing on the geometry and qualitative behavior of spatial spreading from localized initial conditions. The central mathematical object is the asymptotic set of spreading, whose shape, propagation speeds, and local profile are dictated by the structure of the medium and the nonlinear reaction terms. The analysis synthesizes classical results in the homogeneous case, extends them to anisotropic periodic environments, and clarifies the universality of the Freidlin-Gärtner formula for a broad class of reaction terms, not limited to the monostable regime.

Background and Motivation

In homogeneous media, classical results establish that the speed and shape of invasion for localized initial data are determined by the properties of planar traveling fronts. Specifically, invasion is isotropic and characterized by a single critical propagation speed $c^*$. For periodic heterogeneous environments, spatial anisotropy and periodic coefficients lead to direction-dependent spreading speeds and nontrivial, typically non-spherical, shapes for invaded regions.

The transition from the homogeneous to the periodic context introduces new technical challenges: the existence, speed, and profile of traveling fronts become direction-dependent, and the interplay between local geometry and global dynamics becomes richer. These difficulties have motivated the investigation of the generalized Freidlin-Gärtner formula and its implications for both the global extent and the local behavior of solutions.

The Freidlin-Gärtner Formula in Periodic Media

A principal contribution of the paper is a rigorous and unified proof of the Freidlin-Gärtner formula for the asymptotic spreading speed in general reaction-diffusion equations with periodic coefficients, covering monostable, ignition, and (weakly) bistable nonlinearities. For the scalar reaction-diffusion equation

$$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$

with periodic $A(x)$, $q(x)$, and $f(x,u)$, the large-time invaded region associated with compactly supported initial data is described by the asymptotic set of spreading $W$. The set $W$ is determined through the family of direction-dependent critical speeds $c^*(e)$ of pulsating traveling fronts, and the formula reads

$$W = \bigcap_{e \in S^{N-1}} \{ x \in \mathbb{R}^N : x \cdot e < c^*(e) \},$$

where $e$ ranges over unit directions and $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$0 is the minimal (or unique) speed of pulsating traveling fronts in direction $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$1. This construction yields a Wulff shape $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$2, which is convex and generally anisotropic.

Importantly, the formula applies not only to monostable cases (where propagation is universal for positive initial data), but also to ignition and bistable nonlinearities under suitable hypotheses guaranteeing the existence of pulsating fronts in all directions.

Outline of the Proof Strategy

The crux of the generalized proof is to establish simultaneous control on the spreading in all directions, leveraging the geometric properties of expanding sets and the normal velocity determined by the local value $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$3. The argument involves:

• Upper Bound: Utilizing exponential decay properties of both the solution and the pulsating front profiles, one shows that the solution remains negligible outside the expanding set $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$4, with the necessary control given by sliding subsolutions corresponding to the pulsating fronts.
• Lower Bound: For points inside $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$5, a contradiction argument based on first contact points between the solution and suitable sublevel sets (translated plateaus) is deployed. Geometric properties of $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$6 and comparison with whole-space entire solutions (fronts or their compactly supported modifications) ensure that the minimal possible velocity is realized along any admissible direction.

This approach simplifies and unifies previous treatments, capturing invasion dynamics for general periodic reaction-diffusion systems.

Asymptotic Profile: Convergence to Pulsating Fronts

Beyond the global shape of spreading, the paper addresses the local-in-space limiting behavior (asymptotic profile) of solutions. Two primary cases are considered:

• Compactly Supported Initial Data: The $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$7-limit set of translated solutions (in both space and time) contains, in every regular direction, the profiles of pulsating traveling fronts corresponding to the critical velocity of that direction. When the boundary $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$8 is differentiable, the family of possible asymptotic profiles comprises all spatial and temporal translates of the fronts, in addition to the constant states $$\partial_t u - \nabla \cdot (A(x) \nabla u) - q(x) \cdot \nabla u = f(x,u), \qquad x \in {\mathbb{R}}^N, \; t>0,$$9 and $A(x)$0. In corners of $A(x)$1, multiple front profiles may simultaneously arise in the limit.
• General Unbounded Initial Data: For initial data with extended or unbounded support, the invaded set is a Minkowski sum of the initial support and the Wulff shape associated with homogeneous invasion (cf. $A(x)$2 for $A(x)$3). At regular points of $A(x)$4, the asymptotic profile along lines approaching the boundary is again given by the associated pulsating front, validating a local form of the Freidlin-Gärtner formula.

The sliding method, refined to treat asymptotic convergence to the front profiles, exploits regularity and comparison principles to assert that limits of properly translated solutions must match translates of the pulsating front in the supporting direction.

Theoretical and Practical Implications

These results underline the centrality of pulsating fronts: they are not merely special solutions but fundamentally determine both the macroscopic spreading set and the local profile of general solutions at large times. The geometry of $A(x)$5—explicitly linked to the family of front speeds—encodes all information about the global invasion process, while the set of asymptotic (local) profiles is intimately tied to the structure of the medium and the class of nonlinearity.

The presence of corners or singularities in $A(x)$6 implies the possible appearance of multiple front profiles in these directions, and the criteria for global convergence to a single front remain subtle, especially in higher dimensions or for non-compact initial data.

The methodology extends to a variety of nonlinearities and scenarios, with implications for both deterministic propagation in biological, chemical, and physical systems and for stochastic models via the connection to large deviation principles and eigenvalue problems.

Conclusion

The Freidlin-Gärtner formula provides a comprehensive description of both global and local invasion dynamics in reaction-diffusion equations set in periodic media, applicable to monostable, ignition, and bistable nonlinearities. The propagation set and local profile are determined by the speeds and shapes of associated pulsating traveling fronts, encapsulated in an explicit geometric formula. Future research directions include full characterization of convergence (beyond subsequences), extension to non-periodic heterogeneities, and further exploration of the interplay between geometry, nonlinearities, and higher-order effects (2604.27556).