The physical constraints on a new LoBAL QSO at z=4.82

TL;DR

This study characterizes a high-redshift LoBAL QSO at z=4.82, revealing strong outflows and high Eddington ratio, providing insights into early quasar feedback mechanisms and black hole growth.

Contribution

It provides the first detailed physical analysis of a high-redshift LoBAL QSO, including black hole mass, outflow properties, and feedback implications.

Findings

Confirmed the LoBAL nature with new redshift at 4.82

Estimated black hole mass of approximately 2.3 billion solar masses

Derived significant outflow kinetic power exceeding 0.9% of Eddington luminosity

Abstract

Very few low-ionization broad absorption line (LoBAL) QSOs have been found at high redshifts to date. One high-redshift LoBAL QSO, J0122+1216, was recently discovered at the Lijiang 2.4-m Telescope with an initial redshift determination of 4.76. Aiming to investigate its physical properties, we carried out follow-up observations in the optical and near-IR spectroscopy. Near-IR spectra from UKIRT and P200 confirms that it is a LoBAL, with a new redshift determination of $4.82\pm0.01$ based on the \mgii~ emission-line. The new \mgii~ redshift determination reveals strong blueshifts and asymmetry of the high-ionization emission lines. We estimated a black hole mass of $\sim 2.3\times 10^9 M_\odot$ and Eddington ratio of $\sim 1.0$ according to the empirical \mgii-based single-epoch relation and bolometric correction factor. It is possible that strong outflows are the result of an extreme quasar environment driven by the high Eddington ratio. A lower limit on the outflowing kinetic power ($>0.9\% L_{Edd}$) was derived from both emission and absorption lines, indicating these outflows play a significant role in the feedback process to regulate the growth of its black hole as well as host galaxy evolution.

Tables

Table 1. Table 1 Near-IR and optical monitoring observations

Instrument	$\lambda /\Delta \lambda$	Slit	Integration	Observation
Name		Width	Time	Date
		(arcsec)	(seconds)	(UT)
TNG/NICS	500	1.0	800$\times$2	2015-02-13
UKIRT/UIST	350	0.8	240$\times$20	2015-10-30
LJT/YFOSC	320	1.8/2.5	2400-2700	2014 to 2016
P200/TripleSpec	2700	1.0	300$\times$8	2017-02-09

Table 2. Table 2 Magnitudes from multiple surveys

Name	J0122+1216	Reference star
Ra (J2000)	01:22:47.34	01:22:57.8
Dec(J2000)	+12:16:24.0	+12:15:42.0
SDSS-u	23.35 $\pm$0.62	18.88$\pm$0.02
SDSS-g	24.23$\pm$0.41	17.47$\pm$0.01
SDSS-r	22.27$\pm$0.16	16.91$\pm$0.01
SDSS-i	19.36$\pm$0.02	16.73$\pm$0.01
SDSS-z	19.24$\pm$0.06	16.64$\pm$0.01
Y	18.484$\pm$0.036	16.039$\pm$0.006
J	17.979$\pm$0.043	15.751$\pm$0.006
H	17.183$\pm$0.042	15.326$\pm$0.008
K	16.604$\pm$0.042	15.283$\pm$0.012
W1	15.58$\pm$0.053	15.231$\pm$0.044
W2	15.03$\pm$0.095	15.321$\pm$0.121
W3	11.50$\pm$0.165	12.701
W4	8.62	9.184

Table 3. Table 3 Optical monitoring observations

Obs	Ly$\alpha$/Si iv	BI(C iv)	$L^b$	Exp	Ly$\beta$
Date	ratioa		${10}^{45}$	Time	EWc
MJD		($\mathrm{km}{s}^{-1}$)	$\mathrm{erg}{s}^{-1}$	(s)
56954	1.74	15760	6.3	2400	-15.8
57003	1.66	16090	6.7	2400	-8.4
57062	1.61	15050	6.1	2400	-13.3
57228	0.97	-	-	2700	-
57342	1.42	-	-	2400	-
57429	1.91	-	-	2100	-
57750	2.50	-	-	2700	-26.0

Table 4. Table 4 Observational properties derived from the spectra compared with other high- z 𝑧 z BAL QSOs

Name	$z$	type	BI(CIV)	(C iv)a	(Mg ii)	(He i* )	${\alpha }_{\lambda }$b	SDSS-$i$	$M_{i^{\ast }}$d
			$\mathrm{km}{s}^{-1}$	position	REWc	REW
(From Maiolino et al. 2004a)
SDSSJ104845.05+463718.3	6.22	LoBAL	6500$\pm$1100	11000/20000			$-$2.10	22.4	–27.4
SDSSJ104433.04-012502.2	5.78	HiBAL	1950$\pm$250	5000			$-$1.55	21.6	–28.4
SDSSJ075618.14+410408.6	5.08	BAL	270$\pm$60	9000			$-$1.67	20.2	–26.6
SDSSJ160501.21$-$011220.6	4.92	LoBALe	9300$\pm$2000	11000			$-$1.35	19.8	–27.7
(Our Target)
SDSSJ012247.34+121624.0	4.82h	LoBAL	15300$\pm$2000f	12000/21000	14.64	3.54	$-$2.02	19.4	–27.9

Equations

$$M_{\mathrm{BH},\mathrm{vir}}=10^a\times (\frac{\lambda L_{\lambda }}{10^{44}\mathrm{erg}{\mathrm{s}}^{-1}}{)}^b\times (\frac{\mathrm{FWHM}}{\mathrm{km}{\mathrm{s}}^{-1}}{)}^2$$

$$\mathrm{BI}={\int }_{min}^{max}[1-\frac{f(v)}{0.9}]C\mathrm{d}v$$

$$F_{\lambda }=C{F}_{\lambda }^t{\left(\frac{\lambda }{3000\mathring{\mathrm{A}}}\right)}^{(1.62+{\alpha }_{\lambda })}10^{-\frac{A_{3000}}{2.5}\frac{A_{\lambda }}{A_{3000}}},$$

$$A_{\lambda }=k_{\lambda }E(B-V)=\frac{k_{\lambda }{A}_V}{R_V}$$

$$\dot{ϵ}(v)=\frac{1}{2}{\dot{M}}_{\mathrm{out}}^{\mathrm{ion}}{v}^2\approx 1.2\cdot 10^{44}{L}_{45}{v}_{5000}^3{r}_1^{-1}(\mathrm{erg}{\mathrm{s}}^{-1})$$

$$N=\int N_{v}dv=\sum N_v\Delta v=\frac{3.7679\times 10^{14}}{\lambda f_{jk}}\sum {\tau }_v\Delta v$$

