---
title: 'SimBAL: Bayesian Modeling for BAL Quasar Spectra'
url: https://www.emergentmind.com/topics/simbal
type: topic
---

# SimBAL: Bayesian Modeling for BAL Quasar Spectra

Searching arXiv for papers on SimBAL and its later applications.
arxiv_search(query="SimBAL quasar Leighly broad absorption line", max_results=10, sort_by="relevance")
SimBAL is a spectral-synthesis and forward-modeling framework for broad absorption line (BAL) quasar spectroscopy that infers physical outflow properties by comparing synthetic spectra, generated from photoionization calculations, directly to observed data. It was introduced to address regimes in which conventional apparent-optical-depth or line-by-line analyses become unreliable, especially when BAL troughs are broad, blended, saturated, and affected by partial covering [1808.02441]. In SimBAL, grids of ionic column densities computed with Cloudy are mapped into synthetic absorption spectra and fit in a Bayesian framework, typically with MCMC, yielding constraints on ionization parameter, density, column density, kinematics, and inhomogeneous covering structure as functions of velocity [1808.02441]. Subsequent work has extended the method to multiwavelength covering-fraction tests in LoBAL quasars [1811.04174], highly blended FeLoBAL and overlapping-trough systems [2001.07347, 2203.11964], multi-epoch variability analyses [2405.06027], high-redshift feedback studies [2404.12443], and extremely high-velocity outflows [2508.14221].

## 1. Conceptual basis and motivation

SimBAL was introduced as a reversal of the standard BAL-analysis workflow. Rather than first measuring apparent optical depths or fitting individual troughs and only then comparing those measurements to photoionization models, it starts from photoionization models, constructs synthetic spectra directly, and compares those spectra with the observed data [1808.02441]. This design targets a longstanding problem in BAL spectroscopy: the most diagnostic transitions are often blended, partially saturated, and subject to partial covering, so a decomposition into isolated ionic measurements can be ambiguous or impossible.

The original motivation was empirical as well as methodological. Earlier BAL analyses, including approaches based on measured line complexes followed by Cloudy comparison, could yield physically plausible ionic columns yet fail to reproduce the full synthetic spectrum [1808.02441]. SimBAL was designed to avoid that inconsistency by fitting the spectrum itself. This also permits use of constraints from nondetections: the absence of specific transitions enters the likelihood naturally, rather than being discarded as noninformation [1808.02441].

This suggests that SimBAL is best understood not merely as a line-fitting code but as an inference engine that couples photoionization physics, radiative transfer parameterization, and Bayesian calibration. In the published applications, this architecture is repeatedly used where blending is so severe that simpler diagnostics either become lower limits or fail structurally, as in FeLoBALs [2001.07347, 2203.11964] and extremely broad EHVO troughs [2508.14221].

## 2. Core model architecture

At its core, SimBAL uses large grids of ionic column densities generated by Cloudy as a function of the physical state of the gas. In the original implementation, the Cloudy C13.03 grids span \(-4.0 < \log U < 2\), \(2.8 < \log n < 9\), and \(21.5 < \log N_H - \log U < 23.7\), with ionic columns extracted for 179 ground and excited states and combined with a line list of 6267 transitions [1808.02441]. Later FeLoBAL work using Cloudy C17 expanded this to 76488 transitions, 281 ions, and 997 ionic states when excited states are counted, together with a flexible oversampling scheme near the hydrogen ionization front and 619,721 grid points [2001.07347].

The fitted physical parameterization varies somewhat by application, but the common parameter set includes ionization parameter \(\log U\), gas density \(\log n\), the column-density parameter \(\log N_H - \log U\), kinematic parameters, and a covering-fraction parameter \(\log a\) [1808.02441, 2508.14221, 2203.11964]. In some early work, the total column density \(N_H\) is also discussed directly as a derived quantity [1808.02441]. The synthetic spectrum is generated in velocity space, either with Gaussian opacity profiles or, more commonly in later work, with adjacent rectangular “tophat” bins that flexibly represent variation across the trough [1808.02441, 2508.14221, 2404.12443].

The comparison to data is Bayesian. SimBAL generates spectra in real time and compares them with observed spectra via a \(\chi^2\)-based likelihood explored with `emcee` MCMC [1808.02441]. The original paper describes flat priors on the gas parameters to confine sampling to the Cloudy grid and Gaussian priors on line offsets and widths based on inspection [1808.02441]. Posterior chains provide both parameter estimates and derived physical quantities such as radius, mass outflow rate, momentum flux, and kinetic luminosity.

A concise representation of the original logic is:

\[
\text{Cloudy} \rightarrow \text{ionic columns} \rightarrow \text{synthetic spectrum} \rightarrow \text{likelihood} \rightarrow \text{posterior}
\]

as explicitly described in the foundational study [1808.02441].

