REVIEW 6 minor 1 cited by
Asymptotic spreading of interacting species with multiple fronts I: A geometric optics approach
T0 review · 0 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper proves that a slower competing species invading empty habitat can spread at a nonlocally pulled speed, given explicitly in terms of the model parameters, and that the solution develops exactly three invasion fronts with these…
desk verdict Exact spreading speeds for weak competition with a nonlocally pulled front; the main theorem is new and the Hamilton-Jacobi framework works, but the verification of Freidlin's condition needs a fill-in. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the WKB transformation w_ε(t,x) = -ε log u(t/ε, x/ε) and its half-relaxed limits w_* and w^*, which satisfy viscosity sub- and super-solutions of a Hamilton-Jacobi equation whose coefficient switches at the known faster front x=2√dr t. The limits are pinned down by comparison with the explicit solution w1 = max{J1,0}, where J1 is the infimum of the path integral ∫ [|γ̇|²/4 - 1 + a χ_{γ≤c1 s}] ds over curves ending at (t,x). Solving this variational problem gives the front x=c_nlp t, and an exponential decay estimate for u along a line x=ĉt then controls the spreading speed from above via comparison with traveling waves. The machinery is the Evans-Souganidis PDE approach to geometric optics, including half-relaxed limits and Freidlin's condition, which handles coefficients depending on several moving frames.
What would settle it
Run a high-resolution numerical simulation of system (1) with a=0.6, b=0.5, d=1.5, r=1 and initial data u0=χ_{[-1000,0]}, v0=χ_{[-20,0]}, and measure sup{x:u(t,x)>0.4}/t at t=200 and t=400; Theorem 1.3 predicts the ratios converge to c2 = max{c_LLW, c_nlp} ≈ 1.3387, so a limit below that, such as 1.265, would falsify the claimed formula.
Extended reading notes
Core claim
The central discovery is that the spreading speed c2 of the slower species u equals max{c_LLW, c_nlp}, where c_LLW is the classical spreading speed of the coexistence equilibrium into the v-only state, and c_nlp has the explicit formula c_nlp = √dr - √a + (1-a)/(√dr-√a) when √dr ≤ √a + √(1-a), and c_nlp = 2√(1-a) otherwise. The speed c_nlp is characterized by the zero level set of the viscosity solution w1 of the Hamilton-Jacobi equation min{∂t w + |∂x w|^2 + 1 - a χ_{x≤2√dr t}, w} = 0 with initial data 0 on the left half-line and ∞ on the right. In the variational representation of w1, minimizing paths spend time ahead of the moving front, so c_nlp depends on the environment in front of the invasion, making it nonlocally pulled. The paper proves convergence to the four homogeneous equilibria in the regions separated by the three fronts, and obtains the borderline case dr=1 as a limit.
Load-bearing premise
The proof relies on the slower species being exactly absent on a right half-line initially, so that the WKB variable starts as infinity there; if u0 only decays exponentially, the explicit formula for c_nlp no longer follows from the argument.
Editorial extensions
If this is right
- If Theorem 1.3 is correct, the slower species u can spread strictly faster than any traveling wave of the homogeneous system would permit, with c_nlp > c_LLW whenever the faster species is not too fast relative to the competition parameters.
- The asymptotic state of the system is fully classified by three speeds: ahead of c1 t the habitat is empty, between c2 t and c1 t only v persists at density 1, between c3 t and c2 t both species coexist at (k1,k2), and left of c3 t only u persists at density 1.
- When dr=1, the two rightward fronts merge into one front at speed 2, so the invasion goes (1,0) ← (k1,k2) → (0,0) with the coexistence region filling the interval between -c̃_LLW t and 2t.
- The same geometric optics framework is claimed to extend to systems with three or more species and to higher dimensions, because the limiting Hamilton-Jacobi problem accommodates coefficients depending on several distinct front speeds.
Reading between the lines
- If the compact-support hypothesis (H∞) is relaxed to merely exponential decay of u0 on the right, the WKB initial data would be finite rather than infinite, so the limiting Hamilton-Jacobi problem and the formula for c_nlp should change; a testable prediction is that the second front speed becomes a function of the decay rate and approaches c_nlp as the decay rate tends to infinity.
- The algebraic coincidence of c_nlp with the speed found in a related monostable competition case suggests that nonlocally pulled fronts are a general mechanism in reaction-diffusion systems where a fast front modifies the effective growth rate ahead of a slower front, not a peculiarity of this parameter regime.
- A direct numerical test could look for the sharp transition in (4): plotting the measured second-front speed against √dr should show a kink exactly at √dr = √a + √(1-a), with the speed formula switching branches there.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the invasion of two competing species modeled by the Lotka–Volterra competition-diffusion system (1) under the initial condition (H∞): species u initially occupies the left half-line, species v is compactly supported, and the right habitat is empty. The main result (Theorem 1.3) states that for dr>1 the solution develops three fronts with speeds c1=2√(dr), c2=max{c_LLW,c_nlp}, and c3=−\tilde c_LLW, where c_LLW and \tilde c_LLW are the single-front spreading speeds of Lewis et al. and c_nlp is given explicitly by (4); between the fronts the solution converges to (0,0), (0,1), (k1,k2), and (1,0), respectively. The proof combines comparison arguments for c1 and c3 and rough bounds on c2, a WKB transform w_ε=−ε log u_ε with half-relaxed limits leading to the Hamilton–Jacobi equations (33)–(34), a variational solution w1 of the limiting problem whose zero set is {x≤c_nlp t}, and a large-deviation estimate together with wave comparisons (Lemma 2.4) to obtain the upper bound on c2. Theorem 1.4 treats the case dr=1 by a perturbation argument using the dr>1 and dr<1 results.
Significance. The result is significant: it provides the first exact determination of the second spreading speed in the two-species competition invasion problem with separated fronts, resolving a question raised by Shigesada and Kawasaki, and it exhibits a speed c_nlp that is nonlocally determined and can strictly exceed the minimal traveling wave speed. The geometric-optics method, adapted from Freidlin and Evans–Souganidis, is a promising tool for multi-front problems, and the paper contains detailed appendices for the variational calculus and comparison arguments, together with numerical simulations matching the predicted speeds. The main formula is explicit and falsifiable, and the overall derivation is coherent. The only substantive weakness is that one load-bearing verification, Freidlin’s condition in Proposition 3.5(b), is too terse as written; however, the missing argument is local and straightforward, so I do not view it as a threat to the main theorem.
minor comments (6)
- [Proposition 3.5(b), Eq. (41)] The verification of Freidlin’s condition is too terse. In the case c1/2>√a+√(1−a), the boundary is ∂P={x=2√(1−a)t} and the unconstrained minimizer for J1=0 is the straight line lying on ∂P, so (41) is not immediate. Please add the explicit perturbation argument, for example by considering paths γ_ε(s)=2√(1−a)s+ε s(t−s)/t, whose action tends to 0 because the boundary lies strictly below x=c1s and L1 is unchanged there.
- [Lemma 3.4] The sentence “It remains to check that w∗(0,x)=∞ for x>0” appears to refer to the lower half-relaxed limit w_* rather than the upper one, since the proof that follows uses a limit inferior; please correct the symbol to avoid confusion.
- [Section 2 and Proposition 2.1] The notation for maximal and minimal spreading speeds loses the overline/underline distinction in several places, in particular in Proposition 2.1(ii) and in Step 7 of its proof; please disambiguate which of \underline c2 and \bar c2 is used in (11b) and (11c).
- [Lemma 3.2] The line “Furthermore,” immediately before the Proof is an incomplete sentence and should be removed or completed.
- [Proposition 2.1, Step 4] The construction of the compactly supported function \tilde u0 with 0≤\tilde u0≤u0 should be stated explicitly, for instance by taking a cut-off of u0 on a bounded interval inside (−∞,0], rather than merely asserted.
- [Lemma 2.4, Appendix A] The existence of a traveling wave for the perturbed system (71) at any speed above c^δ_LLW is invoked without proof or reference; please add a citation or a short justification.
Circularity Check
No significant circularity: the nonlocally pulled speed c_nlp is derived from an independent HJ variational calculation; the sole self-citation is non-load-bearing.
full rationale
The central new quantity c_nlp is not an input recycled as an output. In Section 3, the WKB transform (28)–(29) and the half-relaxed limits (32) produce the limiting Hamilton–Jacobi equation (7), whose coefficients use only the independently established speed c1=2√dr and the hypothesis (H∞). Proposition 3.5(a) then computes J1 explicitly by minimizing the Lagrangian (39), and Lemma 3.7 identifies w1=max{J1,0} via Freidlin's condition (41); the explicit zero set {(t,x): x≤c_nlp t} is therefore a consequence of the variational formula, not a fitted postulate. Proposition 4.1 and Proposition 4.2 convert this exponential estimate into the upper and lower bounds for c2, with the final value max{cLLW,c_nlp}. The only overlap with prior work by one of the authors is the remark that formula (4) coincides with [18, Theorem 1.1]; the paper does not use [18] to prove the formula or to justify any load-bearing step, so this self-citation does not constitute circularity. External inputs cLLW and c̃LLW are quoted from Lewis et al. and used as benchmarks, not as assumptions about c2. One internal caveat, noted by a careful reader, is that the verification of Freidlin's condition (41) in Proposition 3.5(b) is terse and may be incomplete at the boundary of P; however, a possible proof gap is a correctness risk, not a circularity of the derivation chain.
Assumptions & free parameters
assumptions (5)
- domain assumption Single-front spreading speeds for the competition system exist and satisfy the bounds of Lewis et al. (Theorem 1.1 and Remark 1.2).
- standard math Comparison principle for weakly coupled parabolic systems of competition type.
- standard math Viscosity solution theory for Hamilton-Jacobi equations, including half-relaxed limits, comparison, and the Evans-Souganidis representation formula.
- standard math Continuity of spreading speeds with respect to model parameters (Volpert et al., Ch.3 Theorem 4.2).
- standard math Freidlin's theorem that the viscosity solution of the Hamilton-Jacobi equation equals max(J,0) when J satisfies the Freidlin condition.
Cite this review
Pith. "Pith review of Asymptotic spreading of interacting species with multiple fronts I: A geometric optics approach." pith.science (2026). https://pith.science/paper/ADNI4QHJ
@misc{pith2026190805025,
author = {Pith},
title = {Pith review of: Asymptotic spreading of interacting species with multiple fronts I: A geometric optics approach},
year = {2026},
howpublished = {\url{https://pith.science/paper/ADNI4QHJ}},
note = {Machine review of arXiv:1908.05025}
}
read the original abstract
We establish spreading properties of the Lotka-Volterra competition-diffusion system. When the initial data vanish on a right half-line, we derive the exact spreading speeds and prove the convergence to homogeneous equilibrium states between successive invasion fronts. Our method is inspired by the geometric optics approach for Fisher-KPP equation due to Freidlin, Evans and Souganidis. Our main result settles an open question raised by Shigesada et al. in 1997, and shows that one of the species spreads to the right with a nonlocally pulled front.
Figures
Forward citations
Cited by 1 Pith paper
-
On the logarithmic correction of transition fronts in shifting environments
For Fisher-KPP equations with a piecewise-constant shifting environment, the paper derives the exact logarithmic correction to the front location, extending Bramson's correction to growing domains and moving habitat b...
Reference graph
Works this paper leans on
- [1]
-
[2]
A. Alhasanat and C. Ou, On a conjecture raised by Yuzo Hosono, J. Dynam. Differential Equations , 31 (2019), 287–304
work page 2019
-
[3]
A. Alhasanat and C. Ou, Minimal-speed selection of travelin g waves to the Lotka-Volterra competition model, J. Differential Equations , 26 (2019), 7357–7378
work page 2019
- [4]
-
[5]
G. Barles and B. Perthame, Discontinuous solutions of deter ministic optimal stopping time problems, Model. Math. Anal. Numer. , 21 (1987), 557–579
work page 1987
- [6]
-
[7]
H. Berestycki and J. Fang, Forced waves of the FisherKPP equa tion in a shifting environment, J. Differential Equations , 264 (2018), 2157–2183
work page 2018
-
[8]
H. Berestycki and G. Nadin, Spreading speeds for one-dimens ional monostable reaction- diffusion equations, J. Math. Phys. , 53 (2012), 115619
work page 2012
Show all 43 references
-
[9]
Carr` ere, Spreading speeds for a two-species competitio n-diffusion system, J
C. Carr` ere, Spreading speeds for a two-species competitio n-diffusion system, J. Differential Equations , 264 (2018), 2133–2156
2018
-
[10]
M. B. Davis, Quaternary history and stability of forest communities , In: D.C. W est, H. H. Shugart, D. B. Botkin (editors), Forest Succession: Concep ts and Applications, Springer- Verlag, New York, 1981, pp. 132–153
1981
-
[11]
Y. H. Du and C. H. W u, Spreading with two speeds and mass segreg ation in a diffusive com- petition system with free boundaries, Calc. Var. Partial Differential Equations , 57 (2018), 52
2018
-
[12]
Ducrot, T
A. Ducrot, T. Giletti and H. Matano, Spreading speeds for mul tidimensional reaction-diffusion systems of the prey-predator type, Calc. Var. Partial Differential Equations , 58 (2019), 137
2019
-
[13]
L. C. Evans and P. E. Souganidis, Differential games and repre sentation formulas for solutions of Hamilton-Jacobi-Isaacs equation, Indiana Univ. Math. J. , 33 (1984), 773–797
1984
-
[14]
L. C. Evans and P. E. Souganidis, A PDE approach to geometric o ptics for certain semilinear parabolic equations, Indiana Univ. Math. J. , 38 (1989), 141–172
1989
-
[15]
J. Fang, Y. J. Lou and J. H. W u, Can pathogen spread keep pace wi th its host invasion? SIAM J. Appl. Math. , 74 (2016), 1633–1657
2016
-
[16]
R. A. Fisher, The wave of advance of advantageous genes, Ann. Hum. Genet. , 7 (1937), 355–369
1937
-
[17]
M. I. Freidlin, Limit theorems for large deviations and reac tion-diffusion equation, Ann. Probab., 13 (1985), 639–675
1985
-
[18]
Girardin and K
L. Girardin and K. Y. Lam, Invasion of an empty habitat by two competitors: spreading properties of monostable two-species competition-diffusi on systems, P. Lond. Math. Soc. , 119 (2019), 1279–1335
2019
-
[19]
J. S. Guo and C. H. W u, Dynamics for a two-species competition diffusion model with two free boundaries, Nonlinearity, 28 (2015), 1–27
2015
-
[20]
Hamel and G
F. Hamel and G. Nadin, Spreading properties and complex dyna mics for monostable reaction- diffusion equations, Comm. Partial Differential Equations , 37 (2012), 511–537
2012
-
[21]
Holzer and A
M. Holzer and A. Scheel, Accelerated fronts in a two stage inv asion process, SIAM J. Math. Anal. , 46 (2014), 397–427
2014
-
[22]
X. J. Hou and A. W. Leung, Traveling wave solutions for a compe titive reactiondiffusion system and their asymptotics, Nonlinear Anal. Real World Appl. , 9 (2008), 2196–2213
2008
-
[23]
W. Z. Huang, Problem on minimum wave speed for a LotkaVolterr a reaction-diffusion com- petition model, J. Dynam. Differential Equations , 22 (2010), 285–297
2010
-
[24]
W. Z. Huang and M. A. Han, Non-linear determinacy of minimum w ave speed for a Lotka- Volterra competition model, J. Differential Equations , 251 (2011), 1549–1561
2011
-
[25]
M. Iida, R. Lui and H. Ninomiya, Stacked fronts for cooperati ve systems with equal diffusion coefficients, SIAM J. Math. Anal. , 43 (2011), 1369–1389
2011
-
[26]
A. N. Kolmogorov, I. G. Petrovsky and N. S. Piskunov, ´Etude de l´ equation de la diffusion avec croissance de la quantit´ e de mati´ ere et son applicati on ` a un probl´ eme biologique, Bull. Univ. Moscow, Ser. Internat., Sec. A , 1 (1937), 1-26
1937
-
[27]
M. A. Lewis, B. T. Li and H. F. W einberger, Spreading speed and linear determinacy for two-species competition models, J. Math. Biol. , 45 (2002), 219–233
2002
-
[28]
B. T. Li, Multiple invasion speeds in a two-species integro- difference competition model, J. Math. Biol. , 76 (2018), 1975–2009
2018
-
[29]
B. T. Li, H. F. W einberger and M. A. Lewis, Spreading speeds as slowest wave speeds for cooperative systems, J. Math. Bios. , 196 (2005), 82–98
2005
-
[30]
Liang and X
X. Liang and X. Q. Zhao, Asymptotic speeds of spread and trave ling waves for monotone semiflows with applications, Comm. Pure Appl. Math. , 60 (2007), 1–40
2007
-
[31]
Lin and W
G. Lin and W. T. Li, Asymptotic spreading of competition diffu sion systems: the role of interspecific competitions, European J. Appl. Math. , 23 (2012), 669–689. 32 QIAN LIU, SHUANG LIU AND KING-YEUNG LAM
2012
-
[32]
S. Y. Liu, H. M. Huang and M. X. W ang, Asymptotic spreading of a diffusive competition model with different free boundaries, J. Differential Equations , 266 (2019), 4769–4799
2019
-
[33]
A. J. Majda and P. E. Souganidis, Large-scale front dynamics for turbulent reaction-diffusion equations with separated velocity scales, Nonlinearity , 7 (1994), 1–30
1994
-
[34]
M´ endez, J
V. M´ endez, J. Fort, H. G. Rotstein and S. Fedotov, Speed of re action-diffusion fronts in spatially heterogeneous media, Phys. Rev. E , 68 (2003), 041105
2003
-
[35]
Roques, Y
L. Roques, Y. Hosono, O. Bonnefon and T. Boivin, The effect of c ompetition on the neutral intraspecific diversity of invasive species, J. Math. Biol. , 71 (2015), 465–489
2015
-
[36]
Shigesada and K
N. Shigesada and K. Kawasaki, Biological invasions: Theory and Practice , Oxford University Press, Oxford, 1997
1997
-
[37]
P. E. Souganidis, Front propagation: theory and applications , Lectures Notes in Math., 1660, Fond. CIME/CIME Found. Subser., Springer, Berlin, 1997
1997
-
[38]
M. M. Tang and P. C. Fife, Propagating fronts for competing sp ecies equations with diffusion, Arch. Rational Mech. Anal. , 73 (1980), 69–77
1980
-
[39]
A. I. Volpert, V. A. Volpert and V. A. Volpert, Traveling Wave Solutions of Parabolic Sys- tems, Transl. Math. Monogr., Amer. Math. Soc. 1994
1994
-
[40]
M. X. W ang and Y. Zhang, Note on a two-species competition-di ffusion model with two free boundaries, Nonlinear Anal. , 159 (2017), 458–467
2017
-
[41]
M. X. W ang and Y. Zhang, Dynamics for a diffusive preypredator model with different free boundaries, J. Differential Equations , 264 (2018), 3527–3558
2018
-
[42]
C. H. W u, The minimal habitat size for spreading in a weak comp etition system with two free boundaries, J. Differential Equations , 259 (2015), 873–897
2015
-
[43]
Xin, Front propagation in heterogeneous media, SIAM Rev
J. Xin, Front propagation in heterogeneous media, SIAM Rev. , 42 (2000), 161–230. Received xxxx 20xx; revised xxxx 20xx. E-mail address : liuqian0519@ruc.edu.cn E-mail address : liushuangnqkg@ruc.edu.cn E-mail address : lam.184@osu.edu
2000
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.