$$\dot{ϵ}=\frac{1}{2}{\dot{M}}_{\mathrm{out}}^{\mathrm{ion}}{v}^2=2\pi \mu m_{p}R\Omega N_H{v}^3$$

Full text

††slugcomment: Draft version.11affiliationtext: Yunnan Observatories, Chinese Academy of Sciences, Kunming 650216, China22affiliationtext: Steward Observatory, University of Arizona, Tucson, AZ 85721-0065, USA33affiliationtext: CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, China44affiliationtext: Department of Astronomy & Astrophysics and Institute for Gravitation & the Cosmos, The Pennsylvania State University, 525 Davey Laboratory, University Park, PA 16802, USA55affiliationtext: Hubble Fellow66affiliationtext: Department of Physics, University of Connecticut, 2152 Hillside Rd Unit 3046, Storrs, CT 06269, USA77affiliationtext: Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China88affiliationtext: Department of Astronomy, School of Physics, Peking University, Beijing 100871, China99affiliationtext: Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China1010affiliationtext: Polar Research Institute of China, 451 Jinqiao Road, Shanghai 200136, China1111affiliationtext: Center for Astronomical Mega-science, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing, 100012, P. R. China1212affiliationtext: UKIRT Observatory, Hilo, Hawaii 96720, USA1313affiliationtext: University of the Chinese Academy of Sciences, Beijing 100049, China1414affiliationtext: Department of Physics, Pennsylvania State University, University Park, PA, 16802, USA1515affiliationtext: Key Laboratory for the Structure and Evolution of Celestial Objects, Chinese Academy of Sciences, Kunming 650216, China

1. Introduction

Extreme gaseous outflows appear most conspicuously in quasar spectra as broad absorption lines (BALs), which typically have velocity widths larger than 2000 km s*-1*, reaching velocities up to tenths of the speed of light. Previous studies based on spectral analyses indicated that broad absorption line quasars (BAL QSOs) comprise about $\sim$20% of the quasar population at low and intermediate redshifts (Hewett et al., 2003; Knigge et al., 2008; Reichard et al., 2003b; Scaringi et al., 2009), usually characterized by high-ionization broad absorption lines (HiBALs) (e.g., C iv, Si iv, N v, O vi). A small fraction ($\sim$10%) of the BAL QSOs also shows absorption troughs characterized by low-ionization species (Mg ii$\lambda$2798, Al iii$\lambda$1857) in their spectra. These are called LoBALs (in contrast to HiBALs, which only show high-ionization absorption troughs). An even rarer family of BALs are FeLoBALs, which usually show prominent Fe ii and Fe iii absorption. In general, high-ionization absorption troughs tend to be broader than the troughs associated with low-ionization species (e.g. Filiz Ak et al. 2014), suggesting quasar outflows have a wide range of ionization states. Recently, the discovery of ultra-fast outflows ($v>0.1c$, e.g., Rogerson et al., 2016) presented the phenomenon that the C iv absorption trough overlaps with that of Si iv, and Si iv with N v, thus the outflow redshift is often hard to establish without aid from wide wavelength coverage spectra.

In low-redshift cases, the reddening of normal quasars in general can be well fit by SMC-like extinction curves, while more heavily absorbed AGNs are often best characterized by a flatter extinction curve (e.g., Maiolino et al., 2001, 2004a; Gaskell et al., 2004; Gallerani et al., 2010). In contrast, reddening in BAL QSOs at high redshifts may deviate from the SMC extinction curve, perhaps due to dust absorption associated with an extreme evolutionary phase (Gallerani et al., 2010). Some studies found that HiBAL quasars are on average more reddened than non-BAL quasars. LoBAL quasars are more reddened than HiBAL quasars (e.g., Sprayberry et al., 1992), and FeLoBAL quasars are more reddened than LoBAL quasars (e.g., Reichard et al., 2003b). BAL QSOs are usually characterized by redder continua compared with non-BAL QSOs, which is commonly interpreted as reddening by dust associated with the outflowing gas (Reichard et al. 2003b). BAL QSOs are some of the most reddened objects observable at high redshift, and so they are useful for understanding internal extinction in the early Universe (e.g., Gallerani et al. 2010 )

Originally, the difference between BAL and non-BAL QSOs was mainly ascribed to orientation effects with otherwise similar intrinsic structure (Weymann et al., 1991; Goodrich et al., 1995; Hines et al., 1995; Gallagher et al., 2007). However, one dedicated study on high-redshift BAL QSOs found systematic differences between BAL and non-BAL QSOs with respect to observational properties that should be isotropic, which indicated that these two classes of quasars may be caused by intrinsically different physics to some extent (Gallerani et al., 2010). No matter how many interpretations are proposed, the basic understanding of all the BAL phenomena is that some dense material flows out from the center of the active galactic nucleus (AGN), which can be detected and characterized by observations from IR to X-ray wavelengths. For instance, the general trends that the radio morphology of BAL QSOs is more compact than those of non-BAL QSOs (Becker et al., 2000) and that BAL QSOs have UV to soft X-ray flux ratios 10 - 30 times smaller than unabsorbed quasars (Brandt et al., 2000), greatly advanced our understanding of inner physics of the two kinds of QSOs. Dust-reddened quasars, which are plausibly in an extremely active phase, appear likely to host LoBALs (Urrutia et al. 2009), indicating the LoBAL phenomenon may be linked to evolution, though the orientation effects cannot be ruled out. Some BAL quasars, especially LoBALs, may be cocooned by dust and gas rather than being quasars with a particular line of sight through their disk winds (Becker et al. 2000).

Variability studies give additional insight into the nature of BAL QSOs. For example, Filiz Ak et al. (2013) presented a comprehensive study on the variability of C iv  and Si iv  BALs over multi-year epochs in a large quasar sample, in which they suggested that global changes in ionization state are likely to be the most straightforward mechanism for explaining coordinated variability of multiple C iv  troughs at different velocities. However, some variability of LoBALs may be caused by a BAL structure moving in/out of our line of sight to the UV continuum emitting region (e.g., Hall et al., 2011; Vivek et al., 2014). Variability of BAL trough strengths is relatively common (e.g., Hall et al., 2007; Gibson et al., 2008; Filiz Ak et al., 2013), though only a few cases of dramatic changes ($>$50 %) have been reported (Crenshaw et al., 2000; Lundgren et al., 2007; Hamann et al., 2008; Leighly et al., 2009; Gibson et al., 2010; Filiz Ak et al., 2012; Grier et al., 2015; Rogerson et al., 2016; Vivek et al., 2016). Yet it is difficult to find a single trigger for BAL variability, since it may be caused by either gas motion across the line of sight or changes in the ionization state (e.g., Vivek et al., 2014; Wang et al., 2015).

