---
title: Locally Optimally Emitting Clouds (LOC)
url: https://www.emergentmind.com/topics/locally-optimally-emitting-clouds-loc
type: topic
---

# Locally Optimally Emitting Clouds (LOC)

to=arxiv_search  天天中彩票不ి  手机版天天中彩票?
{"query":"Locally Optimally Emitting Clouds AGN 1907.06669 1912.03382 1408.4397 1310.6402 2001.03625 0911.1173","max_results":10}
Locally Optimally Emitting Clouds (LOC) is a photoionization framework for active galactic nucleus (AGN) emission-line regions in which the broad-line region (BLR) or narrow-line region (NLR) is not treated as a single-zone cloud population, but as a broad distribution of clouds spanning radius, density, and, equivalently, ionizing flux. Its defining principle is local optimal emission: each spectral line is produced most efficiently by the subset of clouds whose local physical conditions maximize that line’s emissivity, so the observed spectrum is the integrated output of an ensemble rather than the signature of a uniquely tuned cloud. In the literature considered here, LOC is used to model BLR line ratios, EUV diagnostics, broad-line variability, changing-look behavior, Fe II anisotropy, and NLR ionization sequences, and it is also reinterpreted physically in terms of clumpy turbulent outflows and stratified disc winds [1907.06669] [1408.4397] [1912.03382] [2001.03625] [0911.1173] [1310.6402].

## 1. Formal definition and mathematical structure

In LOC modeling, the observable line luminosity is obtained by integrating single-cloud emissivities over a cloud ensemble weighted by distribution functions in radius and density. One implementation writes
\[
L_{\rm line} \propto \int\!\!\int F(r)\, f(r)\, g(n)\, dn\, dr,
\]
with empirical choices \(f(r)\propto r^\Gamma\) and \(g(n)\propto n^\beta\) [1907.06669]. An equivalent BLR formulation integrates over hydrogen density and hydrogen-ionizing photon flux,
\[
\int\!\!\int A\,F(n_H,\Phi_H)\,n_H^{\beta_n}\,\Phi_H^{\beta_\Phi}\, d(\log n_H)\,d(\log \Phi_H),
\]
where \(F(n_H,\Phi_H)\) is the photoionization-predicted line flux and the exponents describe the weighting of clouds across parameter space [1408.4397]. For NLR applications, the same idea is expressed as
\[
L_{\rm line}\propto \iint F(r,n_H)\, f(r)\, g(n_H)\, dr\, dn_H,
\]
with \(f(r)\propto r^\gamma\) and \(g(n_H)\propto n_H^\beta\) [1310.6402].

The conceptual content of these expressions is the same across implementations. LOC assumes that different lines are “optimally emitted” in different restricted regions of the density–radius or density–flux plane, so the integrated spectrum is controlled by the cloud distribution functions and by the emissivity contours of the lines themselves. This is why LOC is a selection-effect model rather than a one-zone prescription: the relevant observable is the sum over the ensemble, not the spectrum of any individual cloud [1310.6402].

Several quantities recur in LOC parameterizations. The hydrogen-ionizing photon flux \(\Phi_H\) is used as a proxy for distance from the ionizing source because \(\Phi_H \propto r^{-2}\) in the optically thin limit [1408.4397]. The ionization parameter is commonly written as
\[
U \equiv \frac{\Phi_H}{c\,n_H},
\]
and in one turbulent-outflow treatment the X-ray ionization parameter is
\[
\xi = \frac{L_X}{n r^2},
\]
with \(U \approx \xi/42\) for the adopted NGC 5548 spectral energy distribution [1912.03382].

## 2. Parameter-space coverage and spectroscopic applications

A central strength of LOC is its ability to reproduce multiple lines that peak at different physical conditions. In EUV BLR spectroscopy, CLOUDY grids were computed over \(7 \le \log n_H\ (\mathrm{cm^{-3}}) \le 14\) and \(17 \le \log \Phi_H\ (\mathrm{cm^{-2}\ s^{-1}}) \le 24\), assuming a one-sided 1D slab, constant density, \(N_H = 10^{23}\ \mathrm{cm^{-2}}\), solar metallicity and solar abundance ratios, no cloud shielding of the continuum source by other clouds, purely thermal motions, and a BLR covering factor of 40% [1408.4397]. In that application, the LOC models provide good fits to the measured fluxes, while the single-component models do not. By contrast, the single-component models could fit only two of the six individual AGN spectra, and only after excluding the problematic N \(\lambda\)991 line; the LOC models gave good fits to four of the six individual AGN and to the EUV subset of the 159-object composite spectrum [1408.4397].

