🔭 astrophysics

Mapping the Milky Way's Rotation Curve from Archival 21-cm H I Observations

This undergraduate-level study utilizes archival HI4PI 21-cm data to construct a Galactic rotation curve via the tangent-point method, demonstrating both the utility of 21-cm kinematics and the critical importance of quality control diagnostics when interpreting derived velocities and handling detection limitations.

Original authors: Shibam Ghosh, Joshal Bhokare

Published 2026-10-26

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Shibam Ghosh, Joshal Bhokare

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

To understand our place in the universe, we must first understand the shape of the galaxy we call home. The Milky Way is a vast, spinning disk of stars, gas, and dust, but we live inside it, like a person trying to map a forest while standing in the middle of a tree. We cannot step back to take a photograph of the whole structure. Instead, astronomers must piece together the galaxy's three-dimensional shape and its gravitational pull by observing how things move from our specific vantage point. A key tool for this is neutral atomic hydrogen, a common gas that fills the space between stars. Unlike visible light, which gets blocked by cosmic dust, radio waves emitted by this gas can pass through the dark clouds, revealing the hidden motion of the galaxy's inner regions. By measuring how fast this gas moves toward or away from us, scientists can calculate how the galaxy rotates and, from that, infer how much mass is hidden within.

A recent study by Shibam Ghosh and Joshal Bhokare takes this classic approach and applies it to a massive, publicly available archive of radio data. Their goal was to reconstruct the rotation curve of the inner Milky Way—the relationship between a star's distance from the center and its orbital speed—using a method known as the tangent-point technique. This method relies on a specific geometric trick: when we look toward the inner galaxy, the gas moving fastest in our line of sight is located at the closest point to the galactic center that our view can reach. By measuring the speed of this fastest gas at different angles, researchers can map out the rotation of the galaxy's inner disk. The researchers used data from the HI4PI survey, which covers the entire sky, to build a detailed map of gas speeds and then tested how well a specific set of rules could extract the necessary measurements.

The team began by creating a two-dimensional map that plots the direction we look against the speed of the gas we see. In this map, bright ridges of emission show where large amounts of gas are moving together. To find the rotation speed, the researchers needed to identify the very edge of this emission—the point where the gas speed stops increasing and the signal fades away. They set a strict rule: they would only count a speed as valid if the gas signal was bright enough to be clearly distinguished from the background noise, specifically requiring a brightness level of 20 units on their temperature scale. They also averaged the data over a wide strip of the sky to smooth out small fluctuations. This approach worked well for many directions, successfully tracing the outer boundary of the gas emission and converting those speeds into a rotation curve that stretches from about 1.5 to 7.4 kiloparsecs from the center.

The resulting rotation curve tells a clear story about the inner galaxy. Starting from the innermost sampled region, the orbital speed rises quickly from about 95 kilometers per second. As the distance from the center increases, the speed continues to climb, eventually settling into a range between 230 and 255 kilometers per second for much of the outer part of the measured interval. This high, steady speed is significant because it suggests that the mass of the galaxy is not concentrated in a single central point. If all the mass were packed into the center, the orbital speed would drop off sharply as one moved outward, much like how planets farther from the sun move more slowly. Instead, the gas keeps moving fast, implying that there is a vast amount of mass distributed throughout the disk and beyond, a finding that aligns with the idea of an extended gravitational field.

However, the study is perhaps most valuable for what it reveals about the difficulties of the process, rather than just the final numbers. When the researchers applied their strict brightness rule to the data at high angles, looking toward the outer edges of the inner galaxy, the method failed. In several specific directions, the gas signal was simply too faint to meet the 20-unit threshold, meaning no valid speed could be recorded. Because the data stopped before reaching the angle where the rotation should theoretically slow down to zero, the team could not perform a crucial check to see if their measurements were perfectly accurate. Consequently, they could not apply a correction to fix any small errors in their measurements. The final result shows a rotation curve that looks physically plausible and matches the general behavior expected of a spiral galaxy, but it comes with a clear warning label: the data is incomplete at the edges, and the measurements are sensitive to the specific rules used to define the edge of the signal.

This work serves as a transparent demonstration of how modern astronomy is done. It shows that even with powerful data and clear physical principles, the path to a precise measurement is filled with choices and limitations. The researchers did not hide the gaps in their data or the moments where their method stumbled; instead, they highlighted them to show that a scientific result is only as strong as the diagnostics that support it. By mapping the rotation of the Milky Way and simultaneously exposing the boundaries of their own technique, Ghosh and Bhokare provide a realistic picture of how we learn about the invisible forces that hold our galaxy together. The study confirms that the inner galaxy rotates with a speed that requires a massive, extended distribution of matter, while also reminding us that every measurement is a negotiation between the data we have and the rules we use to interpret it.

Technical Summary: Mapping the Milky Way's Rotation Curve from Archival 21-cm H I Observations

Problem Statement
The Milky Way cannot be observed externally; its three-dimensional structure and gravitational field must be reconstructed from internal observations of angular positions, spectra, and motions. Optical tracers are often obscured by interstellar dust in the Galactic plane. This work addresses the challenge of deriving the inner Galaxy's rotation curve using neutral atomic hydrogen (H I), which penetrates dust via its 21-cm hyperfine line. The specific problem is to construct a reproducible workflow that transforms archival radio spectral data into a kinematic rotation curve while explicitly identifying and reporting the limitations of the chosen analysis parameters.

