HSC-M31 Monitoring Overview

• HSC-M31 monitoring is a multi-mode approach using Subaru’s Hyper Suprime-Cam to study microlensing events, nova demographics, and stellar halo substructures in Andromeda.
• High-cadence microlensing programs apply precise image-difference techniques and strict variability criteria to constrain the abundance of primordial black holes.
• NB515-enabled halo mapping discriminates M31 giants from foreground dwarfs, revealing key features like streams, overdensities, and radial metallicity gradients.

HSC-M31 Monitoring denotes a set of Subaru/Hyper Suprime-Cam observational programs targeting the Andromeda galaxy with distinct but partially overlapping objectives: high-cadence microlensing searches for primordial black holes, seasonal time-domain monitoring for classical and recurrent novae, and wide-field mapping of the stellar halo with the gravity-sensitive NB515 filter. In the published literature, these programs share a common instrumental basis—Subaru 8.2 m with Hyper Suprime-Cam—but differ substantially in cadence, filter choice, source-selection strategy, and statistical interpretation (Niikura et al., 2017). Taken together, they illustrate how a single wide-field imager can support pixel-lensing in crowded bulge fields, transient discovery across the disk, and resolved-star halo archaeology to projected radii of $\sim 120\,\mathrm{kpc}$ (Ogami et al., 2024).

1. Instrumental basis and observing modes

Hyper Suprime-Cam (HSC) is used in the cited M31 programs as a wide-field optical imager on Subaru 8.2 m. In the 2014 microlensing campaign, the telescope/instrument configuration was Subaru 8.2 m plus Hyper Suprime-Cam, with a field-of-view of 1.5° diameter, 104 CCDs in a single pointing covering the entire M31 disk, bulge and halo, a pixel scale of $0.168''\,\mathrm{pixel}^{-1}$, and Sloan-$r$ imaging (Niikura et al., 2017). The observation ran on one dark night, 2014 Nov 23, for a total on-sky duration of $\simeq 7\,\mathrm{hr}$ before M31 set below 30° elevation. It used 194 raw frames of 90 s integration each with 35 s overhead per frame, implying 2 min sampling; six worst-seeing frames with seeing $>1.2''$ were dropped, leaving 188 exposures. No dithering was applied in order to keep stars on the same CCD pixels (Niikura et al., 2017).

A later Subaru/HSC campaign, reanalyzed by Mróz and Udalski, used a closely related high-cadence mode but over multiple epochs: 2014-11-24 in Sloan $r$, 2017-09-20 and all eight nights in 2020 in Sloan $r2$, with 90 s integrations and 30 s readout, again giving 2 min per epoch. The number of usable exposures per night ranged from 44 to 214, and eight “good” nights were used in the candidate search, with total on-sky time of $\simeq 13.5\,\mathrm{hr}$ in 2014/2017 plus 25.8 hr in 2020 (Mróz et al., 31 Mar 2026). The summary of that reanalysis states that HSC has a 1.5 deg² field-of-view, $0.169''$ pixel scale, and 104 science CCDs, and that a single HSC pointing covers essentially the entire M31 disk and inner halo (Mróz et al., 31 Mar 2026). This suggests that the literature contains two slightly different field-of-view characterizations—1.5° diameter and 1.5 deg²—attached to closely related HSC M31 descriptions; both phrasings are part of the published summaries and should therefore be read in context.

For halo work, HSC was used in a different configuration centered on the custom narrow-band NB515 filter. NB515 is centered at $\lambda_c=515\,\mathrm{nm}$ with full-width at half-maximum $0.168''\,\mathrm{pixel}^{-1}$0, designed to straddle the MgH + Mgb absorption complex at $0.168''\,\mathrm{pixel}^{-1}$1–521 nm (Ogami et al., 2024). The halo program consisted of 33 pointings, 5 in 2015 and 28 in 2019, covering $0.168''\,\mathrm{pixel}^{-1}$2 in a roughly circular pattern out to $0.168''\,\mathrm{pixel}^{-1}$3 from M31’s center, corresponding to $0.168''\,\mathrm{pixel}^{-1}$4. Each NB515 field received four exposures of 240 s in most fields, while a few deeper fields received eight or twelve such exposures, under median seeing $0.168''\,\mathrm{pixel}^{-1}$5 with range $0.168''\,\mathrm{pixel}^{-1}$6–$0.168''\,\mathrm{pixel}^{-1}$7 (Ogami et al., 2024).

