Quasar Outflows Research Articles

Radiation-pressure Waves and Multiphase Quasar Outflows

Abstract We report on quasar outflow properties revealed by analyzing more than 60 composite outflow spectra built from ∼60,000 C iv absorption troughs in the SDSS-III/BOSS DR12QBAL catalog. We assess the dependences of the equivalent widths of many outflow metal absorption features on outflow velocity, trough width and position, and quasar magnitude and redshift. The evolution of the equivalent widths of the O vi and N v lines with outflow velocity correlates with that of the mean absorption-line width, the outflow electron density, and the strength of lines arising from collisionally excited metastable states. None of these correlations are found for the other high- or low-ionization species, and different behaviors with trough width are also suggested. We find no dependence on quasar magnitude or redshift in any case. All the observed trends can be reconciled by considering a multiphase stratified outflow structure, where inner regions are colder, denser, and host lower-ionization species. Given the prevalence of radiative acceleration in quasar outflows found by Mas-Ribas & Mauland, we suggest that radiation pressure sweeps up and compresses the outflowing gas outward, creating waves or filaments where the multiphase stratified structure could take form. This scenario is supported by the suggested correlation between electron density and outflow velocity, as well as by the similar behavior observed for the line and line-locking components of the absorption features. We show that this outflow structure is also consistent with other X-ray, radiative transfer, and polarization results, and discuss the implications of our findings for future observational and numerical quasar outflow studies.

AGN radiative feedback in the early growth of massive black holes

ABSTRACT Growing observational evidence confirms the existence of massive black holes ($M_{\rm BH} \sim 10^9 \, \mathrm{M}_{\odot }$), accreting at rates close to the Eddington limit, at very high redshifts ($z \gtrsim 6\!-\!7$) in the early Universe. Recent observations indicate that the host galaxies of the first quasars are chemically evolved systems, containing unexpectedly large amounts of dust. Such a combination of high luminosities and large dust content should form favourable physical conditions for radiative dusty feedback. We explore the impact of the active galactic nucleus (AGN) feedback, driven by radiation pressure on dust, on the early growth of massive black holes. Assuming Eddington-limited exponential black hole growth, we find that the dynamics and energetics of the radiation pressure-driven outflows also follow exponential trends at late times. We obtain modest outflow energetics (with momentum flux $\dot{p} \lesssim L/c$ and kinetic power $\dot{E}_{\rm k} \lesssim 10^{-3} L$), comparable with available observations of quasar-driven outflows at very high redshifts, but significantly lower than typically observed in local quasars and predicted by wind energy-driven models. AGN radiative dusty feedback may thus play an important role in powering galactic outflows in the first quasars in the early Universe.

Ionization Mechanisms in Quasar Outflows

Abstract The various ionization mechanisms at play in active galactic nuclei (AGN) and quasars have been well studied, but relatively little has been done to separately investigate the contributions of these ionization mechanisms within the host galaxy and outflowing components. Using Gemini integral field spectroscopy data, we study the ionization properties of these two components in four nearby (z ≲ 0.2) radio-quiet Type 1 quasars. Emission-line ratios and widths are employed to identify the dominant ionization mechanisms for the host and outflow components in each object. We find that photoionization by the AGN often dominates the ionization of both gaseous components in these systems. In three cases, the outflowing gas is more highly ionized than the gas in the host, indicating that it is more strongly exposed to the ionizing radiation field of the AGN. In two objects, a positive correlation between the line widths and line ratios in the outflowing gas component indicates that shocks with velocities of order 100–500 km s−1 may also be contributing to the ionization and heating of the outflowing gas component. The line ratios in the outflowing gas of one of these two objects also suggest a significant contribution from photoionization by hot, young stars in the portion of the outflow that is closest to star-forming regions in the host galaxy component. The data thus favor photoionization by hot stars in the host galaxy rather than stars formed in the outflow itself.

A high signal-to-noise HST spectrum towards J1009+0713: precise absorption measurements in the CGM of two galaxies

