Pith. sign in

REVIEW 4 major objections 5 minor 51 references

Detailed analysis of multi-line molecular distributions in the Seyfert galaxy NGC 1068: Possible effect of the AGN outflow to the starburst ring

T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A low-variance component of 13 molecular-line maps aligns with the AGN outflow and with broad line widths, suggesting the outflow reaches the starburst ring at 1–2 kpc.

desk verdict Solid PCA application with an honest but unproven outflow-ring claim; deserves review with a mandatory permutation test. read the letter →

arxiv 2506.15999 v1 pith:4NZ5KDLP submitted 2025-06-19 astro-ph.GA

Hiroma Okubo , Toshiki Saito , Shuro Takano , Nario Kuno , Akio Taniguchi , Taku Nakajima , Nanase Harada , Ken Mawatari This is my paper · ORCID
classification astro-ph.GA
keywords principalcomponentanalysisNGC1068molecularlinemappingAGNoutflowstarburstringinterstellarmediumchemistrygalacticfeedback
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper tries to show that a blind statistical decomposition of 13 molecular-line maps can separate the overlapping physical processes in NGC 1068, a galaxy that hosts both an active galactic nucleus and a starburst ring. Its central result is the third principal component computed only on the ring region: negative scores lie along the AGN outflow direction, are carried mainly by CN, C2H, HCN, HCO+, and CS, and the velocity dispersions of those lines anticorrelate with the negative scores. The authors read this as evidence that the AGN outflow is reaching the starburst ring, roughly 1 to 2 kpc from the nucleus, rather than stopping inside the central kiloparsec. If correct, it means low-variance components of multi-line PCA can expose feedback structures that are hard to isolate by eye in complex galaxies.

What carries the argument

The central object is the PC3SB score map: the third principal component computed from 13 standardized integrated-intensity maps of molecular lines on a hexagonal cloud-scale grid covering the starburst ring, an annulus from 750 pc to 2 kpc. Standardization, subtracting each line's mean and dividing by its standard deviation, equalizes the variance across lines so that bright CO emission cannot dominate the decomposition. PCA projects each spatial grid point into a low-dimensional space whose axes are linear combinations of the line intensities, and the load-bearing diagnostic is the anticorrelation between negative PC3SB scores and the velocity dispersions of the negatively loaded lines; that correlation turns a spatial coincidence into a dynamical argument for outflow-driven broadening.

What would settle it

A decisive test would be to measure the kinematics of the southwestern ring at high resolution: if the broad CN, C2H, HCN, and CO line widths in the negative-PC3SB regions are not spatially tied to the outflow cone and do not show outflow-like acceleration along the cone axis, then the feature is more likely tracing bar-end star formation or cloud-cloud collisions than AGN outflow impact. Re-running the PCA with min-max normalization, line-ratio inputs, or no scaling would also settle whether the outflow-aligned pattern is intrinsic or a byproduct of the standardization choice.

Watch

Extended reading notes

Core claim

The central claim is that PC3SB, the third principal component of the starburst-ring maps, which retains about 2% of the variance, isolates a spatial pattern best explained by the AGN outflow interacting with molecular gas in the ring. The negative PC3SB scores fall along the known outflow direction; the lines dominating the negative coefficients are CN, C2H, and HCN, with HCO+ and CS also contributing, exactly the species expected to be enhanced or heated by an outflow; and the velocity dispersions of the negatively loaded lines anticorrelate with the negative scores, with most correlation coefficients more negative than -0.5 and C18O the exception. The paper is careful that this component may contain several physical signals and that it overlaps the bar-end, but it concludes that one coherent interpretation is outflow compression of the ring gas. It also claims that PC1 approximates the H2 column density, that PC2OA separates the circumnuclear disk from the starburst ring, and that PC2SB separates starburst-dominated from shock-dominated gas.

Load-bearing premise