In this paper, we investigate the physical properties of a newly discovered high-redshift LoBAL QSO using optical and near-IR spectroscopic observations. Such distant and powerful outflows are generally believed to exert significant feedback effects on the formation of host galaxies as well as the regulation of supermassive black holes (SMBHs) (e.g., Kauffmann et al., 2000; Richards et al., 2002; Di Matteo et al., 2005; Everett et al., 2005; Reeves et al., 2009; Rupke et al., 2011; Fabian et al., 2012). The study of this newly discovered LoBAL QSO at z $\sim$ 5 has the potential to yield information on whether or not these outflows are energetic enough to provide feedback to its host galaxy. Though LoBAL QSOs are not a large population, their extreme properties may provide more stringent and valuable tests of BAL outflow models and feedback processes than do the more numerous “normal”  BAL QSOs.

This paper is organized as follows. In Section 2, we describe the observations from the Lijiang 2.4m telescope (LJT) and UKIRT and the related data reduction we have performed. Through the analyses of all the spectra collected together, we illustrate the spectral properties of this LoBAL QSO, and estimate its black hole mass, Eddington ratio, extinction and the outflow kinetic power in Section 3. The results and implications are discussed further in Section 4, and summarized in Section 5. A flat cosmology with $H_{0}$ = 70 km s${}^{-1}\cdot$Mpc*-1*, $\Omega_{M}$ = 0.3 and $\Omega_{\Lambda}$ = 0.7 has been adopted unless stated otherwise.

2. Observations and Data Reduction

From new optical-IR selection criteria based on SDSS and WISE photometric data (Wu et al., 2012; Wang et al., 2016), which have been demonstrated to select $z\sim 5$ quasars with both high efficiency and completeness, we discovered a high-redshift BAL QSO, SDSS J012247.34+121624.0 (hereafter J0122+1216, Yi et al. (2015)), with Balnicity Index (BI) $\sim$16000 $\rm km~{}s^{-1}$(possibly a lower limit) using the LJT on 2014 October 24. As BAL QSOs with large BI and high luminosity are rare (Vestergaard et al., 2003; Sulentic et al., 2006; Trump et al., 2006; Gibson et al., 2009), we performed follow-up observations including near-IR spectroscopy with TNG/NICS and UKIRT/UIST, and optical spectroscopic monitoring observations from the Yunnan Faint Object Spectrograph and Camera (YFOSC) mounted on the LJT (see Table 1). All the spectra obtained were first corrected for Galactic extinction (Schlafly & Finkbeiner, 2011) with an average Galactic extinction law ($R_{V}$ = 3.1) and transformed into the quasar rest frame using the redshift determined by the Mg ii emission-line.

2.1. Optically spectrophotometric observations

Optical spectrophotometric monitoring observations of this BAL QSO were obtained with LJT/YFOSC. Among its various imaging and spectroscopic observing modes (Fan et al., 2015), YFOSC provides high-sensitivity, low-resolution spectroscopy, and a quick switch between imaging and spectroscopic modes when observing this faint quasar. Typical FWHM in the YFOSC data vary in the range from 0.8″to 2.0″. Each night the standard calibration data, including bias, sky-flat fields, and internal lamp flats were obtained. Exposures of the Neon and Helium lamps were used for wavelength calibration.

Photometric images of J0122+1216 were obtained before spectroscopic observations to select a nearby non-variable reference star, whose spectra could be obtained simultaneously with those for J0122+1216. We also select reference stars from the Catalina Survey (Drake et al., 2009) in the same field view of YFOSC. As a result, an F/G type star (see Table 2) was chosen as the best candidate according to its little photometric variation between 2006 Sep and 2013 July. The uncertainties of the photometric measurements include the fluctuations due to photon statistics and the scatter (of 0.03 mag or less) in the measurement of the constant-flux dwarf used. The typical total uncertainty is in the range of 0.01-0.05 mag depending on object brightness and the observing conditions. Therefore, our QSO monitoring campaign can offer an efficient and reliable way to search for possible links between emission-line and BAL-trough variability, which may provide clues into the origin of BAL outflows. The photometric magnitudes of J0122+1216 and the reference star obtained from the SDSS, UKIDSS, and WISE databases are listed in Table 2.

During the optical spectroscopic observations, we obtained more than one spectrum for each exposure. To achieve good spectrophotometry, the spectrograph’s focal plane is rotated to an appropriate position angle so that a nearby reference star could be simultaneously located in the slit. This method allows spectrophotometric calibration of the quasar even under non-photometric conditions. The simultaneous observations reduce the effects of atmospheric dispersion at non-parallactic sky angles, as well as light losses due to occasional guiding errors or poor seeing. Exposure times ranged from 2400s to 2700s for each spectroscopic run (see Fig. 1) depending on the transmission and seeing, to aim for a S/N$>$10 pixel*-1* in the continuum. The spectroscopic data were reduced using standard IRAF routines. The images were bias and flat field corrected. The extraction width is typically 10 pixels, and wavelength calibration is applied to data after the spectral extraction. This procedure results in two spectra from each image: one for J0122+1216 and the other for the comparison. The consecutive quasar/star flux ratios are then compared to correct for systematic errors in each observation and to identify and further remove cosmic-ray events. The ratio is usually reproducible to 0.5% - 5% at all wavelengths, and observations with ratios differing from the mean value by more than 5% are discarded. With the spectrophotometric runs carried out with LJT/YFOSC, we got 7 optical spectra of the target, as shown in Table 3. The spectra from LJT/YFOSC cover the wavelengths from 5500 to 9700 Å  at average resolution $R=\lambda/\Delta\lambda\approx$ 320.

To determine variability as well as maximum and minimum velocities of the C iv and Si iv absorption troughs, we smoothed the mean continuum-normalized spectrum using a boxcar smoothing algorithm over three pixels. In order to measure the intrinsic variability of the broad absorption troughs reliably rather than the overall variability of this LoBAL quasar, we made different local fits for each individual spectrum according to the neighboring line-free regions of unabsorbed continua when calculating their Balnicity Index (see more in Section 3.2). To be specific, we first take the unabsorbed parts of the Si iv and C iv emission lines as a reference to recover the continuum in the optical spectra. Moreover, we found that the broad band SED of J0122+1216 is similar to that of a quasar (SDSS J100713.68+285348.4) from Zhou et al. (2010); thus the continuum level could be estimated using the latter as a comparison template to recover the absorption-free spectrum. In addition, a mean spectrum of J0122+1216 and ratio spectra compared to the mean composite are obtained as shown in Fig. 2. At a first glance, the emission/absorption lines appear to vary to some extent during these times. This way of determining variability of emission/absorption lines, however, yields large uncertainties for our lower-S/N spectra and is affected by sky background contamination.