ABSTRACT High signal-to-noise spectra towards background quasars are crucial for uncovering weak absorption in the circumgalactic medium (CGM) of intervening galaxies, such as the diagnostic lines of N v that provide insight to the ionization process of warm gas but typically have low equivalent widths. We present a new spectrum from the Hubble Space Telescope with a signal-to-noise ratio of ∼20–35 towards the quasar SDSS J1009+0713 and analyse absorption systems in the CGM of two L⋆ galaxies close to the line of sight. We identify additional absorption in the CGM of these galaxies that was not reported by the previous lower signal-to-noise spectrum, as well as Milky Way absorbers and quasar outflows from J1009+0713. We measure log (NN v/NO vi) ∼ −1.1 for two CGM absorbers, inconsistent with gas in collisional ionization equilibrium and consistent with a radiatively cooling bulk flow of ∼50–150 km s−1, which could be produced by galactic winds. These column density ratios are also consistent with those found for other L⋆ galaxies and for some gas in the Milky Way’s halo. We place upper limits of log (NN v/NO vi) < −1.8 to −1.2 for other O vi absorbers in the same haloes, which suggests that O vi is produced by different processes in different parts of the CGM, even within the same galactic halo. Together with the kinematically different structure of high- and low-ionization lines, these results indicate there are many components to a single galaxy’s gaseous halo. We find the redshift number density of Ly α forest absorbers and broad Ly α absorbers are consistent with expectations at this redshift.

ERQs are the BOSS of quasar samples: the highest velocity [O iii] quasar outflows

ABSTRACT We investigate extremely red quasars (ERQs), a remarkable population of heavily reddened quasars at redshift z ∼ 2−3 that might be caught during a short-lived ‘blow-out’ phase of quasar/galaxy evolution. We perform a near-IR observational campaign using Keck/NIRSPEC, VLT/X-shooter, and Gemini/GNIRS to measure rest-frame optical spectra of 28 ERQs with median infrared luminosity 〈log L(erg s−1)〉 ∼ 46.2. They exhibit the broadest and most blueshifted [O iii] λ4959,5007 emission lines ever reported, with widths (w90) ranging between 2053 and 7227 km s−1, and maximum outflow speeds (v98) up to 6702 km s−1. ERQs on average have [O iii] outflows velocities about three times larger than those of luminosity-matched blue quasar samples. This discrepancy can be explained by a strong correlation between [O iii] kinematics and i–W3 colour, and not by radio loudness, or higher Eddington ratios. We estimate for these objects that at least 3–5 per cent of their bolometric luminosity is being converted into the kinetic power of the observed wind. Our results reveal that ERQs have the potential to strongly affect the evolution of host galaxies.

VLT/X-Shooter Survey of BAL Quasars: Large Distance Scale and AGN Feedback

Abstract We conducted a survey of quasar outflows using the Very Large Telescope/X-Shooter spectrograph. When choosing the 14 broad absorption line (BAL) and mini-BALs comprising this sample, the data did not cover the S iv and S iv* troughs, whose ratio can be used to determine the distance of the outflows from the central source (R). Therefore, this “Blind Survey” is unbiased toward a particular distance scale. Out of the eight outflows where R can be measured, six have R > 100 pc (spanning the range 100–4500 pc), one has R > 10 pc, and only one (at R < 40 pc) is compatible with a much smaller R scale. At least two of the outflows have a kinetic luminosity greater than 0.5% of their Eddington luminosity, implying that they are able to provide significant active galactic nucleus (AGN) feedback. The outflows span a range of 0 to −10,000 km s−1 in velocity, total column density between 1020.1 and 1022.6 cm−2, ionization parameter ( ) in the range 0.01–1, and electron number density between 103.2 and 105.4 cm−3, with one upper and one lower limit. The results of this survey can be extrapolated to the majority of BAL outflows, implying that most of these outflows are situated far away from the AGN accretion disk, and that a significant portion of them can contribute to AGN feedback processes.

