pith. sign in

arxiv: 2604.19595 · v1 · submitted 2026-04-21 · 🧮 math.AP

Shock wavefronts for parabolic equations with sign-changing diffusivity

Diego Berti , Andrea Corli , Luisa Malaguti This is my paper

Pith reviewed 2026-05-10 01:51 UTC · model grok-4.3

classification 🧮 math.AP
keywords reaction-diffusion equationsign-changing diffusivityshock wavefrontstraveling wavesbistable nonlinearitydiscontinuous solutionspopulation dynamics
0
0 comments X

The pith

Reaction-diffusion models with sign-changing diffusivity admit families of discontinuous shock wavefronts.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines a one-dimensional parabolic equation whose diffusion coefficient switches from positive to negative and back to positive, paired with a bistable reaction whose middle equilibrium sits inside the negative-diffusion interval. Under these conditions the model admits no continuous traveling waves that connect the two stable equilibria. The authors instead construct a family of traveling profiles that contain a single jump discontinuity and determine the admissible propagation speeds for these shocks. The construction supplies explicit conditions on the jump location and speed that ensure the profile connects the equilibria while satisfying the equation in the classical sense on either side of the jump. The same framework is shown to cover a population model that distinguishes isolated and grouped individuals.

Core claim

For the given class of equations, continuous traveling waves connecting the stable states 0 and 1 do not exist, yet there exist shock wavefronts whose profiles possess exactly one jump discontinuity. These profiles are obtained by solving the traveling-wave ODE separately in the regions of positive diffusivity, imposing suitable transmission conditions at the jump, and selecting the speed so that the overall solution advances from one equilibrium to the other.

What carries the argument

Shock wavefront: a traveling-wave profile containing a single jump discontinuity that satisfies the PDE classically away from the discontinuity and an integral balance condition across the jump.

If this is right

  • Each admissible shock wavefront propagates at a speed fixed by the location of its jump and the coefficients of the reaction and diffusion functions.
  • The non-existence of continuous waves follows directly once the middle equilibrium enters the negative-diffusivity interval.
  • The discontinuous profiles can be matched to the population model that distinguishes isolated and grouped individuals, yielding concrete predictions for front speeds in that setting.
  • The jump location must remain inside the negative-diffusivity interval for the construction to close.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Numerical time-dependent simulations initialized near a predicted jump profile could confirm whether the discontinuous wave persists and travels at the analytically predicted speed.
  • The same jump-matching technique might extend to equations with multiple sign changes or to two-dimensional domains with strip-like negative-diffusivity regions.
  • Stability of these shock fronts under small perturbations could be tested by linearizing the evolution operator around the constructed profile.
  • If the negative-diffusivity interval shrinks to zero width, the shock speed is expected to approach the speed of the classical bistable front for the limiting positive-diffusivity equation.

Load-bearing premise

The interior zero of the bistable reaction term must lie inside the interval where the diffusivity is negative.

What would settle it

An explicit continuous traveling-wave solution connecting 0 to 1 for any choice of parameters satisfying the sign-change and bistable-zero conditions would falsify the non-existence claim.

Figures

Figures reproduced from arXiv: 2604.19595 by Andrea Corli, Diego Berti, Luisa Malaguti.

