{"id":"8a1e4643-5423-4c42-ad67-458bfaa45017","arxiv_id":"2501.00506","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"ALMA N2H+ data confirm G11P1-HFS as a small-scale nascent hub-filament system whose gas flows are aligned with gravity along the filaments.","lead":"Using ALMA and JWST observations, this paper confirms that the massive protostar G11P1 sits at the hub of a small, young filament network, and it measures gas flowing inward along the filaments. It offers a rare, relatively feedback-free case for studying how high-mass stars assemble their mass from their surroundings.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The gravity-driven inflow claim rests on a Vg-Fg alignment whose mean angle (41°) is nearly the random expectation (45°); no significance test is reported, so the alignment may be chance.","rationale":"The reader's verdict was CONDITIONAL, and I agree that the paper should be accepted with revisions rather than rejected. However, I identify a different load-bearing weakness than the reader's chosen assumption of N2H+ abundance uncertainty. The most direct threat to the central claim is the statistical significance of the Vg-Fg alignment. The paper's own quantitative statement—mean/median Δθ_{v,g} of 41°/39°—is only slightly below the 45° expected from a uniform random distribution. No significance test, confidence interval, or null-model comparison is provided. Without such a test, the phrase 'Vg and Fg align along the filaments, indicating gravity-driven flows' is an overstatement. This weakness is more fundamental than the abundance issue because it attacks the evidence for gravity-driven inflow even under ideal column-density maps. If the alignment is not significant, the main new conclusion of the paper is unsupported, regardless of whether the column density is accurate to 60% or 5%. The abundance uncertainty, while real and worth propagating, would only change the Fg directions if the N2H+ abundance varies spatially; a constant abundance error would not affect the alignment angle. Thus the statistical test is the first-order check that should be required. The paper does have other merits: it convincingly confirms the small-scale HFS morphology with ALMA N2H+, shows a clear V-shaped velocity profile across the hub, and provides plausible kinematic evidence for inflow. These justify a conditional acceptance with a requested significance test and a toned-down interpretation of the alignment. I therefore keep the verdict CONDITIONAL (UNCHANGED) and partially agree with the reader's framing, since they noticed the 41° value but did not elevate it to the primary concern.","tokens_in":13565,"tokens_out":7279,"duration_ms":72561,"concrete_test":"Perform a Monte Carlo test: randomize the pixel orientations of the V_g vector field (or shuffle the centroid-velocity map) while preserving the spatial masks and vector magnitudes, recompute the mean Δθ_{v,g} for each realization, and compare the observed mean (41°) with the resulting null distribution. If >5% of realizations yield a mean ≤41°, the alignment is not significant. Also compute a bootstrap confidence interval on the mean angle; if 45° is inside the 95% CI, the claim 'align' is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central conclusion that 'V_g and F_g align along the filaments, indicating gravity-driven flows' (Abstract, Sec. 4.2) is supported by the distribution of the converging angle difference Δθ_{v,g} in Appendix A (Fig. A.5b). The reported mean and median are approximately 41° and 39°, respectively. For a uniform distribution of angles in [0°, 90°], the expected mean is 45°. The observed mean is only 4° below random, and the paper provides no uncertainty estimate, no bootstrap, and no comparison to a null distribution. Because the vectors are compared after ignoring their direction (angles folded into [0°, 90°]), the effective signal is further diluted. The claim of alignment is used as the primary evidence for gravity-driven mass inflow, so if the alignment is not statistically significant, the headline result is not established. This concern is independent of the column-density abundance uncertainty; even with a perfect N(H2) map, a 41° mean angle is weak evidence. The paper should report a significance test (e.g., Monte Carlo randomization of the V_g field) and/or focus on the filament spines where alignment is claimed to be stronger.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript uses ALMA N2H+(1-0) observations to confirm the presence of the small-scale (<0.6 pc) hub-filament system G11P1-HFS around the high-mass protostar G11P1, previously identified in JWST images. The authors build a position-position-velocity map with SCOUSEPY, measure on-sky velocity gradients, derive N(H2) with XCLASS, and compute sky-projected gravitational force vectors following Wang et al. (2022). They report a V-shaped VLSR profile toward the hub, steep velocity gradients of about 4.7 and -7 km/s/pc, a decreasing Vg toward the hub, and low relative angles between Vg and Fg vectors, from which they conclude that the system is a nascent HFS with gravity-driven inflows along filaments.","tokens_in":13804,"tokens_out":6525,"duration_ms":69981,"significance":"If the gravity-alignment result holds, the paper provides a rare example of a pre-feedback, small-scale HFS and supports hierarchical, clump-fed accretion scenarios. The strengths are the combination of JWST and ALMA data, the explicit PPV-space description of a wiggled funnel, and the comparison of derived quantities with external benchmarks (Zhou et al. 2023; Padoan et al. 2020). The morphological confirmation of the HFS and the V-shaped velocity profile are defensible from the observations. However, the quantitative evidence for gravity-driven inflow is currently weaker than the abstract and conclusions claim: the reported mean alignment angle is close to the random expectation, and the gravitational force vectors inherit an unpropagated abundance uncertainty. These issues are fixable with additional statistical analysis and robustness tests.","major_comments":[{"comment":"The central claim that \"V_g and F_g align along the filaments, indicating gravity-driven flows\" rests on the distribution of Δθ_v,g in Fig. A.5b. The paper reports mean and median values of approximately 41° and 39°. For uniformly distributed angles in [0°, 90°], the expected mean is 45°, so the observed mean is only about 4° from the random expectation. No uncertainty, sample size, or null-hypothesis test is provided, and the histogram in Fig. A.5b is not shown to have a statistically significant excess at small angles. Please report the standard error or a bootstrap confidence interval for the mean, test the distribution against uniformity with an appropriate test for folded angles, and/or restrict the quantitative comparison to the filament spines where Fig. A.5a visually shows low Δθ_v,g values. Without this, the quantitative support for the abstract's alignment claim is not established.","section":"Appendix A, Fig. A.5; abstract and Sec. 4.2"},{"comment":"The gravitational force vectors F_g are computed from the H2 column density map derived with a fixed abundance ratio X(N2H+/H2)=3×10^-10. The paper explicitly states that the derived N(H2) and mass carry 60-70% uncertainties due to potentially lower N2H+ abundance from chemical differentiation and to uncertainties in X. A constant rescaling of N(H2) would change the magnitude but not the direction of F_g, so the V_g-F_g angle would be unaffected by a global abundance normalization. However, the same passage invokes chemical differentiation, which implies that X could vary spatially and thereby change the direction of F_g and hence Δθ_v,g. Please test the robustness of the alignment to plausible spatial abundance variations, for example by recomputing F_g with a varying-X map or by masking low-column-density pixels, and state explicitly that the alignment angle is invariant only to constant scaling of the column density.","section":"Sec. 3.2, Figs. A.3a and A.5a"}],"minor_comments":[{"comment":"The formula for the mass accretion rate is ambiguous: as written, Ẍ = ΔV_obs M / tan(α), where ΔV_obs is described as a velocity gradient (km/s/pc), omits the filament length, and the tangent factor is not clearly defined. Please state the filament length L explicitly and give the full Kirk et al. (2013) expression, or clarify that ΔV_obs is already the velocity gradient per unit length.","section":"Sec. 4.1"},{"comment":"The abstract quotes a velocity gradient of 5 km/s/pc on the southwest side, while Sec. 3.2 reports 4.7 km/s/pc for the same feature; please harmonize the two numbers.","section":"Abstract and Sec. 3.2"},{"comment":"The histogram legend uses the notation |i,g|, |i,v|, and |v,g|, whereas the text uses Δθ_{i,j}; please define the absolute-value and angle-folding conventions in the caption so the reader can relate the two notations.","section":"Fig. A.5"},{"comment":"The text in Sec. 3.2 states that the low relative angle \"hints\" at the supportive role of gravity, while the abstract and Sec. 4.2 state more decisively that V_g and F_g align and indicate gravity-driven flows; please align the strength of the claims with the quantitative evidence presented.","section":"Sec. 3.2 and Sec. 4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is well suited to A&A Letters, and the identification of G11P1-HFS as a small-scale, pre-feedback system is likely publishable on the strength of the morphological matching and the V-shaped velocity profile. The main risk is the V_g-F_g alignment claim, which is currently overstated relative to a mean angle of 41° with no significance test; I recommend major revision rather than rejection because the deficit is in statistical reporting and robustness testing, not in the fundamental data or the HFS identification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The ALMA N2H+ data do what the title promises: they confirm the G11P1 hub-filament system that JWST only saw morphologically, and they add real kinematic content. The V-shaped VLSR profile and the steep gradients (5 and –7 km/s/pc) on either side of the hub are the strongest results. Those alone make the case that G11P1-HFS is a small-scale, actively accreting structure, independent of the gravitational-vector analysis. The paper is also honest about its main uncertainty source – the N2H+ abundance – and it sensibly excludes F6 as likely outflow-contaminated and flags F5 as tentative. That is good observational discipline.\n\nThe soft spot is the Vg–Fg alignment. The paper reports a mean converging angle of about 41° and calls it evidence of gravity-driven inflow. For a uniform distribution folded into [0°, 90°], the random expectation is 45°, so 41° on its own is not convincing. There is no bootstrap, no Monte Carlo null, and no uncertainty on the mean. The spatial statement that alignment is stronger along filament spines is more interesting, but it is made qualitatively and needs a significance test. Also, the 60–70% column-density uncertainty is not propagated into Fg; if the abundance varies spatially rather than just scaling uniformly, the force directions could shift. This is fixable in revision, not a fatal flaw. The \"wiggled funnel\" language is likewise a qualitative visual impression rather than a quantitative morphological classification, and it leans on Zhou et al. (2023) more than it extends it.\n\nThe paper is a single-source case study, so its value is as a benchmark object for simulations and for comparison with other HFSs, not as a statistical sample. The accretion-rate estimates are two orders of magnitude below the Padoan et al. (2020) values, and the authors attribute this plausibly to missing flux, but it does indicate that quantitative numbers here are uncertain. The citation pattern looks appropriate, and I found no circularity: the velocity gradients and gravitational vectors come from independent moments of the same cube but are not fitted to each other, and the conclusions are checked against external results.\n\nVerdict: a serious referee should see this. The central observational result – a confirmed sub-parsec HFS in a nascent stage with clear gas inflow – is solid and worth publishing. But the gravity-driven interpretation needs a proper significance test, and the propagation of abundance uncertainties into Fg should be addressed. With those revisions, this becomes a useful contribution to the clump-fed accretion literature.","headline":"Solid kinematic confirmation of a sub-parsec hub-filament system, but the gravity-alignment claim needs a significance test before it carries the paper.","tokens_in":14417,"tokens_out":2994,"would_cite":true,"duration_ms":31268,"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":"Using ALMA N2H+ observations, this paper confirms that G11P1-HFS is a small-scale, nascent hub-filament system around the high-mass protostar G11P1, and shows that its dense-gas kinematics are consistent with gravity-driven inflow along…","keywords":["hub-filament system","high-mass star formation","dense gas kinematics","ALMA N2H+ observations","velocity gradients","gravitational inflow","PPV analysis","infrared-dark cloud G11.11-0.12"],"falsifier":"Compute H2 column density independently from the ALMA dust continuum at matched resolution and recompute the gravitational force vectors; if the low relative angle between velocity-gradient and gravitational-force vectors along the filaments does not survive a factor-of-two variation in the assumed abundance or the continuum-based mass map, the gravity-driven inflow interpretation would be falsified.","tokens_in":13363,"feed_emoji":"⭐","tokens_out":6874,"duration_ms":58877,"temperature":0.7,"pith_summary":"This paper aims to establish that G11P1-HFS, a hub-filament system around the high-mass protostar G11P1 in the infrared-dark cloud G11.11-0.12, is a genuine small-scale system (less than 0.6 pc) in a nascent, pre-feedback phase, and that its gas kinematics reveal gravity-driven inflow toward the hub. Using ALMA N2H+ observations, the authors confirm the filamentary morphology first seen in JWST images, measure steep on-sky velocity gradients of about 5 and -7 km/s/pc on opposite sides of the hub, and find that the velocity-gradient vectors align with the gravitational force vectors along the filaments. This alignment is the key evidence that mass assembly is dominated by gravity, supporting clump-fed accretion scenarios for high-mass star formation. The paper also highlights a wiggled funnel morphology in position-position-velocity space, suggesting that sub-filaments and transverse gas flows matter for feeding the hub. If correct, it provides a rare observational snapshot of a high-mass protostar still being built by its parent gas streams.","feed_headline":"Gravity steers gas into a young massive star's hub","feed_subtitle":"JWST-ALMA data reveal a sub-0.6 pc hub-filament system with inflowing gas around G11P1—a rare pre-feedback view.","key_machinery":"The central machinery is the comparison of sky-projected velocity-gradient vectors with sky-projected gravitational force vectors, both derived from the ALMA N2H+ data. Velocity gradients are computed from the centroid-velocity field via second-order central differences of the line-of-sight velocity in right ascension and declination; gravitational force vectors are computed from the H2 column density map, which is obtained by converting N2H+ column densities using a fixed abundance ratio of 3 times 10 to the minus 10. The relative angle between the velocity-gradient and gravitational-force vectors along the filaments is the diagnostic that separates gravity-driven inflow from other mechanisms. A second diagnostic is the V-shaped velocity profile and the radial decrease of the velocity-gradient magnitude toward the hub, which the paper reads as evidence of convergent gas streams onto the protostar.","core_discovery":"On the paper's own terms, the central claim is that ALMA N2H+ observations confirm G11P1-HFS as a small-scale (less than about 0.6 pc) hub-filament system in a nascent phase, with five hub-joining filaments that spatially match the structures traced by JWST. The key kinematic discovery is a V-shaped velocity profile across the hub: velocity gradients of about -7 and +5 km/s/pc on the northeast and southwest sides, respectively, with the gradient magnitude decreasing toward the hub center, which the authors interpret as converging, mass-accreting flows onto the protostar G11P1. In support, the projected gravitational force vectors derived from the H2 column density map are well aligned with the on-sky velocity-gradient vectors along the filament spines, with a mean relative angle of about 41 degrees and median of about 39 degrees, indicating that gravity rather than turbulence drives the inflow. The paper further reports that the system appears as a wiggled funnel in position-position-velocity space, a signature the authors associate with sub-filamentary or transverse gas flows feeding the main hub.","pith_inferences":["If the abundance-uncertainty caveat is set aside, the velocity-gradient versus gravitational-force alignment method could be applied systematically to other nascent hub-filament candidates found in ALMA surveys, testing whether gravity-driven inflow is the norm at sub-parsec scales.","The wiggled-funnel signature might be a generic observational marker for hierarchical, sub-filamentary feeding, and could be searched for in existing ALMA position-position-velocity cubes of other protostellar hubs without requiring new observations.","Comparing G11P1-HFS with the more evolved Mon R2 system suggests an observational sequence in which stellar feedback progressively erases the funnel signature, which could be turned into an evolutionary diagnostic for hub-filament systems.","The missing-flux discrepancy in accretion rates could be checked by combining ALMA data with single-dish N2H+ observations; if the rates rise to model values, the nascent-phase interpretation would be reinforced."],"forward_implications":["G11P1-HFS is confirmed as a real small-scale hub-filament system, strengthening the evidence that hub-filament structures are hierarchical across scales from less than 0.1 pc to more than 10 pc.","The alignment of velocity-gradient and gravitational-force vectors along the filaments implies that gravity, not turbulence or magnetic pressure, drives mass inflow toward the high-mass protostar G11P1.","The decreasing velocity gradient toward the hub and the V-shaped velocity profile are consistent with a mass-accreting hub, supporting clump-fed accretion scenarios such as global hierarchical collapse and inertial inflow.","The wiggled funnel morphology in position-position-velocity space suggests that sub-filaments and transverse gas flows contribute to mass transport into the hub, refining how funnel-inflow signatures should be interpreted.","The estimated filament accretion rates of 0.92 to 1.44 solar masses per million years are about two orders of magnitude below values from models and larger-scale observations, which the paper attributes to missing flux from interferometric observations."],"supporting_citations":[{"why":"Identified G11P1-HFS as a candidate hub-filament system in JWST images, providing the morphological basis that this paper confirms with ALMA.","marker":"Dewangan et al. 2024a"},{"why":"Supplied the ALMA N2H+ (1-0) data cube used for all kinematic and column-density analyses.","marker":"Gieser et al. 2023"},{"why":"Provides the fixed N2H+/H2 abundance ratio of 3e-10 used to convert N2H+ column densities to H2 column densities.","marker":"Caselli et al. 2002"},{"why":"Gives the method for computing sky-projected gravitational force vectors from the column density map.","marker":"Wang et al. 2022"},{"why":"Establishes the V-shaped velocity profile and funnel morphology in PPV space as signatures of mass-accreting hubs.","marker":"Zhou et al. 2023"},{"why":"Provides the accretion-rate formula used to estimate mass inflow along the filaments.","marker":"Kirk et al. 2013"}],"fun_headline_variants":["Nascent hub-filament system shows gravity-fed flows","Gravity not turbulence drives gas into this star hub","ALMA+JWST catch gravity feeding a tiny hub","Funnel-shaped velocity map reveals gravity-fed inflow","Tiny nascent hub-filament system fed by gravity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The gravity-driven inflow conclusion rests on the H2 column density map and its derived gravitational force vectors, which assume a fixed N2H+ abundance ratio of 3e-10; the paper itself notes this introduces 60-70 percent uncertainties in N(H2) and mass, and those uncertainties are not propagated into the velocity-gradient versus gravitational-force alignment.","fun_headline_variants_meta":{"raw":{"variants":["Nascent hub-filament system shows gravity-fed flows","Gravity not turbulence drives gas into this star hub","ALMA+JWST catch gravity feeding a tiny hub","Funnel-shaped velocity map reveals gravity-fed inflow","Tiny nascent hub-filament system fed by gravity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000807,"raw_usage":{"total_tokens":3614,"prompt_tokens":1089,"completion_tokens":2525,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":705,"completion_tokens_details":{"reasoning_tokens":2447}},"tokens_in":705,"tokens_out":2525,"duration_ms":20704,"temperature":1.0,"reasoning_tokens":2447,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:49:08.457103+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute H2 column density independently from the ALMA dust continuum at matched resolution and recompute the gravitational force vectors; if the low relative angle between velocity-gradient and gravitational-force vectors along the filaments does not survive a factor-of-two variation in the assumed abundance or the continuum-based mass map, the gravity-driven inflow interpretation would be falsified.","supporting_citations":[{"cited_title":"J., Myers , P","cited_arxiv_id":null,"evidence_quote":"Provides the fixed N2H+/H2 abundance ratio of 3e-10 used to convert N2H+ column densities to H2 column densities."},{"cited_title":"High-resolution APEX/LAsMA $^{12}$CO and $^{13}$CO (3-2) observation of the G333 giant molecular cloud complex : I. Evidence for gravitational acceleration in hub-filament systems","cited_arxiv_id":"2305.12573","evidence_quote":"Establishes the V-shaped velocity profile and funnel morphology in PPV space as signatures of mass-accreting hubs."}],"review_version":1}