HE0359-3959: an extremely radiating quasar

We present a multiwavelength spectral study of the quasar HE0359-3959, which has been identified as an extreme radiating source at intermediate redshift (z=1.5209). Along the spectral range, the different ionic species give information about the substructures in the broad line region. The presence of a powerful outflow with an extreme blueshifted velocity of $\sim$-6000$\pm$500 km s$^{-1}$ is shown in the CIV{\lambda}1549 emission line. A prominent blueshifted component is also associated with the 1900{\AA} blend, resembling the one observed in CIV{\lambda}1549. We detect a strong contribution of very the low-ionization lines, FeII and Near-Infrared CaII triplet. We find that the physical conditions for the low, intermediate and high-ionization emission lines are different, which indicate that the emission lines are emitted in different zones of the broad line region. The asymmetries shown by the profiles reveal different forces over emitter zones. The high-ionization region is strongly dominated by radiation forces, which also affect the low and intermediate-ionization emitter region, commonly governed by virial motions. These results support the idea that highly radiating sources host a slim disk.


EXTREME POPULATION A SOURCES ALONG THE 4DE1 MAIN SEQUENCE
The 4D Eigenvector 1 (4DE1) parameter space offers a formalism to distinguish and classify type 1 Active Galactic Nuclei (AGN) considering their spectral properties (Sulentic et al., 2000a,b). The Full Width at Half Maximum (FWHM) of Hβ broad component (Hβ BC ), the strength of optical FeII blend at 4,570Å described by the ratio R FeII = I(FeII)/I(Hβ BC ), the velocity shift of the CIVλ1549 profile, and soft X-ray photon index (Ŵ soft ), provide four observationally independent dimensions of the Eigenvector 1. In the 4DE1 optical plane, the type 1 AGN occupy a welldefined sequence, driven mainly by the Eddington ratio, L/L Edd . Along this sequence we observe a variation of the physical parameters and orientation. Then, 4DE1 could be revealing an evolution sequence for type 1 AGN (Sulentic et al., 2000a;Marziani et al., 2010;Zamfir et al., 2010). For more information about the 4DE1 and update of results, see Marziani et al. in this volume. Using the 4DE1 we identify two populations with different spectral features: A and B. Population A has a FWHM (Hβ BC ) ≤4,000 km s −1 . It shows large blue asymmetries in the high-ionization lines like CIVλ1549, and it is majority populated by radio quiet sources. In contrast, population B shows a FWHM (Hβ BC )>4,000 km s −1 and it is mostly composed of radio-loud sources (Sulentic et al., 2002;Zamfir et al., 2010). Each population can be divided into small bins with FWHM (Hβ BC ) = 4,000 km s −1 and R FeII = 0.5, defining subpopulations shown in the Figure 1. In this paper we focus in the subpopulation A3 and A4 (R FeII > 1), which have been identified as highly radiating sources (xA, Marziani and Sulentic, 2014). These kind of sources show high Eddington ratios (L/L Edd > 0.2) probably produced by a slim disk, which is geometrically and optically thick and it could be formed in an advection-dominated accretion flow (Abramowicz et al., 1988;Abramowicz and Straub, 2014) .
We have found selection criteria to identify the xA sources based on the 4DE1 formalism. In the optical region they show a R FeII > 1 (high intensity of FeII) and in the UV range AlIIIλ1860/SiIII]λ1892 ≥ 0.5 and CIII]λ1909/SiIII]λ1892 ≤ 1.0 . Also, they show strong blueshifted components associated with the high ionization lines, for example in CIVλ1549 emission line, indicating the presence of outflows. More details about the xA sources behavior can be found in Martínez-Aldama et al. of this volume.