## 3. Treatment of partial covering and velocity structure

A defining feature of SimBAL is its use of an inhomogeneous, power-law partial-covering formalism. The optical depth is parameterized as

\[
\tau = \tau_{\max} x^a,
\]

where \(x\in(0,1)\) is a fractional surface-area coordinate and \(a\) is the covering-fraction index [1808.02441, 1811.04174, 2405.06027]. The parameter is fit as \(\log a\), and a larger \(a\) corresponds to a smaller effective covering fraction [1811.04174, 2508.14221]. This prescription was adopted because it is mathematically commutative and because it captures the empirical behavior that stronger lines often appear to cover more of the source than weaker lines [1808.02441].

The formalism is central to several SimBAL results. In the SDSS J0850+4451 studies, the model showed that covering fraction decreases strongly with velocity and is often the most strongly varying parameter controlling apparent trough morphology [1808.02441]. In the multiwavelength follow-up, extrapolation of the UV solution to optical and near-infrared data predicted \(\mathrm{Mg\,II}\), \(\mathrm{He\,I^*}\,3889\), and \(\mathrm{He\,I^*}\,10830\) troughs that were too deep, leading to the conclusion that the covering fraction is wavelength dependent and that the UV covering fraction is about 2.5 times higher than in the optical/NIR [1811.04174]. The long-wavelength fit required an upward shift of \(\Delta \log a \approx 0.4\), interpreted geometrically as evidence for a clumpy, structured absorber rather than a smooth screen [1811.04174].

The method also uses velocity-resolved parameterizations to recover internal outflow structure. The original “accordion model” represents the trough as a sequence of adjacent velocity bins; these bins are not interpreted as distinct physical clouds but as a flexible representation of changing outflow conditions with velocity [1808.02441]. In SDSS J0850+4451, the best-performing tophat accordion models improved markedly over single-Gaussian fits and revealed that \(\log U\) and column density increase with velocity while covering fraction decreases [1808.02441]. In the z\(\sim 6.6\) quasar J0923+0402, a 25-bin velocity-resolved tophat accordion fit from \(-32{,}000\) to \(-6{,}800\ \mathrm{km\ s^{-1}}\) showed that the observed LoBAL/high-ionization morphology could be explained without requiring strong velocity-dependent ionization changes; instead, column density and partial covering dominate the profile structure [2404.12443].

A plausible implication is that SimBAL’s velocity binning is less a decomposition into discrete absorbers than a regularized coordinate system for mapping opacity, ionization, and covering structure across the outflow.

## 4. Derived quantities and physical interpretation

SimBAL applications routinely convert fitted parameters into physical radii and energetics using the standard ionization-parameter relation

\[
U=\frac{Q}{4\pi R^2 n c},
\]

with \(Q\) the hydrogen-ionizing photon rate, \(R\) the absorber distance, \(n\) the gas density, and \(c\) the speed of light [1808.02441, 2405.06027, 2508.14221, 2404.12443, 2203.11964]. Once \(U\), \(n\), and \(Q\) are specified or constrained, the absorber distance follows. This is how SimBAL moves from spectral fitting to geometry.

A recurrent correction concerns the conversion from the column-density parameter to a covering-fraction-weighted total column density:

\[
\log N_H = (\log N_H - \log U)+\log U-\log\!\left(1+10^{\log a}\right),
\]

a relation used explicitly in the FeLoBAL survey and the EHVO analysis [2203.11964, 2508.14221]. This correction matters because modifying the covering-fraction structure changes the inferred total column and hence the energetics. In the J0850+4451 multiwavelength covering-fraction study, using the long-wavelength covering fraction lowered the best estimate of the total hydrogen column density by a factor of about 1.6 relative to the UV-only solution [1811.04174].

Mass outflow rate and kinetic luminosity are then computed with the standard thin-shell style relations

\[
\dot M = 8\pi \mu m_p \Omega R N_H v,
\qquad
\dot E_k = \frac{\dot M v^2}{2},
\]

with \(\mu=1.4\) and typically \(\Omega=0.2\) assumed as the global covering fraction [1808.02441, 2203.11964, 2405.06027, 2508.14221, 2001.07347]. Several studies emphasize the strong velocity leverage, effectively \(L_{\rm KE}\propto v^3\) under the adopted scaling [2508.14221, 2203.11964].