Methodology
The study utilizes archival data from the HI4PI all-sky H I survey, which combines the EBHIS (northern) and GASS (southern) surveys. The computational workflow involves the following steps:

Data Processing: FITS spectral cubes are downloaded from the CDS/VizieR archive. The pipeline reconstructs velocity coordinates from FITS headers, applies Astropy WCS for sky coordinates, and extracts a specific longitude interval.
Spatial Averaging: To generate a single spectrum per longitude, pixels within a latitude strip of $|b| \le 2^\circ$ are averaged. This increases signal-to-noise but mixes gas at different vertical heights.
Terminal Velocity Extraction: A longitude-velocity ( $l-v$ ) diagram is constructed. The "terminal velocity" ( $v_{term}$ ) is defined operationally as the velocity edge where the brightness temperature ( $T_B$ ) drops below a fixed absolute threshold of 20 K. The algorithm requires persistence across adjacent channels to avoid noise-induced false detections.
Tangent-Point Geometry: For first-quadrant longitudes ( $0^\circ < l < 90^\circ$ ), the method assumes circular, axisymmetric motion. The Galactocentric radius ( $R_t$ ) is calculated geometrically as $R_t = R_0 \sin l$ , where $R_0 = 8.15$ kpc. The circular speed $\Theta(R_t)$ is derived using the formula $\Theta(R_t) = v_{term} + \Theta_0 \sin l$ , with $\Theta_0 = 236$ km s $^{-1}$ .
Calibration and Diagnostics: The workflow attempts to estimate an "edge-bias" correction by comparing measured terminal velocities near $l \approx 90^\circ$ (where $v_{term}$ should theoretically approach zero) against the geometric expectation. It also generates quality-control overlays comparing extracted points against the raw $l-v$ envelope.

Key Results

Rotation Curve: The derived rotation curve spans Galactocentric radii from approximately 1.5 to 7.4 kpc. It rises from roughly 95 km s $^{-1}$ at the innermost sampled region ( $R \approx 1.5$ kpc) to a plateau of approximately 230–255 km s $^{-1}$ in the intermediate-to-outer measured interval.
Detection Failures: The 20 K threshold and $|b| \le 2^\circ$ averaging configuration failed to yield valid terminal velocities for representative high-longitude sightlines (specifically near $l = 69.9^\circ$ and $84.9^\circ$ ). The outer emission wings in these spectra were too faint to cross the 20 K threshold.
Calibration Outcome: Due to the lack of valid detections in the high-longitude calibration interval (beginning around $l = 65^\circ$ ), no empirical edge-bias correction could be established. Consequently, the applied correction is 0.0 km s $^{-1}$ , which the authors explicitly interpret as "no correction measured" rather than proof of an unbiased edge.
Mass Diagnostics: The inferred enclosed mass (assuming spherical symmetry) grows from $\sim 10^9 M_\odot$ in the inner region to $\sim 10^{11} M_\odot$ near 7 kpc. The curve remains significantly higher than a Keplerian point-mass model, illustrating the need for extended mass distributions.

Key Contributions

Reproducible Undergraduate Investigation: The paper presents a complete, auditable workflow using open archival data (HI4PI) and standard Python libraries (Astropy, NumPy, Matplotlib). It emphasizes transparency by including run records, checksums, and configuration parameters.
Diagnostic Transparency: Unlike studies that present only a final, smoothed rotation curve, this work prioritizes the visibility of failure modes. It explicitly displays spectra where the algorithm fails (NaN results) and overlays extracted points on the raw $l-v$ diagram to show where the method succeeds or deviates from the emission envelope.
Parameter Sensitivity: The study demonstrates that the resulting rotation curve is not an immutable property of the data but depends heavily on analysis choices, specifically the brightness threshold (20 K vs. lower values) and the latitude averaging width ( $2^\circ$ vs. narrower strips).

Significance and Claims
The authors position this work as a pedagogical demonstration of the astrophysical power of 21-cm kinematics and the critical importance of quality control in observational astronomy.

Methodological Validity: The primary success is the demonstration that archival data can be processed to reveal coherent Galactic structure and a plausible inner-Galaxy rotation curve. The diagnostic chain—from $l-v$ maps to individual spectra and threshold crossings—makes the inference process visually auditable.
Limitations as Data: The paper claims that the inability to calibrate the edge bias at high longitudes is a scientifically useful result in itself. It highlights that a "plausible-looking" curve can be generated even when the underlying calibration is incomplete. The authors argue that a scientifically useful pipeline must expose its limitations (such as missing high-longitude detections) as clearly as its successes.
Educational Value: The work serves as a case study in distinguishing between obtaining a visually plausible plot and obtaining a fully validated measurement. It illustrates that high signal-to-noise data (HI4PI) can still yield systematically uncertain derived quantities due to analysis choices (thresholds, averaging) and astrophysical complexities (non-circular motions, the Galactic bar).
Scope: The authors explicitly state this is not a precision Galactic mass model. The results are limited to the inner Galaxy ( $R < 7.4$ kpc) and are subject to the assumptions of axisymmetry and circular motion, which are known to break down near the Galactic center and bar. The comparison with outer-Galaxy stellar kinematics (Eilers et al., 2019) is provided for context only, not as a validation of the H I-derived curve.

Technical Summary: Mapping the Milky Way's Rotation Curve from Archival 21-cm H I Observations

More like this