2. High-cadence microlensing architecture toward M31

The foundational HSC-M31 microlensing search was explicitly designed for dense-cadence detection of short-timescale events from primordial black holes in the halo regions of the Milky Way and M31 (Niikura et al., 2017). In dense M31 fields every pixel contains blended flux, so the analysis operates in the pixel-lensing regime: individual variable sources are detected via difference imaging rather than static source extraction (Niikura et al., 2017). The image-difference pipeline used a reference frame constructed as a coadd of the 10 best-seeing exposures, approximately $0.168''\,\mathrm{pixel}^{-1}$8, and target frames built in two stages: 63 coadds of 3 successive exposures, corresponding to 6 min cadence for initial candidate finding, and then all 188 single-exposure frames at 2 min cadence for light-curve extraction (Niikura et al., 2017).

Kernel-matched subtraction followed Alard & Lupton 1998 in the hscPipe implementation, performed on each “patch” of approximately $0.168''\,\mathrm{pixel}^{-1}$9 iso-latitude tessellation. Difference images were searched for positive or negative PSF-like residuals above a $r$0 threshold over Poisson noise (Niikura et al., 2017). On-site calibrations included bias subtraction and dome-flat division per CCD, plus high-order polynomial background subtraction—10th order over bulge CCDs and 6th order elsewhere—to remove sky and scattered light. The astrometric solution was tied to Pan-STARRS1 reference stars and updated every 11 frames, while photometric zero-points were also referenced to Pan-STARRS1. For PSF modeling, an initial bright-star catalog with $r$1 was used to derive per-CCD PSF models via PSFEx with second-order spatial variation (Niikura et al., 2017).

The light-curve model in difference-flux form was

$r$2

where

$r$3

and

$r$4

These definitions fix the event width in terms of a point-lens microlensing geometry (Niikura et al., 2017).

The later multi-night Subaru/HSC reanalysis retained the same basic physical target—sub-day microlensing in M31—but used an independent difference image analysis pipeline adapted from the OGLE DIA code. In that pipeline, each HSC CCD was divided into four quadrants, then into $r$5-pixel subfields; for each subfield the 15 best-seeing, low-background images from 2014-11-24 were coadded to form a deep reference image. Quadratic WCS mappings to Gaia DR3 achieved approximately 0.2–0.3 pixel RMS per axis, and each target-epoch frame was resampled to the reference grid and photometrically scaled to the same zeropoint. PSF photometry was then performed at each reference-catalog star position on the difference images, with fluxes zeropointed via Pan-STARRS1 DR2 to 0.01–0.03 mag accuracy (Mróz et al., 31 Mar 2026).

3. Event detection, source statistics, and formal sensitivity

In the 2014 microlensing analysis, a master catalog of 15,571 candidate variables was constructed by requiring each difference-image detection to satisfy PSF-shape cuts, a size criterion of 0.75–1.25 times the PSF FWHM, a roundness cut with axis ratio $r$6, and at least 2 detections in the 63 coadds (Niikura et al., 2017). At each of the 15,571 positions, PSF photometry was run on all 188 single-exposure difference images. A local background was estimated from the median in a surrounding $r$7 pixel stamp and subtracted before PSF flux measurement. The empirical noise $r$8 per epoch was derived by performing identical PSF photometry on 1000 random blank positions in each patch and measuring the RMS (Niikura et al., 2017).

The automated selection cuts were applied in sequence. The first, a “bump” requirement, demanded at least 3 consecutive $r$9 points each above $\simeq 7\,\mathrm{hr}$0, leaving 11,703 candidates. The second fitted $\simeq 7\,\mathrm{hr}$1 by minimizing

$\simeq 7\,\mathrm{hr}$2

and required $\simeq 7\,\mathrm{hr}$3 and $\simeq 7\,\mathrm{hr}$4, leaving 227 candidates. The third used the asymmetry metric

$\simeq 7\,\mathrm{hr}$5

with the requirement $\simeq 7\,\mathrm{hr}$6, leaving 146 candidates. A final single-clear-peak requirement, $\simeq 7\,\mathrm{hr}$7, reduced the sample to 66. Visual inspection then removed artifacts near bright stars, chip edges, and one moving object, leaving a single candidate with peak at $\simeq 7\,\mathrm{hr}$8 (Niikura et al., 2017).