The properties of broad absorption line outflows based on a large sample of quasars

Quasar outflows carry mass, momentum and energy into the surrounding environment, and have long been considered a potential key factor in regulating the growth of supermassive black holes and the evolution of their host galaxies. A crucial parameter for understanding the origin of these outflows and measuring their influence on their host galaxies is the distance (R) between the outflow gas and the galaxy center. While R has been measured in a number of individual galaxies, its distribution remains unknown. Here we report the distributions of R and the kinetic luminosities of quasars outflows, using the statistical properties of broad absorption line variability in a sample of 915 quasars from the Sloan Digital Sky Surveys. The mean and standard deviation of the distribution of R are 10^{1.4+/-0.5} parsecs. The typical outflow distance in this sample is tens of parsec, which is beyond the theoretically predicted location (0.01 ~ 0.1 parsecs) where the accretion disc line-driven wind is launched, but is smaller than the scales of most outflows that are derived using the excited state absorption lines. The typical value of the mass-flow rate is of tens to a hundred solar masses per year, or several times the accretion rate. The typical kinetic-to-bolometric luminosity ratio is a few per cent, indicating that outflows are energetic enough to influence the evolution of their host galaxies.

Direct Measurement of Quasar Outflow Wind Acceleration* † ‡

Abstract We search for velocity changes (i.e., acceleration/deceleration) of narrow absorption lines (NALs) that are intrinsic to the quasars using spectra of six bright quasars that have been observed more than once with 8–10 m class telescopes. While variations in line strength and profile are frequently reported (especially in broader absorption lines), definitive evidence for velocity shifts has not been found, with only a few exceptions. Direct velocity shift measurements are valuable constraints on the acceleration mechanisms. In this study, we determine velocity shifts by comparing the absorption profiles of NALs at two epochs separated by more than 10 yr in the observed frame using the cross-correlation function method, and we estimate the uncertainties using Monte Carlo simulations. We do not detect any significant shifts, but we obtain 3σ upper limits on the acceleration of intrinsic NALs (compared to intervening NALs in the same quasars) of ∼0.7 km s−1 yr−1 (∼0.002 cm s−2) on average. We discuss possible scenarios for nondetection of NAL acceleration/deceleration and examine the resulting constraints on the physical conditions in accretion disk winds.

Galaxy and Quasar Fueling Caught in the Act from the Intragroup to the Interstellar Medium

Abstract We report the discovery of six spatially extended (10–100 kpc) line-emitting nebulae in the z ≈ 0.57 galaxy group hosting PKS 0405−123, one of the most luminous quasars at z < 1. The discovery is enabled by the Multi Unit Spectroscopic Explorer and provides tantalizing evidence connecting large-scale gas streams with nuclear activity on scales of <10 proper kpc (pkpc). One of the nebulae exhibits a narrow, filamentary morphology extending over 50 pkpc toward the quasar with narrow internal velocity dispersion (50 ) and is not associated with any detected galaxies, consistent with a cool intragroup medium filament. Two of the nebulae are 10 pkpc north and south of the quasar with tidal-arm–like morphologies. These two nebulae, along with a continuum-emitting arm extending 60 pkpc from the quasar, are signatures of interactions that are expected to redistribute angular momentum in the host interstellar medium (ISM) to facilitate star formation and quasar fueling in the nucleus. The three remaining nebulae are among the largest and most luminous [O iii] emitting “blobs” known (1400–2400 pkpc2) and correspond both kinematically and morphologically to interacting galaxy pairs in the quasar host group, consistent with arising from stripped ISM rather than large-scale quasar outflows. The presence of these large- and small-scale nebulae in the vicinity of a luminous quasar bears significantly on the effect of large-scale environment on galaxy and black hole fueling, providing a natural explanation for the previously known correlation between quasar luminosity and cool circumgalactic medium.