The same study showed why the distributed-cloud approach matters physically. EUV lines such as N IV \(\lambda 765\), O II \(\lambda 833\), and O III \(\lambda 834\) originate primarily from gas with electron temperatures between 37000 K and 55000 K, in BLR clouds with high hydrogen densities (\(n_H \ge 10^{12}\ \mathrm{cm^{-3}}\)) and hydrogen ionizing photon fluxes (\(\Phi_H \ge 10^{22}\ \mathrm{cm^{-2}\ s^{-1}}\)) [1408.4397]. Different lines therefore require access to different regions of the grid. This was one reason single-zone models failed for most objects: even when lines align roughly along constant \(U\), they do not all peak at the same \(U\), and collisionally excited emissivities remain strongly sensitive to temperature, density, and flux [1408.4397].

In the NLR, LOC has been used to fit composite spectra along an AGN ionization sequence. Cloudy v10.0 grids with 7171 cloud models sampled \(17.48 < \log r < 22.0\) and \(2.0 < \log n_H < 8.0\), excluding gas satisfying \(2\log r + \log n_H < 41.5\), corresponding roughly to \(\log U > 0.4\) [1310.6402]. The resulting line-emissivity maps show that different optical diagnostics peak at different radii and densities; for example, \([\mathrm{O\,III}]\,\lambda4363\) peaks at smaller radii than \([\mathrm{O\,III}]\,\lambda5007\), and \([\mathrm{N\,II}]\,\lambda5755\) peaks at smaller radii than \([\mathrm{N\,II}]\,\lambda6584\) [1310.6402]. This is the NLR analogue of the BLR selection effect.

## 3. Physical interpretations: from empirical ensemble to dynamical medium

LOC has often been implemented as a large grid of photoionization calculations without a specific cloud-formation mechanism. A major development is the proposal that the BLR cloud population is the observational imprint of a clumpy turbulent outflow. In that picture, condensations form in thermally unstable zones of an AGN outflow, and the relevant cloud sizes are many orders of magnitude smaller than the global outflow scale \(L_0\), so cloud dynamics can be modeled locally [1912.03382]. The characteristic cooling length is written as
\[
\lambda_{\rm cool} \equiv c_s\, t_{\rm cool},
\]
with the estimate
\[
\lambda_{\rm cool} \approx 3.3\times 10^{10}\, T_5^{3/2}\, n_9^{-1}\, \mathscr{L}_{23}^{-1}\ {\rm cm},
\]
and the essential BLR result is that \(\lambda_{\max} \ll L_\rho\), which justifies local multiphase turbulence simulations [1912.03382].

Using Athena, non-adiabatic hydrodynamics, thermal conduction, heating/cooling, and purely solenoidal forcing, the first 3D local clumpy turbulent outflow simulations for this problem were presented with subsonic turbulent Mach numbers \(M_t \sim 0.05 - 0.75\) and fiducial box size \(L_{\rm box} = 4\,\lambda_{\rm cool}\) [1912.03382]. Those simulations show that condensations form only in a restricted wavenumber interval inside the inertial range,
\[
k_{\min} \le k \le k_F,
\]
where \(k_F\) corresponds to the Field length [1912.03382]. The authors argued that this reproduces the same selection effects characteristic of LOC. A plausible implication is that the broad grid of ionization parameters and densities used in empirical LOC work may arise dynamically from thermal instability plus turbulence rather than being imposed ad hoc.

A related reinterpretation replaces discrete clouds with a continuous, self-shielded, biconical disc wind. In that framework, the wind spans a large portion of the ionizing flux-density \((\phi_H-n_H)\) plane because density varies along streamlines through mass conservation, while ionizing flux is attenuated and reprocessed by the flow itself [2001.03625]. The hydrogen-ionizing photon flux density is
\[
\phi_H = \int_{13.6\mathrm{eV}/h}^{\infty} \frac{4 \pi J_\nu}{h \nu} \, d\nu,
\]
with \(U_H = \phi_H/(n_H c)\), and the authors explicitly stated that the behavior of a line-emitting wind is similar to that of LOC except that the gradients in ionization state and temperature are large-scale and continuous rather than within or between distinct clouds [2001.03625]. They also found that clumpy biconical disc winds can produce BLR-like spectra provided that the wind has a volume filling factor of \(f_V\lesssim0.1\), with the most successful models typically using \(f_V = 0.01\), and that line emission arises almost exclusively from plasma travelling below the escape velocity, implying that “failed winds” are important BLR candidates [2001.03625].

## 4. Responsivity, variability, and broad-line reverberation

