{"id":"43d2e7bc-4155-4ecd-b5b4-13670e1b24e7","arxiv_id":"2504.12381","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The IRAS Vela Shell is a bowl-shaped dust cavity at 353 pc with a mass of 5.1e4 solar masses, likely inflated by an HII region and one to two supernovae, with stellar winds subdominant.","lead":"Astronomers used a new 3D map of interstellar dust to measure the shape and mass of the IRAS Vela Shell, a giant bubble in the southern sky. They found it is a bowl-shaped cavity about 350 parsecs from the Sun, and its expansion was likely driven by hot massive stars and possibly one or two ancient supernova explosions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The single adopted expansion velocity (12 ± 3 km/s) is measured only from cometary globules in the lower bowl, yet is applied to the entire IVS; all dynamical quantities scale with it, and the paper itself flags the mismatch for the upper filament (§4.2).","rationale":"The reader's weakest_assumption identifies exactly the load-bearing point, and the manuscript itself flags it in §4.2, so my read agrees. I do not see a more serious internal inconsistency: the shell-center choice is tested against an independent peak-finding method in Appendix B, the mass estimate is robust to the smoothing-kernel variation, and the CG30 stellar distances place at least one globule on the shell surface. The expansion velocity is therefore the hinge on which the dynamical part of the central claim turns, and it is the least secure condition because the tracer population (cometary globules) samples only the dense lower bowl. A re-analysis of the CO data can settle whether the velocity is representative. The other weaknesses—private-communication values and the unquantified H-alpha correlation—are real but secondary; they would tighten the paper rather than overturn the geometry. Because the reader already conditioned the verdict on this assumption, my read does not change the verdict.","tokens_in":27489,"tokens_out":5578,"duration_ms":64189,"concrete_test":"Re-derive the shell kinematics from the original 12CO data used by Sridharan (1992a) and Rajagopal & Srinivasan (1998), fitting an expanding-shell model separately for the lower bowl and the diffuse upper filament, e.g., in Galactic latitude bins b < -5 deg and b > -5 deg. If the upper-filament expansion velocity differs from 12 ± 3 km/s by more than the quoted 1-sigma uncertainty, recompute the mass-weighted momentum, kinetic energy, N_SN, and t_dyn from Eqs. 6-7 and 17-19; if the mass-weighted momentum drops by more than about 30%, the dynamical conclusions require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's dynamical conclusions (momentum 6.0e5 M_sun km/s, energy 7.1e49 erg, N_SN = 1.7, t_dyn = 1.6 Myr) all pass through p = M_IVS * v_exp with v_exp = 12 ± 3 km/s, a value derived from 12CO radial velocities of cometary globules that trace only the dense lower bowl (Sridharan 1992a; Rajagopal & Srinivasan 1998). Section 4.2 explicitly concedes that 'this velocity range may not perfectly match the top filament due to density differences' and that the same v_exp is nevertheless applied to the entire IVS structure. This is the least secure condition for the central dynamical claim: if the diffuse upper filament is not participating in the same expansion, or expands at a different rate, then the total momentum and the derived supernova count, dynamical age, and HII-vs-SN balance are systematically shifted. The CG30 validation in §5.1 shows only that one globule's gas velocity is consistent with the adopted expansion; it does not independently measure the upper filament. Secondary but relevant: the LyC absorption fraction (21%) and the expected supernova counts (1.16 and 1.21) are private-communication values, so the HII-region and supernova branches are not independently checkable from the published record. The geometric claim—a coherent bowl-shaped cavity at median distance 353 pc—is much better supported and is not threatened by this critique.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the Edenhofer et al. (2024) 3D dust map to reconstruct the three-dimensional geometry of the IRAS Vela Shell (IVS). The authors identify a bowl-shaped cavity below the Galactic plane plus a more diffuse upper filament, at a median distance of about 353 pc, and derive a shell mass of 5.1e4 Msun, a momentum of 6.0e5 Msun km/s, and a kinetic energy of 7.1e49 erg by adopting an expansion velocity of 12 +/- 3 km/s from earlier CO studies of cometary globules. They then compare this momentum with the expected contributions from stellar winds, HII-region pressure, radiation pressure, and supernovae, concluding that winds are subdominant and that one to two supernovae, or an HII region plus radiation, could plausibly drive the shell. A traceback analysis of zeta Puppis, RX J0720.4-3125, and nearby young clusters is used to identify candidate supernova progenitors. The paper also argues for a physical association between the IVS, the Gum Nebula, and the cometary globule system, anchored by the distance and velocity of CG30.","tokens_in":27798,"tokens_out":6103,"duration_ms":62668,"significance":"If the results hold, this is the first 3D geometric model of the IVS and it settles a long-standing distance controversy by placing the shell at roughly 353 pc rather than the older 450 pc estimate. The geometric and mass derivation is carefully done: the authors propagate statistical uncertainty through twelve dust-map samples, vary the smoothing kernel to estimate systematics, compare their boundary-finding method with an independent peak-finding approach, and release their fiducial models and machine-readable astrometry. The CG30 association provides a valuable independent distance/velocity anchor. The dynamical and feedback conclusions are less secure because they scale linearly with a single adopted expansion velocity and because two key quantitative inputs come from private communications, but the paper is honest about these limitations. The work is a useful template for connecting 3D dust mapping with stellar-feedback studies in the local ISM.","major_comments":[{"comment":"The adopted v_exp = 12 +/- 3 km/s is derived from 12CO radial velocities of cometary globules that trace mainly the dense lower bowl, yet the same value is applied to the entire IVS surface, including the diffuse upper filament. The text explicitly concedes that this velocity may not match the top filament, and the CG30 validation in Section 5.1 tests only one globule. Because the momentum p, kinetic energy Ek, supernova count N_SN (Eq. 17), and dynamical age tdyn (Eq. 18) all scale directly with this single velocity, the dynamical conclusions are not robust unless the velocity field of the diffuse component is independently constrained, or unless the results are presented explicitly as conditional on this assumption with a sensitivity analysis that shows how the conclusions change for plausible alternative velocities.","section":"Section 4.2, Eqs. (6)-(7)"},{"comment":"Quantitative values that are central to the feedback budget come from private communications: the 21% Lyman-continuum absorption fraction from McCallum et al. (submitted 2025, private communication) and the expected supernova counts of 1.16 and 1.21 from Swiggum et al. (2024, private communication). These numbers cannot be checked from the published record, but they are used to conclude that HII-region feedback and supernova feedback are consistent with the measured shell momentum. Please include the methodology and resulting values in an appendix, or replace them with quantities from published, citable sources.","section":"Sections 5.2.2 and 5.2.7"},{"comment":"The shell center is identified by visual inspection of the 3D dust map, and the robustness statement that results remain consistent for centers within 15 pc is not documented. Since the radial profiles, shell boundaries, and mass integrals are all measured with respect to this center, please provide a quantitative demonstration of the claimed robustness, for example by showing the variation in total mass, peak radius, and shell thickness over the explored range of centers.","section":"Section 3.1"}],"minor_comments":[{"comment":"The peak-density entry '0.31+11.35-0.28' is an unusual asymmetric uncertainty notation; since the caption states that the values represent 95% percentile ranges, please report the actual percentile values consistently.","section":"Table 1"},{"comment":"The equality between 4*pi*r_i^2*dr/N_pix and (4*pi/3)*(r_{i+1}^3 - r_i^3)/N_pix is exact only in the limit dr -> 0; with the adopted 1 pc sampling it is an approximation, and the text should either state this or use the cell-volume expression directly.","section":"Equation (5)"},{"comment":"The fewer than 10 rays with unrecoverable outer boundaries are said to be rare, but the treatment of these rays in the mass integration is not described; please state whether such rays are excluded or assigned a default outer boundary.","section":"Section 3.1"},{"comment":"The claimed strong spatial correlation between the projected IVS and the H-alpha emission is descriptive rather than quantitative; a simple quantitative metric, such as a contour-overlap fraction or correlation coefficient, would make the association more testable.","section":"Section 5.1"},{"comment":"The statement that less than 0.3% of the gas on the IVS surface is optically thick (column density NH >= 1e21 cm^-2) is not derived in the text; please include the column-density calculation or a specific reference for this threshold.","section":"Section 5.2.2"},{"comment":"There are several typographical and phrasing issues: 'fidicual' for 'fiducial', the title contains 'V ela', and the sentence preceding Table 1 ('the uncertainties are not akin to ...') is difficult to parse. In addition, McCallum et al. is cited both as 'submitted 2025, private communication' and as a published MNRAS entry without volume/page details; please make the citation status consistent.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe headline: this is the first 3D geometric model of the IRAS Vela Shell, from the E24 dust map, and it makes a credible case that the IVS is a coherent bowl-shaped cavity at ~353 pc that belongs to the Gum Nebula complex. If you care about local stellar feedback or 3D dust mapping, this paper is worth engagement.\n\nWhat is genuinely good: the shell geometry work is careful. The boundary finding uses inflection points in the density profiles, they check against a second method (O'Neill et al.'s peak finding), and they propagate both the map samples and the smoothing kernel into the uncertainties. They also release the fiducial and alternate models. The mass derivation is solid. A nice external anchor: six young stars in the cometary globule CG30 sit at 359 pc and land on the shell edge, which independently supports the geometry.\n\nThe soft spots are all in the dynamical interpretation. Momentum, energy, supernova count, and dynamical age pass through one number: the 12±3 km/s expansion velocity measured from CO in the dense lower-bowl globules. The paper says outright in §4.2 that this may not match the upper filament, and applies it anyway. That is an honest caveat, but it means the SN/HII balance and the 1–2 supernovae conclusion are conditional on an untested uniformity assumption. It doesn't threaten the geometric claims.\n\nTwo smaller items: the 21% LyC absorption fraction and the per-cluster expected SN counts (1.16, 1.21) come from private communications, so those numbers are not independently checkable in the published record. And the Hα–dust correlation is stated visually rather than measured. Both are easy to fix: publish the numbers, add a quantitative cross-correlation.\n\nThis should go to peer review, not be desk rejected. The geometry, distance, and mass are the real results and they hold up. The dynamical story is a plausible scaffold, appropriately hedged in the text, but should not be taken as a strong quantitative budget until the expansion velocity is mapped across the whole shell and the private numbers are on the record.","headline":"First 3D geometry of the Vela shell holds up; the dynamical budget leans on one expansion velocity and a few private communication values.","tokens_in":28459,"tokens_out":4151,"would_cite":true,"duration_ms":41584,"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":"The IRAS Vela Shell is a single bowl-shaped cavity at 353 parsecs, and its expansion appears driven by an HII region or one to two supernovae rather than stellar winds.","keywords":["IRAS Vela Shell","Gum Nebula","3D dust mapping","stellar feedback","HII region","supernova remnants","cometary globules","interstellar medium"],"falsifier":"Measure the radial velocity of atomic or molecular gas in the diffuse upper filament of the IVS; if that gas is not expanding at roughly $12\\pm3$ km/s away from the same center, the single-expanding-shell model and its momentum, supernova-count, and age conclusions would need to be revised.","tokens_in":27283,"feed_emoji":"🌌","tokens_out":8261,"duration_ms":74169,"temperature":0.7,"pith_summary":"This paper claims that the IRAS Vela Shell — a conspicuous ring of infrared emission toward the Gum Nebula — is a single coherent three-dimensional cavity, located at a median distance of about 353 parsecs from the Sun rather than the 450 parsecs assumed in older work. Reconstructing the shell from a parsec-resolution 3D dust map, the authors find a dense bowl-shaped lower section below the Galactic plane and a diffuse filament above it, and they place the shell's total mass near $5.1\\times10^{4}$ solar masses. Combining that mass with published expansion velocities of the cometary globules that trace the shell's edge, they obtain a momentum of about $6\\times10^{5}$ solar masses km/s and a kinetic energy of about $7\\times10^{49}$ erg. They argue that the dust-traced cavity and the Gum Nebula's H-$\\alpha$ emission are co-spatial, that stellar winds supply only a few percent of the momentum, and that the remaining budget can be met either by an expanding HII region ionized by $\\gamma^2$ Velorum and $\\zeta$ Puppis or by one to two recent supernovae.","feed_headline":"IRAS Vela Shell is a 353-pc cavity, not a distant ring","feed_subtitle":"3D dust map ties the shell to the Gum Nebula; HII pressure or 1-2 supernovae drive it.","key_machinery":"The load-bearing mechanism is 3D dust mapping itself: the Edenhofer et al. (2024b) map converts the reddening of millions of Gaia-distance stars into a three-dimensional extinction density at parsec-scale resolution. The paper samples that map along 196,608 radial rays from a chosen shell center, converts extinction density to hydrogen volume density with $n_{\\rm H}=1653\\,A'_{\\rm ZGR}\\ \\mathrm{cm}^{-3}$, smooths each profile with a Gaussian kernel, and locates the shell boundaries where the second derivative of the density profile changes sign. Those boundaries define a 3D shell surface that can be integrated along each ray to obtain the mass, and the resulting geometry is combined with literature expansion velocities to build a momentum budget that is compared against stellar winds, HII region thermal and radiation pressure, and supernova momentum injection.","core_discovery":"The central claim is that the IRAS Vela Shell is not a foreground or background projection but a real, approximately spherical but incomplete cavity: a dense bowl about 70 pc below the Galactic plane and a diffuse filament about 70 pc above it, centered near (x,y,z) = (-67,-334,-65) pc in heliocentric Galactic Cartesian coordinates. The paper establishes the shell's geometry by sampling the Edenhofer et al. (2024b) 3D dust map along 196,608 radial rays, converting extinction density to hydrogen volume density, and locating the inner and outer shell boundaries at inflection points of the density profile. It derives a shell mass of $5.1^{+2.4}_{-2.4}\\times10^{4}$ solar masses, a shell momentum of $6.0^{+4.7}_{-3.4}\\times10^{5}$ solar masses km/s, and a dynamical age of $1.6^{+1.5}_{-0.6}$ Myr, and it identifies $\\gamma^2$ Velorum and $\\zeta$ Puppis as currently lying inside the shell. The paper's conclusion is that the shell's expansion is most plausibly powered by a combination of HII region pressure and one to two supernovae, with stellar winds subdominant.","pith_inferences":["The same boundary-from-inflection method could be applied to other H-alpha superbubbles; if bowl-plus-filament shapes are common, many 'shells' in projection are better described as partial cavities, which would change how their masses and ages are derived.","Because the near side of the IVS appears to share a wall with the Local Bubble, a testable prediction is that the two cavities have interacted: the IVS's expansion speed and dust distribution should show a flattened or compressed interface toward the Sun, and the Local Bubble's own history may have set the initial conditions for the IVS's evolution.","The momentum budget treats the HII region as a single pressure source, but the paper's own estimate that only about 21% of LyC photons are absorbed inside the IVS implies that the ionizing stars also photo-evaporate the cometary globules; measuring the globules' mass-loss rates would provide an independent check on whether the HII region can simultaneously drive the shell and erode its fragments.","If the upper filament's expansion velocity could be measured directly and differs from 12 km/s, the single-shell interpretation would split into two structures with possibly separate origins; this is a direct observational test of the paper's main geometric claim."],"forward_implications":["The IVS and the Gum Nebula are physically associated, so the shell's distance of 353 pc revises the distance and size estimates inherited from older work.","The massive stars $\\gamma^2$ Velorum and $\\zeta$ Puppis currently lie inside the shell and can supply enough ionizing photons to sustain an HII region of the IVS's size, making them live candidates for the driver.","If supernovae are the primary driver, only one to two events are required, and the candidate clusters identified in the paper would be expected to have produced about 1.2 supernovae in the past 3 Myr, matching the requirement.","The dynamical age of roughly 1.6 Myr is much younger than the 10-20 Myr proposed earlier, implying that the IVS records a recent feedback event rather than the older cluster formation episode.","The cometary globule CG30, with its embedded young stars at 359 pc, lies on the shell surface, directly connecting the shell's dust geometry to the expanding globule system."],"supporting_citations":[{"why":"Supplies the 3D dust map and its 12 posterior samples, the primary data from which the shell geometry and mass are derived.","marker":"Edenhofer et al. (2024b)"},{"why":"Provides the CO radial velocities of cometary globules from which the adopted expansion velocity of 12 ± 3 km/s is taken.","marker":"Sridharan (1992a)"},{"why":"Independent measurement of the IVS expansion velocity from CO observations of the whole region, used to set the velocity prior.","marker":"Rajagopal & Srinivasan (1998)"},{"why":"Established the IRAS Vela Shell as a named structure and provided the first projection-based size and distance estimates this paper revises.","marker":"Sahu & Sahu (1992)"},{"why":"Identifies the six young stars embedded in CG30 whose mean distance places the globule on the shell surface.","marker":"Yep & White (2020b)"},{"why":"Provides the young cluster catalog and expected supernova counts used to identify candidate sources within the past 3 Myr.","marker":"Swiggum et al. (2024)"},{"why":"Supplies the momentum-per-supernova sampling and shell-expansion formalism used for the supernova count and dynamical age.","marker":"Bialy et al. (2021)"},{"why":"The earlier 10-20 Myr supernova scenario that this paper's 1.6 Myr dynamical age directly challenges.","marker":"Cantat-Gaudin et al. (2019)"}],"fun_headline_variants":["IRAS Vela Shell is a bowl-shaped cavity, not a ring","3D map ties IRAS Vela Shell to Gum Nebula feedback","Shell momentum suggests HII region or 1-2 supernovae","Gamma2 Vel and zeta Pup sit inside IRAS Vela Shell","Dust map shows IRAS Vela Shell is a real 3D cavity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The shell expansion velocity of $12\\pm3$ km/s, measured from CO radial velocities of cometary globules in the lower bowl, is applied to the entire shell surface including the diffuse upper filament; all momentum, energy, supernova-count, and age estimates scale linearly with this velocity.","fun_headline_variants_meta":{"raw":{"variants":["IRAS Vela Shell is a bowl-shaped cavity, not a ring","3D map ties IRAS Vela Shell to Gum Nebula feedback","Shell momentum suggests HII region or 1-2 supernovae","Gamma2 Vel and zeta Pup sit inside IRAS Vela Shell","Dust map shows IRAS Vela Shell is a real 3D cavity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00079,"raw_usage":{"total_tokens":3596,"prompt_tokens":1171,"completion_tokens":2425,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":787,"completion_tokens_details":{"reasoning_tokens":2327}},"tokens_in":787,"tokens_out":2425,"duration_ms":15749,"temperature":1.0,"reasoning_tokens":2327,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:32:53.602993+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the radial velocity of atomic or molecular gas in the diffuse upper filament of the IVS; if that gas is not expanding at roughly $12\\pm3$ km/s away from the same center, the single-expanding-shell model and its momentum, supernova-count, and age conclusions would need to be revised.","supporting_citations":[{"cited_title":"1998, Journal of Astrophysics and Astronomy, 19, 79, doi: 10.1007/BF02714912","cited_arxiv_id":null,"evidence_quote":"Independent measurement of the IVS expansion velocity from CO observations of the whole region, used to set the velocity prior."}],"review_version":1}