On the structure and energetics of quasar broad absorption-line outflows

Quasar accretion-disk outflows might play an important role in galaxy evolution, but they are notoriously difficult to study due to line saturation and blending problems in the Ly-alpha forest. We circumvent these problems by constructing median composite spectra of diverse broad absorption lines (BALs) and `mini-BALs' in SDSS-III BOSS quasars at redshifts 2.3 < z < 3.5. Sorting by CIV 1549,1551 absorption-line strength with AlIII 1855,1863 as an additional indicator of low ionisations (LoBALs) we find the following: (1) Deeper and broader BALs are accompanied by weaker HeII 1640 emission lines, consistent with softer ionising spectra producing more effective radiative acceleration. (2) PV 1118,1128 absorption is present with resolved ~1:1 depth ratios in all composites from mini-BALs to strong BALs indicating that line saturation, large total column densities log N_H(cm^-2) > 22.7, and large ionisation parameters log U > -0.5 are common. (3) Different observed depths in saturated lines identify inhomogeneous partial covering on spatial scales <0.006 pc, where weak/low-abundance transitions like PV form in small high-column density clumps while stronger/broader lines like CIV form in larger volumes. (4) Excited-state SiIV* 1073 and CIII* 1176 lines in BAL outflows indicate typical densities n_e > 3 x 10^5 cm^-3 and maximum radial distances R < 23pc from the quasars. (5) For reasonable actual distances, the median BAL outflow has minimum kinetic energy L_K/L > 0.005(R/1.2pc) sufficient (by some estimates) for feedback to galaxy evolution. (6) LoBAL quasars have the largest median outflow column densities, highest velocities, and weakest HeII 1640 emission in our study; they appear to be at one extreme in a distribution of quasar properties where softer ionising spectra drive more powerful outflows.

Distance, Energy, and Variability of Quasar Outflows: Two HST/COS Epochs of LBQS 1206+1052∗

Abstract We analyze new HST/COS spectra for two quasar absorption outflows seen in the quasi-stellar object LBQS 1206+1052. These data cover, for the first time, absorption troughs from S iv, Si ii, and P v. From the ratio of the S iv* to S iv column densities, we measure the electron number density of the higher-velocity (−1400 km s−1, v1400) outflow to be cm−3 and constrain the lower-velocity (−730 km s−1, v700) outflow to log(n e ) &gt; 5.3 cm−3. The n e associated with the higher-velocity outflow is an order of magnitude larger than reported in prior work. We find that the previous measurement was unreliable since it was based on density-sensitive absorption troughs that were likely saturated. Using photoionization models, we determine the best χ 2-minimization fit for the ionization parameter and hydrogen column density of the higher-velocity outflow: log( and log( cm−2, respectively. We calculate from U H and n e a distance of pc from the central source to the outflow. Using an SED attenuated by the v700 outflow yields a two-phase photoionization solution for the v1400 outflow, separated by a . Otherwise, the resultant distance, mass flux, and kinetic luminosity are similar to the unattenuated case. However, the attenuated analysis has significant uncertainties due to a lack of constraints on the v700 outflow in 2017.

Prevalence of radio jets associated with quasar outflows and feedback

AbstractWe have identified that radio jets are commonly associated with “radiative mode” feedback in quasars. By performing a systematic multi-wavelength study of z &lt; 0.2 quasars, we have found that 70–80% of our sample of ‘radio-quiet’ type 2 quasars, which host kpc-scale ionized gas outflows, exhibit radio jet structures. Here, we discuss our results on the pilot sample of 10 objects that combine high resolution (∼ 0.25 - 1 arcsec) radio imaging at 1-7 GHz with optical IFU observations. Our results demonstrate that it is extremely common for jets to be spatially and kinematically linked to kpc-scale ionized gas kinematics in such quasars. Therefore, radio jets may be an important driver of outflows during ‘radiative mode’ feedback, apparently blurring the lines between the traditional divisions of feedback modes.

Emergence and Variability of Broad Absorption Line Quasar Outflows

Abstract We isolate a set of quasars that exhibit emergent C iv broad absorption lines (BALs) in their spectra by comparing spectra in the Sloan Digital Sky Survey (SDSS) Data Release 7 and the SDSS/BOSS Data Releases 9 and 10. After visually defining a set of emergent BALs, follow-up observations were obtained with the Gemini Observatory for 105 quasars. We find an emergence rate consistent with the previously reported disappearance rate of BAL quasars given the relative numbers of non-BAL and BAL quasars in the SDSS. We find that candidate newly emerged BALs are preferentially drawn from among BALs with smaller balnicity indices, shallower depths, larger velocities, and smaller widths. Within two rest-frame years (average) after a BAL has emerged, we find it equally likely to continue increasing in equivalent width in an observation 6 months later (average) as it is to start decreasing. From the time separations between our observations, we conclude that the coherence timescale of BALs is less than 100 rest-frame days. We observe coordinated variability among pairs of troughs in the same quasar, likely due to clouds at different velocities responding to the same changes in ionizing flux, and the coordination is stronger if the velocity separation between the two troughs is smaller. We speculate that the latter effect may be due to clouds having on average lower densities at higher velocities owing to mass conservation in an accelerating flow, causing the absorbing gas in those clouds to respond on different timescales to the same ionizing flux variations.

Quenching star formation with quasar outflows launched by trapped IR radiation

We present cosmological radiation-hydrodynamic simulations, performed with the code Ramses-RT, of radiatively-driven outflows in a massive quasar host halo at $z = 6$. Our simulations include both single- and multi-scattered radiation pressure on dust from a quasar and are compared against simulations performed with thermal feedback. For radiation pressure-driving, we show that there is a critical quasar luminosity above which a galactic outflow is launched, set by the equilibrium of gravitational and radiation forces. While this critical luminosity is unrealistically high in the single-scattering limit for plausible black hole masses, it is in line with a $\approx 3 \times 10^9 \, \rm M_\odot$ black hole accreting at its Eddington limit, if infrared (IR) multi-scattering radiation pressure is included. The outflows are fast ($v \, \gtrsim \, 1000 \, \rm km \, s^{-1}$) and strongly mass-loaded with peak mass outflow rates $\approx 10^3 - 10^4 \, \rm M_\odot \, yr^{-1}$, but short-lived ($< 10 \, \rm Myr$). Outflowing material is multi-phase, though predominantly composed of cool gas, forming via a thermal instability in the shocked swept-up component. Radiation pressure- and thermally-driven outflows both affect their host galaxies significantly, but in different, complementary ways. Thermally-driven outflows couple more efficiently to diffuse halo gas, generating more powerful, hotter and more volume-filling outflows. IR radiation, through its ability to penetrate dense gas via diffusion, is more efficient at ejecting gas from the bulge. The combination of gas ejection through outflows with internal pressurisation by trapped IR radiation leads to a complete shut down of star formation in the bulge. We hence argue that radiation pressure-driven feedback may be an important ingredient in regulating star formation in compact starbursts, especially during the quasar's `obscured' phase.

