{"id":"2c46f319-fb0e-41d7-9642-936600ab666e","arxiv_id":"2607.25824","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Solar wind protons experience substantial perpendicular, but not parallel, heating from 0.05 to 1 au, reducing the expected adiabatic growth of temperature anisotropy.","lead":"Using Parker Solar Probe and Solar Orbiter data, this study finds that solar wind protons are heated in the direction perpendicular to the magnetic field between 0.05 and 1 astronomical units, while their parallel temperature shows no clear deviation from adiabatic expansion. The result explains why proton temperature anisotropy grows more slowly than adiabatic theory predicts and how firehose instabilities eventually limit it.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The attribution of the observed C_perp rise to perpendicular heating assumes that heat-flux and non-gyrotropic terms in Eq. (A.12) are negligible; this is not demonstrated, so the central 'heating' claim is not yet secured.","rationale":"The reader's weakest_assumption correctly identifies the Lagrangian reduction and the neglected heat-flux/non-gyrotropic terms as the key risk. I focus specifically on the neglected terms because, even if the radial sequence were perfectly Lagrangian, the measured C_perp increase could be caused by non-heating processes such as heat-flux divergence or non-gyrotropic stresses. This is more directly load-bearing for the abstract's use of the word 'heating' and for the quantitative heating-rate estimates in Fig. 2. However, the reader's concern about stream mixing is also valid and independent; I do not downgrade the overall verdict because the qualitative evidence for non-adiabatic perpendicular evolution is strong and corroborated by previous studies (e.g., Mozer et al. 2023; Marsch et al. 1983). The paper is CONDITIONAL as judged by the reader; my concern reinforces the need for the requested revisions rather than changing the verdict. A concrete heat-flux and non-gyrotropy check is feasible with existing data and would settle whether the 'heating' terminology is justified.","tokens_in":16337,"tokens_out":7342,"duration_ms":73953,"concrete_test":"Compute from the PSP and SolO 3D proton distributions the heat-flux tensor q and non-gyrotropic pressure P_ng for the same radial bins as Fig. 1, and evaluate the terms in Eqs. (A.9)–(A.10), specifically (1/2 ∇·q):bb, Tr(1/2 ∇·q), and the P_ng contributions. If these are collectively smaller than ~10% of u_p,R d/dR ln C_perp over 0.05–0.25 au and 0.3–1 au, the heating interpretation stands; if they are comparable or larger, the inferred Q_perp is not a pure heating rate and the central claim requires rephrasing. A decisive positive control would be to run the same calculation on a synthetic dataset where the true heating rate is known.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the solar wind is heated perpendicularly rests on interpreting the radial increase of ln C_perp (Fig. 1b) as a Lagrangian time derivative d ln C_perp/dt = u_p,R d/dR ln C_perp, and then equating this with Q_perp/P_perp in Eq. (A.12). This is only valid if all other terms in the energy equation — heat-flux divergence, non-gyrotropic pressure, non-ideal electric fields, and collisions — are negligible. The paper does not estimate or bound these terms, despite PSP and SolO measuring the full 3D proton distribution from which the heat-flux tensor and non-gyrotropic pressure can be computed. If, for example, the parallel heat-flux divergence is a significant fraction of the inferred Q_perp, then the 'heating rates' in Fig. 2 and the abstract's 'perpendicular heating' overstate actual energy input, and the comparison to turbulent dissipation rates (Sec. 5) is not clean. Moreover, the radial sequence mixes different plasma parcels; the per-bin percentile splitting (Appendix B) does not tag the same parcel, so the rise could reflect a changing population mix. Both issues make the identification of the observed invariant growth as 'heating' load-bearing and currently unverified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript combines Parker Solar Probe (0.05–0.25 au) and Solar Orbiter (0.3–1 au) proton moments to study the radial evolution of the two CGL adiabatic invariants C_parallel = T_parallel (B/n)^2 and C_perp = T_perp/B (Eq. 1). Splitting the data into fast and slow wind using 3–33% speed percentiles per radial bin, the authors report that ln C_parallel is approximately constant with distance while ln C_perp increases for both populations (Fig. 1). Interpreting this increase through d/dt = u_R d/dR and Eq. (3), they convert the radial slope into perpendicular heating-rate densities Q_perp and per-mass rates epsilon_perp (Fig. 2), finding rates that decrease with distance and are larger in the fast wind. They then argue that this perpendicular heating slows the double-adiabatic decrease of T_perp/T_parallel and shapes the plasma distribution in the beta_parallel–T_perp/T_parallel plane, with the slow wind approaching the firehose thresholds and the fast wind staying near the proton-cyclotron/mirror conditions (Figs. 3–5). The central claim is that solar-wind protons experience significant perpendicular heating throughout 0.05–1 au while remaining adiabatic in the parallel direction.","tokens_in":16660,"tokens_out":6758,"duration_ms":64256,"significance":"If the heating attribution is correct, the paper extends evidence for non-double-adiabatic proton behavior to 0.05 au, bridges PSP and Solar Orbiter observations, and provides quantitative heating rates that can be compared with turbulent dissipation estimates. The main strengths are the large multi-spacecraft dataset, the explicit percentile-based population separation that accounts for wind acceleration, and the use of two fitting forms to show that the qualitative rise of ln C_perp and flatness of ln C_parallel are robust. The paper is important, but the central inference depends on assumptions that are not demonstrated: the neglected heat-flux, non-gyrotropic, and population-mixing effects, and the quantitative rates lack a statistical uncertainty budget and are affected by a large PSP/SolO offset. These issues are addressable and do not, in my view, invalidate the qualitative result, but they are load-bearing for the abstract's 'significant average heating' statement and for the comparison with turbulence rates.","major_comments":[{"comment":"The central inference equates u_R d ln C_perp/dR with Q_perp/P_perp and then calls Q_perp a 'perpendicular heating rate.' But the text states that Q_perp includes heat-flux divergence, non-ideal electric fields, and non-gyrotropic pressure-tensor contributions. Equation (A.10) explicitly contains (1/2 ∇·q):bb and Tr(P_ng·∇u + ...), and none of these terms is estimated or bounded in the manuscript. If, for example, the heat-flux divergence is a non-negligible fraction of the advective term, the inferred 'heating' is not local energy dissipation. Since both PSP and SolO measure full 3D distribution functions, the authors should provide at least order-of-magnitude estimates of the heat-flux and agyrotropy terms, or alternatively use a more cautious term such as 'effective non-adiabatic rate' and avoid the direct comparison with turbulent dissipation in Section 5 until the neglected terms ar","section":"Appendix A / Eq. (A.12); Section 3, Fig. 2"},{"comment":"The reduction d/dt = u_R ∂_R treats the radial sequence of binned measurements as a Lagrangian time derivative along a single steady flow. The data are not Lagrangian: PSP and SolO sample different plasma parcels at different times, and Appendix B shows that the selected wind populations accelerate by about 150–200 km/s between 0.05 and 1 au. The 3–33% percentile split is applied independently to each radial bin, so it does not track the same parcel; a radial change in the source mix within the slow/fast categories could produce a rise in ln C_perp without any heating. This is load-bearing because Eq. (3) is only valid if the same population is followed and if the non-advective terms are negligible. The authors should test for population-mixing effects (e.g., identify and track individual stream intervals, back-map plasma parcels, or compare with a constant-speed-threshold split) or soft","section":"Section 3 / Appendix B"},{"comment":"The quantitative heating rates carry no statistical uncertainties. The linear and power-law fits are applied to binned means without reporting confidence intervals, goodness-of-fit, or sensitivity to the chosen bin size. More importantly, the PSP- and SolO-based rates differ by a factor of about 5–10 in epsilon_perp at R ≈ 0.25–0.3 au, and the text attributes this to calibration without quantification. As a result, the combined radial trend quoted in Section 5 (epsilon from about 10^5 to 5×10^3 J s^-1 kg^-1) is not established at a stated precision. The authors should propagate fit uncertainties and either calibrate the PSP/SolO offset using the overlap or present the two spacecraft's rates separately rather than as one continuous trend. This is needed before the rates can be compared quantitatively with turbulent dissipation estimates.","section":"Section 3 / Fig. 2"}],"minor_comments":[{"comment":"The sentence stating SolO data are taken at 'R ∈ [0.05, 0.25] au' appears to be a typo; the rest of the paper consistently uses R ∈ [0.3, 1] au for SolO. Please correct.","section":"Section 2"},{"comment":"In the sentence 'both C_p⊥ and C_p⊥∥ are conserved', the second symbol should be C_p∥, not C_p⊥∥.","section":"Section 2"},{"comment":"The adiabatic CGL curves in Figs. 3–5 are computed from power-law exponents s(B) and s(n) fitted to the same dataset. The table note already cautions about interpretation, but the main text should state explicitly that these are data-derived baselines, not independent theoretical predictions, so that readers do not mistake the comparison for a test of double-adiabatic theory against an external model.","section":"Table 1"},{"comment":"The marginal-stability curves depend on the assumptions of bi-Maxwellian protons, beta_e = 1, and gamma = 10^-3 Omega_ci. The text mentions some limitations, but the figures would benefit from a sentence in the captions reiterating that the thresholds are not exact boundaries for the observed non-Maxwellian distributions.","section":"Section 4 / Figs. 4–5"}],"recommendation":"major_revision","confidential_remarks":"The paper is worth publishing if the authors can secure the central attribution of the ln C_perp rise to perpendicular heating. The qualitative pattern is likely robust, but the abstract's 'significant average heating' statement currently rests on unquantified neglect of heat-flux and agyrotropy terms and on a Lagrangian interpretation of non-Lagrangian measurements. The missing uncertainty budget and the large PSP/SolO offset also prevent the quantitative rates from being used as they are. I therefore recommend major revision rather than rejection: the authors have access to the data needed to estimate the neglected terms and to quantify the population-mixing and calibration issues. There is no novelty concern relative to Mozer et al. (2023) because the present work adds a systematic fast/slow wind separation and the full 0.05–1 au range; however, the overlap discontinuity must be handled transparently."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe key thing to know: this is a well-executed, consolidating study. It confirms what Helios and early PSP already suggested — perpendicular proton heating, no parallel heating, firehose-regulated anisotropy — but extends it to 0.05 au with a large PSP dataset and a coordinated SolO comparison, and it gives radial heating-rate estimates for fast and slow wind separately. The central measurement (ln C_perp rising with distance, ln C_par flat) is robust across binning and population choices. The authors are careful with FOV limitations, use QTN densities for PSP, and split wind populations by percentile to handle acceleration. That is good work.\n\nThe new contribution is mainly quantitative: heating rates Q_perp and epsilon_perp from 0.05 to 1 au, split by speed, and the comparison to turbulent dissipation rates. The rates are plausible and consistent with Mozer et al. and Kontar et al., though an order of magnitude below Bowen et al. The anisotropy evolution in the beta-T plane is a nice summary, and the discussion of instability thresholds (including alpha effects) is honest.\n\nSoft spots: the heating rates have no propagated statistical uncertainties — just two fit choices. The PSP/SolO offset in anisotropy (factor ~1.3) and in heating rates (factor 5-10) is acknowledged but not corrected, so the combined radial profile is qualitative across the joint. The more foundational issue: the authors identify d ln C_perp/dt with u d/dR and attribute the residual to Q_perp, assuming heat-flux divergence, non-gyrotropic pressure, and collisions are negligible. They do not estimate those terms, even though PSP and SolO measure the full 3D distributions. That is a real caveat, and the stress-test is right to flag it. But I don't think it breaks the paper: the same pattern appears in Helios, the parallel invariant is flat, and the heating explanation is consistent with the observed anisotropy evolution. The heat-flux term would have to be large and anisotropic to produce this. Still, a referee should ask them to bound it, or at least discuss why it's small.\n\nNo code or processed data is released, which is a reproducibility annoyance.\n\nBottom line: this is a useful reference paper for heliophysics. It deserves a serious referee — conditional acceptance with requests for uncertainty treatment and a clearer discussion of the neglected terms. I'd cite it for the heating-rate profiles.","headline":"Useful consolidation of perpendicular proton heating with new quantitative radial profiles; core claim holds, but rates need uncertainty and heat-flux caveats.","tokens_in":17189,"tokens_out":3042,"would_cite":true,"duration_ms":30946,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Solar wind protons are heated perpendicular to the magnetic field from 0.05 to 1 au, and this heating substantially slows the growth of temperature anisotropy predicted by double-adiabatic expansion.","keywords":[],"falsifier":"A measurement that follows the same plasma parcel or streams as they move outward (for instance, using a radial-aligned pair of spacecraft, or tracking features in the distribution function) that shows C_⊥ not increasing with distance would refute the claim of genuine perpendicular heating. Alternatively, a full evaluation of the heat-flux term in the C_⊥ evolution equation that shows it accounts for most of the observed dC_⊥/dR would falsify the attribution to heating.","tokens_in":16216,"feed_emoji":"","tokens_out":2023,"duration_ms":23386,"temperature":0.7,"pith_summary":"This paper tries to establish that, across the entire inner heliosphere sampled by Parker Solar Probe and Solar Orbiter, solar wind protons receive non-adiabatic heating in the direction perpendicular to the local magnetic field, while the parallel direction behaves adiabatically. The perpendicular heating is strong enough to significantly reduce the temperature anisotropy that would develop if the plasma expanded double-adiabatically. Despite this heating, slower wind streams still develop a parallel-dominated anisotropy that becomes constrained by kinetic firehose instabilities. If correct, the solar wind is not double-adiabatic even close to the Sun, and perpendicular energy input is a major control on the plasma's stability evolution.","feed_headline":"","feed_subtitle":"","key_machinery":"The central objects are the two adiabatic invariants C_∥ = T_∥(B/n)² and C_⊥ = T_⊥/B, whose log-time derivatives equal the non-adiabatic heating rates normalized by pressure: (1/2)d/dt ln C_∥ = Q_∥/P_∥ and d/dt ln C_⊥ = Q_⊥/P_⊥. Under the double-adiabatic (CGL) approximation these invariants are exactly conserved; any measured secular change in them with heliocentric distance directly signals non-adiabatic heating or cooling. The paper exploits this by converting the radial sequence of binned spacecraft measurements into a temporal derivative via d/dt = u_R d/dR, allowing it to quantify the heating rates from the observed slopes of ln C_∥ and ln C_⊥.","core_discovery":"By monitoring the radial evolution of the two Chew-Goldberger-Low adiabatic invariants, C_∥ = T_∥(B/n)² and C_⊥ = T_⊥/B, in Parker Solar Probe (0.05–0.25 au) and Solar Orbiter (0.3–1 au) data, the paper finds that C_⊥ increases with distance for both slow and fast solar wind, definitively indicating non-adiabatic perpendicular proton heating. C_∥, by contrast, shows no clear deviation from a constant value, indicating no average parallel heating or cooling. The perpendicular heating rates are higher in the faster wind and decrease with distance, but across all distances the heating is sufficient to make the observed decrease of T_⊥/T_∥ slower than the double-adiabatic prediction. Consequentl","pith_inferences":["The paper's finding that T_⊥/T_∥ ≳ 1 for both wind populations inside 0.075 au, extrapolated sunward, implies that the high corona likely also has T_⊥ > T_∥; this is a testable constraint for coronal heating models, though it rests on the same radial-sequence assumption.","The authors leave open the physical mechanism; a natural testable extension is to sort the same data by Alfvénicity (as they suggest) to see whether the perpendicular heating is concentrated in Alfvénic slow wind, which would link the heating to switchbacks or coherent turbulent structures.","The jump in C_⊥ and heating rates between PSP and SolO data near 0.3 au is attributed to calibration; if a cross-calibration or a better-aligned radial dataset were available, the true radial profile of C_⊥ might be smoother, and the inferred total heating rate could be revised."],"forward_implications":["If the perpendicular heating is real, then standard double-adiabatic (CGL) models of the inner heliosphere are inadequate even close to the Sun; models must include a perpendicular energy source.","The reduced rate of anisotropy growth means the onset of firehose and other kinetic instabilities is delayed or weakened relative to adiabatic predictions, changing where in the solar wind unstable conditions first appear.","Because the heating rates are higher in fast wind than slow wind (by a factor of ~5 in specific heating rate), the mechanism responsible likely scales with wind speed or turbulence amplitude.","The observed heating rate near 0.05 au (~10⁵ J s⁻¹ kg⁻¹) matches independent estimates of turbulent energy dissipation, suggesting that turbulence is a viable and possibly dominant heating source at these distances."],"fun_headline_variants":["Solar wind protons heat sideways, not along field lines","Perpendicular proton heating slows solar wind anisotropy growth","Proton heating reduces anisotropy but firehose instability persists","Radial proton heating: perpendicular wins, parallel stays adiabatic","Solar wind proton heating dims temperature anisotropy, not firehose"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The identification of the observed increase in C_⊥ with distance as 'heating' assumes that the binned radial sequence of measurements represents the Lagrangian evolution of the same plasma population (d/dt = u_R d/dR), and that neglected heat-flux and non-gyrotropic pressure terms are small—an assumption that is stated but not demonstrated.","fun_headline_variants_meta":{"raw":{"variants":["Solar wind protons heat sideways, not along field lines","Perpendicular proton heating slows solar wind anisotropy growth","Proton heating reduces anisotropy but firehose instability persists","Radial proton heating: perpendicular wins, parallel stays adiabatic","Solar wind proton heating dims temperature anisotropy, not firehose"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000167,"raw_usage":{"total_tokens":1073,"prompt_tokens":699,"completion_tokens":374,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":443,"completion_tokens_details":{"reasoning_tokens":293}},"tokens_in":443,"tokens_out":374,"duration_ms":4125,"temperature":1.0,"reasoning_tokens":293,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T01:21:25.742998+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement that follows the same plasma parcel or streams as they move outward (for instance, using a radial-aligned pair of spacecraft, or tracking features in the distribution function) that shows C_⊥ not increasing with distance would refute the claim of genuine perpendicular heating. Alternatively, a full evaluation of the heat-flux term in the C_⊥ evolution equation that shows it accounts for most of the observed dC_⊥/dR would falsify the attribution to heating.","supporting_citations":[],"review_version":1}