HE0359-3959: An Extreme xA Source
In our extreme luminosity Hamburg-ESO sample (Marziani et al., 2009;Sulentic et al., 2017), we have identified four cases of highly radiating quasars that show an extreme behavior, i.e., a high Eddington ratio and a strong blue asymmetry [c(1/2)<-4,000 km s −1 ; centroid at half intensity] in the CIVλ1549 profile (Sulentic et al., 2017). The most extreme case corresponds to the quasars HE0359-3959, with z = 1.5209, FIGURE 1 | 4DE1 Optical Plane reproduced from Marziani and Sulentic (2014). Gray points correspond to the sample of 470 bright low-z QSOs from Zamfir et al. (2010). The plane is divided in bins according to Sulentic et al. (2002). Extreme accretor population A sources (xA) are located in A3 and A4 bins. The black dot indicates the position of 1 Zw 1, the prototype of low-z xA sources. And, the red dot marks the location of HE0359-3959, an extreme xA source with high-z. log(L bol ) = 47.6 erg s −1 and a R FeII = 1.12. It is cataloged as an A3 source (see Figure 1).
The aim of this paper is to analyze the spectral behavior of an extreme xA source, the quasar HE0359-3959. We performed multicomponent fits in a wide spectral range: UV, optical and Near-Infrared (Section 2); which gives us information about the dynamics and the physical conditions of the broad line region (BLR) (Section 3). In Section 4, we summarize the main results of our work.

Observations and Data Reduction
Ultraviolet (UV), optical, and Near-Infrared spectra were observed with the Very Large Telescope (VLT-ESO). Optical and Near-Infrared spectra were obtained with the Infrared Spectrometer And Array Camera (ISAAC; decommissioned in 2013) using a slit of 0.6 ′′ . The near-infrared spectrum was observed in 2010 in the K band with a total exposure time of 1,120 s. The optical spectrum was observed in 2004 in the J band with a total exposure time of 3,600 s. For the ultraviolet spectrum we used the Focal Reducer and low dispersion Spectrograph (FORS1) and a slit of 1.0 ′′ with a total exposure time of 1,440 s. It was observed in 2008. The data reduction was done using the IRAF package. The procedures followed are explained in Marziani et al. (2009)

Multicomponent Fits
We perform multicomponent fits using SPECFIT, an IRAF routine (Kriss, 1994) to get the information of the most important emission lines. In each spectral range we fit a local continuum. The FWHM of all the broad components (BC) for Hβ, AlIIIλ1860, SiIII]λ1892, CIVλ1549, and SiIVλ1397 was taken equal. In the Figure 2, we present the multicomponents fits after continuum subtraction, for the CIVλ1549 and Ca II triplet range. The rest of the fits will be shown in an upcoming paper.

Multiwavelength Analysis
Low-ionization lines (LIL) have an ionization potential (IP) 20 eV. The Hβ line is the prototype of LIL. In population A3 and A4 sources Hβ has associated a blueshifted component (Bachev et al., 2004). In the case of HE0359-3959, the blueshifted component has a contribution to the total flux of ∼9%, and shows a centroid a half intensity of c(1/2)≈−500 ± 70 km s −1 .
The FeII (IP∼16 eV) has an important contribution in the optical and near-infrared regions. To reproduce it we used the templates modeled by Marziani et al. (2009) and Garcia-Rissmann et al. (2012) for the optical and near-infrared ranges, respectively. Several works have found (Persson, 1988;Ferland and Persson, 1989;Joly, 1989;Dultzin-Hacyan et al., 1999;Martínez-Aldama et al., 2015) a close relationship between the FeII and the NIR Ca II λ8498, λ 8542, and λ8662 Å triplet. This relation is very well appreciable in this object: as well as the optical FeII is strong, the NIR Ca II triplet also is. It is the first time where we observe the Ca II triplet lines isolated at high redshift. Strong intensities of both ions imply an extremely low-ionization degree (U < 10 −2 ; U: ionization parameter) and a high density (n H ∼ 10 11-13 cm −3 ) (Baldwin et al., 2004;Matsuoka et al., 2007;Martínez-Aldama et al., 2015).
In the UV region, the 1900Å blend is formed by two intermediate-ionization lines (IIL; IP ∼ 20-40 eV), AlIIIλ1860 and SiIII]λ1892, which are accompanied by CIII]λ1909 and some FeIII transitions. In this blend we appreciate a blueshifted component. This component should be most likely associated with AlIIIλ1860. Respect to AlIIIλ1860, the blueshifted component has a contribution of the total profile of 60%. The centroid a half intensity is c(1/2)≈−3,200 ± 250 km s −1 , which indicates the presence of an outflow generated by radiation forces presented in the intermediate-ionization lines . On the other hand, considering the high intensity of AlIIIλ1860, SiIII]λ1892, Ca II and FeII, it could suggest a possible chemical enrichment of the BLR (Juarez et al., 2009).