The sensitivity analysis combined lensing geometry, halo modeling, efficiency simulations, and source-count estimation. The Einstein radius and timescale were written as

$\simeq 7\,\mathrm{hr}$9

and

$>1.2''$0

The optical depth for one M31 star was

$>1.2''$1

which is independent of $>1.2''$2 since $>1.2''$3; for an NFW halo of the Milky Way plus M31, with virial masses $>1.2''$4 and $>1.2''$5, one finds $>1.2''$6 (Niikura et al., 2017).

The differential event rate per star, using the Griest 1991 formalism, was given as

$>1.2''$7

with

$>1.2''$8

The integration yields $>1.2''$9 peaking at $r$0–1 hr for $r$1 (Niikura et al., 2017).

Detection efficiency $r$2 was estimated with Monte Carlo light-curve injection. For each $r$3, $r$4 simulated microlensing curves were generated, sampled at the real $r$5, had epoch-by-epoch Gaussian noise $r$6 added, and were passed through all selection cuts. The resulting efficiency was approximately 70–60% at $r$7–24 for $r$8–3 hr, falling to $r$9–30% at $r2$0–26; fake-star image injections with GalSim confirmed the Monte Carlo approach to within 10% (Niikura et al., 2017). The number of monitored source stars, $r2$1, was bounded conservatively at $r2$2 from peaks in the reference coadd and estimated more fully at $r2$3 down to $r2$4 by overlap with the HST/PHAT luminosity function and the transformation

$r2$5

The expected number of events was then

$r2$6

With one candidate, the Poisson 95% CL upper bound was $r2$7, leading to upper limits on $r2$8 as a function of mass (Niikura et al., 2017).

4. Primordial-black-hole constraints and their later reassessment

The original HSC-Andromeda study concluded that, given simultaneous monitoring of tens of millions of stars in M31, many microlensing events would be expected if light primordial black holes constituted a significant fraction of dark matter, but only a single candidate event was identified. This translated into the most stringent upper bounds on the abundance of primordial black holes in the mass range $r2$9 (Niikura et al., 2017). The key performance metrics reported for that campaign were a limiting magnitude per 90 s of $\simeq 13.5\,\mathrm{hr}$0 for a point source, typical image quality of $\simeq 13.5\,\mathrm{hr}$1 FWHM, 2 min sampling, optimal sensitivity to $\simeq 13.5\,\mathrm{hr}$2–3 hr, and an effective stellar sample reaching up to $\simeq 13.5\,\mathrm{hr}$3 sources depending on magnitude cut (Niikura et al., 2017). The final constraint was summarized as follows: for monochromatic primordial black holes, the mass fraction $\simeq 13.5\,\mathrm{hr}$4 must satisfy $\simeq 13.5\,\mathrm{hr}$5 in the range $\simeq 13.5\,\mathrm{hr}$6–$\simeq 13.5\,\mathrm{hr}$7, thereby closing a previously unconstrained lunar-mass window (Niikura et al., 2017).

The same analysis also reported two low-mass degradations of sensitivity. Finite source size, assuming $\simeq 13.5\,\mathrm{hr}$8 for most main-sequence targets, reduces magnification when $\simeq 13.5\,\mathrm{hr}$9, weakening bounds for $0.169''$0. Wave optics, comparing $0.169''$1 with $0.169''$2, further suppresses small-mass sensitivity at $0.169''$3 (Niikura et al., 2017).

A later controversy arose when a separate preprint by Sugiyama et al., summarized in the 2026 reanalysis, reported twelve candidates for short-timescale microlensing events and attributed them to a large population of planetary-mass primordial black holes. Mróz and Udalski reanalyzed the Subaru data using an independent difference image analysis photometric pipeline and concluded that all twelve candidates exhibited asymmetric light curves and/or variability on multiple nights of Subaru observations. Their classification assigned ten objects to RR Lyrae stars, one to an eclipsing binary, and one to an unclassified variable star; no compelling evidence for short-timescale microlensing events remained (Mróz et al., 31 Mar 2026).