The load-bearing premise is that standardizing each line map by variance is the physically appropriate scaling; when the maps are instead rescaled to a fixed 0–1 range, the southwest–northeast contrast in PC3SB weakens and the analogous component in the overall region becomes dominated by the bar-end rather than the outflow.

Editorial extensions

If this is right

  • AGN outflows can influence molecular gas at radii of 1 to 2 kpc, extending the known scale of such interactions beyond the central kiloparsec.
  • Principal components that retain only about 2% of the variance can still carry a physically interpretable signal, so truncating at PC1 or PC2 would discard the outflow signature.
  • Multi-line PCA can separate overlapping components of column density, chemical state, and feedback in galaxies with both an AGN and a starburst ring.
  • The anticorrelation between a low-variance PC score and line velocity dispersion offers a transferable kinematic test for outflow-cloud interactions in other galaxies.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, the prominence of the outflow-aligned feature under standardization suggests it may partly encode dynamic-range differences between bright CO and faint dense-gas tracers; re-running the analysis on line ratios or with explicit intensity scaling would separate that effect from an intrinsic ISM structure.
  • A natural prediction is that other nearby Seyfert galaxies with millimeter-wave line maps will show low-variance components aligned with their radio or X-ray outflow axes, loaded on CN, C2H, and HCN, and correlated with broad line widths.
  • If the interaction is real, the southwestern ring should show independent corroborating signs, such as enhanced dust temperature, SiO or HNCO shock emission, or elevated CN/CO and C2H/CO ratios, which existing or future observations could test.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. This paper applies principal component analysis (PCA) to ALMA Band 3 integrated intensity maps of 13 molecular lines toward the Seyfert galaxy NGC 1068, separately for an overall central region within about 2 kpc and for the starburst ring between 750 pc and 2 kpc. The authors standardize each line map to zero mean and unit variance, use the scikit-learn PCA implementation, and assign coefficient errors by Monte Carlo perturbation of the noise. They interpret PC1 as an approximate H2 column density tracer, PC2 as a density/chemical contrast between the circumnuclear disk and the starburst ring, and, as the central claim, PC3SB negative scores as a possible signature of interaction between the AGN outflow and gas in the starburst ring. This interpretation is supported by the spatial coincidence of negative scores with the AGN outflow direction, by the loading of CN, C2H, HCN, HCO+, CS, and CO isotopologues on those scores, and by an anti-correlation between negative scores and the velocity dispersions of several lines (Table 5). An appendix repeats the analysis with min-max normalization and reports notable differences in the PC3OA and PC3SB maps.

Significance. If the central claim holds, it would extend the influence of the AGN outflow in NGC 1068 from the central kiloparsec to the starburst ring at about 1–2 kpc, and it would demonstrate that low-variance PCA components of multi-line maps can reveal physically meaningful structure in a galaxy hosting both an AGN and a starburst ring. The paper has clear strengths: the data processing is described in detail, the analysis uses a uniform PHANGS-ALMA pipeline with a common grid and mask, the Monte Carlo error bars on the PCA coefficients are a useful addition, and the PC2OA interpretation is checked against lines not used in the PCA (Table 4). The velocity-dispersion test in Table 5 is a genuine external check. However, the central claim currently rests on a component that explains only about 2% of the variance, whose spatial pattern is sensitive to the scaling choice, and for which no null model or quantitative alignment test is provided.

