Biconical Outflow Model in Astrophysics

• Biconical outflow models are defined by two oppositely directed cones with specific geometric parameters such as half-opening angle, radial extents, and inclination.
• The models incorporate various velocity laws—constant, accelerating, or decelerating—to explain observed double-peaked emission-line profiles and projection effects.
• They are applied to estimate mass outflow rates, kinetic energy, and feedback strengths in systems ranging from AGNs and starbursts to galactic superwinds.

A biconical outflow model describes astrophysical winds or jet structures consisting of two oppositely directed conical regions—bicones—originating from a central engine (such as an active galactic nucleus, nuclear starburst, or young stellar object) and propagating into the surrounding interstellar or circumgalactic medium. These models are essential for interpreting the morphology, kinematics, mass, and feedback effects of outflows observed in emission and absorption lines across diverse astrophysical contexts, from galactic superwinds to AGN narrow-line regions (NLRs) and galactic-scale jets.

1. Geometric Definitions and Core Model Structure

A standard biconical outflow model consists of two cones, each characterized by a half-opening angle (θ), an axis orientation (position angle and inclination relative to the observer), and a spatial extent defined by minimum and maximum radii (r_in, r_out) from the nucleus. The apex of both cones coincides with the central source.

• Mathematical boundaries: For an axisymmetric bicone along the z-axis:

$$x^2 + y^2 \leq (z \tan\theta)^2, \quad r_{\rm in} \leq \sqrt{x^2 + y^2 + z^2} \leq r_{\rm out}$$

• Orientation: The model provides the capability to describe arbitrary inclinations (i) with respect to the line of sight, enabling synthetic observables as functions of projection effects (Nevin et al., 2017, Bae et al., 2016).
• Variants: Hollow bicones (with inner and outer opening angles), filled bicones, and nested or asymmetric bicones are all well explored in the literature (Nevin et al., 2017, Bae et al., 2016).

A summary of observationally derived geometric parameters from recent studies is presented below:

Object/System	Half-Opening Angle (θ)	Radial Extent	Inclination (i)	Reference
NGC 7172 (AGN; neon lines)	60°	0.5–0.7 kpc	~70° from LOS	(Muñoz et al., 2024)
F08572+3915 NW (ULIRG, H₂)	25° ± 5°	0.3–1.4 kpc	~45° (adopted estimate)	(Dan et al., 2024)
Milky Way Fermi Bubbles	55°	2–6.5 kpc	perpendicular to disk	(1412.14801612.01578)
3C 191 (quasar, [O III])	30–50°	5–11 kpc	Not directly measured	(Zhao et al., 26 Feb 2025)
PSOJ183+05 (z~6.4 QSO, [CII])	60°	0.3–4.8 kpc	Perpendicular to disk	(Bischetti et al., 21 Apr 2025)
Cygnus A (NLR, JWST)	46–54°	1.2–1.6 kpc	~60°	(Ogle et al., 10 Feb 2025)

Opening angles, inclinations, and physical sizes vary widely but are tightly constrained in systems with spatially resolved IFU or high-resolution spectroscopy.

2. Kinematic Prescriptions and Line-of-Sight Velocity Mapping

Biconical outflow models incorporate parametric or analytic radial velocity laws to describe the flow of gas within the cones. Canonical forms include:

• Constant velocity: $v(r) = v_{\rm out} = \text{const}$, a robust approximation for outflows lacking strong acceleration or deceleration over the observed range (Fox et al., 2014, Zhuoqi et al., 16 Mar 2025).
• Acceleration/deceleration: $v(r) = v_{\max} (1 - e^{-r/r_{\rm acc}})$ for an acceleration phase, or piecewise linear acceleration-deceleration with a turnover radius $r_t$ (Nevin et al., 2017, Dan et al., 2024).
• Radial projection: The observed line-of-sight velocity at position $r$ and polar angle $\theta$ (relative to the cone axis and LOS) is

$$v_{\rm LOS}(r, \theta, \phi) = v(r) \cdot \mathbf{\hat{r}} \cdot \mathbf{\hat{z}}$$

or,

$$v_{\rm LOS} = v(r) [\sin\theta \sin\phi \sin i + \cos\theta \cos i]$$

enabling synthetic velocity field and profile predictions as functions of geometry (Nevin et al., 2017, Bae et al., 2016, Ogle et al., 10 Feb 2025).