The physical interpretations vary by object class. In the original SDSS J0850+4451 analysis, SimBAL placed the outflow at 1–3 pc, found total \(\log N_H = 22.9\) for solar metallicity and 22.4 for \(Z=3Z_\odot\), and derived a mass outflow rate of 17–28 \(M_\odot\,\mathrm{yr^{-1}}\) with kinetic power 0.8–0.9% of \(L_{\rm Bol}\) [1808.02441]. In the overlapping-trough FeLoBAL SDSS J1352+4239, SimBAL inferred \(\log N_H = 23.22 \pm 0.05\), radius \(\sim 10\) pc, \(\log \dot M = 3.5 \pm 0.04\ M_\odot\,\mathrm{yr^{-1}}\), and \(\log L_{KE} = 48.1 \pm 0.04\ \mathrm{erg\,s^{-1}}\), with \(L_{KE}/L_{Bol}\sim 1\) [2001.07347]. In the high-redshift quasar J0923+0402, SimBAL delivered \(\log U = -0.32^{+0.01}_{-0.02}\), \(\log N_H = 22.57^{+0.01}_{-0.02}\), and a BAL radius bracketed as \(3 \lesssim R_{\rm BAL} \lesssim 210\) pc [2404.12443].

## 5. Observational domains and representative applications

SimBAL has been applied across a range of BAL phenomenology. The table summarizes representative studies documented in the literature.

| Object / study | SimBAL role | Principal result |
|---|---|---|
| SDSS J085053.12+445122.5 I [1808.02441] | Foundational HST/COS UV spectral synthesis | Velocity-resolved \(\log U\), \(N_H\), density, and covering-fraction structure; outflow at 1–3 pc |
| SDSS J085053.12+445122.5 II [1811.04174] | UV solution extrapolated to optical/NIR | UV covering fraction about 2.5 times higher than optical/NIR; clumpy structured absorber |
| SDSS J135246.37+423923.5 [2001.07347] | Overlapping-trough FeLoBAL forward modeling | Extremely fast FeLoBAL outflow; scattered light and radiation filtering required |
| 50 low-z FeLoBAL quasars [2203.11964] | First systematic sample analysis | Wide range of \(\log U\), density, and radius; 18% with feedback-level outflows |
| WPVS 007 [2405.06027] | First simultaneous multi-epoch SimBAL analysis | Variability primarily driven by changes in covering fraction |
| J0923+0402 at \(z\sim6.6\) [2404.12443] | BAL fit in multi-phase feedback study | Powerful compact LoBAL wind on \(\lesssim210\) pc scale |
| J164653.72+243942.2 [2508.14221] | EHVO best-estimate forward modeling with multi-epoch tying | Variability best explained by changing \(\log a\); pc-scale extremely high-velocity outflow |

The foundational application to SDSS J0850+4451 established the method’s capacity to recover velocity-resolved physical conditions from HST/COS far-UV data and showed that simple single-Gaussian opacity profiles were inadequate [1808.02441]. The follow-up study extended the analysis into rest-frame optical and near-IR bands, especially He I* \(\lambda10830\), to exploit the fact that the continuum-emitting region is much larger at longer wavelengths [1811.04174]. Because the 1 \(\mu\)m emission region was estimated to be about a factor of 12 larger in radius and about 140 in area than the 1100 Å region, the weaker long-wavelength absorption implied that absorber uniformity depends on source scale [1811.04174].

In FeLoBAL work, SimBAL has been particularly important because near-UV spectra can contain thousands of overlapping Fe II transitions. The single-object analysis of SDSS J1352+4239 showed that the code could simultaneously fit continuum, emission, absorption, reddening, and a scattered-light component in a heavily absorbed overlapping-trough quasar [2001.07347]. The broader survey of 50 low-redshift FeLoBAL quasars then used SimBAL uniformly to derive ranges of \(-4\lesssim\log U\lesssim 1.2\), \(2.8\lesssim\log n\lesssim8\), and \(0\lesssim\log R\lesssim4.4\) pc, finding no evidence for disk winds at \(R\ll0.01\) pc in that sample [2203.11964].

Time-domain applications represent another expansion. In WPVS 007, the first simultaneous multi-epoch SimBAL analysis fit pairs of epochs while allowing only one parameter class to vary and concluded that covering-fraction changes dominate the variability [2405.06027]. The EHVO quasar J1646 used a similar multi-epoch tying strategy across DR5 and DR9, and again the best model was the one in which only \(\log a\) changed, with \(\log U\) and \(N_H\) held fixed [2508.14221].

## 6. Scientific inferences enabled by SimBAL

Several scientific themes recur across SimBAL studies. The first is that covering fraction is often the dominant control on observed trough shape. In SDSS J0850+4451, the covering fraction decreases strongly with velocity [1808.02441]. In WPVS 007, BAL variability is determined primarily by changes in covering fraction rather than ionization [2405.06027]. In J1646, the extreme two-epoch variability is likewise best explained by changes in \(\log a\), not by changes in ionization or total column density [2508.14221]. These convergent results suggest that absorber substructure and projected source coverage are central observables rather than nuisance parameters.