major comments (4)
  1. [Appendix, Figs. 13(d), 14] The robustness of the central claim to the choice of scaling is not established. The appendix states that with min-max normalization the PC3SB negative-score contrast between the southwestern and northeastern parts of the ring becomes less pronounced, and that the PC3OA negative scores primarily extract the bar-end rather than the AGN outflow (Figs. 13d and 14). This is load-bearing because standardization is an arbitrary preprocessing choice; the outflow-aligned PC3SB feature may reflect the scaling rather than an intrinsic ISM structure. Please quantify the southwest–northeast contrast under both scalings and under additional choices (for example, median/RMS normalization), give a physical justification for the chosen scaling, or substantially temper the central claim.
  2. [Sec. 2.3, Fig. 9(a)] No null model is presented for PC3SB. The Monte Carlo analysis in Sec. 2.3 perturbs each pixel by its noise and shows that the coefficient signs are stable against random noise, but it does not test whether a principal component with only about 2% of the variance and with a spatially coherent pattern would arise from 13 independent maps with the same spatial autocorrelation. The SB ring region contains only 482 hexagons smoothed to 150 pc (Sec. 2.1), so the effective number of independent samples is much smaller than 482. Please add a permutation or surrogate test that destroys the spatial correspondence between the line maps while preserving each map’s autocorrelation, and report the null distributions of the PC3 variance fraction, the overlap of negative scores with the outflow cone, and the Table 5 correlations.
  3. [Sec. 4.2.3, Fig. 9(d)] The spatial alignment of the negative scores with the AGN outflow is assessed only visually. Because this alignment is the primary spatial evidence for the interaction claim, a quantitative statistic is needed: for example, the fraction of negative-score pixels inside the outflow cone compared with a random control, or a two-dimensional cross-correlation with an outflow mask. The acknowledged overlap with the bar-end region should also be separated or explicitly modeled, since the bar-end is a plausible alternative origin for the negative scores.
  4. [Table 5, Sec. 4.2.3] The anti-correlation between PC3SB negative scores and line velocity dispersions is reported without uncertainties and without correction for spatial autocorrelation. With correlation coefficients between about –0.26 and –0.65 across nine lines, and with multiple comparisons, these values are not sufficient to establish the claimed effect by themselves. Please provide bootstrap or jackknife confidence intervals, account for the effective number of independent spatial samples, and report the correlation over all ring pixels rather than only over the negative-score-selected pixels, to avoid selection effects.
minor comments (5)
  1. [Figure captions 2, 4, 5] The captions contain the typo “usd” for “used”; please correct this.
  2. [Table 4 heading] The word “Talbe” in the table heading should be “Table”.
  3. [Table 5 header] The column header for Table 5 is garbled: the labels for the two CN hyperfine components and for C2H are not cleanly aligned with the data columns, so the reader cannot tell which transition corresponds to which coefficient. Please reformat the table.
  4. [Sec. 4.2.3] The text refers to “C2H (N=11/2-01/2)”, but Table 5 lists simply “C2H”; please use consistent notation for all lines.
  5. [Abstract and Sec. 4.2.2] The phrase “have a possibility to reconstruct” is awkward and should be rewritten as “may reconstruct”; also, “collisions between molecular molecular clouds” in Sec. 4.2.2 repeats “molecular”.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: PC3SB interpretation is a descriptive PCA result cross-checked against external tracer calibrations and independent velocity-dispersion data.

full rationale

The paper is descriptive rather than derivational: it applies PCA to 13 standardized integrated-intensity maps (Section 2.2) and interprets the resulting score maps using tracer assignments cited from prior literature. The central PC3SB claim—that the negative-score pattern is consistent with AGN outflow interaction—is not forced by construction. The negative scores are a linear combination of the input intensities, while the supporting velocity-dispersion anti-correlation in Section 4.2.3 and Table 5 uses measurements not entering the PCA, so the check is external to the derived component. The AGN outflow geometry is anchored to independent references (Das et al. 2006; Mingozzi et al. 2019). The self-citation to Saito et al. (2022b), which shares co-authors, is used as a prior chemical calibration for CN/C2H enhancement by UV/AGN outflow, but the same association is also supported by García-Burillo et al. (2017), and it is not the sole evidence for the interaction claim. The scaling sensitivity shown in the Appendix weakens the robustness of low-variance PCs, but that is a statistical concern, not a circularity in the derivation. No fitted parameter is renamed as a prediction, and no quantity is defined in terms of the conclusion it is used to support.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new physical entities or fitted physical constants. Its central claim rests instead on data-processing choices such as S/N thresholds, grid scale and standardization, plus literature tracer assignments and the linearity and independence assumptions of PCA.