• Line decomposition: In spatially integrated spectra, biconical outflows commonly produce double-peaked line profiles (blueshifted and redshifted peaks), with Gaussian decomposition used to extract centroids, widths, and flux ratios (Cheng et al., 8 Nov 2025, Bizyaev et al., 2022).
• Monte Carlo inversion: Parameter estimation via Markov Chain Monte Carlo (MCMC) is employed to fit model parameters to slit or IFU velocity data (Nevin et al., 2017).

3. Physical Quantities: Outflow Mass, Rate, and Energetics

Conversion of observed luminosities and velocities to mass outflow rates (Ṁ), kinetic power ($\dot{E}_{\rm kin}$), and momentum flux (Ṗ) employs conical geometry and gas diagnostics:

• Mass outflow rate:

$$\dot{M}_{\rm out} = \Omega\, \frac{M_{\rm out}\, v_{\max}}{R}$$

where Ω is the fractional solid angle ($\Omega=1-\cos\theta$ per cone), $M_{\rm out}$ is gas mass, $v_{\max}$ is bulk velocity, and $R$ is the characteristic radius (Bischetti et al., 21 Apr 2025, Dan et al., 2024).

• Ionized and molecular mass: Derived using line fluxes (e.g., [O III], [C II], H$_2$) and electron or molecular density measurements, typically through emission line diagnostics or fine-structure ratios (Muñoz et al., 2024, Dan et al., 2024, Zhao et al., 26 Feb 2025).
• Kinetic energy and momentum:

$$\dot{E}_{\rm kin} = \frac{1}{2} \dot{M}_{\rm out} v_{\max}^2 \qquad \dot{p} = \dot{M}_{\rm out} v_{\max}$$

allowing for direct assessment of the feedback potential relative to AGN bolometric luminosity (Nevin et al., 2017).

Observed outflow rates and energetics display wide dispersion across systems:

System	Ṁ (M$_\odot$ yr$^{-1}$)	$v_{\max}$ (km/s)	$\dot{E}_{\rm kin}$ (erg/s)	Phase	Reference
NGC 7172 (AGN)	$\sim0.03$	400–600	$10^{39}$–$10^{40}$	Ionized	(Muñoz et al., 2024)
F08572+3915 NW (ULIRG)	$0.16$	1100	$6\times10^{40}$	H$_2$	(Dan et al., 2024)
3C 191 (Quasar)	$9.5$–13.4	800–975	$2.6$–$3.7\times10^{42}$	[O III]	(Zhao et al., 26 Feb 2025)
Milky Way Fermi Bubbles	$>0.2$	1000–1300	$>6\times10^{55}$ (total kinetic E)	UV/ionized	(Bordoloi et al., 2016)
PSOJ183+05 (z~6.4 QSO)	930	790	Not stated	[C II]	(Bischetti et al., 21 Apr 2025)
Cygnus A (NLR, Fe II)	40	150	Not stated	Ionized	(Ogle et al., 10 Feb 2025)

4. Observational Diagnostics and Model Selection

Biconical outflow model parameters are constrained through multi-faceted observations and quantitative fitting:

• Spatially resolved velocity and dispersion fields: IFU spectroscopy of emission lines ([O III], H$\alpha$, [C II], H$_2$) enables mapping of kinematic bicones and direct measurement of opening angles and extents (Zhao et al., 26 Feb 2025, Bischetti et al., 21 Apr 2025, Ogle et al., 10 Feb 2025).
• Spectral decomposition: Profile fitting (often with multi-Gaussian models) gives velocities, dispersions, and flux ratios of blue/red components. F-test statistics assess the number of components required (Cheng et al., 8 Nov 2025).
• Morphology-kinematics mapping: Shell shapes and limb-brightened cones in emission/absorption are used to differentiate between rotational and outflow models. For example, asymmetric double-peaked lines or global blueshifts/faint red wings favor biconical outflows over rotating disks or dual AGN (Nevin et al., 2017, Cheng et al., 8 Nov 2025).
• Polarization and radiative transfer: In starburst-driven winds, polarization mapping of scattered nuclear emission by biconical dust outflows constrains cone inclination and deceleration (Yoshida et al., 2010).
• Line ratio diagnostics: Electron density, temperature, and phase structure are constrained via line ratios (e.g., [Ne V]14/24 μm), crucial for mass and energetic calculations (Muñoz et al., 2024, Ogle et al., 10 Feb 2025).
• Absorption-line profile synthesis: Analytical biconical models predict the shape and depth of absorption/emission features in resonant lines, enabling fitting via MCMC for geometry and kinematics (Carr et al., 2018).
• Comparison to alternative models: In most systems, outflow models are tested against rotation-dominated and dual-nucleus scenarios, and are selected based on reproduction of asymmetry and velocity amplitude (Nevin et al., 2017, Cheng et al., 8 Nov 2025, Zhuoqi et al., 16 Mar 2025).

