{"id":"35af072f-9519-4f9d-8733-d84683e6d56c","arxiv_id":"2607.11880","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Empirical velocity-space diffusion coefficients from PSP proton VDFs match scale-dependent intermittent stochastic heating and rule out parallel-ICW resonant heating and helicity-barrier SH in a sub-Alfvénic stream.","lead":"Parker Solar Probe data were inverted to measure how solar-wind protons diffuse in velocity space near the Sun. The measured heating matches intermittent stochastic heating, not cyclotron resonance or helicity-barrier models.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"The inversion of the guiding-center equation presupposes a pure Fokker–Planck diffusion operator; residual non-diffusive terms can bias D⊥⊥ and the velocity-space peak of dQ_emp that underpins the SH discrimination.","rationale":"The paper’s strongest claim is a genuine methodological advance: an empirical velocity-space heating rate that can be compared directly with analytic SH and RH expressions without the intermediate LET assumptions of B25. The reader’s weakest-assumption (stable radial parcel) is real and correctly identified, but it is secondary; even with perfect connectivity the inversion still embeds the pure-diffusion ansatz of Eq. 1. That ansatz is what converts the observed VDF evolution into a D⊥⊥ whose peak location and amplitude then select intermittent SH over J25 and ICW. Because the paper also presents the mere existence of a finite D_ii as evidence for a Fokker–Planck process, the circularity is load-bearing for the discrimination claim. The proposed synthetic-residual test is concrete, uses the same data pipeline, and would settle whether the peak location is an artifact of the operator form. The verdict therefore remains CONDITIONAL (method and stream-specific discrimination accepted once the robustness check and multi-interval statistics are supplied), but the primary soft spot is the diffusion-operator assumption rather than solely the connectivity premise. No contradiction or pure fitting is present, so REJECT is unwarranted.","tokens_in":18177,"tokens_out":766,"duration_ms":7103,"concrete_test":"Inject controlled non-diffusive residuals (e.g., a small CGL-violating source term or a synthetic non-gyrotropic moment scaled to 5–10 % of the observed ∂f/∂r) into the left-hand side of Eq. 1 for the same 116 windows, re-invert for D⊥⊥, recompute dQ_emp, and check whether the peak remains within 0.2 v_th,⊥ of 1.1 v_th,⊥ and whether the intermittent-SH match still holds at the reported magnitude. If either fails, the discrimination is not robust.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim rests on the empirical dQ_emp (Eq. 2) peaking near 1.1 v_th,⊥ and matching intermittent scale-dependent SH while disfavoring J25 and ∥-ICW RH (Fig. 3). That dQ_emp is obtained by solving the steady-state gyrotropic guiding-center equation (Eq. 1) for the three diffusion coefficients via SVD after RBF interpolation of SPANi VDFs. The right-hand side of Eq. 1 is written exclusively as a Fokker–Planck operator; any residual non-diffusive contributions (finite-Larmor-radius corrections, non-gyrotropic moments, unmodeled large-scale forces, or residual FOV/instrumental systematics) are absorbed into the inverted D_ii. Because the same operator form is later cited as “evidence of a Fokker–Planck like process,” the match to SH is not fully independent of the modeling assumption. The reader correctly flags the single-stream radial-evolution premise, but the deeper load-bearing issue is that the inversion itself can force a diffusion-like solution even when the true evolution contains non-diffusive terms, thereby locking the location and magnitude of the dQ_emp peak that discriminates among heating models.","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.5","summary":"The manuscript introduces an inversion of the steady-state gyrotropic proton guiding-center equation (Eq. 1) against RBF-interpolated SPANi VDFs to obtain empirical velocity-space diffusion coefficients D_ii(v_perp, v_parallel) in a carefully selected two-hour low-beta, highly imbalanced sub-Alfvenic PSP Encounter-8 interval. From the measured D_perpperp the authors form a fully kinetic perpendicular heating rate dQ_emp (Eq. 2) and show that it peaks near 1.1 v_th,perp and matches both the magnitude and velocity-space location of the intermittent, scale-dependent stochastic-heating (SH) expression, while the rms SH, the J25 helicity-barrier SH, and the quasilinear parallel-ICW resonant heating rates neither peak in the same region nor reach the required amplitude (Fig. 3). The work therefore claims both a novel observational methodology for constraining collisionless heating and direct evidence for a Fokker-Planck-like diffusive process dominated by intermittent SH in the near-Sun wind.","tokens_in":18468,"tokens_out":1385,"duration_ms":16818,"significance":"If the central comparison holds, the paper supplies a genuinely new, phase-space-resolved diagnostic that can be applied to other PSP intervals, multi-spacecraft data sets, and laboratory plasmas, thereby moving beyond scalar cascade-rate or moment-based estimates. The explicit discrimination among three SH formulations and ICW heating on the same velocity grid, together with the demonstration that intermittency is required for magnitude agreement, constitutes a falsifiable advance over prior work (including Bowen et al. 2025). The machine-readable inversion pipeline and the use of centered finite differences that conserve particle number are methodological strengths that the community can reuse.","major_comments":[{"comment":"Data section and Eq. (1): the entire inversion rests on the premise that the two-hour interval samples a single, radially evolving plasma parcel from a stable coronal-hole source so that successive SPANi VDFs can be used to evaluate the advection and CGL terms while gravity and the ambipolar electric field are neglected. The connectivity argument is qualitative; no quantitative test (e.g., ballistic mapping uncertainty, non-radial flow residuals, or comparison with a second independent stream) is provided. If the premise fails, the inverted D_perpperp and the location of the dQ_emp peak that discriminates among models are systematically biased.","section":"Data section; Eq. (1)"},{"comment":"Eq. (1) and Discussion: the right-hand side is written exclusively as a Fokker-Planck diffusion operator. Any residual non-diffusive contributions (finite-Larmor-radius corrections, residual non-gyrotropy, unmodeled large-scale forces, or FOV/instrumental systematics) are absorbed into the inverted D_ii. The subsequent claim of “evidence of a Fokker-Planck like process” is therefore not fully independent of the modeling assumption that was imposed a priori. A quantitative residual analysis or a controlled synthetic-data test demonstrating that non-diffusive terms do not shift the dQ_emp peak is required before the discrimination in Fig. 3 can be regarded as robust.","section":"Eq. (1); Discussion"},{"comment":"The analysis is performed on a single, carefully chosen two-hour interval that yields only 116 independent D_ii measurements. While the interval is well-resolved and FOV-clean, the central claim that intermittent SH is the dominant mechanism in the sub-Alfvenic wind cannot be generalized from one stream. At minimum the authors should demonstrate that the same peak location and magnitude ordering persist in at least one additional independent sub-Alfvenic interval, or clearly reframe the result as a case study.","section":"Results; Fig. 3"},{"comment":"End Matter, Eqs. (6)–(8): the SH coefficients c1 = 0.75 and c2 = 0.34 are taken from earlier test-particle work and held fixed. The paper notes that raising c1 by a factor of ~4 would bring the J25 rate up to the observed magnitude, yet still leaves the peak at the wrong velocity. Because the empirical D_perpperp are now available, a direct least-squares fit for c1 and c2 (or a Bayesian posterior) on the same velocity grid would remove this free-parameter ambiguity and strengthen the claim that only the intermittent scale-dependent form is viable.","section":"End Matter; Eqs. (6)–(8)"}],"minor_comments":[{"comment":"Fig. 3 normalizes all curves to the maximum of dQ_emp; absolute heating rates (or a second panel with absolute units) would allow direct comparison with the LET rates of Bowen et al. 2025 and with independent cascade-rate estimates.","section":"Fig. 3"},{"comment":"The flux-tube area A(r) ~ R^6/(R^4+6) is adopted without sensitivity tests; a brief check with the simpler A ~ r^2 would quantify the impact on the left-hand side of Eq. (1).","section":"Eq. (1)"},{"comment":"Notation for the intermittent versus rms SH rates (D_SH,int vs D_SH,rms) is introduced only in the text surrounding Fig. 3; a short table or explicit definitions earlier would improve readability.","section":"Comparison with Stochastic..."},{"comment":"The secondary peak at ~3 v_th,perp in the intermittent SH curve is attributed to a breakdown of the SH model; a quantitative estimate of the velocity at which the chaotic-drift assumption fails would make this statement more precise.","section":"Fig. 3 caption / Discussion"}],"recommendation":"major_revision","confidential_remarks":"The methodological novelty is real and the comparison in Fig. 3 is carefully executed; the main risk is over-generalization from a single stream plus the circularity inherent in assuming a pure diffusion operator. If the authors can add even one additional interval and a residual/synthetic test, the paper would be a strong PRL candidate. The citation pattern is appropriate and does not appear to suppress competing work."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real advance here is the inversion itself. They take successive SPANi 3-D VDFs, RBF-interpolate them, form the left-hand side of the steady gyrotropic guiding-center equation (advection + CGL), and solve for the three diffusion coefficients by SVD. That gives an empirical D⊥⊥(v⊥,v∥) and a fully kinetic perpendicular heating rate without presupposing the SH formula. The subsequent comparison (Fig. 3) is clean: the empirical dQ peaks near 1.1 vth,⊥ and matches the intermittent, scale-dependent SH expression in both location and amplitude; the J25 helicity-barrier version and the quasilinear ∥-ICW rate do neither. That is a genuine observational constraint that prior moment- or spectrum-based estimates could not deliver.\n\nThe data handling looks careful—FOV metric, centered differences that conserve density, consistent velocity grids, error shading on the empirical curve. The theory side cites the right papers (Chandran, Klein & Chandran, Mallet, Isenberg, Bowen et al. 2025, Johnston et al.) and uses the standard c1/c2 values from test-particle work. No obvious citation padding or missing counter-literature.\n\nSoft spots are real but proportionate. Everything rests on one two-hour, low-β, highly imbalanced interval chosen for stable coronal-hole connectivity and good FOV; if the parcel is not purely radial or if non-diffusive terms (FLR, residual non-gyrotropy, unmodeled large-scale forces) leak into the RHS, the inverted D absorbs them. That makes the later claim of “evidence for a Fokker–Planck-like process” mildly circular—the operator form was assumed. Flux-tube area is a model, gravity and ambipolar field are dropped, and there is no public code. None of these sink the central discrimination for this stream; they just mean the result is a strong first demonstration rather than a settled statistical fact.\n\nThis is for people who care about kinetic heating diagnostics in the solar wind or who want an observational test of SH versus cyclotron resonance. The method is worth having. I would send it to peer review without hesitation; a good referee will push for multi-interval statistics and a clearer separation between assumed operator and claimed evidence, but the paper already earns that attention.","headline":"New inversion of the guiding-center equation yields empirical D⊥⊥ and dQ that cleanly favor intermittent scale-dependent SH over J25 and ∥-ICW RH in one sub-Alfvénic stream, with the usual single-interval and pure-diffusion caveats.","tokens_in":19198,"tokens_out":661,"would_cite":true,"duration_ms":14965,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Direct inversion of solar-wind proton distributions shows intermittent stochastic heating matches the observed velocity-space heating rate.","keywords":["stochastic heating","velocity-space diffusion","Parker Solar Probe","collisionless plasma","solar wind","guiding-center equation","ion cyclotron waves","intermittency"],"falsifier":"Repeat the identical inversion on another well-resolved sub-Alfvénic stream whose source connectivity is independently verified; if the empirical dQ_⊥ peak moves away from ~1.1 v_th,⊥ or no longer matches the intermittent stochastic-heating curve, the claimed identification fails.","tokens_in":19009,"feed_emoji":"☀️","tokens_out":681,"duration_ms":5721,"temperature":0.7,"pith_summary":"Collisionless plasmas such as the near-Sun solar wind cool more slowly than adiabatic expansion allows and show strong perpendicular temperature anisotropies that require continuous preferential heating. This paper asks whether that heating is stochastic (random kicks that break magnetic-moment conservation) or resonant with parallel ion-cyclotron waves. Using Parker Solar Probe measurements of three-dimensional proton velocity distributions in a low-beta, highly imbalanced sub-Alfvénic stream, the authors invert the guiding-center equation to extract the perpendicular velocity-space diffusion coefficient and the associated phase-space heating rate without assuming a particular mechanism. The measured heating rate peaks near 1.1 thermal speeds and has a magnitude that matches only the scale-dependent stochastic-heating formula once intermittent large-amplitude fluctuations are retained. Competing expressions that incorporate a helicity barrier or quasilinear cyclotron damping peak in the wrong region of velocity space and fall short in amplitude. The result supplies direct kinetic evidence for a Fokker–Planck-like diffusive process and a practical method for identifying heating channels in any collisionless plasma that can be observed with full distribution functions.","feed_headline":"Solar-wind protons reveal intermittent stochastic heating","feed_subtitle":"Direct diffusion measurements match theory only when large fluctuations are kept","key_machinery":"Inversion of the steady-state proton guiding-center equation (advection plus CGL terms equal velocity-space diffusion) via singular-value decomposition on nine neighboring distribution-function samples, yielding D_⊥⊥(v_⊥,v_∥) and, after integration by parts, the fully kinetic heating rate dQ_⊥(v_⊥).","core_discovery":"When the proton guiding-center equation is inverted on successive three-dimensional SPANi distributions, the empirically recovered perpendicular heating rate dQ_⊥ peaks near 1.1 v_th,⊥ and equals the magnitude predicted by scale-dependent stochastic heating only when intermittency is retained; helicity-barrier stochastic heating and parallel-ion-cyclotron resonant heating neither peak at the same velocity nor reach the required amplitude.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Direct diffusion coefficients reveal intermittent stochastic heating","Proton VDF inversions match scale-dependent stochastic heating rates","Empirical phase-space rates confirm SH over RH in solar-wind stream","Intermittency-retained SH peaks match measured perpendicular heating","Guiding-center inversion isolates Fokker-Planck diffusion near Sun"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The two-hour interval must be a single, radially evolving plasma parcel from a stable coronal-hole source so that the observed change in the distribution can be treated as a pure radial derivative with large-scale forces neglected.","fun_headline_variants_meta":{"raw":{"variants":["Direct diffusion coefficients reveal intermittent stochastic heating","Proton VDF inversions match scale-dependent stochastic heating rates","Empirical phase-space rates confirm SH over RH in solar-wind stream","Intermittency-retained SH peaks match measured perpendicular heating","Guiding-center inversion isolates Fokker-Planck diffusion near Sun"]},"model":"grok-4.5","effort":"low","cost_usd":0.003528,"raw_usage":{"total_tokens":1210,"prompt_tokens":885,"num_sources_used":0,"completion_tokens":86,"cost_in_usd_ticks":35280000,"prompt_tokens_details":{"text_tokens":885,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":239,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":885,"tokens_out":86,"duration_ms":2530,"temperature":1.0,"reasoning_tokens":239,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T02:27:20.308676+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Repeat the identical inversion on another well-resolved sub-Alfvénic stream whose source connectivity is independently verified; if the empirical dQ_⊥ peak moves away from ~1.1 v_th,⊥ or no longer matches the intermittent stochastic-heating curve, the claimed identification fails.","supporting_citations":[],"review_version":1}