free parameters (4)
  • S/N threshold for CO mask = 4
    Grids with CO(1-0) S/N <= 4 are removed before PCA (Section 2.1); this choice changes the sample and therefore the PCs.
  • Minimum detected grid count = 60
    Lines must have S/N > 4 in at least 60 hexagons to enter the PCA (Section 2.1); this is an ad hoc sample selection that removes weak lines.
  • Integration velocity range = +/-300 km/s
    Maps integrate +/-300 km/s around 1116 km/s (Section 2.1), which sets which emission enters every map.
  • Spatial grid scale = 150 pc
    Data are regridded to 150 pc hexagonal cells (Section 2.1); the resolution determines the effective number of independent samples and the PCA outcome.
assumptions (4)
  • domain assumption PCA assumes linear relationships and independent samples; each hexagon is treated as one molecular cloud with no spatial correlation.
    Section 2.2 states this assumption; if nonlinearity or correlated pixels dominate, PC3 can mix or invent features.
  • domain assumption Standardization (zero mean, unit variance per line) is the appropriate scaling for comparing lines.
    This choice is central to PC3SB; the appendix shows min-max normalization weakens the outflow-aligned negative scores.
  • domain assumption Velocity dispersion broadening of molecular lines indicates interaction or shock rather than rotation, optical depth, or hyperfine structure.
    Section 4.2.3 uses anticorrelation with velocity dispersion to support the outflow interaction; C2H hyperfine structure is noted but not fully corrected.
  • domain assumption Literature tracer assignments, including CN and C2H enhancement by AGN outflow and CH3OH as a shock tracer, are correct.
    Table 3 and Section 4 rely on these assignments, some from prior work by the same group (Saito et al. 2022b).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Detailed analysis of multi-line molecular distributions in the Seyfert galaxy NGC 1068: Possible effect of the AGN outflow to the starburst ring." pith.science (2026). https://pith.science/paper/4NZ5KDLP

@misc{pith2026250615999,
  author       = {Pith},
  title        = {Pith review of: Detailed analysis of multi-line molecular distributions in the Seyfert galaxy NGC 1068: Possible effect of the AGN outflow to the starburst ring},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4NZ5KDLP}},
  note         = {Machine review of arXiv:2506.15999}
}
abstract

