Table of Contents
Fetching ...

OutLines: Modeling Spectral Lines from Winds, Bubbles, and Outflows

Sophia R. Flury

TL;DR

OutLines delivers a physically driven framework for modeling emission and absorption lines from winds, bubbles, and outflows, addressing limitations of empirical and multi-Gaussian approaches. By employing the Sobolev approximation, multiple velocity-field choices (including CAK β-law, AccPlaw, and VelPlaw), and diverse radial density profiles across spherical, filled-cone, and hollow-cone geometries (with disks and cavities), the code links spectral features to wind physics and geometry. The authors validate the method on four astrophysical cases—an H II region knot, a Green Pea-like super star cluster, a starburst galaxy, and an AGN—deriving physically meaningful outflow properties and insights into feedback and LyC escape. As an openly available tool, OutLines enables robust, physics-based inference from large spectroscopic surveys, with potential to inform studies of baryon cycling, galaxy evolution, and AGN-driven feedback in the era of WEAVE-LOFAR and 4MOST/WAVES.

Abstract

Common methods for studying the kinematics and geometry of outflowing gas rely on modeling emission and absorption lines in integrated spectra using methods that are not physically motivated, including empirical quantiles or fitting multiple Gaussian or Voigt profiles. Such methods are not always consistent with the interpretation of these features and, as a result, miss key underlying physics and can even lead to inaccurate interpretations of observations. To address this problem, we present the publicly available python code OutLines, which provides astrophysical models of spectral emission and absorption line profiles produced by outflows in a variety of environments. The OutLines code accounts for differences in parameterization of the velocity field and density profile while allowing for different outflow geometries, making OutLines versatile and useful for a wide variety of astrophysical phenomena. We demonstrate the wide applicability of OutLines by using the code to model line profiles in an H II region knot, super star clusters, a starburst galaxy, and an AGN. In each of these contexts, we illustrate how OutLines can illuminate key underlying physics in ways that improve our scientific understanding and address important open questions in astronomy, including the key mechanisms in the baryon cycle, the evolution of H II regions and galaxies, and even Lyman continuum escape. OutLines will be a critical resource as massively multiplexed spectroscopic surveys like WEAVE-LOFAR and 4MOST/WAVES come online, providing the means to probe feedback kinematics with deeper, higher resolution spectroscopy for unprecedented large samples of galaxies.

OutLines: Modeling Spectral Lines from Winds, Bubbles, and Outflows

TL;DR

OutLines delivers a physically driven framework for modeling emission and absorption lines from winds, bubbles, and outflows, addressing limitations of empirical and multi-Gaussian approaches. By employing the Sobolev approximation, multiple velocity-field choices (including CAK β-law, AccPlaw, and VelPlaw), and diverse radial density profiles across spherical, filled-cone, and hollow-cone geometries (with disks and cavities), the code links spectral features to wind physics and geometry. The authors validate the method on four astrophysical cases—an H II region knot, a Green Pea-like super star cluster, a starburst galaxy, and an AGN—deriving physically meaningful outflow properties and insights into feedback and LyC escape. As an openly available tool, OutLines enables robust, physics-based inference from large spectroscopic surveys, with potential to inform studies of baryon cycling, galaxy evolution, and AGN-driven feedback in the era of WEAVE-LOFAR and 4MOST/WAVES.

Abstract

Common methods for studying the kinematics and geometry of outflowing gas rely on modeling emission and absorption lines in integrated spectra using methods that are not physically motivated, including empirical quantiles or fitting multiple Gaussian or Voigt profiles. Such methods are not always consistent with the interpretation of these features and, as a result, miss key underlying physics and can even lead to inaccurate interpretations of observations. To address this problem, we present the publicly available python code OutLines, which provides astrophysical models of spectral emission and absorption line profiles produced by outflows in a variety of environments. The OutLines code accounts for differences in parameterization of the velocity field and density profile while allowing for different outflow geometries, making OutLines versatile and useful for a wide variety of astrophysical phenomena. We demonstrate the wide applicability of OutLines by using the code to model line profiles in an H II region knot, super star clusters, a starburst galaxy, and an AGN. In each of these contexts, we illustrate how OutLines can illuminate key underlying physics in ways that improve our scientific understanding and address important open questions in astronomy, including the key mechanisms in the baryon cycle, the evolution of H II regions and galaxies, and even Lyman continuum escape. OutLines will be a critical resource as massively multiplexed spectroscopic surveys like WEAVE-LOFAR and 4MOST/WAVES come online, providing the means to probe feedback kinematics with deeper, higher resolution spectroscopy for unprecedented large samples of galaxies.

Paper Structure

This paper contains 44 sections, 49 equations, 28 figures, 6 tables.

Figures (28)

  • Figure 1: Radial velocity fields included in OutLines and given by Equations \ref{['eqn:betaCAK']}-\ref{['eqn:velplaw']}. Acceleration from rest at the base $R_0$ of the wind ($x=1$) proceeds until the wind reachs the terminal velocity $v_\infty$ ($w[x]=1$). Velocity fields include the 1986AA...164...86P$\beta$ law approximation to CAK theory 1975ApJ...195..157C (blue, Equation \ref{['eqn:betaCAK']}), the 2010ApJ...717..289S power law acceleration (green, Equation \ref{['eqn:accplaw']}), and the power law velocity field (purple, Equation \ref{['eqn:velplaw']}). Dash length increases with increasing $\beta$.
  • Figure 2: Nebular emission ( top) and resonant absorption ( bottom) line profiles for each of the velocity fields included in OutLines assuming $\beta=1.5$ (CAK theory), $\beta=1$ (velocity power law), and $\beta=2$ (acceleration power law) for the OutLines defaults.
  • Figure 3: Radial density profiles included in OutLines and given by Equations \ref{['eqn:denplaw']}-\ref{['eqn:fred']}. Pulse-like density profiles, including the normal ('Normal', green, Eqn \ref{['eqn:norm']}), log normal ('LogNormal', yellow, Eqn \ref{['eqn:lognorm']}), and fast-rise exponential decay ('FRED', red, Eqn \ref{['eqn:fred']}) distributions, represent physical scenarios such as blast waves, expanding shells, and single burst events. For visualization, the pulse-like profiles are arbitrarily set to maximize density at $r=3R_0$. Continuous-like density profiles, including power law ('PowerLaw', purple, Eqn \ref{['eqn:denplaw']}), exponential ('Exponential', orange, Eqn \ref{['eqn:expon']}), and double power law ('DoublePowerLaw', blue, Eqn \ref{['eqn:dplaw']}) distributions, represent an integral over many pulses and will always maximize density at the base of the outflow where $r=R_0$.
  • Figure 4: Nebular emission ( top) and absorption ( bottom) line profiles predicted by OutLines for fixed density, column density, and velocity field for different density distributions assuming a spherical geometry.
  • Figure 5: Nebular emission ( top) and absorption ( bottom) line profiles predicted by OutLines for a series of expanding bubbles, shells, or other pulse-like episodes in cases of unresolved (orange) or resolved (blue) shells with (dashed) and without (solid) damping to mimic the effects of dissipation. In all cases, density, column density, velocity field, and terminal velocity are fixed, and a spherical geometry is assumed.
  • ...and 23 more figures