{"id":"f2e73faf-669a-418c-834d-b3c85abc3c1e","arxiv_id":"2505.03236","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Compact symmetric objects show low-amplitude optical variability and a bluer-when-brighter trend, with weaker amplitude than blazars, consistent with a jet-origin interpretation.","lead":"Researchers used five years of Zwicky Transient Facility data to measure optical brightness changes in 38 compact symmetric objects, young galaxies with small radio jets. They found these objects vary less than blazars, whose jets point almost at us, and that both types show a trend of becoming bluer when brighter.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Host-galaxy light is not subtracted or bounded; the 2–3× CSO/blazar Fvar gap and the BWB trend used to infer beaming can be produced by constant red starlight diluting the nuclear signal.","rationale":"I read the paper as a first systematic optical variability census of a bona fide CSO sample, with a secondary comparative claim that the CSO/blazar Fvar contrast and common BWB behaviour are explained by viewing angle and relativistic beaming. The census itself is supported by public ZTF data, a well-defined source sample from Kiehlmann et al. (2024a), and standard Fvar and chi-square machinery; I do not see an internal inconsistency in the descriptive statistics. The load-bearing weak point is precisely the step from observed Fvar and colour slopes to the beaming interpretation. That step requires the assumption that the CSO light curves are dominated by the active nucleus, but the paper's own sample is 31/38 galaxies and the matched sub-sample includes 6 galaxies. A constant red host dilutes the observed fractional variability and simultaneously creates a blue-when-brighter colour trend, so both observable signatures used to argue for weak beaming are degenerate with host contamination. The paper does not quantify, model, or subtract host light, and its Section 4 attributes the Fvar difference entirely to Doppler boosting. This is a missing control, not a matter of disagreement with the community consensus on CSO orientation. The concrete test is straightforward and decisive: if host-corrected CSO variability remains suppressed and the BWB slopes persist, the physical interpretation is supported; if not, the central comparison is an artifact of photometric aperture contamination. Since the reader's conditional verdict already targets this same assumption, I see no need to move the verdict; it should remain conditional pending that test.","tokens_in":22347,"tokens_out":3858,"duration_ms":42560,"concrete_test":"For the 9 matched CSOs and ideally all 38, fit PSF + Sérsic host models to SDSS or PanSTARRS images (or use ZTF difference imaging or a published nuclear-to-total flux ratio) to estimate the host flux in g and r. Use the best-fit host flux to correct each epoch's ZTF flux, then recompute Fvar (Eq. 3) and the g−r versus g slopes (Table 5). If host-corrected CSO Fvar remains ~2–3× below the blazar mean and most CSOs still show a significant negative slope, the beaming interpretation survives; if the gap shrinks toward unity or the BWB slopes become flat or positive, host dilution is the dominant explanation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central physical claim (Section 4 and Conclusions item 4) is that the lower CSO variability is due to weak Doppler beaming. The evidence is the Fvar ratio and the BWB colour trend, but both are computed on total ZTF flux. For a constant host galaxy of flux H and a variable nucleus of flux N(t), Eq. (3) yields Fvar_obs ≈ Fvar_nuc × N/(N+H); the colour-magnitude relation also gains a blue-when-brighter component whenever a red, constant host is added to a variable blue nucleus, even if the nucleus itself has no wavelength-dependent colour change. 31 of 38 CSOs are classified as galaxies, and the matched sub-sample still contains six galaxies, while the blazar comparison objects are point-like. The paper nowhere measures or subtracts the host contribution, so the entire 2–3× amplitude gap and the apparent BWB behaviour can be explained by dilution without any difference in jet beaming. The descriptive census—that CSOs show low-amplitude optical variability—is likely robust; the attribution-specific interpretation that beaming sets the CSO/blazar contrast is not isolated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a ZTF-based study of long-term optical variability of 38 compact symmetric objects (CSOs) from the Kiehlmann et al. (2024a) bona fide catalog. The authors compute chi-square variability classifications and fractional variability amplitudes (Fvar) in g, r, and i bands, and compare a redshift- and magnitude-matched sub-sample of 9 CSOs with 5 FSRQs and 12 BL Lacs. They report low-amplitude variability in CSOs (mean Fvar roughly 0.084, 0.070, and 0.066 in g, r, and i), a 2-3x larger Fvar in blazars, and a bluer-when-brighter (BWB) colour trend in both classes. They interpret these results as evidence that CSO optical emission is jet-dominated and that the smaller amplitude in CSOs reflects weaker Doppler beaming due to larger viewing angles.","tokens_in":22587,"tokens_out":8448,"duration_ms":77393,"significance":"The descriptive census is useful: this appears to be the first systematic multi-band optical variability characterization of a bona fide CSO sample, and it uses standard, reproducible recipes on public ZTF data. The detection of variability in roughly half of the CSOs in all three bands and the BWB trend are interesting empirical results that will be a useful reference for future CSO work. However, the beaming interpretation is not uniquely supported because host-galaxy starlight can mimic both the amplitude suppression and the colour trend. If the authors can bound or remove the host contribution, the comparison would become a valuable constraint on orientation-based unification. As it stands, the central comparative claim is weakened.","major_comments":[{"comment":"The attribution of the lower CSO variability to weaker Doppler beaming is not isolated from host-galaxy dilution. In Section 2.1, 31 of 38 CSOs are classified as galaxies, and the matched sub-sample used for Table 4 includes six galaxies, whereas the blazar comparison objects are point-like. For a constant host flux H and a variable nuclear flux N(t), the fractional variability amplitude defined in Eq. (3) satisfies Fvar_obs = Fvar_nuc * <N>/(<N>+H); therefore a 2-3x suppression can be produced by host starlight alone, with no difference in beaming. Similarly, adding a red constant host to a variable blue nucleus produces a BWB colour-magnitude correlation even if the nuclear colour is constant. The manuscript does not measure or bound H/N, so the beaming interpretation stated in Section 4 and Conclusions items 4 and 5 is one possibility but not a demonstrated one. I ask for a host-subtraction or host-fraction analysis (for example, surface-brightness decomposition of the ZTF/SDSS images, or an upper limit on H/N derived from quasar-dominated CSOs) and a re-comparison on host-corrected nuclear fluxes.","section":"Section 4 / Conclusions item 4"},{"comment":"The claim that blazars are two to three times more variable than CSOs is not supported by a statistical test. The comparison in Table 4 involves only 7 CSOs, 5 FSRQs, and 10 BL Lacs, yet the paper reports only mean values and errors. Please add a two-sample test per band (e.g., Kolmogorov-Smirnov or Mann-Whitney) on the Fvar distributions, and also test whether the ordering Fvar_g > Fvar_r > Fvar_i is significant within each class rather than inferred from the means.","section":"Section 3.1.2 / Table 4"},{"comment":"The colour-magnitude slopes reported in Table 5 are not consistent with the 'similar BWB behaviour' narrative. The CSO slopes are typically 0.8-1.5 mag per mag, while the blazar slopes are mostly below 0.5 mag per mag. A steep CSO slope is a natural signature of dilution by a red constant host, so this quantitative difference needs to be fitted with the host model and not simply grouped with blazars under the same jet interpretation. Please report the host-corrected slopes or explain the difference explicitly.","section":"Section 3.2 / Table 5"}],"minor_comments":[{"comment":"The matching criteria are stated as r-band magnitude within ±0.5 mag, the abstract says similar g-band magnitudes, and Table 1 lists only g; please clarify which quantity was matched and list the r magnitudes used.","section":"Section 2.2"},{"comment":"The panel labelled J1158+1022 corresponds to the source J1158+2450 in Table 1; fix the label.","section":"Figure 5"},{"comment":"After requiring variability in all three bands, the matched sample in Table 4 has 7 CSOs and 10 BL Lacs, while Section 2.2 introduces 9 CSOs and 12 BL Lacs; state the reduction explicitly.","section":"Section 3.1.2 / Table 4"},{"comment":"The quasar/galaxy classification via absolute B magnitude uses total SDSS magnitudes and therefore includes host starlight; this should be stated as a caveat when interpreting Table 3.","section":"Section 2.1"}],"recommendation":"major_revision","confidential_remarks":"The host-galaxy dilution issue is the main obstacle; if the authors can quantify H/N and show that the Fvar gap persists for quasar-dominated CSOs, I would support acceptance. The statistical tests and the colour-slope difference in Table 5 should also be addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper delivers something the field didn't have—a systematic, bona fide-sample measurement of long-term optical variability in CSOs. The ZTF-based Fvar values and detection fractions are useful and will be cited. But the paper's central comparative claim—that the 2–3x lower variability in CSOs relative to blazars is due to weaker Doppler beaming—is not supported as it stands, because host-galaxy light was never subtracted or bounded. That is a real gap, not a minor quibble.\n\nWhat is good: first systematic optical variability census of 38 bona fide CSOs with 5-year ZTF light curves; standard, reproducible methods (chi-square, Fvar, colour-magnitude fits); clear tabulations; appropriate comparison blazar sub-sample matched in redshift and brightness. The descriptive result—CSOs do vary at low amplitude (mean Fvar 0.08, 0.07, 0.06 in g, r, i)—is solid and new.\n\nWhere it gets shaky: Section 4 and Conclusions items 4–5 attribute the CSO/blazar Fvar gap and the BWB trend to beaming. But 31 of 38 CSOs are classified as galaxies, and the matched sub-sample still includes six galaxies. Blazars are point-like; CSO photometry includes unresolved host starlight. For a constant red host plus variable blue nucleus, the observed Fvar is suppressed by N/(N+H), and the colour-magnitude relation acquires a BWB component. So both the amplitude gap and the BWB signal can arise without any difference in intrinsic beaming. The paper even reports CSO quasars (brighter, less host-dominated) varying more than galaxies—consistent with the dilution story. The authors never measure or model H. Until they do, the beaming interpretation is not isolated. The low-amplitude variability result stands; the explanation does not.\n\nMinor: the matched sample is small (9 CSOs, 5 FSRQs, 12 BL Lacs), and the colour-magnitude slopes are not formally compared between populations.\n\nBottom line: worthwhile, citable for the census; the physical interpretation needs re-analysis with host-galaxy correction or at least a quantitative dilution bound. Send to a good referee.","headline":"Useful first census of CSO optical variability; the beaming interpretation is not yet isolated from host-galaxy dilution.","tokens_in":23101,"tokens_out":2850,"would_cite":true,"duration_ms":28654,"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":"The paper claims that compact symmetric objects show low-amplitude optical variability—about two to three times weaker than blazars—and that the difference traces to weaker relativistic beaming, not to a different emission mechanism.","keywords":["compact symmetric objects","AGN optical variability","relativistic jets","blazars","bluer-when-brighter","fractional variability amplitude","Doppler beaming","Zwicky Transient Facility"],"falsifier":"Measure the host-galaxy light in the ZTF images of the nine matched CSOs by fitting and subtracting a galaxy profile, then recompute $F_{\\mathrm{var}}$ from the nuclear flux alone. If the host-subtracted amplitudes rise to blazar levels, the variability gap is an artifact of steady starlight rather than weaker beaming; if the gap persists, the viewing-angle interpretation is supported.","tokens_in":22202,"feed_emoji":"🌌","tokens_out":12719,"duration_ms":112474,"temperature":0.7,"pith_summary":"This paper tries to establish that compact symmetric objects (CSOs), the young, small radio galaxies whose jets are thought to lie at large angles to the line of sight, are low-amplitude optical variables whose light comes from the relativistic jet rather than from the host galaxy alone. Using five years of Zwicky Transient Facility data in g, r, and i bands for 38 bona fide CSOs and a matched sub-sample of 17 blazars, it reports mean fractional variability amplitudes of $F_{\\mathrm{var}} \\approx 0.084$, $0.070$, and $0.066$ in g, r, and i—roughly two to three times below the blazar values. It also reports that both populations turn bluer as they brighten. The paper concludes that the shared bluer-when-brighter behaviour points to the same jet-dominated emission mechanism, while the smaller CSO amplitude follows because their jets are seen at large angles and are only weakly Doppler-boosted. If correct, this makes CSOs a useful low-beaming control population for jet-physics studies.","feed_headline":"Compact symmetric objects flicker far less than blazars","feed_subtitle":"Five years of ZTF data show their optical variability is 2–3 times weaker yet follows the same bluer-when-brighter trend.","key_machinery":"The main measurement is the fractional variability amplitude $F_{\\mathrm{var}} = \\sqrt{(S^2 - \\bar{f}_{\\mathrm{err}}^2)/\\bar{f}^2}$, which compares the excess variance of a light curve to its mean flux and serves as the paper's yardstick for how variable a source is. The second piece is the matched-sample comparison, in which CSOs, FSRQs, and BL Lacs are paired by redshift and g-band brightness so that any variability difference is attributed to orientation rather than distance or luminosity. The third is the colour-magnitude diagram of near-simultaneous g-r colour against g magnitude, fit with a weighted least-squares line, whose slope reveals the bluer-when-brighter behaviour. The physical mechanism invoked to explain the amplitude gap is relativistic Doppler beaming tied to viewing angle.","core_discovery":"The paper's central discovery is that a bona fide sample of 38 compact symmetric objects varies optically on month-to-year timescales with low amplitude: roughly 76% of those with usable g data, 87% of those with r data, and 78% of those with i data are variable, with 21 sources variable in all three bands. For the nine CSOs matched by redshift and g-band brightness to 5 flat-spectrum radio quasars and 12 BL Lac objects, the same analysis shows the CSO variability to be about two to three times weaker, with mean $F_{\\mathrm{var}}$ values of $0.081\\pm0.006$, $0.078\\pm0.001$, and $0.065\\pm0.002$, against $0.236\\pm0.003$, $0.188\\pm0.001$, and $0.136\\pm0.004$ for FSRQs and $0.251\\pm0.001$, $0.228\\pm0.001$, and $0.185\\pm0.001$ for BL Lacs. In addition, both CSOs and blazars show a bluer-when-brighter colour-magnitude trend, with 20 of the 25 CSOs that have detectable colour variations showing it. The paper interprets these results as evidence that the optical emission of CSOs is jet-produced and that the variability gap is a beaming effect: blazar jets point nearly at us, so perturbations are Doppler-amplified, while CSO jets lie at large viewing angles and receive only modest amplification.","pith_inferences":["If host-galaxy starlight substantially dilutes CSO variability, the beaming interpretation would need revision; a clean test is to recompute $F_{\\mathrm{var}}$ from host-subtracted images of the matched CSOs.","Because CSOs show bluer-when-brighter colour changes at amplitudes near the noise floor, they offer a low-Doppler-boost regime in which particle-injection and cooling models can be compared without the extreme amplification seen in blazars.","Extending the matching scheme to rest-frame bands and k-corrected luminosities could sharpen the colour-slope comparison and reveal whether the bluer-when-brighter slope itself depends on viewing angle.","A larger CSO sample could test a quantitative prediction: if beaming drives variability, the fractional variability amplitude should scale with radio core prominence, a proxy for viewing angle, across the population."],"forward_implications":["CSOs join the class of jet-dominated, optically variable AGN rather than being quiet young radio sources in the optical band.","The shared bluer-when-brighter trend in CSOs and blazars implies that the same basic jet process, likely particle acceleration in shocks, operates across a wide range of jet orientations and Doppler factors.","The tentatively higher variability of CSO quasars over CSO galaxies in the r band supports orientation-based unification and motivates a larger sample to test it.","The measured $F_{\\mathrm{var}}$ ratios between blazars and CSOs provide a route to estimate relative Doppler-boosting factors once host-galaxy light is accounted for.","Longer or higher-cadence monitoring of CSOs can connect their optical flickering to their known gamma-ray detections and test whether the jets remain active at small scales."],"supporting_citations":[{"why":"Defines the bona fide CSO catalogue that supplies the 38 sources and their optical and radio classifications.","marker":"Kiehlmann et al (2024a)"},{"why":"Describes the ZTF survey whose public database supplied the multi-band optical light curves.","marker":"Graham et al 2019"},{"why":"Documents the ZTF camera and observing system that produced the g, r, and i time series.","marker":"Bellm et al, 2019"},{"why":"Provides the fractional variability amplitude estimator and its uncertainty formula used as the main variability metric.","marker":"Vaughan et al, 2003"},{"why":"Supports the same $F_{\\mathrm{var}}$ definition and the practice of classifying sources as variable, probable, or non-variable.","marker":"Sesar et al, 2007"},{"why":"Supplies the ROMA-BZCAT catalogue from which the FSRQ and BL Lac comparison objects were selected.","marker":"Massaro et al, 2015"},{"why":"Gives the Galactic extinction corrections applied to the magnitudes before computing fluxes and colours.","marker":"Schlafly and Finkbeiner (2011)"},{"why":"Provides the shock-in-jet model invoked to explain jet-driven optical flux variations.","marker":"Marscher and Gear, 1985a"},{"why":"Explains how changes in the Doppler factor can produce a bluer-when-brighter colour trend.","marker":"Villata et al, 2004"}],"fun_headline_variants":["CSOs flicker less than blazars, but share the same color trend","Why CSOs are optically steadier: beaming, not absence of jets","Optical twinkle of CSOs is 2-3 times weaker than blazars","Same bluer-when-brighter behavior, but CSOs are less variable","Jets at an angle: CSOs show muted optical variability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion that the smaller CSO variability is caused by weaker Doppler beaming assumes that steady host-galaxy starlight does not dilute the measured variability, even though 31 of the 38 CSOs are galaxies while the comparison blazars are point-like.","fun_headline_variants_meta":{"raw":{"variants":["CSOs flicker less than blazars, but share the same color trend","Why CSOs are optically steadier: beaming, not absence of jets","Optical twinkle of CSOs is 2-3 times weaker than blazars","Same bluer-when-brighter behavior, but CSOs are less variable","Jets at an angle: CSOs show muted optical variability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000538,"raw_usage":{"total_tokens":2667,"prompt_tokens":1114,"completion_tokens":1553,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":730,"completion_tokens_details":{"reasoning_tokens":1450}},"tokens_in":730,"tokens_out":1553,"duration_ms":16032,"temperature":1.0,"reasoning_tokens":1450,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:56:31.416414+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the host-galaxy light in the ZTF images of the nine matched CSOs by fitting and subtracting a galaxy profile, then recompute $F_{\\mathrm{var}}$ from the nuclear flux alone. If the host-subtracted amplitudes rise to blazar levels, the variability gap is an artifact of steady starlight rather than weaker beaming; if the gap persists, the viewing-angle interpretation is supported.","supporting_citations":[{"cited_title":"Exploring the Variable Sky with the Sloan Digital Sky Survey","cited_arxiv_id":"0704.0655","evidence_quote":"Supports the same $F_{\\mathrm{var}}$ definition and the practice of classifying sources as variable, probable, or non-variable."}],"review_version":1}