LOC has been used to explain why broad Mg II behaves differently from broad Balmer lines in quasars. In one fiducial implementation, the cloud distribution adopted \(\Gamma=-1.1\), \(\beta=-1\), solar abundance, \(N_{\rm H}=10^{23}\,\mathrm{cm^{-2}}\), and a global covering factor of 50%, with a fiducial quasar of \(M_{\rm BH}=10^{8.5}M_\odot\), \(L_{3000\AA}\sim10^{44-45}\,\mathrm{erg\,s^{-1}}\), inner boundary \(R_{\rm in}=10^{16.5}\) cm, and outer boundary \(R_{\rm out}=10^{18}\) cm, or about 0.3 pc [1907.06669]. In that model, Mg II-emitting gas is on average more distant from the ionizing source than the H\(\alpha\)/H\(\beta\) gas, and responds with a lower amplitude to continuum variations [1907.06669].

The physical distinction between Mg II and the Balmer lines is central. H\(\alpha\) and H\(\beta\) are recombination lines, while Mg II is dominated by collisional excitation and has a low excitation energy of 4.4 eV [1907.06669]. For typical BLR conditions, \(T_e\sim10^4\) K and \(n_e\sim10^{10}\,\mathrm{cm^{-3}}\), the Balmer recombination timescale is
\[
\tau_{\rm rec} = (n_e\alpha)^{-1} \sim 0.1\left(10^{10}\,\mathrm{cm^{-3}}/n_e\right)\,\mathrm{hr},
\]
while the Mg II collisional timescale is \(\sim 0.01\) s [1907.06669]. The line response is therefore not limited by atomic timescales; it is set by where in the BLR the line forms. One-cloud calculations showed that Mg II emission becomes significant only at large enough column density, peaking around \(N_{\rm H}\sim10^{23}\,\mathrm{cm^{-2}}\), and dropping again at very high columns where clouds become too optically thick [1907.06669].

Responsivity was formalized as
\[
\eta=\frac{d\log F(r)}{d\log \Phi(H)}\propto -0.5\,\frac{d\log F(r)}{d\log r},
\]
so that if \(F(r)\propto r^\gamma\), then \(\eta=-\gamma/2\) [1907.06669]. In the fiducial model, when the continuum drops by 1 dex, Mg II luminosity falls by only \(\sim0.45\) dex, compared with about 0.6 dex for H\(\alpha\) and 0.7 dex for H\(\beta\) [1907.06669]. This lower responsivity, combined with a larger average formation radius, naturally dilutes and slows Mg II variability. If the BLR is truncated at \(\sim 0.3\) pc, most of the Mg II flux is emitted near that outer boundary, so the line does not display strong breathing; depending on \(R_{\rm out}\), the same LOC framework can produce fully breathing, partially breathing, or no-breathing behavior [1907.06669].

These results have direct implications for reverberation mapping. Broad Mg II lags are intrinsically harder to measure because the line varies less strongly, and because the Mg II-emitting region is somewhat farther out, longer monitoring baselines are needed [1907.06669]. Even so, a measured Mg II lag can still be used to infer a BLR size and black hole mass through
\[
M_{\rm BH} = \frac{f\,\Delta V^2 R}{G}.
\]
The same model suggests, however, that Mg II may not have a strong intrinsic size-luminosity relation for an individual quasar, because its emission can remain tied to the outer truncation radius rather than tracking luminosity in the same way as H\(\beta\) [1907.06669].

Changing-look behavior is interpreted similarly. Simulations of a continuum decline from \(L_{3000\AA}=10^{45}\) to \(10^{43}\,\mathrm{erg\,s^{-1}}\) showed Balmer lines fading first while Mg II remains detectable over a wider luminosity range [1907.06669]. This provides a natural LOC explanation for the persistence of broad Mg II in changing-look quasars defined on H\(\alpha\)/H\(\beta\), and for the rare population of broad Mg II emitters in otherwise normal galaxy spectra [1907.06669].

## 5. Low-ionization clouds, column density, and anisotropic emission

LOC has also been extended by adding dynamical filtering to the low-ionization part of parameter space. For Fe II-emitting gas in quasars, the relevant clouds occupy the low-ionization region of the ionizing flux–density plane, but the decisive parameter is column density because radiation pressure can expel low-column clouds [0911.1173]. The force multiplier, defined as the ratio of total gas opacity to electron scattering opacity, is \(\sim 10^3 - 10^4\) in Fe II-emitting gas [0911.1173]. Since observed systems typically have \(L/L_{Edd}\sim 10^{-1}\), whereas optically thin low-ionization gas would require \(L/L_{Edd} \sim 10^{-3} - 10^{-4}\) for infall, the conclusion is that infalling Fe II clouds must be sufficiently thick that radiation acts primarily on an illuminated surface layer rather than the entire cloud [0911.1173].

