{"id":"f4a2bd6b-fe08-4a87-8bc3-c07a94db0ea0","arxiv_id":"2501.04986","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"Using a CR transport code with prescribed parallel and perpendicular diffusion laws, the authors find power-law energy-density profiles, but the headline anisotropy scaling is recovered from the input rather than derived.","lead":"This paper simulates cosmic-ray electron transport around a pulsar wind nebula using prescribed magnetic-turbulence diffusion laws, and reports power-law scalings for the electron energy density. A reader may care because the scalings, if real, would connect gamma-ray halo observations to magnetic turbulence parameters, but the key anisotropy law is an input, not a discovery.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed emergent scalings are prescribed inputs: Appendix A Step 2 hard-wires chi = MA^alpha into CRIPTIC, and the UCR~MA^-6 result is entangled with an inverted effective-radius transformation; the central 'intrinsic property' claim is unsupported.","rationale":"The reader correctly identified the hard-wired input chi = MA^alpha as the weakest assumption. The stress-test confirms this: the simulation measures a diffusion anisotropy ratio that is fixed by the code's sub-grid parameters, so the claim that D_perp/D_parallel ~ MA^alpha 'maintains its features' over a long evolution is circular. The additional load-bearing concern is the effective-radius transformation in Sect. 4.1. The paper defines r = sqrt((x/sqrt(chi))^2 + y^2 + z^2), but the anisotropic diffusion equation with D_perp = chi*D_parallel requires either scaling the parallel coordinate by sqrt(chi) or scaling the perpendicular coordinates by 1/sqrt(chi) to reduce to isotropic diffusion. The paper's transformation is inverted, so the radial profiles and the volume averages used for UCR ~ MA^-6 are computed in a nonphysical metric. This matters because the abstract's headline result is that the MA^-6 power law is an intrinsic property independent of energy and losses. If the Green's function analysis with the correct radius does not reproduce the slope, the claim collapses; even if it does reproduce the slope, the relation is not intrinsic but a derived consequence of the assumed input. Either way, the central scientific advance is not supported as presented. The radiative-loss and spectral-index effects on morphology are plausibly valid, but those are secondary to the paper's main claim. The recommendation is REJECT because the central emergent-scaling claim is unsupported and the effective-radius transformation is internally inconsistent.","tokens_in":18088,"tokens_out":5319,"duration_ms":54059,"concrete_test":"Compute the no-loss, continuous-point-source Green's function solution of Eq. (1) with uniform D_perp = chi*D_parallel, and evaluate the averaged UCR(MA) using the correct isotropic radius R_eff = sqrt(x^2 + (y^2+z^2)/chi). If the predicted slope is still approximately -6, the scaling is an analytic consequence of the prescribed input rather than an emergent numerical result; if the slope changes, the reported MA^-6 scaling is an artifact of the inverted effective radius.","verdict_should_be":"REJECT","load_bearing_attack":"The central claim that D_perp/D_parallel ~ MA^alpha and UCR ~ MA^-6 survive as intrinsic properties is not established, because the first relation is put in by hand and the second is entangled with an inverted coordinate transformation. In Appendix A Step 2, CRIPTIC is initialized with chi = D_perp/D_parallel = MA^alpha, with alpha = 4 sub-Alfvenic and alpha = 3 super-Alfvenic. These sub-grid diffusion coefficients are fixed inputs; the simulation never evolves the turbulent magnetic field. Recovering D_perp/D_parallel ~ MA^4.0 and MA^3.0 in Fig. 8(c) is therefore tautological rather than a numerical discovery, and no independent test-particle MHD calculation is performed. Separately, Sect. 4.1 defines the effective radius as r = sqrt((x/sqrt(chi))^2 + y^2 + z^2). For Eq. (1) with D_perp = chi D_parallel, the transformation that makes transport isotropic is either X = x*sqrt(chi), giving R_eff^2 = chi*x^2 + y^2 + z^2, or scaling the perpendicular coordinates by 1/sqrt(chi), giving R_eff^2 = x^2 + (y^2+z^2)/chi. The paper's r rescales the parallel coordinate in the opposite direction, so the reported UCR(r) profiles and the volume average leading to UCR ~ MA^-6.0+/-0.5 are computed in a distorted metric. The MA^-6 scaling may therefore be a coordinate artifact of the inverted radius rather than a physical relation. This directly undermines the abstract's claim that the MA^-6 law is independent of energy and radiative losses.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents CRIPTIC simulations of high-energy electron transport in a pulsar-wind-nebula environment, using prescribed anisotropic diffusion coefficients D_∥ = D0(E/mec^2)^δ and D_⊥/D_∥ = M_A^α with α = 4 (sub-Alfvénic) and α = 3 (super-Alfvénic). Three groups of runs vary radiative losses, the turbulence spectral index γ, and the Alfvén Mach number M_A. The authors report power-law radial and spectral energy-density profiles (UCR ∝ r^-6/5, UCR ∝ 1/E without losses; steeper with losses), M_A- and viewing-angle-dependent morphologies, slow-to-normal diffusion transitions, an averaged relation UCR ∝ M_A^-6.0±0.5, and persistence of the input diffusion scalings over 3 kyr. They interpret the last two as intrinsic properties connecting CR transport to MHD turbulence.","tokens_in":18489,"tokens_out":10615,"duration_ms":103618,"significance":"If the headline scalings were genuinely emergent, the paper would be valuable: TeV halo morphology and the UCR–M_A relation could then be used as observational probes of turbulence magnetization and spectral index. The paper has strengths: it uses a documented transport code, runs a large parameter sweep (90 simulations), includes realistic electron cooling (synchrotron plus inverse Compton in the Klein-Nishina regime), and provides clear diagnostics of morphology and anisotropy. The qualitative findings that radiative losses suppress anisotropy, steepen radial profiles, and produce sub-diffusive behavior for PeV electrons are plausible and potentially useful. However, the two central claims are not supported by the simulation design: the D_⊥/D_∥ vs. M_A scaling is an input rather than a discovery, and the UCR ∝ M_A^-6 law is entangled with an inverted coordinate transformation. The paper's contribution is therefore substantially reduced, and the abstract's 'intrinsic property' conclusion is not justified.","major_comments":[{"comment":"Appendix A Step 2 hard-wires χ = D_⊥/D_∥ = M_A^α as a CRIPTIC input (α = 4 sub-Alfvénic, α = 3 super-Alfvénic), and the turbulent magnetic field is never evolved. Section 4.4 and Fig. 8(c) then report that the measured D_⊥/D_∥ follows M_A^4.0 and M_A^3.0 and call this a 'numerical finding' and an 'inherent property.' Because the particle deviations used to compute D_∥ and D_⊥ are generated by the very same prescribed diffusion tensor, the agreement is a consistency check of the code, not an independent confirmation. The same objection applies to the persistence of D_∥ ∼ E^δ: Eq. (2) is the prescribed parallel diffusion law. The abstract's claim that the distribution law is an intrinsic property therefore rests on circular reasoning.","section":"§2.2, §4.4, Appendix A Step 2"},{"comment":"The effective radius is defined as r² = x²/χ + y² + z² with χ = M_A^α. For Eq. (1) with D_⊥ = χD_∥, the change of variables that makes the spatial operator isotropic is X = x, Y = y/√χ, Z = z/√χ, giving R_iso² = x² + (y² + z²)/χ; scaling the parallel coordinate x by 1/√χ instead is the inverse transformation. Consequently, the UCR(r) profiles in Figs. 1–2 and the volume-averaged UCR in Fig. 8(a) are computed in a distorted metric, and the reported UCR ∝ M_A^-6.0±0.5 law (Summary item 4) is not established as a physical scaling. The radial power-law indices UCR ∝ r^-6/5 and r^-11/5 are likewise affected by this coordinate choice.","section":"§4.1, Eq. (1), footnote 3"},{"comment":"The averaging procedure that produces UCR for each M_A is not defined, and the M_A^-6.0±0.5 scaling is not derived from the transport solution. For a fixed-time anisotropic Gaussian profile with D_⊥ = χD_∥ and χ = M_A^α, the volume-averaged density scales roughly as χ^-1 = M_A^-α, i.e., M_A^-4 (sub-Alfvénic) and M_A^-3 (super-Alfvénic); the additional steepening to M_A^-6 is unexplained and may be an artifact of the inverted effective radius in §4.1. As it stands, the abstract's statement that the M_A^-6 law is independent of energy and radiative losses is not supported.","section":"§4.4, Fig. 8(a)"}],"minor_comments":[{"comment":"The code name is spelled 'CRIPTIC' in the main text but 'CRIPPTIC' in Appendix A; please standardize the spelling.","section":"§3, Appendix A"},{"comment":"There are grammatical slips such as 'from which can see that' (Sect. 4.1.1) and 'fast/normal-diffusion diffusion' (Sect. 4.3); these should be corrected.","section":"§4.1.1, §4.3"},{"comment":"The footnote justifies the coordinate change by reference to a standard Green's function solution but does not show the transformation; please provide an explicit derivation from Eq. (1), since the current definition is directly related to the major issue above.","section":"§4.1, footnote 3"},{"comment":"Please specify exactly how UCR is averaged (over which volume or radial range, and at which time) and report fit uncertainties for the power-law indices β and ε; the error bars shown are standard deviations of UCR, not uncertainties on the fitted slopes.","section":"§4.4, Fig. 8(a)"},{"comment":"When presenting Eq. (3), the text should distinguish more clearly between theoretically predicted scalings and numerically verified scalings from the cited test-particle simulations, so that the input status of this relation in the present work is transparent.","section":"§2.2"}],"recommendation":"reject","confidential_remarks":"The central claims are circular (the D_⊥/D_∥ vs. M_A scaling is an input) and the UCR–M_A result is entangled with an inverted effective-radius transformation. I do not see a revision within the current scope that would establish the headline 'intrinsic property' and UCR ∝ M_A^-6 results; establishing them would require either evolving the MHD turbulence or reframing the paper as a parameter study of a prescribed anisotropic diffusion model, with the UCR–M_A and 'intrinsic property' claims removed or heavily qualified. A significantly revised version focusing on the loss-modified spectral evolution and sub-diffusion behavior could be reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the central results are largely predetermined by the input assumptions. Appendix A Step 2 hard-wires chi = M_A^alpha into CRIPTIC, with alpha = 4 sub-Alfvenic and 3 super-Alfvenic. Recovering D_perp/D_parallel ~ M_A^alpha in Fig. 8(c) is therefore a consistency check, not a numerical discovery. The UCR ~ M_A^-6 relation also follows directly from that input: for alpha = 4, the volume of the rescaled isotropic diffusion region scales as (chi D_parallel)^{3/2}, giving UCR ~ chi^{-3/2} ~ M_A^{-6}. So that claim is not independent either.\n\nThere is a more concrete technical problem. The effective radius in Sec. 4.1 is defined as r^2 = x^2/chi + y^2 + z^2. For D_perp = chi D_parallel, the coordinate change that makes the Fokker-Planck operator isotropic is either X = sqrt(chi) x (giving R_eff^2 = chi x^2 + y^2 + z^2) or Y = y/sqrt(chi), Z = z/sqrt(chi) (giving x^2 + (y^2+z^2)/chi). The paper rescales the parallel coordinate by 1/sqrt(chi), which is the inverse of either correct transformation. So the radial profiles (r^-6/5, r^-11/5) and the averaged UCR(M_A) are computed in a distorted metric. This undermines the quantitative power-law claims.\n\nWhat the paper does well: it is transparent about the input scalings and reviews the theoretical background clearly. The qualitative effects of radiative losses—steeper radial profile and sub-diffusion in the PeV band—are sensible. The morphology maps showing viewing-angle dependence and the transition from sub- to super-Alfvenic anisotropy are consistent with earlier anisotropic diffusion work and could serve as useful illustrations.\n\nThe soft spots are serious, but the paper is not incoherent. The radiative-loss and morphology results stand on their own, and the error in the effective radius is fixable. The main issue is framing: the authors present properties of their input model as emergent findings. Also, because the turbulent field is never evolved self-consistently, the paper cannot test MHD scaling laws; it only probes the consequences of assuming them.\n\nWho is this for? Someone wanting a worked parameter study of anisotropic diffusion with radiative losses, or a classroom example of why input assumptions matter. The quantitative results should not be cited until the effective radius is corrected and the circularity claims are reframed.\n\nRecommendation: send to peer review, not desk reject. A good referee will catch the inverted radius and force a careful reanalysis and reframing. The paper has enough substance to be salvageable.","headline":"The paper's headline scalings are mostly hard-wired inputs, and the effective radius transformation appears inverted, so the main discovery claims do not survive scrutiny; still, the parameter study is substantive enough to warrant a referee rather than a desk reject.","tokens_in":19048,"tokens_out":5783,"would_cite":false,"duration_ms":53679,"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":"In a long numerical evolution of cosmic-ray electrons in MHD turbulence, the paper finds that the anisotropic diffusion scalings persist and the energy density falls as $M_A^{-6}$.","keywords":["cosmic-ray transport","anisotropic diffusion","magnetohydrodynamic (MHD) turbulence","Alfvén Mach number","pulsar wind nebula","TeV halos","radiative losses","Fokker-Planck transport"],"falsifier":"Run a collisionless test-particle simulation in which the MHD turbulent field is evolved self-consistently, measuring $D_\\parallel$ and $D_\\perp$ for 1 GeV to 10 PeV electrons in sub- and super-Alfvénic regimes with radiative cooling; if the recovered exponent of $D_\\perp/D_\\parallel$ vs $M_A$ departs systematically from $\\alpha=4,3$, or the recovered $U_{CR}$-vs-$M_A$ slope departs from $-6$, the claimed universality fails. Observationally, compare halo brightness around pulsars with independently estimated $B_0$ and check whether the inferred $M_A$ agrees with the measured magnetization.","tokens_in":17889,"feed_emoji":"🌌","tokens_out":11055,"duration_ms":99234,"temperature":0.7,"pith_summary":"The paper tries to establish that anisotropic diffusion—fast along a mean magnetic field, slow across it—is an intrinsic, long-lived description of high-energy cosmic-ray transport near pulsar wind nebulae. Using a numerical Fokker-Planck transport simulation with sub-grid diffusion coefficients, the authors evolve injected electron populations in a turbulent magnetized medium, treating the magnetization parameter $M_A$ and the turbulence spectral index $\\gamma$ as free parameters. They report that the two input diffusion laws remain intact after a 3 kyr evolution: $D_\\parallel \\propto E^{2-\\gamma}$ and $D_\\perp/D_\\parallel \\propto M_A^\\alpha$ with $\\alpha=4$ sub-Alfvénic and $\\alpha=3$ super-Alfvénic. They also find that the volume-averaged cosmic-ray energy density obeys $U_{CR} \\propto M_A^{-6.0\\pm0.5}$, independent of particle energy and radiative losses, while losses convert normal PeV diffusion into sub-diffusion. If these results hold, the morphology and brightness of TeV halos become observable probes of the magnetization and turbulence spectrum of the surrounding interstellar medium.","feed_headline":"Cosmic-ray halo density falls as M_A to the -6th power","feed_subtitle":"Long-running simulations show the scaling survives energy changes and radiative losses, linking halo shape to turbulence.","key_machinery":"The machinery is the ratio law $D_\\perp/D_\\parallel \\approx M_A^\\alpha$, with $\\alpha=4$ for sub-Alfvénic and $\\alpha=3$ for super-Alfvénic turbulence, joined to the energy law $D_\\parallel = D_0(E/m_e c^2)^\\delta$ with $\\delta=2-\\gamma$. The paper feeds these laws into a Fokker-Planck transport equation through an anisotropic effective radius $r=\\sqrt{(x/\\sqrt{\\chi})^2+y^2+z^2}$ with $\\chi=M_A^\\alpha$, evolves injected electrons for 3 kyr under synchrotron and inverse-Compton cooling, and then re-extracts $D_\\parallel$ and $D_\\perp$ from ensemble-averaged squared displacements. The argument stands on the fact that the returned exponents match the input exponents, and that the scan over $M_A$ yields an emergent $U_{CR}\\propto M_A^{-6.0\\pm0.5}$.","core_discovery":"The authors' central claim is that the anisotropic-diffusion relation $D_\\perp/D_\\parallel \\sim M_A^\\alpha$ and the energy scaling $D_\\parallel \\sim E^{2-\\gamma}$ are preserved through a long cosmic-ray evolution, so they can be treated as intrinsic properties of transport in turbulent magnetic fields. On top of that, they report a new volume-averaged power law $U_{CR} \\sim M_A^{-6.0\\pm0.5}$ that is independent of particle energy and of radiative losses, with the spatial profile following $U_{CR}\\propto r^{-6/5}$ near the source and steepening to $r^{-11/5}$ when PeV electrons cool, while radiative cooling turns PeV transport from normal diffusion into sub-diffusion with $d^2\\propto t^{4/5}$. They further show that halo morphology is oval when viewed across the mean field and round when viewed along it, with the anisotropy elongation controlled by $M_A$ through the exponent $\\eta=M_A^{\\alpha/2}$. Thus the paper converts TeV halo shape and brightness into readouts of turbulence magnetization and spectral index.","pith_inferences":["A self-consistent test-particle simulation that actually evolves the turbulent field, rather than prescribing the diffusion scalings, would show whether $\\alpha=4,3$ and the $M_A^{-6}$ law are emergent or merely inherited from the input assumptions; the paper itself does not evolve the turbulence.","The energy independence of the $M_A$ scaling suggests that the magnetization could be inferred from ratios of halo brightness at two energies rather than from absolute fluxes, reducing model dependence in future observations.","Because the study uses electrons and only synchrotron and inverse-Compton cooling, whether the $M_A^{-6}$ correlation survives for protons with hadronic losses is a testable next step the paper leaves open.","Observational maps of several pulsar halos with independently estimated magnetic-field strengths could test the predicted $U_{CR}$-versus-$M_A$ slope directly, since $M_A$ enters through the mean field $B_0$."],"forward_implications":["The recovered $D_\\perp/D_\\parallel \\sim M_A^{4}$ sub-Alfvénic and $\\sim M_A^{3}$ super-Alfvénic relations, independent of energy and radiative losses, justify using these laws in interpretive models of TeV halos.","The averaged energy density $U_{CR}\\propto M_A^{-6.0\\pm0.5}$ gives a quantitative link between halo brightness and the Alfvén Mach number of the surrounding medium.","The near-source spatial profile $U_{CR}\\propto r^{-6/5}$, steepened to $r^{-11/5}$ by losses, provides a template for fitting observed halo surface-brightness profiles.","Radiative losses convert PeV transport from normal diffusion to sub-diffusion with $d^2\\propto t^{4/5}$, so high-energy halo morphology encodes cooling as well as turbulence.","The anisotropy exponent $\\eta=M_A^{\\alpha/2}$ connects the oval aspect ratio of a halo to the magnetization, allowing viewing angle and turbulence properties to be separated with multi-angle observations."],"supporting_citations":[{"why":"Supplies the transport code and the analytic Green's-function solutions used to evolve the cosmic-ray distributions in the three simulation groups.","marker":"Krumholz et al. 2022"},{"why":"Provides the analytical prediction $D_\\perp/D_\\parallel\\approx M_A^\\alpha$ with $\\alpha=4$ and $\\alpha=3$ that the simulations inherit and test.","marker":"Lazarian & Yan 2014"},{"why":"Establishes the $M_A^4$ dependence for perpendicular diffusion in sub-Alfvénic turbulence, the exponent adopted for the sub-Alfvénic regime.","marker":"Yan & Lazarian 2008"},{"why":"Numerical test-particle simulation supporting the sub-Alfvénic perpendicular-diffusion scaling used as input.","marker":"Xu & Yan 2013"},{"why":"Numerical confirmation of the predicted $\\alpha=4$ and $\\alpha=3$ scalings in sub- and super-Alfvénic regimes.","marker":"Maiti et al. 2022"},{"why":"Reports the TeV halos around Geminga and Monogem with a slow diffusion coefficient, the observation the anisotropic-diffusion model is built to explain.","marker":"Abeysekara et al. 2017"},{"why":"Prior model explaining slow diffusion with anisotropic diffusion and small Alfvén Mach number and viewing angle; the current paper widens this parameter space.","marker":"Liu et al. 2019"},{"why":"Constrains anisotropic diffusion models to $M_A\\lesssim0.3$ and viewing angle $\\psi\\lesssim5^\\circ$; the paper tests whether slow diffusion can occur outside that narrow space.","marker":"De La Torre Luque et al. 2022"},{"why":"Source of the empirical energy-dependent parallel diffusion relation $D_\\parallel = D_0(E/m_e c^2)^\\delta$ used as the baseline law.","marker":"Seo & Ptuskin 1994"}],"fun_headline_variants":["Cosmic-ray halo density follows M_A^-6 law independent of losses","Magnetization alone sets cosmic-ray halo density as -6 power","Halo density falls with magnetization to the -6, energy-independent","Turbulent magnetization imprints -6 scaling on cosmic-ray halo density","Cosmic-ray halos: density tied to turbulence via -6 power law"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two input scalings—parallel diffusion growing as energy to a power and the perpendicular-to-parallel ratio growing as the fourth or third power of the magnetization—are correct descriptions of real turbulent transport, since the simulations impose those scalings as sub-grid inputs and never evolve the turbulent magnetic field itself.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic-ray halo density follows M_A^-6 law independent of losses","Magnetization alone sets cosmic-ray halo density as -6 power","Halo density falls with magnetization to the -6, energy-independent","Turbulent magnetization imprints -6 scaling on cosmic-ray halo density","Cosmic-ray halos: density tied to turbulence via -6 power law"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001048,"raw_usage":{"total_tokens":4458,"prompt_tokens":1057,"completion_tokens":3401,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":673,"completion_tokens_details":{"reasoning_tokens":3306}},"tokens_in":673,"tokens_out":3401,"duration_ms":24842,"temperature":1.0,"reasoning_tokens":3306,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:22:02.401586+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a collisionless test-particle simulation in which the MHD turbulent field is evolved self-consistently, measuring $D_\\parallel$ and $D_\\perp$ for 1 GeV to 10 PeV electrons in sub- and super-Alfvénic regimes with radiative cooling; if the recovered exponent of $D_\\perp/D_\\parallel$ vs $M_A$ departs systematically from $\\alpha=4,3$, or the recovered $U_{CR}$-vs-$M_A$ slope departs from $-6$, the claimed universality fails. Observationally, compare halo brightness around pulsars with independently estimated $B_0$ and check whether the inferred $M_A$ agrees with the measured magnetization.","supporting_citations":[{"cited_title":"& Yan , H","cited_arxiv_id":null,"evidence_quote":"Provides the analytical prediction $D_\\perp/D_\\parallel\\approx M_A^\\alpha$ with $\\alpha=4$ and $\\alpha=3$ that the simulations inherit and test."},{"cited_title":"& Lazarian , A","cited_arxiv_id":null,"evidence_quote":"Establishes the $M_A^4$ dependence for perpendicular diffusion in sub-Alfvénic turbulence, the exponent adopted for the sub-Alfvénic regime."},{"cited_title":"2022, , 926, 94","cited_arxiv_id":null,"evidence_quote":"Numerical confirmation of the predicted $\\alpha=4$ and $\\alpha=3$ scalings in sub- and super-Alfvénic regimes."},{"cited_title":"2019, , 123, 221103","cited_arxiv_id":null,"evidence_quote":"Prior model explaining slow diffusion with anisotropic diffusion and small Alfvén Mach number and viewing angle; the current paper widens this parameter space."},{"cited_title":"2022, , 106, 123033","cited_arxiv_id":null,"evidence_quote":"Constrains anisotropic diffusion models to $M_A\\lesssim0.3$ and viewing angle $\\psi\\lesssim5^\\circ$; the paper tests whether slow diffusion can occur outside that narrow space."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the empirical energy-dependent parallel diffusion relation $D_\\parallel = D_0(E/m_e c^2)^\\delta$ used as the baseline law."}],"review_version":1}