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 $^{-1}$0 km s$^{-1}$1 and $^{-1}$2 km s$^{-1}$3, suggesting a rotating disk inclination $^{-1}$4 over a spatial extent $^{-1}$5 kpc.

No optical, UV, NIR, or IR counterpart was visible in Pan-STARRS1, 2MASS, GALEX, or WISE/AKARI data at stringent limits; using $^{-1}$6 mag and $^{-1}$7 mag, the corresponding upper bound for the stellar mass was $^{-1}$8 and $^{-1}$9 mag, yielding $\int S_\text{HI} dv = 0.424 \pm 0.020$0 (Xu et al., 2023). The dynamical mass, $\int S_\text{HI} dv = 0.424 \pm 0.020$1, implied a baryon-to-dynamical mass ratio $\int S_\text{HI} dv = 0.424 \pm 0.020$2, with $\int S_\text{HI} dv = 0.424 \pm 0.020$3. The environment is isolated: no massive neighbors within $\int S_\text{HI} dv = 0.424 \pm 0.020$4 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 ($\int S_\text{HI} dv = 0.424 \pm 0.020$5-band limit $\int S_\text{HI} dv = 0.424 \pm 0.020$6 mag arcsec$\int S_\text{HI} dv = 0.424 \pm 0.020$7) and GALEX NUV (limiting $\int S_\text{HI} dv = 0.424 \pm 0.020$8 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 $\int S_\text{HI} dv = 0.424 \pm 0.020$9 mag arcsec$^{-1}$0, four magnitudes deeper than Pan-STARRS1 (Mitrašinović et al., 31 Dec 2025):

• Detected a LSB optical counterpart with centroid ($^{-1}$1, $^{-1}$2)
• Central $^{-1}$3-band surface brightness $^{-1}$4 mag arcsec$^{-1}$5
• Half-light radius $^{-1}$6 kpc
• Offset from HI peak much less than the FAST beam

Spectroscopic observations (6 m BTA, SCORPIO-2; 3650–7300 Å, $^{-1}$7) detected H$^{-1}$8 emission at $^{-1}$9 km s$W_{50} = 38.9 \pm 1.9$0, 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 $W_{50} = 38.9 \pm 1.9$1-band luminosity $W_{50} = 38.9 \pm 1.9$2
• Mean color $W_{50} = 38.9 \pm 1.9$3
• Stellar mass $W_{50} = 38.9 \pm 1.9$4 (using Kroupa IMF, Into & Portinari relation)
• HI mass (from FAST detection at matched distance): $W_{50} = 38.9 \pm 1.9$5
• Gas-to-stellar mass ratio $W_{50} = 38.9 \pm 1.9$6

Comparisons with limits from Pan-STARRS1 ($W_{50} = 38.9 \pm 1.9$7) 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 $W_{50} = 38.9 \pm 1.9$8 hr on-source, provided a synthesized beam of $W_{50} = 38.9 \pm 1.9$9, corresponding to $^{-1}$0 kpc at $^{-1}$1 Mpc (Šiljeg et al., 18 Jan 2026). Key findings include:

• Refined HI centroid within $^{-1}$2 of the optical counterpart identified in Pan-STARRS1 stacking at $^{-1}$3 mag
• HI line width $^{-1}$4 km s$^{-1}$5; $^{-1}$6 km s$^{-1}$7
• HI mass $^{-1}$8 at $^{-1}$9 Mpc
• Stellar mass (from aperture photometry and mass-to-light ratios): $W_{20} = 51.1 \pm 2.0$0
• $W_{20} = 51.1 \pm 2.0$1, confirming an extremely gas-dominated, LSB dwarf

Rotation velocities, after correcting for inclination and line-broadening, yield $W_{20} = 51.1 \pm 2.0$2 km s$W_{20} = 51.1 \pm 2.0$3, 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$W_{20} = 51.1 \pm 2.0$4 locus of HI-selected dwarfs on the $W_{20} = 51.1 \pm 2.0$5–$W_{20} = 51.1 \pm 2.0$6 relation, but at low $W_{20} = 51.1 \pm 2.0$7 and high $W_{20} = 51.1 \pm 2.0$8 (Šiljeg et al., 18 Jan 2026)
• On the baryonic Tully–Fisher relation, $W_{20} = 51.1 \pm 2.0$9, its $^{-1}$0 km s$^{-1}$1 yields $^{-1}$2 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-$^{-1}$3 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$^{-1}$4), 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 $^{-1}$5 mag arcsec$^{-1}$6 (Š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	$^{-1}$7 km s$^{-1}$8	(Xu et al., 2023)
Distance	$^{-1}$9 Mpc	(Xu et al., 2023)
HI mass	$D = 28.8 \pm 2.9$0	(Xu et al., 2023, Šiljeg et al., 18 Jan 2026)
Stellar mass	$D = 28.8 \pm 2.9$1	(Mitrašinović et al., 31 Dec 2025, Šiljeg et al., 18 Jan 2026)
Central $D = 28.8 \pm 2.9$2	$D = 28.8 \pm 2.9$3 mag arcsec$D = 28.8 \pm 2.9$4	(Mitrašinović et al., 31 Dec 2025)
$D = 28.8 \pm 2.9$5	$D = 28.8 \pm 2.9$6	(Mitrašinović et al., 31 Dec 2025, Šiljeg et al., 18 Jan 2026)
$D = 28.8 \pm 2.9$7	$D = 28.8 \pm 2.9$8 km s$D = 28.8 \pm 2.9$9	(Š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 $M_{\rm HI} = 8.3 \times 10^7\,M_\odot$0 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).