Research papers | Newest first

Maps and simulations of the nearby Milky Way

This page collects six papers on ionized gas, three-dimensional maps, and the local Galactic magnetic field, with a few figures from each.

Scroll to the latest paper Milky Way image: ESA/Gaia/DPAC, Stefan Payne-Wardenaar · CC BY-SA 3.0 IGO

2026

H I density and radial velocity

A velocity-resolved 3D map of atomic hydrogen within 1.25 kpc

McCallum et al. (2026), Monthly Notices of the Royal Astronomical Society

Read the paper

The 21-cm sky gives us velocity but not distance. This reconstruction separates nearby emission from distant Galactic structure and places the local atomic gas into a three-dimensional volume around the Sun.

Top-down comparison of reconstructed atomic hydrogen and dust column density around the Sun
Reconstructed H I and dust column density. The same clouds and cavities appear independently in both tracers.
Observed, reconstructed, distant, and local components of the 21-centimetre sky
Decomposing the observed sky into local gas inside 1.25 kpc and the more distant Milky Way.
Top-down density-weighted radial velocity of reconstructed atomic hydrogen
Density-weighted line-of-sight velocity shows both Galactic rotation and coherent local motions.
Integrated atomic-hydrogen density in the Galactic midplane and several height intervals Atomic hydrogen density and residual velocity across several heights above and below the Galactic plane
H I density and residual radial velocity in matching height intervals above and below the Galactic midplane.

2026

Galactic magnetic field

The radial component of the local Galactic magnetic field in 3D

McCallum et al. (2026), Monthly Notices of the Royal Astronomical Society

Read the paper

This reconstruction combines my electron-density map with the Hutschenreuter et al. Faraday-rotation sky to infer how the line-of-sight magnetic field runs through the local volume. The animation follows samples from the posterior, showing both the recovered structure and where the reconstruction remains uncertain.

Samples of the reconstructed radial magnetic field. Positive Br points towards the observer; negative Br points away.

2025

Multiwavelength ionized gas

A 3D, multiwavelength view of the Milky Way's local ionized gas

McCallum, Wood, Benjamin, Krishnarao & Vandenbroucke (2025), MNRAS

Read the paper

Using the Edenhofer 3D dust map and a census of nearby O stars, this simulation follows ionizing photons through the local ISM and predicts Hα, [S II], and [O III] emission. The lines respond differently to temperature and ionization, making it easier to compare where massive stars ionize gas and where recent supernova shocks affect the emission.

Face-on false-colour view of local ionized gas with named H II regions
Face-on emission-line map: Hα in red, [S II] in green, and [O III] in blue. The labels identify several familiar nearby nebulae in the same coordinate system.
Face-on evolution of H alpha, sulphur, and oxygen emission in a radiation-hydrodynamics simulation
A separate face-on radiation-hydrodynamics simulation, shown in the same Hα, [S II], and [O III] colour channels as the static map.

Stellar environment

The surrounding gas changes how far a star's photons travel

Ionizing luminosity is not the only factor. ζ Puppis and the Bajamar star emit nearly the same number of ionizing photons, but their environments give them very different reach. In this model, ζ Puppis influences a volume 827 times larger, showing how strongly a star's environment controls its impact on the wider ISM.

Face-on and edge-on comparison of the ionized volumes reached by zeta Ophiuchi, the Bajamar star, and zeta Puppis
Where Lyman-continuum photons from three individual stars are absorbed, shown face-on and edge-on against the surrounding gas.

The same paper, followed through time

Birth cloud, H II region, supernova, remnant

A typical massive star first lights and disperses its birth cloud. Its supernova then flashes as an [O III]-bright shell, cools into [S II], and finally folds back into the turbulent interstellar medium.

Six time steps following the environment of a massive star from before birth through its supernova remnant
Six moments from 19.1 to 25.0 Myr, in chronological order from upper left to lower right.

2025

The local Hα sky

The Hα sky in three dimensions

McCallum et al. (2025), Monthly Notices of the Royal Astronomical Society: Letters

Read the paper

Starting from a 3D dust map and the nearby O-star population, the model predicts the Hα sky from the Sun. Comparing the result with WHAM provides a check on how well the mapped gas and nearby stars reproduce the observed emission.

All-sky comparison of observed WHAM H alpha, simulated H alpha, dust-scattered H alpha, and the simulation quality factor
The observed WHAM sky, the simulated Hα sky, the dust-scattered contribution, and the difference between the model and observations.
Observed and simulated H alpha views of the Orion-Eridanus superbubble
Orion-Eridanus in the observed WHAM sky and in the simulated sky.
Top-down simulated H alpha and neutral hydrogen maps around the Sun
Top-down Hα and neutral hydrogen connect named sky regions to their physical surroundings.

2024

Time-dependent metal ionization

Persistent collisionally excited emission from diffuse ionized gas

McCallum, Wood, Benjamin, Krishnarao & Vandenbroucke (2024), MNRAS

Read the paper

Oxygen does not move instantly between ionization states as the gas heats and cools. The result is an onion-skin structure around hot bubbles, with high ions surviving in thin gas long after a simple equilibrium model would erase them.

Slices through five oxygen ion fractions and gas temperature in a tall-box galaxy simulation
O I through O V+ and temperature after 25 Myr. Successive ions trace different layers around hot, low-density structures.

2024

High-altitude diffuse ionized gas

The persistence of high-altitude non-equilibrium ionized gas

McCallum et al. (2024), Monthly Notices of the Royal Astronomical Society

Read the paper

In low-density gas, recombination can take millions of years. Following that delay directly produces a stable, extended ionized layer; forcing the gas into equilibrium instead makes the layer flicker with the short lives of its stars.

Chronological sequence of projected gas column density during a galactic inflow and outflow cycle
The first inflow-outflow cycle: gas leaves the midplane, reaches several kiloparsecs, and returns.
Time-dependent radiation-hydrodynamics simulation of diffuse ionized gas above a galactic midplane
The evolving tall-box simulation. The ionized layer survives between generations of short-lived massive stars.
Time sequence comparing non-equilibrium and equilibrium ionized hydrogen column density
A direct comparison: non-equilibrium gas stays extended and stable, while the equilibrium model flickers.