{"id":"9c770b32-7f4f-4daa-afff-d48a6f3817cf","arxiv_id":"2506.07589","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"At parsec scales, FRII radio galaxies with efficient and inefficient accretion have comparable radio luminosity, and the optical-radio classification separation is already visible.","lead":"Astronomers measured the radio brightness of 3CR radio galaxies with very long baseline interferometry, reaching scales of a few parsecs. They found that the usual separation between optical-radio classes already exists at the base of the jet, which favors the idea that the central engine, not the environment, sets the radio class.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The equality claim rests on the VLBI subsample being class-unbiased, yet the Fundamental Radio Catalog's flux-dependent selection is never modeled, and the paper's parent-sample checks do not test selection on parsec-scale core flux.","rationale":"I read the paper in good faith: it is a useful, clearly presented comparison of public VLBI data, and the parent-sample histograms and nonparametric tests are reasonable first-order checks. However, the strongest claim—that FRII RGs with different accretion rates already have comparable radio luminosity on parsec scales—requires that the VLBI-detected subsample be representative of the parent 3CR classes in exactly the quantity being compared, L8GHz at parsec scales. The Fundamental Radio Catalog has a selection function dominated by VLBI detectability and astrometric usability, which depends on compact correlated flux density. That is not the same quantity as L178MHz or L[OIII], so the paper's checks cannot exclude class-dependent truncation. In particular, if FRII-LERGs with low parsec-scale core luminosity are less likely to enter the catalog, their observed L8GHz distribution will be biased upward and could mimic equality with FRII-HERGs. The small class subsamples and the requirement that both MW and KS p-values fall below 0.05 only aggravate the problem: the tests reject a null only when the evidence is quite strong, so a p-value above threshold is not evidence of equivalence. This is why the conditional verdict already given by the reader is appropriate: the engine-scenario interpretation should not be accepted as established until the selection function is modeled or a direct check of the missing sources is performed. Credit should go to the authors for using a public catalog, showing the parent-overlap histograms, and including compactness information; the concern is not that they acted improperly, but that a specific unsupported assumption is load-bearing. The proposed logistic-regression plus inverse-probability-weighting check, or a targeted VLBA snapshot of the 24 non-detected sources, would settle whether the null result is physical or an artifact. Until then, the central assertion is best treated as a tantalizing but not fully demonstrated result.","tokens_in":11268,"tokens_out":8259,"duration_ms":103974,"concrete_test":"Use the B10 parent sample of 79 sources: assemble for every source an arcsecond-scale core flux from VLA images, the literature, or archive cutouts, and fit a logistic model for RFC detection as a function of class, redshift, L178MHz, L[OIII], and core flux. Then recompute the FRII-LERG versus FRII-HERG L8GHz comparison with inverse-probability-of-detection weights, or equivalently treat the observed fluxes as truncated at the RFC detection threshold. If the weighted median offset between the two classes moves by more than 0.3 dex, or if either p-value crosses 0.05, selection bias is the likely driver of the claimed equality; if the weighted result is unchanged and detection probability is class-independent, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 builds the 3CR-VLBI sample by cross-matching 79 B10 3CR RGs with the Fundamental Radio Catalog (Petrov & Kovalev 2025), obtaining 55 associations. Because the catalog is a heterogeneous reanalysis of 1,088 VLBI experiments and not a complete survey, a source appears only if it was targeted, observed, and detected in at least one experiment. The central claim of Section 3—that FRII-LERGs and FRII-HERGs have comparable 8 GHz parsec-scale luminosity—therefore presupposes that RFC detection is independent of radio-optical class at fixed redshift and total radio power. The paper's checks do not establish this. Figure 1 compares L178MHz and L[OIII] histograms between the parent and VLBI subsamples, but L178MHz traces kiloparsec-scale lobes and L[OIII] traces accretion, not the parsec-scale core flux that governs RFC inclusion; Figure 2's compactness C8GHz is conditional on detection, so sources with cores too faint to detect never enter the sample. With 55 objects split into class subsamples of roughly 12-30 sources, the nonparametric tests have limited power to detect a selection difference, and no uncertainty or detection-completeness model is reported. If FRII-LERGs with weaker parsec-scale cores are preferentially omitted, their observed L8GHz distribution shifts upward and the reported pMW=0.3/pKS=0.4 null result could be a selection artifact rather than evidence for the engine scenario. The additional α=0 extrapolation for eight sources is a smaller issue, but it further underscores that the luminosity values carry unquantified systematic uncertainty.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses VLBI data from the Fundamental Radio Catalog (Petrov & Kovalev 2025) to measure parsec-scale 8 GHz radio luminosities of 55 3CR radio galaxies with optical and radio classifications from Buttiglione et al. (2009, 2010, 2011). The authors compare the luminosity distributions of FRI-LERGs, FRII-LERGs, and FRII-HERGs using Mann-Whitney and Kolmogorov-Smirnov tests. They find no statistically significant difference between FRII-LERGs and FRII-HERGs (p_MW = 0.3, p_KS = 0.4) and interpret this as evidence for the 'engine scenario,' in which the accretion regime and central engine properties determine parsec-scale jet power, against the evolutionary and environmental scenarios.","tokens_in":11569,"tokens_out":3696,"duration_ms":43493,"significance":"If the central claim is correct—that FRII-LERGs and FRII-HERGs have comparable parsec-scale radio luminosities—the result would place the differentiation of radio-optical classes at the jet base, before kiloparsec-scale environmental interactions, and would strengthen the engine scenario for jet launching. The paper's strengths include the use of a well-defined parent sample (79 sources with optical classifications), an independent VLBI catalog, and a first-order compactness check for orientation and resolution biases. However, the main inference rests on a null result obtained with a small sample and a selection function that is not fully modeled, so the conclusion is currently not as robust as the text suggests.","major_comments":[{"comment":"The claim that 'FRII RGs with different accretion rates already have comparable radio luminosity on parsec scales' (final sentence of Section 3) is not established by the reported p-values. Failing to reject the null hypothesis with subsamples of roughly 12-30 objects does not demonstrate that the distributions are comparable; it only indicates that the test lacked sensitivity to the offset present in the data. The paper should include a power analysis or a simulation-based estimate of the minimum median luminosity offset that could be detected at the adopted 5% threshold, and should state the physical offset that would be ruled out. Without this, the null result is uninformative and the wording 'comparable' overreaches.","section":"Section 3"},{"comment":"The Fundamental Radio Catalog is a heterogeneous compilation of 1,088 VLBI experiments, not a complete survey, so inclusion in the 3CR-VLBI subsample depends on whether a source was targeted, detected, and resolved at parsec scales. The parent-sample checks shown in Figure 1 compare L178MHz and L[OIII], which are proxies for kpc-scale lobe power and accretion rate, respectively, but not for the parsec-scale core flux that governs RFC selection. The compactness distributions in Figure 2 (left panel) are conditional on detection and cannot reveal the properties of sources whose cores were too faint or too resolved to enter the catalog. If FRII-LERGs with weaker parsec-scale cores are preferentially missing, the observed similarity at 8 GHz could be a selection artifact. The authors should either model the selection function (e.g., using survival analysis with upper limits from other VLBI surveys or comparing with a complete VLBI sample of 3CR sources) or quantitatively assess the magnitude of a plausible selection bias that would erase a real luminosity offset.","section":"Section 2"},{"comment":"The sentence 'We note that no redshift-related biases are present' is asserted without any supporting test or figure. Given that 3CR is a flux-limited sample and the parent and VLBI subsamples have different sizes, the redshift distributions of FRII-LERGs and FRII-HERGs should be compared directly (e.g., with a two-sample KS test) to justify this claim, especially because luminosity at a fixed frequency is redshift-dependent through distance and K-correction.","section":"Section 2"}],"minor_comments":[{"comment":"For eight sources the 8 GHz luminosity is extrapolated from the nearest band assuming a flat spectral index (alpha = 0). The uncertainty introduced by this assumption is not propagated into the statistical tests; a robustness check excluding these sources (or varying alpha within plausible bounds, e.g., -0.5 to 0.5) would strengthen the results.","section":"Appendix A, Table A.1"},{"comment":"Unresolved flux densities are available for only 44 of 55 objects, so compactness C8GHz is undefined for 20% of the sample. The impact of excluding these objects from the compactness comparisons should be mentioned explicitly.","section":"Section 2"},{"comment":"The histograms in the right panel use a fixed bin size of 0.5 dex; with subsample sizes of 12-30, the choice of binning can affect the visual impression. Showing individual luminosity values (e.g., as a strip chart or bean plot) would provide a more transparent view of the distributions.","section":"Figure 2"},{"comment":"The abstract states that 'the separation between FRII classes with different accretion rates occurs near the central engine,' but the body finds no statistically significant separation in luminosity between FRII-LERGs and FRII-HERGs. This phrasing is ambiguous and could be misread as claiming a separation where none is found; it should be clarified that the similarity itself is the finding.","section":"Abstract"},{"comment":"The reference to Edge et al. (1959) as 'Geochim. Cosmochim. Acta' appears to be incorrect; the Third Cambridge Catalogue of Radio Sources was published in Monthly Notices of the Royal Astronomical Society. Please verify and correct this citation.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is concise and addresses a timely question in jet/accretion physics. The main concern is that the central inference—comparable parsec-scale luminosities—is built on a null result without a power analysis or a model of the VLBI selection function. I believe the paper is potentially publishable after the authors quantify the sensitivity of their tests and address the selection bias in a more quantitative way. I also note the possibly incorrect Edge et al. (1959) reference, which may indicate a broader citation-verification issue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this paper does something genuinely useful — it assembles the 3CR radio galaxies with optical-radio classifications and checks whether their parsec-scale 8 GHz luminosities separate by class, using the Fundamental Radio Catalog. That specific comparison is new, and it is a sensible, low-cost use of public VLBI data. The sample construction is careful: cross-matching with a 5-arcminute radius, visual inspection of VLA cutouts, and the histograms in Fig. 1 show the VLBI subset broadly tracks the parent sample in L[OIII] and L178MHz. The nonparametric tests are appropriate for the sample sizes.\n\nThe soft spots are real. The central claim that FRII-LERGs and FRII-HERGs have comparable parsec-scale luminosity rests on failing to reject the null with 55 objects split into small subsamples. There is no power analysis, no equivalence test, no quantification of what luminosity offset the data could actually rule out. The abstract even says \"separation between FRII classes with different accretion rates occurs near the central engine,\" which reads like the opposite of what the tests show; the FRII-LERG and FRII-HERG distributions overlap. That wording needs to be fixed.\n\nMore importantly, the selection function of the Fundamental Radio Catalog is never modeled. The catalog is a heterogeneous collection of VLBI experiments; a source enters only if it was targeted and detected. The parent-sample checks compare L178MHz and L[OIII], but those do not control for parsec-scale core flux, which is exactly what governs inclusion in a VLBI catalog. The compactness distribution in Fig. 2 is conditional on detection, so it cannot reveal faint-core sources that were missed. If the missing FRII-LERGs tend to have weaker parsec-scale cores, the observed \"equality\" is a selection artifact. The alpha=0 extrapolation for eight sources is a minor extra worry, but it reinforces that the luminosities carry unquantified systematic uncertainty.\n\nThe theoretical framing is reasonable: the engine scenario is a legitimate alternative, and the BZ/BP discussion is fine. But the data as presented do not distinguish between scenarios as strongly as claimed. A null result with an unmodeled selection function cannot carry that weight.\n\nWho is this for? People working on AGN jet-accretion connections and radio-galaxy populations will want to know about the comparison, and the paper deserves a serious referee — it is a legitimate observational test with a clear question. But it needs substantial revision: add a power analysis or equivalence bounds, model or at least discuss the RFC selection function, report luminosity uncertainties, and tighten the abstract. I would send it to review rather than desk reject, but I would expect major changes before it is solid.\n\nRecommendation: engage with it, but treat the central conclusion as provisional.","headline":"A genuinely useful archival VLBI comparison whose central null result is oversold; needs a power analysis and a harder look at catalog selection before favoring the engine scenario.","tokens_in":12131,"tokens_out":3665,"would_cite":false,"duration_ms":40421,"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 radio-optical class of a radio galaxy — FRI versus FRII, high- versus low-excitation — is already visible in the parsec-scale jet, before it interacts with the kiloparsec-scale environment.","keywords":["radio galaxies","Fanaroff-Riley classification","high-excitation radio galaxies","low-excitation radio galaxies","very long baseline interferometry","parsec-scale jets","accretion modes","3CR sample"],"falsifier":"Run a complete 8 GHz VLBI survey of all 79 3CR radio galaxies at z < 0.3 without preselection by compactness; if the 24 sources currently missing from the VLBI subsample are included and the FRII-LERG and FRII-HERG luminosity distributions then differ significantly, the paper's central claim would be falsified.","tokens_in":11031,"feed_emoji":"📡","tokens_out":8975,"duration_ms":101654,"temperature":0.7,"pith_summary":"This paper tests three competing explanations for why some radio galaxies have edge-brightened FRII morphology yet inefficient, low-excitation accretion. It does so by measuring the 8 GHz radio luminosity of the jets at parsec scales using very long baseline interferometry, for 55 galaxies from the 3CR sample with known optical-radio classes. The key result is that the classes keep the same relative ordering in the parsec-scale luminosity–[OIII] plane as they do in the kiloparsec-scale 178 MHz plane, and FRII galaxies with efficient accretion (HERGs) and inefficient accretion (LERGs) have statistically indistinguishable 8 GHz luminosities. A sympathetic reading is that the paper establishes that the engine, not the environment or later fuel depletion, sets the radio-optical class: the distinction is already present within about 10 pc of the black hole.","feed_headline":"Radio galaxy class is set at the jet base, 3CR VLBI data show","feed_subtitle":"In 55 3CR galaxies, 8 GHz core luminosities match kiloparsec classes, pointing to the jet engine itself.","key_machinery":"The load-bearing comparison is between two radio luminosities of the same objects at different spatial scales: the parsec-scale 8 GHz luminosity $L_{8\\,\\mathrm{GHz}}$ measured by VLBI and the kiloparsec-scale 178 MHz luminosity from the 3CR survey, both plotted against the [OIII]$\\lambda5007$ luminosity that serves as an accretion-rate proxy. The paper also uses the compactness parameter $C_{8\\,\\mathrm{GHz}} = f_{\\mathrm{un}}/f_{\\mathrm{tot}}$, the fraction of the 8 GHz flux that is unresolved by VLBI, to check that the observed subsample is not biased toward a narrow range of orientations. Because the milliarcsecond beam corresponds to a few parsecs at the sample's redshifts, the 8 GHz luminosity reads the jet state near its launch region, so the matching class ordering at 8 GHz and 178 MHz becomes evidence about the central engine rather than about large-scale propagation.","core_discovery":"The discovery is that the separation between radio-optical classes is already evident at milliarcsecond scales, before the jet has propagated far enough to interact with the kiloparsec-scale medium. For 52 of the 55 sources, the 8 GHz VLBI luminosity maps emission from regions smaller than 10 pc (uncorrected for projection) in 75% of cases, close to the Bondi radius. Plotted against [OIII] luminosity, the classes arrange themselves at 8 GHz just as they do at 178 MHz: FRIs faint in both, FRII-HERGs bright in both, and FRII-LERGs intermediate. The Mann-Whitney and Kolmogorov-Smirnov tests give $p_0 = 0.3$ and $p_0 = 0.4$ for the FRII-LERG versus FRII-HERG luminosity comparison, so the paper finds no statistical evidence that different accretion regimes produce different parsec-scale jet luminosities. The authors conclude that FRII-LERGs are not faded former HERGs and that the environment does not create the FRI/FRII difference; instead, the differentiation happens during the initial ejection phase, driven by the accretion-disk type and black hole properties.","pith_inferences":["A direct testable extension is to measure jet collimation profiles within the first few tens of Schwarzschild radii for these same objects: if the engine scenario is correct, FRII-LERGs and FRII-HERGs should show different expansion shapes even where their 8 GHz luminosities overlap.","The luminosity overlap at 8 GHz does not by itself prove equal kinetic power; high-frequency VLBI polarization and spectral-index maps could test whether FRII-LERG jets carry the same magnetic flux density as FRII-HERGs.","Applying the same 8 GHz versus [OIII] comparison to complete flux-limited samples at millijansky levels would show whether the engine scenario governs the whole radio-galaxy population or only the bright 3CR subset."],"forward_implications":["FRII-LERGs are unlikely to be former FRII-HERGs whose fuel has run out: a faded accretion flow would leave a weak parsec-scale jet, but these jets are as luminous as those of actively accreting FRII-HERGs.","The kiloparsec-scale environment cannot be the main agent that separates FRI from FRII and HERG from LERG, because the classes are already separated at the jet base.","The engine scenario is the remaining explanation: black hole mass, spin, and magnetic flux, together with the accretion-disk configuration, determine the radio-optical class during the initial ejection phase.","The same relative ordering at 8 GHz and 178 MHz means the parsec-scale core luminosity carries information about the long-term jet power, which can be used where extended lobe emission is not detectable."],"supporting_citations":[{"why":"Provides the VLBI positions and flux densities from the Fundamental Radio Catalog used to compute the parsec-scale 8 GHz luminosities.","marker":"Petrov & Kovalev 2025"},{"why":"Supplies the optical and radio classifications and the [OIII] and 178 MHz luminosities for the 3CR sample.","marker":"Buttiglione et al. 2009"},{"why":"Completes the optical-radio classification of the 3CR sources used in this analysis.","marker":"Buttiglione et al. 2010"},{"why":"Extends the spectroscopic classification that defines the HERG and LERG classes.","marker":"Buttiglione et al. 2011"},{"why":"The X-ray study that proposed the evolutionary and environmental scenarios the paper sets out to test.","marker":"Macconi et al. 2020"},{"why":"The earlier millijansky study that introduced the engine scenario by showing similar inefficient accretion can produce different radio morphologies and luminosities.","marker":"Grandi et al. 2021"},{"why":"Defines the FRI/FRII morphological classification adopted throughout the paper.","marker":"Fanaroff & Riley 1974"},{"why":"Introduces the high- and low-excitation optical classification (HERG and LERG) used to separate accretion regimes.","marker":"Jackson & Rawlings 1997"},{"why":"Defines the compactness parameter $C_{8\\,\\mathrm{GHz}} = f_{\\mathrm{un}}/f_{\\mathrm{tot}}$ used to check orientation bias in the VLBI subsample.","marker":"Popkov et al. 2021"},{"why":"Linking accretion mode to VLBI jet expansion profiles supports the interpretation that class differences emerge at the jet base.","marker":"Boccardi et al. 2021"}],"fun_headline_variants":["Radio class fixed at jet base, parsec-scale data show","Class separation already at parsec scale in 3CR galaxies","Jet engine sets radio class, VLBI data reveal","FRII differences visible before kiloparsec scales","Accretion type imprinted on parsec-scale jet luminosity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes the 55 galaxies with high-resolution radio observations fairly represent all 79 galaxies: if compact or bright jets were more likely to be observed, the apparent similarity of the two FRII classes could be an artifact of which sources got observed.","fun_headline_variants_meta":{"raw":{"variants":["Radio class fixed at jet base, parsec-scale data show","Class separation already at parsec scale in 3CR galaxies","Jet engine sets radio class, VLBI data reveal","FRII differences visible before kiloparsec scales","Accretion type imprinted on parsec-scale jet luminosity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00068,"raw_usage":{"total_tokens":3169,"prompt_tokens":1103,"completion_tokens":2066,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":719,"completion_tokens_details":{"reasoning_tokens":1987}},"tokens_in":719,"tokens_out":2066,"duration_ms":16759,"temperature":1.0,"reasoning_tokens":1987,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:29:40.150644+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a complete 8 GHz VLBI survey of all 79 3CR radio galaxies at z < 0.3 without preselection by compactness; if the 24 sources currently missing from the VLBI subsample are included and the FRII-LERG and FRII-HERG luminosity distributions then differ significantly, the paper's central claim would be falsified.","supporting_citations":[{"cited_title":"2009, A&A, 495, 1033","cited_arxiv_id":null,"evidence_quote":"Supplies the optical and radio classifications and the [OIII] and 178 MHz luminosities for the 3CR sample."},{"cited_title":"2011, A&A, 525, A28+","cited_arxiv_id":null,"evidence_quote":"Extends the spectroscopic classification that defines the HERG and LERG classes."},{"cited_title":"2020, MNRAS, 493, 4355","cited_arxiv_id":null,"evidence_quote":"The X-ray study that proposed the evolutionary and environmental scenarios the paper sets out to test."},{"cited_title":"2021, ApJ, 911, 17","cited_arxiv_id":null,"evidence_quote":"The earlier millijansky study that introduced the engine scenario by showing similar inefficient accretion can produce different radio morphologies and luminosities."},{"cited_title":"& Rawlings, S","cited_arxiv_id":null,"evidence_quote":"Introduces the high- and low-excitation optical classification (HERG and LERG) used to separate accretion regimes."}],"review_version":1}