The 2026 summary is explicit that the temporal and spatial distributions of the Subaru candidates were inconsistent with expectations for microlensing events, and that the results were in clear tension with previous searches toward the Magellanic Clouds, such as OGLE (Mróz et al., 31 Mar 2026). It also gives example limit translations: for $0.169''$4, Sugiyama et al. predicted $0.169''$5, so requiring $0.169''$6 at 95% CL gives $0.169''$7; at $0.169''$8, $0.169''$9 implies $\lambda_c=515\,\mathrm{nm}$0 (Mróz et al., 31 Mar 2026). Mróz and Udalski state that these Subaru-M31 limits are weaker than the 95% upper limits from OGLE-LMC/SMC, reported as $\lambda_c=515\,\mathrm{nm}$1 for $\lambda_c=515\,\mathrm{nm}$2, and also weaker than the earlier Niikura et al. bound (Mróz et al., 31 Mar 2026).

A common misconception is that very high cadence alone is sufficient for robust sub-day microlensing discovery in M31. The reanalysis argues otherwise: pulsating variables with sub-day periods, especially RR Lyrae, can mimic very short microlensing bumps if only a single night is used. In that account, strict symmetry cuts, robust multi-epoch checks, and external variable-star cross-matches are necessary components of the inference pipeline (Mróz et al., 31 Mar 2026).

5. Nova monitoring as a complementary HSC-M31 time-domain program

A distinct but closely related use of HSC-M31 monitoring concerns nova demographics. Shafter and Hornoch’s 2025 analysis of the recurrent nova population in M31 extends earlier work by Shafter et al. (2015) and is summarized with explicit recommendations for an HSC M31 monitoring program (Shafter et al., 19 Apr 2026). The historical nova discovery rate is reported as $\lambda_c=515\,\mathrm{nm}$3–10 yr$\lambda_c=515\,\mathrm{nm}$4 for 1909–1950, rising through the CCD era and reaching a plateau at $\lambda_c=515\,\mathrm{nm}$5 during 2005–2025; the extrapolated rate for 2021–2025 is $\lambda_c=515\,\mathrm{nm}$6 (Shafter et al., 19 Apr 2026). The Milky Way nova rate is given as $\lambda_c=515\,\mathrm{nm}$7–12 yr$\lambda_c=515\,\mathrm{nm}$8, approximately one third that of M31 (Shafter et al., 19 Apr 2026).

For confirmed recurrent novae, the summary quotes a sample size of 22 confirmed recurrent novae and 79 total recurrent-nova eruptions by mid 2025. Inter-eruption intervals range from 2.4 yr to 88.1 yr, with median observed recurrence time $\lambda_c=515\,\mathrm{nm}$9, 25th percentile $0.168''\,\mathrm{pixel}^{-1}$00, 75th percentile $0.168''\,\mathrm{pixel}^{-1}$01, and a fraction $0.168''\,\mathrm{pixel}^{-1}$02 with $0.168''\,\mathrm{pixel}^{-1}$03, shorter than U Sco (Shafter et al., 19 Apr 2026). The summary notes that no simple analytic fit is reported and that small-number statistics dominate.

Spatially, all novae, including classical and recurrent systems, are described as following the $0.168''\,\mathrm{pixel}^{-1}$04-band starlight profile of M31, encompassing the bulge and inner disk. The cumulative fractions within isophotal major-axis radius $0.168''\,\mathrm{pixel}^{-1}$05 are tabulated as follows (Shafter et al., 19 Apr 2026):

Radius $0.168''\,\mathrm{pixel}^{-1}$06 (′)	$0.168''\,\mathrm{pixel}^{-1}$07	$0.168''\,\mathrm{pixel}^{-1}$08
5′	0.30	0.28
10′	0.50	0.47
20′	0.75	0.78
30′	0.90	0.93

The Kolmogorov-Smirnov test gives $0.168''\,\mathrm{pixel}^{-1}$09, indicating no significant difference between classical and recurrent novae in this spatial comparison (Shafter et al., 19 Apr 2026).

The same summary gives a general M31 maximum-magnitude/rate-of-decline relation from Clark et al. (2024),

$0.168''\,\mathrm{pixel}^{-1}$10

with $0.168''\,\mathrm{pixel}^{-1}$11 and $0.168''\,\mathrm{pixel}^{-1}$12 for classical novae, and also states the empirical form