Physical Properties of HE0359-3959
In order to study the physical properties of the quasar HE0359-3959, we built a grid of photoionization simulations using the CLOUDY code (Ferland et al., 1998(Ferland et al., , 2013. For our simulations we considerer a Mattews and Ferland continuum (Mathews and Ferland, 1987), a plane-parallel geometry, a metallicity 5Z⊙ with an overabundance of Al and Si with respect to carbon (by a factor of three), and a column density of N c = 10 23 cm −2 . See Negrete et al. (2012) for more details. Our simulations span the density range 7.00 ≤ log (n H ) ≤ 14.00 and −4.5 ≤log (U) ≤0.00 for the ionization parameter, in intervals of 0.25 dex. More details about the CLOUDY simulations can be found in Negrete et al. (2014). Using the UV lines, we define three groups of diagnostic ratios:: • The flux ratio AlIIIλ1860/SiIII]λ1892 is a useful density diagnostic. • The flux ratio SiIVλ1397/SiIII]λ1892 for the ionization parameter. • The flux ratio CIVλ1549/SiIVλ1397 is mainly sensitive to the relatives abundances of C and Si.
In Figure 3 is shown the result of the simulations. We obtained that the flux ratios are intersected in log(n H ) = 12.32 cm −3 and log(U) = −2.95. Compared to not highly radiating AGNs (Negrete et al., 2013), this source shows a high density and a low ionization parameter, which marks a different behavior in the BLR, probably causing by the slim disk hosted in these kind of sources. Taking into account the high intensity of AlIIIλ1860, FeII, and Ca II we conclude that effectively the lowionization emitter zone has a high density and low-ionization parameter. Negrete et al. (2012) proposed a new method to determine the size of the BLR (r BLR ) and the black hole mass (M BH ) based on the product n H · U and independently of redshift. This method gives similar results to the obtained from the classical methods such as reverberation mapping at low-z (Negrete et al., 2014). Knowing the product of n H · U obtained from the CLOUDY simulations, we compute the size of the BLR (r BLR ) and considering the FWHM of the broad components as the velocity dispersion, we can get the black hole mass (M BH ) and the Eddington ratio. The size of the BLR is log(r BLR ) = 18.37 ± 0.04 cm and the black hole mass is log(M BH ) = 9.52 ± 0.41 M ⊙ . These values are in agreement with the ones found for a large xA sample at high-redshift (Martínez-Aldama et al., in preparation).
The Eddington ratio for this source is L/L Edd = 0.74 ± 0.11. Considering that it shows a c(1/2)∼−6,000 ± 500 km s −1 for CIVλ1549, we confirm the directly proportional relation between c(1/2) and L/L Edd . Indicating that L/L Edd could be the driver of the outflows (Sulentic et al., 2017).

CONCLUSIONS
The information given by the multiwavelength analysis indicates that in HE0359-3959 there is coexistence of substructures in the broad line region. Low and intermediate-ionization regions, where Hβ, AlIIIλ1860 and SiIII]λ1892 are emitted, are dense (n H ∼10 11-12 cm −3 ) and optically thick (U ∼10 −2.5 ). They are mainly governed by virial motions and the presence of a blueshifted component indicates the influence of radiation forces. On the other hand, according to Marziani et al. (2010) the highionization region is less dense (n H ∼10 10 cm −3 , U ∼10 −1 ), pointing out a difference with the physical conditions shown by the low and intermediate-ionization lines.
High ionization lines are dominated by strong radiation forces, producing outflows in high-ionization lines like CIVλ1549, HeIIλ1640, and SiIVλ1397. The high Eddington ratio value suggests the presence of a slim optically thick disk which could be related to the extreme outflow properties observed in HE0359-3959. The presence of strong outflows has been related with the co-evolution of the active galactic nuclei and the host galaxy.