pith. machine review for the scientific record. sign in

arxiv: 2512.13954 · v1 · submitted 2025-12-15 · 🌌 astro-ph.CO · astro-ph.GA

Recognition: no theorem link

Catalogs of optically-selected clusters and photometric luminous red galaxies from the Hyper Suprime-Cam Subaru Strategic Program final year dataset

Authors on Pith no claims yet

Pith reviewed 2026-05-16 21:24 UTC · model grok-4.3

classification 🌌 astro-ph.CO astro-ph.GA
keywords galaxy clustersluminous red galaxiesphotometric redshiftsweak gravitational lensingcluster catalogslarge scale structureSubaru Hyper Suprime-Cam
0
0 comments X

The pith

The final Hyper Suprime-Cam dataset yields catalogs of more than 10,000 optically selected galaxy clusters and about 6 million photometric luminous red galaxies.

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

This paper releases new catalogs of galaxy clusters and luminous red galaxies selected from optical imaging in the Hyper Suprime-Cam Subaru Strategic Program. The cluster catalog covers approximately 1200 square degrees and includes over 10,000 clusters with richness above 15 extending to redshift 1.38, containing more than 600,000 probable red member galaxies. Photometric redshifts for clusters achieve scatter better than 0.01 across a wide range, while the LRG catalog reaches redshift 1.25 with photo-z scatter under 0.02 between 0.4 and 1.0. Stacked weak lensing signals from the clusters are measured out to redshift 1 to refine the mass-richness relation. These large samples support detailed investigations of galaxy evolution in clusters and large-scale structure in cosmology.

Core claim

We construct samples of optically-selected clusters and photometric luminous red galaxies from the Hyper Suprime-Cam Subaru Strategic Program final year dataset covering ∼1200 deg². The cluster catalogs extend out to the redshift of 1.38 and contain more than 10000 clusters with richness larger than 15, where the richness is defined to be a membership probability weighted number of galaxies above the stellar masses of approximately 10^{10.2}M_⊙. The total number of probable red cluster member galaxies in these clusters are more than 6×10^5. Photometric redshifts of the clusters are shown to be precise with the scatter better than ∼0.01 for a wide redshift range. We detect stacked weak lensed

What carries the argument

Richness, defined as the membership probability weighted number of galaxies above a stellar mass threshold of about 10^{10.2} solar masses, which selects cluster members and serves as a mass proxy.

Load-bearing premise

Color-based selection and photometric redshift estimates remain reliable with minimal bias or contamination across the full redshift range up to 1.38.

What would settle it

Spectroscopic observations of a subsample of high-redshift clusters revealing that the photometric redshift scatter exceeds 0.01 or that the stacked weak lensing signal is inconsistent with the reported mass-richness relation.

Figures

Figures reproduced from arXiv: 2512.13954 by Atsushi J. Nishizawa, Hung-Yu Jian, I-Non Chiu, Lihwai Lin, Masamune Oguri, Momoka Fujikawa, Naomi Ota, Nobuhiro Okabe, Rhythm Shimakawa, Tadayuki Kodama, Yen-Ting Lin, Yoshiki Toba.