Fortunately, one feasible way to alleviate uncertainties from continuum fitting is to study the flux ratios between our target and the comparison star from simultaneous spectroscopic observations during different epochs. When using the spectra of reference star to recover the continuum level for all the spectra of J0122+1216 obtained from the monitoring campaign, the variability of the C iv absorption troughs is not detected (less than 10%) , which is consistent with most low-redshift BAL QSOs previously studied. Our flux calibration is robust due to simultaneous observations for the non-variable star within the same slit. This also minimizes any possible spectral variability that may be caused either by system fluctuations or weather changes.

We then determined trough velocities by including only regions with flux below 90% of the estimated continuum level (i.e., where the normalized flux density is $<0.9$). A maximum velocity of $\sim 29000$ $\rm km~{}s^{-1}$  was set to avoid interference by broad lines since outflows at such high velocities are very rare in the UV band. It is clear that both C iv  and Si iv  have similar velocity structures, which is particularly characterized by their two deep absorption troughs. For N v , a shallow absorption trough seems to be imprinted at $\sim 2000$ $\rm km~{}s^{-1}$, but it is impossible to obtain other absorption features at higher velocities due to the heavily absorbed Ly $\alpha$  forest. In addition, the Si iv and C iv  BALs have a similar trough velocity width of 23000 $\rm km~{}s^{-1}$. The maximum and minimum velocities of the Si iv and C iv troughs in each epoch were measured separately, and were consistent over all epochs within the wavelength calibration uncertainty.

2.2. Near-IR spectroscopic observations

With director’s discretionary time from the TNG observatory, the first near-IR spectrum was obtained by TNG/NICS on 2015 Feb 13 through service mode. Unfortunately, the S/N of the raw spectrum was too low to be analyzed due to poor seeing and transparency plus large airmass effects. Successful near-IR spectroscopic observations were carried out on 2015 Oct 30 in Hawaii using UKIRT/UIST with an ABBA dither pattern, and each exposure was set to 240 seconds to subtract optimally the sky emission lines. Based on the analysis of the performance of UIST and the fact that this target is faint in the JHK bands, we adopted a 0.8″slit and HK grism, which covers 1.395$\mu$m to $2.506\mu$m with average spectral resolution of $R\sim$ 350. Bias, flat-field and standard-star observations were also taken before and after exposures of the target. The IR observations had excellent seeing ($\sim 0.4$ ″  in H-band) and low humidity ($\tau\sim 0.05$). The data-reduction was automatically performed using the UIST pipeline based on Starlink packages.

Because the spectral resolution from UKIRT/UIST is low, another high-resolution near-IR spectrum covering the JHK bands was obtained from P200/TripleSpec on Feb 9, 2017 for the purpose of investigating further details. In order to obtain a higher S/N spectrum of J0122+1216, we took 10 Fowler depths to optimize the readout noise after each single exposure using the ABBA mode. The total exposure time was initially assigned to be 5400s, but we then reduced the time to $300\times 8$s when noticing the S/N of each single spectrum started to weaken. This is probably due to the average seeing becoming larger than 1.5″  and clouds starting to increase. The data reduction of all spectra was performed using Spextool packages provided by the Palomar Observatory website. The S/N of the combined spectra from P200/TripleSpec was low; however, this spectrum can be used for determining the redshift associated with that from UKIRT/UIST when noticing they have similar Mg ii  emission/absorption profiles (see Fig. 4).

With the aid of the near-IR spectrum obtained by UKIRT/UIST and P200/TripleSpec, this quasar is confirmed as a LoBAL QSO with the largest BI(C iv) $\sim 16000$ $\rm km~{}s^{-1}$  at $z>3$, matching more tentative observations in the discovery paper. Some low-ionization absorption lines (red dashed lines in Fig. 5) such as Mg ii  and He i* were also found with greater than 2$\sigma$ significance, while only tentative evidence for other lines like Fe II is found (i.e., they are found at less than 2$\sigma$ significance). These low-ionization absorption lines, in turn, provide additional diagnostics to study outflow properties.

3. Spectral properties and analyses

Through a comprehensive analysis of C iv $\lambda$1549, previous studies revealed that blueshifted emission lines associated with outflows in radio quiet quasars are more common and larger amplitude at high luminosity (e.g., Wang et al., 2011; Marziani et al., 2016). Due to the blueshifts of C iv  and other high-ionization emission lines, the low-ionization Mg ii  emission-line is more useful for estimating the systemic redshift (Richards et al., 2002; Hewett et al., 2010; Shen et al., 2016). This argument has been further supported by an observed symmetry of the Mg ii emission-line. The Mg ii line in our near-IR spectrum is consistent with being symmetric, leading us to conclude that it is a reliable indicator of the systemic redshift of J0122+1216.

The initial redshift of $z=4.76$ was determined by Yi et al. 2015 based on optical spectra using the Ly$\alpha$, Ly$\beta$, C iv and Si iv emission lines matched to the average composite template, while our new redshift determined from the Mg ii emission-line in the near-IR spectroscopy is 4.82$\pm 0.01$. With the latter redshift, the unabsorbed emission line peaks of Ly$\beta$, C iv and Si iv exhibit blueward shifts, which are commonly seen among low-redshift BAL QSOs (e.g., Richards et al., 2002; Wang et al., 2011), but most of them are rarely with such a dramatic blueshift. This scenario, in fact, can be further corroborated with the evidence of significant variability among these emission lines (see Table 3), which is discussed more in the following sections.

Up to now, there is no widely accepted definition to classify LoBAL or FeLoBAL QSOs with respect to the width and depth of the Mg ii or Fe ii absorption lines (e.g., Trump et al., 2006; Gibson et al., 2009; Zhang et al., 2010). At wavelengths shortward of 14500 Å, the increasing atmospheric opacity makes detailed line identifications difficult. With this caveat in mind, the near-IR spectrum has potential Fe ii absorption features (e.g., the UV62 and UV63 components in the Fe absorption atlas and shallow Fe ii UV1 absorption at $\sim$ 2600 Å ; see Hall et al. 2002). Alternatively, it may be experiencing a transition phase between FeLoBAL and LoBAL QSOs characterized by the emergence or disappearance of Fe ii/ Fe iii absorption lines. The possibility of blended Fe emission/absorption lines makes the measurement of its continuum more uncertain in this spectral region.

