{"id":"6fbb9878-5a13-4650-b70b-d9ff6dd346cf","arxiv_id":"2508.19386","paper_version":7,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"An axisymmetric HDG solver for stellar oscillation equations with gravity and differential rotation is validated by matching observed solar mode-splitting coefficients.","lead":"This paper implements and validates a computer solver for solar oscillations that includes both gravity and differential rotation, using the HDG numerical method. It reproduces observed rotation-induced shifts in acoustic mode frequencies, a capability needed for local helioseismology to probe the Sun's interior.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the central a-coefficient validation is strongly supported, and the flagged f_ghost inconsistency is a negligible second-order effect.","rationale":"The paper's central claim is that HDG reproduces observed rotation-induced mode splitting. This is tested two ways: (1) against an analytic formula for a simplified rotation profile, where a1/a3/a5 agree to within 0.03 nHz; (2) against global-mode observations for the solar profile, agreeing within 0.01–0.03 nHz, well inside observational errors. The analytic test is independent of the observations and exercises the same rotation terms, so it provides strong evidence that the rotation operator is implemented correctly. Combined with the no-rotation power spectrum matching GYRE, the methodology is well supported. The one disclosed modeling approximation (Remark 7) is that the spherical Model-S background plus rotation does not satisfy rotational hydrostatics, leaving a ghost force f_ghost=ρ0Ω×Ω×x that is not implemented. However, for the validated high-degree p-mode (ℓ=85, n=7, ~3 mHz), this is a centrifugal correction of order (Ω/ω)^2 ~ 2×10^-8 relative to the frequency, implying a frequency bias ~10^-4 nHz—orders of magnitude below the quoted agreement and the observational uncertainties. The 3% deviation near the center is irrelevant because the mode sensitivity is near the surface. I also considered whether the a-coefficient match is circular because the rotation profile is derived from the same observed splittings. It is a consistency check rather than a new prediction, but it still validates the forward solver, and the independent analytic profile test eliminates the risk of tuned artifacts. Overall, no load-bearing concern was identified; the reader's flagged f_ghost issue is a legitimate but quantitatively negligible approximation for the claimed result.","tokens_in":947,"tokens_out":866,"duration_ms":114407,"concrete_test":"Recompute Table 1 for the solar rotation profile with a self-consistently centrifugally distorted background (or with the omitted f_ghost term included) and compare a1, a3, a5 for ℓ=85, n=7. If any a-coefficient shifts by more than 0.05 nHz, the f_ghost omission would be load-bearing; otherwise the current verdict stands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I find no load-bearing objection to the central claim. The strongest validation—a-coefficients matching observations to within 0.01–0.03 nHz for ℓ=85, n=7—is corroborated by an independent analytic check for the simplified rotation profile and by GYRE agreement in the no-rotation case. The one flagged inconsistency (Remark 7: using the non-rotating Model-S background with rotation implies f_ghost=ρ0Ω×Ω×x≠0, which is not implemented) is a real modeling approximation, but it is second-order in Ω/ω. For the validated high-degree mode at ~3 mHz, the centrifugal correction changes eigenfrequencies by about (Ω/ω)^2 times the frequency, i.e. ~10^-4 nHz—far below the quoted agreement and the 0.05–0.10 nHz observational errors. The 3% deviation near the center is irrelevant because high-degree p-modes do not sample the deep interior. The paper discloses this limitation and the stabilization choices, and the numerical method is independently supported by the no-rotation GYRE comparison and by the analytic simplified-profile a-coefficients. No assumption was found that, if wrong, would overturn the accepted verdict.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Hybridizable Discontinuous Galerkin (HDG) solver for the time-harmonic equations of linear stellar oscillations in axisymmetric backgrounds with differential rotation and gravity. The working equations are a first-order system in Lagrangian displacement and Eulerian pressure perturbation, equivalent to Galbrun-type formulations but with lower regularity requirements. The authors derive the azimuthal-mode decomposition, give a full HDG discretization with a stabilization strategy, and study the role of the buoyancy frequency in creating elliptic/hyperbolic regions. They then synthesize helioseismic observables (power spectra and cross-covariances) from the computed Green's kernels. Without rotation, the synthetic power spectrum matches GYRE eigenfrequencies; with a simplified rotation profile, fitted a-coefficients agree with analytic values; with the observed solar rotation profile of Larson & Schou, the simulated a1, a3, a5 agree with observed values to within roughly 0.01-0.03 nHz. The central claim is that the solver can accurately compute Green's kernels for realistic solar backgrounds including gravity and differential rotation.","tokens_in":60084,"tokens_out":19999,"duration_ms":222412,"significance":"If the results hold, the paper provides a validated, open-source 2.5D forward solver for local helioseismology with gravity and rotation, a substantial step beyond earlier no-gravity or spherically symmetric tools. The strongest points are the independent no-rotation validation against GYRE, the analytic simplified-rotation check, and the excellent agreement with observed rotation-induced mode splitting. The detailed empirical study of HDG stabilization in the presence of a non-zero buoyancy frequency is also valuable. Main caveats are the non-equilibrium background approximation (Remark 7), the partly circular nature of the solar-profile comparison, and the fact that the rotation validation is performed for a single (ℓ,n) mode; none of these appears to invalidate the quantitative agreement, but they should be addressed explicitly.","major_comments":[{"comment":"The rotating validation uses the non-rotating Model-S background together with Ω≠0, so f_ghost=ρ0Ω×Ω×x is nonzero and not implemented. This is disclosed, and it is not fatal: the matched mode is ℓ=85, n=7 at ~3 mHz, a high-degree p-mode concentrated near the surface, while the hydrostatic deviation is largest in the deep interior; the fitted a1,a3,a5 are first order in Ω, whereas f_ghost is O(Ω^2) and affects even coefficients. Still, the paper should state this scaling and give a quantitative bound on the effect. In addition, because profile [39] is inferred from those same observed a-coefficients, the solar-profile comparison is partly a consistency check; the independent analytic simplified-profile validation in Table 1 should be presented as the primary quantitative evidence.","section":"§3.1 (Remark 7), §7.3.2, Table 1"}],"minor_comments":[{"comment":"The sentence 'ellipticity is lost when N^2>0 and for frequencies above N' appears to reverse the definitions in (3.49) and §3.5: the hyperbolic region is ω^2<N^2. Please correct the wording.","section":"§1"},{"comment":"The cross-covariance definition omits a complex conjugate. The later derivation using the addition theorem implicitly assumes conjugation of spherical harmonics; without a stated convention, (7.22) is ambiguous. Please clarify the Fourier convention and the role of conjugation.","section":"§7.1.2, Eqs. (7.17)-(7.22)"},{"comment":"The stabilization scaling α_tune is a free parameter with empirically chosen values (α_tune=±10^6 in (6.5)). A brief statement on the sensitivity of the §7.3 a-coefficient results to this parameter, or a note that the source geometry used for the observables is stable across τa/τb, would improve reproducibility.","section":"§5.5, §6.3"},{"comment":"Typo: 'hystrostatic' should be 'hydrostatic'. Also, the 'at most 3%' deviation should be defined more precisely (which norm? at which radius?), since it is an important caveat for the rotating background.","section":"Remark 7"},{"comment":"State explicitly that the observed a-coefficients quoted from [39] are for ℓ=85, n=7, and note any mode-selection details; currently the reader must infer this from Table 1 and the surrounding text.","section":"§7.3.1, Eq. (7.44)"},{"comment":"Please specify whether γ_att=10 μHz and 2 μHz are angular frequencies or ordinary frequencies, since the equations use angular frequency ω. This affects how readers interpret the attenuation relative to 3 mHz and 0.2 mHz.","section":"§6.1"},{"comment":"The caption and legend should define 'Approx', 'Fit', 'Max' and the green/gray curves more explicitly; the current caption relies heavily on the body text.","section":"Figure 16"}],"recommendation":"minor_revision","confidential_remarks":"The paper is a solid numerical-methods contribution with unusually strong validation. The main caveats—the non-implemented f_ghost term and the partly circular solar-profile comparison—are disclosed and are, in my reading, second-order for the claimed 0.03 nHz agreement, but they deserve explicit quantitative comment. The one-mode rotation validation is narrow but acceptable given the independent analytic simplified-rotation check. I would not require new simulations, though a numerical or analytic estimate of the f_ghost effect would strengthen the central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know about this paper is that it actually delivers what the abstract promises: an HDG solver for the axisymmetric oscillation equations with both gravity and differential rotation, with validation that goes beyond self-consistency. The synthetic power spectrum matches GYRE eigenfrequencies for the no-rotation case, and the simulated a-coefficients for rotation match both the analytic simplified-profile values and the observed solar values within 0.01–0.03 nHz. For a numerical methods paper in local helioseismology, that is strong evidence the physics and the discretization are both right.\n\nWhat's genuinely new: previous axisymmetric work in [24] ignored gravity. Here the inclusion of gravity changes the character of the operator via the buoyancy frequency N, producing elliptic/hyperbolic regions that matter for stabilization. The authors analyze that, design the HDG stabilization accordingly, and show empirically that a naive stabilization choice gives artifacts. The derivations are detailed and careful, and Remark 7 openly discloses the main modeling approximation: the non-rotating Model-S background is used together with a rotation profile, so the implied ghost force f_ghost = rho0 Omega x Omega x x is not balanced. That's a real inconsistency, but for the high-degree p-modes used in the validation it is second-order — the centrifugal correction to eigenfrequencies is about (Omega/omega)^2, far below the quoted agreement. The 3% deviation near the center is irrelevant because those modes don't sample the deep interior. So this isn't a load-bearing flaw.\n\nThe soft spots are minor. The stabilization parameter alpha_tune is empirically tuned, so the method needs case-by-case care, but the paper says this and studies it. The a-coefficients are extracted via Lorentzian fitting, which is standard, and the observational comparison uses a rotation profile inferred from mode splitting, so there is some circularity in the final validation; the analytic simplified-profile check breaks that circularity. Citation patterns are mostly the group's own previous work, which is appropriate given the method lineage.\n\nWho should read this: helioseismologists who need Green's kernels for inversions, and numerical analysts interested in HDG for mixed-type operators. It deserves a serious referee. I'd send it to review.","headline":"A careful, well-validated forward-modeling paper for axisymmetric solar oscillations; the f_ghost approximation is real but negligible for the modes tested, and the paper deserves serious refereeing.","tokens_in":60535,"tokens_out":1795,"would_cite":true,"duration_ms":20266,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N30","85-08"],"pacs":["96.60.Ly","02.70.Dh"],"model":"deepseek-v4-flash","headline":"The paper claims its HDG solver for axisymmetric solar oscillations—with gravity, attenuation, and differential rotation—produces power spectra whose rotation-induced splittings match observed a-coefficients to within 0.03 nHz.","keywords":["helioseismology","stellar oscillations","differential rotation","Hybridizable Discontinuous Galerkin","mode splitting","buoyancy frequency","Green's functions","power spectra"],"falsifier":"Run the same pipeline for the lowest harmonic degrees (ℓ ≲ 20), whose acoustic paths reach the deep interior where the neglected ghost force peaks near 3% of the hydrostatic balance; a systematic departure of the fitted a-coefficients from observed values, beyond the few-hundredths-of-nHz accuracy obtained at ℓ = 85, would show the background inconsistency matters. A complementary check: rebuild the background to satisfy rotational hydrostatics (as done for polytropic rotating stars) and test whether the ℓ = 85 coefficients shift by more than 0.03 nHz.","tokens_in":59670,"feed_emoji":"☀️","tokens_out":17858,"duration_ms":157703,"temperature":0.7,"pith_summary":"This paper builds a numerical solver for the equations that govern linear adiabatic oscillations of the Sun in a realistic background that includes gravity, acoustic attenuation, and differential rotation, and shows that the simulated wavefields reproduce observed helioseismic data. The authors solve the oscillation equations in first-order form (displacement and pressure perturbation as unknowns), decompose the axisymmetric problem into azimuthal modes, and discretize each modal problem on the meridional plane with a Hybridizable Discontinuous Galerkin (HDG) method. The key quantitative result: for the observed solar rotation profile, the synthetic power spectrum gives frequency-splitting a-coefficients a1/2π = 442.86 nHz, a3/2π = 22.17 nHz, a5/2π = −3.22 nHz, matching measured values 442.85±0.05, 22.15±0.08, −3.21±0.10 nHz. If correct, the solver supplies accurate Green's kernels for local helioseismology—forward modeling of power spectra and cross-covariances that spherical eigensolvers cannot provide once rotation breaks spherical symmetry. The paper also shows that gravity changes the mathematical nature of the wave operator (elliptic versus hyperbolic regions set by the buoyancy frequency), which makes the choice of HDG stabilization load-bearing for accuracy.","feed_headline":"Solar-wave solver matches observed rotation splitting to 0.02 nHz","feed_subtitle":"Computed a-coefficients land within 0.03 nHz of observed values, a key check for helioseismic forward modeling.","key_machinery":"The load-bearing object is the first-order system A u + β1 w + ∇w = g, ∇·u + β2·u + ϱw = h, with u the rescaled Lagrangian displacement and w the pressure perturbation, obtained from the Gough–Thompson equations via a Liouville change of variables. Axisymmetry reduces it to azimuthal modal problems on the meridional half-disk, whose operator relates to a scalar PDE with leading term ∇m·A_m^{-1}∇m. The buoyancy frequency N decides the local type of that operator: elliptic where ω²>N², hyperbolic where ω²<N² (at zero attenuation). The HDG method couples cells through a face trace λ of w and a stabilization parameter τ = |A^{-1}β1·n| + iω α_tune |n^T A^{-1} n|, whose imaginary part receives its","core_discovery":"This paper claims a Hybridizable Discontinuous Galerkin (HDG) scheme for the first-order stellar oscillation equations computes Green's kernels for axisymmetric solar backgrounds with gravity, acoustic attenuation, and differential rotation. Working in the Gough–Thompson formulation (displacement and pressure perturbation as unknowns), the axisymmetric problem splits into independent 2D modal problems on the meridional half-disk. The central validation: for the observed solar rotation profile, the synthetic power spectrum yields a1/2π = 442.86 nHz, a3/2π = 22.17 nHz, a5/2π = −3.22 nHz, matching observed values 442.85±0.05, 22.15±0.08, −3.21±0.10 nHz.","pith_inferences":["The strongest untested consequence of the background inconsistency sits in the deep interior: modes with low degree ℓ reach depths where the neglected ghost force peaks near 3% of the hydrostatic balance. If such modes were simulated, a systematic departure from observed a-coefficients beyond the 0.03 nHz accuracy seen at ℓ = 85 would finger the reference-state error rather than the numerical meth","The space-like/time-like stabilization rule for the hyperbolic region is a general recipe for HDG on first-order systems whose type changes across the domain; it could be transferred to any stratified-wave problem with an effective buoyancy frequency (atmospheric or oceanic internal waves) and tested on idealized profiles before committing to stellar backgrounds.","Because the HDG solve yields the full displacement vector, synthetic observables could move from the pressure proxy to line-of-sight velocity with projection and radiative transfer; the fidelity gain is not quantified in this paper, but the solver architecture makes it a direct next step."],"forward_implications":["Rotation-induced splitting becomes a quantitative validation: the synthetic power spectrum for the observed solar rotation profile yields a1/2π = 442.86 nHz, a3/2π = 22.17 nHz, a5/2π = −3.22 nHz, within 0.02–0.03 nHz of measured values—below the observational noise level.","Green's kernels computed this way can serve as forward models for local helioseismic inversions (time-distance and holography) on axisymmetric backgrounds with realistic flows, which spherically symmetric eigensolver-based approaches cannot supply once rotation breaks spherical symmetry.","Without rotation, the power-spectrum ridges coincide with GYRE eigenfrequencies, cross-checking the Green's-function approach against an independent solver and reproducing the discrete mode structure.","Because the HDG solve provides the full 3D displacement field, synthetic Doppler signals can be improved by including projection and radiative-transfer effects rather than using the pressure perturbation as a proxy.","The same axisymmetric pipeline applies to other rotating bodies—rapidly rotating stars and gas giants—and, at low frequencies, to inertial modes whose sensitivity reaches the deep interior."],"supporting_citations":[{"why":"Supplies the Gough–Thompson form of the oscillation equations (displacement and pressure perturbation unknowns) from which the paper's first-order working system is derived.","marker":"[31]"},{"why":"Provides the Model-S solar background (density, pressure, sound speed, gravity potential) used for all numerical experiments.","marker":"[14]"},{"why":"Supplies the observed solar rotation profile and the measured a-coefficients against which the synthetic power spectrum is validated.","marker":"[39]"},{"why":"GYRE's eigenfrequencies and eigenfunctions provide the independent no-rotation comparison and the Coriolis correction in the analytic a-coefficient estimate.","marker":"[63]"},{"why":"Gives the simplified rotation profile for the analytic comparison and the earlier gravity-free axisymmetric formulation this solver extends.","marker":"[24]"},{"why":"Builds the C2-smooth representation of the solar background coefficients that the numerical implementation requires.","marker":"[19]"},{"why":"The unified hybridization framework on which the HDG discretization of the modal problem is based.","marker":"[15]"},{"why":"Carries the HDG framework (local problems, numerical traces, static condensation) that the paper adapts to the oscillation system.","marker":"[50]"},{"why":"Defines the a-coefficient expansion basis and the rotation-coefficient relation (including the Coriolis correction) used to interpret the fitted splittings.","marker":"[59]"}],"fun_headline_variants":["HDG solar-wave solver nails rotation splitting within 0.02 nHz","New solar oscillation model matches observed rotation to 0.02 nHz","Solar interior model reproduces rotation-induced frequency shifts","HDG method computes solar Green's kernels, matches observed a-coefficients"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The observed rotation is overlaid on the non-rotating, spherically symmetric Model-S background, so the reference state does not satisfy rotational hydrostatics and the solver neglects the implied ghost force f_ghost = ρ0 Ω×Ω×x (up to 3% of the balance, largest near the centre, per Remark 7); the validated match to observed mode splitting depends on that inconsistency being negligible where the compared modes live.","fun_headline_variants_meta":{"raw":{"variants":["HDG solar-wave solver nails rotation splitting within 0.02 nHz","New solar oscillation model matches observed rotation to 0.02 nHz","Solar interior model reproduces rotation-induced frequency shifts","HDG method computes solar Green's kernels, matches observed a-coefficients"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000168,"raw_usage":{"total_tokens":1147,"prompt_tokens":844,"completion_tokens":303,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":227}},"tokens_in":588,"tokens_out":303,"duration_ms":3236,"temperature":1.0,"reasoning_tokens":227,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:48:14.842767+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same pipeline for the lowest harmonic degrees (ℓ ≲ 20), whose acoustic paths reach the deep interior where the neglected ghost force peaks near 3% of the hydrostatic balance; a systematic departure of the fitted a-coefficients from observed values, beyond the few-hundredths-of-nHz accuracy obtained at ℓ = 85, would show the background inconsistency matters. A complementary check: rebuild the background to satisfy rotational hydrostatics (as done for polytropic rotating stars) and test whether the ℓ = 85 coefficients shift by more than 0.03 nHz.","supporting_citations":[{"cited_title":"Howe , Solar interior rotation and its variation , Living Reviews in Solar Physics, 6 (2009), pp","cited_arxiv_id":null,"evidence_quote":"Supplies the Gough–Thompson form of the oscillation equations (displacement and pressure perturbation unknowns) from which the paper's first-order working system is derived."},{"cited_title":"Cockburn, J","cited_arxiv_id":null,"evidence_quote":"Provides the Model-S solar background (density, pressure, sound speed, gravity potential) used for all numerical experiments."},{"cited_title":"Lynden-Bell and J","cited_arxiv_id":null,"evidence_quote":"Supplies the observed solar rotation profile and the measured a-coefficients against which the synthetic power spectrum is validated."},{"cited_title":"Gizon, A","cited_arxiv_id":null,"evidence_quote":"Gives the simplified rotation profile for the analytic comparison and the earlier gravity-free axisymmetric formulation this solver extends."},{"cited_title":"Fournier, N","cited_arxiv_id":null,"evidence_quote":"Builds the C2-smooth representation of the solar background coefficients that the numerical implementation requires."},{"cited_title":"Rouxelin , M \\'e thodes num \\'e riques mixtes condens \\'e es pour l’ \\'e tude de la propagation des ondes acoustiques en \\'e coulement","cited_arxiv_id":null,"evidence_quote":"Carries the HDG framework (local problems, numerical traces, static condensation) that the paper adapts to the oscillation system."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the a-coefficient expansion basis and the rotation-coefficient relation (including the Coriolis correction) used to interpret the fitted splittings."}],"review_version":1}