{"id":"bb13c588-9e6f-4527-92dc-f04291916a39","arxiv_id":"2506.12974","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":8,"one_line_summary":"HS2 has been detected for the first time in a cold dark cloud, TMC-1, with an abundance of about 6e-11 relative to H2, inferred from five hyperfine lines.","lead":"Astronomers report the first detection of HS2, a molecule with two sulfur atoms, in the cold dark cloud TMC-1. The finding adds a new constraint to the long-standing puzzle of where sulfur hides in star-forming clouds.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reported detection significance is not supported by Table 1: all five HS2 lines are ~3σ, and the paper does not demonstrate consistency with the three nondetected hyperfine components.","rationale":"The paper's central claim is the first detection of HS2 in TMC-1, so the evidence establishing the reality of the lines is the load-bearing element. The most concrete problem is the mismatch between the significance asserted in Section 3 and the Gaussian-fit errors in Table 1: the table implies every line is only ~3σ, so the detection is weaker than advertised. This also contradicts the reader's weakest-assumption summary, which took the weakest line to be 3.5σ. The additional issue of the three nondetected hyperfine components is important because a self-consistent LTE identification should also explain why predicted lines are absent; the paper's blanket statement about higher noise is not quantified. Both issues are checkable from the existing data and model. A first detection by a single group based on five ~3σ lines should at least be vetted before being accepted at high confidence; hence the verdict should be conditional on these checks rather than unchanged. I agree with the reader only partially: line identification is indeed critical, but the specific vulnerability is the overstated significance and the missing nondetection analysis, not primarily the catalog frequencies.","tokens_in":10975,"tokens_out":22646,"duration_ms":230450,"concrete_test":"Re-measure the QUIJOTE spectra at the five detected and three nondetected HS2 frequencies: compute the line-free rms, the Gaussian-fit S/N per line, and the predicted T_A from the MADEX LTE model for all eight hyperfine components using the fitted Trot and column density. Confirm (1) whether any detected line has peak or integrated S/N below 4σ when the same noise measure is used as in the text; (2) whether any nondetected line has a predicted peak above 3σ times its local rms. If either answer is yes, the detection significance or the LTE identification needs to be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 3 claims four detected HS2 lines above 4σ and one at 3.5σ, but the 1σ Gaussian-fit errors in Table 1 give peak S/N of 0.60/0.22=2.7, 1.02/0.28=3.6, 0.47/0.17=2.8, 1.41/0.50=2.8, and 1.00/0.36=2.8; the integrated intensities give 2.8-3.7σ. Hence, by the errors in Table 1, no line individually meets the stated 4σ threshold and the weakest line is 2.7σ, not 3.5σ. The paper does not define the noise measure used for the >4σ claim. If the quoted fit errors are the true uncertainties, the first-detection claim depends on five marginal (~3σ) features. The self-consistency of the LTE identification is also not demonstrated: the three nondetected hyperfine components are dismissed as being in regions of higher noise without quantitative noise estimates or predicted intensities from the fitted model, so a predicted line could be present but missed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the first detection of the HS2 radical toward the cold dark cloud TMC-1 using the QUIJOTE Q-band survey. Five hyperfine components of the 20,2-10,1 and 30,3-20,2 rotational transitions are identified with the help of CDMS/MADEX spectroscopic data. An LTE model with Trot = 7 K and an assumed source size of 80 arcsec yields a column density of (5.7 ± 1.1) × 10^11 cm^-2 and an abundance of (5.7 ± 1.1) × 10^-11 relative to H2. The authors also run Nautilus chemical models with the KIDA network and attribute the main formation of HS2 to dissociative recombination of H2S2+ and H3S2+.","tokens_in":11231,"tokens_out":4690,"duration_ms":50976,"significance":"If the detection is secure, it would be the first identification of a molecule containing more than one sulfur atom in a cold dark cloud, and it would add a new observational constraint on sulfur chemistry in dense, quiescent gas. The comparison with HSO and the use of standard, publicly available tools (CDMS/MADEX and Nautilus/KIDA) are strengths. However, the detection claim currently rests on a significance statement that is contradicted by the tabulated line parameters, and the self-consistency of the hyperfine pattern is not demonstrated quantitatively. These issues are load-bearing for the central claim, so the paper needs major revision.","major_comments":[{"comment":"The claimed detection significance is not supported by the uncertainties in Table 1. From the tabulated peak antenna temperatures and their 1σ errors, the signal-to-noise ratios are 0.60/0.22 = 2.7, 1.02/0.28 = 3.6, 0.47/0.17 = 2.8, 1.41/0.50 = 2.8, and 1.00/0.36 = 2.8. The integrated values give 3.0, 3.7, 3.0, 2.8, and 3.3. None of the five lines reaches 4σ, and the weakest is 2.7σ, not the quoted 3.5σ. The statement in Section 3 that four lines are detected above 4σ and one at 3.5σ is therefore inconsistent with the table. Please define the precise noise measure and smoothing used for the significance estimate, or revise the detection claim. If the Table 1 errors are the true uncertainties, a statistical combination of the five lines (and the nondetected components) should be presented to justify a secure identification.","section":"Section 3, Table 1"},{"comment":"The three nondetected hyperfine components (20,2-10,1, J = 3/2-1/2, F = 1-0; 30,3-20,2, J = 7/2-5/2, F = 3-2; and 30,3-20,2, J = 5/2-3/2, F = 2-1) are dismissed as being weaker or falling in noisier regions, but no quantitative evidence is provided. To establish that the five detected lines form a self-consistent hyperfine pattern, give the local rms noise at each nondetected frequency, the expected peak intensity from the best-fit LTE model, and the resulting upper limit or detection. Without these numbers, the possibility that one or more of the five features are unrelated chance fluctuations or lines of a different carrier is not excluded.","section":"Section 3"},{"comment":"The quoted 20% uncertainty in the HS2 column density appears to reflect only the dispersion of the line fits and does not propagate the assumed source size of 80 arcsec, the assumed H2 column density of 10^22 cm^-2, or the range of plausible rotational temperatures. Since the abundance is directly compared with HSO and with the Horsehead value, please specify how the 20% error was computed and report the sensitivity of N(HS2) to the source-size and temperature assumptions. This does not affect the reality of the lines, but it is essential for the quantitative abundance comparison that is a central result of the paper.","section":"Section 3, LTE fit"}],"minor_comments":[{"comment":"There are several typographical errors, including 'dest uction' in Figure B.2, 'effotts' in Section 4.1, and the corrupted text 't/greaterorsimilar104 yr' in Section 4.2. Please correct these and check the final PDF for symbol rendering.","section":"Throughout"},{"comment":"The header for the integrated intensity column is difficult to parse; consider writing it as '∫ T_A^* dv' and using consistent units in the column caption.","section":"Table 1"},{"comment":"The sentence describing the FFTS spectral resolution says 'a spectral resolution of 38.15 kHz (~0.27 km/s) was also used' but the verb is missing; please rephrase.","section":"Section 2"}],"recommendation":"major_revision","confidential_remarks":"The discrepancy between the claimed >4σ detections and the ~3σ values in Table 1 is serious. If the authors cannot provide a valid noise definition that justifies the stated significances, or a statistically combined detection that accounts for all hyperfine components, I would not support publication of this as a first-detection Letter. The paper may be salvageable, but the central claim needs to be reworked and re-reviewed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the headline: the paper reports a genuinely new result — HS2 in a cold dark cloud, the first sulfur-sulfur molecule in TMC-1 — and the LTE modeling and chemical network analysis are useful. But the detection significance is overstated. Section 3 says four lines are above 4σ and one at 3.5σ; Table 1's own Gaussian errors give peak S/N between 2.7 and 3.6, and integrated intensities between 2.8 and 3.7. No line individually meets 4σ. That matters for a first detection. The stress-test note is correct on this point.\n\nWhat the paper does well: the observations come from a deep, well-calibrated survey; the line velocities and widths are consistent across five hyperfine components of two rotational transitions; the LTE fit with T_rot=7K is plausible and matches other sulfur species; the column density uncertainty is honestly stated at 20%; the chemical model uses public code and a standard network, and the authors explicitly flag missing surface and neutral reactions. The comparison with HSO and the sulfur depletion problem gives the result context.\n\nThe soft spots, in proportion: the significance issue is the main one. The nondetected hyperfine components are dismissed without quantitative noise estimates or predicted intensities, so a reader can't check whether a real line could be hidden. The assumed source size of 80 arcsec is reasonable but unverified, and the LTE assumption rests on a low T_rot without collisional rates. These are standard for the field, but the significance point is not a standard quibble — it's a number that should be corrected or re-derived.\n\nBottom line: the detection is likely real — five lines at consistent velocities from a known molecule don't usually all appear by chance — but the paper doesn't demonstrate that statistically. A referee should ask for a combined significance estimate (e.g., product of probabilities across lines given the known hyperfine pattern) and a table of predicted vs. observed intensities for all components, including the nondetected ones.\n\nWho it's for: astrochemists working on sulfur chemistry and dark cloud surveys. It deserves peer review, but with the significance analysis fixed. The central claim can probably survive, but not as it stands.","headline":"HS2 in TMC-1 is a plausible first detection, but the paper's >4σ claim doesn't survive Table 1; the lines are ~3σ, so the detection needs a proper combined significance.","tokens_in":11798,"tokens_out":2348,"would_cite":false,"duration_ms":24007,"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":"This paper reports the first detection of the HS2 radical in the cold dark cloud TMC-1, with an abundance about seven times that of its oxygen analogue HSO.","keywords":["HS2","TMC-1","cold dark cloud","sulphur chemistry","interstellar molecules","radio spectral line survey","dissociative recombination","sulphur depletion"],"falsifier":"Observe the three hyperfine components of HS2 that this survey did not detect ($2_{0,2}-1_{0,1}$, $J=3/2-1/2$, $F=1-0$, and the $F=3-2$ and $F=2-1$ components of $3_{0,3}-2_{0,2}$) with deeper integration: if the assignment is correct they must appear at their predicted relative intensities and at the same source velocity as the five detected lines, near 5.8–6.1 km s$^{-1}$.","tokens_in":10789,"feed_emoji":"🔭","tokens_out":10886,"duration_ms":108559,"temperature":0.7,"pith_summary":"This paper reports the first detection of the sulfur radical HS2 in a cold dark cloud, at the cyanopolyyne peak of TMC-1. The evidence consists of five faint radio lines, below 1.5 mK in antenna temperature, assigned to fine and hyperfine components of the $2_{0,2}-1_{0,1}$ and $3_{0,3}-2_{0,2}$ rotational transitions of HS2. A local thermodynamic equilibrium model with a rotational temperature of 7 K reproduces the observed profiles and gives a column density of $(5.7\\pm1.1)\\times10^{11}$ cm$^{-2}$, hence an abundance of $(5.7\\pm1.1)\\times10^{-11}$ relative to H$_2$, about seven times the abundance of the oxygen-containing counterpart HSO. A gas-grain chemical model reproduces this abundance and points to dissociative recombination of H$_3$S$_2^+$ and H$_2$S$_2^+$ as the main formation route. If correct, this is the first observation of a molecule with more than one sulfur atom in this type of source, and it gives a new handle on the long-standing question of where sulfur hides in dense clouds.","feed_headline":"First two-sulfur molecule detected in a cold dark cloud","feed_subtitle":"The HS2 radical turns up in TMC-1 at seven times the abundance of its oxygen twin HSO.","key_machinery":"The object that carries the claim is the HS2 radical, a bent molecule with an unpaired electron, which produces fine and hyperfine line structure that gives each rotational transition a characteristic multiplet. The argument works by matching five observed lines to the predicted multiplet positions and relative intensities of the $2_{0,2}-1_{0,1}$ and $3_{0,3}-2_{0,2}$ transitions, then using an LTE excitation model with rotational temperature and column density as free parameters to reproduce the full set of profiles. Because the upper levels lie at only 2.2 and 4.5 K, the derived column density is insensitive to the exact temperature as long as it is not below about 4.5 K. A three-phase gas-grain chemical model (gas, surface, mantle) with a standard astrochemical network supplies the formation and destruction rates used to interpret the observed abundance.","core_discovery":"On the paper's own terms, the discovery is that HS2 exists in the cold gas of TMC-1 with observable emission. Five lines, four above $4\\sigma$ and one at $3.5\\sigma$, match catalog rest frequencies and hyperfine patterns of HS2 and are not blended with known spectral features. Fitting all five line profiles with a local thermodynamic equilibrium model at $T_{\\rm rot}=7$ K yields $N({\\rm HS}_2)=(5.7\\pm1.1)\\times10^{11}$ cm$^{-2}$, and adopting $N({\\rm H}_2)=10^{22}$ cm$^{-2}$ gives a fractional abundance of $(5.7\\pm1.1)\\times10^{-11}$. The chemical model identifies dissociative recombination of H$_3$S$_2^+$ and H$_2$S$_2^+$ as the dominant gas-phase formation channels, with destruction by proton transfer from H$^+$, H$_3^+$, HCO$^+$, and H$_3$O$^+$, and finds no significant grain-surface contribution in the adopted network. The paper concludes that this is the first detection of a species containing more than one sulfur atom in a cold dark cloud.","pith_inferences":["A natural extension not developed in the paper is a direct search for the precursor ions H$_3$S$_2^+$ and H$_2$S$_2^+$ in the same Q-band data; detecting them would test the proposed formation pathway rather than only its product.","The chemical network omits the grain-surface reaction S$_2$ + H $\\rightarrow$ HS$_2$ while treating H$_2$ + S$_2$ as barrierless; updating that choice could raise or lower the predicted HS$_2$ abundance at early times and is a concrete way to test the assumed surface chemistry.","If the HS2/HSO ratio reflects a general chalcogen pattern, cold clouds should also show measurable differences between HO2 and HSO-type chemistries; better abundance limits on HO2 in TMC-1 would show whether sulfur and oxygen really diverge in this way."],"forward_implications":["HS2 becomes the first two-sulfur molecule found in a cold dark cloud, so searches for related di-sulfur species such as S2, H2S2, and HSSH in TMC-1 become well motivated.","The HS2/HSO abundance ratio of about seven in TMC-1 means the sulfur analogue is not suppressed relative to its oxygen counterpart, a constraint for chemical networks that treat S and O similarly.","The model's ability to reproduce the observed abundance with a depleted initial sulfur abundance of $1.5\\times10^{-6}$ supports the idea that HS2 forms mainly via gas-phase ion-electron recombination rather than requiring new grain-surface reactions.","Because the modeled HS2 abundance depends strongly on cloud density and on the sulfur depletion factor, future detections in other dark clouds could be used as probes of sulfur depletion.","The very low upper-level energies of the detected transitions imply that HS2 emission traces the cold envelope of TMC-1, consistent with other sulfur molecules in the same source."],"supporting_citations":[{"why":"Defines the Q-band line survey and the calibration and folding procedures from which the five detected lines are drawn.","marker":"Cernicharo et al. (2021a)"},{"why":"Supplies the catalog spectrum of HS2 used to predict the line positions and hyperfine patterns.","marker":"Müller et al. (2001)"},{"why":"Documents the updated database from which the HS2 rest frequencies used in the fit are taken.","marker":"Endres et al. (2016)"},{"why":"Describes the excitation code used to model the line profiles under LTE and derive the column density.","marker":"Cernicharo (2012)"},{"why":"Reports the previous HS2 detection in the Horsehead PDR, the benchmark against which this cold-cloud detection is compared.","marker":"Fuente et al. (2017)"},{"why":"Provides the HSO abundance in TMC-1 used for the HS2/HSO abundance comparison.","marker":"Marcelino et al. (2023)"},{"why":"Gives the H2 column density of TMC-1 adopted to convert the HS2 column density to a fractional abundance.","marker":"Cernicharo & Guelin (1987)"},{"why":"Describes the three-phase gas-grain chemical model used to compute formation and destruction rates for HS2.","marker":"Ruaud et al. (2016)"},{"why":"Supplies the updated astrochemical network used inside the gas-grain model.","marker":"Wakelam et al. (2024)"}],"fun_headline_variants":["First multi-sulfur molecule found in a cold dark cloud","HS2 detected in TMC-1: a second sulfur atom joins in","TMC-1 reveals HS2, a sulfur pair in cold space","HS2 in TMC-1 outshines its oxygen twin by 7x","First sighting of HS2 in a dark cloud's cold gas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The identification hinges on the catalog rest frequencies and hyperfine intensities of HS2 and on the five lines being genuine rather than chance coincidences with unknown carriers; the weakest detected line is only at $3.5\\sigma$.","fun_headline_variants_meta":{"raw":{"variants":["First multi-sulfur molecule found in a cold dark cloud","HS2 detected in TMC-1: a second sulfur atom joins in","TMC-1 reveals HS2, a sulfur pair in cold space","HS2 in TMC-1 outshines its oxygen twin by 7x","First sighting of HS2 in a dark cloud's cold gas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000322,"raw_usage":{"total_tokens":1846,"prompt_tokens":1014,"completion_tokens":832,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":746}},"tokens_in":630,"tokens_out":832,"duration_ms":8175,"temperature":1.0,"reasoning_tokens":746,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:37:10.650249+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the three hyperfine components of HS2 that this survey did not detect ($2_{0,2}-1_{0,1}$, $J=3/2-1/2$, $F=1-0$, and the $F=3-2$ and $F=2-1$ components of $3_{0,3}-2_{0,2}$) with deeper integration: if the assignment is correct they must appear at their predicted relative intensities and at the same source velocity as the five detected lines, near 5.8–6.1 km s$^{-1}$.","supporting_citations":[{"cited_title":"2012, in EAS Publications Series, Vol","cited_arxiv_id":null,"evidence_quote":"Describes the excitation code used to model the line profiles under LTE and derive the column density."},{"cited_title":"R., Pety , J., et al","cited_arxiv_id":null,"evidence_quote":"Reports the previous HS2 detection in the Horsehead PDR, the benchmark against which this cold-cloud detection is compared."},{"cited_title":"2023, , 674, L13","cited_arxiv_id":null,"evidence_quote":"Provides the HSO abundance in TMC-1 used for the HS2/HSO abundance comparison."},{"cited_title":"& Guelin , M","cited_arxiv_id":null,"evidence_quote":"Gives the H2 column density of TMC-1 adopted to convert the HS2 column density to a fractional abundance."}],"review_version":1}