pith. sign in

arxiv: 2506.18741 · v2 · submitted 2025-06-23 · 🧮 math.AP · math.PR

Free boundary regularity and well-posedness of physical solutions to the supercooled Stefan problem

Sebastian Munoz This is my paper

Pith reviewed 2026-05-19 08:03 UTC · model grok-4.3

classification 🧮 math.AP math.PR
keywords supercooled Stefan problemfree boundary regularityphysical solutionswell-posednessweighted obstacle problemjump discontinuities
0
0 comments X

The pith

Physical solutions to the supercooled Stefan problem have a free boundary that is C^1 in space and infinitely smooth except at isolated jump times.

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

This paper proves regularity and uniqueness results for physical solutions to the supercooled Stefan problem when the initial temperature is only assumed integrable. The free boundary is shown to be continuously differentiable in the space variable and infinitely differentiable away from a closed countable set of times, with no accumulation of those jump times. Short-time uniqueness is shown to imply global uniqueness, which extends well-posedness to initial data beyond previous regimes and resolves open questions from earlier works. A sympathetic reader would care because the results describe interface evolution in supercooled systems under minimal data assumptions.

Core claim

Physical solutions to the supercooled Stefan problem with only integrable initial temperature have a free boundary that is C^1 as a function of space and C^infty outside of a closed, countable set. The set of positive times when a jump occurs cannot have accumulation points. Short-time uniqueness of physical solutions implies global uniqueness, allowing uniqueness for very general initial data.

What carries the argument

A weighted obstacle problem satisfied by the physical solutions, used to establish regularity and non-degeneracy estimates and to classify the free boundary points.

If this is right

  • The free boundary is C^1 in space under only integrability assumptions on initial temperature.
  • Jump times in the free boundary cannot accumulate.
  • Short-time uniqueness of physical solutions yields global uniqueness.
  • Uniqueness holds for initial data outside the scope of prior well-posedness results.

Where Pith is reading between the lines

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

  • The weighted obstacle problem technique might extend regularity results to other parabolic free boundary problems with jumps.
  • The explicit classification of the countable set could support numerical methods that track singularity times in simulations.
  • Backward propagation of oscillation may apply to long-time analysis in related solidification models.
  • These uniqueness implications could guide stability studies for supercooled processes in materials science.

Load-bearing premise

The solutions are physical solutions that satisfy the weighted obstacle problem derived from the supercooled Stefan formulation.

What would settle it

A physical solution with integrable initial data whose free boundary fails to be C^1 in space or whose jump times accumulate would disprove the regularity and non-accumulation claims.

Figures

Figures reproduced from arXiv: 2506.18741 by Sebastian Munoz.

