{"id":"1bd9433d-183b-41db-a219-c7d226aa4ce3","arxiv_id":"2508.20215","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Comet 103P/Hartley 2's peak brightness continued fading in 2023 at about 0.59 magnitudes per orbit, implying its active area has shrunk by roughly an order of magnitude since 1991.","lead":"Citizen astronomers' 2023 observations show comet Hartley 2 keeps getting dimmer, losing about 42% of its brightness at every return. The data suggest the well-studied 'hyperactive' comet may be running out of easily vaporized ice.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 2 reduced magnitudes are not reproducible from Eq. (4) and the tabulated mGmin, rh, Δ, so the ΔGmin=0.59 slope is not yet quantitatively supported.","rationale":"The reader's conditional verdict is appropriate. I identify a more internal, checkable weakness than the cross-apparition passband comparison: the headline slope is fit to reduced magnitudes that do not reproduce from the paper's own equations and tables. The four-point trend may survive after recomputation, and prior apparitions independently show fading, so outright rejection is not warranted. But the specific quantitative claims—ΔGmin=0.59±0.11, 42% flux loss per return, active fraction declining by an order of magnitude—are not yet supported as written. The proposed concrete test would settle whether this is a typographical issue in Table 2 or a deeper reduction error.","tokens_in":19811,"tokens_out":21293,"duration_ms":223724,"concrete_test":"Recompute Table 2: for each apparition take the binned upper-envelope photometry and the actual rh, Δ at the quoted peak date, form G = m − 5log10(rhΔ), and refit the four points. If the resulting per-orbit slope differs from 0.59 by more than ~0.15 mag (or if applying Eq. (4) to the published table columns fails to reproduce the published Gmin values), the 42%-per-return and active-fraction claims must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The quantitative central claim—ΔGmin=0.59±0.11 mag/orbit, a 42% flux loss per return—is fit to the four Gmin values in Table 2. Those values do not follow from the paper's own reduction equation. Applying Eq. (4) to the tabulated mGmin and the tabulated rh and Δ gives: 1991: 6.83 − 5log10(0.95×0.719) = 7.66 (table: 7.28); 1997: 7.52 − 5log10(1.03×0.700) = 8.23 (table: 7.84); 2010: 4.47 − 5log10(1.06×0.695) = 5.13 (table: 8.92); 2023: 8.34 − 5log10(1.06×0.693) = 9.01 (table: 10.24). All four entries disagree by ≥0.38 mag, and 2010 disagrees by 3.8 mag. Also, plugging the 2023 fit parameters from Table 1 into Eq. (3) at perihelion (rh=1.06, Δ=0.693) gives m≈7.1, not the tabulated mGmin=8.34. The Gmin column may have been computed with different (peak-date) distances, but then the table's Δ column and the text's claim that the table gives perihelion distances are wrong. Until the reduced magnitudes are recomputed from the original binned photometry and the correct ephemeris, the fitted slope and the derived active-fraction decline are not trustworthy; the qualitative fading may survive, but the headline numbers do not.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a 2023 observing campaign of comet 103P/Hartley 2 by Unistellar/AFA citizen astronomers, presents unfiltered broadband photometry, and fits an upper-envelope lightcurve to the 2023 data. The fitted minimum apparent magnitude is normalized to a reduced magnitude G_min = 10.24 ± 0.47. Combining this value with COBS-based fits for the 1991, 1997, and 2010 apparitions, the authors derive a linear secular fading of ΔG_min = 0.59 ± 0.11 mag per apparition, corresponding to a ~42% flux loss per orbit. From this they estimate that the active fraction has dropped by an order of magnitude since 1991 and that the comet may no longer be hyperactive. The paper interprets the trend as progressive volatile depletion and projects a nuclear-magnitude endpoint around 2112.","tokens_in":20282,"tokens_out":14099,"duration_ms":134460,"significance":"If the quantitative claim is correct, the paper provides a valuable, citizen-science-enabled confirmation of secular fading in a well-studied JFC on human-observable timescales. The public data table and the explicit comparison with the independent narrowband result of Schleicher & Bair (2024) are strengths. However, the central reduction chain has serious internal inconsistencies as printed: the G_min values in Table 2 do not follow from Eq. (4) with the quantities given, and the Table 1 fit parameters do not reproduce the tabulated m_G,min values. Because the headline slope, flux-loss fraction, and active-fraction statements are arithmetic on those tabulated values, the quantitative claims are not currently supported. The qualitative fading trend is plausible and consistent with prior work, but the paper's headline numbers need to be re-derived transparently before publication.","major_comments":[{"comment":"The reduced magnitudes in Table 2 are not reproducible from Eq. (4) and the tabulated perihelion distances. For 2023: 8.34 − 5log10(1.06×0.693) = 9.01, not 10.24. For 2010: 4.47 − 5log10(1.06×0.695) = 5.13, not 8.92. For 1991: 6.83 − 5log10(0.95×0.719) = 7.66, not 7.28. The 2023 and 2010 values are approximately recovered if one uses the geocentric distance at the stated “Peak Date” (e.g., Δ≈0.39 AU on 2023 Oct 2; Δ≈0.12 AU on 2010 Oct 22) rather than the tabulated perihelion Δ. If that is the intended procedure, the table caption and column headings are wrong, and the 1991/1997 entries must be recomputed and shown. Since Fig. 4 and the ΔG_min = 0.59±0.11 fit are based directly on these four G_min values, the slope and all derived quantities (flux ratio f, active-fraction decline, “no longer hyperactive”) are not supported as printed.","section":"§4.1, Table 2, Eq. (4)"},{"comment":"The fitted lightcurve parameters are inconsistent with the dataset they are supposed to describe. For 2023, using Table 1 (HG=9.85, k_pre=17.89) in Eq. (3) at the October 2 peak (rh≈1.07, Δ≈0.39) gives m_G ≈ 9.85 + 0.15 − 7.32 ≈ 2.7 mag, which is ~5.6 mag brighter than the tabulated m_G,min = 8.34 and ~4.9 mag brighter than the brightest Table 3 point near that date (7.60). Similar large discrepancies occur for 1991 and 2010 when the perihelion or perigee geometry is used. This suggests a sign error in Eq. (3), a mismatch between the quantity actually fitted (apparent vs. reduced magnitude), or an error in Tables 1–2. The authors need to supply a worked reduction example or fit diagnostics (residuals, fitted curves over data) so the chain from raw magnitudes to G_min and ΔG_min is checkable.","section":"§3, Eq. (3), Table 1"},{"comment":"The quantitative slope assumes that the unfiltered 2023 magnitudes (described as “most comparable to the Gaia G filter”) are directly comparable to COBS broadband magnitudes from 1991, 1997, and 2010, with no passband or zero-point correction. A systematic offset of only ~0.2–0.3 mag between apparitions would change ΔG_min by 0.05–0.08 mag/orbit and the derived “order of magnitude” active-fraction decline by a factor of roughly 1.5–2. The formal 1σ error of ±0.11 mag does not include such systematics. I request a sensitivity analysis that shifts the 2023 (or each historical) point by a plausible passband/zero-point range and recomputes the slope, or an explicit anchoring of the 2023 photometry to a common photometric system (e.g., comparison with contemporaneous COBS or narrowband magnitudes on the same nights).","section":"§2.1/§4.1"},{"comment":"The active-fraction estimates factive(2017)≈0.7 and factive(2023)≈0.4 are not independent measurements: they are obtained by applying the same fitted slope that is under question, and they also assume factive ∝ flux and constant nucleus size. This is acceptable arithmetic, but should be phrased as a model projection rather than a new observational constraint. The paper's wording (“we can estimate active fractions directly from the measured flux ratios”) is stronger than justified. Please label these as derived from the assumed trend.","section":"§4.2"}],"minor_comments":[{"comment":"The caption says “at perihelion” and gives rh and Δ at perihelion, but the G_min values apparently use peak-date distances. Clarify both in the caption and in the column headers (e.g., “r_h at peak” vs. “Δ at peak”), and state explicitly which distances enter Eq. (4).","section":"Table 2"},{"comment":"The phrase “minimum magnitude in each bin” should be “minimum numerical magnitude (i.e., brightest)” to avoid ambiguity; this is standard but should be stated because the paper is aimed at a broad citizen-science audience.","section":"§3"},{"comment":"In the second bullet, “The minimum reduced magnitude, m_G,min” should be G_min (the symbol m_G,min is used for the apparent minimum magnitude elsewhere). This typo also appears in the Conclusions bullet list.","section":"§5"},{"comment":"Several entries have no reported uncertainty or photometric radius (e.g., many AFA rows). Please indicate whether these were intentionally omitted because errors were not recorded, and give at least a representative description of how the AFA pipeline assigned errors.","section":"Table 3"},{"comment":"The color transformation from V to G is cited to Jordi et al. (2010), but the equation as written appears to have no unit/zero-point explanation. A one-line sanity check would help: with B−V=0.75, V−R=0.43, the offset is −0.0176 −0.005 −0.074 = −0.097, so G≈V−0.10; confirm this is the intended sign and magnitude for the cited transformation.","section":"Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The paper has a valuable dataset and a plausible scientific conclusion, but the quantitative reduction chain is internally inconsistent in ways that are easily demonstrated with the paper's own equations and tables. I do not believe this requires rejection: the authors can likely fix the issue by recomputing the fits, adding a worked example, and correcting Tables 1–2, and then the trend may well stand. However, until the G_min values are reproducible, the headline numbers (ΔG_min=0.59±0.11, 42% flux loss, active fraction decline, no-longer-hyperactive claim) should not be accepted. I recommend major revision with the specific request to provide a reproducible reduction trail from Table 3 to Fig. 4. I found no evidence of fabrication or authorship problems; the main risk is a presentation/reduction-code issue, not scientific misconduct."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new 2023 observations and the citizen-science effort behind them are genuinely useful, and the data are laid out in a transparent table. The qualitative claim that the comet continues to fade is almost certainly right—it was already in the literature through 2010, and your Table 3 does show the comet fainter than in 2010. But the quantitative central claim, ΔGmin = 0.59 ± 0.11 mag/orbit and everything built on it, does not survive contact with the paper's own arithmetic.\n\nThe stress-test check is correct. Applying Eq. (4) to the tabulated mGmin, rh, and Δ gives Gmin values that do not match Table 2: 1991 gives 7.66 instead of 7.28, 2010 gives 5.13 instead of 8.92, and 2023 gives 9.01 instead of 10.24. The 2010 discrepancy is 3.8 mag—not a rounding issue. On top of that, plugging the 2023 fit parameters (HG = 9.85, kpre = 17.89) into Eq. (3) at perihelion gives m ≈ 7.1, not the reported mGmin = 8.34. So the table and the fit disagree with each other. That is not a minor typo; the slope and the derived flux ratio, active-fraction decline, and 'no longer hyperactive' conclusion all rest on those numbers. The paper needs a full recomputation of the reduced magnitudes from the original binned photometry and a re-fit.\n\nThe passband comparability issue (unfiltered 2023 vs. heterogeneous COBS data) is real but secondary. The arithmetic is the load-bearing problem. I also note the 2023 fit parameters look poorly constrained—k values around 15–18 with large errors, and the implied reduced magnitude at the fit's minimum is inconsistent with the table.\n\nWhat is good: the paper gives a complete photometric table, makes a good-faith effort to handle a crowded field, and cites prior work honestly. The qualitative fading trend is consistent with Knight & Schleicher (2013) and with the independent narrow-band result of Schleicher & Bair (2024). That is credit where it is earned.\n\nThis paper should be sent to peer review—the dataset is useful and the question matters—but the referee should be asked to verify the reduction, recompute the Gmin values, and revisit the slope before the quantitative conclusions can be accepted. I would not cite the numbers in their current form.","headline":"The 2023 Hartley 2 photometry is valuable, but the paper's headline numbers do not reproduce from its own equations.","tokens_in":21261,"tokens_out":5294,"would_cite":false,"duration_ms":54365,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Comet Hartley 2 continues a steady fade, losing about 42% of its brightness each six-and-a-half-year return.","keywords":["103P/Hartley 2","comet activity decline","Jupiter-family comet","reduced magnitude","photometry","volatile depletion","citizen science","hyperactivity"],"falsifier":"A decisive check is to re-reduce the 2023 images using only frames whose surface-brightness profile shows no star-contamination spikes and a fixed aperture radius, then compare that Gmin with a same-passband series from the next apparition around 2030; if the 2023 Gmin moves by more than about 0.2 mag under the stricter reduction, or the next apparition does not come in about 0.6 mag fainter, the claimed per-return fading rate is not secure.","tokens_in":19735,"feed_emoji":"☄️","tokens_out":7776,"duration_ms":75706,"temperature":0.7,"pith_summary":"This paper uses photometry from a worldwide citizen-astronomer network to measure Comet 103P/Hartley 2's peak brightness during its 2023 apparition. It finds the comet's minimum reduced magnitude Gmin = 10.24 ± 0.47, and when this is compared with fits to the 1991, 1997, and 2010 apparitions, the peak fades by 0.59 ± 0.11 mag per orbit, meaning roughly 42% of the flux is lost each return. The authors argue this is a real secular decline in activity, not an effect of orbital geometry, and that it implies the comet's active surface fraction has dropped by about an order of magnitude since 1991. They conclude Hartley 2 may no longer qualify as hyperactive and, if the fading continues, could reach bare-nucleus brightness around the year 2112 after at least 13 more active apparitions.","feed_headline":"Hartley 2 fades 42% per return, 2023 data confirm","feed_subtitle":"The once-hyperactive comet has lost ~15x its 1991 peak brightness and may no longer earn the label.","key_machinery":"The load-bearing device is the reduced magnitude G = mG − 5 log10(rhΔ), which removes distance effects, applied to the upper envelope of each apparition's light curve. Brightness is fitted with mG = HG + 5 log10(rh) + k log10(Δ), with separate pre- and post-perihelion slopes; the minimum reduced magnitude from each fit is then regressed against apparition number to obtain the per-return decline. The flux ratio f = 10^(−0.4ΔGmin) and the assumption that active fraction scales linearly with coma flux convert the photometric trend into an activity trend.","core_discovery":"The paper's central claim is that 103P/Hartley 2's activity is still declining on a human-observable timescale. Combining 2023 observations with archived photometry from earlier apparitions, the authors fit each apparition's light curve and track the minimum reduced magnitude Gmin—the brightness the comet would have if placed 1 AU from both the Sun and Earth. They find Gmin rose from 7.28 in 1991 to 10.24 in 2023, a linear trend of 0.59 ± 0.11 mag per 6.48-year orbit. Converting this to flux, each return delivers only 58% ± 6% of the previous apparition's peak outgassing, so peak activity is now lower than in 1991 by a factor of about 15, with a 1σ range of 9–26. Because coma brightness is t","pith_inferences":["Beyond the paper's claims: if the trend is real, Hartley 2 could become a live example of a comet crossing from hyperactive to ordinary and eventually dormant, giving observers a rare chance to watch the end stages of a Jupiter-family comet over just a few decades.","A cleaner test than fitting heterogeneous archived magnitudes would be a single-passband, aperture-matched photometric program across the next apparition around 2030, plus contemporaneous narrowband water-production measurements; the prediction is another ~0.6 mag fainter Gmin and a commensurate drop in gas production.","The upper-envelope approach assumes the brightest points are uncontaminated; if any of those points are background-star blends, Gmin would be biased. Checking the 2023 frames for star-free windows against the surface-brightness spike filter would tighten the result.","The active-fraction numbers assume constant nucleus size and a linear flux–production scaling over five orbits; both are approximations, so the 'no longer hyperactive' conclusion should be read as a quantitative estimate with systematic uncertainty comparable to the reported statistical error."],"forward_implications":["Hartley 2's peak brightness in 2023 was about 15 times lower than in 1991, with a 1σ range of roughly 9–26 times lower.","The comet's effective active surface fraction dropped from about 1.17 at the 2010 perihelion to roughly 0.4 in 2023, so it no longer meets the definition of a hyperactive comet.","If the linear trend continues, Hartley 2 will reach bare-nucleus magnitude around 2112 ± 18, with about 13 more active orbits as a lower limit; activity more likely fades asymptotically.","The broadband decline of 42% ± 6% per orbit agrees with the roughly 40% per-orbit drop in water production reported previously, supporting progressive volatile depletion rather than observing geometry as the cause.","Because the unfiltered flux includes a gas contribution of roughly 15–25%, dust mass-loss estimates from this photometry are upper limits."],"supporting_citations":[{"why":"Established the prior per-apparition decline in water production (~40%) that the new broadband decline is compared against.","marker":"Knight & Schleicher (2013)"},{"why":"Supplies the archived photometry used to fit the 1991, 1997, and 2010 apparition light curves.","marker":"COBS (2024)"},{"why":"Provides the 2010 active-fraction value (~1.17) used as the anchor for estimating later active fractions.","marker":"Lisse et al. (2009)"},{"why":"Independent narrowband photometry that also found the secular decrease continuing through 2023.","marker":"Schleicher & Bair (2024)"},{"why":"Documents the icy-chunk hyperactivity mechanism that the declining-active-fraction conclusion modifies.","marker":"A'Hearn et al. (2011)"},{"why":"Supplies the color transformation used to convert the nucleus V-band magnitude to the Gaia G band for the 2112 projection.","marker":"Jordi et al. (2010)"},{"why":"Provides the nucleus absolute magnitude and colors used in the nucleus-brightness projection.","marker":"Li et al. (2013)"}],"fun_headline_variants":["Hartley 2 fades 42% each return, may lose 'hyperactive' status","Comet's peak brightness drops 42% per orbit, 2023 data show","Once-hyperactive comet now 15x dimmer than in 1991","Brightness of Hartley 2 falls 42% per apparition, study finds","Comet 103P: activity decline points to volatile exhaustion"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The trend assumes that unfiltered 2023 magnitudes, whose effective passband is closest to Gaia G, are directly comparable to the heterogeneous archived magnitudes from 1991, 1997, and 2010; a systematic passband or zero-point offset between apparitions would change the fitted ΔGmin and the active-fraction conclusion.","fun_headline_variants_meta":{"raw":{"variants":["Hartley 2 fades 42% each return, may lose 'hyperactive' status","Comet's peak brightness drops 42% per orbit, 2023 data show","Once-hyperactive comet now 15x dimmer than in 1991","Brightness of Hartley 2 falls 42% per apparition, study finds","Comet 103P: activity decline points to volatile exhaustion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000256,"raw_usage":{"total_tokens":1444,"prompt_tokens":806,"completion_tokens":638,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":531}},"tokens_in":550,"tokens_out":638,"duration_ms":6190,"temperature":1.0,"reasoning_tokens":531,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:13:43.782048+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is to re-reduce the 2023 images using only frames whose surface-brightness profile shows no star-contamination spikes and a fixed aperture radius, then compare that Gmin with a same-passband series from the next apparition around 2030; if the 2023 Gmin moves by more than about 0.2 mag under the stricter reduction, or the next apparition does not come in about 0.6 mag fainter, the claimed per-return fading rate is not secure.","supporting_citations":[{"cited_title":"M., Fernandez, Y","cited_arxiv_id":null,"evidence_quote":"Provides the 2010 active-fraction value (~1.17) used as the anchor for estimating later active fractions."},{"cited_title":"2024, 56, 314.06","cited_arxiv_id":null,"evidence_quote":"Independent narrowband photometry that also found the secular decrease continuing through 2023."}],"review_version":1}