$0.168''\,\mathrm{pixel}^{-1}$13

again for classical novae only (Shafter et al., 19 Apr 2026). Recurrent novae occupy the “faint-and-fast” corner, with typical $0.168''\,\mathrm{pixel}^{-1}$14 and $0.168''\,\mathrm{pixel}^{-1}$15 days, whereas classical novae have $0.168''\,\mathrm{pixel}^{-1}$16 and $0.168''\,\mathrm{pixel}^{-1}$17–50 days (Shafter et al., 19 Apr 2026). Recurrent novae lie systematically $0.168''\,\mathrm{pixel}^{-1}$18 mag below the classical-nova MMRD at a given $0.168''\,\mathrm{pixel}^{-1}$19 (Shafter et al., 19 Apr 2026).

For an HSC observing season with limiting magnitude $0.168''\,\mathrm{pixel}^{-1}$20, the summary gives the distance modulus of M31 as $0.168''\,\mathrm{pixel}^{-1}$21, implying $0.168''\,\mathrm{pixel}^{-1}$22, and states that all known M31 novae peak at $0.168''\,\mathrm{pixel}^{-1}$23, or $0.168''\,\mathrm{pixel}^{-1}$24, hence 100% detectable above $0.168''\,\mathrm{pixel}^{-1}$25 (Shafter et al., 19 Apr 2026). Detection completeness depends mainly on decline rate: classical novae with $0.168''\,\mathrm{pixel}^{-1}$26–100 d are recovered at $0.168''\,\mathrm{pixel}^{-1}$27 even with a 7 d cadence, whereas recurrent novae with $0.168''\,\mathrm{pixel}^{-1}$28 d require cadence $0.168''\,\mathrm{pixel}^{-1}$29–3 d to ensure $0.168''\,\mathrm{pixel}^{-1}$30 detection of peak (Shafter et al., 19 Apr 2026). The expected HSC detections per year are quoted as $0.168''\,\mathrm{pixel}^{-1}$31 classical novae and $0.168''\,\mathrm{pixel}^{-1}$32 recurrent novae; for seasons of $0.168''\,\mathrm{pixel}^{-1}$33 months, these scale to $0.168''\,\mathrm{pixel}^{-1}$34 classical novae and $0.168''\,\mathrm{pixel}^{-1}$35–2 recurrent novae (Shafter et al., 19 Apr 2026).

This suggests that HSC-M31 monitoring is naturally bifurcated by cadence. For classical novae, $0.168''\,\mathrm{pixel}^{-1}$36 d is sufficient, whereas recurrent novae require $0.168''\,\mathrm{pixel}^{-1}$37–3 d and especially careful image registration below $0.168''\,\mathrm{pixel}^{-1}$38 to identify positional coincidences of repeated eruptions. The published recommendations accordingly emphasize difference-imaging in the crowded bulge, uniform imaging out to $0.168''\,\mathrm{pixel}^{-1}$39 major-axis radius, deeper bulge coverage within 5′, and multi-band follow-up in $0.168''\,\mathrm{pixel}^{-1}$40 and a bluer band (Shafter et al., 19 Apr 2026).

6. Halo-resolved monitoring with NB515 and implications for field design

The halo-oriented HSC-M31 program described by Komiyama and collaborators uses resolved-star selection rather than transient detection, but it is still a monitoring architecture in the broader sense of repeated, calibrated wide-field HSC observations of M31’s outskirts (Ogami et al., 2024). Its central methodological contribution is the use of NB515 to separate M31 halo giants from Milky Way foreground dwarfs with approximately 90% accuracy (Ogami et al., 2024).

The broadband and narrow-band catalogs were produced with hscPipe 6.7, applying bias/dark subtraction, flat-fielding, sky subtraction, cosmic-ray masking, and per-CCD astrometric and photometric calibration against Pan-STARRS1. Point-spread-function photometry yielded calibrated catalogs in $0.168''\,\mathrm{pixel}^{-1}$41, $0.168''\,\mathrm{pixel}^{-1}$42, and NB515, all corrected for Galactic extinction using the Schlegel et al. (1998) maps and a Fitzpatrick (1999) $0.168''\,\mathrm{pixel}^{-1}$43 law:

$0.168''\,\mathrm{pixel}^{-1}$44

$0.168''\,\mathrm{pixel}^{-1}$45