Evidence that 50% of BALQSO Outflows Are Situated at Least 100 pc from the Central Source

Abstract The most robust way for determining the distance of quasar absorption outflows is the use of troughs from ionic excited states. The column density ratio between the excited and resonance states yields the outflow number density. Combined with a knowledge of the outflow’s ionization parameter, a distance from the central source (R) can be determined. Here we report results from two surveys targeting outflows that show troughs from S iv. One survey includes 1091 SDSS and BOSS quasar spectra, and the other includes higher-quality spectra of 13 quasars observed with the Very Large Telescope. Our S iv samples include 38 broad absorption line (BAL) outflows and four mini-BAL outflows. The S iv is formed in the same physical region of the outflow as the canonical outflow-identifying species C iv. Our results show that S iv absorption is only detected in 25% of C iv BAL outflows. The smaller detection fraction is due to the higher total column density (N H) needed to detect S iv absorption. Since R empirically anticorrelates with N H, the results of these surveys can be extrapolated to C iv quasar outflows with lower N H as well. We find that at least 50% of quasar outflows are at distances larger than 100 pc from the central source, and at least 12% are at distances larger than 1000 pc. These results have profound implications for the study of the origin and acceleration mechanism of quasar outflows and their effects on the host galaxy.

Study of Outflows in Luminous Quasars at Redshift ∼0.4–0.8

