{"id":"68d561aa-4320-4756-9bfb-a5e8beb30e0a","arxiv_id":"2509.01834","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"PRIMA's FIRESS/FTM spectrometer is predicted to measure D/H in bright comets in ~6 hours, enabling comet and meteorite isotope samples of comparable size for the first time.","lead":"This paper tests whether the proposed PRIMA space telescope could measure the deuterium-to-hydrogen ratio in comet water, a key fingerprint for the origin of Earth's oceans. Its model predicts PRIMA can detect the needed HDO lines in bright comets in about six hours, opening the door to a large sample of comet isotope measurements.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 2 HDO line fluxes are model-only and unvalidated; a factor-2 excitation error turns the 6 h detection claim into ~25 h and shrinks the accessible comet sample.","rationale":"The reader's weakest-assumption already identifies the non-LTE model flux as the crux. My read agrees and sharpens it: the model predictions for the specific FIR HDO lines are not calibrated against any observation, and the two parameters most likely to affect them (HDO photolysis rate and outflow velocity) are fixed at assumed values. The 184.44 µm line has upper energy 144 K, so it is not anchored by the only good cometary HDO dataset, the Herschel/HIFI 509 GHz ground-state line. Since the paper's quantitative claims scale as the inverse square of these fluxes, even a factor-2 model uncertainty is material. The rest of the paper is careful: line-blending analysis uses conservative high abundances, line-to-continuum estimates are grounded in PACS data, and the instrument sensitivity is the stated CBE. The conditional verdict remains appropriate; no change to the reader's recommendation is needed, but the paper would be materially strengthened by a validation/error-budget section on the model fluxes.","tokens_in":16910,"tokens_out":14235,"duration_ms":166340,"concrete_test":"For a comet with a measured D/H and a Herschel/HIFI detection of the 509 GHz HDO line (e.g., 103P/Hartley 2), run the same non-LTE code used for Table 2 and require it to reproduce the 509 GHz line within its error bars. Then read off the predicted 184.44 µm flux and repeat with the HDO photolysis rate varied over the range allowed by available photodissociation cross-section data (e.g., Huebner & Mukherjee 2015) and outflow velocity varied from 0.6 to 1.0 km/s. If the 184.44 µm flux varies by more than a factor of 1.5 across this parameter range, the paper should quote an error budget on Table 2 before the 6 h / 'up to a dozen' feasibility claim is accepted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The feasibility argument (Section 6) is built on the four non-LTE HDO line fluxes in Table 2. These fluxes are produced by a model (Section 3.1) that fixes the outflow velocity at 0.8 km/s and assumes the HDO photolysis rate equals water's, but no error bars are quoted and no cometary observation validates the two preferred lines. The 184.44 µm and 234.64 µm lines have upper-level energies 144.4 K and 83.6 K, whereas the only space-based HDO line with a good cometary calibration, Herschel/HIFI 509 GHz, has upper energy 46.8 K; it therefore constrains a different excitation regime. Because integration time scales as flux^-2, a factor-2 overprediction in the model changes the headline 6 h at S/N=5 to ~25 h, and the Section 4 estimate of 'up to a dozen' comets shrinks correspondingly. The model is not circular or internally inconsistent—it is an established code—but the missing validation and missing uncertainty on the predicted fluxes are the load-bearing gap in the paper's central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a feasibility study for measuring the D/H ratio in cometary water with the proposed PRIMA space mission and its FIRESS/FTM instrument. Using a non-LTE excitation model for a reference comet (Q = 2e28 s^-1, r_h = Δ = 1 au, D/H = 2 VSMOW), it predicts HDO line fluxes at 125.85, 130.84, 184.44, and 234.64 μm, and finds that the 184.44 and 234.64 μm lines can be detected at S/N = 5 in about 6 hours. The paper then examines two confounding factors: blending with other coma species (using PSG LTE models with conservative abundances) and line-to-continuum ratio (using Herschel/PACS continuum measurements), concluding that neither is a blocking issue. It further estimates that up to a dozen comets with water production rate FOM = 1e28 s^-1 should be accessible during the 5-year PRIMA primary mission after accounting for the Sun avoidance angle, and that radial studies of a Garradd-like comet are possible out to ~2.7 au. The central claim is that PRIMA can build a comet D/H sample comparable in size to the meteorite sample, far exceeding the current four accurate space-based measurements.","tokens_in":17176,"tokens_out":6568,"duration_ms":72712,"significance":"If the detection-time estimates and target counts are robust, the proposed observations would provide a step-change in comet D/H statistics, enabling tests of reservoir differences, hyperactivity correlations, and the origin of Earth's water. The paper has clear strengths: it adopts a well-specified reference model, uses instrument current-best-estimate sensitivities, propagates the two main observational confounders (line blending and line-to-continuum) with conservative assumptions, and presents the scaling transparently. The non-LTE model is an established code (Cordiner et al. 2022), and the D/H ratio is used as an input, not derived, so there is no internal circularity. The main weakness is that the predicted HDO line fluxes are unvalidated for the specific transitions used, with no quoted uncertainty; because integration time scales as flux^-2, the headline 6-hour detection time and the accessible-comet sample are sensitive to a factor-2 excitation error.","major_comments":[{"comment":"The central detection-time estimates rest entirely on the four non-LTE HDO line fluxes in Table 2. The model fixes the outflow velocity at 0.8 km/s and assumes the HDO photolysis rate equals that of H2O, but no error bars are given and no cometary observation is shown validating the 184.44 and 234.64 μm lines. The only space-based HDO line with a good cometary calibration, Herschel/HIFI 509 GHz, has an upper-level energy of 46.8 K, whereas the 184.44 and 234.64 μm lines have Eu = 144.4 K and 83.6 K, respectively, so it samples a different excitation regime. Because integration time scales as flux^-2, a factor-2 overprediction of the model changes the 6-hour claim to ~25 hours and shrinks the accessible sample correspondingly. I request a quantitative uncertainty estimate for the predicted fluxes (e.g., from model parameter variations or a comparison with multi-transition observations of","section":"Section 3.1, Table 2"},{"comment":"The accessible-comet count is derived from the Origins Report (Ref. 23) and a 1/3 solar-elongation factor, but it is not explicitly tied to the detection model. The text states that 17±2 comets are accessible at FOM = 2e28 and that the number increases to about 36 at FOM = 1e28, then reduces to 'up to a dozen' after Sun-avoidance. However, the Table 2 detection times are computed for FOM = 2e28, and the conclusion only states in passing that 'integrations ~4 times longer should be feasible.' Please state explicitly which FOM threshold is required for a usable D/H measurement (including S/N, line choice, and continuum noise), and propagate the uncertainty in the comet flux distribution through to the expected number of actual measurements. As written, the target-count claim is stronger than the error budget presented.","section":"Section 4"},{"comment":"The line-to-continuum analysis uses Herschel/PACS continuum fluxes scaled as S ~ Q Δ^-1 r_h^-0.5, and the FIRESS/FTM sensitivity shown in Figure 2 is quoted 'w/o shot noise from the continuum.' The conclusion that the line-to-continuum ratio is not a confounding factor for the 184.44 and 234.64 μm lines assumes that the CBE sensitivity still applies when the continuum contributes at the 4-10% level relative to the line. Please state explicitly whether the quoted sensitivities include photon noise from the expected continuum level at these wavelengths, and discuss the systematic uncertainty in the continuum scaling exponent used to compare comets with different dust-to-gas ratios.","section":"Section 3.3, Table 5"}],"minor_comments":[{"comment":"The first sentence of Section 6 says 'at a distance r_h = Δ = 1' without units; it should be '1 au'.","section":"Section 6"},{"comment":"There is a typo: 'water production production rate' should be 'water production rate'.","section":"Section 3.1"},{"comment":"The text says 'line blending ... is not a confounding factor for D/H measurements in comets with PRIMA,' but the 234.64 μm line is one channel away from a bright methanol line and requires careful subtraction. This statement should be qualified, since the methanol contamination is a systematic that must be modeled, not simply absent.","section":"Section 3.2"},{"comment":"The solar-elongation correction is derived from Oort Cloud comets observed in 2000-2020. Please state whether this correction is applied equally to Jupiter-family/Kuiper belt comets, whose orbital and brightness distributions may differ.","section":"Section 4"},{"comment":"The ASTHROS row reports the HDO 588.65 μm line flux in mK km/s, with a 5σ 1h sensitivity of 14.9 mK km/s. For consistency, clarify that this sensitivity is the Herschel/HIFI measured value, not a projection for ASTHROS, since the text later says a future heterodyne instrument would improve on HIFI by a factor of ~2.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"This is a mission-feasibility paper rather than a discovery paper, and the authors are well-positioned given their involvement in Herschel/HIFI comet observations and the SUBLIME coma model. The central claim is plausible but hinges on unvalidated non-LTE line fluxes; the requested uncertainty analysis is well within the manuscript's scope. The self-citations (Refs. 3, 30, 38) are appropriate and do not indicate a circular or inflated argument. The paper fits the journal's scope, and with the sensitivity analysis added, it could become a solid reference for PRIMA's comet-science case."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a well-structured feasibility study for cometary D/H measurements with PRIMA/FIRESS. What is actually new: quantitative non-LTE line flux predictions for the four strongest HDO lines in the PRIMA range, detection-time estimates, a PSG-based blending analysis, line-to-continuum ratios derived from Herschel PACS data, and a corrected accessible-comet count that factors in the telescope Sun-avoidance angle. That last correction (36 down to about a dozen) is an honest adjustment to the earlier Origins Report numbers, and the paper deserves credit for it.\n\nThe arguments for the 184.44 and 234.64 um lines as the best targets are sensible. The blending analysis uses conservative assumptions (highest observed molecular abundances, 40K LTE), and the methanol-subtraction strategy for 234.64 um is plausible. The line-to-continuum ratios, 4-10% at 180 um, are grounded in real PACS observations and look robust.\n\nThe soft spot is the non-LTE flux scale. Table 2 fluxes come from a model with fixed outflow velocity (0.8 km/s) and HDO photolysis rate equal to water's, but no uncertainty is given, and the two preferred lines (upper energies 144 K and 84 K) are not validated against the one well-calibrated cometary HDO line (Herschel 509 GHz, 47 K), which samples a different excitation regime. A factor-2 error in predicted flux changes the headline 6 h detection to ~25 h and shrinks the accessible sample accordingly. That does not kill the mission concept—the brightest comets are still detectable—but it does mean the 'up to a dozen' target count is soft. The authors should add a validation check using existing HDO observations and propagate an uncertainty range through the detection-time and target-count tables. The accessible-comet estimate also depends on launch date and on extrapolating elongation statistics from 2000-2020; the paper acknowledges launch date but doesn't quantify it.\n\nThe citation pattern is appropriate; self-citations are to the relevant Herschel papers and the Origins Report, and the D/H value is an input, not a derived result. The paper is clearly written and the logic is internally consistent.\n\nRecommendation: send it to peer review. The core feasibility claim likely holds, but a referee should push for uncertainty quantification and a model-validation step before the headline numbers are taken at face value.","headline":"A careful PRIMA feasibility study with honest target-count corrections; the detection-time claims rest on unvalidated model fluxes, so referee for revision rather than desk reject.","tokens_in":17721,"tokens_out":2909,"would_cite":false,"duration_ms":31563,"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":"This paper argues that PRIMA's far-infrared spectrometer can measure the D/H ratio in cometary water by detecting HDO lines at 184 and 235 microns, reaching 5-sigma detections in about 6 hours per comet.","keywords":["comets","water D/H ratio","far-infrared spectroscopy","isotopic composition","HDO lines","PRIMA mission","cometary coma excitation","origin of Earth's oceans"],"falsifier":"Measure the HDO 184.44 µm line with FIRESS/FTM toward a comet with water production about 2x10^28 s^-1 at 1 au; if the integrated flux is below 8.23x10^-20 W/m^2 by more than the model uncertainty, the stated 6-hour detection time fails. A high-resolution check on the same comet's HDO 509 GHz line would test the non-LTE excitation model independently.","tokens_in":1982,"feed_emoji":"☄️","tokens_out":5019,"duration_ms":122342,"temperature":0.7,"pith_summary":"This paper tries to establish that the far-infrared space telescope PRIMA can measure the deuterium-to-hydrogen ratio in cometary water for a dozen comets, using lines of the rare isotopologue HDO at 184 and 235 microns. If true, the cometary D/H sample would grow from four accurate space-based measurements to one large enough to compare Oort-cloud and Kuiper-belt reservoirs and to search for correlations with hyperactivity. Such a sample would provide quantitative constraints on where Earth's water came from and whether cometary ice preserves interstellar chemistry or was reprocessed in the protoplanetary disk.","feed_headline":"PRIMA could measure comet water D/H in about 6 hours","feed_subtitle":"That would grow the accurate comet-water sample from 4 objects to a dozen, testing Earth-ocean origins.","key_machinery":"The mechanism is the optically thin HDO line. Since HDO is a rare isotopologue, its far-infrared lines escape the coma and their integrated flux scales with the D/H ratio. A non-LTE excitation code predicts the flux for a given water production rate, outflow velocity, and radiation field; an LTE spectral simulator then checks for contamination by other coma molecules. Together these tools translate a proposed observation into an integration time and an accessible target count.","core_discovery":"The paper claims that PRIMA's FIRESS/FTM can measure the D/H ratio in cometary water by detecting HDO lines at 184.44 and 234.64 µm. For a reference comet with water production 2x10^28 s^-1 at 1 au and D/H twice VSMOW, a non-LTE model gives line fluxes of 8.23x10^-20 and 1.30x10^-19 W/m^2, requiring about 6.5 and 5.7 hours for 5-sigma detections. Blending and line-to-continuum analyses show these lines are usable. With up to a dozen comets accessible during the 5-year mission, PRIMA could expand the accurate comet D/H sample from four to a statistically meaningful set.","pith_inferences":["The same method likely applies to any sufficiently bright icy body, including main-belt comets and Centaurs, even though the paper does not quantify those targets.","The simultaneous oxygen-isotope lines could produce a three-isotope plot for comets analogous to meteorite oxygen-isotope diagrams, a stronger classification tool than D/H alone.","The pre-launch line-flux predictions could be checked against archival far-infrared observations of earlier comets, giving an early test of the excitation model before PRIMA flies.","If on-orbit sensitivity is a factor of two worse than the current best estimate, the six-hour integrations grow to roughly a day and the accessible sample shrinks accordingly."],"forward_implications":["A sample of up to a dozen comets becomes feasible, roughly tripling the number of accurate space-based D/H measurements.","Simultaneous HDO, H2-18O, and H2-17O lines provide multiple isotopic ratios per comet for cross-checks.","For a bright, Garradd-like comet, D/H can be measured from 1 to about 2.7 au, testing whether coma measurements reflect the bulk nucleus or sublimation fractionation.","Comparing Oort-cloud and Kuiper-belt D/H distributions tests whether comets formed in place or were captured from other stars.","Correlating D/H with hyperactivity could identify which physical class of comet delivered Earth's water."],"supporting_citations":[{"why":"Defines the PRIMA mission and the FIRESS/FTM instrument on which the D/H survey is based.","marker":"[7]"},{"why":"Provides the FIRESS/FTM sensitivity curve used for all detection-time calculations.","marker":"[32]"},{"why":"Establishes the HDO and H2-18O line method for D/H measurement and supplies a heterodyne sensitivity baseline.","marker":"[2]"},{"why":"Provides a D/H and oxygen-isotope measurement in a bright Oort-cloud comet used as the heliocentric-distance case.","marker":"[3]"},{"why":"The non-LTE cometary coma model used to compute the HDO line fluxes that drive the detection-time estimates.","marker":"[38]"},{"why":"Sets the HDO photolysis rate assumption used in the excitation model.","marker":"[39]"},{"why":"Generates the LTE coma spectra used to test whether target HDO lines are blended with other molecular lines.","marker":"[41]"},{"why":"Source of the comet feasibility statistics used to estimate the number of accessible targets.","marker":"[23]"},{"why":"Revised isotopic analysis that justifies sampling well-mixed coma gas as representative of the bulk nucleus.","marker":"[29]"},{"why":"Defines the hyperactive-comet science case and provides previous terrestrial D/H values for that class.","marker":"[30]"}],"fun_headline_variants":["PRIMA to triple comet water D/H sample, test ocean origins","Comet water D/H: PRIMA expands sample from 4 to 12","PRIMA's 6-hour comet water D/H detections could triple sample","PRIMA: 6 hours per comet water D/H, 3x more data"],"cache_read_input_tokens":19456,"weakest_assumption_plain":"The whole detection-time and target-count scaling rests on the non-LTE model's predicted HDO line flux; if actual HDO excitation is weaker, the accessible sample shrinks.","fun_headline_variants_meta":{"raw":{"variants":["PRIMA to triple comet water D/H sample, test ocean origins","Comet water D/H: PRIMA expands sample from 4 to 12","PRIMA's 6-hour comet water D/H detections could triple sample","PRIMA: 6 hours per comet water D/H, 3x more data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000627,"raw_usage":{"total_tokens":2702,"prompt_tokens":671,"completion_tokens":2031,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":415,"completion_tokens_details":{"reasoning_tokens":1956}},"tokens_in":415,"tokens_out":2031,"duration_ms":14494,"temperature":1.0,"reasoning_tokens":1956,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T12:08:16.605852+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the HDO 184.44 µm line with FIRESS/FTM toward a comet with water production about 2x10^28 s^-1 at 1 au; if the integrated flux is below 8.23x10^-20 W/m^2 by more than the model uncertainty, the stated 6-hour detection time fails. A high-resolution check on the same comet's HDO 509 GHz line would test the non-LTE excitation model independently.","supporting_citations":[{"cited_title":"Glenn and et al","cited_arxiv_id":null,"evidence_quote":"Defines the PRIMA mission and the FIRESS/FTM instrument on which the D/H survey is based."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the FIRESS/FTM sensitivity curve used for all detection-time calculations."},{"cited_title":"Hartogh , D","cited_arxiv_id":null,"evidence_quote":"Establishes the HDO and H2-18O line method for D/H measurement and supplies a heterodyne sensitivity baseline."},{"cited_title":"Bockel \\'e e-Morvan , N","cited_arxiv_id":null,"evidence_quote":"Provides a D/H and oxygen-isotope measurement in a bright Oort-cloud comet used as the heliocentric-distance case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The non-LTE cometary coma model used to compute the HDO line fluxes that drive the detection-time estimates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets the HDO photolysis rate assumption used in the excitation model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Generates the LTE coma spectra used to test whether target HDO lines are blended with other molecular lines."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Revised isotopic analysis that justifies sampling well-mixed coma gas as representative of the bulk nucleus."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the hyperactive-comet science case and provides previous terrestrial D/H values for that class."}],"review_version":1}