Long-time dynamics of Vlasov--Hartree systems across the Coulomb threshold

Mengyi Xie; Wenrui Huang

Long-time dynamics of Vlasov--Hartree systems across the Coulomb threshold

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

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 →

arxiv 2608.04402 v2 pith:4M64EEMJ submitted 2026-08-05 math.AP

Wenrui Huang , Mengyi Xie This is my paper · ORCID

classification math.AP MSC 35Q8335Q5535B40

keywords Vlasov–HartreesystemCoulombinteractionmodifiedscatteringlong-timeasymptoticsBose–Fermimixtureglobalwell-posednessGalileanvectorfieldslogarithmicphasecorrection

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper studies the Vlasov–Hartree system in three dimensions with Coulomb interaction, a mean-field model of a Bose–Fermi mixture in which the bosonic density generates a force on the fermions and the fermionic density generates a potential for the bosons. For sufficiently small, regular, and localized initial data it proves global well-posedness and, more importantly, identifies the exact long-time asymptotics. The central claim is modified scattering with the coupling visible at leading order: the fermionic distribution converges along characteristics shifted by $\log(t) F_\infty(v)$, where $F_\infty$ is built from the limiting bosonic profile, while the bosonic wave acquires the phase $e^{-i \log(t)\phi_\infty(x/t)}$, with $\phi_\infty$ built from the limiting fermionic density. A sympathetic reading is that the two components never decouple: each species' memory of the other survives in the asymptotic data at order one.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Dear colleague,

The short version: this paper proves the first modified scattering theorem for the 3D Coulomb Vlasov–Hartree system, and the crossed log corrections (fermionic profile determines the bosonic phase, bosonic profile determines the fermionic drift) are genuinely new. The main bootstrap is detailed and internally coherent. But Theorem 1.1 as stated has a real regularity hole: assumption (1.7) does not imply f0 in C^1, while the local theory in Proposition 3.1 requires C^1 data, and the proof supplies no approximation argument. This is patchable in a paragraph (add f0 in C^1 to (1.7), or weaken the conclusion), and the asymptotic machinery itself only uses L^2∩L^∞ derivative bounds, so the main result is not in danger.

What is new: the coupled asymptotic pair (f∞, h∞) with mutually determined log corrections. This is not present in the Vlasov–Poisson modified scattering literature (Ionescu–Pausader–Wang–Widmayer et al.) or in Hartree long-range theory. The physical-space method from the companion [25] is adapted and mostly reproduced; the paper's own propositions (3.2–3.10, 4.1–4.7) carry the estimates, so [25] functions more as a methodological source than an unstated black box. Still, the authors should confirm which estimates are imported verbatim, since an unpublished dependency should be flagged.

Soft spots, in proportion:

1. The regularity gap (the stress-test has it right). A non-C^1 f0 with only weak L^2∩L^∞ derivative bounds will not become C^1; Vlasov transport preserves the singularity. The theorem's proof as written does not produce f in C^1. Fix: add f0 in C^1 to (1.7), or state the conclusion in W^{1,∞}∩W^{1,2} for the Vlasov component. Minor fix, but it must be made.

2. The dependence on [25]. The paper calls it "the driving force." A referee cannot fully audit that input if key estimates are only proved there. The present text seems to prove the needed Hartree asymptotics in Section 4, but the authors should state explicitly which lemmas are original here and which are quoted from [25]. A clarity issue, not a correctness crisis.

3. Small technical points: the abstract has a typo ("SCA TTERING"), and the transition from the local continuation criterion to the global logarithmic growth is terse but fine. I did not check every Lorentz-space constant; the estimates look consistent, and I found no contradiction in the main bootstrap.

Who it is for: anyone working on kinetic–dispersive coupling, modified scattering, or Vlasov–Poisson/Hartree asymptotics. It deserves a serious referee. I would send it out, with a request to patch the regularity statement and clarify the [25] dependence.

Recommendation: accept with revision, after the authors fix the f0 regularity hypothesis and state the provenance of the [25] estimates.

Referee Report

3 major / 5 minor

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

0 steps flagged · score 2.0 of 10

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 0 free parameters · 6 assumptions · 0 invented entities

No numbers are fitted to data and no new physical objects are introduced. The proof rests on standard analytic estimates plus explicit smallness and regularity assumptions. The only non-public input is the unpublished companion paper [25], cited as a source of the physical-space method but largely reproduced in the present proof.

assumptions (6)

standard math Strichartz estimates, including endpoint Strichartz, for the three-dimensional Schrodinger propagator (Keel-Tao).
Used in Section 2.3 and in the local well-posedness iteration to control u in mixed space-time norms.
standard math Lorentz-space Holder, convolution, Riesz potential, and Riesz transform estimates (O'Neil, Stein, Grafakos), such as (2.3), (2.6)-(2.8).
These estimates are the main analytic tools for bounding Coulomb potentials and fields in Sections 2 and 3.
standard math Liouville theorem: divergence-free Hamiltonian vector fields preserve phase-space measure, so transport along characteristics preserves L^p norms (Lemma 2.3).
Used throughout the Vlasov component analysis, including bootstrap bounds for weighted moments and derivatives of gamma.
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.
Used in Section 3.2 to define the effective fields phi_R, E_R, phi_r^0, and E_r^0 and to derive sharp decay.
domain assumption Initial data smallness and localization (1.7) with epsilon_0 below a universal threshold.
The theorem is conditional on this small-data regime; no explicit admissible value of epsilon_0 is given.
domain assumption C^1 regularity of f0 for the local well-posedness theory in Proposition 3.1.
Proposition 3.1 explicitly assumes f0 in C^1, while Theorem 1.1 states C^1 solutions under assumption (1.7) only. This mismatch is a gap in the statement.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Long-time dynamics of Vlasov--Hartree systems across the Coulomb threshold." pith.science (2026). https://pith.science/paper/4M64EEMJ

@misc{pith2026260804402,
  author       = {Pith},
  title        = {Pith review of: Long-time dynamics of Vlasov--Hartree systems across the Coulomb threshold},
  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 and inverse-power interactions. For small, regular, localized data, we prove global well-posedness and determine long-time dynamics. The main result concerns the Coulomb case, where the coupling persists at leading order through reciprocal asymptotic corrections: the fermionic distribution scatters along logarithmically corrected free characteristics set by the asymptotic bosonic profile, while the bosonic wave gains a logarithmic phase set by the asymptotic fermionic density. We also treat longer- and shorter-range interactions. Our proof combines Hamiltonian methods for the Vlasov equation with purely physical-space vector-field methods for dispersive equations. For stronger long-range interactions, we introduce a symplectic transformation and a regularized effective-phase cancellation, which remove leading nonintegrable interactions while preserving Hamiltonian structure and avoiding derivative loss.

Discussion (0). Continue with ORCID to comment.

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.