Abstract We perform a systematic study of outflow in the narrow-line region (NLR) of active galactic nuclei (AGNs) at z ∼ 0.4–0.8 based upon a large sample of ∼900 quasars at z ∼ 0.4–0.8. The sample is extracted from the Sloan Digital Sky Survey by mainly requiring that (1) the g-band magnitude is brighter than 19 mag and (2) the [O iii]λ5007 emission line has a signal-to-noise ratio larger than 30. Profiles of multiple emission lines are modeled by a sum of several Gaussian functions. The spectral analysis allows us to identify (1) a prevalence of both [O iii]λ5007 line blue asymmetry and bulk velocity blueshift of both [Ne iii]λ3869 and [Ne iii]λ3426 lines, when the [O ii]λ3727 line is used as a reference. The velocity offset of the [O iii]λ5007 line is, however, distributed around a value of zero, except for a few outliers. (2) Not only the significant [O iii]λ5007 line asymmetry, but also the large bulk velocity offsets of [Ne iii]λ3869 and [Ne v]λ3426 emission lines tend to occur in the objects with high L/L Edd, which is considerably consistent with the conclusions based on local AGNs. With three M BH estimation methods, the significance level of the trend is found to be better than 2.9σ, 3.2σ, and 1.8σ for [O iii], [Ne iii], and [Ne v], respectively. After excluding the role of radio jets, the revealed dependence of NLR gas outflow on L/L Edd allows us to argue that the pressure caused by the wind/radiation launched/emitted from a central supermassive black hole is the most likely origin of the outflow in these distant quasars, which implies that the outflow in luminous AGNs up to z ∼ 1 have the same origin.

Quasar outflow energetics from broad absorption line variability

Quasar outflows have long been recognized as potential contributors to the co-evolution between supermassive black holes (SMBHs) and their host galaxies. The role of outflows in AGN feedback processes can be better understood by placing observational constraints on wind locations and kinetic energies. We utilize broad absorption line (BAL) variability to investigate the properties of a sample of 71 BAL quasars with P$\thinspace$V broad absorption. The presence of P$\thinspace$V BALs indicates that other BALs like C$\thinspace$IV are saturated, such that variability in those lines favours clouds crossing the line of sight. We use these constraints with measurements of BAL variability to estimate outflow locations and energetics. Our data set consists of multiple-epoch spectra from the Sloan Digital Sky Survey and MDM Observatory. We detect significant (4$\sigma$) BAL variations from 10 quasars in our sample over rest frame time-scales between < 0.2-3.8 yr. Our derived distances for the 10 variable outflows are nominally < 1-10 pc from the SMBH using the transverse-motion scenario, and < 100-1000 pc from the central source using ionization-change considerations. These distances, in combination with the estimated high outflow column densities (i.e. $N_{\textrm{H}}$ > 10$^{22}$ cm$^{-2}$), yield outflow kinetic luminosities between ~ 0.001-1 times the bolometric luminosity of the quasar, indicating that many absorber energies within our sample are viable for AGN feedback.

Ionized Gas Outflows from the MAGNUM Survey: NGC 1365 and NGC 4945