The minimum column density for infall is written as
\[
N\left( {\rm{H}} \right)_{\rm infall}  \ge {f_{\rm thick}\over \sigma_T } \frac{L}{L_{Edd}}
\simeq 1.5 \times 10^{24}f_{\rm thick}\frac{L}{L_{Edd}}{\rm\,cm^{ - 2}},
\]
with \(f_{\rm thick}\sim 1\) [0911.1173]. In the standard cloud, the ionized surface layer has \(N(\mathrm{H}) \approx 5 \times 10^{21}\) before the ionization front is reached, and the cloud becomes neutral and Fe II-producing deeper inside [0911.1173]. This produces a physically filtered LOC ensemble: clouds with too small a column are pushed out, while only the high-column subset contributes to the infalling Fe II component.

Anisotropy is then unavoidable. The observed spectrum is dominated by the shielded face of the infalling clouds rather than a symmetric distribution of emitters [0911.1173]. In this geometry, optical Fe II emission is nearly isotropic, whereas UV Fe II is predominantly inwardly beamed. Once the column exceeds the ionization-front threshold, around \(N(\mathrm{H}) > 10^{21.7}\), Fe II emission becomes strongly inwardly beamed; as the cloud column increases further into the \(10^{22}-10^{23}\) range, the Fe II/H\(\beta\) ratio increases, approaching an asymptotic value of about 3 for sufficiently large \(N(\mathrm{H})\) or \(L/L_{Edd}\) [0911.1173]. The paper interpreted this as a physical driver for Eigenvector 1: cloud column density acts as the hidden variable coupling \(L/L_{\rm Edd}\), Fe II strength, and spectral differences [0911.1173].

Within LOC language, this is a notable conceptual shift. The model preserves the idea that different lines arise from locally optimal conditions, but it adds a dynamical survival criterion to the cloud ensemble. This suggests that, at least for low-ionization species, the cloud distribution sampled by the spectrum may be shaped jointly by photoionization optimality and radiation-pressure selection.

## 6. NLR sequences, unresolved issues, and nomenclature

In the NLR, LOC has been used to interpret systematic changes in ionization state across AGN composite spectra. After integrating over a wide range of radii and densities, the models indicate that the radial extent of the NLR is the major parameter determining the position of high- to moderate-ionization AGN along the sequence [1310.6402]. The best dust-free model kept the density weighting nearly fixed at \(\beta = -1.4\) while varying the radial weighting \(\gamma\) from \(-0.75\) at the high-ionization end through \(-0.5\), \(-0.25\), and \(0.0\) to \(1.0\) at the low-ionization end [1310.6402]. The physical interpretation was explicit: higher-ionization AGN contain optimally emitting clouds more concentrated toward the central continuum source, while lower-ionization AGN have more extended NLRs [1310.6402].

The same analysis inferred that ionizing luminosity is anticorrelated with NLR ionization level for the AGN sequence selected by mean-field independent component analysis, and, for a fiducial covering factor \(\Omega = 0.4\), characteristic sizes up to roughly 7 kpc were obtained for the lowest-ionization subset [1310.6402]. The authors also considered whether the sequence might instead be a mixing curve of star formation and AGN emission, but argued that while many galaxies do show such composite behavior, the AGN sequence they isolated appears to be a special set of objects with negligible star-formation excitation [1310.6402]. A possible interpretation deserving further exploration was that the ionization sequence might be an age sequence in which lower-ionization objects are older and have systematically cleared out their central regions by radiation pressure [1310.6402].

LOC does not eliminate all difficulties. In the NLR application, the observed increases of \([\mathrm{O\,III}]\,4363/5007\) and \([\mathrm{N\,II}]\,5755/6584\) toward lower ionization were not reproduced, and the temperature-sensitive auroral-line ratios remained a significant unresolved problem [1310.6402]. In the BLR EUV application, N \(\lambda\)991 was systematically underpredicted, with the most likely explanation being a combination of line-fitting and deblending uncertainties rather than a definitive abundance anomaly [1408.4397]. These shortcomings delimit the explanatory range of current LOC implementations without undermining the core selection-effect picture.

A separate source of confusion is terminological rather than physical. The acronym “LOC” is also used for a line radiative transfer program, explicitly “line transfer with OpenCL,” which is a deterministic non-LTE ray-tracing code for 1D and 3D interstellar-medium line transfer and is not related to the astrophysical concept of Locally Optimally Emitting Clouds [2009.12609]. In AGN emission-line theory, by contrast, LOC denotes the cloud-ensemble framework summarized above.

Source: https://www.emergentmind.com/topics/locally-optimally-emitting-clouds-loc