$0.168''\,\mathrm{pixel}^{-1}$46

The 50% completeness is 23.21 mag in NB515, while PAndAS $0.168''\,\mathrm{pixel}^{-1}$47 and $0.168''\,\mathrm{pixel}^{-1}$48 reach 50% completeness at $0.168''\,\mathrm{pixel}^{-1}$49 mag and $0.168''\,\mathrm{pixel}^{-1}$50 mag (Ogami et al., 2024).

Foreground rejection is based on the fact that dwarfs with $0.168''\,\mathrm{pixel}^{-1}$51–5 exhibit stronger MgH+Mgb absorption than giants with $0.168''\,\mathrm{pixel}^{-1}$52, making dwarfs redder in $0.168''\,\mathrm{pixel}^{-1}$53. The analysis defines

$0.168''\,\mathrm{pixel}^{-1}$54

where $0.168''\,\mathrm{pixel}^{-1}$55 is the dwarf ridge line. A dwarf-likelihood $0.168''\,\mathrm{pixel}^{-1}$56 is computed from a Gaussian model in $0.168''\,\mathrm{pixel}^{-1}$57, with $0.168''\,\mathrm{pixel}^{-1}$58, and latitude-dependent foreground density is modeled by

$0.168''\,\mathrm{pixel}^{-1}$59

leading to the M31 RGB membership probability

$0.168''\,\mathrm{pixel}^{-1}$60

Inside a Dartmouth-isochrone RGB box, the NB515-RGB sample is then defined by $0.168''\,\mathrm{pixel}^{-1}$61, and comparison with Keck/DEIMOS SPLASH classifications indicates $0.168''\,\mathrm{pixel}^{-1}$62 purity (Ogami et al., 2024).

The resulting NRGB map recovers the Giant Southern Stream, Eastern and Western shell fans, Streams C and D, and the North-Western stream, and identifies three new overdensities at $0.168''\,\mathrm{pixel}^{-1}$63: the Metal-Rich Cloud at $0.168''\,\mathrm{pixel}^{-1}$64 with $0.168''\,\mathrm{pixel}^{-1}$65 to $0.168''\,\mathrm{pixel}^{-1}$66, the SE Stream at $0.168''\,\mathrm{pixel}^{-1}$67 with $0.168''\,\mathrm{pixel}^{-1}$68 to $0.168''\,\mathrm{pixel}^{-1}$69, and the Metal-Poor Cloud at $0.168''\,\mathrm{pixel}^{-1}$70 with $0.168''\,\mathrm{pixel}^{-1}$71 to $0.168''\,\mathrm{pixel}^{-1}$72 (Ogami et al., 2024). Distances were derived by forward-modeling the color-magnitude diagram with a 4-parameter power-law-plus-plateau luminosity function,

$0.168''\,\mathrm{pixel}^{-1}$73

combined with a Gaussian metallicity model and Markov Chain Monte Carlo posterior estimation (Ogami et al., 2024). For the Giant Southern Stream, seven contiguous subfields yielded distances from $0.168''\,\mathrm{pixel}^{-1}$74 at $0.168''\,\mathrm{pixel}^{-1}$75 to $0.168''\,\mathrm{pixel}^{-1}$76 at $0.168''\,\mathrm{pixel}^{-1}$77, corresponding to a line-of-sight gradient of $0.168''\,\mathrm{pixel}^{-1}$78 (Ogami et al., 2024).

Globally, the ensemble metallicity distribution has $0.168''\,\mathrm{pixel}^{-1}$79, median $0.168''\,\mathrm{pixel}^{-1}$80, and $0.168''\,\mathrm{pixel}^{-1}$81, with a linear radial trend

$0.168''\,\mathrm{pixel}^{-1}$82

The $0.168''\,\mathrm{pixel}^{-1}$83-band surface-brightness profile obeys power laws $0.168''\,\mathrm{pixel}^{-1}$84 with $0.168''\,\mathrm{pixel}^{-1}$85 for the metal-poor population, $0.168''\,\mathrm{pixel}^{-1}$86 for the metal-rich population, and $0.168''\,\mathrm{pixel}^{-1}$87 for the full sample (Ogami et al., 2024). The abstract further emphasizes that the photometric metallicity distribution is spatially non-uniform with nonmonotonic trends with radius, suggesting insufficient time to dynamically homogenize the accreted populations (Ogami et al., 2024).

