{"id":"de81761a-5ffa-4ab1-b66f-80a5d79da854","arxiv_id":"2507.22107","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Young stars in Orion, Upper Sco, Taurus, and Perseus carry kinematic imprints of turbulence, supernova energy injection, anisotropic expansion, and magnetic fields.","lead":"This paper uses full 3D positions and 3D velocities of young stars in four nearby star-forming regions to read the leftover imprints of turbulence, supernova explosions, gravity, and magnetic fields. It finds anisotropic motions and supernova-like features in several groups, and in Perseus the stars move mostly perpendicular to the local magnetic field.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"VSF bumps attributed to supernovae are not robust to substructure and group-size artifacts; the Taurus+Perseus ~100 pc peak likely reflects bulk relative motion between the two groups rather than an expanding shell.","rationale":"The reader's conditional verdict already identifies the VSF-bump-to-supernova interpretation as the weakest assumption, and our read agrees. The concern is real but does not overturn the paper: the central claim is supported by independent lines of evidence, including anisotropic expansion profiles in Orion BCD, λ Ori, and Upper Sco (Section 4.2), and the magnetic-field alignment in Perseus (Section 4.3, p ≈ 0.01 against a Monte Carlo null). The paper is transparent about many limitations (Section 5) and explicitly flags the Perseus and Upper Sco peaks as possible artifacts. Our added emphasis is that the Taurus+Perseus ~100 pc peak is even more likely to be a cross-group bulk-motion artifact than a physical shell, and that a null test is straightforward. This does not change the CONDITIONAL verdict; it strengthens the case for requesting such a test in revision. Credit is due for the use of two independent magnetic-field maps, a Monte Carlo null for the field alignment, and consistency with prior expansion measurements—these are genuine checks that support parts of the central claim even if the SN-bump evidence is weakened.","tokens_in":21658,"tokens_out":6618,"duration_ms":82106,"concrete_test":"Recompute the Taurus+Perseus 6D VSF after subtracting each group's median velocity vector from its stars before computing pair velocity differences. If the ~100 pc peak disappears, it is a bulk-motion artifact of combining two groups, not a supernova shell. Repeat this subgroup-mean subtraction for the λ Ori and Perseus VSFs to see whether their peaks persist; persisting peaks would support the shell interpretation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim's supernova-injection evidence leans on interpreting VSF bumps as SN shells (Section 4.2.2). This is insecure because several bumps occur at scales where the VSF is dominated by group geometry rather than a physical shell. The paper concedes for Perseus that the ~80 pc peak 'could instead be an artifact, due to this ℓ approaching the size of the group' and for Upper Sco that the peak sits at the largest probed ℓ. Most tellingly, the combined Taurus+Perseus VSF (Figure 5) has a ~100 pc peak that is expected from cross-pairs between two spatially separated groups moving with different bulk velocities; the resulting velocity offset appears as a bump at the inter-group separation and yields a spurious age estimate of ℓ/⟨|δv|⟩. No null test is provided to distinguish an expanding shell (radial velocity differences aligned with pair separations) from substructure or bulk relative motion. If these bumps are artifacts, the 'multiple generations of supernovae' conclusion loses much of its support, although the anisotropic expansion and magnetic-field results remain.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses Gaia DR3 astrometry and APOGEE DR17 radial velocities to construct 6D kinematics of young stars in Orion (split into Orion BCD, ONC, Orion A, and lambda Ori), Upper Sco, Taurus, and Perseus. It computes 6D and 4D first-order velocity structure functions, radial expansion profiles from iteratively defined group centers, and compares stellar proper motions to plane-of-sky magnetic field orientations from Planck and the VGT. The main claims are that the VSFs are broadly consistent with Larson's relation at intermediate scales, that several VSF bumps indicate local energy injection from supernovae, that several groups show clear anisotropic expansion, and that Perseus stars move preferentially perpendicular to the local magnetic field. The paper concludes that young stars retain kinematic memory of turbulence, gravity, supernovae, and magnetic fields, i.e., more than just the turbulent state of their natal clouds.","tokens_in":21814,"tokens_out":5166,"duration_ms":61857,"significance":"If the VSF-bump interpretations hold, the paper would provide a notable observational demonstration of multiple physical processes imprinted in the 6D kinematics of young stars, extending the earlier Ha21/Ha22 analyses. The study has clear strengths: the supernova age estimates use a parameter-free scaling (ell/<|delta v|>), uncertainties are propagated through Monte Carlo resampling and MCMC fits, two independent magnetic field maps are used for the alignment analysis, and the authors explicitly discuss several limitations and biases. The anisotropic expansion profiles and the Perseus magnetic-field alignment result are potentially valuable even if the supernova interpretation of some VSF bumps is weakened by geometric artifacts.","major_comments":[{"comment":"The combined Taurus+Perseus VSF is computed from LSR-frame velocities without subtracting each group's bulk velocity. Cross-pairs between the two spatially separated groups contribute at separations comparable to the inter-group distance, so a bump at ell ~ 100 pc with <|delta v|> ~ 20 km/s can be produced by the systemic velocity difference alone, independent of any expanding shell. The age estimate of ~5 Myr from ell/<|delta v|> is therefore not evidence for the Per-Tau shell unless a null test (e.g., randomizing velocities within each group, or subtracting group mean velocities before computing the VSF) shows that the bump survives. This is load-bearing for the conclusion that 'multiple generations of supernovae and their effect on stellar kinematics are observed in the combined VSFs of Taurus and Perseus.'","section":"Section 4.2.2 and Figure 5"},{"comment":"The paper itself states that the Perseus ~80 pc peak 'could instead be an artifact, due to this ell approaching the size of the group' and that the Upper Sco peak sits at the largest probed ell; these are precisely the scales where sparse pair counts and group-edge effects dominate the VSF. The derived supernova ages of 4-6 Myr for Perseus and ~4 Myr for Upper Sco are therefore not robust, and these cases should either be explicitly excluded from the supernova-injection claim or supported by a jackknife or edge-removal test. Without such a test, the robust supernova evidence reduces to Orion BCD and lambda Ori only.","section":"Section 4.2.2, Perseus and Upper Sco"},{"comment":"The rescaling of proper-motion errors to the radial-velocity error width is ad hoc, and the paper later concedes that it cannot determine whether the heightened x-direction 4D VSF amplitudes are physical or due to the intrinsic line-of-sight velocity scatter. Because the 4D VSF amplitude comparison underpins the 'mild anisotropy' statement, the systematic uncertainty from propagating the reported errors with and without the rescaling should be quantified, and the anisotropy claim should be explicitly restricted to slopes and features if the amplitude information cannot be separated from the measurement noise.","section":"Section 3.1 and Section 5"}],"minor_comments":[{"comment":"There is a typo: 'we preform random sampling' should be 'we perform random sampling'; also in Section 4.1 'Komolgorov' should be 'Kolmogorov'.","section":"Section 3.1"},{"comment":"'Plank Polarization' should be 'Planck Polarization'.","section":"Figure 2 caption"},{"comment":"The cut described as 'proper motion > 60 mas/yr' needs clarification: as written it removes high proper-motion stars, which appears to contradict the intended purpose of removing contaminants; an absolute-value or inequality direction should be stated explicitly.","section":"Section 2.1"},{"comment":"The 'factor of 1.023 for unit conversion' used when inverting expansion slopes is not defined; please state the unit conversion (e.g., 1 pc/Myr = 1.023 km/s) explicitly.","section":"Section 4.2.2"},{"comment":"The radial-velocity error cut is described in the text as 1.5 times the MAD, but Figure 7 shows a fixed red dashed line; clarify that the plotted line is the group-specific 1.5 MAD threshold.","section":"Section 2.1 and Figure 7"}],"recommendation":"major_revision","confidential_remarks":"The central uncertainty is whether the Taurus+Perseus VSF peak and some other bumps survive a bulk-motion/geometry null test; if they do not, the 'multiple generations of supernovae' conclusion should be substantially softened. The paper's breadth is a strength, but the current presentation overstates the robustness of several supernova interpretations relative to the admitted group-size and sampling artifacts."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, here is my read on Velguth et al. (arXiv:2507.22107). It is a solid follow-up to Ha21/Ha22 on the same star-forming regions, with three genuinely new pieces: the 4D VSFs split by galactic Cartesian direction, the 3D expansion profiles from MCMC fits, and the proper-motion vs magnetic-field alignment in Taurus and Perseus. The last of these is the most interesting. The Perseus result (p ~ 0.01, preferential perpendicular motion) is a first, and the analysis looks careful, with two independent field maps (Planck and VGT) and a random-position null test. The anisotropic expansion in lambda Ori and Orion BCD also comes through cleanly and matches earlier work.\n\nThe soft spot is the supernova interpretation of VSF bumps. The paper itself concedes for Perseus that the ~80 pc peak could be an artifact from approaching the group size, and for Upper Sco the peak sits at the largest probed ell. The combined Taurus+Perseus VSF (Figure 5) is even shakier: the ~100 pc peak is exactly what you get from cross-pairs between two spatially separated groups with different bulk velocities. The paper mentions under-sampling and the dip at the peak but does not test the null that this is just group separation rather than a shell. So the conclusion that 'multiple generations of supernovae are observed in the combined VSFs' overreaches. The Orion and lambda Ori bumps are more plausible, but even there, substructure can produce peaks.\n\nAlso, the proper-motion error scaling in Section 3.1 is ad hoc, and the paper admits it cannot distinguish a physical x-direction VSF amplitude from line-of-sight measurement scatter. That is honest, but it means the 4D anisotropy result rests on slope and feature differences, not amplitudes, which is fine where slopes differ (lambda Ori, Taurus).\n\nOverall: the central idea that young stars retain kinematic imprints beyond turbulence is plausible and supported by the expansion and B-field lines. The SN-injection component is the least secure. The paper deserves peer review, and I would send it out, but the SN interpretation needs a null test against substructure and relative motion before publication. The B-field and expansion work alone justify referee time.","headline":"A useful follow-up with a strong B-field result and clean anisotropic expansion profiles, but the supernova-sourced VSF bumps need a substructure null test before they carry the weight the conclusions place on them.","tokens_in":22468,"tokens_out":2418,"would_cite":true,"duration_ms":26760,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Young stars in nearby star-forming regions retain kinematic memories of turbulence, gravity, supernova explosions, and magnetic fields, so their 6D motions can probe processes gas observations alone cannot.","keywords":["velocity structure functions","young stellar kinematics","star-forming regions","ISM turbulence","supernova feedback","magnetic field alignment","Gaia DR3 astrometry","Larson's relation"],"falsifier":"A decisive test would be to build null catalogs with the same number of stars, the same spatial distribution, and the same velocity noise but no expanding shell, run the identical VSF binning, and check whether bumps at ~60 pc, ~100 pc, and ~20 pc survive; if they do, the supernova-injection reading loses its evidence. Alternatively, a larger sample that fills in the poorly sampled large-$\\ell$ bins (for example, from future Gaia data releases) should preserve the bumps if they are physical shells.","tokens_in":21384,"feed_emoji":"🌌","tokens_out":11403,"duration_ms":113071,"temperature":0.7,"pith_summary":"The paper asks what young stars remember about the gas they were born from, using Gaia DR3 astrometry and APOGEE radial velocities to build full 6D position-plus-velocity maps of four nearby star-forming regions: Orion, Upper Sco, Taurus, and Perseus. It argues that the stars' kinematics are not a pure turbulent cascade: first-order velocity structure functions — the mean velocity difference between star pairs versus their separation — show bumps that the authors read as supernova-injected shells, radial profiles reveal anisotropic expansion from gravity and shell dynamics, and in Perseus the stars' proper motions are preferentially perpendicular to the local magnetic field. The central claim is that young stars retain the imprints of multiple overlapping physical processes rather than just the turbulent state of their natal clouds, so their 6D motions can serve as tracers of the interstellar medium in ways gas observations alone cannot. If right, this gives a new way to date past supernovae from the scale of the velocity bumps and to connect star formation to cloud-scale dynamics.","feed_headline":"Young stars remember supernovae, turbulence, and magnetic fields","feed_subtitle":"Stellar motions in four nearby clouds reveal supernova shells, anisotropic expansion, and magnetic alignment","key_machinery":"The central tool is the first-order velocity structure function (VSF), $\\langle|\\delta v|\\rangle(\\ell)$, the mean absolute velocity difference between pairs of young stars as a function of their three-dimensional separation $\\ell$. Different physical processes predict different slopes and features: Kolmogorov turbulence gives $\\langle|\\delta v|\\rangle \\propto \\ell^{1/3}$, Burgers-like compressible turbulence gives $\\ell^{1/2}$, free expansion gives $\\ell^1$, and a relaxed cluster gives slope 0. Computing the VSF in 4D (one velocity component at a time) as well as 6D isolates anisotropic motion, while radial expansion profiles (outward velocity versus radius from the group center, fit with MCMC linear models) separate bulk expansion from turbulent scaling. The magnetic-field analysis compares each star's proper-motion position angle with the local plane-of-sky magnetic field orientation from Planck polarization and from the velocity gradient technique, testing the resulting $\\delta\\theta$ distribution against a random distribution with a Kolmogorov-Smirnov test.","core_discovery":"Working in galactic Cartesian coordinates, the paper computes 6D and 4D first-order velocity structure functions (VSFs) for Orion (split into λ Ori, ONC, Orion A, and Orion BCD), Upper Sco, Taurus, and Perseus. At intermediate separations (roughly 10 to 100 pc) the VSFs mostly follow Larson's relation, showing that a turbulent signal is retained; the anisotropy between the x, y, and z components is mild except in λ Ori and Taurus. Several VSFs show bumps or turnovers that the authors interpret as local energy injection from supernovae: a ~60 pc peak in λ Ori (~6 Myr), a ~70 pc bump in Orion (~7 Myr), a ~30 pc turnover in Upper Sco that is contaminated by older stars and poorly sampled, a ~20 pc peak in Taurus, and a ~100 pc peak in the combined Taurus and Perseus sample (~5 Myr, consistent with the Per-Tau shell from Bialy et al. 2021). A ~80 pc peak in Perseus is flagged by the paper itself as possibly an artifact of the group's size. Radial velocity profiles show clear anisotropic expansion in Orion BCD, λ Ori, and Taurus, which the authors attribute to a mix of gravitational collapse, subgroup motions, and expanding supernova shells. In Perseus, the proper motions of young stars are preferentially perpendicular to the plane-of-sky magnetic field measured by both Planck polarization and the velocity gradient technique, while Taurus shows no significant alignment; the authors hypothesize that cloud collapse dragged the field into this configuration. The paper concludes that each group shows multiple, overlapping causes and that young stars retain strong kinematic memories of the conditions before their formation.","pith_inferences":["Beyond the paper: if the VSF bumps are genuine shells, age-binning the same groups should make the young-population peaks sharper and reveal the underlying Kolmogorov slope in the newest stars; the paper notes this as future work.","Beyond the paper: the perpendicular motion in Perseus could be checked against high-resolution field maps in other filaments; a prediction is that YSO proper motions in strongly magnetized, sheet-collapsed clouds will also run perpendicular to the field.","Beyond the paper: the combined Taurus-Perseus VSF peak at ~100 pc offers a kinematic age for the Per-Tau shell that could be compared with stellar traceback; tracing the 3D velocities of shell members back in time should place the expansion center near the proposed shell.","Beyond the paper: a possible artifact test is to inject synthetic expanding shells into mock catalogs with the same group sizes and membership selection and ask whether the observed bump amplitudes and locations are recovered; this would calibrate how much of the peaks could be sampling geometry."],"forward_implications":["If the VSF-bump interpretation is right, supernova shell ages can be estimated directly from young-star kinematics via $\\ell/\\langle|\\delta v|\\rangle$, giving a stellar-based complement to gas-based shell dating.","The anisotropic expansion seen in Orion BCD, λ Ori, and Taurus means that bulk motions, not just turbulence, must be included when interpreting velocity statistics of young associations.","The Perseus perpendicular alignment, if real, implies that magnetic fields shape the kinematics stars inherit in at least some clouds, making stellar proper motions a new observational probe of magnetic cloud structure.","The combination of turbulence, gravity, supernovae, and field effects in every group implies that no single physical process can be inferred from a group's velocity structure alone; multi-probe analyses are required.","Because the paper recovers Larson-consistent scaling at intermediate scales, young stars can be used as 3D tracers of ISM turbulence in regions where gas line-of-sight velocities are the only gas tracer."],"supporting_citations":[{"why":"Introduced the VSF method for young stars and reported Kolmogorov-like scaling, providing the baseline this paper extends with 6D and 4D analyses.","marker":"Ha et al. (2021, 2022)"},{"why":"Establishes the size-velocity dispersion relation used as the turbulent baseline against which the observed VSFs are compared.","marker":"Larson (1981)"},{"why":"Provides the $\\ell^{1/3}$ velocity scaling that defines the turbulent signature the paper seeks to identify and deviate from.","marker":"Kolmogorov (1941)"},{"why":"Supplies evidence that Orion subgroups are expanding and that a supernova may have driven the expansion, anchoring the Orion BCD and λ Ori interpretations.","marker":"Kounkel et al. (2018)"},{"why":"Gives the ~6 Myr age estimate for the Orion supernova that the paper's ~7 Myr VSF-based age is compared with.","marker":"Kounkel (2020)"},{"why":"Identifies the Per-Tau shell and Tau Ring that the combined Taurus-Perseus VSF and the Taurus ~20 pc peak are interpreted as signatures of.","marker":"Bialy et al. (2021)"},{"why":"Reports the Ophiuchus supernova ~2 Myr ago against which the Upper Sco VSF turnover is discussed.","marker":"Neuhäuser et al. (2020)"},{"why":"Provides the expansion-profile methodology and the λ Ori expansion measurement that the paper's anisotropic expansion results are compared with.","marker":"Armstrong & Tan (2024)"},{"why":"Shows the 3D magnetic field in Perseus is perpendicular to the star-forming clouds, supporting the interpretation of the observed proper-motion alignment.","marker":"Tahani et al. (2022)"},{"why":"Provides the polarized dust emission maps used as one of the two magnetic field orientation tracers.","marker":"Planck Collaboration et al. (2020)"}],"fun_headline_variants":["6D stellar motions expose supernova energy injection and magnetic alignment","Young stars retain kinematic memories of turbulence and supernova shocks","Nearby star motions reveal supernova shells and magnetic field effects","Stellar kinematics in four clouds show turbulent and supernova signatures","Gaia and APOGEE reveal supernova and magnetic imprints in young stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim depends on treating a bump or turnover in the velocity structure function at a given separation as the imprint of an expanding supernova shell, rather than as the group's own internal substructure or the under-sampling of separations nearly as large as the group.","fun_headline_variants_meta":{"raw":{"variants":["6D stellar motions expose supernova energy injection and magnetic alignment","Young stars retain kinematic memories of turbulence and supernova shocks","Nearby star motions reveal supernova shells and magnetic field effects","Stellar kinematics in four clouds show turbulent and supernova signatures","Gaia and APOGEE reveal supernova and magnetic imprints in young stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00065,"raw_usage":{"total_tokens":3115,"prompt_tokens":1210,"completion_tokens":1905,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":826,"completion_tokens_details":{"reasoning_tokens":1816}},"tokens_in":826,"tokens_out":1905,"duration_ms":15860,"temperature":1.0,"reasoning_tokens":1816,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T12:03:26.341855+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be to build null catalogs with the same number of stars, the same spatial distribution, and the same velocity noise but no expanding shell, run the identical VSF binning, and check whether bumps at ~60 pc, ~100 pc, and ~20 pc survive; if they do, the supernova-injection reading loses its evidence. Alternatively, a larger sample that fills in the poorly sampled large-$\\ell$ bins (for example, from future Gaia data releases) should preserve the bumps if they are physical shells.","supporting_citations":[{"cited_title":"Expansion Kinematics of Young Clusters. I. Lambda Ori","cited_arxiv_id":"2407.11845","evidence_quote":"Provides the expansion-profile methodology and the λ Ori expansion measurement that the paper's anisotropic expansion results are compared with."}],"review_version":1}