Figure 1
Figure 1. Figure 1: Distribution of CAMIRA clusters on the sky. Filled red circles show positions of clusters in the Wide layer, while filled blue squares show posi￾tions of clusters in the Deep+UD layer. The bright star mask is applied. Alt text: The plot shows how CAMIRA clusters are distributed on the sky. choices of parameters used in the CAMIRA algorithm. As a byproduct, photometric LRG catalogs are also constructed taki… view at source ↗
Figure 2
Figure 2. Figure 2: Distribution of CAMIRA clusters in the redshift-richness plane. Contours indicates the density of clusters. Top and right panels show the histograms of the cluster redshift zcl and the richness Nˆmem, respectively. Alt text: In the main plot, the x axis shows the redshift from 0.1 to 1.38, and the y axis shows the richness from 15 to 200. Points show how CAMIRA clusters are distributed in this plane. 0.2 0… view at source ↗
Figure 3
Figure 3. Figure 3: Upper: The comparison between cluster photometric redshifts zcl and spectroscopic redshifts of BCGs zBCG,spec. Values of the bias δz, the scatter σz, and the outlier fraction of the redshift residual (zcl − zBCG,spec)/(1 + zBCG,spec) that are derived with 4σ clipping are also shown. Lower: The bias δz, the scatter σz of the redshift residual as a function of redshift. Alt text: the x and y axes show the re… view at source ↗
Figure 4
Figure 4. Figure 4: Examples of gri-composite color images of low- (left), mid- (middle), and high-redshift (right) CAMIRA clusters. The size of each image is ∼ 5. ′ 5. Alt text: concentrations of red galaxies are shown. 20 100 Nˆmem 101 102 103 104 L X E ( z ) − 2 [1042 erg/s] z < 0.2 0.2 ≤ z < 0.5 0.5 ≤ z [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: Histogram of positional offsets Roff , which are defined by physical transverse distances, between CAMIRA cluster centers and eRASS1 X￾ray cluster centers. The solid line shows the best-fitting model assuming a two-component model. Alt text: The x axis shows the positional offset from 0 to 1 megaparsec. 3.3 Stacked weak lensing analysis The mass-richness relation of HSC-SSP CAMIRA clusters has been constra… view at source ↗
Figure 7
Figure 7. Figure 7: Lensing profiles measured from the three-year HSC-SSP shear catalog. Different panels show results for different redshift bins. In each panel, we plot profiles for different richness bins with different symbols. Shaded regions indicate ranges of the radii that are excluded from the analysis (see text for more details). Solid lines indicate model predictions with the best-fit parameters listed in [PITH_FUL… view at source ↗
Figure 9
Figure 9. Figure 9: The color density map of photometric LRGs in the zLRG–M∗ plane, where zLRG and M∗ refer to photometric redshifts and stellar masses of LRGs, respectively. The density is calculated by a kernel density estima￾tion using Gaussian kernels. Due to the large number of photometric LRGs for Wide, here we show the density map for Deep+UD as they are similar between Wide and Deep+UD. Top and right panels show the h… view at source ↗
Figure 10
Figure 10. Figure 10: Similar to [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Example of lensing profiles Σ(R) (upper) calculated by equa￾tion (A13) and ∆Σ(r) (lower) calculated by equation (A14). We assume the halo mass M200m = 1 × 1014h −1M⊙, the redshift z = 0.5, the mis￾centering parameters fcen = 0.6 and σoff = 0.3h −1Mpc, and the cosmo￾logical parameters shown in Sec. 3.3. Note that we show lensing profiles in the physical length scale. Dashed, dash-dotted, and dotted lines s… view at source ↗
Figure 11
Figure 11. Figure 11: Appendix. References Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2025, Phys. Rev. D, 112, 083535 Aguado, D. S., Ahumada, R., Almeida, A., et al. 2019, ApJS, 240, 23 Aihara, H., Armstrong, R., Bickerton, S., et al. 2018a, PASJ, 70, S8 Aihara, H., Arimoto, N., Armstrong, R., et al. 2018b, PASJ, 70, S4 Aihara, H., AlSayyad, Y., Ando, M., et al. 2019, PASJ, 71, 114 —. 2022, PASJ, 74, 247 Alarcon, A., Ga… view at source ↗
read the original abstract

We construct samples of optically-selected clusters and photometric luminous red galaxies (LRGs) from the Hyper Suprime-Cam Subaru Strategic Program final year dataset covering $\sim 1200$~deg$^2$. The cluster catalogs extend out to the redshift of $1.38$ and contain more than 10000 clusters with richness larger than $15$, where the richness is defined to be a membership probability weighted number of galaxies above the stellar masses of approximately $10^{10.2}M_\odot$. The total number of probable red cluster member galaxies in these clusters are more than $6\times 10^5$. Photometric redshifts of the clusters are shown to be precise with the scatter better than $\sim 0.01$ for a wide redshift range. We detect stacked weak lensing signals of clusters out to the redshift of $1$, and use them to update constraints on the mass-richness relation. Our catalog of about 6 million photometric LRGs extend out to the redshift of $1.25$, and have the scatter of the photometric redshift better than $\sim 0.02$ for the redshift range between $0.4$ and $1.0$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript constructs catalogs of optically-selected galaxy clusters and photometric luminous red galaxies from the Hyper Suprime-Cam Subaru Strategic Program final year dataset covering ~1200 deg². The cluster catalog contains more than 10,000 clusters with richness >15 (defined as membership-probability-weighted galaxy counts above ~10^10.2 M_⊙) extending to z=1.38, reports photo-z scatter better than ~0.01, detects stacked weak lensing signals to z=1 to update the mass-richness relation, and provides a catalog of ~6 million LRGs to z=1.25 with photo-z scatter better than ~0.02 for 0.4<z<1.0.

Significance. If the photo-z and membership calibrations hold, the catalogs represent a substantial increase in high-redshift optically selected clusters and LRGs over a wide area, enabling improved cosmological constraints and galaxy evolution studies; the weak-lensing mass calibration update is a clear strength for the cluster sample.

major comments (1)
  1. [photometric redshift validation and cluster membership section] The section describing cluster photo-z performance and membership probabilities: the reported scatter better than ~0.01 up to z=1.38 and the resulting richness values (>10k clusters, >6e5 members) are derived from color-based selection and photo-z weighting, but the manuscript shows only the scatter on the full sample without an independent spectroscopic validation set that is both (a) separate from the training set and (b) representative of the cluster environment at z≈1.2–1.38. This directly affects the reliability of the high-redshift richness and stacked lensing mass-richness relation.
minor comments (2)
  1. Clarify in the text whether the ~1200 deg² coverage represents the full final-year footprint after masking or a specific subset, and provide the exact effective area used for the cluster and LRG samples.
  2. [weak lensing section] In the weak-lensing results, specify the exact redshift bins used for the stacked signals and report the signal-to-noise or detection significance per bin to allow direct assessment of the mass-richness update.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough review and valuable comments on our manuscript. We address the major comment point by point below, clarifying our validation approach and indicating the revisions we will make.

read point-by-point responses
  1. Referee: [photometric redshift validation and cluster membership section] The section describing cluster photo-z performance and membership probabilities: the reported scatter better than ~0.01 up to z=1.38 and the resulting richness values (>10k clusters, >6e5 members) are derived from color-based selection and photo-z weighting, but the manuscript shows only the scatter on the full sample without an independent spectroscopic validation set that is both (a) separate from the training set and (b) representative of the cluster environment at z≈1.2–1.38. This directly affects the reliability of the high-redshift richness and stacked lensing mass-richness relation.

    Authors: We appreciate the referee pointing out the need for clearer validation details. The photo-z scatter is indeed reported for the full photometric sample, but this includes a cross-validation where the spectroscopic training set is split, with the validation portion held out and not used in training. This held-out set is representative of the general galaxy population. For cluster-specific performance, we have performed checks using spectroscopic redshifts of confirmed cluster members from overlapping surveys at z<1, showing consistent scatter. At z=1.2-1.38, while a large independent spectroscopic sample of cluster galaxies is not available (as noted in the literature due to observational challenges), the photo-z method relies on colors that are calibrated across environments, and the richness definition uses probability weighting to mitigate outliers. The stacked weak lensing detection to z=1 provides supporting evidence for the catalog reliability. In the revised manuscript, we will add a dedicated subsection detailing the training/validation split, present any available high-z spec comparisons, and discuss the potential impact on high-redshift results. This will strengthen the presentation without altering the core results. revision: partial

Circularity Check

0 steps flagged

Observational catalog paper with no load-bearing circular derivations

full rationale

This is a data release and catalog construction paper. Cluster richness is defined directly as a membership-probability-weighted count of galaxies above a stellar-mass threshold, and photometric redshifts are computed from standard color-based methods. No equations or results are presented as 'predictions' that reduce to the same fitted inputs by construction. The stacked weak-lensing mass-richness update uses external shear measurements and does not loop back to the selection function. Self-citations exist for prior HSC work but are not load-bearing for the central catalog numbers. Overall circularity is minimal (score 1).

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

As a data catalog paper the central claim rests on observational processing pipelines rather than new theoretical constructs; no free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption Standard assumptions in photometric redshift estimation and optical cluster finding algorithms hold across the surveyed redshift range
    The reported photo-z scatters and richness definition rely on established color-based methods without new validation shown in the abstract.

pith-pipeline@v0.9.0 · 5573 in / 1231 out tokens · 54402 ms · 2026-05-16T21:24:01.337478+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages · 2 internal anchors

  1. [1]

    Publications of the Astronomical Society of Japan, 2024,00, 1–13 doi: 10.1093/pasj/xxx000 Advance access publication date: Catalogs of optically-selected clusters and photometric lu- minous red galaxies from the Hyper Suprime-Cam Subaru Strategic Program final year dataset Masamune OGURI, 1,2,∗ Yen-Ting LIN, 3 Nobuhiro OKABE, 4,5,6 Naomi OTA, 7 I-Non CHIU...

  2. [2]

    We detect stacked weak lensing signals of clusters out to the redshift of1, and use them to update constraints on the mass-richness relation

    Photometric redshifts of the clusters are shown to be precise with the scatter better than∼0.01for a wide redshift range. We detect stacked weak lensing signals of clusters out to the redshift of1, and use them to update constraints on the mass-richness relation. Our catalog of about 6 million photometric LRGs extend out to the redshift of1.25, and have t...

  3. [3]

    Catalogs of optically-selected clusters and photometric luminous red galaxies from the Hyper Suprime-Cam Subaru Strategic Program final year dataset

    Astronomical Society of Japan. arXiv:2512.13954v1 [astro-ph.CO] 15 Dec 2025 2Publications of the Astronomical Society of Japan(2024), Vol. 00, No. 0 ing weak gravitational lensing (see Oguri & Miyazaki 2025, for a review). Clusters of galaxies are efficiently identified from concentra- tions of red galaxies (e.g., Gladders & Yee 2000; Koester et al. 2007;...

  4. [4]

    (2018a) construct an optically-selected cluster catalog out to the redshift of1.1from the first-year (S16A) HSC-SSP dataset covering∼232deg

    that makes use of the stellar pop- ulation synthesis model of Bruzual & Charlot (2003) to model the red-sequence, Oguri et al. (2018a) construct an optically-selected cluster catalog out to the redshift of1.1from the first-year (S16A) HSC-SSP dataset covering∼232deg

  5. [5]

    2019; Chiu et al

    The mass-richness relation of CAMIRA clusters are calibrated by weak gravitational lens- ing (Murata et al. 2019; Chiu et al. 2020b) and clustering (Chiu et al. 2020a). Follow-up observations with X-ray (Ota et al

  6. [6]

    2021; Kitayama et al

    and the Sunyaev-Zel’dovich effect (Okabe et al. 2021; Kitayama et al

  7. [7]

    As the HSC-SSP survey progresses, the CAMIRA cluster catalog is continuously updated (e.g., Willis et al

    confirm that the hot gas is indeed associated with rich CAMIRA clusters. As the HSC-SSP survey progresses, the CAMIRA cluster catalog is continuously updated (e.g., Willis et al. 2021; Oguri et al. 2021; Klein et al. 2022; Hashiguchi et al. 2023; Toba et al. 2024; Chen et al. 2024). Wen & Han (2021) construct an independent cluster catalog from the HSC-SS...

  8. [8]

    We study the properties of our cluster catalogs, including the mis-centering and the mass- richness relation. In addition to the cluster catalogs, we construct catalogs of photometric luminous red galaxies (LRGs), which are byproducts of the CAMIRA algorithm to fit galaxy colors and red- shifts. The photometric LRG catalog from the first-year (S16A) HSC-S...

  9. [9]

    2 Data and method 2.1 HSC-SSP final year dataset The HSC-SSP is a wide-field optical imaging survey with the Hyper Suprime-Cam (HSC; Miyazaki et al

    Catalogs will be made available at https://github.com/oguri/cluster_catalogs. 2 Data and method 2.1 HSC-SSP final year dataset The HSC-SSP is a wide-field optical imaging survey with the Hyper Suprime-Cam (HSC; Miyazaki et al. 2018; Komiyama et al. 2018; Furusawa et al

  10. [10]

    The HSC-SSP observations are conduced in 2014–2021 and the data are reduced with hscPipe (Bosch et al

    over more than 1,000 deg2 with a limiting magnitude ofi lim ∼26at5σfor point sources with a2 ′′ diameter aperture. The HSC-SSP observations are conduced in 2014–2021 and the data are reduced with hscPipe (Bosch et al. 2018). The HSC-SSP data are made publicly available through a series of data releases (Aihara et al. 2018a, 2019, 2022). The HSC-SSP final ...

  11. [11]

    Additionally we remove galaxies withz-bandcmodelmagnitude error larger than0.1

    brighter thanz-band magnitudes of24. Additionally we remove galaxies withz-bandcmodelmagnitude error larger than0.1. As in Oguri et al. (2018a), we usecmodelmagni- tudes for magnitudes of galaxies, while the point spread function (PSF)-matched aperture photometry with the target PSF size of 1.′′3and the aperture of1. ′′5in diameter to avoid the deblending...

  12. [12]

    2.1), using a matching separation tolerance of1 ′′

    We cross match these spectroscopic galaxies with photometric galaxies in the HSC-SSP with the bright star mask (see Sec. 2.1), using a matching separation tolerance of1 ′′. We remove any spec- troscopic galaxies with measurement error larger than0.01(1 + zspec). The numbers of spectroscopic galaxies that are cross- matched are 2925976 for Wide and 172930 ...

  13. [13]

    The number of the richest clusters declines with increasing redshift, reflecting a rapid evo- lution of the number density of massive dark matter halos

    It is found that many clusters of galaxies are identified out toz∼1. The number of the richest clusters declines with increasing redshift, reflecting a rapid evo- lution of the number density of massive dark matter halos. We note that the selection of cluster member galaxies used for com- puting the richness, which correspond to the stellar mass limit ofM...

  14. [14]

    3.2 Comparison with an X-ray cluster catalog We cross match the CAMIRA cluster catalog with SRG/eROSITA All-Sky Survey first cluster catalog (hereafter eRASS1; Bulbul et al. 2024). The matching condition is the redshift difference of |∆z| ≤0.05and the physical transverse distance ofR≤1h −1Mpc, where the transverse distance is computed assuming the redshif...

  15. [15]

    2018a; Okabe et al

    We find a positive correlation between these two quantities, which is consistent with previous analysis of X-ray properties of CAMIRA clusters (e.g., Oguri et al. 2018a; Okabe et al. 2019; Willis et al. 2021; Ota et al. 2023). The relatively large scatter mainly reflects the large scatter of the richness-mass relation. A detailed analysis of the richness-...

  16. [16]

    We derive differential surface mass density profiles∆Σ(R) for subsamples of CAMIRA clusters following Medezinski et al

    that covers ∼430deg 2 of the HSC-SSP footprint with the raw number density of22.9arcmin −2. We derive differential surface mass density profiles∆Σ(R) for subsamples of CAMIRA clusters following Medezinski et al. (2018). Errors on the measurements include intrinsic ellipticities of galaxies. We use theP(z)cut method, for which an integrated probability of ...

  17. [17]

    In Table 2, we summarize the median and 68% credible interval of each model parameter after marginalizing over the other parameters

    From the MCMC analysis, we find that our model fits the ob- served lensing profiles reasonably well with the minimumχ 2 of χmin/dof = 182.4/182. In Table 2, we summarize the median and 68% credible interval of each model parameter after marginalizing over the other parameters. We find that model parameters of the richness-mass relation are broadly consist...

  18. [18]

    have to be taken into account as emphasized in Murata et al. (2019). We compare observed lensing profiles with our best-fit model in Figure

  19. [19]

    Given the richness-mass relation constrained from observed lensing profiles, we can derive the mass-richness relation using the Bayes theorem (Murata et al. 2019). P(lnM|N, z) = P(lnN|M, z)P(lnM|z)R dlnM P(lnN|M, z)P(lnM|z) ,(8) whereP(lnN|M, z)is defined in equation (2) and we use the halo mass function Publications of the Astronomical Society of Japan(2...

  20. [20]

    Alt text: In each panel, the x axis shows the distance from the cluster center from 0.1 to 10 megaparsec

    In sub-panels, we also show residuals (i.e., observed minus model predicted values) normalized by errors. Alt text: In each panel, the x axis shows the distance from the cluster center from 0.1 to 10 megaparsec. The y axis shows the differential surface mass distribution. 8Publications of the Astronomical Society of Japan(2024), Vol. 00, No. 0 15 20 30 40...

  21. [21]

    2018; Sunayama et al

    However, in order to place accurate constraints on cosmological parameters, the projection effect arising from the asphericity of clusters as well as neighboring galaxies have to be carefully taken into account (e.g., Osato et al. 2018; Sunayama et al. 2020; Wu et al. 2022; Lee et al. 2025). For instance, Sunayama et al. (2020) show that optically-selecte...

  22. [22]

    photometric LRG

    the method described in Sec. 2.3.2 The photometric LRG catalogs contain 6705294 galaxies from the Wide layer after applying for the bright star mask, and 119955 galaxies from the Deep+UD layer after applying for the bright star mask. We note that the photomet- ric LRG catalogs are constructed independently of the CAMIRA cluster catalogs presented in Sec. ...

  23. [23]

    This work was supported by JSPS KAKENHI Grant Numbers JP25H00662, JP25H00672, JP25H01551, JP25K01026, JP25K01044, JP24H00002, JP24K00684, JP23H00108, JP23K22537, JP22K21349

    Acknowledgments We thank an anonymous referee for useful suggestions. This work was supported by JSPS KAKENHI Grant Numbers JP25H00662, JP25H00672, JP25H01551, JP25K01026, JP25K01044, JP24H00002, JP24K00684, JP23H00108, JP23K22537, JP22K21349. YTL acknowledge support from National Science and Technology Council of Taiwan under grant NSTC 114-2112-M-001-03...

  24. [24]

    (2019), here we take the so-called halo model approach to compute cluster lensing profiles, largely following Oguri & Hamana (2011)

    has been used for the previous analysis of the mass-richness relation in Murata et al. (2019), here we take the so-called halo model approach to compute cluster lensing profiles, largely following Oguri & Hamana (2011). In this approach, contributions from one-halo (1h) and two-halo (2h) terms are considered separately. For the one-halo term, we adopt the...

  25. [25]

    The parameterρ s andr s are connected with the halo mass and the concentration parameter ignoring the truncation (see Oguri & Hamana 2011)

    This density profile is smoothly truncated atr≃r t. The parameterρ s andr s are connected with the halo mass and the concentration parameter ignoring the truncation (see Oguri & Hamana 2011). We introduce the dimensionless truncation radius defined byτ 200m =r t/r200m, and assumeτ 200m = 2.5through- out the paper based on the comparison with ray-tracing s...

  26. [26]

    (2020) and also used inDarkEmulator(Nishimichi et al

    implemented by Fang et al. (2020) and also used inDarkEmulator(Nishimichi et al

  27. [27]

    for fast calculations of the inverse Fourier transform. Following Oguri & Hamana (2011) and Oguri & Takada (2011), the Fourier transform of the two-halo term is computed as ˜Σ2h(k) = ¯ρm0bh(M)Pm(k/(1 +z);z),(A12) 10□2 10□1 100 101 102 R [h□1Mpc] 1011 1012 1013 1014 1015 Σ(R) [ h M⊙Mpc□2] total 1h center 1h mis-center 2h 10□2 10□1 100 101 102 R [h□1Mpc] 10...

  28. [28]

    Data Release 1 of the Dark Energy Spectroscopic Instrument

    Appendix. References Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2025, Phys. Rev. D, 112, 083535 Aguado, D. S., Ahumada, R., Almeida, A., et al. 2019, ApJS, 240, 23 Aihara, H., Armstrong, R., Bickerton, S., et al. 2018a, PASJ, 70, S8 Aihara, H., Arimoto, N., Armstrong, R., et al. 2018b, PASJ, 70, S4 Aihara, H., AlSayyad, Y ., Ando, M., et al. 2019, ...