{"id":"76e3d3d2-fd70-436c-97cb-54f7a1a0f61c","arxiv_id":"2506.14218","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"PAH emission appears in 64% of 124 Herbig Ae/Be stars in re-reduced Spitzer spectra, peaks around 7000-11000 K, and the 6.2 micron band shifts blueward with stellar temperature.","lead":"This paper uses reprocessed Spitzer telescope spectra of 124 young, medium-weight stars to sort their surrounding dust and gas by chemical fingerprint. It finds which stars are most likely to show emission from large organic molecules, and shows that one emission band shifts with stellar temperature, giving future James Webb observations a target shortlist.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"6.2 µm–Teff shift lacks quantitative uncertainties and robustness checks; low-resolution fitting and nebular contamination could produce the trend.","rationale":"The paper is a competent survey with public data and reproducible methods, so I read the central claim in its strongest form: a temperature-driven blueward shift of the 6.2 µm PAH band, based on the largest HAeBe sample tested to date. For that claim to hold, two conditions must be met: the measured centroids must be accurate relative to the shift, and the shift must not be a proxy for another variable. Both conditions are currently unchecked. Sec. 3.6 does not report errors on fitted peak positions or a correlation statistic; Figure 5 is presented without a significance test. Given the low spectral resolution, centroid uncertainties may be comparable to the 0.03–0.09 µm Class A/B separation. The continuum subtraction is described only as 'cubic spline', with no anchor points or stability analysis; this is not an accusation of sloppiness, but the robustness of the fits is not demonstrated. In addition, the authors themselves note that Herbig Be stars show mixed Class A/B profiles and that reflection nebulae contribute Class A features. Since Herbig Be stars populate the high-TeF end, the correlation might track the presence of nebulosity rather than stellar temperature. The 6.0 µm non-shift is subject to the same deblending systematics; a fixed or weakly constrained 6.0 Gaussian could trivially produce a constant centroid. A robustness re-fit of the public spectra would settle these questions. For these reasons the conditional verdict is appropriate; my stress-test does not change it.","tokens_in":19941,"tokens_out":6560,"duration_ms":78024,"concrete_test":"Refit all 50 spectra with an automated, reproducible continuum method—e.g., PAHFIT or fixed spline anchor points at 5.0, 5.5, 6.5, and 6.9 µm, with two alternative anchor sets—and compute bootstrap uncertainties on the 6.2 and 6.0 µm centroids with Gaussian centroids left free. Report the Spearman correlation between 6.2 centroid and Teff with a 95% confidence interval for (i) all 50 sources and (ii) the subsample excluding HBe stars and known reflection nebulae. If the confidence interval includes zero, or if the correlation disappears after excluding nebulous sources, the claimed temperature-driven shift is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. 3.6 reports the central claim—a blueward shift of the 6.2 µm PAH peak with Teff—using only a scatter plot. No peak-position uncertainties, fitted slope, or correlation coefficient are given, and Figure 5 is not accompanied by a significance test. This matters because Spitzer IRS low-resolution spectra at 6.2 µm have R≈60–127 (Δλ≈0.05–0.10 µm), on the same order as the quoted Class A–B separation (6.19–6.28 µm). The centroids come from a cubic-spline continuum and double-Gaussian deblending whose anchor points and initial conditions are not specified in Sec. 3.6. If hotter stars have slightly different continuum shapes, or if the deblending drifts when the 6.0 µm component is weak, the apparent Teff trend could be a fitting artifact. The paper itself states that Herbig Be stars show a mix of Class A and B 'irrespective of Teff' and that associated reflection nebulae contribute Class A PAHs; this is a Teff-correlated environmental confound that is never removed. The same deblending affects the claimed 6.0 µm non-shift. Without a robustness test and significance estimate, the largest-sample claim is not yet supported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a mid-infrared spectroscopic survey of 124 Herbig Ae/Be stars using uniformly processed Spitzer IRS spectra from the CASSISjuice database. The authors visually classify the spectra into five categories (PAH-only, silicate-only, PAH+silicate, no feature, and HACs), compute detection frequencies as functions of effective temperature and the 2–24 micron spectral index, and report that the 6.2 micron PAH peak shifts blueward with increasing Teff while the 6.0 micron feature remains roughly constant. They also measure the 11.0/11.2 micron PAH ratio in a subset of stars and interpret the results as evidence for Teff-driven PAH processing and distinct molecular origins for the 6.0 micron band.","tokens_in":20172,"tokens_out":4372,"duration_ms":46146,"significance":"If the 6.2 micron–Teff correlation is quantitatively robust, this would be the largest sample yet used to test that peak shift and would provide valuable statistical evidence for temperature-driven PAH processing in intermediate-mass pre-main-sequence stars. The paper also offers a useful uniformly reduced sample and classification of the public SSHC catalog. However, the central quantitative claims are currently supported only by scatter plots and visual classifications, without per-point uncertainties, significance tests, or robustness checks against the fitting procedure. The paper's value as a catalog and classification resource is real, but the headline astrophysical conclusions need substantial additional analysis before they can be accepted.","major_comments":[{"comment":"The central claim that the 6.2 micron PAH peak shifts blueward with Teff is supported only by a scatter plot; no per-source peak-position uncertainties, fitted slope, correlation coefficient, or significance test are reported. At Spitzer IRS low resolution (R approximately 60–127, corresponding to Delta-lambda approximately 0.05–0.10 micron), the measurement uncertainty is comparable to the quoted Class A–B separation of 6.19–6.28 microns, so the claim requires a quantitative statistic (e.g., Spearman rank correlation with a bootstrap confidence interval) before it can be regarded as established.","section":"Section 3.6, Figure 5"},{"comment":"The continuum subtraction and deblending procedure is underspecified: the cubic-spline anchor points are not listed, the initial conditions and number of Gaussians in the double-Gaussian fit are not given, and no robustness test is presented. Because hotter stars may have different underlying continuum shapes, the apparent Teff trend could be a fitting artifact; the authors should test stability against alternative continuum definitions, anchor choices, and deblending assumptions, and report the resulting dispersion in fitted peak wavelengths.","section":"Section 3.6 and Figure A.1"},{"comment":"The aggregate Teff trend is vulnerable to a known environmental confound acknowledged in the same section: Herbig Be stars show a mix of Class A and Class B PAH profiles 'irrespective of Teff,' and associated reflection nebulae contribute Class A PAHs. Since Be stars populate the high-TeFF end of the sample, the blueward shift could reflect nebular contamination rather than stellar-temperature-driven processing. The authors should analyze the Herbig Ae and Herbig Be subsamples separately, or explicitly test whether the trend persists after removing stars with associated nebulosity.","section":"Section 3.6"},{"comment":"The PAH and silicate detection-frequency histograms are plotted without error bars or confidence intervals. With 10 bins of 12 stars each, binomial counting uncertainties are substantial (about +/-14 percent for a 50% detection fraction), so the claimed 'sweet spot' around 7000–11000 K and the high-temperature decline need a significance test or confidence intervals before they can be considered robust.","section":"Section 3.3, Figure 3"},{"comment":"The claim that the 6.0 micron feature does not shift is based only on a 3xMAD visual criterion, with no correlation coefficient or significance test relating the 6.0 micron peak position to the 6.2 micron peak position or to Teff. Since both peaks are extracted from the same double-Gaussian fit, correlated fitting errors are possible, and the 'distinct origin' conclusion requires a quantitative test.","section":"Section 3.6, Figure 5 right panel"}],"minor_comments":[{"comment":"The reported 11.0/11.2 micron ratios are internally inconsistent: Section 3.8 gives a range of 0.03–0.12 and cites HD 97300 and PDS 144N values of 0.16 and 0.14, while the Conclusion states a range of 0.046–0.12. These numbers should be reconciled.","section":"Section 3.8 and Section 4"},{"comment":"The caption describes the left panel as a transition 'from Class B/C (circumstellar) to Class A (nebular),' but the text defines Class C as low-mass T Tauri disk profiles; please clarify whether Class C profiles occur in the HAeBe sample or only in the comparison literature.","section":"Figure 5 caption"},{"comment":"There are several typographical issues: 'cubic spine' in Section 3.6, 'spilne anchor points' in Section 3.8, 'PSD 144N' in the Figure 4 caption, and duplicated 'stretching stretching' in Section 3.7.","section":"Throughout"},{"comment":"The Pearson correlation matrix is reported without p-values or corresponding sample sizes for each ratio pair; adding significance levels would help the reader assess the strength of the 0.9–0.99 correlations.","section":"Section 3.7, Figure 6"},{"comment":"Capitalization of 'CASSISjuice' is inconsistent (CASSISjuice vs CASSISJuice); please unify the spelling throughout the manuscript.","section":"Section 2.2"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the manuscript's main novelty is the 6.2 micron–Teff shift and its interpretation, but this central claim currently lacks quantitative uncertainties, significance testing, and robustness checks. The data are public and the required analysis is well within the paper's scope, so I recommend major revision rather than rejection. The classification catalog and uniform reduction are a useful community resource that should be preserved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a useful, workmanlike survey: 124 HAeBe stars with uniformly processed Spitzer IRS spectra, clearly classified, with public data and reproducible methods. The genuinely new items are the first systematic 11.0/11.2 µm ratio measurement in a HAeBe sample, the largest-sample test of the 6.2 µm peak shift with Teff (50 sources), and an empirical target-selection rule (n2-24 > -1) that will be handy for JWST proposals. The PAH detection frequency trends and the flared-disk preference match what Acke et al. 2010 and Seok & Li 2017 already showed, but consolidating them on a uniform sample is still worth having.\n\nThe paper's main quantitative claims, though, are under-supported. The 6.2 µm blueward shift with Teff is the headline result, but Fig. 5 is just a scatter plot: no peak-position uncertainties, no fitted slope or correlation coefficient, no significance test. That matters because at Spitzer IRS low-resolution sampling (R≈60–120) the wavelength bin is comparable to the quoted 6.19–6.28 µm Class A/B separation. The paper also never tests the stability of the fitted peak positions against alternative continuum anchor points or deblending choices, which is exactly the robustness check that would make the trend convincing. The authors do honestly note that Herbig Be stars show a mix of Class A and B irrespective of Teff and that reflection nebulae contribute Class A PAHs, but that is also a Teff-correlated environmental confound that no attempt is made to remove. So the largest-sample claim is not yet fully supported.\n\nThere are also some internal inconsistencies: the 11.0/11.2 analysis mentions 20 resolved sources in one place and 21 in another, and the reported ratio range in the text (0.03–0.12) differs from the conclusion (0.046–0.12), while two reference stars are quoted at 0.16 and 0.14. These need to be reconciled.\n\nNone of this is fatal. The survey is a solid basis for future work, and the public data let anyone re-check the numbers. A careful revision that adds uncertainties, significance tests, and a continuum-choice robustness test would substantially raise its value.\n\nI'd send it to peer review and would probably cite the n2-24 selection rule and detection statistics in my own work. It is not a paradigm changer, but it is a legitimate contribution.","headline":"Solid, useful HAeBe PAH survey with public data; the headline 6.2 µm–Teff shift is plausible but under-quantified at the current spectral resolution.","tokens_in":20812,"tokens_out":2567,"would_cite":true,"duration_ms":24136,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In 50 Herbig Ae/Be stars, the 6.2-micron PAH band shifts blueward as stellar temperature rises, and this is the largest sample yet used to test that peak shift.","keywords":["Herbig Ae/Be stars","polycyclic aromatic hydrocarbons","mid-infrared spectroscopy","Spitzer IRS","PAH peak shift","circumstellar disks","pre-main-sequence stars","PAH ionization"],"falsifier":"Re-fit the same 50 spectra with an independent continuum definition, such as a local linear baseline or a spectral decomposition that does not rely on spline anchor points, and test whether the 6.2 micron peak still moves blueward with effective temperature. If the trend disappears or weakens when the continuum choice changes, the claimed shift would be an artifact of the fitting procedure.","tokens_in":19719,"feed_emoji":"🔭","tokens_out":7218,"duration_ms":71344,"temperature":0.7,"pith_summary":"This paper presents a mid-infrared survey of 124 Herbig Ae/Be stars using newly processed Spitzer/IRS spectra, classifying each by its molecular and dust features into PAH-only, silicate-only, both, featureless, and HAC groups. It finds PAH emission in 64% of the sample and silicates in 50%, with detections peaking in the effective-temperature range 7000–11000 K. The central result is that the 6.2 micron PAH band, detected in 50 stars, shifts blueward as stellar effective temperature increases, and this is the largest sample yet used to test that peak shift. The adjacent 6.0 micron feature does not shift with it, which the authors take as evidence that the two bands arise from different molecular vibrations: C=O carbonyl or olefinic C=C stretching rather than aromatic C–C stretching. If correct, the shift maps how ultraviolet processing reshapes the PAH population in planet-forming disks and gives JWST a concrete target to check with spatially resolved spectroscopy.","feed_headline":"Hotter Herbig stars shift the 6.2-micron PAH band blueward","feed_subtitle":"The 6.0-micron band does not follow, hinting at a distinct carrier — a result JWST can test directly.","key_machinery":"The load-bearing object is the 6.2 micron PAH emission complex, produced by aromatic C–C stretching vibrations, together with the weaker 6.0 micron band blended on its blue side. The analysis isolates these bands by subtracting a cubic-spline continuum and deblending the two peaks with a double-Gaussian model, then compares the measured 6.2 micron peak position with stellar effective temperature and assigns it to Class A (6.19–6.235 microns), B (6.235–6.28 microns), or C (6.28–6.32 microns) profiles. The argument that carries the paper's main claim is the observed blueward progression of this peak across a 50-star sample, with the 6.0 micron feature serving as a counterexample that does not follow the same trend.","core_discovery":"The paper's central claim is that in Herbig Ae/Be stars the peak wavelength of the 6.2 micron PAH emission band is not fixed but moves blueward, from Class B/C toward Class A profiles, as the effective temperature of the central star rises, with 50 sources providing the largest sample yet used to test this correlation. The authors interpret the shift as progressive UV-driven chemical processing: hotter stars destroy aliphatic sidegroups, dehydrogenate PAHs, and fragment smaller molecules, leaving a more compact, stable aromatic population that emits at shorter wavelengths. The weaker 6.0 micron band, by contrast, stays near a median of 6.02 microns regardless of temperature, implying a distinct carrier that does not respond to the same processing. The paper also finds that PAH detections cluster in stars with spectral index n2−24 > −1 and flared disk morphologies, and reports what it describes as the first systematic measurement of the 11.0/11.2 micron PAH ratio in a Herbig Ae/Be sample, spanning 0.03–0.12 and indicating a range of ionization conditions.","pith_inferences":["If the 6.2 micron shift is real, JWST spectroscopy of the same targets should reproduce the blueward progression within spatially resolved disks, confirming it is an intrinsic processing trend rather than a sample-level artifact.","The stationary 6.0 micron feature offers a partial control: a continuum-fitting artifact that shifted both bands together would be ruled out, though an artifact affecting only the 6.2 micron blend remains possible.","A prospective test would observe the roughly half of known Herbig Ae/Be stars not yet covered by Spitzer, selecting targets with n2−24 > −1, and check whether PAH detection rates there exceed those in lower-index stars as predicted.","The 11.0/11.2 ratios overlapping PDR values suggest that disk PAHs can be as ionized as those in photodissociation regions, implying heavy processing; spatially resolved mid-infrared observations could separate disk from nebular PAH populations to test this."],"forward_implications":["Future PAH surveys of intermediate-mass pre-main-sequence stars can prioritize targets with n2−24 > −1 and flared morphologies, where PAH detection rates are highest.","PAH emission is most frequently detected in the 7000–11000 K range and declines at higher temperatures, consistent with ultraviolet destruction or chemical processing of PAHs in hotter environments.","The blueward 6.2 micron shift with stellar temperature supports a picture in which hotter radiation fields leave a smaller, more processed PAH population.","Because the 6.0 micron feature stays near 6.02 microns while 6.2 shifts, the two bands trace distinct molecular carriers, a distinction future high-resolution JWST spectra can test directly.","The measured 11.0/11.2 micron ratios, spanning 0.03–0.12, show that PAH ionization in Herbig Ae/Be disks overlaps PDR-like conditions and reaches higher fractions in flared disks such as PDS 144N."],"supporting_citations":[{"why":"Established PAH peak shifts with effective temperature in Herbig Ae stars and provides the earlier sample that this work extends to Herbig Be stars.","marker":"Acke et al. 2010"},{"why":"Defined the Class A/B/C PAH band profile classification and attributed the 6.2 and 6.0 micron features to specific molecular vibrations.","marker":"Peeters et al. 2002"},{"why":"Supplied the cubic-spline continuum subtraction method and earlier PAH detection statistics for Herbig stars.","marker":"Seok & Li 2017"},{"why":"Built the Spitzer Spectral HAeBe Catalog that is the foundation of this sample and identified the HAC and C60 sources used here.","marker":"Arun et al. 2023"},{"why":"Described the CASSIS pipeline and data products used to process the Spitzer/IRS spectra.","marker":"Lebouteiller et al. 2015"},{"why":"Provided DFT modeling showing how PAH size and ionization shift the 7.7 micron complex, the interpretive basis for the 6.2 micron shift.","marker":"Maragkoudakis et al. 2023"},{"why":"Supplied the Orion PDR JWST template used to compare the 6.0 micron feature's peak position.","marker":"Chown et al. 2024"},{"why":"Introduced the Group I/II flared versus self-shadowed disk classification used to link PAH detection to disk morphology.","marker":"Meeus et al. 2001"},{"why":"Defined the wavelength integration ranges used for PAH band ratios including 11.0/11.2 microns.","marker":"Zhang et al. 2024"}],"fun_headline_variants":["Hotter Herbig stars shift 6.2-micron PAH emission blueward","6.0-micron PAH stays fixed as 6.2-micron shifts blueward","Large sample of Herbig stars shows 6.2-micron PAH blue shift","6.0-micron PAH band implies distinct carrier in Herbig stars","First 11.0/11.2 PAH ratio survey in Herbig Ae/Be stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The measured 6.2 micron peak wavelengths are not systematically biased by the cubic-spline continuum subtraction and double-Gaussian deblending; if hotter stars have differently shaped dust continua that tilt the fitted peak, the blueward shift with temperature could be an artifact.","fun_headline_variants_meta":{"raw":{"variants":["Hotter Herbig stars shift 6.2-micron PAH emission blueward","6.0-micron PAH stays fixed as 6.2-micron shifts blueward","Large sample of Herbig stars shows 6.2-micron PAH blue shift","6.0-micron PAH band implies distinct carrier in Herbig stars","First 11.0/11.2 PAH ratio survey in Herbig Ae/Be stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002956,"raw_usage":{"total_tokens":11290,"prompt_tokens":1089,"completion_tokens":10201,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":705,"completion_tokens_details":{"reasoning_tokens":10082}},"tokens_in":705,"tokens_out":10201,"duration_ms":59552,"temperature":1.0,"reasoning_tokens":10082,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:18:17.393102+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the same 50 spectra with an independent continuum definition, such as a local linear baseline or a spectral decomposition that does not rely on spline anchor points, and test whether the 6.2 micron peak still moves blueward with effective temperature. If the trend disappears or weakens when the continuum choice changes, the claimed shift would be an artifact of the fitting procedure.","supporting_citations":[],"review_version":1}