Adopting the method from Liu et al. (2015), we quantitatively analyzed the reliability of different absorption features in the near-IR spectrum. All three of the low-ionization absorption lines have velocity widths larger than 1600 $\rm km~{}s^{-1}$. Among them, the Mg ii and He i* absorption lines are detected at 3.2$\sigma$ and 2.5$\sigma$ significance, respectively, while Fe ii was only detected at 0.9$\sigma$ significance and should thus be treated with caution. However, we argue that the Fe ii  absorption feature is likely to be real, given that the Fe ii absorption trough peaks lie at the same position in the common velocity absorption system with respect to the Mg ii , He i* and Fe ii outflowing ions (Fig. 6). These absorption lines have similar detached velocities ($\sim 8000$ $\rm km~{}s^{-1}$) and velocity widths ( $\sim 2000$ $\rm km~{}s^{-1}$). Such observational properties would be useful for setting physical constraints with respect to the origin of outflows. However, a more detailed study of these ions would require spectra of higher S/N and resolution.

To determine whether this newly discovered LoBAL QSO is typical when compared to the general BAL QSO population, we compare its properties with those of the BAL QSO sample of Gibson et al. (2009). J0122+1216 has obvious emission/absorption blueshifts that are commonly seen among the BAL population, but J0122+1216 itself is less reddened than typical BAL QSOs. Moreover, it is clear to see the steep onset of the C iv and Ly$\alpha$ absorptions, which may be caused by inherently similar geormetry and physics of absorption. The velocity width of the C iv trough is not unusual, but is significantly larger than 90% of C iv BAL troughs (Gibson et al., 2009). Though the velocity of the C iv trough falls within the normal range of the BAL QSO population, the maximum velocity of the trough falls toward the high end of the distribution. In general, neither the emission nor its absorption features are strong outliers in the observational properties examined, though they do fall in less common regions with respect to the BI index, the velocity width and range.

3.1. Estimated BH mass based on the Mg ii line

The single-epoch Mg ii  estimator derived from reverberation mapping has been applied to the study of large samples of AGN at $1.5\leq z\leq 5$, and it is also preferentially adopted for high-luminosity cases at high redshifts in recent near-IR spectroscopic studies (e.g., Jiang et al., 2007; Yi et al., 2014; Wang et al., 2015; Zuo et al., 2015; Wu et al., 2015). Furthermore, the emission-line width of Mg ii  has been demonstrated to correlate with that of H$\beta$ in single-epoch spectra (e.g., Salviander et al., 2007; Shen et al., 2008; McGill et al., 2008; Wang et al., 2009; Shen et al., 2016), so it has become a preferred estimator of BH masses from intermediate redshifts to the most-distant quasars (e.g., Shen et al., 2011; Willott et al., 2010; Trakhtenbrot et al., 2012; Zuo et al., 2015).

The empirical single-epoch relation usually adopted to estimate SMBH masses from the Mg II line and its associated continuum at 3000 Å is as follows:

[TABLE]

This relation, in fact, is derived from a virial assumption of gaseous in the broad emission line region. Therefore, both $a$ and $b$ are determined by the virial coefficient $f$. This factor $f$ could be calibrated by an independent way (such as the empirical $M-\sigma$ relation) to a higher confidence level. Therefore, finding reliable measurements for the Mg ii  line and the $f$ factor would be critically important to estimate the BH mass. However, different fits among these broad emission lines (Wang et al., 2009), different $f$ factors between classical bulges and pseudobulges (Ho & Kim, 2014) , and a potentially new $R_{BLR}-L$ relation depending upon accretion rates (Du et al., 2016), should be taken into account when estimating the BH mass, which in turn, make it hard to mitigate intrinsic uncertainties of the single-epoch spectral relation. For simplicity, $a=0.86$ and $b=0.5$ are adopted in this paper respectively, according to the study from Vestergaard et al. (2009).

Though the Mg ii absorption line appears on the blueward part of the emission line, the measurement of FWHM should be reliable from the symmetric profile fit of the emission line. Here, we calculated the 3000Å  monochromatic luminosity ($L_{\rm 3000\AA~{}}\sim 5.81\times 10^{46}$ $\rm erg~{}s^{-1}$) based on the analysis of the near-IR spectrum after correction for Galactic extinction (Schlafly & Finkbeiner, 2011). We then assume a bolometric correction factor of $\sim 5.15\times L_{\rm 3000\AA~{}}$ (Richards et al., 2006) to obtain the bolometric luminosity after subtracting the UV Fe II contribution. Noticing random atmospheric changes may lead to a large uncertainty for near-IR spectral flux calibration, an alternative method based on the luminosity relation ($L_{\rm 3000\AA~{}}\sim 1.62L_{\rm 5100\AA~{}}$, see Trakhtenbrot et al., 2012) is also taken into account. Since Mg ii BALs are more frequently detected in quasars with narrower H$\beta$ and weaker [O iii] emission lines, stronger UV and optical Fe ii multiplets, and higher luminosity (e.g., Zhang et al., 2010), we neglect the [O iii] emission contribution on the marginal blueward end of the W1 bandpass. This is also consistent with the anti-correlation between Fe ii and [O iii] strength (Boroson & Green, 1992). Then $L_{\rm 3000\AA~{}}\sim 6.19\times 10^{46}$ $\rm erg~{}s^{-1}$  was derived. Using the empirical single-epoch relation with FWHM of the Mg ii line $\sim$3600 km s*-1* and the 3000 Å  monochromatic luminosity, the BH mass is estimated to be $\sim 2.3\times 10^{9}M_{\odot}$. The Eddington ratio, a non-dimensional physical parameter defined by the bolometric luminosity over its Eddington luminosity, is $\sim 1.0$ based on the flux calibration and spectral fit (see Fig. 7).

A BAL QSO such as this with extreme outflows might have additional velocity structure even in Mg ii  in addition to the gas motions dominated by gravitation that are assumed to provide the measure of the black hole mass (Plotkin et al., 2015). We also note the possibility of systematic error in the Mg ii estimator for very high luminosities that is assumed to be the same as that for the low-luminosity QSOs that have been reverberation mapped. In such a case, the inherent BH mass uncertainty of J0122+1216 would be higher than that of the single-epoch spectral relation ($\sim 0.3$ dex) according to previous studies at low and intermediate redshifts (e.g., Vestergaard et al., 2009; Wang et al., 2009; Shen et al., 2011).

3.2. Determining the BI(Mg ii)

To calculate the BI, we first normalize the spectrum by dividing by the continuum fit to the observed spectrum in less-absorbed spectral regions around 1700 Å  and to the red of 3500 Å  in the rest frame. Combining the continuum derived from optical and near-IR spectra provides the framework to analyze absorption troughs quantitatively for both the Mg ii and C iv lines.

For the LoBAL QSO J0122+1216, with the detached Mg ii  absorption trough far away from its emission peak, we adopt the traditional absorption BI to be consistent with the discovery paper, which is defined as follows (Weymann et al., 1991):

[TABLE]

