J0139+4328: Dark Galaxy Candidate Revisited

• Dark Galaxy Candidate J0139+4328 is an HI-rich system initially lacking detectable starlight until deep imaging revealed a faint, low-surface-brightness counterpart.
• Multi-instrument observations—including FAST surveys, deep optical imaging, and VLA mapping—establish its gas-dominated, rotating disk nature with a high HI-to-stellar mass ratio.
• This case underscores the need for combined HI and ultra-deep optical surveys to distinguish between truly star-free halos and faint stellar populations in dwarf galaxies.

A dark galaxy candidate is an HI-rich system with no detectable stellar component or emission at wavelengths tracing stars or dust. J0139+4328, also designated FAST J0139+4328 and AGC 249460 (in ALFALFA), has played a central role in the modern re-examination of the “dark galaxy” hypothesis through deep single-dish HI surveys, multi-wavelength optical/IR imaging, interferometric follow-up, and spectroscopic confirmation. Once considered a prototypical dark galaxy due to the absence of an optical counterpart in major surveys and the presence of a rotationally supported HI disk, J0139+4328 is now classified as an extremely gas-rich, low-surface-brightness (LSB) dwarf galaxy, whose stellar population lies below the detection thresholds of most all-sky imaging. This entry summarizes the discovery, multi-instrument characterization, controversy, and current astrophysical status of this object.

1. HI Detection and Initial Dark Galaxy Classification

The first identification of J0139+4328 as a dark galaxy candidate occurred through the Five-hundred-meter Aperture Spherical radio Telescope (FAST) HI survey (Xu et al., 2023). The source presented hallmark features of a disk galaxy:

• Double-horned HI profile and S-shaped position–velocity (PV) diagram structure
• Heliocentric velocity $V_\text{hel} = 2464.4 \pm 0.8$ km s$^{-1}$
• Integrated flux $\int S_\text{HI} dv = 0.424 \pm 0.020$ Jy km s$^{-1}$
• Fitted profile widths $W_{50} = 38.9 \pm 1.9$ km s$^{-1}$, $W_{20} = 51.1 \pm 2.0$ km s$^{-1}$

A Hubble-flow or Cosmicflows-3 distance of $D = 28.8 \pm 2.9$ Mpc implies $M_{\rm HI} = 8.3 \times 10^7\,M_\odot$. Tilted-ring modeling yielded $V_\text{rot}=26.9\pm9.4$ km s$^{-1}$ and $\sigma_v=7.8\pm1.4$ km s$^{-1}$, suggesting a rotating disk inclination $i = 27^\circ \pm 12^\circ$ over a spatial extent $R_\text{eff}=23.9\pm2.4$ kpc.

No optical, UV, NIR, or IR counterpart was visible in Pan-STARRS1, 2MASS, GALEX, or WISE/AKARI data at stringent limits; using $m_r > 21.7$ mag and $m_g > 22.0$ mag, the corresponding upper bound for the stellar mass was $M_\star < 6.9 \times 10^5 M_\odot$ and $M_B > -10$ mag, yielding $M_\text{HI}/L_B \gtrsim 59\,M_\odot/L_\odot$ (Xu et al., 2023). The dynamical mass, $M_\text{dyn} = (5.1 \pm 2.8) \times 10^9\,M_\odot$, implied a baryon-to-dynamical mass ratio $M_\text{dyn}/M_\text{bary} = 47 \pm 27$, with $M_\text{bary} \simeq 1.1 \times 10^8\,M_\odot$. The environment is isolated: no massive neighbors within $100$ kpc.

These properties led to the classification as the first bona-fide isolated dark dwarf galaxy in the nearby universe, consistent with cold dark matter models predicting HI-rich, star-free halos (sometimes called RELHICs) (Kwon et al., 4 Jun 2025, Xu et al., 2023).

2. Optical Non-Detection, Deep Imaging, and Stellar Counterpart Discovery

