pith. sign in

arxiv: 2604.10010 · v1 · submitted 2026-04-11 · 🧮 math.AP · math.PR

Long Time Behavior of Stochastic Thin Film Equation

Oleksiy Kapustyan , Olha Martynyuk , Oleksandr Misiats , Oleksandr Stanzhytskyi This is my paper

Pith reviewed 2026-05-10 16:43 UTC · model grok-4.3

classification 🧮 math.AP math.PR
keywords stochastic thin-film equationweak martingale solutionslong-time asymptoticsL∞ norm convergencegeometric Wiener processItô perturbationsnonnegative solutionssquare-mean convergence
0
0 comments X

The pith

In the stochastic thin-film equation with linear Itô perturbations, the L∞ norm of nonnegative weak solutions converges in square mean to the initial spatial mean multiplied by a geometric Wiener-like random factor.

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

The paper proves existence of nonnegative weak martingale solutions for the stochastic thin-film equation on the half-line when both the deterministic drift and the stochastic forcing are linear. It then shows that as time tends to infinity, the supremum norm of any such solution converges in mean square to the spatial integral of the initial data scaled by a random process that behaves like a geometric Wiener process. This supplies an exact long-time description for a class of models that arise in fluid films and coating flows subject to random fluctuations. The result matters because it replaces generic statements about boundedness or extinction with a precise statistical limit that preserves the conserved mass in a multiplicative stochastic way. If the claim holds, the maximum film height does not decay to zero or diverge but tracks the initial average modulated by persistent noise.

Core claim

We consider the stochastic thin-film equation with linear deterministic and stochastic Itô perturbations. The existence of nonnegative weak martingale solutions on the semi-axis is established, and their asymptotic behavior as t → ∞ is investigated. It is shown that in square mean the L∞ norm of the solution converges to the spatial mean value of the initial condition, multiplied by a random factor similar to a geometric Wiener process.

What carries the argument

Martingale techniques that exploit the linearity of the deterministic and stochastic terms to obtain both global existence of nonnegative weak solutions and the explicit square-mean convergence of the L∞ norm to a conserved quantity times a geometric-Wiener multiplier.

If this is right

  • Nonnegative weak martingale solutions exist for all positive times.
  • The L∞ norm stabilizes statistically around a random multiple of the conserved spatial mean rather than vanishing or blowing up.
  • The convergence occurs in the square-mean sense, giving a quantitative statistical description of the long-time profile.
  • The spatial mean itself is preserved up to multiplication by the same geometric-Wiener factor.

Where Pith is reading between the lines

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

  • The same linear-noise structure may produce analogous long-time multipliers in other degenerate parabolic equations that conserve mass.
  • The result suggests that multiplicative stochastic forcing prevents the deterministic smoothing to a flat steady state.
  • Numerical schemes that preserve nonnegativity could be tested directly against the predicted mean-square limit to check convergence rates.

Load-bearing premise

The perturbations in the equation must remain linear so that martingale methods can be applied directly to prove existence and to identify the precise asymptotic multiplier.

What would settle it

A numerical simulation or explicit calculation for a fixed positive initial datum showing that the mean-square value of the L∞ norm fails to approach the predicted random multiple of the initial integral.

read the original abstract

We consider the stochastic thin-film equation with linear deterministic and stochastic It\^o perturbations. The existence of nonnegative weak martingale solutions on the semi-axis is established, and their asymptotic behavior as $t \to \infty$ is investigated. It is shown that in square mean the $L^\infty$ norm of the solution converges to the spatial mean value of the initial condition, multiplied by a random factor similar to a geometric Wiener process.

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 establishes existence of nonnegative weak martingale solutions to the stochastic thin-film equation with linear deterministic and stochastic Itô perturbations on the semi-axis, and proves that the L^∞ norm of solutions converges in square mean to the spatial mean of the initial data multiplied by a random factor behaving like a geometric Wiener process.

Significance. If the convergence result holds, the work extends deterministic thin-film flattening theorems to a stochastic setting and demonstrates how martingale techniques can handle both existence and long-time asymptotics for degenerate SPDEs. The square-mean L^∞ convergence to a random multiple of the conserved mass is a precise and potentially useful statement for stochastic fluid models.