$f(v)$ is the continuum-normalized spectral flux at a velocity $v$ (in $\rm km~{}s^{-1}$) in the rest frame. The dimensionless value $C$, is set to 1 at velocity width more than 2000 km $\rm s^{-1}$ from the start of a continuous trough with flux density less than 90% of the continuum, and $C=0$ elsewhere. In addition, we define velocities flowing outward along our line of sight to be negative, which are distinguished from the QSO emission rest frame (see Fig 8). Here, the 3000 km s*-1* red limit and 10000 km s*-1* blue limit are set to avoid absorption lines unrelated to the Mg ii  ion. Further, the Mg ii  velocity components are not seen at higher velocities. Then, the BI(Mg ii) is calculated to be $\sim$ 1350 $\rm km~{}s^{-1}$, with the main uncertainty being the determination of the intrinsic continuum.

3.3. Broad-Band Spectral Energy Distribution

Extensive studies of interstellar attenuation from reddening and scattering have been done, but most of them are focused on the local universe. Generally, the SMC extinction curve is steeper than those of the LMC, the MW, and starburst galaxies (Calzetti et al., 1994). However, some individual studies showed that extinction curves for LoBAL QSOs may be steeper than the SMC extinction curve shortward of 2000 Å  (e.g., Hall et al., 2002). In this scenario, a smaller average grain size is needed to produce steeper extinction curves, since small particles cannot scatter efficiently at long wavelengths. In contrast, the flat extinction curves in the circumnuclear regions of Seyfert galaxies could be interpreted as being produced from dust dominated by very large grains.

To match the observed spectra, different quasar templates were adopted mainly by adjusting the continuum slope ($\alpha_{\lambda}$), the absolute extinction ($A_{3000}$), and by applying some empirical and theoretical extinction curves. The average intrinsic slope ($\alpha_{\lambda}$) of BAL QSOs is about $-2.02$ at low redshifts (Reichard et al., 2003a), so we took an unreddened QSO template with this slope to fit the continuum sections at $\lambda>3000$ Å, which would be less affected by intrinsic extinction. In order to match the composite quasar spectrum to the input spectra, we allow for changes in the slope and shape of the continuum between the input quasar spectrum and the template quasar spectrum, and the spectrum of the LoBAL is fitted with the following expression:

[TABLE]

where $C$ is a normalization constant, $F_{\lambda}^{t}$ is a quasar template spectrum, $\alpha_{\lambda}$ is the intrinsic slope of the unreddened spectrum, and $A_{\lambda}/A_{3000}$ is the extinction curve normalized at 3000 Å, which is derived from empirical and theoretical extinction curves (e.g., Maiolino et al., 2004a; Gallerani et al., 2010). The extinction at a wavelength $\lambda$ is related to the color excess $E(B-V)$ and to the reddening curve by:

[TABLE]

Here, we adopt the template compiled by Reichard et al. (2003a), obtained using 892 quasars classified as non-BAL, whose continuum is described by a power-law with average index $\alpha_{\lambda,t}\sim-1.6$ ($F_{\lambda}^{t}\propto\lambda^{\alpha_{\lambda,t}}$). We change the slope of the template with the term $(\lambda/3000)^{1.62+\alpha_{\lambda}}$. The redward part of the spectra at $\lambda_{rest}>3000$ Å  is the least-extinguished region and so is used to anchor the reddening fit. However, we did not include $2500<\lambda<3000$ Å  in the fit since it has the prominent Fe ii bump and absorptions as shown in Figure 7. Meanwhile, we exclude the spectral format edges ($<2500$ Å  and $>4300$ Å ) from the reddening fit. The continuum emission peak between Si iv and Ly $\alpha$ was also excluded because it falls into the strong atmospheric absorption region around $\sim$7600 Å . Adopting an intrinsic spectrum with slope of $-2.02$ and an absolute extinction at 3000 Å  of 0.35 mag, the blue-part continuum slope of J0122+1216 is estimated to be $\sim-1.6$, which is very close to that of the average composite spectral slope. In principle, the slope of the optical spectrum, from the Mg ii emission-line profile redward in the observed-frame near-IR spectrum to rest wavelengths larger than 3500 Å , could be used to visually check the consistency of our fit to the whole observed spectra (see the cyan dashed line in Fig. 9). However, for a given extinction curve $f_{\lambda}=A_{\lambda}/A_{3000}$, any other extinction curve of the form $f_{\lambda}^{\prime}=k(f_{\lambda}-1)+1$, can fit the observed spectrum equally well by simply changing the absolute extinction $A_{3000}$. Thus, we cannot eliminate a degeneracy in the slope of the extinction curve without knowledge of the intrinsic luminosity.

On the other hand, adopting a lower uncertainty factor of $\sim$2 from the single-epoch spectral relation, the maximum BH mass and intrinsic bolometric luminosity are $\sim 4.6\times 10^{9}M_{\odot}$ and $\sim 6.3\times 10^{47}$ $\rm erg~{}s^{-1}$  constrained by the conventional Eddington limit, which corresponds to a monochromatic luminosity of $1.2\times 10^{47}$ $\rm erg~{}s^{-1}$  at 3000 Å . According to these values, the maximum absolute extinction $A_{3000}$ should not be larger than $\sim 0.72$ mag giving an intrinsic continuum slope not steeper than $-2.8$. Very few QSOs have been found with such a blue UV slope, and a new Eddington limit ($4.4\times 10^{47}$ $\rm erg~{}s^{-1}$) could be derived after one iteration using the monochromatic 3000 Å  luminosity and relation (1). Thus we adopt the values of intrinsic continuum slope of $-2.02$ and 0.35 mag extinction at 3000 Å  as reasonable and matching the observed spectra, application of which yields the final spectral fit shown in Figure 9.

The assumptions discussed above were made for a reasonable case since it is hard to determine the exact reddening when considering the degeneracy between $\alpha_{\lambda}$ and $A_{V}$. Among the repeated fits of the intrinsic and observed spectra, the SMC extinction curve seems to deviate substantially from that of the LoBAL QSO J0122+1216 when considering the fact that it will make the whole spectrum redder. However, a detailed investigation of the extinction properties of J0122+1216 would require higher S/N spectral data obtained in simultaneous observations and wider wavelength coverages, which is beyond the scope of this paper.

3.4. Physical constraints on the outflows

AGN feedback seems to play an important role in controlling the co-evolution between central SMBHs and host galaxies over a wide range of cosmological scales, and thus it is essential to understand the detailed structure and physics of outflows (e.g., Di Matteo et al., 2005). Generally, both emission and absorption lines can be efficient diagnostics of the geometry and dynamics of quasars. To be specific, the emission lines are usually used to determine the black-hole mass or even to derive a global covering factor of outflows; the absorption lines, which are associated with outflows in the line of sight, are good for determining the total hydrogen column density ($N_{H}$) and the velocity distribution of the outflows (Liu et al., 2015).