For future HSC-M31 monitoring, this halo study recommends a dual-filter strategy using a surface-gravity-sensitive narrow band such as NB515 plus two broad bands bracketing the red-giant branch, such as $0.168''\,\mathrm{pixel}^{-1}$88 and $0.168''\,\mathrm{pixel}^{-1}$89, reaching $0.168''\,\mathrm{pixel}^{-1}$90 at $0.168''\,\mathrm{pixel}^{-1}$91 in both NB515 and $0.168''\,\mathrm{pixel}^{-1}$92, corresponding to approximately $0.168''\,\mathrm{pixel}^{-1}$93 s NB515 exposures per field under $0.168''\,\mathrm{pixel}^{-1}$94 seeing (Ogami et al., 2024). A plausible implication is that halo-resolved monitoring and time-domain monitoring are not independent design problems: the same wide-field HSC footprint can be optimized simultaneously for substructure mapping, foreground control, and transient localization.

7. Methodological synthesis and recurring technical lessons

Across the microlensing, nova, and halo programs, several technical themes recur. First, precise external calibration is fundamental. The microlensing studies tied astrometry and photometric zero-points to Pan-STARRS1, with the later reanalysis also using Gaia DR3 for WCS solutions (Niikura et al., 2017, Mróz et al., 31 Mar 2026). The halo study likewise calibrated against Pan-STARRS1 and then imposed explicit extinction corrections (Ogami et al., 2024). Second, crowded-field inference in M31 depends heavily on image subtraction or probabilistic source classification rather than naive static-source photometry. The 2014 microlensing campaign used kernel-matched subtraction in hscPipe (Niikura et al., 2017); the 2026 reanalysis adapted the OGLE DIA code (Mróz et al., 31 Mar 2026); the nova-monitoring recommendations explicitly state that difference-imaging methods such as Alard & Lupton 1998 are strongly recommended in the bulge (Shafter et al., 19 Apr 2026).

Third, cadence must be matched to the phenomenon. The 2 min cadence of the original HSC microlensing search was chosen to resolve events as short as a few minutes and to optimize sensitivity to $0.168''\,\mathrm{pixel}^{-1}$95–3 hr (Niikura et al., 2017). For nova science, by contrast, $0.168''\,\mathrm{pixel}^{-1}$96 d is adequate for classical novae, while recurrent novae require $0.168''\,\mathrm{pixel}^{-1}$97–3 d to catch their fast evolution (Shafter et al., 19 Apr 2026). For halo mapping, monitoring is effectively replaced by depth, areal coverage, and filter complement, with four 240 s NB515 exposures per field under sub-arcsecond seeing (Ogami et al., 2024).

Fourth, variable-star rejection is a central issue wherever sub-day transients are sought. Mróz and Udalski argue that a symmetry test such as $0.168''\,\mathrm{pixel}^{-1}$98 was too loose to reject steep “sawtooth” pulsators and that multi-night tests must be enforced without subjective rescues for poor subtraction (Mróz et al., 31 Mar 2026). Their recommendations include at least 3–4 independent seasons, no correlated residuals on any other night at the same position above 2–3$0.168''\,\mathrm{pixel}^{-1}$99, stricter asymmetry cuts such as $r$00, multi-band follow-up in $r$01, and cross-matching against deep variable-star catalogs such as PS1 and Gaia variability flags (Mróz et al., 31 Mar 2026). This suggests that the main epistemic risk in HSC-M31 microlensing is not raw sensitivity but classification error in a crowded, variable-rich field.

Finally, HSC-M31 monitoring exemplifies a broader convergence of wide-field time-domain astronomy and resolved stellar-population studies. The same single-pointing coverage that supports pixel-lensing of tens of millions of stars in the bulge and disk (Niikura et al., 2017) also supports nova census work over the bulge and inner disk (Shafter et al., 19 Apr 2026), while the larger NB515 program extends the HSC-M31 framework into halo substructure mapping over $r$02 (Ogami et al., 2024). The published record therefore presents HSC-M31 monitoring not as one homogeneous survey, but as a family of technically linked observing modes whose scientific outputs range from dark-matter limits to recurrent-nova demographics and the chemodynamical structure of Andromeda’s stellar halo.