Figure 1
Figure 1. Figure 1: Left: plot of the potential P. Right: plots of the diffusivity D (solid line) and the source-sink term g (dashed line). Assumption (P) makes (1.1) a forward-backward-forward equation; we allow P ′ = D to possibly vanish at either 0 or 1. We recall that a source term g satisfying (g) is called bistable, and g is called monostable if it is positive in (0, 1) and vanishes at 0 and 1. Equation (1.1) is known t… view at source ↗
Figure 2
Figure 2. Figure 2: Above: the case P(1) > P(0) and the interval I. Middle: the case P(1) = P(0). Below: the case P(1) = P(γ) = P(0). where ξs is the jump point. All terms appearing in the previous estimate exist and the condition requires that the quantity P ′ (φ) + cφ is continuous across the shock (see Proposi￾tion 3.1). Thus, it closely resembles the classical Rankine–Hugoniot condition for shocks in the hyperbolic settin… view at source ↗
Figure 3
Figure 3. Figure 3: The profile φ when P(β) ≤ P(0) < P(1) (left) and P(0) < P(1) ≤ P(α) (right). (φℓ , φr) fulfilling (1.7) (see [31, Proposition 2]). However, such a pair may fail to exist if P(α) > P(1) as we will show in Proposition 2.2 (see also Remark 2.1 and [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Condition (1.7) states that the area of the yellow region equals that of the blue [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Plots of P. Red points denote α and β, the blue point is the point ( 2 3 , 0), black points are the extrema of the intervals I and η(I). On the left: Di = 35, Dg = 8, so that ω = 1/3. On the right, Di = 32, Dg = 5 so that ω = 2/3. Proof. We rewrite formula (2.1) by exploiting (2.3); we obtain D(u) = 3(Di − Dg)  u 2 − (α + β)u + αβ = (Di − Dg) [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Plots of the functions η and P, D. On the left column the plot of η, on the right the plot of D vs. P, where the red points denote P(β) and P(α). In the top line the case 0 < ω < 1 2 , in the bottom line the case √ 1 3 < ω < 1. Notation and values of Di , Dg are as in [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The values of φℓ for which condition (2.10) is satisfied are shown in orange. The values of φℓ that minimize the corresponding speed are shown in blue and yellow. Finally, we provide an estimate of the closed interval where the speeds vary. The result makes use of (5.8) (see Proposition 5.3), derived in the general case. Proposition 2.4. Consider (2.1)-(2.2), with Di > 4Dg > 0 and (2.4). Then, c ∗ (φℓ) sat… view at source ↗
Figure 8
Figure 8. Figure 8: Pictorial representation of the results in Propositions 4.1 and 4.2. The arrows [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The traveling waves ˆφ and ˇφ. Lemma 5.2. Assume c ∈ R. Let φ be defined as in (5.2) for some ξs in R. If both conditions in (3.7) are satisfied, then φ is the unique profile with speed c of a shock wavefront with a single jump from φr to φℓ located at ξs. Proof. Let φ be defined as in (5.2) for some c, ξs ∈ R. Assume conditions (3.7). If (φℓ , φr) = (0, 1), then (3.7)1 implies P(1) = P(0) and (3.7)2 yield… view at source ↗
Figure 10
Figure 10. Figure 10: The geometric meaning of condition (5.1). The speed [PITH_FULL_IMAGE:figures/full_fig_p029_10.png] view at source ↗
read the original abstract

We consider a reaction-diffusion equation in a one-dimensional space, where the diffusion coefficient changes sign from positive to negative and back to positive. The reaction term is bistable, with its interior zero located in the region where the diffusivity is negative. The model does not admit continuous wavefronts, i.e., continuous traveling waves that connect the steady states $0$ and $1$. We prove the existence of a family of shock wavefronts, that is, wavefronts with profiles exhibiting a jump discontinuity. We investigate the properties of these profiles and their propagation speeds. Finally, we apply the results to a recently proposed model describing the movement of a population composed of both isolated and grouped individuals.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript considers a reaction-diffusion equation in one space dimension with a diffusivity coefficient that changes sign from positive to negative and back to positive. The reaction term is bistable, with its interior zero located inside the negative-diffusivity interval. The paper proves that no continuous traveling waves connecting the equilibria 0 and 1 exist, establishes the existence of a family of discontinuous shock wavefronts, analyzes the properties of these profiles together with their propagation speeds, and applies the results to a population model with isolated and grouped individuals.

Significance. If the central claims hold, the work introduces a new class of discontinuous traveling-wave solutions for parabolic equations with sign-changing diffusion. The geometric condition that the bistable zero lies in the negative-diffusivity region supplies a clear criterion separating the non-existence of continuous waves from the existence of shocks. The explicit construction and speed characterization, together with the concrete application, add practical value. The rigorous existence proof for this family of solutions is a notable strength.

major comments (2)
  1. [Hypotheses and main theorems] The location of the interior zero of the bistable reaction term inside the negative-diffusivity interval is load-bearing for both the non-existence of continuous wavefronts and the construction of admissible jumps. The manuscript must state the precise parameter ranges guaranteeing this location and verify that the condition remains satisfied throughout the family of solutions and in the population-model application; otherwise both the exclusion argument and the jump-matching conditions lose their justification.
  2. [Existence proof for shock wavefronts] The derivation of the jump conditions from the weak form of the PDE across the discontinuity must be carried out explicitly, showing that the selected speed satisfies the integrated balance and that the resulting profile is consistent with the sign-changing diffusivity on either side of the jump.
minor comments (2)
  1. Notation for the diffusivity function and the reaction term should be made uniform across all sections to avoid ambiguity when referring to the positive and negative intervals.
  2. All figures illustrating the profiles and the sign-changing diffusivity should include clear labels for the jump location and the intervals of positive/negative diffusion.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The suggestions help clarify the hypotheses and strengthen the presentation of the existence proof. We address each major comment below and indicate the revisions made.

read point-by-point responses

  1. Referee: [Hypotheses and main theorems] The location of the interior zero of the bistable reaction term inside the negative-diffusivity interval is load-bearing for both the non-existence of continuous wavefronts and the construction of admissible jumps. The manuscript must state the precise parameter ranges guaranteeing this location and verify that the condition remains satisfied throughout the family of solutions and in the population-model application; otherwise both the exclusion argument and the jump-matching conditions lose their justification.

    Authors: We agree that explicit parameter ranges strengthen the justification. The manuscript already assumes in the statement of Theorem 1.1 and the setup of Section 2 that the interior zero lies in the negative-diffusivity interval. In the revision we add Subsection 2.3, which specifies the precise ranges on the diffusivity coefficients and the bistable nonlinearity parameters that guarantee this location. We further verify that the family of shock profiles constructed in Theorem 3.1 satisfies the condition uniformly (the jump is placed outside the negative interval), and that the parameter choices in the population-model application of Section 5 keep the zero inside the negative region. These additions make the load-bearing role of the geometric condition fully transparent. revision: yes

  2. Referee: [Existence proof for shock wavefronts] The derivation of the jump conditions from the weak form of the PDE across the discontinuity must be carried out explicitly, showing that the selected speed satisfies the integrated balance and that the resulting profile is consistent with the sign-changing diffusivity on either side of the jump.

    Authors: We concur that an explicit derivation improves rigor. In the revised Section 3.2 we now derive the jump conditions directly from the weak formulation: we integrate the PDE across an arbitrary interval containing the discontinuity, obtain the Rankine-Hugoniot-type balance that determines the admissible speed c, and verify that the resulting profile is consistent with positive diffusivity on both sides while the jump location avoids the negative-diffusivity interval. This calculation confirms that the constructed family satisfies the integrated balance and is admissible. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the existence proof for shock wavefronts

full rationale

The paper is a self-contained mathematical existence proof for discontinuous traveling-wave solutions (shock wavefronts) to a reaction-diffusion PDE with sign-changing diffusivity. The key hypothesis—that the interior zero of the bistable nonlinearity lies strictly inside the negative-diffusivity interval—is introduced as an explicit modeling assumption and is used once to rule out continuous profiles and again to construct admissible jumps; it is not derived from or equivalent to the claimed solutions. No parameters are fitted to data and then relabeled as predictions, no ansatz is smuggled via self-citation, and the derivation proceeds by standard phase-plane analysis, energy estimates, and weak-form matching conditions that remain independent of the final existence statement. The application to a population model is a straightforward extension rather than a load-bearing justification.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The result rests on standard existence and comparison principles for ODEs and parabolic PDEs together with the explicit sign conditions on D and f; no free parameters or invented entities are introduced in the abstract.

axioms (1)
  • standard math Standard comparison and maximum principles for parabolic equations hold under the given regularity assumptions on D and f.
    Invoked implicitly to rule out continuous waves and to construct discontinuous ones.

pith-pipeline@v0.9.0 · 5406 in / 1231 out tokens · 49518 ms · 2026-05-10T01:51:51.313885+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [1]

    D. G. Aronson, H. F. Weinberger. Multidimensional nonlinear diffusion arising in population genetics.Adv. in Math., 30(1):33–76, 1978

    work page 1978

  2. [2]

    Berti, A

    D. Berti, A. Corli, L. Malaguti. Uniqueness and nonuniqueness of fronts for degenerate diffusion- convection reaction equations.Electron. J. Qual. Theory Differ. Equ., No. 66, pp. 1–34, 2020

    work page 2020

  3. [3]

    Berti, A

    D. Berti, A. Corli, L. Malaguti. Wavefronts for degenerate diffusion-convection reaction equa- tions with sign-changing diffusivity.Discrete Contin. Dyn. Syst., 41(12):6023–6046, 2021

    work page 2021

  4. [4]

    Berti, A

    D. Berti, A. Corli, L. Malaguti. Diffusion-convection reaction equations with sign-changing diffusivity and bistable reaction term.Nonlinear Anal. Real World Appl., 67:Paper No. 103579, 29, 2022

    work page 2022

  5. [5]

    Berti, A

    D. Berti, A. Corli, L. Malaguti. The role of convection in the existence of wavefronts for biased movements.Math. Methods Appl. Sci., 47(1):491–515, 2024

    work page 2024

  6. [6]

    Bertsch, R

    M. Bertsch, R. Dal Passo. Hyperbolic phenomena in a strongly degenerate parabolic equation. Arch. Rational Mech. Anal., 117(4):349–387, 1992

    work page 1992

  7. [7]

    Bressan.Hyperbolic systems of conservation laws

    A. Bressan.Hyperbolic systems of conservation laws. Oxford University Press, 2000

    work page 2000

  8. [8]

    Brezis.Functional analysis, Sobolev spaces and partial differential equations

    H. Brezis.Functional analysis, Sobolev spaces and partial differential equations. Universitext. Springer, New York, 2011

    work page 2011

  9. [9]

    Calvo, J

    J. Calvo, J. Campos, V. Caselles, O. S´ anchez, J. Soler. Flux-saturated porous media equations and applications.EMS Surv. Math. Sci., 2(1):131–218, 2015

    work page 2015

  10. [10]

    Calvo, J

    J. Calvo, J. Campos, V. Caselles, O. S´ anchez, J. Soler. Qualitative behaviour for flux-saturated mechanisms: travelling waves, waiting time and smoothing effects.J. Eur. Math. Soc. (JEMS), 19(2):441–472, 2017. 29

    work page 2017

  11. [11]

    Campos, A

    J. Campos, A. Corli, L. Malaguti. Saturated fronts in crowds dynamics.Adv. Nonlinear Stud., 21(2):303–326, 2021

    work page 2021

  12. [12]

    Caselles

    V. Caselles. On the entropy conditions for some flux limited diffusion equations.J. Differential Equations, 250(8):3311–3348, 2011

    work page 2011

  13. [13]

    Caselles

    V. Caselles. Flux limited generalized porous media diffusion equations.Publ. Mat., 57(1):155– 217, 2013

    work page 2013

  14. [14]

    Chertock, A

    A. Chertock, A. Kurganov, P. Rosenau. On degenerate saturated-diffusion equations with con- vection.Nonlinearity, 18(2):609–630, 2005

    work page 2005

  15. [15]

    R. M. Colombo, A. Corli. Continuous dependence in conservation laws with phase transitions. SIAM J. Math. Anal., 31(1):34–62, 1999

    work page 1999

  16. [16]

    C. M. Dafermos.Hyperbolic conservation laws in continuum physics. Springer-Verlag, Berlin, third edition, 2010

    work page 2010

  17. [17]

    Dr´ abek, S

    P. Dr´ abek, S. Jung, E. Ko, M. Zahradn´ ıkov´ a. Traveling waves for monostable reaction-diffusion- convection equations with discontinuous density-dependent coefficients.J. Math. Anal. Appl., 539(1, part 1):Paper No. 128481, 26, 2024

    work page 2024

  18. [18]

    Dr´ abek, M

    P. Dr´ abek, M. Zahradn´ ıkov´ a. Bistable equation with discontinuous density dependent diffusion with degenerations and singularities.Electron. J. Qual. Theory Differ. Equ., pages Paper No. 61, 16, 2021

    work page 2021

  19. [19]

    Dr´ abek, M

    P. Dr´ abek, M. Zahradn´ ıkov´ a. Traveling waves for generalized Fisher-Kolmogorov equation with discontinuous density dependent diffusion.Math. Methods Appl. Sci., 46(11):12064–12086, 2023

    work page 2023

  20. [20]

    Ferracuti, C

    L. Ferracuti, C. Marcelli, F. Papalini. Travelling waves in some reaction-diffusion-aggregation models.Adv. Dyn. Syst. Appl., 4(1):19–33, 2009

    work page 2009

  21. [21]

    P. C. Fife.Mathematical aspects of reacting and diffusing systems, volume 28 ofLecture Notes in Biomathematics. Springer-Verlag, Berlin-New York, 1979

    work page 1979

  22. [22]

    B. H. Gilding, R. Kersner.Travelling waves in nonlinear diffusion-convection reaction. Birkh¨ auser Verlag, Basel, 2004

    work page 2004

  23. [23]

    Goodman, A

    J. Goodman, A. Kurganov, P. Rosenau. Breakdown in Burgers-type equations with saturating dissipation fluxes.Nonlinearity, 12(2):247–268, 1999

    work page 1999

  24. [24]

    Guarnotta, C

    U. Guarnotta, C. Marcelli. Traveling waves for highly degenerate and singular reaction-diffusion- advection equations with discontinuous coefficients.J. Differential Equations, 460:Paper No. 114075, 26, 2026

    work page 2026

  25. [25]

    S. T. Johnston, R. E. Baker, S. D. McElwain, M. J. Simpson. Co-operation, competition and crowding: a discrete framework linking Allee kinetic, nonlinear diffusion, shocks and sharp- fronted travelling waves.Sci. Rep., 7:42134, 2017

    work page 2017

  26. [26]

    Kurganov, P

    A. Kurganov, P. Rosenau. Effects of a saturating dissipation in Burgers-type equations.Comm. Pure Appl. Math., 50(8):753–771, 1997

    work page 1997

  27. [27]

    Kuzmin, S

    M. Kuzmin, S. Ruggerini. Front propagation in diffusion-aggregation models with bi-stable reaction.Discrete Contin. Dyn. Syst. Ser. B, 16(3):819–833, 2011

    work page 2011

  28. [28]

    Y. Li, P. van Heijster, R. Marangell, M. J. Simpson. Travelling wave solutions in a negative nonlinear diffusion-reaction model.J. Math. Biol., 81(6-7):1495–1522, 2020

    work page 2020

  29. [29]

    Y. Li, P. van Heijster, M. J. Simpson, M. Wechselberger. Shock-fronted travelling waves in a reaction–diffusion model with nonlinear forward–backward–forward diffusion.Physica D, 423, 2021. 30

    work page 2021

  30. [30]

    Miller, A

    T. Miller, A. K. Y. Tam, R. Marangell, M. Wechselberger, B. H. Bradshaw-Hajek. Analytic shock-fronted solutions to a reaction-diffusion equation with negative diffusivity.Stud. Appl. Math., 153(1):Paper No. e12685, 23, 2024

    work page 2024

  31. [31]

    Miller, A

    T. Miller, A. K. Y. Tam, R. Marangell, M. Wechselberger, B. H. Bradshaw-Hajek. Shock selection in reaction–diffusion equations with partially negative diffusivity using nonlinear reg- ularisation.Phys. D, 474:Paper No. 134561, 2025

    work page 2025

  32. [32]

    Schwartz.Th´ eorie des distributions

    L. Schwartz.Th´ eorie des distributions. Tome I. Hermann & Cie, Paris, 1950

    work page 1950

  33. [33]

    T. P. Witelski. Shocks in nonlinear diffusion.Appl. Math. Lett., 8(5):27–32, 1995. 31

    work page 1995