3.4.1 The high-ionization C iv emission-line

Since the velocity field and gas density distribution cannot be determined according to the analyses of the optical and near-IR spectra available, an alternative way to estimate the kinetic power of J0122+1216 is necessary. The dramatic blueshifts of high-ionization lines as shown in Figure 5 provide evidence that the C iv line region may be dominated by outflows rather than the gravitational virial motions prevalent in most normal quasars at low redshifts (Wang et al., 2011), or at least, these blueshifted lines are partially caused by outflows. Therefore, a simplified method based on an outflow model was adopted to obtain an order of magnitude estimate of the kinetic power associated with the C iv emitting outflows under several assumptions. We start by considering the fraction of the flux in the C iv blue component that is above the expected projected escape velocity at $r\sim 1500$ $r_{g}$. This radius ($\sim 10$ times larger than that of the C iv broad emission line region) is set to make sure it will cover the relevant broad emission line regions.

The outflowing kinetic power, at the velocity of $v$ in units of 5000 km s*-1* and the integrated C iv line luminosity in units of $\rm erg~{}s^{-1}$  (Marziani et al., 2016), is:

[TABLE]

where $r^{-1}_{\rm 1}$ is the inner radius of the C iv emission-line outflowing region (here $r_{1}=1500r_{g}\sim 2.8$ pc is assumed). First, we measure the integrated C iv  line luminosity ($\sim 4.3\times 10^{45}$ $\rm erg~{}s^{-1}$) after subtraction of the continuum. From the blueshifted side of the C iv  line ($\sim$5000 $\rm km~{}s^{-1}$ ) which is likely driven by outflows, the kinetic power associated with the C iv  line can be roughly estimated to be $\sim$4.0 $\times 10^{45}$ $\rm erg~{}s^{-1}$  ($\sim$1.2% LEdd). We note that inner radius of $1500r_{g}$ and an electron density of $10^{9}$ cm*-3* are assumed to derive equation (5), but both of them are likely linked to lower limits when compared with physical properties of the C iv  emission-line region. Thus, adopting this simplified approach, the outflow kinetic power may be underestimated to some extent.

3.4.2 Analysis based on absorption lines

As a comparison, we also follow a conventional procedure to analyze the C iv , He i3889 and Mg ii absorption troughs, which is helpful to study outflowing material along our line of sight. For an ion giving rise to an unsaturated absorption line, the optical depth at velocity $v$ is directly proportional to the column density corresponding to $v$ (here we use $N_{v}$ to represent this physical quantity). Then, the total column density $N$, in units of cm-2*, can be approximated by

[TABLE]

where $f_{jk}$ is the oscillator strength corresponding to a specific ion, $\tau_{v}$ is the optical depth at velocity $v$, and $\Delta v$ is the velocity width element in the center of $v$. Using Equation 6, we derived the He i*  column density to be $4.5\times 10^{14}$ cm*-2*, adopting the oscillator strength of 0.064, a velocity width of 2000 $\rm km~{}s^{-1}$  and an apparent optical depth of 0.15.

The lower limit on the total hydrogen column density ($N_{H}$), which tightly depends on the ionization states and abundances, can also be estimated to some extent. Because there are still substantial uncertainties that depend on results from photoionization simulations, we perform only a qualitative analysis here. First, the He i*  absorption line multiplets can be used as probes to obtain the real ion column density as opposed to a lower limit available from the apparent column density (Arav et al., 2001). Then, the column density of helium in the metastable state ($4.5\times 10^{14}$ cm*-2*) derived from the He i*($\lambda 3889$) absorption line can yield a conservative lower limit on the ionized helium column density and hence the ionized hydrogen column density ($N_{HII}\sim 10^{21}$ cm*-2* ). Furthermore, a lower limit on the hydrogen column density ($N_{H}=N_{HI}+N_{HII}>10^{22.5}$ cm*-2*) can be obtained according to the detailed photoionization simulations for reasonable parameter ranges (e.g. Leighly et al., 2011). The column density estimated above, in fact, is also corroborated by other observational features such as a LoBAL with Fe ii  absorption that would require high column densities (e.g., Mudd et al., 2016).

With a lower limit on the total hydrogen column density ($3.2\times 10^{22}$ cm*-2*), the fraction of outflowing kinetic power can be somewhat constrained. The relation of outflowing kinetic power is expressed as follows:

[TABLE]

where $\mu m_{p}$=1.4$m_{p}$ is the mean mass per proton in the absorber, and $\Omega$ is the global covering factor.

For the reliably detected Mg ii BALs of J0122+1216 at $z=4.82$, a mean absorbed flux-weighted centroid velocity of 8000 $\rm km~{}s^{-1}$, an assumed value of $R=60$ pc (e.g., Ji et al., 2015; Liu et al., 2015), and a global covering factor of $\Omega=0.1$ were adopted according to investigations of LoBAL QSOs at $z<1$ (e.g., Trump et al., 2006; Zhang et al., 2010; Liu et al., 2015). The lower limit on kinetic power is estimated to be 2.8 $\times 10^{45}$ $\rm erg~{}s^{-1}$  ($\sim$0.9% of LEdd) using Equation (7). Here, we adopted a lower value of $\Omega$, which is mainly based on optical spectroscopic statistics, and thus could be underestimated due to obscuration effects.

4. Discussion

Understanding the origin of outflows and feedbacks are fundamental to our understanding of various observational phenomena of BALs. Some authors have proposed that the different line types represent different manifestations of a single outflow phenomenon viewed at different angles (e.g, Elvis, 2000). However, BAL QSOs only appear in a small fraction of the observed spectra, indicating that either they are only seen in particular directions with respect to the axis of the accretion disk or they are observed during particular phases of the BAL QSO life (Farrah et al. 2007, and references therein). Thus, the possibility still remains that a combination of evolutionary and orientation effects can explain the separation of BAL QSOs and non-BAL QSOs. For example, the structure of the outflows could change either with cosmic time or an orientation dependence, but it is difficult or impossible to favor one of them based on observations of individual objects.

A redshift dependence of the broad absorption line quasar fraction has been found by Allen et al.(2011), implying that an orientation effect alone is not sufficient to explain the presence of BAL troughs in some but not all quasar spectra. The outflows giving rise to these broad absorption lines once had been thought to originate in accretion disk instabilities and are accelerated to velocities of up to $0.2c$ in the quasar rest-frame (Hamann et al., 2013) by magnetic effects (e.g., Everett et al., 2005) or thermal winds (e.g., Giustini et al., 2012) from the accretion disk. For the LoBAL J0122+1216, the idea of powerful radiation pressure driven by the central SMBH seems to be plausible in terms of the outflows appearing among a fairly large fraction of luminous QSOs (e.g., Arav et al., 1995; Proga et al., 2000; Chelouche & Netzer, 2001), especially noticing its high luminosity and Eddington ratio.

