{"id":"b592a80c-35c4-4602-aa6c-2e44006af5a1","arxiv_id":"2411.16372","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A spiral pattern in the stellar velocity dispersion map of NGC 4030 correlates with star-forming regions, supporting the idea that stars inherit the motion of their parent gas and are later heated.","lead":"Using MUSE observations of the spiral galaxy NGC 4030, astronomers detected a spiral pattern in the velocity dispersion of its stars, aligned with regions of active star formation. The result suggests newborn stars inherit the calm motion of the gas they form from, then gradually heat up over time, offering a new way to study how galaxy disks evolve.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The grand-design sigma spiral may be an artifact of KINEMETRY harmonic subtraction; a null-model or independent-smooth-subtraction test is required before the detection is secure.","rationale":"The reader's conditional verdict already identifies the KINEMETRY subtraction as a key assumption, and the present stress-test agrees that this is the most load-bearing step: the paper's headline discovery is a spiral pattern in the sigma residual map, and all of the maps used to define and compare the spiral arms pass through the same harmonic-subtraction procedure. The internal consistency checks with different templates and fitters show that the measured sigma maps are robust, but they do not test whether the subtraction of a smooth elliptical model can generate a coherent spiral-shaped residual of the observed amplitude. The EW(H-alpha)-based radial profiles are a valuable independent confirmation of the local anti-correlation between sigma and star formation, yet they do not establish that the residual spiral itself is intrinsic. A null-model or independent-smooth-fit test would directly settle this concern, and it is a reasonable requirement before the strong claim of a 'striking grand-design spiral pattern' is accepted. Since the reader's conditional verdict already requires such validation, no change to the verdict is recommended.","tokens_in":5,"tokens_out":5652,"duration_ms":130246,"concrete_test":"Re-run the analysis replacing the KINEMETRY subtraction with an independent smooth model, e.g., a 2D spline or Gaussian process fitted to the observed sigma map after masking the manually selected spiral-arm spaxels. If the residual spiral disappears or drops below the claimed 2-8 km/s amplitude, the pattern is an artifact of the KINEMETRY harmonic subtraction rather than an intrinsic stellar kinematic feature.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central detection is a residual spiral pattern in the stellar velocity dispersion map after subtracting a KINEMETRY harmonic model. The paper does not specify the harmonic order, the radial range, or the allowed variation of the ellipse parameters used in this subtraction, and it does not validate the subtraction on a null test. Because the spiral arms are themselves identified in a residual map produced by the same algorithm (HST F450W after KINEMETRY subtraction), the sigma residual pattern may reflect systematic misfit between the true, non-elliptical sigma distribution and the smooth elliptical harmonic model. The 2-8 km/s arm/inter-arm residual is only slightly larger than the reported 1.7 km/s median uncertainty, so a modeling residual of a few km/s could produce the same signature. The paper's consistency checks (multiple templates, fitters, nebular subtraction) validate the input sigma maps, but not the subtraction step. The EW(H-alpha)-selected radial profiles of observed sigma provide an important independent check of the sigma-age correlation, but they do not by themselves establish that the residual spiral pattern is intrinsic. If the KINEMETRY subtraction artificially creates a coherent spiral residual, the paper's 'striking grand-design spiral pattern in sigma' claim collapses, even though some of the underlying local correlations may survive.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyses MUSE IFU data of the grand-design spiral galaxy NGC 4030, deriving stellar velocity dispersion and luminosity-weighted age maps with multiple fitting methods (PPXF and BAYES-LOSVD; IUS and MILES templates). After subtracting a smooth KINEMETRY harmonic model from the observed maps, the authors report a grand-design spiral pattern in the stellar velocity-dispersion residual map, with 2-8 km/s lower sigma inside the spiral arms, where H-alpha and EW(H-alpha) are enhanced and stellar ages are 0.25-1.25 Gyr younger than in inter-arm regions. They interpret this as evidence that stars inherit the low velocity dispersion of the progenitor cold molecular gas and are subsequently heated by gravitational interactions, comparing the resulting age-velocity relation with cosmological zoom-in simulations (Bird et al. 2021; Martig et al. 2014).","tokens_in":10484,"tokens_out":7653,"duration_ms":67658,"significance":"If the detection is robust, it is an important new observable: a large-scale correlation between stellar kinematics, stellar age, and ongoing star formation in an external disk galaxy, with direct bearing on the age-velocity relation and disk-heating mechanisms. The paper's strength is the extensive cross-validation of the kinematic maps: different template libraries, fitting codes, numbers of Gauss-Hermite moments, and with/without nebular-continuum subtraction all reproduce the pattern (Fig. A.1). The EW(H-alpha)-selected radial profiles (Fig. 3b,c) provide an independent check on the sigma-age correlation that does not rely on the manually defined spiral-arm masks. The comparison with simulations is suggestive, though the use of luminosity-weighted ages and a constant sigma offset in the g106 comparison weakens quantitative conclusions.","major_comments":[{"comment":"The KINEMETRY harmonic subtraction is the step that defines the residual maps used to measure the spiral pattern and its amplitude. The manuscript does not report the harmonic order, the radial range, or the allowed variation of the ellipse parameters, and it does not provide a null test of the subtraction. Because the residual amplitude (2-8 km/s) is only slightly larger than the median per-spaxel uncertainty (1.7 km/s), a smooth-model misfit of a few km/s could in principle generate a coherent spiral-like residual in a non-elliptical sigma map. The consistency checks in Fig. A.1 validate the input sigma maps, not the subtraction step. Please specify the KINEMETRY settings and demonstrate that the residual spiral is robust to the subtraction method, e.g., by comparing with an alternative smooth baseline (azimuthal median, spline, or unsharp mask) or by injecting a known spiral pattern into a smooth map.","section":"Section 2 / Fig. A.5"},{"comment":"The paper quotes a median sigma uncertainty of 1.7 km/s and residual amplitudes of 2-8 km/s, but no uncertainties are given for the radial profiles, and the statistical significance of the arm/inter-arm difference is not quantified. The authors should report the uncertainty on each radial-bin value, including spatial covariance, and a significance level for the spiral pattern, ideally from Monte Carlo realizations that propagate the full analysis chain.","section":"Section 3, Fig. 1c"},{"comment":"The core data processing relies on an unpublished pipeline described only as 'Breda et al., 2024b (in prep.)', which is not listed in the references. This makes the results difficult to reproduce. The paper should provide a detailed account of the pipeline's steps in the appendix, include a reference to a public repository or companion paper, or state explicitly which steps are new versus standard.","section":"Section 2 / Breda et al. 2024b"},{"comment":"The AVR comparison with the simulation g106 involves 'adjusting sigma by summing a constant value'. Since this is an unconstrained parameter, the apparent agreement with the S-shape should be tested for sensitivity to the offset; alternatively, the offset should be physically motivated. Without this, the comparison remains speculative and should be framed accordingly.","section":"Section 4 / Fig. 3a"}],"minor_comments":[{"comment":"There are several typographical and formatting errors: 'EWH(α)' for 'EW(Hα)', 'week nebular emission' for 'weak nebular emission', '0.25 Gr' for '0.25 Gyr' in the Fig. 2 caption, and 'van van Donkelaar' in the Introduction. These should be corrected.","section":"Throughout"},{"comment":"The statement about the full-map radial profile 'showing however lower values, due to the inclusion of the inner regions where sigma is systematically lower' is unclear; please clarify whether the full-map values are lower than the outside-arm values at all radii and why this is expected.","section":"Fig. 1c"},{"comment":"The reference list contains entries not cited in the text (e.g., Seabroke & Gilmore 2007) and formatting inconsistencies (e.g., 'PASJ, 44, 601Lea17'); please harmonize with the journal style.","section":"References"},{"comment":"The abstract introduces 'A VR' with a space and uses the vague term 'Post-processing'; please state explicitly that the pattern is detected after subtracting a smooth baseline model and define the AVR abbreviation at first use.","section":"Abstract"},{"comment":"The Fig. A.1 caption states that 'The sigma_spiral features are evident in all cases', but the main text emphasizes the residual map. Please state explicitly whether the spiral is visible in the raw sigma maps before KINEMETRY subtraction, and if so, show a direct comparison between the raw and residual maps in the main text.","section":"Fig. A.1"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the KINEMETRY subtraction. I would ask the authors to add a null test and full details of the subtraction before acceptance. If the spiral is indeed present in the raw sigma maps (as the Fig. A.1 caption suggests), the paper could be accepted after a moderate revision that presents this evidence more prominently and quantifies the significance. The unpublished pipeline reference is also a reproducibility concern that the editor may wish to see addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper is worth reading and worth refereeing, but the headline claim — a grand-design spiral pattern in the stellar velocity dispersion — is not yet nailed. The authors have done something genuinely new: using high-quality MUSE data of NGC 4030, they see coherent 2–8 km/s drops in sigma along the spiral arms, and these drops track younger ages and higher EW(Ha). They check this with multiple template sets (IUS, MILES), two fitters (pPXF, Bayes-LOSVD), and with/without nebular-continuum subtraction. That is a solid set of internal consistency checks, and the correlation between low sigma, young luminosity-weighted ages, and star formation activity is convincing in the radial-profile EW(Ha) selection, which does not depend on the manually drawn arms.\n\nThe soft spot is the residual map itself. The pattern is defined as the sigma map minus a KINEMETRY harmonic model. The paper never states the harmonic order, the radial range, or the allowed variation of ellipse parameters. More importantly, there is no null test: subtract the same harmonic model from an artificially smooth or scrambled sigma map and see if coherent spiral residuals appear. Because the spiral arms are drawn on an HST residual map produced by the same algorithm, the sigma residual pattern could in principle be a misfit between the true, non-elliptical sigma distribution and the smooth elliptical model. A few km/s of systematic misfit is exactly the size of the claimed signal. So the stress-test concern is legitimate and lands on the central claim: the 'striking grand-design spiral pattern' is not yet demonstrated.\n\nThat said, the underlying physics correlation does not collapse. The EW(Ha)-selected radial profiles are independent of KINEMETRY and of the arm selection, and they show the same anti-correlation between sigma and star formation. So even if the grand-design spiral is an artifact, the age–sigma–star formation connection in this galaxy probably survives.\n\nThe AVR interpretation is qualitative. The comparison to Bird et al. and Martig et al. is reasonable but loose — they shift a constant to align with g106, and the paper admits the comparison is speculative. The S-shape is interesting but depends on FADO luminosity-weighted ages, which have their own biases. Not a fatal issue, but the paper should be more careful about claiming support from simulations.\n\nBottom line: this is a single-galaxy result with a low-amplitude signal and an unvalidated subtraction step. It deserves peer review, not because the claim is proven, but because it is a plausible new observable and the internal consistency checks are above average. The referee should ask for a null test and full disclosure of the KINEMETRY settings. I would not cite it as a secure detection yet, but I would bring it to a reading group for discussion.","headline":"A plausible but unproven detection; the sigma spiral needs a null test before it can be called secure, but the age–sigma–star formation correlation is solid.","tokens_in":11132,"tokens_out":3495,"would_cite":false,"duration_ms":30945,"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":"A grand-design spiral pattern in the stellar velocity dispersion of NGC 4030 shows stars forming in spiral arms are born colder and heat up as they age.","keywords":["galaxies: spiral","galaxies: evolution","galaxies: kinematics and dynamics","age-velocity relation","integral field spectroscopy","disk heating","velocity dispersion","NGC 4030"],"falsifier":"A decisive check would be to run the same harmonic-expansion subtraction on a velocity dispersion map of a galaxy with no spiral structure, or on a synthetic smooth galaxy with comparable noise; if a similar spiral residual appears, the pattern is an artifact. Alternatively, mapping the cold molecular gas velocity dispersion in NGC 4030 at comparable resolution should reveal a matching spiral of low-dispersion gas if the inheritance picture is correct.","tokens_in":10026,"feed_emoji":"🌀","tokens_out":6881,"duration_ms":58441,"temperature":0.7,"pith_summary":"This letter reports a grand-design spiral pattern in the stellar velocity dispersion map of the nearby spiral galaxy NGC 4030, seen after subtracting a smooth model of the galaxy's main body. Stars inside the spiral arms are colder (lower $\\sigma_\\star$ by 2–8 km/s), younger, and coincident with HII regions, so the pattern traces active star formation. The paper interprets this as evidence that stars inherit the velocity dispersion of the cold molecular gas they form from, so the age–velocity relation records both the turbulence of the gas at birth and later gravitational heating. A sympathetic reader would care because the same residual-map technique could expose disk-heating physics in other galaxies, and the data provide a sharp target for cosmological simulations.","feed_headline":"Stars are born cold in NGC 4030's spiral arms","feed_subtitle":"Young stars in the arms move 2–8 km/s slower than older stars, linking gas birth to disk heating.","key_machinery":"The analysis rests on harmonic-expansion modeling and subtraction: a two-dimensional map of each observable (velocity dispersion, age, H$\\alpha$) is fitted along best-fitting ellipses with a smooth multi-harmonic model, and the residual map retains only non-axisymmetric structure. Full spectral fitting of the integral-field data cube, with nebular-continuum subtraction and Voronoi binning, supplies the stellar kinematics and luminosity-weighted ages; selection of progressively higher EW(H$\\alpha$) spaxels provides an independent check that the arm/inter-arm contrasts are not artifacts of the manually drawn arm regions. The harmonic-expansion residual is the load-bearing object because it is what reveals the spiral pattern in $\\sigma_\\star$ in the first place.","core_discovery":"The central claim is that NGC 4030 shows a large-scale, grand-design spiral pattern in the residual stellar velocity dispersion map, with $\\Delta\\sigma_\\star$ systematically 2–8 km/s lower inside the arms than in the surrounding disk, matching regions of younger luminosity-weighted ages (0.25–1.25 Gyr younger) and enhanced H$\\alpha$ emission. This spatial coincidence is read as star formation from cold gas: newly born stars keep the low velocity dispersion of their progenitor gas, while older populations have been heated by cumulative gravitational interactions. The galaxy's age–velocity relation is S-shaped rather than a smooth power law, with a wide spread of $\\sigma_\\star$ at 7–8 Gyr ages that the authors attribute to a rapid formation phase across the sampled region. The paper argues that the observed AVR is consistent with high-resolution cosmological zoom-in simulations where birth dispersion and post-birth heating each contribute roughly half of the present-day dispersion.","pith_inferences":["Applying the same harmonic-expansion subtraction to other integral-field galaxy samples could reveal whether such sigma spirals are common or rare; the paper only demonstrates one clear case.","Luminosity-weighted ages dilute old stellar components, so the true age–sigma separation between arms and inter-arm regions may be larger than the 0.25–1.25 Gyr reported; mass-weighted stellar ages would test this.","If the inheritance picture is right, high-resolution CO observations of NGC 4030 should show a matching spiral pattern in the molecular gas velocity dispersion, a directly testable prediction.","The S-shape at 7–8 Gyr may encode a genuine formation epoch; comparing the location of that knee with simulated merger timescales could date the disk's assembly."],"forward_implications":["If the interpretation is right, spiral arms in this galaxy are sites where stars are currently being born from cold gas, and their present-day velocity dispersion measures the gas dispersion at birth rather than just accumulated heating.","The S-shaped age–velocity relation implies that disk heating models must include both a birth-dispersion term set by ISM turbulence and a post-birth heating term, rather than a single power law.","Residual velocity-dispersion maps can serve as a tracer of recent star formation and spiral structure even when photometric tracers are faint.","The measured 2–8 km/s contrast sets a quantitative benchmark that cosmological zoom-in simulations of Milky Way-mass galaxies should reproduce.","EW(H$\\alpha$)-selected profiles confirm the arm/inter-arm trends, so the result is not tied to the particular arm boundaries chosen."],"supporting_citations":[{"why":"Supplies the integral-field data cube and the galaxy parameters (distance, stellar mass, star-formation rate, effective radius) defining the analyzed system.","marker":"Erroz-Ferrer et al. 2019"},{"why":"The spectral synthesis approach used to derive stellar ages from the observed spectra.","marker":"Gomes & Papaderos 2017"},{"why":"The Voronoi binning method that sets the signal-to-noise of the spectra before kinematic and age extraction.","marker":"Cappellari & Copin 2003"},{"why":"The harmonic-expansion method used to model and subtract the underlying map structure, generating the residuals that reveal the spiral pattern.","marker":"Krajnović et al. 2006"},{"why":"The kinematic extraction code that measures the stellar velocity dispersion from the binned spectra.","marker":"Cappellari 2017"},{"why":"Provides the high-resolution cosmological zoom-in simulation whose age–velocity relation is compared to NGC 4030's, offering the birth-plus-heating interpretive framework.","marker":"Bird et al. 2021"},{"why":"The prior integral-field study of spiral-arm and disk kinematics that this paper extends by finding clear arm/inter-arm velocity-dispersion differences.","marker":"Pessa et al. 2023"},{"why":"Simulated age–velocity relations, especially for the galaxy resembling NGC 4030, whose S-shape supports the combined birth-heating interpretation.","marker":"Martig, Minchev & Flynn 2014a,b"}],"fun_headline_variants":["Spiral pattern in stellar velocity exposes disk heating","Young stars move slower in NGC 4030's spiral arms","S-shaped age-velocity relation in NGC 4030 hints at disk heating","Velocity spirals reveal stellar birth temperatures in NGC 4030","Cold birth, hot later: NGC 4030's spiral velocity signature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the smooth harmonic-expansion model subtracted from the velocity dispersion map does not itself create the spiral pattern; if the residual 2–8 km/s signal is an artifact of the subtraction, the claimed correlation between low dispersion, young age, and star formation collapses.","fun_headline_variants_meta":{"raw":{"variants":["Spiral pattern in stellar velocity exposes disk heating","Young stars move slower in NGC 4030's spiral arms","S-shaped age-velocity relation in NGC 4030 hints at disk heating","Velocity spirals reveal stellar birth temperatures in NGC 4030","Cold birth, hot later: NGC 4030's spiral velocity signature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001444,"raw_usage":{"total_tokens":5834,"prompt_tokens":977,"completion_tokens":4857,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":4769}},"tokens_in":593,"tokens_out":4857,"duration_ms":31854,"temperature":1.0,"reasoning_tokens":4769,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:11:35.110129+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check would be to run the same harmonic-expansion subtraction on a velocity dispersion map of a galaxy with no spiral structure, or on a synthetic smooth galaxy with comparable noise; if a similar spiral residual appears, the pattern is an artifact. Alternatively, mapping the cold molecular gas velocity dispersion in NGC 4030 at comparable resolution should reveal a matching spiral of low-dispersion gas if the inheritance picture is correct.","supporting_citations":[{"cited_title":"M., den Brok M., Onodera M., Brinchmann J., Marino R","cited_arxiv_id":null,"evidence_quote":"Supplies the integral-field data cube and the galaxy parameters (distance, stellar mass, star-formation rate, effective radius) defining the analyzed system."},{"cited_title":"& Copin Y., 2003, MNRAS, 342, 345","cited_arxiv_id":null,"evidence_quote":"The Voronoi binning method that sets the signal-to-noise of the spectra before kinematic and age extraction."},{"cited_title":"C., Loebman S","cited_arxiv_id":null,"evidence_quote":"Provides the high-resolution cosmological zoom-in simulation whose age–velocity relation is compared to NGC 4030's, offering the birth-plus-heating interpretive framework."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The prior integral-field study of spiral-arm and disk kinematics that this paper extends by finding clear arm/inter-arm velocity-dispersion differences."}],"review_version":1}