Figure 4.1
Figure 4.1. Figure 4.1: Short time construction of γ1(t). The zero curve z˜ − acts as a barrier that precludes z − from hitting the boundary at t = t0. Since u is smooth away from the boundary, t0 ≤ t < t1. If we had ux(t, x1 + M(t1 − t)) < 0, then t = t0 and we may simply take γ1(t) ≡ x1 + M(t1 − t). We may assume then that ux(t, x1 + M(t1 − t)) = 0. We define the curve γ1 in two steps: first on the interval [t, t1], and then … view at source ↗
Figure 4
Figure 4. Figure 4: ) [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗
Figure 4.2
Figure 4.2. Figure 4.2: Correspondence between jumps of Λ at t = t1, t2, . . . , tN and oscillations of u at t = t0. The curve x = γ(t) acts as a barrier keeping all the curves in a compact region. We can now prove the main result of this section by iterating Lemma 4.7. Proof of Theorem 1.3. Proceeding by contradiction, assume that t∞ ∈ (0, ∞) is an accumulation point for the set of discontinuities of Λ, and fix t0 < t∞. By Lem… view at source ↗
read the original abstract

We study the regularity and well-posedness of physical solutions to the supercooled Stefan problem. Assuming only that the initial temperature is integrable, we prove that the free boundary, known to have jump discontinuities as a function of the time variable, is $C^1$ as a function of the space variable, and is $C^{\infty}$ outside of a closed, countable set, which we describe explicitly. We also prove that, as conjectured in arXiv:1902.05174, the set of positive times when a jump occurs cannot have accumulation points. In addition, we prove that short-time uniqueness of physical solutions implies global uniqueness, which allows us to obtain uniqueness for very general initial data that fall outside the scope of the current well-posedness regime. In particular, we answer two questions left open in arXiv:1811.12356, arXiv:2302.13097, regarding the global uniqueness of solutions. We proceed by deriving a weighted obstacle problem satisfied by the solutions, which we exploit to establish regularity and non-degeneracy estimates and to classify the free boundary points. We also establish a backward propagation of oscillation property, which allows us to control the occurrence of future jumps in terms of the past oscillation of the solution.

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 claims that, for physical solutions to the supercooled Stefan problem with merely integrable initial temperature, the free boundary is C^1 as a function of the space variable and C^infty outside an explicitly described closed countable set. It further asserts that jump times cannot accumulate, that short-time uniqueness implies global uniqueness (resolving open questions from arXiv:1811.12356 and arXiv:2302.13097), and that these conclusions follow from deriving a weighted obstacle problem satisfied by the solutions together with a backward oscillation propagation property.

Significance. If the results hold, the work advances free-boundary regularity theory for the supercooled Stefan problem by obtaining C^1 spatial regularity and smoothness off a countable set under minimal L^1 data assumptions. The introduction of the weighted obstacle problem and the explicit classification of free-boundary points constitute concrete technical contributions that may apply to related parabolic obstacle problems. The global-uniqueness implication from short-time uniqueness is a useful structural observation.

major comments (2)
  1. [§2] §2 (derivation of the weighted obstacle problem): With initial temperature only in L^1, the temperature field is a priori merely integrable in space-time. The passage from the weak Stefan condition to the precise weighted obstacle formulation (involving specific test-function choices or integration by parts localized near the free boundary) is not obviously justified; additional approximation or density arguments appear necessary to avoid tacitly assuming continuity or boundedness that is recovered only later.
  2. [§4] §4 (non-degeneracy and classification): The C^1 regularity, C^infty regularity off the countable set, and the non-accumulation of jumps all rest on the weighted obstacle problem and the associated non-degeneracy estimates. If the derivation in §2 requires extra regularity not available from L^1 data alone, these estimates do not directly apply to the stated solution class.
minor comments (2)
  1. [Throughout] Notation for the free boundary and the weighted obstacle problem should be introduced once and used uniformly; occasional shifts between different symbols for the same object reduce readability.
  2. [§5] A short remark clarifying how the explicit description of the closed countable set follows from the oscillation propagation would help readers trace the countability argument.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for the positive assessment of the significance of our results. We address the two major comments below, clarifying the technical justifications and indicating revisions to strengthen the exposition.

read point-by-point responses

  1. Referee: [§2] §2 (derivation of the weighted obstacle problem): With initial temperature only in L^1, the temperature field is a priori merely integrable in space-time. The passage from the weak Stefan condition to the precise weighted obstacle formulation (involving specific test-function choices or integration by parts localized near the free boundary) is not obviously justified; additional approximation or density arguments appear necessary to avoid tacitly assuming continuity or boundedness that is recovered only later.

    Authors: We agree that the derivation merits additional detail to be fully rigorous under L^1 data. In §2 we begin from the distributional form of the Stefan condition satisfied by physical solutions and approximate by mollification of the initial temperature. For the smoothed data the weighted obstacle problem follows by direct integration by parts against localized test functions. The limit passage is justified by the stability of physical solutions under L^1 convergence together with the fact that the free boundary has locally finite perimeter (already established from the weak formulation). This procedure uses only integrability and avoids any a priori continuity or boundedness. We will revise §2 to include an explicit subsection spelling out the approximation scheme, the choice of test functions, and the passage to the limit. revision: yes

  2. Referee: [§4] §4 (non-degeneracy and classification): The C^1 regularity, C^infty regularity off the countable set, and the non-accumulation of jumps all rest on the weighted obstacle problem and the associated non-degeneracy estimates. If the derivation in §2 requires extra regularity not available from L^1 data alone, these estimates do not directly apply to the stated solution class.

    Authors: Because the weighted obstacle problem is established for the L^1 class via the approximation argument outlined above, the non-degeneracy estimates and point classification in §4 apply directly. These estimates are obtained by testing the obstacle formulation with suitable cut-off functions and using only the L^1 integrability to control oscillations; the backward propagation of oscillation is likewise derived from the variational inequality without presupposing higher regularity. The C^1 and C^infty conclusions are obtained after the formulation is in place, so there is no circularity. We will add a short clarifying paragraph in §4 that explicitly ties the estimates back to the L^1 derivation in §2. revision: partial

Circularity Check

0 steps flagged

No circularity; derivation self-contained via new obstacle problem

full rationale

The paper derives a weighted obstacle problem directly from the supercooled Stefan formulation under L1 initial data and then applies it to obtain the stated regularity, non-degeneracy, and uniqueness results. No quoted step reduces a prediction or central claim to a fitted input, self-citation chain, or definitional tautology. Cited prior works address open questions or background but do not supply the target regularity statements or uniqueness theorems used here. The logical order (weak formulation to obstacle problem to estimates) is presented as forward and independent of the final conclusions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The claims rest on the integrability assumption for initial data and on the existence of physical solutions that satisfy the derived weighted obstacle problem; no free parameters or new entities are introduced.

axioms (1)
  • domain assumption Initial temperature belongs to L^1
    Stated explicitly as the sole assumption on the data for all regularity and uniqueness statements.

pith-pipeline@v0.9.0 · 5754 in / 1185 out tokens · 35931 ms · 2026-05-19T08:03:41.161591+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

28 extracted references · 28 canonical work pages · 1 internal anchor

  1. [1]

    Blanchet, J

    A. Blanchet, J. Dolbeault, R. Monneau.On the continuity of the time derivative of the solution to the parabolic obstacle problem with variable coefficients.J. Math. Pures Appl. (9) 85 (2006), no. 3, 371–414

    work page 2006

  2. [2]

    L. A. Caffarelli,The obstacle problem revisited.J. Fourier Anal. Appl. 4 (1998), no. 4-5, 383–402

    work page 1998

  3. [3]

    L. A. Caffarelli, A. Petrosyan, H. Shahgholian.Regularity of a free boundary in parabolic potential theory. J. Amer. Math. Soc. 17 (2004), no. 4, 827–869

    work page 2004

  4. [4]

    Nonlinearity 28 (2015), no

    Qualitative properties of solutions for the noisy integrate and fire model in computational neuroscience. Nonlinearity 28 (2015), no. 9, 3365–3388

    work page 2015

  5. [5]

    S. Choi, I. C. Kim, Y.-H. Kim, Existence for the Supercooled Stefan Problem in General Dimensions. arXiv:2402.17154 (2024)

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  6. [6]

    Cuchiero, S

    C. Cuchiero, S. Rigger, S. Svaluto-Ferro,Propagation of minimality in the supercooled Stefan problem. Ann. Appl. Probab. 33 (2023), no. 2, 1388–1418

    work page 2023

  7. [7]

    Delarue, S

    F. Delarue, S. Nadtochiy, M. Shkolnikov,Global solutions to the supercooled Stefan problem with blow-ups: regularity and uniqueness,Prob. Math. Phys. 3 (2022), no. 1, 171–213

    work page 2022

  8. [8]

    Delarue, J

    F. Delarue, J. Inglis, S. Rubenthaler, E. Tanré,Particle systems with a singular mean-field self-excitation. Application to neuronal networks.Stochastic Process. Appl. 125 (2015), no. 6, 2451–2492. 24

    work page 2015

  9. [9]

    Durrett,Probability—theory and examples.Camb

    R. Durrett,Probability—theory and examples.Camb. Ser. Stat. Probab. Math., 49

    work page

  10. [10]

    Duvaut,Résolution d’un problème de Stefan (fusion d’un bloc de glace à zéro degré).C

    G. Duvaut,Résolution d’un problème de Stefan (fusion d’un bloc de glace à zéro degré).C. R. Acad. Sci. Paris Sér. A-B 276 (1973), A1461–A1463

    work page 1973

  11. [11]

    Nassib, Y.-H

    G. Nassib, Y.-H. Kim, A. Z. Palmer,PDE methods for optimal Skorokhod embeddings.Calc. Var. Partial Differential Equations 58 (2019), no. 3, Paper No. 113, 31 pp

    work page 2019

  12. [12]

    Hambly, S

    B. Hambly, S. Ledger, A. Søjmark.A McKean-Vlasov equation with positive feedback and blow-ups.Ann. Appl. Probab. 29 (2019), no. 4, 2338–2373

    work page 2019

  13. [13]

    Hambly, A

    B. Hambly, A. Søjmark.An SPDE model for systemic risk with endogenous contagion.Finance Stoch. 23 (2019), no. 3, 535–594

    work page 2019

  14. [14]

    Figalli, X

    A. Figalli, X. Ros-Oton, J. Serra.The singular set in the Stefan problem.J. Amer. Math. Soc. 37 (2024), no. 2, 305–389

    work page 2024

  15. [15]

    Inglis, D

    J. Inglis, D. Talay.Mean-field limit of a stochastic particle system smoothly interacting through threshold hitting-times and applications to neural networks with dendritic component.SIAMJ.Math.Anal.47(2015), no. 5, 3884–3916

    work page 2015

  16. [16]

    B. F. Jones, Jr.,A Class of Singular Integrals.Amer. J. Math. 86 (1964), 441–462

    work page 1964

  17. [17]

    I. C. Kim, Y.-H. Kim,The Stefan problem and free targets of optimal Brownian martingale transport. Ann. Appl. Probab. 34 (2024), no. 2, 2364–2414

    work page 2024

  18. [18]

    O. A. Ladyženskaja, V. A. Solonnikov, N. N. Ural’ceva.Linear and quasilinear equations of parabolic type. Translated from the Russian by S. Smith Transl. Math. Monogr., Vol. 23, American Mathematical Society, Providence, RI, 1968. xi+648 pp

    work page 1968

  19. [19]

    Ledger, A

    S. Ledger, A. Søjmark.Uniqueness for contagious McKean-Vlasov systems in the weak feedback regime. Bull. Lond. Math. Soc. 52 (2020), no. 3, 448–463

    work page 2020

  20. [20]

    Ledger, A

    S. Ledger, A. Søjmark.At the mercy of the common noise: blow-ups in a conditional McKean-Vlasov problem. Electron. J. Probab. 26 (2021), Paper No. 35, 39 pp

    work page 2021

  21. [21]

    G. M. Lieberman, Second order parabolic differential equations.World Scientific Publishing Co., Inc., River Edge, NJ, 1996. xii+439 pp

    work page 1996

  22. [22]

    V. P. Mikha˘ ılov,Partial differential equations.Translated from the Russian by P. C. Sinha “Mir”, Moscow; distributed by Imported Publications, Chicago, IL, 1978. 397 pp

    work page 1978

  23. [23]

    Nadtochiy, M

    S. Nadtochiy, M. Shkolnikov,Particle systems with singular interaction through hitting times: application in systemic risk modeling.Ann. Appl. Probab. 29 (2019), no. 1, 89–129

    work page 2019

  24. [24]

    Nadtochiy, M

    S. Nadtochiy, M. Shkolnikov,Mean field systems on networks, with singular interaction through hitting times. Ann. Probab. 48 (2020), no. 3, 1520–1556

    work page 2020

  25. [25]

    Revuz, M

    D. Revuz, M. Yor, Continuous martingales and Brownian motion. Third edition. Grundlehren Math. Wiss., 293 [Fundamental Principles of Mathematical Sciences] Springer-Verlag, Berlin, 1999. xiv+602 pp

    work page 1999

  26. [26]

    Sherman,A general one-phase Stefan problem.Quart

    B. Sherman,A general one-phase Stefan problem.Quart. Appl. Math. 28 (1970), 377–382. MR 0282082

    work page 1970

  27. [27]

    Mustapha, M

    S. Mustapha, M. Shkolnikov, Well-posedness of the supercooled Stefan problem with oscillatory initial conditions. Electron. J. Probab. 29 (2024), Paper No. 193, 21 pp

    work page 2024

  28. [28]

    Wang,On the Regularity Theory of Fully Nonlinear Parabolic Equations I.Comm

    L. Wang,On the Regularity Theory of Fully Nonlinear Parabolic Equations I.Comm. Pure Appl. Math. 45 (1992), 27–76. 25

    work page 1992