Boroson (2002) found that BAL QSOs on average have higher Eddington ratios and accretion rates than those of non-BAL QSOs in a small sample of BAL QSOs, and a similar conclusion has been corroborated by recent studies (e.g., Ganguly et al., 2007; Zubovas et al., 2013). Moreover, some studies suggest that BAL QSOs are redder and/or more luminous than other quasars (Reichard et al., 2003a; Trump et al., 2006; Gibson et al., 2009). The high Eddington ratio of the LoBAL QSO J0122+1216, in fact, fits into this scenario of violently accreting systems having stronger outflows. Further, observational evidence to support the luminosity - outflow connection has been found between the blueshift and asymmetry of the C iv profile and the Eddington ratio (e.g.,Laor & Brandt, 2002; Wang et al., 2011). It is possible that high Eddington ratios, even though not the main driver in the feedback processes of BAL QSOs, may control the formation of extreme outflows with sub-relativistic components under extreme circumnuclear conditions. Though the variability of Ly$\alpha$ and Ly$\beta$ emission lines has been detected, no corresponding pattern with the variability of C iv or Si iv absorption troughs has been found. Thus, due to the effects of low S/N spectra and a short time span of monitoring, we cannot establish a relation between the high-ionization broad emission/absorption lines to pin down the origin of circumnuclear outflows.

According to theoretical models, at least 0.5% of $L_{\rm Edd}$ is required to provide substantial feedback to the host galaxy in the form of outflowing material (e.g., Di Matteo et al., 2005; Hopkins et al., 2005, 2010). The lower limit on kinetic power estimated from both emission and absorption lines of J0122+1216 are already larger than that required for significant AGN feedback. Due to the uncertainties of estimating the kinetic power of outflows, which are typically one order of magnitude or even more, we caution that the method adopted in this paper may be more uncertain given all the assumptions adopted. However, a lower limit on the outflowing kinetic power of J0122+1216 can be safely set through the analyses mentioned above. Thus, we find supporting observational evidence that the circumnuclear outflows of J0122+1216 could contribute significantly to the origin of galactic-scale feedback effects, based on the measurements of the BH mass and bolometric luminosity.

5. Conclusion

The quasi-ubiquitous outflows in quasars may play an important role on sub-parsec as well as kilo-parsec scales in controlling the growth of the central black hole, the evolution of the host galaxy, and the chemical enrichment of the intergalactic medium according to observational and theoretical predictions (e.g., Moe et al., 2009; Ostriker et al., 2010), yet more observational evidence is needed, especially for high-redshift cases.

There are various methods being adopted to learn about the evolution of quasars, however, studying absorption lines from the outflows of BAL QSOs is the most direct and efficient way to learn their geometrical structures and chemical abundances, kinematics and extinction, or even more fundamental physics in terms of the accretion disk, the powering mechanism, the feedback and growth of SMBH processes etc. in the early Universe. In this paper, the follow-up investigations of this newly discovered LoBAL QSO have been made from an observational view, and we briefly conclude our findings and suggested implications as follows:

1)Based on near-IR spectra with Mg ii absorption/emission lines, J0122+1216 was identified to be a LoBAL QSO, and its redshift was determined to be 4.82 instead of 4.76. With this systemic redshift, the high-ionization UV emission lines present large blueshifts, which could be directly associated with extreme activity of circumnuclear outflows.

2)During 16 months of monitoring, the variability of C iv absorption troughs is small after spectrophotometric calibration, while it is conspicuous for the Ly$\alpha$ and Ly$\beta$ emission lines, which again, could be an observational signature characterizing extreme outflows.

3)The Mg ii-based BH mass of this LoBAL QSO is about 2.3 $\times 10^{9}$ M⊙ according to the empirical single-epoch spectral relation. However, the intrinsic deviation from the relation could be larger than the typical 0.3 dex, in terms of these dramatic outflows launching from the high-luminosity region at such a redshift.

4)We find observational evidence that the Eddington ratio of this LoBAL QSO is very close to or even beyond the conventional limit, which may become the main driver of its sub-relativistic outflows under extreme circumnuclear environments.

Outflows are generally believed to be fundamental to understand the process of AGN feedback, so a better understanding of their geometries and kinematics could shed light on the physical basis of BAL QSOs. Considering the fact that large variabilities of BALs usually happen during a relatively long time and the typical life time of outflows is longer than several tens of years in the rest-frame for most BAL QSOs, it is meaningful to put this newly discovered LoBAL QSO in a long-term monitoring campaign for further study. In addition, some observational findings, (e.g., X-ray emission from BAL AGNs is, on average, much weaker than non-BAL AGNs, and their circumnuclear regions tend to be more compact and polarized), are of great interest to advance our understanding of BAL QSOs. Therefore, more extensive follow-up observations would be helpful to confirm these findings and further reveal intrinsic physics of the LoBAL QSO J0122+1216.

We acknowledge stimulating discussions with Shaohua Zhang, Wenjuan Liu and Xiaobo Dong. We acknowledge the support of the staff of the Lijiang 2.4m telescope (LJT). Funding for the telescope has been provided by CAS and the People’s Government of Yunnan Province. W. Yi thanks the support from the West Light Foundation of The Chinese Academy of Sciences (Y6XB016001). W. Yi also thanks the financial support from the program of China Scholarships Council No. 201604910001 for his postdoctoral study at the Pennsylvania State University. CJG and WNB acknowledge support from NSF grants AST-1516784 and AST-1517113. The United Kingdom Infrared Telescope (UKIRT) is supported by NASA and operated under an agreement among the University of Hawaii, the University of Arizona, and Lockheed Martin Advanced Technology Center; operations are enabled through the cooperation of the East Asian Observatory. This research uses data obtained through the Telescope Access Program (TAP), which has been funded by the National Astronomical Observatories of China, the Chinese Academy of Sciences (the Strategic Priority Research Program ”The Emergence of Cosmological Structures” Grant No. XDB09000000), and the Special Fund for Astronomy from the Ministry of Finance. Observations obtained with the Hale Telescope at Palomar Observatory were obtained as part of an agreement between the National Astronomical Observatories, Chinese Academy of Sciences, and the California Institute of Technology. We used the quasar spectral fitting code kindly provided by Yue Shen. We also acknowledge the assistance from Antonio Magazzu, Emilio Molinari and Thomas Augusteijn at the TNG and NOT observatories. Facilities: The Lijiang 2.4m Telescope (LJT), YFOSC, TNG, NICS, UKIRT, UIST, P200, TripleSpec.