We apply principal component analysis (PCA) to the integrated intensity maps of 13 molecular lines of the nearby type-2 Seyfert galaxy NGC 1068 obtained by Atacama Large Millimeter/sub-millimeter Array (ALMA) to objectively visualize the features of its center, (1) within a radius of about 2 kpc ($\sim$ 27".5; hereafter the "overall region") and (2) the ring shaped starburst region between 750 pc ($\sim$ 10") and 2 kpc ($\sim$ 27".5) of the galaxy (hereafter the "SB ring region"). PCA is a powerful unsupervised machine learning technique that extracts key information through dimensionality reduction. The PCA results for the overall region have a possibility to reconstruct a map representing the approximate H$_2$ column density and difference of volume density and/or chemical composition between the circumnuclear disk (CND) and the starburst ring (SB ring). Additionally, the PCA results for the SB ring region have a possibility to reconstruct a map representing the approximate H$_2$ column density and distinction between starburst dominated region and shock dominated region. Furthermore, the PCA results for the SB ring region indicate a possible interaction between the Active Galactic Nucleus (AGN) outflow and gas in the SB ring. Although further investigation is required, we suggest that the AGN outflow interacts with gas in the SB ring, as this feature is consistent with the direction of the AGN outflow and is contributed by CN, C$_2$H and HCN, which are known to be enhanced by the AGN outflow. These results demonstrate that PCA can effectively extract features even for galaxies with complex structures, such as AGN + SB ring. This study also implies that PCA has the potential to uncover previously unrecognized phenomena by visualizing latent structures in multi-line data.

Figures

Figures reproduced from arXiv: 2506.15999 by the authors.

Figure 1
Figure 1. The 13CO(1–0) standardized integrated intensity maps of (a) overall region and (b) SB ring region of NGC 1068 used for PCA. The two dashed lines intersecting the AGN position (Roy et al. 1998) indicate the approximate outer boundaries of the ionized gas cones (Mingozzi et al. 2019). The central coordinates of this image are (α, δ)J2000 = (2h 42m40s .7132, -0◦ 0 ′ 47′′.655). The black circle marks the field of view (… view at source ↗
Figure 2
Figure 2. The standardized integrated intensity maps usd for the PCA [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. As same as Figure 2 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (11 more)
Figure 4
Figure 4. Figure 4: The standardized integrated intensity maps usd for the PCA [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: As same as Figure 4 [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: An image showing an example of PCA in three [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 7
Figure 7. Figure 7: PCA results in overall region: (a) contribution rate (CR) of each PC and cumulative contribution rate (CCR), (b) PC1 score map, [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]
Figure 8
Figure 8. Figure 8: PCAOA coefficients. Alt text: The three subfigures by PCA. From left to right, they represent PC1 coefficients, PC2 coefficients, and PC3 coefficients for overall region. tensity maps where the S/N of the CO (J=1-0) integrated inten￾sity map is ≤ 4. This step helps min…
Figure 9
Figure 9. Figure 9: These show PCA results in SB ring region: (a) contribution rate (CR) of each PC and cumulative contribution rate (CCR), (b) [PITH_FULL_IMAGE:figures/full_fig_p010_9.png]
Figure 10
Figure 10. Figure 10: PCASB coefficients. Alt text: The three subfigures by PCA. From left to right, they represent PC1 coefficients, PC2 coefficients, and PC3 coefficients for starburst ring region. 3 Results [PITH_FULL_IMAGE:figures/full_fig_p011_10.png]
Figure 11
Figure 11. Figure 11: The PCAOA results by normalized integrated intensity maps. Alt text: The results applied PCA for overall region by normalized integrated intensity maps labeled from (a) to (d). (a) shows contribution rate and cumulative contribution rate. From (b) to (d) show PC1 scor…
Figure 12
Figure 12. Figure 12: The PCAOA coefficients. Alt text: The three subfigures by PCA by normalized integrated intensity maps. From left to right, they represent PC1 coefficients, PC2 coefficients, and PC3 coefficients for overall region. within each fine structure line is nearly identical (…
Figure 13
Figure 13. Figure 13: The PCASB results by normalized integrated intensity maps. Alt text: The results applied PCA for starburst ring region by normalized integrated intensity maps labeled from (a) to (d). (a) shows contribution rate and cumulative contribution rate. From (b) to (d) show P…
Figure 14
Figure 14. Figure 14: The PCASB coefficients. Alt text: The three subfigures by PCA by normalized integrated intensity maps. From left to right, they represent PC1 coefficients, PC2 coefficients, and PC3 coefficients for starburst ring region. Acknowledgments This paper makes use of the fo…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

51 extracted references · 49 canonical work pages

  1. [1]

    and Neri, R., 2012, A & A, 537A, 44

    Henkel, C., Costagliola, F. and Neri, R., 2012, A & A, 537A, 44

    work page 2012

  2. [2]

    Aladro, R., et al., 2013, A & A, 549, 39

    work page 2013

  3. [3]

    P., Tollerud, E

    Aladro, R., et al., 2015, A & A, 579, 101 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33 Astropy Collaboration, Price-Whelan, A. M., Sip ˝ocz, B. M., et al. 2018, AJ, 156, 123

    work page 2015

  4. [4]

    Boger, G. I. and Sternberg, A., 2005, ApJ, 632, 1

    work page 2005

  5. [5]

    M., Kraemer, S

    Das, V ., Crenshaw, D. M., Kraemer, S. B., & Deo, R. P. 2006, AJ, 132, 620

    work page 2006

  6. [6]

    R., Rigopoulou, D

    Donnan, F. R., Rigopoulou, D. and García, I., 2024, MNRAS, 532, 75

    work page 2024

  7. [7]

    P., Schlemmer, S., Schilke, P., Stutzki, J., & Müller, H

    Endres, C. P., Schlemmer, S., Schilke, P., Stutzki, J., & Müller, H. S. 2016, JMoSp, 327, 95

    work page 2016

  8. [8]

    R., & Pineau des Forêts, G

    Flower, D. R., & Pineau des Forêts, G. 2012, MNRAS, 421, 2786

    work page 2012

Show all 51 references
  1. [9]

    and Bachiller, R., 1993, A & A, 276, 473

    Fuente, A., Martin-Pintado, J., Cernicharo, J. and Bachiller, R., 1993, A & A, 276, 473

    1993

  2. [10]

    2014, A &A, 567, A125 García-Burillo, S., et al., 2017, A & A, 608, 56

    Garay, G., Lizano, S., 1999, PASP, 111, 1049 García-Burillo, S., et al., 2010, A & A, 519, A2 García-Burillo, S., Combes, F., Usero, A., et al. 2014, A &A, 567, A125 García-Burillo, S., et al., 2017, A & A, 608, 56

    1999

  3. [11]

    Ginard, D., et al., 2015, A & A, 578, 49

    2015

  4. [12]

    2019, radio- astro-tools/spectralcube: Release v0.4.5, v0.4, Zenodo, doi: 10.5281/zenodo.35586145

    Ginsburg, A., Koch, E., Robitaille, T., et al. 2019, radio- astro-tools/spectralcube: Release v0.4.5, v0.4, Zenodo, doi: 10.5281/zenodo.35586145

    2019 doi

  5. [13]

    and Hirota, A., 2019, ApJ, 884, 2

    Harada, N., Sakamoto, K., Martín, S., Watanabe, Y ., Aladro, R. and Hirota, A., 2019, ApJ, 884, 2

    2019

  6. [14]

    Harada, N., et al., 2024, ApJ, 271, 38

    2024

  7. [15]

    and Wakelam, V ., 2010, ApJ, 721, 1570

    Harada, N., Herbst, E. and Wakelam, V ., 2010, ApJ, 721, 1570

    2010

  8. [16]

    and Herbst, E., 2013, ApJ, 765, 108

    Harada, N., Thompson, T, A. and Herbst, E., 2013, ApJ, 765, 108

    2013

  9. [17]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Natur, 585, 357

    2020

  10. [18]

    Y ., et al., 2023, A & A, 675, 151

    Huang, K. Y ., et al., 2023, A & A, 675, 151

    2023

  11. [19]

    Y ., Viti, S., 2023, Faraday Discussions, 245, 181

    Huang, K. Y ., Viti, S., 2023, Faraday Discussions, 245, 181

    2023

  12. [20]

    2009, AJ, 137, 3581

    Imanishi, M., Nakanishi, K., Tamura, Y ., & Peng, C.-H. 2009, AJ, 137, 3581

    2009

  13. [21]

    2013, Publications of the Astronomical Society of Japan, 65, 100

    Izumi, T., et al. 2013, Publications of the Astronomical Society of Japan, 65, 100

    2013

  14. [22]

    Menten, Karl M., 2017, A & A, 605L, 5

    2017

  15. [23]

    and Dalgarno, A., 1996, A & A, 306, 21 Martín, S., et al., 2021, A & A, 656, 46

    Lepp, S. and Dalgarno, A., 1996, A & A, 306, 21 Martín, S., et al., 2021, A & A, 656, 46

    1996

  16. [24]

    S., and Turner, J

    Meier, D. S., and Turner, J. L., 2005, ApJ, 618, 259

    2005

  17. [25]

    and Spaans, M., 2005, A & A, 436, 397

    Meijerink, R. and Spaans, M., 2005, A & A, 436, 397

    2005

  18. [26]

    Mingozzi, M., Cresci, G., Venturi, G., et al., 2019, A & A, 622, A146

    2019

  19. [27]

    P., Waters, B., Schiebel, D., Young, W., & Golap, K

    McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in ASP Conf. Ser. 376, Astronomical Data Analysis Software and Systems XVI, ed. R. A. Shaw, F. Hill, & D. J. Bell (San Francisco, CA: ASP),127

    2007

  20. [28]

    Nakajima, Y ., et al., 2024, ApJ, 974, 243

    2024

  21. [29]

    Nakajima, T., et al., 2015, PASJ, 67, 8

    2015

  22. [30]

    and Herbst, E., 2018, PASJ, 70, 7

    Nakajima, T., Takano, S., Kohno, K., Harada, N. and Herbst, E., 2018, PASJ, 70, 7

    2018

  23. [31]

    Nakajima, T., et al., 2023, ApJ, 955, 1

    2023

  24. [32]

    Okoda, Y ., Oya, Y ., Abe, S., Komaki, A., Watanabe, Y ., and Yamamoto, S., 2021, 923, 168

    2021

  25. [33]

    1901, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 2, 559

    Pearson, K. 1901, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 2, 559

    1901

  26. [34]

    2011, JMLR, 12, 2825

    Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, JMLR, 12, 2825

    2011

  27. [35]

    J., Tafalla, M., Gueth, F., Bachiller, R

    Rico-Villas, F., et al., 2021, MNRAS, 502, 3021 Rodríguez-Fernández, N. J., Tafalla, M., Gueth, F., Bachiller, R. 2010, A & A, 516, 98

    2021

  28. [36]

    L., Colbert, E

    Roy, A. L., Colbert, E. J. M., Wilson, A. S., & Ulvestad, J. S., 1998, ApJ, 504, 147

    1998

  29. [37]

    2017, ApJ, 834, 6

    Saito, T., et al. 2017, ApJ, 834, 6

    2017

  30. [38]

    2022, ApJL, 927, L32

    Saito, T., et al. 2022, ApJL, 927, L32

    2022

  31. [39]

    Saito, T., et al., 2022, ApJ, 935, 2

    2022

  32. [40]

    S., Wiedner, M

    Sakamoto, K., Aalto, S., Evans, A. S., Wiedner, M. C., and Wilner, D. J. 2010, ApJ, 725, L228 Sánchez-García, M., 2022, A & A, 660, 83

    2010

  33. [41]

    Scourfield, M., et al., 2020, MNRAS, 496, 4

    2020

  34. [42]

    Z., Matthews, K., Carico, D

    Scoville, N. Z., Matthews, K., Carico, D. P. and Sanders, D. B., 1988, ApJ, 327, 61

    1988

  35. [43]

    Takano, S., et al., 2014, PASJ, 66, 4

    2014

  36. [44]

    2018, ApJS, 236, 40

    Tanaka, K., Nagai, M., Kamegai, K., Iino, T., and Sakai, T. 2018, ApJS, 236, 40

    2018

  37. [45]

    Tanaka, K., et al., 2024, ApJ, 961, 18

    2024

  38. [46]

    Tosaki, T., et al., 2017, PASJ, 69, 18

    2017

  39. [47]

    Tsai, M., Hwang, C., Matsushita, S., Baker, A. J. and Espada, D., 2012, ApJ, 746, 129

    2012

  40. [48]

    2020, NatMe, 17, 261

    Virtanen, P., Gommers, R., Oliphant, T., et al. 2020, NatMe, 17, 261

    2020

  41. [49]

    Wang, J., Zhang, Zhi-Yu, Qiu, J., Shi, Y ., Zhang, Jiangshui and Fang, M., 2014, ApJ, 796, 57

    2014

  42. [50]

    and Kouchi, A., 2003, ApJ, 588, 121

    Watanabe, N., Shiraki, T. and Kouchi, A., 2003, ApJ, 588, 121

    2003

  43. [51]

    Xie, J., et al., 2023, ApJ, 949, 89

    2023

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.