5. Model Variants, Extensions, and Physical Regimes

Advanced biconical outflow models incorporate further complexity:

• Hollow and nested bicones: Allow modeling of physically vacant interior cones or nested conical structures, explaining multi-component kinematic data (Bae et al., 2016, Nevin et al., 2017).
• Density and velocity gradients: Non-uniform profiles ($n_e(r) \propto r^{-\beta}$, $v(r) \propto r^{\alpha}$) influence projected emission, mass rates, and kinetic flux calculations (Dan et al., 2024, Bae et al., 2016, Muñoz et al., 2024).
• Rotation and jet-driven flows: Inclusion of rotation (e.g., spiral streamlines in Cygnus A) or jet interaction (bullets/streamers) adds diagnostic power for interpreting features not explained by simple radial outflow (Ogle et al., 10 Feb 2025).
• Thin shell and external interaction: In models of jet/ambient interaction (e.g., Mira, starburst superwinds), the momentum-balance at the shell determines the observed bow shock and cavity shape (López-Cámara et al., 2011).
• Impact of dust and disk extinction: Forward modeling includes dust attenuation, with bicone–disk overlap producing asymmetric or nested profiles (Bae et al., 2016, Nevin et al., 2017).
• Time dependence and feedback: Outflow ages are constrained by dividing observed spatial scales by bulk velocity, e.g., Fermi Bubbles (age ≈ 2.5–9 Myr) (Fox et al., 2014, Bordoloi et al., 2016).

6. Model Constraints, Feedback, and Astrophysical Implications

Large-sample statistical analyses and model–data comparison put quantitative limits on the prevalence, velocity range, opening angle, and feedback potential of biconical outflows:

• AGN-driven winds in type 2 AGNs: Monte Carlo model grids and MCMC sampling (incorporating disk and outflow components, opening angle, and extinction) recover launching speeds of $v_{\max}\sim 250$–350 km/s (bulk) to $\sim1000$–1500 km/s (strongest $\sim$2–5%) and typical outer half-opening angles $\theta_{\text{out}}\sim 30$–40°, rarely exceeding 50° (Kim et al., 11 Nov 2025, Bae et al., 2016).
• Energetic efficiency: In moderate-luminosity AGNs, kinetic energy injection typically remains 0.01%–0.1% of $L_{\mathrm{bol}}$, below thresholds required for effective feedback or quenching, though individual systems reach higher values (Nevin et al., 2017, Zhao et al., 26 Feb 2025).
• Feedback nature: There is robust evidence for both positive (outflow-triggered star formation in the CGM) and negative feedback (gas mass-loading, global quenching), varying by phase, system, and spatial scale (Muñoz et al., 2024, Bischetti et al., 21 Apr 2025).
• Limitations from spatial resolution and seeing: Statistical size–luminosity correlations (e.g., $R_{\rm out}\propto L_{\rm [O III]}^{0.30}$) must consider seeing effects, which can bias outflow radius determination by a factor of 2–5 in unresolved objects (Kim et al., 11 Nov 2025).

7. Summary Table: Representative Biconical Outflow Model Parameters

Model Feature	Range / Formula	Typical Value(s)	Key Reference(s)
Half-opening angle	θ	25–60°, up to 120° (full)	(Dan et al., 2024, Muñoz et al., 2024, Bischetti et al., 21 Apr 2025)
Outflow velocity	$v_{\rm max}$, $v(r)$	100–1300 km/s (bulk); up to 3000	(Dan et al., 2024, Bordoloi et al., 2016)
Mass-outflow rate	$\dot{M}$	0.03–930 $M_\odot$/yr	(Muñoz et al., 2024, Bischetti et al., 21 Apr 2025)
Kinetic power	$\dot{E}_{\rm kin}$	$10^{39}$–$10^{42}$ erg/s	(Muñoz et al., 2024, Zhao et al., 26 Feb 2025)
Size (outer radius)	$r_{\rm out}$	0.5–10 kpc	(Dan et al., 2024, Zhao et al., 26 Feb 2025)

The biconical outflow framework provides a robust, physically motivated, and mathematically explicit model that underpins the interpretation of emission and absorption line features arising from nuclear and galactic-scale feedback, allowing derivation of kinematic, geometric, and energetic parameters essential for understanding galaxy evolution and the baryon cycle (Nevin et al., 2017, Bae et al., 2016, Muñoz et al., 2024, Dan et al., 2024, Kim et al., 11 Nov 2025).