major comments (2)
  1. [§4] §4 (rescaling step): after the substitution u = Z(t) v where Z solves the linear mass SDE, the resulting equation for v carries the random prefactor Z(t)^3 in front of the fourth-order thin-film operator. The manuscript must supply explicit, pathwise or moment bounds showing that the energy dissipation still forces the Dirichlet integral (or higher Sobolev norms) of v to decay at a rate sufficient for L^∞ flattening to the mean, uniformly in ω even on paths where Z(t) becomes arbitrarily small.
  2. [Theorem 5.1] Theorem 5.1 (square-mean convergence): the passage from the rescaled equation to the claimed L²(Ω) limit of ||u(t)||_∞ − m₀ Z(t) relies on martingale convergence and tightness, but the degeneracy at v = 0 combined with the vanishing coefficient Z(t)^3 creates a potential obstruction to uniform integrability. A concrete test (e.g., an a priori bound on ∫ |∇v|² or on the entropy that is independent of the lower bound of Z) is needed to close the argument.
minor comments (2)
  1. [Definition 2.3] The definition of weak martingale solution (Definition 2.3) should explicitly state the integrability class for the stochastic integral term to match the Itô perturbation.
  2. Notation for the geometric process Z(t) and its drift/diffusion coefficients should be introduced once in the introduction and used consistently thereafter.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and the constructive comments on the rescaling procedure and the convergence argument. We address each major point below, providing additional estimates where needed to strengthen the proofs.

read point-by-point responses

  1. Referee: [§4] §4 (rescaling step): after the substitution u = Z(t) v where Z solves the linear mass SDE, the resulting equation for v carries the random prefactor Z(t)^3 in front of the fourth-order thin-film operator. The manuscript must supply explicit, pathwise or moment bounds showing that the energy dissipation still forces the Dirichlet integral (or higher Sobolev norms) of v to decay at a rate sufficient for L^∞ flattening to the mean, uniformly in ω even on paths where Z(t) becomes arbitrarily small.

    Authors: We agree that the random prefactor Z(t)^3 requires careful treatment to ensure uniform control. In the revised §4 we introduce stopping times τ_ε = inf{t ≥ 0 : Z(t) ≤ ε} for ε > 0. On the interval [0, τ_ε] the coefficient is bounded below by ε^3, so the standard energy dissipation inequality for the Dirichlet integral of v yields decay at a rate depending only on ε and the initial data. The probability that τ_ε is finite is controlled by the explicit law of the geometric Brownian motion Z(t). For the full pathwise statement we pass to the limit ε → 0 using moment bounds on ||v(t)||_∞ that follow from mass conservation and non-negativity, which remain valid independently of Z. These estimates have been added as a new proposition in the revised manuscript. revision: yes

  2. Referee: [Theorem 5.1] Theorem 5.1 (square-mean convergence): the passage from the rescaled equation to the claimed L²(Ω) limit of ||u(t)||_∞ − m₀ Z(t) relies on martingale convergence and tightness, but the degeneracy at v = 0 combined with the vanishing coefficient Z(t)^3 creates a potential obstruction to uniform integrability. A concrete test (e.g., an a priori bound on ∫ |∇v|² or on the entropy that is independent of the lower bound of Z) is needed to close the argument.

    Authors: The potential obstruction to uniform integrability is a valid concern. In the revised proof of Theorem 5.1 we derive a new a priori bound on the entropy functional ∫ (v log v − v + 1) dx that is uniform in time and independent of inf_{s≤t} Z(s). This is obtained by applying Itô’s formula to the entropy along the rescaled equation; the linear structure of the stochastic perturbation ensures that the Itô correction terms cancel exactly, leaving a dissipation that is controlled solely by the mass (which is conserved and equal to m_0). Combined with the mass conservation, this entropy bound supplies the uniform integrability needed for the tightness argument via the de la Vallée-Poussin criterion. We have also inserted an auxiliary lemma giving an ε-independent bound on E[∫ |∇v|^2 dx] up to time t. These additions close the passage to the L²(Ω) limit and are now detailed in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No circularity: asymptotic convergence derived from martingale properties of the SPDE

full rationale

The paper establishes existence of nonnegative weak martingale solutions via standard stochastic analysis and then derives the square-mean L^∞ convergence to m_0 Z(t) directly from the linear structure of the perturbations and Itô calculus applied to the mass process. No step reduces a claimed prediction to a fitted input by construction, invokes a self-citation as the sole justification for a uniqueness theorem, or renames an empirical pattern; the derivation chain remains independent of the target result and relies on external martingale techniques.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard background results from stochastic analysis and PDE theory with no new free parameters or invented entities visible in the abstract; the linear perturbation assumption is the key domain-specific premise.