Initial confirmation efforts in surveys such as DESI Legacy Imaging ($r$-band limit $\sim 28.5$ mag arcsec$^{-2}$) and GALEX NUV (limiting $m_\mathrm{NUV} \sim 24$ AB mag) failed to reveal any stellar component (Kwon et al., 4 Jun 2025). However, follow-up with dedicated, long-exposure imaging using the 1.4 m Milanković and 0.6 m Nedeljković telescopes reached $B \sim 28$ mag arcsec$^{-2}$, four magnitudes deeper than Pan-STARRS1 (Mitrašinović et al., 31 Dec 2025):

• Detected a LSB optical counterpart with centroid ($$\alpha=01^{\mathrm h}39^{\mathrm m}29.6^{\mathrm s}$$, $\delta=+43\degree28'29.8''$)
• Central $B$-band surface brightness $\mu_{0,B} \approx 24.1$ mag arcsec$^{-2}$
• Half-light radius $r_\mathrm{eff} = 5.24'' = 0.73$ kpc
• Offset from HI peak much less than the FAST beam

Spectroscopic observations (6 m BTA, SCORPIO-2; 3650–7300 Å, $R\simeq1000$) detected H$\alpha$ emission at $V_\text{hel} = 2492 \pm 25$ km s$^{-1}$, closely matching the HI velocity and spatially coincident with the luminous centroid, thereby unambiguously confirming the HI–stellar association.

3. Stellar and Gas Content Quantification

Direct integration of the observed surface-brightness profile and color-dependent mass-to-light ratio yielded (Mitrašinović et al., 31 Dec 2025):

• Total $B$-band luminosity $L_B = (1.11 \pm 0.23) \times 10^7\, L_{\odot,B}$
• Mean color $(B-V) \sim 0.48 \pm 0.18$
• Stellar mass $M_* = (7.2 \pm 3.7) \times 10^6\, M_\odot$ (using Kroupa IMF, Into & Portinari relation)
• HI mass (from FAST detection at matched distance): $M_\text{HI} = 8.3 \pm 1.7 \times 10^7\, M_\odot$
• Gas-to-stellar mass ratio $M_\text{HI}/M_* = 11.5 \pm 6.4$

Comparisons with limits from Pan-STARRS1 ($M_* < 6.9\times10^5\,M_\odot$) demonstrate that previous surveys were simply not deep enough to detect the LSB disk; the true gas-to-star ratio was an order of magnitude lower than initially inferred.

4. Confirmatory Interferometric HI Imaging and Precise Galaxy Characterization

High-resolution HI mapping with the Karl G. Jansky Very Large Array (VLA), combining D and C configurations for $20.5$ hr on-source, provided a synthesized beam of $21.8'' \times 15.5''$, corresponding to $\sim 3$ kpc at $D \sim 31$ Mpc (Šiljeg et al., 18 Jan 2026). Key findings include:

• Refined HI centroid within $1''$ of the optical counterpart identified in Pan-STARRS1 stacking at $g\approx20.6$ mag
• HI line width $W_{50}=28.5 \pm 1.2$ km s$^{-1}$; $W_{20}=43.4\pm1.9$ km s$^{-1}$
• HI mass $M_\text{HI}=5.9\pm2.2 \times 10^7\,M_\odot$ at $D=31.0$ Mpc
• Stellar mass (from aperture photometry and mass-to-light ratios): $M_* = 3.2^{+6.4}_{-2.0}\times10^6\,M_\odot$
• $M_\text{HI}/M_* = 18$, confirming an extremely gas-dominated, LSB dwarf

Rotation velocities, after correcting for inclination and line-broadening, yield $V_\text{rot} = 12.3 \pm 2.8$ km s$^{-1}$, with substantial uncertainties due to warp, inclination, and marginal spatial resolution; nonetheless, morphology and kinematics are fully consistent with a very low-mass rotating disk.

5. Placement in Scaling Relations and Astrophysical Context

With the revised stellar and HI content, J0139+4328 occupies an extreme position among known galaxies:

• Within the 1$\sigma$ locus of HI-selected dwarfs on the $M_\text{HI}$–$M_*$ relation, but at low $M_*$ and high $M_\text{HI}/M_*$ (Šiljeg et al., 18 Jan 2026)
• On the baryonic Tully–Fisher relation, $$M_\text{bar} = M_* + 1.36 M_\text{HI} \approx 8.3 \times 10^7\,M_\odot$$, its $V_\text{rot}\sim12$ km s$^{-1}$ yields $M_\text{bar}$ substantially above the canonical relation, but well within the observed dispersion for low-mass, gas-rich dwarfs and UDG analogs where kinematic uncertainties are severe (Šiljeg et al., 18 Jan 2026, Mitrašinović et al., 31 Dec 2025)

No evidence persists for a star-free halo; instead, J0139+4328 exemplifies the ultra-faint, high-$M_\text{HI}/M_*$ regime accessible with deep HI and optical synergy.

6. Lessons for Dark Galaxy Searches and Future Prospects

The progression from dark galaxy candidate to confirmed LSB dwarf in J0139+4328 highlights methodological and astrophysical challenges:

• Single-dish HI centroids (e.g., from FAST, ALFALFA) have large positional uncertainties (20–30$''$), often exceeding the characteristic size of faint LSB disks.
• Reliable identification of low-mass, HI-dominated dwarfs requires both high-resolution interferometric HI imaging and deep, matched-depth optical (or near-IR) imaging to $\mu \gtrsim 28$ mag arcsec$^{-2}$ (Šiljeg et al., 18 Jan 2026).
• Multi-band photometry, robust mass-to-light calculations, and optical spectroscopy are essential to confirm physical association.
• The theoretical interpretation of isolated HI clouds with high gas fractions must avoid hasty “dark galaxy” claims unless stellar counterparts are excluded to extremely low surface-brightness limits.
• J0139+4328 conforms to the predicted properties of ΛCDM star-free halos only if prior optical limits are accepted; deeper imaging demonstrates these systems can harbor faint, yet detectable stellar populations.

A plausible implication is that the “dark galaxy” phenomenon in current extragalactic HI surveys reflects limitations in optical imaging depth and centroid assignment, not the presence of truly star-free HI halos (Mitrašinović et al., 31 Dec 2025, Šiljeg et al., 18 Jan 2026).

7. Summary Table: Key Observational Properties

Property	Value / Limit	Reference
Heliocentric velocity	$2464.4\pm0.8$ km s$^{-1}$	(Xu et al., 2023)
Distance	$28.8 \pm 2.9$ Mpc	(Xu et al., 2023)
HI mass	$(5.9\mbox{--}8.3)\times10^7\,M_\odot$	(Xu et al., 2023, Šiljeg et al., 18 Jan 2026)
Stellar mass	$3.2^{+6.4}_{-2.0}\mbox{--}7.2\pm3.7\times10^6\,M_\odot$	(Mitrašinović et al., 31 Dec 2025, Šiljeg et al., 18 Jan 2026)
Central $\mu_B$	$24.1$ mag arcsec$^{-2}$	(Mitrašinović et al., 31 Dec 2025)
$M_\text{HI}/M_*$	$11.5 \pm 6.4\mbox{--} 18$	(Mitrašinović et al., 31 Dec 2025, Šiljeg et al., 18 Jan 2026)
$V_\text{rot}$	$12.3 \pm 2.8$ km s$^{-1}$	(Šiljeg et al., 18 Jan 2026)
Baryonic Tully–Fisher	Consistent within scatter for gas-rich dwarfs	(Šiljeg et al., 18 Jan 2026)
Environment	Isolated; no massive neighbors within $100$ kpc	(Xu et al., 2023)

J0139+4328 remains among the most gas-rich dwarf galaxies known, and its reclassification underscores the importance of coordinated HI and optical observations in the extragalactic low-mass regime (Kwon et al., 4 Jun 2025, Mitrašinović et al., 31 Dec 2025, Šiljeg et al., 18 Jan 2026, Xu et al., 2023).