{"id":"d9ca452a-0f6e-4475-99f0-b2b2dda07c8a","arxiv_id":"2506.20468","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"JWST spectra of two PDRs show aliphatic carbon bonds are destroyed at cloud surfaces while aromatic bonds survive, and deuterium substitution is more efficient in aliphatic than aromatic carbon bonds.","lead":"This paper presents JWST spectroscopic maps of two nearby star-forming clouds and lists the hundreds of gas and dust spectral features seen across their surfaces. It shows that the dusty outflow from the Horsehead cloud edge contains only aromatic carbon while aliphatic (chain-like) carbon bonds are destroyed, and that deuterium substitution happens more readily in aliphatic than aromatic carbon bonds.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative D/H ratios rest on PAH-specific theoretical band strengths and on 4.4 μm measurements that overlap the 13CO2 ice feature; neither uncertainty is propagated into the abstract's headline numbers.","rationale":"I read the paper as primarily a spectral inventory plus a set of new PDR conclusions: outflow composition, aliphatic processing, and deuterium fractionation. The qualitative statements—3.3 μm only in the Horsehead outflow, 3.4/3.3 decline across the front, and no 4.4 μm feature in the Horsehead—are supported by the displayed spectra and spatial maps, and the MIRI F770W outflow caveat is honestly stated. The weakest link is the quantitative D/H conversion, because the abstract's main numerical claim inherits two unquantified systematics. First, the adopted intrinsic band strengths are for PAH molecules with aliphatic side chains; the paper deliberately avoids identifying the 3–5 μm carriers with PAHs, so applying these values to a-C(:H)-type carriers is an unsupported extrapolation. Second, the 4.4 μm aromatic C-D band overlaps the 13CO2 ice feature at ~4.39 μm in the very regions (NGC 7023 MOL, DF3) used for the shielded aromatic D/H points; no correction or caveat is given in Table 6. Both issues are addressable and do not necessarily overturn the qualitative trends, but they directly affect the abstract numbers. This is essentially the same core assumption the reader flagged, though I add the 13CO2/4.4 μm overlap as a concrete internal inconsistency. I therefore see no reason to change the CONDITIONAL verdict.","tokens_in":37282,"tokens_out":6106,"duration_ms":68484,"concrete_test":"Recompute Table 6, columns 4 and 7, replacing the Yang et al. (2020)/Yang & Li (2023b) PAH-based intrinsic band strengths with an independent set appropriate to disordered carbonaceous carriers (e.g., laboratory-measured a-C:H/a-C:D band strengths), and, before measuring the 4.4 μm feature, subtract a fitted 13CO2 ice absorption profile from the NGC 7023 MOL and DF3 spectra. Then check whether the aliphatic-over-aromatic ND/NH ordering and the quoted 0.1–0.3 versus 0.03 ranges survive; if not, the abstract's quantitative deuterium-fraction claims require revision or rephrasing as limits.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—ND/NH ~ 0.1–0.3 for aliphatic C-D bonds versus ~0.03 for aromatic C-D bonds—is obtained in Sect. 4.5.3 by converting observed 3.4/4.7 and 3.3/4.4 μm intensity ratios with intrinsic C-D/C-H band strengths from Yang et al. (2020) and Yang & Li (2023b). Those calculations assume an underlying PAH structure with aliphatic side groups, yet Sect. 1 of this paper explicitly declines to identify the carriers as PAHs and leaves open disordered a-C(:H) carbonaceous material. If the true carriers are not PAH-like, the adopted band strengths do not apply, and the absolute D/H values—and possibly the claimed aliphatic-over-aromatic preference—change. Additionally, the 4.4 μm aromatic C-D band is centered at 4.38 μm in the NGC 7023 MOL and DF3 templates, while Sect. 4.4 reports 13CO2 ice absorption at ~4.39 μm in exactly those templates; the Table 6 extraction does not describe separating these components, so the shielded-region aromatic D/H points carry a specific, unmodeled blend risk. Neither the carrier-model uncertainty nor the ice-blend risk is propagated into the abstract's headline numbers.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents an overview of JWST GTO 1192 NIRSpec and MIRI IFU observations of the Horsehead and NGC 7023 NW PDRs, including data reduction, cross-calibration, region extraction, line identification, and an inventory of atomic, H2, CH+, CO, and ice features. The analysis focuses on the 3-5 μm carbonaceous emission: it confirms the entrainment of 3.3 μm aromatic carriers in the Horsehead photo-evaporative outflow, finds no 3.4 μm aliphatic emission there, and interprets this as complete removal of aliphatic bonds. It also reports detections of the 4.4 and 4.7 μm features attributed to aromatic and aliphatic C-D stretching modes and derives ND/NH ~0.1-0.3 for aliphatics versus ~0.03 for aromatics using theoretical band strengths from Yang et al. (2020) and Yang & Li (2023b).","tokens_in":37557,"tokens_out":6548,"duration_ms":71079,"significance":"The dataset and the public line and spectrum tables are a valuable community resource; the cross-calibration against JWST imaging is careful, and the per-spaxel treatment of the 3.3/3.4 μm bands demonstrates spatial resolution of hydrocarbon processing at the ~100 au scale. If the aliphatic destruction and deuterium fractionation claims hold, they provide new, spatially resolved constraints on carbonaceous dust evolution in PDRs and on deuterium incorporation into C-H bonds. The paper is appropriately cautious in the body about carrier identification (Sect. 1) and about the plateau origin, but the abstract's headline numbers inherit external model dependencies and blend uncertainties that are not fully propagated.","major_comments":[{"comment":"The abstract's statement that the outflow is 'purely aromatic' and that there is 'complete removal of aliphatic bonds' is based on scaling the DF1 spectrum by a hand-selected factor of 0.28 to match the 3.3 μm feature and then visually asserting that 3.4 μm emission is absent. No quantitative upper limit on I(3.4) in the H ii region is provided, and the uncertainty in the 0.28 scaling factor is not propagated. A formal non-detection limit, together with a discussion of whether size-selective entrainment could produce the same 3.3-only appearance, is needed before 'complete removal' can be stated.","section":"§4.5.1, Fig. 15"},{"comment":"The ND/NH values in Table 6 are converted from observed intensity ratios using intrinsic C-D/C-H band strengths from Yang et al. (2020) and Yang & Li (2023b), which assume an underlying PAH structure with aliphatic side groups. Section 1 explicitly declines to identify the carriers as PAHs and leaves disordered a-C(:H) materials as an open possibility. If the true carriers are not PAH-like, the adopted band strengths do not apply, and the absolute D/H ratios, and possibly the claimed aliphatic-over-aromatic preference, change. The paper should either justify the PAH-based band strengths for these specific PDR carriers or present the D/H numbers as model-dependent and temper the abstract accordingly.","section":"§4.5.3, Table 6, Sect. 1"},{"comment":"The 4.4 μm aromatic C-D feature used for the NGC 7023 MOL and DF3 entries in Table 6 (central wavelength ~4.38 μm, width ~0.1 μm) overlaps the 13CO2 ice absorption at ~4.39 μm identified in exactly those regions in Sect. 4.4. The extraction described in Sect. 4.5.3 does not state how this blend was separated, so the shielded-region aromatic ND/NH values (0.028 and 0.027) may be contaminated by the ice feature. A decomposition of the 13CO2 component, or an explicit uncertainty term for this blend, must be added before those points can support the conclusion that aromatic deuteration is less efficient.","section":"§4.4, §4.5.3, Fig. 12, Table 6"},{"comment":"The suggestion of a reduced aromatic deuteration efficiency in the Horsehead relative to NGC 7023 rests on a non-detection of the 4.4 μm feature in all Horsehead regions, combined with an assumed detection threshold. Because the 4.4 μm feature is weak even in NGC 7023 (0.11-0.48 in the units of Table 6, with 1σ errors of 0.02-0.13), the Horsehead absence should be reported as a quantitative upper limit with an explicitly stated detection threshold rather than a qualitative efficiency difference.","section":"§4.5.3"}],"minor_comments":[{"comment":"The term 'Disassociation Regions' should read 'Dissociation Regions'.","section":"Abstract"},{"comment":"The caption begins with 'Extractions regions on images'; this should be 'Extraction regions on images'.","section":"Fig. 1 caption"},{"comment":"The label 'HIIb' in the Horsehead sections is typographically inconsistent with the 'H ii' notation used throughout the text and should be unified.","section":"Tables 2 and 3"},{"comment":"The text states that the aliphatic fraction is roughly constant at ~0.15-0.25 in both objects, but Table 6 contains Horsehead DF2 at 0.274±0.049 and NGC 7023 DF3 at 0.233±0.041; the quoted range should be reconciled with these values or explicitly described as a typical range.","section":"§4.5.2, Table 6"},{"comment":"The symbol Cλ is used both for extinction in Eq. (2) and for continuum flux in Eq. (A.1); a different symbol for one of the two quantities would remove ambiguity.","section":"Eq. (2)"}],"recommendation":"major_revision","confidential_remarks":"The D/H derivation is not circular: it depends on external theoretical band strengths, not on the paper's own outputs. The main revision burden is on quantifying the two headline claims (the aliphatic-free outflow and the D/H ratios), including the 13CO2 ice blend in the 4.4 μm extraction. The manuscript is otherwise well within the scope of A&A and the underlying dataset is genuinely valuable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Best quick take: this is a data-rich paper that will get used. The JWST IFU mosaics of the Horsehead and NGC 7023 NW PDR deliver ~100 AU resolution on the 3.3/3.4 um aliphatic-aromatic transition, a clean detection of the 4.4 and 4.7 um deuterated bands, and a thorough line inventory. The line lists and filter contribution tables are exactly the kind of community resource that makes follow-up work easier. Credit where due: the per-spaxel maps showing the 3.4 um emission dropping to zero just inside the PDR front while 3.3 um persists, and the absence of 3.4 um in the outflow, are a clear visual demonstration of preferential aliphatic destruction. That part is convincing.\n\nThe soft spots are real but addressable. First, 'complete removal of aliphatic bonds' is an inference from a non-detection after scaling DF1 by a hand-picked factor of 0.28. That scaling is reasonable, and the spaxel gradient supports it, but a null detection is a null detection; the abstract should say 'no aliphatic emission detected' rather than claiming complete removal as fact. Second, the ND/NH ratios in Table 6 are converted using intrinsic band strengths from Yang et al. (2020) and Yang & Li (2023b), which assume PAH-like carriers with aliphatic side groups. The paper itself in Section 1 declines to commit to PAHs and leaves open a-C(:H) disorder. That is an internal tension; the quoted D/H numbers (0.1-0.3 aliphatic, 0.03 aromatic) are therefore model-dependent, and the error bars shown are statistical only. The systematic uncertainty from band strengths should be stated prominently and reflected in the headline numbers or removed from the abstract. Third, in MOL and DF3, the 4.4 um aromatic C-D band sits right on the 13CO2 ice absorption at 4.39 um; the paper never discusses that blend in Section 4.5.3. That could bias those two points.\n\nNone of this kills the paper. The qualitative conclusions look solid, and the dataset is a significant step forward. The right fix is calibration in the abstract and a paragraph acknowledging the systematic and blend issues. This deserves peer review; a good referee will push for those caveats, but the core result should survive.","headline":"New JWST PDR data of real value, but the deuterium fractions and 'complete removal' claim carry more model dependence than the abstract admits.","tokens_in":38158,"tokens_out":5237,"would_cite":true,"duration_ms":52199,"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":"JWST spectra of two star-forming clouds show deuterium preferentially replaces hydrogen in aliphatic carbon bonds, at 3–10 times the aromatic rate.","keywords":["photo-dissociation regions","JWST spectroscopy","carbonaceous dust","aliphatic C-H bonds","aromatic C-H bonds","deuterium fractionation","Horsehead nebula","NGC 7023"],"falsifier":"A laboratory measurement of the intrinsic C–D stretch band strength for the actual aliphatic carrier (for example hydrogenated amorphous carbon with alkane side chains) at interstellar temperatures, combined with a high-resolution spectrum that separates the 4.7 µm feature from CO ro-vibrational lines and CO ice absorption, would settle whether the reported $N_D/N_H \\sim 0.1$–$0.3$ is real or an artifact of band-strength and blending assumptions.","tokens_in":37096,"feed_emoji":"🔭","tokens_out":6488,"duration_ms":61030,"temperature":0.7,"pith_summary":"This paper presents a spatially resolved spectral inventory of two benchmark photo-dissociation regions, the Horsehead nebula and the NW filament of NGC 7023, using JWST IFU spectroscopy from 0.97 to 28 µm at physical resolutions down to about 50 AU. Its central new result is an asymmetry in how carbonaceous dust responds to ultraviolet radiation: aliphatic (chain-like) carbon–hydrogen bonds are destroyed faster than aromatic (ring-like) bonds at the cloud edge, yet where aliphatic bonds survive, deuterium substitution into them is far more efficient, with $N_D/N_H \\sim 0.1$–$0.3$ versus $\\sim 0.03$ for aromatic bonds. The paper also confirms spectroscopically that the photo-evaporative outflow from the Horsehead PDR carries only aromatic carriers, with no detectable aliphatic 3.4 µm emission. If correct, these findings give observers a direct infrared tracer of deuterium fractionation in dust and a sharp spatial diagnostic of UV processing across PDR fronts.","feed_headline":"Deuterium locks into aliphatic dust bonds 3-10x more than aromatic","feed_subtitle":"New JWST maps show UV strips chain-like carbon bonds first; deuterium marks the survivors.","key_machinery":"The load-bearing object is the 3–5 µm spectral complex of carbonaceous emission: the aromatic C–H stretch at 3.3 µm, the aliphatic C–H stretch at 3.4 µm with its red plateau, and their deuterated counterparts at roughly 4.4 µm (aromatic C–D) and 4.7 µm (aliphatic C–D). The argument proceeds by decomposing region-averaged and per-spaxel spectra into a polynomial continuum plus Drude profiles, mapping the 3.4/3.3 ratio across the front at roughly 100 AU resolution, and converting the observed 4.4/3.3 and 4.7/3.4 intensity ratios into $N_D/N_H$ using theoretical intrinsic band strengths for PAH-like carriers with aliphatic side groups.","core_discovery":"The paper establishes, from JWST NIRSpec and MIRI IFU spectra of the Horsehead and NGC 7023 PDRs, that aliphatic and aromatic carbonaceous carriers respond differently to the UV field and that deuterium uptake marks the difference. At the Horsehead front, the 3.3 µm aromatic C–H feature continues into the H ii region while the 3.4 µm aliphatic feature drops to zero just inside the front, showing that the outflow is purely aromatic and that exposure to even moderate UV removes aliphatic bonds. In the more shielded regions of both objects, the aliphatic-to-aromatic bond ratio is roughly constant at $\\sim 0.2$, but the deuterated fractions differ systematically: $N_D/N_H \\sim 0.1$–$0.3$ for aliphatic bonds versus $\\sim 0.03$ for aromatic bonds, and no 4.4 µm aromatic C–D feature is detected in the Horsehead at all. The absence of gas-phase HD lines, together with the high deuterium fraction in the carbonaceous bonds, is presented as consistent with the energetics by which deuterium replaces H in C–H bonds and catalyzes H$_2$ formation.","pith_inferences":["If the aliphatic C–D bond is as robust as the energetics suggest, deuterated aliphatic carriers could survive farther into H ii regions than their hydrogenated counterparts, so the 4.7 µm band may trace processed dust even where the 3.4 µm band has vanished.","The absence of the 4.4 µm aromatic C–D feature in the Horsehead, at a column where NGC 7023 shows it, hints that aromatic deuteration depends on the hardness of the radiation field or on the formation history of the carriers; a test would be to observe a PDR with an intermediate UV hardness.","The reported correlation of the 3.4–3.6 µm plateau with the aliphatic rather than the aromatic band, opposite to an earlier study, suggests the plateau carrier may be a separate aliphatic-rich population; separating its spatial profile from the main 3.4 µm band would settle the assignment.","The same 4.7 µm spectral window contains CO ro-vibrational lines and CO ice absorption, so future higher-resolution observations that resolve the C–D band from CO would directly test whether part of the inferred aliphatic deuterium fraction is actually CO emission."],"forward_implications":["Aliphatic C–H bonds act as a sensitive UV dosimeter: their disappearance at the PDR front marks the transition from shielded to unshielded material before the aromatic emission fades.","The 3.4 µm aliphatic feature can be used to isolate the photo-evaporative component of PDR outflows; its absence identifies outflows that have been stripped of aliphatic carriers.","Deuterium fractionation in carbonaceous dust is observable in the infrared: the 4.4 and 4.7 µm C–D bands offer a way to map $N_D/N_H$ in PDRs and, by extension, to test models of interstellar deuterium chemistry.","The measured line-contribution fractions in the JWST NIRCam and MIRI filters allow imaging-only studies to correct for line contamination in PDRs with similar physical conditions.","Ice formation (H$_2$O, CO$_2$, CO) begins in the deeper, more shielded regions of NGC 7023, providing a spatial link between the destruction of aliphatic bonds at the front and the freeze-out of volatiles behind it."],"supporting_citations":[{"why":"Imaging discovery of the photo-evaporative dust outflow in the Horsehead that this paper confirms spectroscopically and uses to define the PDR front and scaling factors.","marker":"Abergel et al. (2024)"},{"why":"Supplies the theoretical intrinsic band strengths used to convert observed 4.4/3.3 and 4.7/3.4 µm intensity ratios into $N_D/N_H$.","marker":"Yang et al. (2020)"},{"why":"Provides the intrinsic band strengths for the aromatic and aliphatic C–H and C–D stretches and the assumption of PAH carriers with aliphatic side groups.","marker":"Yang & Li (2023b)"},{"why":"Gives the energetics argument that deuterium replaces H in C–H bonds and catalyzes H$_2$ formation, used to explain the high deuterium fraction and the absence of HD lines.","marker":"Draine (2006)"},{"why":"Laboratory results cited for the greater UV stability of aromatic relative to aliphatic structures.","marker":"Sandford et al. (2013)"},{"why":"Laboratory and theoretical evidence that aromatic structures are stable relative to aliphatic ones under UV irradiation.","marker":"Marciniak et al. (2021)"},{"why":"Attribution of the 3.3 and 3.4 µm bands to aromatic and aliphatic C–H stretches, the foundation of the band inventory.","marker":"Tielens (2008)"},{"why":"Supplies the PDRS4ALL line list used for line identification and the NIRSpec flux scale comparison.","marker":"Peeters et al. (2024)"}],"fun_headline_variants":["JWST finds deuterium favors aliphatic dust bonds 10x over aromatic","Aliphatic bonds vanish in UV outflow, deuterium marks survivors","Deuterium clues: JWST shows aromatic dust survives UV, aliphatic does not"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The $N_D/N_H$ ratios in Table 6 assume that the 4.4 and 4.7 µm bands are C–D stretching modes and that the theoretical intrinsic band strengths, computed for PAH-like carriers with aliphatic side groups, apply to the actual interstellar carriers; if those band strengths or identifications are systematically wrong, the quantitative ratios change or collapse.","fun_headline_variants_meta":{"raw":{"variants":["JWST finds deuterium favors aliphatic dust bonds 10x over aromatic","Aliphatic bonds vanish in UV outflow, deuterium marks survivors","Deuterium clues: JWST shows aromatic dust survives UV, aliphatic does not"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000575,"raw_usage":{"total_tokens":2835,"prompt_tokens":1185,"completion_tokens":1650,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":801,"completion_tokens_details":{"reasoning_tokens":1586}},"tokens_in":801,"tokens_out":1650,"duration_ms":11182,"temperature":1.0,"reasoning_tokens":1586,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:47:43.124411+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A laboratory measurement of the intrinsic C–D stretch band strength for the actual aliphatic carrier (for example hydrogenated amorphous carbon with alkane side chains) at interstellar temperatures, combined with a high-resolution spectrum that separates the 4.7 µm feature from CO ro-vibrational lines and CO ice absorption, would settle whether the reported $N_D/N_H \\sim 0.1$–$0.3$ is real or an artifact of band-strength and blending assumptions.","supporting_citations":[{"cited_title":"D., et al","cited_arxiv_id":null,"evidence_quote":"Imaging discovery of the photo-evaporative dust outflow in the Horsehead that this paper confirms spectroscopically and uses to define the PDR front and scaling factors."},{"cited_title":"J., Li , A., & Glaser , R","cited_arxiv_id":null,"evidence_quote":"Supplies the theoretical intrinsic band strengths used to convert observed 4.4/3.3 and 4.7/3.4 µm intensity ratios into $N_D/N_H$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the energetics argument that deuterium replaces H in C–H bonds and catalyzes H$_2$ formation, used to explain the high deuterium fraction and the absence of HD lines."}],"review_version":1}