REVIEW 3 major objections 5 minor 30 references
Long-Time Dynamics and Modified Scattering for the Coulomb Vlasov--Hartree System
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper claims that the Coulomb Vlasov–Hartree system, a classical–quantum model of a Bose–Fermi mixture, has global smooth solutions for small localized data and that its long-time behavior is modified scattering in which the…
desk verdict First modified scattering for the coupled Coulomb Vlasov–Hartree system, with a genuine but easily patched regularity gap in the main theorem. 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 proof is carried out in the common self-similar variable $y=x/t$. For the Coulomb potential it uses an exact shell decomposition $1/|z|=C_0\int_0^\infty \chi(z/R)\,dR/R^2$, which lets the physical field $E(t,x)$ be compared with an effective field $E^0(t,x/t)$ whose evolution can be controlled through the transport equation for the free-transport profile $\gamma(t,x,v)=f(t,x+tv,v)$. On the Schrödinger side the central objects are the Galilean vector field $G=ix+2t\nabla$ and the renormalized profile $h(t,y)=t^{3/2}e^{i|y|^2t/4}u(t,ty)e^{i\theta(t,y)}$, where $\theta$ removes the growing Coulomb phase. The key structural fact is that both components propagate in the same variable: convergence of $h$ to $h_\infty$ gives the limiting force $F_\infty=-\nabla(V*|h_\infty|^2)$ that enters the Vlasov characteristics, and convergence of $\gamma$ gives the density $\rho_\infty$ that defines the phase $\phi_\infty=V*\rho_\infty$ in the Hartree law.
What would settle it
Solve (1.1) numerically in three dimensions for small Gaussian data with $V(x)=\lambda|x|^{-1}$ and check whether $\tilde\gamma(t,x,v)=\gamma(t,x-\log(t)F_\infty(v),v)$ converges in $L^\infty_{x,v}$ with the stated $t^{-\delta_0}$ rate and whether the rescaled Schrödinger profile, after removing the phase $e^{-i\log(t)\phi_\infty(x/t)}$, converges in $L^\infty_x$; failure of either for arbitrarily small data would refute Theorem 1.1. Alternatively, finding data satisfying (1.7) but not $C^1$ for which the system has no strong solution would refute the theorem's regularity statement.
Extended reading notes
Core claim
With $V(x)=\lambda|x|^{-1}$, the paper proves (Theorem 1.1) that for initial data satisfying the smallness condition (1.7) there is a unique global strong solution, and that as $t\to\infty$ the free-transport profile obeys $f(t, x+tv-\log(t)F_\infty(v), v) = f_\infty(x,v) + O_{L^\infty_{x,v}}(t^{-\delta_0})$, while the Hartree wave obeys $u(t,x) = t^{-3/2} e^{-i|x|^2/(4t)} e^{-i\log(t)\phi_\infty(x/t)} h_\infty(x/t) + O_{L^\infty_x}(t^{-3/2-\delta_0})$. Here $\rho_\infty(v)=\int f_\infty^2(x,v)\,dx$, $\phi_\infty = V*\rho_\infty$, and $F_\infty = -\nabla(V*|h_\infty|^2)$, so the asymptotic fermionic state sets the bosonic phase and the asymptotic bosonic state bends the fermionic trajectories. The same coupled argument, with integrable time weights, yields ordinary scattering for inverse-power interactions $V(x)=\lambda|x|^{-\alpha}$, $1<\alpha<3/2$, with rates $t^{1-\alpha}$.
Load-bearing premise
The whole long-time picture rests on the small-data bootstrap closing with the logarithmic growth rates stated in Proposition 3.10, on the validity of the unpublished companion's physical-space Hartree estimates [25] used as a black box, and on smallness in (1.7) implying the $C^1$ regularity that the local well-posedness argument requires.
Editorial extensions
If this is right
- The asymptotic state of the system is genuinely coupled: $f_\infty$ and $h_\infty$ cannot be read off from the free evolution of the initial data, because each carries a leading-order imprint of the other.
- The fermionic sector satisfies a modified scattering law along the corrected characteristics $x+tv-\log(t)F_\infty(v)$, so the long-time force on the fermions is entirely described by the limiting bosonic profile.
- The bosonic sector obeys $u(t,x)\approx t^{-3/2}e^{-i|x|^2/(4t)}e^{-i\log(t)\phi_\infty(x/t)}h_\infty(x/t)$, meaning the fermionic density acts as a logarithmic phase clock on the condensate.
- For $1<\alpha<3/2$, the same bootstrap yields ordinary (unmodified) scattering with the explicit rates (1.11)–(1.12), confirming that $\alpha=3/2$ is the borderline for long-range behavior in this coupled setting.
- Conservation laws imply the constructed asymptotic profiles preserve the initial fermionic and bosonic masses, so the scattering map, once defined, will be mass-preserving.
Reading between the lines
- A natural next step, not taken in the paper, is to construct the modified wave operator: Theorem 1.1 gives the asymptotic pair $(f_\infty,h_\infty)$ that a direct scattering map would need to hit, and the effective fields $F_\infty,\phi_\infty$ are exactly the data such an operator would invert.
- The common self-similar variable $y=x/t$ suggests the method may carry over to other kinetic–dispersive systems with borderline long-range interactions, such as Vlasov–Nordström or Vlasov–Schrödinger with gravitational coupling, where a similar two-sided log correction should appear.
- The restriction $\alpha<3/2$ is technical (an endpoint Sobolev embedding); a sharper embedding or a different norm might push the short-range result up to $\alpha=2$, and the failure at $\alpha=3/2$ would then be a genuine threshold rather than an artifact.
- One testable consequence of the coupling is that the asymptotic phase $\phi_\infty$ should equal the Coulomb potential of the limiting fermionic spatial density, equation (4.24); measuring the phase difference in a numerical two-species simulation would verify the self-consistent closure of the asymptotic data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the three-dimensional Vlasov--Hartree system with Coulomb interaction, in which the Vlasov force is generated by the bosonic density |u|^2 and the Hartree potential by the fermionic density \int f^2 dv. The main theorem (Theorem 1.1) asserts that sufficiently small, localized, and regular data produce a unique global strong solution and that both components exhibit modified scattering: f(t, x+tv-\log t F_\infty(v), v) -> f_\infty in L^\infty, while u(t,x) is asymptotic to t^{-3/2} e^{-i|x|^2/4t} e^{-i\log t \phi_\infty(x/t)} h_\infty(x/t). Theorem 1.2 states a short-range analogue with ordinary scattering for interactions |x|^{-\alpha}, 1<\alpha<3/2. The proof couples a bootstrap for the Vlasov free-transport profile \gamma with Galilean vector-field estimates for the Hartree component and derives the asymptotic profiles by constructing limits rather than assuming them.
Significance. If the main theorem is correct, this is the first sharp long-time asymptotic description for this coupled kinetic-dispersive model, and it shows that the two components interact at leading order through logarithmic phase and characteristic corrections. The proof is largely self-contained and contains parameter-free constructions: the asymptotic profiles are limits of the evolution, not fitted objects, and the Lorentz-space estimates used are mostly standard. The extension to short-range interactions is a useful additional result. However, the regularity gap in Theorem 1.1, the sign/phase inconsistency in Section 4.2, and the incomplete energy estimates in Section 3.2 must be fixed before the claims are fully supported.
major comments (3)
- [Section 1, Theorem 1.1; Section 3.1, Proposition 3.1] The hypotheses of Theorem 1.1 do not imply the regularity required by the local well-posedness theory. Assumption (1.7) controls \nabla_{x,v} f0 only in L^2 \cap L^\infty and with a weight in L^\infty; this does not imply that f0 is C^1, whereas Proposition 3.1 requires f0 \in C^1 and yields f \in C^1 on a short interval. Since the proof neither approximates f0 by C^1 data nor supplies a separate local existence argument under (1.7), the global C^1 statement in Theorem 1.1 is not justified. This is a genuine gap in the theorem as stated, although it is localized: adding f0 \in C^1 to (1.7), or weakening the conclusion for f to W^{1,\infty} \cap W^{1,2} in (x,v), would close it.
- [Section 4.2, Eqs. (4.3)-(4.7) and the paragraph preceding (1.10)] The proof of the Hartree asymptotic is inconsistent with the definition of h_\infty and with the sign of the phase in (1.10). Since h(t,y)=g(t,y)e^{i\theta(t,y)} and \theta(t,y)=l(y)+\log(t)\phi_\infty(y)+O(t^{-1/40}), the correct relation is g(t,y)=h(t,y)e^{-i\theta(t,y)}, so g converges to h_\infty(y)e^{-i(l(y)+\log(t)\phi_\infty(y))}. The displayed identity in the text instead writes g(t,y)-h_\infty(y)e^{+i\lambda(\log(t)\phi_\infty(y)+l(y))} and concludes convergence to the plus-phase profile. Consequently, as written, the proof yields a different phase than the e^{-i\log(t)\phi_\infty} in (1.10), and the l(y) phase is also missing from (1.10). This is fixable by redefining the final h_\infty to absorb e^{-il(y)} and by correcting the signs throughout Section 4.2, but the current text does not support the stated law.
- [Section 3.2, Lemmas 3.5 and 3.7] The energy estimates for Gu and G^2u omit the potential term proportional to \phi. From (i\partial_t-\Delta)Gu=\phi Gu+2t(\nabla\phi)u, the standard identity gives \partial_t\|Gu\|_{L^2} \le \|\phi\|_{L^\infty}\|Gu\|_{L^2}+2t\|\nabla\phi\|_{L^\infty}\|u\|_{L^2}. The lemma states only the second term. The missing term is handled by Gr\"onwall using Corollary 3.4 (\|\phi\|_{L^\infty} \lesssim \varepsilon^2\langle t\rangle^{-1}), which inserts an additional factor of t^{O(\varepsilon^2)} or \log^{O(\varepsilon^2)}; for \varepsilon small this is compatible with the bootstrap, but the estimate as printed is not justified. The same omission occurs in the G^2 equation in Lemma 3.7. Please provide the corrected inequality and verify that the log powers absorb the extra factor.
minor comments (5)
- [Eq. (4.16)] The definition of F_\infty(y) has a typo: it integrates with dy and uses h_\infty(y) inside a convolution in x. The intended formula is F_\infty(y)=\lambda\int (x-y)/|x-y|^3 |h_\infty(x)|^2 dx, as used in Corollary 4.4.
- [Section 1.2 and references] The paper should state explicitly which estimates from the unpublished companion [25] are actually used as inputs. The proof appears self-contained, but Section 1.2 describes [25] as the driving force; the reader needs a precise dependency list to audit the argument.
- [Title, Lemma 2.1, and Lemma 3.7] There are several typographical errors: 'SCA TTERING' in the title, 'gives' in Lemma 2.1(2), and 'femionic' in Lemma 3.7.
- [Section 3.3, Lemma 3.9] The chain-rule identity '\nabla_x\gamma(t,x-tv,v)=t^{-1}\nabla_v\gamma(t,x-tv,v)-t^{-1}\nabla_v[\gamma(t,x-tv,v)]' is confusing; please rewrite it with the intended variable change and clarify how integration by parts in v is used.
- [Section 4.2 and Proposition 4.3] In the final Hartree step, the rate t^{-1/80} is stated without specifying the choice of m \in (3/2,2) and \kappa in Proposition 4.3; please spell out the choice so the reader can verify the exponent.
Circularity Check
No significant circularity: the asymptotic profiles are constructed as limits rather than assumed, and the only self-referential element is a non-load-bearing citation to the authors' own forthcoming work [25]; the regularity gap under (1.7) is a correctness issue, not circularity.
full rationale
Inspection of the derivation chain shows that the asymptotic objects in Theorem 1.1 are not inputs but constructed limits. In Section 4.1, h(t)=g(t)e^{i\theta(t)} is introduced and Proposition 4.1 proves convergence to h_infty using only the bootstrap bounds on G^2u from Proposition 3.10; Corollary 4.4 then defines F_infty by (4.16) and proves t^2F(t,tv) -> F_infty. No parameter is fitted and no asymptotic profile is inserted by hand; the logarithmic corrections are derived from the time-integrated potential. Similarly, Lemma 4.5 constructs phi_infty as the limit of the effective potential t*phi^0(t,a), and Lemma 4.7 verifies phi_infty = V*rho_infty after convergence, so the coupling identities in (1.8) are proven consequences, not assumptions. The Section 3 bootstrap closes by strict improvement of the norms; assumptions (3.9) are auxiliary and are later consequences, so no circularity is present. The only self-referential element is the citation of the authors' own forthcoming work [25] in Section 1.2, described as 'the driving force of our paper.' However, no displayed estimate from [25] is invoked in Sections 3-4; the needed physical-space Hartree bounds are derived directly, so this citation is not load-bearing. Separately, the theorem as stated has a regularity gap: assumption (1.7) does not imply f0 in C^1, which Proposition 3.1 requires, and no approximation argument is supplied. That gap affects correctness of the stated regularity, not circularity of the derivation.
Assumptions & free parameters
assumptions (6)
- standard math Strichartz estimates, including endpoint Strichartz, for the three-dimensional Schrodinger propagator (Keel-Tao).
- standard math Lorentz-space Holder, convolution, Riesz potential, and Riesz transform estimates (O'Neil, Stein, Grafakos), such as (2.3), (2.6)-(2.8).
- standard math Liouville theorem: divergence-free Hamiltonian vector fields preserve phase-space measure, so transport along characteristics preserves L^p norms (Lemma 2.3).
- standard math The Coulomb kernel shell representation with a compactly supported chi and the identity C0/|z| = integral_0^infty chi(z/R) dR/R^2.
- domain assumption Initial data smallness and localization (1.7) with epsilon_0 below a universal threshold.
- domain assumption C^1 regularity of f0 for the local well-posedness theory in Proposition 3.1.
Cite this review
Pith. "Pith review of Long-Time Dynamics and Modified Scattering for the Coulomb Vlasov--Hartree System." pith.science (2026). https://pith.science/paper/4M64EEMJ
@misc{pith2026260804402,
author = {Pith},
title = {Pith review of: Long-Time Dynamics and Modified Scattering for the Coulomb Vlasov--Hartree System},
year = {2026},
howpublished = {\url{https://pith.science/paper/4M64EEMJ}},
note = {Machine review of arXiv:2608.04402}
}
read the original abstract
We study the three-dimensional Vlasov--Hartree system with Coulomb interaction, a classical--quantum model for Bose--Fermi mixtures. For sufficiently small, regular, and localized initial data, we prove global well-posedness and determine the precise long-time dynamics. The coupling remains visible at leading order: the Vlasov distribution scatters along logarithmically corrected free characteristics determined by the asymptotic bosonic profile, while the Hartree wave acquires a logarithmic phase correction generated by the asymptotic fermionic density. Our proof combines techniques from kinetic theory with vector-field methods for dispersive equations and also extends to shorter-range inverse-power interactions.
Reference graph
Works this paper leans on
-
[25]
B. Pausader and M. Xie,A Fourier-Free New Method for the 3D Hartree Equation: Full Space and half-space with Robin Condition,in preparation
-
[1]
C. Bardos and P. Degond,Global existence for the Vlasov–Poisson equation in 3 space variables with small initial data, Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire2(1985), no. 2, 101–118
work page 1985
-
[2]
L. Bigorgne and R. Velozo Ruiz,Late-time asymptotics of small data solutions for the Vlasov–Poisson system,Nonlinearity39(2026), no. 4, Paper No. 045007
work page 2026
-
[3]
E. C´ ardenas, J. K. Miller and N. Pavlovi´ c,On the effective dynamics of Bose–Fermi mixtures,Ars Inven. Anal. 2025, Paper No. 5, 54 pp
work page 2025
-
[4]
Global Existence and Time Decay for the Vlasov-Hartree System
K. Cavanagh,Global existence and time decay for the Vlasov–Hartree system,arXiv:2607.04444
-
[5]
S.-H. Choi and S. Kwon,Modified scattering for the Vlasov–Poisson system, Nonlinearity29(2016), 2755– 2774. MODIFIED SCATTERING FOR THE COULOMB VLASOV–HARTREE SYSTEM 29
work page 2016
- [6]
-
[7]
J. Ginibre and T. Ozawa,Long range scattering for non-linear Schr¨ odinger and Hartree equations in space dimensionn≥2,Comm. Math. Phys.151(1993), no. 3, 619–645
work page 1993
Show all 30 references
-
[8]
R. T. Glassey,The Cauchy problem in kinetic theory, SIAM, Philadelphia, PA, 1996; MR1379589
1996
-
[9]
Grafakos,Classical Fourier analysis, third edition, Graduate Texts in Mathematics, 249, Springer, New York, 2014; MR3243734
L. Grafakos,Classical Fourier analysis, third edition, Graduate Texts in Mathematics, 249, Springer, New York, 2014; MR3243734
2014
-
[10]
Hayashi and P
N. Hayashi and P. I. Naumkin,Asymptotics for large time of solutions to the nonlinear Schr¨ odinger and Hartree equations,Amer. J. Math.120(1998), no. 2, 369–389
1998
-
[11]
Hayashi, P
N. Hayashi, P. I. Naumkin and T. Ozawa,Scattering theory for the Hartree equation,SIAM J. Math. Anal. 29(1998), no. 5, 1256–1267
1998
-
[12]
Hong and S
Y. Hong and S. Pankavich,Modified scattering for the Vlasov–Riesz system with long-range interactions, arXiv:2604.04256
-
[13]
Huang and H
W. Huang and H. Kwon,Scattering of the Vlasov–Riesz system in the three dimensions,Ann. Henri Poincar´ e (2026)
2026
-
[14]
Huang and H
W. Huang and H. Kwon,Scattering map for the Vlasov–Poisson system with a repulsive harmonic potential, arXiv:2602.21344
-
[15]
Huang, B
W. Huang, B. Pausader and M. Suzuki,The Vlasov–Poisson system with a perfectly conducting wall: Convex domains,Comm. Math. Phys., to appear arXiv:2412.13434
-
[16]
Iacobelli, S
M. Iacobelli, S. Rossi, and K. Widmayer,On the stability of vacuum in the screened Vlasov–Poisson equation, Journal of the London Mathematical Society113(2026), no. 1, e70426
2026
-
[17]
A. D. Ionescu, B. Pausader, X. Wang and K. Widmayer,On the asymptotic behavior of solutions to the Vlasov–Poisson system,Int. Math. Res. Not. IMRN2022(2022), no. 12, 8865–8889
2022
-
[18]
Kato and F
J. Kato and F. Pusateri,A new proof of long-range scattering for critical nonlinear Schr¨ odinger equations, Differential Integral Equations24(2011), no. 9–10, 923–940
2011
-
[19]
Keel and T
M. Keel and T. C. Tao, Endpoint Strichartz estimates, Amer. J. Math.120(1998), no. 5, 955–980
1998
-
[20]
Kepka and K
B. Kepka and K. Widmayer,Modified scattering dynamics in the Vlasov–Poisson equation near an attractive point mass,arXiv:2511.04363
-
[21]
Lions and B
P.-L. Lions and B. Perthame,Propagation of moments and regularity for the 3-dimensional Vlasov–Poisson system, Invent. Math.105(1991), 415–430
1991
-
[22]
O’Neil, Convolution operators andL(p, q) spaces, Duke Math
R. O’Neil, Convolution operators andL(p, q) spaces, Duke Math. J.30(1963), 129–142; MR0146673
1963
-
[23]
Pausader and K
B. Pausader and K. Widmayer,Stability of a point charge for the Vlasov–Poisson system: The radial case, Comm. Math. Phys.385(2021), no. 3, 1741–1769
2021
-
[24]
Pausader, K
B. Pausader, K. Widmayer and J. Yang,Stability of a point charge for the repulsive Vlasov–Poisson system, J. Eur. Math. Soc.28(2026), no. 7, 2751–2848
2026
-
[26]
Peetre, Espaces d’interpolation et th´ eor` eme de Soboleff, Ann
J. Peetre, Espaces d’interpolation et th´ eor` eme de Soboleff, Ann. Inst. Fourier (Grenoble)16(1966), fasc. 1, 279–317; MR0221282
1966
-
[27]
Pfaffelmoser,Global classical solutions of the Vlasov–Poisson system in three dimensions for general initial data, J
K. Pfaffelmoser,Global classical solutions of the Vlasov–Poisson system in three dimensions for general initial data, J. Differential Equations95(1992), 281–303
1992
-
[28]
Schlue and M
V. Schlue and M. Taylor,Inverse modified scattering and polyhomogeneous expansions for the Vlasov– Poisson system, Nonlinearity38(2025), no. 9, 095019
2025
-
[29]
E. M. Stein,Singular integrals and differentiability properties of functions, Princeton Mathematical Series, No. 30, Princeton Univ. Press, Princeton, NJ, 1970
1970
-
[30]
Van Hoose,Modified scattering for the Hartree nonlinear Schr¨ odinger equation, arXiv preprint arXiv:2407.18411 (2024)
T. Van Hoose,Modified scattering for the Hartree nonlinear Schr¨ odinger equation, arXiv preprint arXiv:2407.18411 (2024). Brown University Email address:wenrui huang@brown.edu Yale University Email address:mengyi.xie@yale.edu
2024 arXiv
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.