{"id":"af46602c-1caa-44e3-856f-4c828242934c","arxiv_id":"2506.19784","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Anisotropic viscous-resistive compressible MHD weak solutions converge, up to subsequences, to weak solutions of reduced MHD in T3 and R3, with identified initial data and the constraint bϱ+B=0.","lead":"Global weak solutions of strongly anisotropic, strongly magnetized compressible MHD are proven to converge to the reduced MHD model that fusion and astrophysics codes use. This is the first rigorous justification of the limit in the weak solution setting, in both periodic domains and the whole space.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 1's uniform L^kappa bound on varrho^varepsilon fails for 3/2<gamma<2: the energy scaling in (33) only gives an epsilon^{gamma-2} blow-up, so the weak compactness asserted in Theorems 1-2 is not established for this range.","rationale":"The reader identified the scaling in Section 3.1 as the weakest assumption and accepted the paper with moderate confidence. My stress-test found a different, more serious problem: the uniform bound on the rescaled density varrho^varepsilon, which is the basis for all subsequent weak compactness, is false for the parameter range 3/2<gamma<2 under the paper's own energy estimates. The energy in (33) controls the relative entropy with the prefactor epsilon^{-2(gamma-1)}. For gamma<2, Lemma 11's convexity inequalities imply only the diverging bounds recorded above. An explicit admissible initial datum (rho_0=1+epsilon^{gamma-1}f, v_0=0, B_0=0) satisfies assumption (27) yet makes varrho_0 unbounded in L^gamma. Since Theorem 1 asserts weak-* convergence of varrho^varepsilon in L^infty_loc(R_+;L^kappa) with kappa=min(2,gamma), the statement cannot hold for this data, independent of the evolution. The same defect propagates into the whole-space result through Lemma 6. This is not a matter of a missing technical detail in a single step; it invalidates the main theorems for a nontrivial interval of gamma. I therefore recommend rejecting the paper in its current form. A revised version restricting to gamma>=2, or adding conditions on the initial data that restore boundedness of varrho^varepsilon, could be reconsidered, but the central claim as stated is not supported.","tokens_in":51570,"tokens_out":47649,"duration_ms":483821,"concrete_test":"Recompute the bound in Lemma 1 for gamma=1.8 on T^3. Choose a smooth mean-zero f, set rho^varepsilon_0=1+epsilon^{0.8}f, v^varepsilon_0=0, B^varepsilon_0=0. Verify that the relative-entropy term in (27) is a gamma(gamma-1)/2 ||f||_{L^2}^2 (bounded), while ||(rho^varepsilon_0-rho^varepsilon_0)/epsilon||_{L^gamma}=epsilon^{-0.16}||f||_{L^gamma}, which diverges. Then use Lemma 11 and inequality (30)-(33) to derive the sharp bound on the large-density part: it is epsilon^{gamma-2}, not epsilon^{(2/gamma)-1}. If the derived bound diverges, Lemma 1's estimate (42) is wrong and the claimed weak compactness of varrho^varepsilon fails for this admissible initial data.","verdict_should_be":"REJECT","load_bearing_attack":"The proof of Lemma 1 claims, for 3/2<gamma<2, that varrho^varepsilon=(rho^varepsilon-rho^varepsilon)/varepsilon is uniformly bounded in L^infty_loc(R_+;L^kappa) with kappa=gamma (see (42)). This does not follow from the stated energy. Energy inequality (30)-(33) controls Pi_2(rho^varepsilon) = a/epsilon^{2(gamma-1)}[(rho^varepsilon)^gamma - gamma rho^varepsilon (rho^varepsilon)^{gamma-1} + (gamma-1)(rho^varepsilon)^gamma]. For gamma<2, Lemma 11 gives a lower bound A >= nu_3 |rho^varepsilon-rho^varepsilon|^gamma on large density regions, so energy only implies ||varrho^varepsilon 1_{rho^varepsilon>=R}||_{L^gamma}^gamma <= C epsilon^{gamma-2}, which diverges as epsilon->0. Similarly, on bounded-density regions the quadratic lower bound yields ||varrho^varepsilon 1_{rho^varepsilon<R}||_{L^2}^2 <= C epsilon^{2gamma-4}, which also diverges. The estimate (42) appears to assume an epsilon^2 energy rate that (33) does not provide for gamma<2. This is not a cosmetic issue: take gamma in (3/2,2) and initial data rho^varepsilon_0=1+epsilon^{gamma-1}f with f smooth, mean zero, v_0=0, B_0=0. Then (27) holds with relative entropy of order O(1), while varrho^varepsilon_0=epsilon^{gamma-2}f has L^gamma norm epsilon^{gamma-2}||f||_{L^gamma}->infty. No subsequence can converge weak-* in L^infty_loc(R_+;L^gamma), contradicting the stated conclusion of Theorem 1 at t=0. The same problem enters Theorem 2 via Lemma 6. Since Theorems 1 and 2 are the central claim and they are stated for all gamma>3/2, the gamma in (3/2,2) range is not covered by the proof as written.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies, for the anisotropic compressible MHD system (3)-(4) with strong vertical magnetic field, the singular limit of global weak solutions as the aspect-ratio and time-scale parameters go to zero. The authors claim that, up to subsequences, these solutions converge to weak solutions of the reduced MHD system (12)-(17), with perpendicular initial data given by transverse Leray projections and with the algebraic relation bϱ+B=0 for positive times. The proofs are carried out separately in the periodic domain T^3 and in the whole space R^3, using energy estimates, fast magnetosonic wave filtering by a unitary group, Simon compactness, and the Lions-Masmoudi compactness lemma. The paper is written for adiabatic exponents γ>3/2 and includes both prepared and unprepared initial data, with explicit discussion of the time boundary layer.","tokens_in":51938,"tokens_out":9364,"duration_ms":99922,"significance":"If the main theorem is correct, the paper is a substantial contribution: it gives the first rigorous weak-solution level justification of the RMHD model from compressible MHD, in both periodic and whole-space settings, with general unprepared data and with the correct projection operator appearing at t=0. The proof is detailed, carefully structured, and built on appropriate and well-tested tools (Hu-Wang existence theory, Lions-Masmoudi compactness, Simon's lemma, and an explicit unitary group for transverse fast magnetosonic waves). The authors are also transparent about limitations, such as the lack of anisotropic Strichartz estimates. However, as detailed below, the stated range γ>3/2 is not supported by the argument; the density-corrector estimate fails for 3/2<γ<2, so the theorems as stated are false in that range. The paper would be sound if the statement were restricted to γ≥2, unless a new density estimate is found.","major_comments":[{"comment":"The claimed uniform L^κ bound for ϱε fails for 3/2<γ<2. From the energy inequality (30)-(33) and Lemma 11, the correct consequences are ‖ϱε 1_{ρε≤R}‖_{L∞L^2}^2 ≲ ε^{2γ-4} and ‖ϱε 1_{ρε>R}‖_{L∞L^γ}^γ ≲ ε^{γ-2}; both right-hand sides diverge as ε→0 for γ<2, whereas (42) asserts boundedness and uses the exponent ε^{2/γ-1}. This estimate is the only mechanism producing the weak limit ϱ in L∞L^κ and consequently the relation bϱ+B=0 (Lemma 4, point 3, and Section 4.7). The gap is load-bearing: with initial data ρ0ε=1+ε^{γ-1}f (f smooth with zero mean), v0ε=0, B0ε=0, assumption (27) is satisfied with O(1) energy, but ϱ0ε=ε^{γ-2}f diverges in L^γ for γ<2, so the weak-* convergence claimed in Theorem 1 cannot hold at t=0.","section":"§4.1, Lemma 1, estimate (42)"},{"comment":"The same issue appears in the whole-space proof. The bound ‖ρε−1‖_{L∞L^γ_loc}≤Cε stated for 1<γ<2 is stronger than what (68) together with Lemma 11 yields; the correct local bound is O(ε^{2(γ−1)/γ}), which tends to zero but is weaker than ε. Consequently the uniform bound for ϱε in L^κ_loc in the fifth assertion of Lemma 6 also fails for γ<2, and Theorem 2 inherits exactly the same obstruction as Theorem 1. A different density estimate, or an explicit restriction to γ≥2, is needed before the passage to the limit is justified.","section":"§5.1, Lemma 6, third and fifth assertions"},{"comment":"Because the density-corrector compactness is used in several load-bearing places (the filtered profile in Lemma 4, the pressure terms πε2, the identification of initial data in (17), and the derivation of the parallel equations (16)), the gap is not local. The counterexample at t=0 described above shows that the theorems as stated are false for 3/2<γ<2, not merely unproved. The authors should either prove a genuinely new uniform estimate for ϱε in this range or explicitly restrict the main theorems to γ≥2, with all subsequent statements adjusted accordingly.","section":"Theorems 1 and 2"}],"minor_comments":[{"comment":"The word 'ANISTROPIC' in the title should be 'ANISOTROPIC'.","section":"Title"},{"comment":"In the whole-space section, the reference 'from (27)' should be 'from (36)', and the regularity lines in (37) and Theorem 2 contain T^3 where R^3 is intended (e.g., H^1(T^3) and L^2(T^3)).","section":"§2.4 and Theorem 2"},{"comment":"The notation H^1_0(T^3) is used for the mean-zero subspace of \\.H^1(T^3); since H^1_0 usually denotes zero trace, a different symbol would avoid confusion.","section":"§4.3, proof of Lemma 3"},{"comment":"The symbol κ is overloaded: it is defined both as κ=min{2,γ} and later as κ=max{1/2,3/(2γ)}; please use a distinct letter for the second quantity.","section":"§5.3, Lemma 9"},{"comment":"There are several typos (e.g., 'seventhies', 'ubiquitus', 'regim', 'anistotropic') and inconsistent uses of the pair (ρ,ρ̄) in the initial-data discussion; a careful proofreading pass is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The γ<2 issue is a genuine counterexample at t=0, not a missing line in an otherwise correct proof. The theorem as stated is false for 3/2<γ<2, which includes the physically common value γ=5/3. If the authors can prove a new estimate, the result would be very strong; if not, a restriction to γ≥2 would make the proof sound but would remove a physically relevant range. I would ask the authors to address this head-on before publication, and would not accept the paper with the current statement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know first: the stress-test note is right. The paper's central claim, weak compactness of the density corrector varrho^epsilon in L^kappa for all gamma>3/2, does not hold for 3/2<gamma<2. Lemma 1 asserts a bound of order epsilon^{2/gamma-1} for the L^gamma norm of varrho^epsilon on large-density regions, but the energy only gives epsilon^{2(gamma-1)} for the L^gamma distance between rho^epsilon and its mean, which translates to epsilon^{1-2/gamma} for varrho^epsilon. That exponent is negative in the range 3/2<gamma<2, so the bound blows up. Concretely, take rho^epsilon_0 = 1 + epsilon^{gamma-1} f with f nonzero, mean zero, and smooth. The energy (27) is satisfied with O(1) relative entropy, while varrho^epsilon_0 = epsilon^{gamma-2} f diverges in L^gamma. So the weak convergence of varrho^epsilon at t=0 fails, and Theorem 1, as stated, is false. Theorem 2 inherits the same problem through Lemma 6. This is not a cosmetic issue: the compactness of varrho^epsilon is used in the filtering argument and in identifying the limit system.\n\nThe paper nevertheless has real substance. It is the first attempt I know to pass from global weak solutions of anisotropic compressible MHD to RMHD, in both T^3 and R^3, with unprepared data. The unitary-group filtering of perpendicular fast magnetosonic waves, the cancellations that extract the parallel dynamics, and the careful treatment of the whole-space case show serious work. The authors also flag the limitations honestly: Q_perp v^epsilon_perp only converges weakly, and |B^epsilon|^2 does not converge in general.\n\nOther soft spots are minor by comparison: the proof of Lemma 10 is very intricate, and Appendix Lemmas 11 and 12 leave nontrivial details to the reader. There are also typos.\n\nThe load-bearing flaw is the gamma range. If the theorems are restricted to gamma >= 2, the proof appears to go through, and that remains a substantial result. For gamma in (3/2,2), either extra assumptions on initial data or a weaker notion of convergence for varrho^epsilon is needed.\n\nMy recommendation: send this to a serious referee, but alert them to the exponent error in Lemma 1. The paper should not be accepted as is; it needs major revision or a narrowed statement. I would not cite the current version until the gamma question is resolved.","headline":"A serious paper with a real gap in the gamma range: the density corrector bound fails for 3/2<gamma<2, so the main theorems are not established as stated.","tokens_in":52559,"tokens_out":7208,"would_cite":false,"duration_ms":69969,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B25","35Q35","35Q85","76W05","35D30"],"pacs":[],"model":"deepseek-v4-flash","headline":"MHD weak solutions converge to RMHD under strong anisotropy","keywords":["reduced magnetohydrodynamics","singular limits","weak solutions","anisotropic plasmas","compressible MHD","low Mach number limit","fast magnetosonic waves","large magnetic field"],"falsifier":"A decisive numerical or analytical test takes unprepared data in $\\mathbb{T}^3$ with $\\nabla_\\perp\\cdot v_{0\\perp}\\neq 0$ and $B_{0\\parallel}+b\\varrho_0\\neq 0$, solves the linearized version of (3) for small $\\varepsilon$, and checks that after the fast magnetosonic transient the perpendicular velocity has projected to $P_\\perp v_{0\\perp}$ and the parallel magnetic field has relaxed to $(B_{0\\parallel}-\\varrho_0)/c$; persistent $O(1)$ deviations from those projected initial values, or any surviving component of $B_\\parallel + b\\varrho$, would contradict the theorem.","tokens_in":51287,"feed_emoji":"🧲","tokens_out":7933,"duration_ms":76428,"temperature":0.7,"pith_summary":"The paper proves that global weak solutions of viscous and resistive barotropic compressible magnetohydrodynamics pass, in the singular limit of strong magnetization and spatial anisotropy, to weak solutions of the reduced magnetohydrodynamic (RMHD) equations used in fusion and astrophysics. The convergence is established in the periodic box $\\mathbb{T}^3$ and in the whole space $\\mathbb{R}^3$, up to subsequences, with the density converging strongly to $1$ and the perpendicular magnetic field, the solenoidal perpendicular velocity, and the parallel fields converging in the stated strong or weak topologies. The limit is only perpendicularly incompressible: $\\nabla_\\perp \\cdot v_\\perp = 0$ and $\\nabla_\\perp \\cdot B_\\perp = 0$, while the parallel direction keeps compressible information through the algebraic relation $b\\varrho + B_\\parallel = 0$. The proof works for unprepared initial data because a time boundary layer at $t=0$ projects the initial perpendicular velocity and modifies the initial parallel magnetic field.","feed_headline":"MHD weak solutions converge to RMHD under strong anisotropy","feed_subtitle":"Rigorous limit on the torus and in all of space: plasma becomes perpendicularly incompressible but keeps parallel compressibility.","key_machinery":"The load-bearing object is the fast-oscillatory unitary group $S(\\tau) = \\exp(\\tau L)$ generated by the skew-adjoint operator $L$ acting on the pair $U^\\varepsilon = (b_\\varepsilon \\varrho^\\varepsilon + B_\\parallel^\\varepsilon,\\, Q_\\perp(\\rho^\\varepsilon v_\\perp^\\varepsilon))$, where $P_\\perp$ and $Q_\\perp$ are the perpendicular divergence-free and gradient projections. The operator $L$ encodes the transverse fast magnetosonic waves, whose perpendicular speeds are $\\pm\\sqrt{b+1}$; wrapping $U^\\varepsilon$ in $S(-t/\\varepsilon)$ filters those waves and leaves a compact filtered profile. In the parallel direction the proof uses an algebraic cancellation: the combination $B_\\parallel^\\varepsilon - \\varrho^\\varepsilon/\\rho^\\varepsilon$ evolves with no $\\varepsilon^{-1}$ singularity, which together with a product compactness criterion yields the linear parallel system and the relation $b\\varrho + B_\\parallel = 0$.","core_discovery":"The central claim is a two-part convergence statement. First, in both $\\mathbb{T}^3$ and $\\mathbb{R}^3$ and for adiabatic exponent $\\gamma > 3/2$, any sequence of global weak solutions to the penalized anisotropic MHD system (3)-(4) converges, up to a subsequence, to $(1,\\varrho,v,B)$ where the perpendicular pair $(v_\\perp,B_\\perp)$ solves the closed nonlinear RMHD system (12)-(13) with initial data $(P_\\perp B_{0\\perp}, P_\\perp v_{0\\perp})$, the solenoidal part of the original data in the perpendicular plane. Second, the parallel components $(v_\\parallel, B_\\parallel)$ solve the linear transport system (16) with $c = 1 + 1/b$ and initial data $B_\\parallel(0) = (B_{0\\parallel} - \\varrho_0)/c$, and the first-order density is slaved to the parallel magnetic field by $b\\varrho + B_\\parallel = 0$. For unprepared data the same limit emerges because a time boundary layer at $t=0$ removes the irrotational part $Q_\\perp v_{0\\perp}$ and redefines $B_\\parallel(0)$.","pith_inferences":["A natural test is to vary the scaling ratios, e.g. $\\mu_\\perp \\sim \\varepsilon^2$ or $\\eta_\\perp \\sim \\varepsilon^2$; the paper's limit equations depend on the chosen ratios, so one would expect different reduced models or failure of RMHD beyond the threshold.","The proof's separation of perpendicular and parallel mechanisms suggests the same strategy can be carried to extended MHD with extra Hall or electron-pressure terms, where the fast transverse waves would have modified speeds.","A Fourier computation in the periodic case shows that resonant modes with equal perpendicular wavenumbers survive the filtering; this indicates that the perpendicular filtered component feeds the pressure term, which could be probed numerically by computing $\\pi_3$ from two-point correlations.","An anisotropic Strichartz estimate for the transverse wave equation, if available, would likely upgrade the weak convergence of $Q_\\perp v_\\perp$ to strong convergence in the whole space, settling an open question noted in the paper."],"forward_implications":["The RMHD system is the true weak limit of dissipative compressible MHD in the strongly anisotropic, strongly magnetized regime; it is not only a formal asymptotic model.","In this limit the plasma flow is incompressible in the plane perpendicular to the background field, while the parallel direction retains compressibility through the relation $b\\varrho + B_\\parallel = 0$.","Unprepared initial data are admissible: fast magnetosonic oscillations are filtered by a time boundary layer that effectively replaces $v_{0\\perp}$ by $P_\\perp v_{0\\perp}$ and $B_{0\\parallel}$ by $(B_{0\\parallel}-\\varrho_0)/c$.","The same passage to the limit gives an alternative construction of global weak solutions of the RMHD system by letting $\\varepsilon \\to 0$.","In the whole space the proof manages without isotropic Strichartz estimates, using a truncation-and-mollification argument adapted to the anisotropic wave equation, which is the reason the result covers $\\mathbb{R}^3$ as well."],"supporting_citations":[{"why":"supplies the global existence of weak solutions to the three-dimensional compressible MHD system (1), i.e. the solutions whose limit is studied.","marker":"[23]"},{"why":"provides the filtered-profile and unitary-group method for the low-Mach-number limit that the paper adapts to perpendicular magnetosonic waves.","marker":"[33]"},{"why":"supplies the product compactness criterion and weak-solution energy framework used to pass to the limit in the nonlinear and parallel terms.","marker":"[32]"},{"why":"introduces the slow/fast classification and the unitary group for fast singular limits used to filter fast oscillations.","marker":"[42]"},{"why":"gives the classical existence theory for weak solutions of the incompressible MHD system at the base of the RMHD limit system.","marker":"[41]"},{"why":"supplies the classical existence of weak solutions to the limiting incompressible MHD equations used in the RMHD well-posedness discussion.","marker":"[14]"}],"fun_headline_variants":["MHD meets RMHD: rigorous anisotropic limit","Anisotropy drives MHD to RMHD in strong field","Plasma limit: only perpendicular incompressible","MHD weak solutions converge to RMHD system"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on the exact anisotropic ordering that fixes the perpendicular and parallel length scales in ratio $\\varepsilon$ and scales the viscosities and resistivities as $\\mu_\\perp = \\varepsilon \\mu_\\perp^\\varepsilon$, $\\mu_\\parallel = \\mu_\\parallel^\\varepsilon/\\varepsilon$, and similarly for $\\lambda$ and $\\eta$, with $\\varepsilon \\to 0$, so if a plasma obeys different anisotropy scalings the RMHD limit, including which Laplacians survive, can change or fail.","fun_headline_variants_meta":{"raw":{"variants":["MHD meets RMHD: rigorous anisotropic limit","Anisotropy drives MHD to RMHD in strong field","Plasma limit: only perpendicular incompressible","MHD weak solutions converge to RMHD system"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000584,"raw_usage":{"total_tokens":2796,"prompt_tokens":1047,"completion_tokens":1749,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":1687}},"tokens_in":663,"tokens_out":1749,"duration_ms":12860,"temperature":1.0,"reasoning_tokens":1687,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:25:15.247444+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive numerical or analytical test takes unprepared data in $\\mathbb{T}^3$ with $\\nabla_\\perp\\cdot v_{0\\perp}\\neq 0$ and $B_{0\\parallel}+b\\varrho_0\\neq 0$, solves the linearized version of (3) for small $\\varepsilon$, and checks that after the fast magnetosonic transient the perpendicular velocity has projected to $P_\\perp v_{0\\perp}$ and the parallel magnetic field has relaxed to $(B_{0\\parallel}-\\varrho_0)/c$; persistent $O(1)$ deviations from those projected initial values, or any surviving component of $B_\\parallel + b\\varrho$, would contradict the theorem.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the global existence of weak solutions to the three-dimensional compressible MHD system (1), i.e. the solutions whose limit is studied."},{"cited_title":"Lions, N","cited_arxiv_id":null,"evidence_quote":"provides the filtered-profile and unitary-group method for the low-Mach-number limit that the paper adapts to perpendicular magnetosonic waves."},{"cited_title":"Lions, Mathematical topics in fluid dynamics","cited_arxiv_id":null,"evidence_quote":"supplies the product compactness criterion and weak-solution energy framework used to pass to the limit in the nonlinear and parallel terms."},{"cited_title":"Schochet, Fast singular limits of hyperbolic PDEs, J","cited_arxiv_id":null,"evidence_quote":"introduces the slow/fast classification and the unitary group for fast singular limits used to filter fast oscillations."},{"cited_title":"Sermange, R","cited_arxiv_id":null,"evidence_quote":"gives the classical existence theory for weak solutions of the incompressible MHD system at the base of the RMHD limit system."},{"cited_title":"Duvaut, J.-L","cited_arxiv_id":null,"evidence_quote":"supplies the classical existence of weak solutions to the limiting incompressible MHD equations used in the RMHD well-posedness discussion."}],"review_version":2}