A second theme is that multi-line forward modeling yields systematically stronger physical constraints than AOD baselines. The J1646 paper explicitly contrasts SimBAL with AOD: AOD fits only continuum-normalized absorption, does not model emission-line fill-in, relies mostly on C IV, and returns conservative lower limits, whereas SimBAL fits continuum, emission lines, and absorption together, uses multiple ions, and recovers hidden opacity beneath emission peaks [2508.14221]. In that case, AOD gives \(\log N_H > 21.24\) and \(>20.91\) in the two epochs, while SimBAL returns \(\log N_H = 21.79 \pm 0.06\) and \(21.63^{+0.05}_{-0.06}\) [2508.14221].

A third theme is that SimBAL can expose nontrivial outflow geometry. The multiwavelength J0850+4451 result rejects the “giant screen” picture because the absorber does not fully cover even the compact UV continuum source [1811.04174]. The favored interpretation is a distribution of small clumps, either internally structured or clustered, that more effectively cover the smaller UV source than the larger optical/NIR source [1811.04174]. In J0923+0402, the balance of high- and low-ionization BAL features is explained without invoking strong ionization stratification, instead pointing to velocity-dependent opacity and covering structure [2404.12443].

A fourth theme is feedback energetics. In the original J0850+4451 study, the momentum flux was consistent with \(L_{\rm Bol}/c\), and the kinetic luminosity was 0.8–0.9% of \(L_{\rm Bol}\) [1808.02441]. The FeLoBAL survey adopted \(L_{KE}/L_{Bol}>0.005\) as a feedback threshold and found that 8 of 55 outflowing BAL components meet it, corresponding to roughly 18% of the full sample [2203.11964]. The high-redshift J0923+0402 analysis concluded that the BAL wind carries \(0.02 \lesssim \dot E_{\rm BAL}/L_{\rm Bol} \lesssim 1.0\), consistent with efficient AGN feedback [2404.12443]. J1646 showed that even when column density is lower than typical BAL values, very high velocity can dominate the kinetic budget, yielding \(\log L_{KE}\sim 46.5\)–47.2 for the EHVO component [2508.14221].

## 7. Limitations, caveats, and related terminology

SimBAL’s published applications also identify recurrent limitations. Density is frequently weakly constrained unless a density-sensitive transition is present. In the original SDSS J0850+4451 analysis, \(\ion{C}{3}^*\lambda1175\) provided the density constraint needed for the radius estimate [1808.02441]. In J1646 and J0923+0402, the spectra lack such diagnostics, so the density is fixed at \(\log n = 6\) as a modeling assumption rather than being measured directly [2508.14221, 2404.12443]. Radius and energetics in those cases are therefore degenerate with the adopted density.

Global covering fraction \(\Omega\) is another systematic. Many studies adopt \(\Omega=0.2\), motivated by BAL fractions, but explicitly note that mass outflow rates and kinetic luminosities scale with this assumption [1808.02441, 2203.11964, 2404.12443, 2508.14221]. Metallicity and SED choices also matter. In SDSS J0850+4451, a hard SED required substantially larger columns and gave a poorer fit than a soft quasar-like SED, while enhanced metallicity \(Z=3Z_\odot\) improved the fit and lowered the inferred total column [1808.02441]. In J0923+0402, the authors note that higher metallicity would lower the inferred column density and weaken the outflow energetics [2404.12443].

Model flexibility can itself be a source of ambiguity. The FeLoBAL survey introduced spectral PCA for emission-line modeling, priors on PCA coefficients to prevent the continuum model from spuriously fitting away absorption, and modified partial-covering schemes for especially opaque overlapping-trough objects [2203.11964]. In SDSS J1352+4239, SimBAL had to include anomalous reddening and a scattered-light component with a best-fit scattering fraction of about \(29 \pm 0.5\%\) to reproduce the residual flux in the troughs [2001.07347]. Those extensions show both the adaptability of the framework and the dependence of inference on the adopted composite source model.

A separate terminological caveat concerns names. In quasar spectroscopy, “SimBAL” refers to the spectral-synthesis BAL code introduced by Leighly and collaborators [1808.02441]. This should be distinguished from “SimBaL” in low-\(\ell\) CMB analysis, where the term denotes a simulation-based likelihood for Planck temperature-polarization spectra rather than a BAL-quasar spectral-synthesis code [2103.14378]. The shared spelling reflects nomenclature overlap, not methodological continuity.

Overall, SimBAL occupies a specific niche in quasar-outflow research: it is a forward-modeling BAL inference framework optimized for spectra where line blending, partial covering, and non-black saturation render simpler techniques incomplete. The literature to date shows its utility in deriving velocity-resolved physical conditions, testing absorber geometry across wavelength and time, and quantifying outflow energetics from pc to sub-kpc scales across LoBAL, FeLoBAL, and EHVO phenomenology [1808.02441, 1811.04174, 2001.07347, 2203.11964, 2405.06027, 2404.12443, 2508.14221].

Source: https://www.emergentmind.com/topics/simbal