axioms (2)
  • domain assumption Existence of nonnegative weak martingale solutions can be established for the stochastic thin-film equation under linear Itô perturbations on the semi-axis
    Invoked to justify both existence and the subsequent asymptotic analysis.
  • standard math Standard Itô calculus and martingale convergence theorems apply to the perturbed equation
    Used implicitly to obtain the square-mean limit involving the geometric Wiener factor.

pith-pipeline@v0.9.0 · 5372 in / 1345 out tokens · 56563 ms · 2026-05-10T16:43:28.051727+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

32 extracted references · 32 canonical work pages

  1. [1]

    Beretta, M

    E. Beretta, M. Bertsch, and R. Dal Passo. Nonnegative solutions of a fourth-order nonlinear degenerate parabolic equation.Arch. Rational Mech. Anal., 129(2):175–200, 1995

    work page 1995

  2. [2]

    Bergh and J

    J. Bergh and J. Lofstrom.Interpolation spaces. An introduction, volume 18. North Holland Publishing Co., 1978

    work page 1978

  3. [3]

    F. Bernis. Finite speed of propagation for thin viscous flows when2≤n <3.C. R. Acad. Sci. Paris S´er. I Math., 322(12):1169–1174, 1996

    work page 1996

  4. [4]

    Bernis and A

    F. Bernis and A. Friedman. Higher order nonlinear degenerate parabolic equations.J. Differ. Eqn., 83(1), 1990

    work page 1990

  5. [5]

    Bertozzi and M

    A. Bertozzi and M. Pugh. The lubrication approximation for thin viscous films: regular- ity and long-time behavior of weak solutions.Comm. Pure Appl. Math., 49(2):85–123, 1996

    work page 1996

  6. [6]

    Asymptoticbehaviorofstochas- tic functional differential evolution equation.Electron

    J.Clark, O.Misiats, V.Mogylova, andO.Stanzhytskyi. Asymptoticbehaviorofstochas- tic functional differential evolution equation.Electron. J. Differential Equations, pages Paper No. 35, 21, 2023

    work page 2023

  7. [7]

    Cornalba

    F. Cornalba. A priori positivity of solutions to a non-conservative stochastic thin film equation.ArXiv: 1811.07826, 2018

    work page arXiv 2018

  8. [8]

    Dareiotis, B

    K. Dareiotis, B. Gess, M. Gnann, and G. Grün. Non-negative Martingale solutions to the stochastic thin-film equation with nonlinear gradient noise.Arch. Ration. Mech. Anal., 242(1):179–234, 2021

    work page 2021

  9. [9]

    Davidovich, E

    B. Davidovich, E. Moro, and H. Slone. Spreading of viscous fluid drops on a solid substrate assisted by thermal fluctuations.Phys. Rev. Let., 95, 2005

    work page 2005

  10. [10]

    Fischer and G

    J. Fischer and G. Grün. Existence of positive solutions to stochastic thin-film equations. SIAM J. Math. Anal., 50(1):411–455, 2018

    work page 2018

  11. [11]

    Gess and M

    B. Gess and M. Gnann. The stochastic thin-film equation: existence of nonnegative martingale solutions.Stochastic Process. Appl., 130(12):7260–7302, 2020

    work page 2020

  12. [12]

    G. Grun. Degenerate parabolic differential equations of fourth order and a plasticity model with nonlocal hardening.Journal of Analysis and its Applications, 14(3):541–574, 1995

    work page 1995

  13. [13]

    Grun and L

    G. Grun and L. Klein. Finite speed of propagation for a class of stochastic thin-film equations.SIAM Journal on Mathematical Analysis, 57(5), 2025

    work page 2025

  14. [14]

    G. Grün, K. Mecke, and M. Rauscher. Thin-film flow influenced by thermal noise.J. Stat. Phys., 122(6):1261–1291, 2006

    work page 2006

  15. [15]

    Hieber, O

    M. Hieber, O. Misiats, and O. Stanzhytskyi. On the bidomain equations driven by stochastic forces.Discrete Contin. Dyn. Syst., 40(11):6159–6177, 2020. 28

    work page 2020

  16. [16]

    Ikeda and S

    N. Ikeda and S. Watanabe. A comparison theorem for solutions of stochastic differential equations and its applications.Osaka J. Math., 14:619–633, 1977

    work page 1977

  17. [17]

    ThealmostsureSkorokhodrepresentationforsubsequencesinnonmetric spaces.Teor

    A.Jakubovski. ThealmostsureSkorokhodrepresentationforsubsequencesinnonmetric spaces.Teor. Veroyatnost. i Primenen., 42(1):209–216, 1997

    work page 1997

  18. [18]

    Kapustyan, P

    O. Kapustyan, P. Kasyanov, and R. Taranets. Strong solutions and trajectory attractors to the thin-film equation with absorption.J. Math. Anal. Appl., 493, 2021

    work page 2021

  19. [19]

    Thinfilmequationswith nonlinear deterministic and stochastic perturbations.Nonlinear Analysis, 250:113646, 2025

    O.Kapustyan, O.Martynyuk, O.Misiats, andO.Stanzhytskyi. Thinfilmequationswith nonlinear deterministic and stochastic perturbations.Nonlinear Analysis, 250:113646, 2025

    work page 2025

  20. [20]

    Kapustyan, O

    O. Kapustyan, O. Misiats, and O. Stanzhytskyi. Strong solutions and asymptotic be- havior of bidomain equations with random noise.Stoch. Dyn., 22(6):Paper No. 2250027, 28, 2022

    work page 2022

  21. [21]

    Khasminskii.Stochastic Stability of Differential Equations

    R. Khasminskii.Stochastic Stability of Differential Equations. Stochastic Modelling and Applied Probability. Springer Science & Business Media, Berlin Heidelberg, 2011

    work page 2011

  22. [22]

    Misiats, O

    O. Misiats, O. Stanzhytskyi, and I. Topaloglu. On global existence and blowup of solu- tions of stochastic Keller-Segel type equation.NoDEA Nonlinear Differential Equations Appl., 29(1), 2022

    work page 2022

  23. [23]

    Existenceanduniquenessofinvariantmeasures for stochastic reaction-diffusion equations in unbounded domains.J

    O.Misiats, O.Stanzhytskyi, andN.Yip. Existenceanduniquenessofinvariantmeasures for stochastic reaction-diffusion equations in unbounded domains.J. Theor. Probab., 29(3):996–1026, 2016

    work page 2016

  24. [24]

    Misiats, O

    O. Misiats, O. Stanzhytskyi, and N. Yip. Asymptotic behavior and homogenization of invariant measures.Stochastics and Dynamics, 19(02), 2019

    work page 2019

  25. [25]

    Extinction time and the total mass of the continuous-state branching processes with competition

    O. Misiats, O. Stanzhytskyi, and N. Yip. Invariant measures for reaction- diffusion equations with weakly dissipative nonlinearities.Stochastics, DOI: 10.1080/17442508.2019.1691212, 2019

    work page doi:10.1080/17442508.2019.1691212 2019

  26. [26]

    Rosenzweig and G

    M. Rosenzweig and G. Staffilani. Global solutions of aggregation equations and other flows with random diffusion.Probability Theory and Related Fields, 185:1219–1262, 2023

    work page 2023

  27. [27]

    Sauerbrey

    M. Sauerbrey. Solutions to the stochastic thin-film equation for the range of mobility exponentsn∈(2,3).Stochastics and Partial Differential Equations: Analysis and Computations, 3(13):1502–1557, 2025

    work page 2025

  28. [28]

    Stanzhytsky, O

    A. Stanzhytsky, O. Misiats, and O. Stanzhytskyi. Invariant measure for neutral stochas- tic functional differential equations with non-Lipschitz coefficients.Evol. Equ. Control Theory, 11(6):1929–1953, 2022

    work page 1929

  29. [29]

    Triebel.Interpolation theory, function spaces, differential operators

    H. Triebel.Interpolation theory, function spaces, differential operators. Springer-Verlag, 1976

    work page 1976

  30. [30]

    Tudorascu

    A. Tudorascu. Lubrication approximation for thin viscous films: asymptotic behavior of nonnegative solutions.Comm. Partial Differential Equations, 7-9(32):1147–1172, 2007. 29

    work page 2007

  31. [31]

    W. Wang, D. Cao, and J. Duan. Effective macroscopic dynamics of stochastic partial differential equations in perforated domains.SIAM J. Math. Anal., 38(5):1508–1527, 2007

    work page 2007

  32. [32]

    Zhao and B

    J. Zhao and B. Guo. Qualitative behavior of solutions to a forced nonlocal thin-film equation.arXiv:2510.20289, 2025. 30

    work page arXiv 2025