{"id":"e912c491-583b-4968-954f-e68ed6adaeca","arxiv_id":"2411.09099","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"VLBI radio jet position angles preferentially align with the projected minor axes of their host galaxies, with p <= 0.01 in 5,853 DESI Legacy Survey matches that have spectroscopic redshifts.","lead":"Astronomers compared the directions of radio jets from thousands of supermassive black holes with the shapes of their host galaxies. They find a weak but statistically significant tendency for the jets to point along the shortest axis of galaxies that have measured distances.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The signal appears only in the spec-z subsample, and the permutation null preserves marginal PA distributions but not selection correlations with redshift or sky position; this selection channel is the load-bearing untested assumption.","rationale":"I agree with the reader's weakest assumption. The strongest evidence that this is load-bearing is the internal contrast in Fig. 3: the full DESI LS sample shows no alignment, and the signal appears only in the spectroscopic-redshift rows. A physical interpretation based on resolution or redshift is plausible, but the current permutation test cannot exclude the alternative that spec-z availability or shape-quality selection correlates both PAs with a common variable such as declination or redshift. Since the same data and code can implement a stratified permutation null, this concern is directly testable. I do not see an internal inconsistency in the statistical construction, and the paper is transparent about its cuts and systematics. The proposed test is cheap and decisive; until it is run, or an equivalent check is provided, the CONDITIONAL verdict remains appropriate. The reader's verdict should not change because the same condition is already identified, and no new fatal flaw has emerged.","tokens_in":32355,"tokens_out":8973,"duration_ms":116081,"concrete_test":"Recompute the DESI LS p-values for the spec-z and spec-z/good-case cells using a conditional permutation null in which jet PAs are shuffled only within narrow bins of declination (e.g. 10-degree bins), spectroscopic redshift or z-band magnitude (as a redshift proxy), and optical semi-minor axis or ellipticity, keeping the selected sample fixed. If the p-values remain below 0.01 under this stratified null, the selection-bias explanation is ruled out; if they rise above 0.05, the reported alignment is a selection artifact rather than an astrophysical signal.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the contrast between Fig. 3a-c (p5 ~ 0.5-0.6, no alignment) and Fig. 3d-l (p <= 0.01) in DESI LS: the signal is introduced by requiring a measured spectroscopic redshift and, in the right column, a good optical shape. The Q statistic and reshuffling bands correctly account for non-uniform PA marginals by shuffling the observed PA lists within the selected sample. They do not account for selection that is jointly dependent on a third variable. Spec-z availability from OCARS is not random: it correlates with redshift, magnitude, optical spectral properties, and survey footprint. If, within the spec-z/good-case subsample, the jet PA has a residual systematic dependence on declination or redshift (the paper itself discusses a known VLBI north-south elongation effect) and the optical PA has any similar dependence from PSF or survey systematics, the two PA lists can become correlated through that common variable even with no intrinsic jet-galaxy alignment. The reshuffling null breaks the PA-PA correlation in the observed sample but also destroys PA-declination and PA-redshift correlations, so it cannot distinguish these cases. The paper does not report a declination-stratified or redshift-stratified null for the spec-z/good-case cells, nor a multiple-testing correction across the 12 panels and the additional survey, binning, and magnitude cuts. The headline p <= 0.01 is therefore not yet shown to be robust against this selection channel.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript cross-matches roughly 6273 VLBI-detected AGN with optical survey catalogues to compare parsec-scale jet position angles with the projected kpc-scale shapes of their host galaxies. After applying resolution and quality cuts, the authors use a Monte Carlo reshuffling null, KS tests, and robustness checks on binning and magnitude cuts to search for an excess of small Delta PA between the VLBI jet and the optical minor axis. They find no alignment in the full DESI LS sample (Fig. 3a-c, p5 ~ 0.55-0.62), but report a weak signal with p <= 0.01 in the subset with measured spectroscopic redshifts (Fig. 3d-l). The same analysis is applied to a combined sample and to individual surveys, and an EAGLE projection model is used to argue that the observed Delta PA distribution is consistent with a ~33-degree scatter about a close 3D spin-minor-axis relation. The central claim is that the pc-scale jet direction is statistically tied to the kpc-scale projected minor axis of the host galaxy over three orders of magnitude in scale, at least for sources with well-measured shapes and spectroscopic redshifts.","tokens_in":32723,"tokens_out":8329,"duration_ms":88679,"significance":"If the result is robust, it would be an interesting and rare observational link between the innermost AGN jet direction and the large-scale stellar morphology of the host galaxy, with implications for merger-driven galaxy evolution and SMBH-galaxy co-evolution. The study is based on a large VLBI sample, uses independent jet and optical catalogues, and ships its cross-matching code and data publicly. The statistical machinery is carefully presented: the reshuffling null preserves the marginal PA distributions, the binning and magnitude checks are useful, and the EAGLE projection exercise is a sensible way to illustrate projection effects. The main weakness is that the headline signal appears only after selecting sources with spectroscopic redshifts and good shapes, and the analysis does not yet rule out the possibility that the signal is introduced by selection correlations with redshift, declination, or survey systematics rather than by an intrinsic alignment. The cross-survey comparison is also weaker than the DESI LS result alone, so the robustness of the central claim is not yet established.","major_comments":[{"comment":"The headline signal is introduced by the spectroscopic-redshift selection: Fig. 3a-c are consistent with the null (p5 = 0.62, 0.58, 0.55; pks ~ 0.09-0.11), while Fig. 3d-l reach p5 <= 0.01. The Q-statistic null reshuffles the two PA lists within the selected subsample (Eqs. 6-7), which preserves the marginal PA distributions but removes any joint dependence of jet PA and optical PA on a third variable. Spectroscopic-redshift availability in OCARS is not random: it correlates with redshift, magnitude, and survey footprint, and the paper itself documents a VLBI north-south elongation effect and non-uniform optical PAs in the section 'Underlying distribution of the galaxy and jet PAs'. A residual common dependence of the two PAs on declination or redshift within the spec-z/good-case subsample would therefore produce a Delta PA excess under the null. Please provide declination- and redshift-stratified null tests (or an equivalent conditional analysis) to close this selection channel.","section":"Fig. 3 and Methods 'Quantifying the alignment signal' (Eqs. 6-7)"},{"comment":"No multiple-testing correction is applied across the 12 panels of Fig. 3, the alternative bin numbers in M7, the magnitude cuts in M8, and the survey-by-survey and combined analyses in M3-M6. With the top row null and p-values of order 0.01 in the spec-z rows, a look-elsewhere effect is a plausible alternative to an astrophysical alignment. The authors should report the effective number of independent tests or demonstrate a pre-specified hierarchy of cuts in which the signal is present; the qualitative statement that 'the overall picture' is unchanged by binning does not provide a quantitative correction.","section":"Fig. 3 and Supplementary Figs. M3-M8"},{"comment":"The cross-survey comparison does not currently support a universal signal. DESI LS spec-z rows give p5 ~ 0.005-0.01, but the combined sample in M3 gives much weaker p-values (e.g., M3d p5 = 0.05, M3l p5 = 0.08), DES and SkyMapper in M5-M6 are consistent with the null in most rows, and the only independent significant detection is SDSS at z < 0.1 with N = 52 (M4j-l, p5 = 0.006). The text explains that DESI LS is used because it is less affected by systematics, but a quantitative error budget for why the other surveys lose the signal is needed before the survey dependence can be attributed to systematics rather than to the same selection channel that may be at work in DESI LS.","section":"Supplementary Figs. M3-M6 and 'Survey by survey analysis'"}],"minor_comments":[{"comment":"The conclusion states that the jet is 'oriented perpendicularly to the projected optical shape of the host galaxy with a p-value < 0.05' without the qualifiers 'in sources with measured spectroscopic redshifts' and 'weak'; please align this wording with the abstract and with the null results in the top row of Fig. 3.","section":"Conclusion"},{"comment":"The phrase 'perpendicularly to the projected optical shape' is ambiguous: the measured Delta PA is relative to the projected minor axis, so the signal is an alignment with the minor axis (equivalently a perpendicular orientation to the galaxy disk plane), not a perpendicular orientation to the projected shape as a whole. Consider rewording for clarity.","section":"Abstract and Fig. 1"},{"comment":"The Fig. 3 caption repeats the sentence 'Each panel corresponds to a specific selection of sources' twice, and the KiDS row of Table 1 appears to be missing a column value; the table should be formatted consistently.","section":"Fig. 3 caption and Table 1"},{"comment":"pks values of 0.000 (e.g., Fig. 3g) should be reported as p < 0.001, given that the p-value resolution is set by 1000 reshufflings.","section":"Fig. 3 and Fig. M7"},{"comment":"The EAGLE scatter amplitude epsilon = 0.33 is chosen after the fact to match the observed 68% interval; the main text does acknowledge this, but the Fig. 4 caption should explicitly mark panels (c) and (d) as illustrative rather than as a quantitative fit to the data.","section":"Fig. 4"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth taking seriously, but the headline result is not yet nailed down. The genuinely new thing here is scale: roughly 5,800 DESI LS cross-matches with parsec-scale VLBI jet position angles, compared against earlier work with a handful of pairs or kpc-scale radio emission. The multi-survey comparison is a real step up, and the authors are unusually transparent about systematics—PA uncertainties, seeing cuts, binning checks, magnitude cuts, even the non-uniform optical PA distributions in SkyMapper. The reshuffling null and KS tests are sensible, and shipping the cross-match data and code is good practice. Give credit where it is due: this is the first statistical look at jet–minor-axis alignment at VLBI resolution, and the analysis is mostly careful.\n\nThe soft spot is the one the reader flagged, and I think it lands. The signal is absent in the full DESI LS sample and appears only once you require a spectroscopic redshift and, in the right column, a well-measured shape. The permutation null correctly breaks PA–PA correlations within the selected sample, but it also breaks any correlations between PA and declination or redshift. Spec-z availability is not random—it tracks redshift, magnitude, spectral type, and survey footprint—and the paper itself notes a known north-south VLBI PA effect. If jet PA and optical PA both have weak systematics tied to sky position or redshift, the reshuffling null will happily call the induced correlation a detection. The paper does not report a declination-stratified or redshift-stratified null for the spec-z/good-case cells, and it does not correct for the multiple redshift, quality, and binning cuts examined across twelve panels and several surveys. That is a genuine gap, not a manufactured one.\n\nThat said, the flaw is in the significance claim, not the method. The underlying observable—an excess of small ΔPA in the low-redshift, spectroscopically confirmed subset—is worth checking with a better null. The EAGLE projection stuff is illustrative rather than evidential; the epsilon=0.33 scatter is tuned to reproduce the observed median and is not independent confirmation. I would not treat it as support for the detection.\n\nNet: a solid, reproducible observational study with a load-bearing selection assumption that remains untested. The paper deserves referee time and should be published after the authors either add the stratified nulls and multiple-testing checks or soften the claim to a candidate signal. If the spec-z selection effect is real, the p<=0.01 will not survive; if it is not, this becomes a useful anchor for the field. Worth a serious referee, but I would not cite it as an established alignment yet.","headline":"A careful large-sample search for a jet–host alignment that finds a real candidate signal, but the headline p-value appears only after selection cuts that the null does not yet control for.","tokens_in":33240,"tokens_out":1122,"would_cite":false,"duration_ms":17243,"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 finds that parsec-scale AGN jets appear statistically perpendicular to the projected minor axes of their host galaxies, connecting scales from sub-parsec to kiloparsec.","keywords":["AGN jets","VLBI","host galaxy alignment","galaxy morphology","supermassive black holes","position angles","radio galaxies"],"falsifier":"Compute the same $\\Delta PA$ histogram for sources that have only photometric redshifts while applying the identical shape-quality cuts; if the perpendicular alignment is a real astrophysical effect it should survive, whereas if it disappears the reported signal is a selection artifact tied to spectroscopic-redshift availability.","tokens_in":32110,"feed_emoji":"🔭","tokens_out":8257,"duration_ms":82548,"temperature":0.7,"pith_summary":"The paper sets out to test whether the direction of the relativistic jet launched near a supermassive black hole carries any information about the optical shape of its host galaxy. The authors compare VLBI jet position angles for more than six thousand AGN against galaxy shapes from several optical surveys, using the survey least affected by systematics as the main sample. In the subset of sources with measured spectroscopic redshifts and reliable shape measurements, they find a weak but significant excess of jets oriented perpendicular to the projected minor axis of the galaxy, with $p \\leq 0.01$ in the primary sample. No alignment appears when the full sample is used without the redshift and quality cuts. If the effect is real, it implies that the black-hole–disk–jet system and the kiloparsec-scale stellar body of the galaxy are not independently oriented.","feed_headline":"Parsec-scale jets align with host galaxies' minor axes","feed_subtitle":"VLBI jets tend to point perpendicular to their host galaxy's short axis, hinting at black hole–galaxy co-evolution.","key_machinery":"The central object is the angular separation $\\Delta PA = |PA_{\\rm jet} - PA_{\\rm minor}|$ wrapped into $[0^\\circ, 90^\\circ]$, measuring how close the jet is to perpendicular to the galaxy's projected minor axis. Jet position angles come from VLBI images with formal errors; optical minor axes come from elliptical model fits to galaxy images. The null hypothesis is built by reshuffling the jet and optical position-angle distributions 1000 times, producing a reference band and a statistic $Q$; a Kolmogorov–Smirnov test against a uniform distribution is used as a binning-independent check. The load-bearing cuts are a semi-minor axis larger than $1.3$ arcsec, position-angle errors below $22.5^\\circ$, survey-specific 'good case' morphological flags, and the presence of a spectroscopic redshift. A projection simulation using a cosmological hydrodynamical simulation shows that projection alone cannot produce the observed peak unless the stellar spin is closely related to the galaxy plane.","core_discovery":"The central claim is that the parsec-scale jet direction is statistically tied to the kiloparsec-scale optical shape of its host galaxy: jets preferentially point along the projected minor axis of the optical image, i.e., perpendicular to the major axis. The evidence comes from the distribution of the angle difference $\\Delta PA$ between the jet position angle and the closest minor-axis position angle, which peaks toward small $\\Delta PA$ only after restricting to sources with spectroscopic redshifts and well-measured, non-spherical shapes. The reported significance is $p \\leq 0.01$ for the primary optical survey sample (and $p < 0.05$ in the conclusion), with consistent results from two-bin, five-bin, and Kolmogorov–Smirnov tests. The paper interprets the absence of signal in the unrestricted sample as a resolution and selection effect rather than as evidence against a physical connection.","pith_inferences":["Because several redshift, quality, and binning choices were examined, the reported $p \\leq 0.01$ likely overstates the evidence unless corrected for multiple comparisons; re-running the same statistic on an independent sample would settle the matter.","If the perpendicular alignment is real, it suggests that the jet's launch direction is not determined purely by local accretion-disk dynamics but is correlated with the galaxy-scale angular momentum distribution, with consequences for feedback and merger models.","A clean test would be to split the sample by host morphology and radio luminosity; the merger-based interpretation predicts the alignment should be strongest in elliptical galaxies formed by major mergers and weaker in disk galaxies."],"forward_implications":["A genuine pc–kpc alignment would connect sub-parsec jet physics to kiloparsec-scale stellar structure, supporting the idea that supermassive black holes and their host galaxies co-evolve.","The signal's appearance only in the spectroscopic-redshift, well-shaped subsample implies that source resolution and shape reliability control the detectability of the alignment in current data.","The projection simulations indicate that the observed peak is stronger than projection effects alone can produce, so a physical relation between the black-hole spin axis and the galaxy's angular momentum is the natural explanation.","Larger samples from upcoming wide-area radio and optical surveys should sharpen the significance and test whether the alignment persists at higher redshift and in different galaxy types."],"supporting_citations":[{"why":"Supplies the VLBI jet position angles and formal errors for 9220 AGN, the jet dataset used throughout.","marker":"[46]"},{"why":"Provides the optical galaxy shapes and morphological fits that constitute the primary comparison sample.","marker":"[56]"},{"why":"Earlier kpc-scale radio–optical orientation comparison whose perpendicular result is extended here to parsec scales.","marker":"[15]"},{"why":"Documents a north-south elongation systematic in VLBI jet position angles that the analysis must account for.","marker":"[65]"},{"why":"Supplies the cosmological simulation whose galaxies are used to model projection effects.","marker":"[22]"},{"why":"Companion simulation calibration paper used in generating the projected galaxy sample.","marker":"[23]"},{"why":"Provides the catalogue of spectroscopic redshifts used to define the redshift-restricted subsample.","marker":"[47]"},{"why":"Earlier VLBI-host-galaxy alignment study with a small sample that this work builds on.","marker":"[6]"}],"fun_headline_variants":["Jets align with galaxy minor axes","Black hole jets point along galaxy short axes","AGN jets favor host minor axis","Parsec jets track kpc galaxy shape","Jets and galaxies: perpendicular tie"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the subsample with measured spectroscopic redshifts and reliable shapes is not biased in a way that manufactures the angle correlation, and the analysis does not correct for the number of redshift, quality, and binning choices examined.","fun_headline_variants_meta":{"raw":{"variants":["Jets align with galaxy minor axes","Black hole jets point along galaxy short axes","AGN jets favor host minor axis","Parsec jets track kpc galaxy shape","Jets and galaxies: perpendicular tie"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000161,"raw_usage":{"total_tokens":1244,"prompt_tokens":965,"completion_tokens":279,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":581,"completion_tokens_details":{"reasoning_tokens":216}},"tokens_in":581,"tokens_out":279,"duration_ms":3980,"temperature":1.0,"reasoning_tokens":216,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:02:00.543950+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same $\\Delta PA$ histogram for sources that have only photometric redshifts while applying the identical shape-quality cuts; if the perpendicular alignment is a real astrophysical effect it should survive, whereas if it disappears the reported signal is a selection artifact tied to spectroscopic-redshift availability.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier VLBI-host-galaxy alignment study with a small sample that this work builds on."}],"review_version":1}