AGN feedback, acting through strong outflows accelerated in the nuclear region of AGN hosts, is invoked as a key ingredient for galaxy evolution by many models to explain the observed BH-galaxy scaling relations. Recently, some direct observational evidence of radiative mode feedback in action has been finally found in quasars at $z$>1.5. However, it is not possible to study outflows in quasars at those redshifts on small scales ($\lesssim$100 pc), as spatial information is limited by angular resolution. This is instead feasible in nearby active galaxies, which are ideal laboratories to explore outflow structure and properties, as well as the effects of AGN on their host galaxies. In this proceeding we present preliminary results from the MAGNUM survey, which comprises nearby Seyfert galaxies observed with the integral field spectrograph VLT/MUSE. We focus on two sources, NGC 1365 and NGC 4945, that exhibit double conical outflows extending on distances >1 kpc. We disentangle the dominant contributions to ionization of the various gas components observed in the central $\sim$5.3 kpc of NGC 1365. An attempt to infer outflow 3D structure in NGC 4945 is made via simple kinematic modeling, suggesting a hollow cone geometry.

The origin of fast molecular outflows in quasars: molecule formation in AGN-driven galactic winds

We explore the origin of fast molecular outflows that have been observed in Active Galactic Nuclei (AGN). Previous numerical studies have shown that it is difficult to create such an outflow by accelerating existing molecular clouds in the host galaxy, as the clouds will be destroyed before they can reach the high velocities that are observed. In this work, we consider an alternative scenario where molecules form in-situ within the AGN outflow. We present a series of hydro-chemical simulations of an isotropic AGN wind interacting with a uniform medium. We follow the time-dependent chemistry of 157 species, including 20 molecules, to determine whether molecules can form rapidly enough to produce the observed molecular outflows. We find H$_2$ outflow rates up to 140 M$_\odot$ yr$^{-1}$, which is sensitive to density, AGN luminosity, and metallicity. We compute emission and absorption lines of CO, OH and warm (a few hundred K) H$_2$ from the simulations in post-processing. The CO-derived outflow rates and OH absorption strengths at solar metallicity agree with observations, although the maximum line of sight velocities from the model CO spectra are a factor $\approx$2 lower than is observed. We derive a CO (1-0) to H$_2$ conversion factor of $\alpha_{\rm{CO} (1-0)}$ = 0.13 M$_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, 6 times lower than is commonly assumed in observations of such systems. We find strong emission from the mid-infrared lines of H$_2$. The mass of H$_2$ traced by this infrared emission is within a few per cent of the total H$_2$ mass. This H$_2$ emission may be observable by JWST.

Narrow absorption lines complex I: one form of broad absorption line

We discover that some of the broad absorption lines (BALs) are actually a complex of narrow absorption lines (NALs). As a pilot study of this type of BAL, we show this discovery through a typical example in this paper. Utilizing the two-epoch observations of J002710.06-094435.3 (hereafter J0027-0944) from the Sloan Digital Sky Survey (SDSS), we find that each of the C IV and Si IV BAL troughs contains at least four NAL doublets. By resolving the Si IV BAL into multiple NALs, we present the following main results and conclusions. First, all these NALs show coordinated variations between the two-epoch SDSS observations, suggesting that they all originate in the quasar outflow, and that their variations are due to global changes in the ionization condition of the absorbing gas. Secondly, a BAL consisting of a number of NAL components indicates that this type of BAL is basically the same as the intrinsic NAL, which tends to support the inclination model rather than the evolution model. Thirdly, although both the C IV and Si IV BALs originate from the same clumpy substructures of the outflow, they show different profile shapes: multiple absorption troughs for the Si IV BAL in a wider velocity range, while P-Cygni for the C IV BAL in a narrower velocity range. This can be interpreted by the substantial differences in fine structure and oscillator strength between the Si IV {\lambda}{\lambda}1393, 1402 and C IV {\lambda}{\lambda}1548, 1551 doublets. Based on the above conclusions, we consider that the decomposition of a BAL into NALs can serve as a way to resolve the clumpy structure for outflows, and it can be used to learn more about characteristics of the clumpy structure and to test the outflow model, when utilizing high-resolution spectra and photoionization model.