Column-density estimation from the equivalent widths of absorption doublets
Abstract
Na I D absorption is widely used to estimate the dust extinction toward Galactic and extragalactic sources, but the inferred column densities are unreliable once the lines saturate. We present an analytical framework for estimating column densities from the equivalent widths (EWs) of absorption-line doublets. The method requires only that the two doublet members be measured separately -- the velocity structure need not be resolved -- so a resolving power λ/Δλ≳103 suffices for Na I D, two orders of magnitude below velocity-resolved methods. For the Gaussian curve of growth, the classical inversion reduces to a universal saturation correction depending only on the doublet ratio, R=EW2/EW1: at fixed R, the inferred column density is linear in the measured EW. The formulation yields the exact Gaussian inversion, Nκ̃, which provides the Doppler parameter; a closed-form second-order approximation; and a strict lower bound, Nmin, valid for an arbitrary velocity structure. We recommend reporting Nmin: an optimal, assumption-free lower bound, within a few per cent of Nκ̃ for R≳1.4 and at most a factor $3.5$ below it at R=1.1. The ratio Nκ̃/N2nd provides a per-object saturation diagnostic. We validate the framework against published profile-fitting Na I D column densities toward Type Ia supernovae, using only the integrated EWs. An explicit closed-form formula converts the two EWs into an extinction estimate through the observed log N(Na I)--AV relation. A public implementation with full uncertainty propagation accompanies the paper.
Versions
Related
- Where the Forest Goes Dark: Halo-centered Characterization of High Column Density Systems with IllustrisTNG
- Emergence of high-mass stars in complex fiber networks (EMERGE) VII. Multi-Tracer Column Density Synthesis
- A measurement-based approach to equivalent power flux density calculation, to quantify the impact of satellite mega-constellations on low-frequency radio astronomy
- Disc-Based Estimation of the Local Dark Matter Density and Velocity Distribution in IllustrisTNG50
- SPHEREx 3.3 Micron Aromatic Emission in the Magellanic Clouds: Separating Dust Column, Excitation, and Environmental Suppression
BibTeX
@misc{badash2026columndensity,
title = {{Column-density estimation from the equivalent widths of absorption doublets}},
author = {Avshalom Badash and Doron Kushnir},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2608.04088},
eprint = {2608.04088},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.04088},
}
RIS
TY - UNPB
TI - Column-density estimation from the equivalent widths of absorption doublets
AU - Avshalom Badash
AU - Doron Kushnir
PY - 2026
DA - 2026/10/30
PB - arXiv
DO - 10.48550/arxiv.2608.04088
UR - https://arxiv.org/abs/2608.04088
AB - Na I D absorption is widely used to estimate the dust extinction toward Galactic and extragalactic sources, but the inferred column densities are unreliable once the lines saturate. We present an analytical framework for estimating column densities from the equivalent widths (EWs) of absorption-line doublets. The method requires only that the two doublet members be measured separately -- the velocity structure need not be resolved -- so a resolving power λ/Δλ≳10^3 suffices for Na I D, two orders of magnitude below velocity-resolved methods. For the Gaussian curve of growth, the classical inversion reduces to a universal saturation correction depending only on the doublet ratio, R=EW_2/EW_1: at fixed R, the inferred column density is linear in the measured EW. The formulation yields the exact Gaussian inversion, N_κ̃, which provides the Doppler parameter; a closed-form second-order approximation; and a strict lower bound, N_min, valid for an arbitrary velocity structure. We recommend reporting N_min: an optimal, assumption-free lower bound, within a few per cent of N_κ̃ for R≳1.4 and at most a factor $3.5$ below it at R=1.1. The ratio N_κ̃/N_2nd provides a per-object saturation diagnostic. We validate the framework against published profile-fitting Na I D column densities toward Type Ia supernovae, using only the integrated EWs. An explicit closed-form formula converts the two EWs into an extinction estimate through the observed log N(Na I)--A_V relation. A public implementation with full uncertainty propagation accompanies the paper.
ER -