Issue |
A&A
Volume 647, March 2021
|
|
---|---|---|
Article Number | A133 | |
Number of page(s) | 8 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202039375 | |
Published online | 23 March 2021 |
Detection of an ionized gas outflow in the extreme UV-luminous star-forming galaxy BOSS-EUVLG1 at z = 2.47
1
Centro de Astrobiología (CSIC-INTA), Carretera de Ajalvir, 28850 Torrejón de Ardoz, Madrid, Spain
e-mail: javier.alvarez@cab.inta-csic.es
2
International Associate, Cosmic Dawn Center (DAWN), Copenhagen, Denmark
3
Instituto de Astrofísica de Canarias, C/Vía Láctea, s/n, 38205 San Cristóbal de La Laguna, Tenerife, Spain
4
Universidad de La Laguna, Dpto. Astrofísica, 38206 San Cristóbal de La Laguna, Tenerife, Spain
Received:
9
September
2020
Accepted:
28
January
2021
BOSS-EUVLG1 is the most ultraviolet (UV) and Lyα luminous galaxy to be going through a very active starburst phase detected thus far in the Universe. It is forming stars at a rate of 955 ± 118 M⊙ yr−1. We report the detection of a broad Hα component carrying 25% of the total Hα flux. The broad Hα line traces a fast and massive ionized gas outflow characterized by a total mass, log(Mout[M⊙]), of 7.94 ± 0.15, along with an outflowing velocity (Vout) of 573 ± 151 km s−1 and an outflowing mass rate (Ṁout) of 44 ± 20 M⊙ yr−1. The presence of the outflow in BOSS-EUVLG1 is also supported by the identification of blueshifted UV absorption lines in low and high ionization states. The energy involved in the Hα outflow can be explained by the ongoing star formation, without the need for an active galactic nucleus to be included in the scenario. The derived low mass-loading factor (η = 0.05 ± 0.03) indicates that, although it is massive, this phase of the outflow cannot be relevant for the quenching of the star formation, namely, the negative feedback. In addition, only a small fraction (≤15%) of the ionized outflowing material with velocities above 372 km s−1 has the capacity to escape the gravitational potential and to enrich the surrounding circumgalactic medium at distances above several tens of kpc. The ionized phase of the outflow does not carry sufficient mass or energy to play a relevant role in the evolution of the host galaxy nor in the enrichment of the intergalactic medium. As predicted by some recent simulations, other phases of the outflow could be responsible for most of the outflow energy and mass in the form of hot X-ray emitting gas. The expected emission of the extended X-ray emitting halo associated with the outflow in BOSS-EUVLG1 and similar galaxies could be detected with the future ATHENA X-ray observatory, however, there are no methods at present that would assist in their spatial resolution.
Key words: galaxies: starburst / galaxies: high-redshift / Galaxy: kinematics and dynamics / galaxies: ISM / intergalactic medium
© ESO 2021
1. Introduction
Cosmological simulations exploring the formation and evolution of galaxies in the early universe invoke outflows to explain the observed stellar mass function of galaxies. While outflows produced by luminous accreting black holes are favoured for the high-mass end, outflows associated with intense starbursts are claimed for the low-mass end (e.g., Somerville & Davé 2015; Naab & Ostriker 2017; Nelson et al. 2019; Mitchell et al. 2020; Valentini et al. 2020). These outflows do not only control the growth of stellar mass in galaxies across redshifts, they are also responsible for the metal enrichment of the circum- and intergalactic medium (Kim et al. 2020).
Observationally, outflows in galaxies are a universal phenomenon that have been detected at both low (e.g., Arribas et al. 2014; Cicone et al. 2016; Concas et al. 2019) and high redshifts (e.g., Freeman et al. 2019; Förster Schreiber et al. 2019; Swinbank et al. 2019; Finnerty et al. 2020), covering all types of galaxies, ranging from those featuring an active galactic nucleus (AGN; e.g., Harrison et al. 2014, 2018; King & Pounds 2015; Fiore et al. 2017) to those featuring starbursts (e.g., Heckman et al. 1990; Veilleux et al. 2005, 2020). The complex, intrinsic multi-phase nature (molecular, neutral, and ionized) of galactic outflows has been traced at all wavelengths from the X-rays (Strickland et al. 2004), to ultraviolet (Steidel et al. 2010), optical (Arribas et al. 2014; Cazzoli et al. 2016; Cicone et al. 2016; Concas et al. 2019), near-infrared (Hill & Zakamska 2014; Emonts et al. 2017), far-infrared (Veilleux et al. 2013), and millimeter wavelengths (Cicone et al. 2014; Pereira-Santaella et al. 2018; Spilker et al. 2020a,b; Veilleux et al. 2020).
Most of the studies of ionized outflows at high-z have focused on galaxies located on the main sequence (MS) of star-forming galaxies (SFG), namely, galaxies with average star formation rates (SFRs) of several to a few tens of M⊙ yr−1 (Förster Schreiber et al. 2019; Freeman et al. 2019; Swinbank et al. 2019). Ionized gas outflows appear to be frequent at the peak of the cosmic star formation history (i.e. z ∼ 1 − 3). As concluded in the recent KMOS3D Hα survey (Förster Schreiber et al. 2019), about a third of the non-AGN dominated SFGs show evidence of ionized gas outflows with masses (Mout) of 107 − 108 M⊙, outflowing mass rates (Ṁout) of few tens M⊙ yr−1, and low mass-loading factors (η = Ṁout/SFR) of ∼0.1−0.2. Similar mass outflowing rates and mass-loading factors have been measured in MS star-forming galaxies at z ∼ 1 (Swinbank et al. 2019). These studies concluded that the warm ionized phase of the outflow is not sufficient to explain low-mass galaxy formation efficiency, additionally asserting that cold molecular and hot X-ray emitting phases are expected to dominate the mass, energy, and momentum of the outflows as observed in low-z galaxies (e.g., Veilleux et al. 2020; Strickland et al. 2004).
Studies of outflows in high-z galaxies with extreme star formation, that is, with SFRs of several hundreds to thousand M⊙ yr−1, are far more limited. A recent study of massive, dusty star-forming galaxies (DSFG) at z > 4 reported the detection of massive molecular outflows in a large fraction (73%) of the galaxies (Spilker et al. 2020a). These outflows are characterized by high velocities (Vmax ∼ 430 − 1200 km s−1), mass outflowing rates of several M⊙ yr−1, and mass-loading factors just below one (Spilker et al. 2020b). However, the presence of the (partially) ionized outflowing component was not detected in those sources with a high quality [CII]158 μm line and with clear molecular outflow, as traced by the OH119 μm absorption line. High-velocity outflows have also been detected in z > 2 Lyman break galaxies (LBGs, Steidel et al. 2010; Sugahara et al. 2019). These LBGs show blueshifted low- and high-ionization ultraviolet (UV) interstellar medium (ISM) absorption lines tracing outflows with maximum velocities of about 800 km s−1. Therefore, high-velocity outflows appear to be common in both infrared (i.e., dusty) and ultraviolet (i.e., dust-poor) luminous star-forming galaxies at high redshifts. However, the multi-phase structure and their impact on the formation and evolution of galaxies is far from clear. Detailed multi-wavelength studies of individual prototypes, or small samples of galaxies, are required to quantify and establish the relevance of each of the main outflow components (cold dense molecular, warm (partially-) ionized, or hot diffuse X-rays) in the stellar buildup of galaxies and in the metal enrichment of the intergalactic medium (IGM).
Very recently, an extremely UV luminous galaxy, BOSS-EUVLG1, at a redshift of 2.47, was discovered in the study of Marques-Chaves et al. (2020a; hereafter MC20) within the extended Baryon Oscillation Spectroscopic Survey (Abolfathi et al. 2018). BOSS-EUVLG1 is a compact (size ≤ 1.2 kpc) and young (∼4−5 Myr) starburst forming stars at a rate of ∼1000 M⊙ yr−1, with low metallicity and dust content (12 + log(O/H) = 8.13 ± 0.19, E(B − V)≃0.07). BOSS-EUVLG1 is the most UV-luminous galaxy identified so far in the universe and can be considered as the prototype of a new and unique type of extremely UV luminous galaxies (EUVLGs) characterized as having extremely high UV continuum and Lyα luminosities (MUV < −23.5; log(LLyα[erg s−1]) > 43.5), and unobscured SFR above 450 M⊙ yr−1. In addition, BOSS-EUVLG1 has a specific star formation (sSFR) ∼106 ± 15 Gyr−1, a factor 30 times higher than 1010 M⊙ main-sequence galaxies at redshifts 2−2.5 (Whitaker et al. 2015) but in the range of the Hα-excess galaxies at redshifts ∼4−5 that were recently identified in the Spitzer UltraVISTA ultra-deep stripes survey (SMUVS survey; Caputi et al. 2017). An important fraction (15%) of the SMUVS galaxies in the 9.4 ≤ log(Mstellar[M⊙]) ≤ 11.0 mass range appear to be in a starburst phase with SFR of hundreds to thousands M⊙ yr−1 and sSFR well above 24 Gyr−1, and up to 120 Gyr−1. The SMUVS galaxies in the starburst phase account for more than 50% of the cosmic SFR density at redshifts ∼3.9−4.9. Therefore, extreme UV luminous galaxies such as BOSS-EUVLG1 could represent the short lived phase of this class of Hα-excess galaxies, and could be important contributors to the SFR density at intermediate and high redshifts. BOSS-EUVLG1, therefore, brings on the opportunity to investigate the presence and properties of ionized gas outflows for the first time in this class of extreme star-forming galaxies.
In this article, we report the detection of the ionized phase of the outflow in BOSS-EUVLG1 using the Hα emission line. Section 2 describes the observations and calibrations of the Hα spectrum. Section 3 presents the analysis of Hα emission and rest-frame UV absorption lines. The results of the analysis and interpretation as outflowing material are presented in Sect. 4, including the overall kinematic and mass (Sect. 4.1), the outflowing mass rate and mass-loading factor (Sect. 4.2), the nature of the energy source (Sect. 4.3), and the mass escape fraction into the intergalactic medium (Sect. 4.4). The results are discussed in Sect. 5, focusing in the comparison with low- and high-z SFGs (Sect. 5.1), the multi-phase nature of the outflows (Sect. 5.2), and the possibility of detecting these outflows with the future X-ray observatory Advanced Telescope for High ENergy Astrophysics (ATHENA; Sect. 5.3). Finally, Sect. 6 present the summary and outlook. We adopted H0 = 67.7 km s−1 Mpc−1 and Ωm = 0.307 (Planck Collaboration XIII 2016) as cosmological parameters and the Salpeter initial mass function (IMF; Salpeter 1955) throughout this paper. A factor of 1.66 (Newman et al. 2012) is used to convert the SFRs and stellar masses (Mstellar) from a Chabrier (2003) to a Salpeter (1955) IMF when using some published data.
2. Observations
Near-infrared (Near-IR) long-slit spectroscopy of BOSS-EUVLG1 was obtained with the Espectrógrafo Multiobjeto Infra-Rojo (EMIR) instrument on the 10.4 m Gran Telescopio Canarias (GTC). Observations were taken in two service mode nights – 2020 February 10 and May 8 – in bright Moon and subarcsec seeing (≃0.8″ full width half maximum) conditions, as part of the GTC program GTCMULTIPLE2E-19A (PI: R. Marques-Chaves). The observations used the Kgrism with a dispersion of 1.71 Å pix−1 and a 0.8″-wide long-slit, providing a spectral resolution of 82 km s−1 at ∼2.25 μm. The slit was centered on a bright (K = 17.43 ± 0.02, AB) reference star (SDSS J122039.27+084225.0, ≃25″ SW of BOSS-EUVLG1) and oriented so as to encompass BOSS-EUVLG1, with a sky position angle of PA = 58.8° (north toward east). The total exposure time on-source was 6400 s, split into 40 individual exposures of 160 s each, using an ABBA nodding pattern with 6″ offset along the slit.
Data reduction was performed using the EMIR pipeline (PYEMIR)1. Near-IR 2D spectra were flux calibrated using the HIP 59174 telluric standard star observed in the second night. To account for slit-losses, the fluxes of BOSS-EUVLG1 and the reference star were matched to those obtained from photometry. The 1D spectrum of Hα was extracted applying an optimal extraction algorithm (Horne 1986) in an aperture of 3″. Figure 1 shows the final 2D calibrated and 1D extracted spectrum of BOSS-EUVLG1.
![]() |
Fig. 1. Observed Hα emission of BOSS-EUVLG1. Upper panel: near-IR 2D calibrated spectrum centered on the Hα emission. Bottom panel: observed Hα emission line profile of BOSS-EUVLG1, together with the two-components Gaussian fit and its residuals. Gray line: 1D observed Hα spectrum. Black line: best two-component Gaussian Hα fitted model. Red and blue lines: narrow and broad Hα components. Green line: [N II]6548,6583 Å lines. Gray dotted line: illustration of observed spectral line-spread function, derived from a sky line at 2.25 μm, normalized to the peak of the Hα emission, and centered within a velocity of −1200 km s−1. The χ2, calculated from −1000 to 750 km s−1, is also included. |
3. Data analysis
The Hα emission of BOSS-EUVLG1 exhibits an asymmetric line profile characterized by a main component with a blue-shifted wing (Fig. 1). We implement a two-component Gaussian fit, using a Python version of the MPFIT routine (Markwardt 2009)2, to characterize the Hα emission. The [N II]6548,6583 Å doublet and a continuum emission, together with the instrumental broadening, are also included in the fit. Figure 1 shows the results of the two-components fit with the narrow and a blue-shifted broad emission of Hα. Additionally, we test the possibility to fit the profile with just one Gaussian component but it fails to reproduce its asymmetry, yielding to a χ2 50% larger than using the two-component fit. The errors on the parameters are estimated using Monte Carlo. The noise of the spectrum is measured as the root mean square (rms) of the continuum surrounding the Hα+[NII] lines. This noise is used to generate new spectra (3000), where a random Gaussian noise with a sigma equal to the rms is added to the original spectrum before the lines are fitted again.
The Hα narrow component has a flux (FN) of 1.43 ± 0.16 × 10−15 erg s−1 cm−2 and a full width at half maximum (FWHMN) of 241 ± 11 km s−1. The broad component is characterized by an Hα flux (FB) of 4.9 ± 1.7 × 10−16 erg s−1 cm−2, a width (FWHMB) of 511 ± 145 km s−1, and a velocity offset (Vpeak) of −139 ± 87 km s−1 (i.e., blueshift) with respect to the narrow component. In addition, the [N II]6583 Å emission appears as a weak line with a measured flux of (5.9 ± 1.5)×10−17 erg s−1 cm−2 for the best fit. These results are consistent, within the uncertainties, with the slight asymmetries in the low-resolution spectrum of the [O III] 4959,5007 Å lines (R ∼ 735, MC20). Higher spectral resolution is required to constrain the properties of the outflow traced by the high ionization [O III] 4959,5007 Å lines.
Evidence for outflowing material in BOSS-EUVLG1 comes from the analysis of the UV rest-frame absorption line spectrum taken with the GTC (see MC20 for a complete description of the observed spectrum). We use several UV ISM absorption lines in a variety of ionization states, from partial neutral gas (Si II 1260 Å, C II 1334 Å, and Si II 1526 Å) to highly ionized species (Si IV 1393,1402 Å). The rest-frame UV spectrum of BOSS-EUVLG1 is first normalized using spectral windows that are free of absorption and emission features identified by Rix et al. (2004). Next, we derive the average ISM absorption line profiles by performing a mean of the individual absorption lines in the velocity plane. Left panel of Fig. 2 shows the resulting average profiles using low- and high-ionization UV ISM absorption lines and the comparison with the Hα line. As seen in this figure, the centroids of low- and high-ionization ISM absorption lines are clearly blueshifted with respect to the systemic redshift by Δvlow-ion = −305 ± 75 km s−1 and Δvhigh-ion = −415 ± 75 km s−1, respectively. The errors account for the uncertainties of the spectral bin size and the standard deviation of the centroid of each individual ISM line. Due to the low spectral resolution of the rest-frame UV spectrum (R ∼ 700), these lines appear almost unresolved so that the detailed characterization of the outflowing gas from ISM absorption lines (i.e., Ṁout, Vmax, column density, covering fraction, etc.), is not possible as it requires a better spectral resolution. However, they do present further evidence for the existence of outflows in this galaxy and, therefore, lend their support to the interpretation of the broad Hα component as outflowing material.
![]() |
Fig. 2. View of the UV absorption lines and the UV peak velocity of BOSS-EUVLG1. Left panel: comparison between the Hα emission and the UV ISM absorption lines. Upper left panel: normalized Hα spectrum (black line), together with the narrow (red line) and broad (blue line) Hα components. Bottom left panel: UV ISM absorption lines for two different ionization states. Orange line: combination of UV ISM absorption lines of partially neutral elements (Si II 1260 Å, C II 1334 Å, and Si II 1526 Å). Green line: combination of UV ISM absorption lines of high ionization elements (Si IV 1393,1402 Å). Both lines are normalized in a spectral windows clean of absorption and emission features. The peak velocities of both UV absorption ionization states and the broad Hα component are shown as color-coded vertical lines. Right panel: peak velocity of the UV low-ionization ISM absorption lines as a function of the UV absolute magnitude. Yellow star: BOSS-EUVLG1. Blue dots: Lyman alpha emitters (Patrício et al. 2016; Marques-Chaves et al. 2020b). Red squares: Lyman break galaxies (Pettini et al. 2000; Quider et al. 2009, 2010; Dessauges-Zavadsky et al. 2010; Marques-Chaves et al. 2018). |
4. Results
4.1. Kinematics and mass of the ionized outflow
The narrow component of Hα is assigned to the emission in the host galaxy while the broad component is identified with outflowing ionized gas. The main properties of this outflowing material (i.e., outflowing mass, velocity of the bulk of the outflow and its maximal velocity) are derived from the flux, offset velocity, and line profile of the broad Hα line component, assuming the electron density is known. For BOSS-EUVLG1, the electron density is derived from the low-ionization [O II] 3729,3727 Å ratio (0.68 ± 0.21, MC20) giving a density of 1500 cm−3 (Kewley et al. 2019). This value provides the mean electron density of the ionized gas in the galaxy. Since no direct measurement of the electron density in the outflow is available within our data, this value of the density will also be assumed for outflowing material in all the subsequent analysis and derivation of the physical quantities. This is a valid assumption. The median electron density of the outflows traced by the low ionization [SII]6717,6731 Å lines in a large sample of low-z U/LIRGs is a factor of 1.23 larger than the density of the ionized gas in these galaxies (Arribas et al. 2014).
The peak velocity (Vpeak) of the broad Hα component is blueshifted by −139 ± 87 km s−1. High-spatial resolution zoom-in simulations of the ionized gas velocity field traced by the Hα emission line in z ∼ 2 non-interacting galaxies predict the presence of a broad Hα line component in two-third of the simulated galaxies and with blueshifted velocities of 110 ± 20 km s−1 (Ceverino et al. 2016). This result is explained as a direct consequence of the extended structure of the outflowing material, coming from the near-side of the galaxy, above the plane of the galaxy (Ceverino et al. 2016). These simulations cover galaxies with stellar masses 1−4 × 1010 M⊙, similar to that of BOSS-EUVLG1, but with SFR ∼ 20−60 M⊙ yr−1, namely, factors of 16−50 lower than the SFR derived for BOSS-EUVLG1. A similar spatial and velocity structure can be expected in BOSS-EUVLG1 where, due to its higher SFR, the outflowing material could be launched to distances further out, and up to several kpc above the plane of the galaxy.
The observed Hα outflowing velocity is substantially lower than that measured using the UV absorption lines (Δvlow-ion = −305 ± 75 km s−1 and Δvhigh-ion = −415 ± 75 km s−1), but it is in agreement with the typical blueshifted velocity observed in Lyman Break Galaxies at intermediate redshifts (−164 ± 16 km s−1, Steidel et al. 2010). Such differences are expected because the outflowing gas is stratified, that is, regions of different ionization state (neutral versus ionized gas, low- versus high-ionization ISM) show different outflowing velocities (e.g., Erb 2015; Rupke 2018; de la Cruz et al. 2020).
The kinematics of the Hα outflow is characterized by the peak velocity offset and the outflow velocity (Vout = Vpeak + 2 × σB). For a FWHMB of 511 ± 145 km s−1 and a Vpeak of −139 ± 87 km s−1, the Vout correspond to 573 ± 151 km s−1. Following Colina et al. (1991), the total outflowing ionized mass (Mout ∝ L(Hα)/ne) corresponds to log(Mout[M⊙]) = 7.94 ± 0.15, derived from the extinction corrected (E(B − V) = 0.07, MC20) broad Hα luminosity. A more conservative estimate of the total outflowing ionized mass can be derived by assuming only the ionized gas with velocities that are well above the expected dispersion or rotational velocities derived from the narrow Hα component. Considering velocities only above 241 km s−1 (FWHMN) and the offset velocity of the outflowing component, the lower limit to the outflowing mass would correspond to about a third of the total mass. Therefore, the outflowing ionized mass would be in the mass range of 7.42 < log(Mout[M⊙]) < 7.94.
4.2. Outflowing mass rate and mass-loading factor
The outflowing ionized mass rate is given as K × Mout × Vout × size−1, where K is a factor that depends on the geometry, with values between 1 and 3 for an elongated and spherical filled outflow, respectively (Lutz et al. 2020). Here, we adopt a value of one as the spatially resolved outflows in nearby galaxies support in general elongated geometries (Pereira-Santaella et al. 2018) and the limited spatial resolution here does not provide information on the geometry of the outflow. For a Vout of 573 ± 151 km s−1 and an outflow size comparable to the UV-continuum source (i.e., 1.2 kpc, MC20), the outflowing mass rate corresponds to 44 ± 20 M⊙ yr−1.
The SFR derived from the observed Hα flux corrected by the internal extinction (E(B − V) = 0.07, MC20) corresponds to 955 ± 118 M⊙ yr−1 for a Salpeter IMF (Salpeter 1955) and a SFR = 7.9 × 10−42 × L(Hα) (Kennicutt 1998). The corresponding mass-loading factor (η = Ṁout/SFR) is 0.05 ± 0.03. This low value implies that while the amount of ionized gas in the outflow is significant, it does not represent a substantial amount relative to the ongoing star formation process that is taking place. Therefore, the ionized warm phase of the outflow, traced by the Hα emission line, appears irrelevant in the quenching of the extreme on-going star formation in BOSS-EUVLG1. However, other phases could serve as the dominant phases of the outflow, carrying most of its mass and energy, therefore acting as the major contributors to the removal of gas and quenching of star formation (see discussion in Sect. 5.2).
4.3. Energy source of the outflow
For young starbursts, the rates of outflowing mass, kinetic energy, and momentum liberated by supernovae explosions and stellar winds are a function of the SFR. Following Veilleux et al. (2005), these quantities correspond to 248 M⊙ yr−1, 6.7 ± 0.8 × 1044 erg s−1, and 4.8 ± 0.6 × 1036 dynes for a SFR of 955 M⊙ yr−1, respectively. The measured rate of kinetic energy () and momentum (Ṁout × Vout) of the Hα outflow corresponds to 4.6 ± 3.2 × 1042 erg s−1 and 1.6 ± 0.8 × 1035 dynes, respectively. The ionized outflow represents in BOSS-EUVLG1 a minor fraction of the outflowing mass (18%), kinetic energy (0.7%), and momentum (3.3%) rates liberated by the starburst. There is therefore no need to invoke the presence of an AGN in the nucleus to explain the energetics of the warm ionized outflow. However, unless the starburst in BOSS-EUVLG1 is extremely inefficient in converting the available energy into outflows in the surrounding interstellar medium, contributors to the outflow could still be present and locked in other phases such as in the hot, X-ray emitting gas (see Sect. 5.2 for further discussion).
4.4. Ionized mass escape fraction and IGM metal enrichment
Independent of the energy source, gas outflows in galaxies are relevant, not only as potential mechanisms of star formation quenching, but also as one of the main mechanisms for expelling gas from galaxies into the intergalactic medium, with subsequent metal enrichment (Somerville & Davé 2015). An estimate of the fraction of the outflowing material in BOSS-EUVLG1 that is capable of escaping the gravitational potential of the galaxy and enrich the IGM can be derived assuming only the outflowing material with velocities above the escape velocity. For a galaxy represented by an isothermal sphere truncated to a maximum radius (Rmax), the escape velocity for a dynamical mass (Mdyn) within a radius (R) is given as (e.g., Arribas et al. 2014 and references therein). BOSS-EUVLG1 has a dynamical mass of 1.8 × 1010 M⊙, derived for a R of 1.2 kpc and velocity dispersion from the narrow Hα component of 102 km s−1 (Bellocchi et al. 2013). For a Rmax of 10 kpc, the escape velocity corresponds to 372 km s−1.
The velocity of the outflow (Vout = 573 km s−1) is 1.5 times larger than the escape velocity. Thus, only a small fraction of the total ionized outflowing material (15%, i.e., ∼1.3 × 107 M⊙), with velocities well above Vesc, would potentially be able to escape and enrich the surrounding IGM at distances larger than 10 kpc. Most of the outflowing material would not escape the gravity of the host galaxy and would be likely to fall back into the system, providing the raw material for subsequent delayed star formation episodes, as in low-z infrared luminous star-forming galaxies (Arribas et al. 2014; Ceverino et al. 2016). We note, however, that Hα only traces the warm (104 K) ionized gas component of the outflow. In fact, the peak velocities of the blueshifted low- and high-ionization UV ISM absorption lines (Δvlow-ion ∼ −305 km s−1 and Δvhigh-ion ∼ −415 km s−1) are factors of 2 and 3 larger than the corresponding Hα (Vpeak = −139 ± 85 km s−1). This already demonstrates the existence of outflowing material that is rich in metals, with velocities higher than those of the Hα emitting gas – velocities that are high enough to escape into the IGM and contribute to its metal enrichment.
5. Discussion
5.1. Ionized Outflows in BOSS-EUVLG1, low- and intermediate-z star-forming galaxies
BOSS-EUVLG1 is an extreme starburst galaxy with a star-forming rate of 955 ± 118 M⊙ yr−1, which can be considered an ultraluminous galaxy based on its UV luminosity (i.e., L1500 Å ∼ 1012 L⊙). Its SFR is similar to those measured in high-z dust-enshrouded extreme starbursts (Reuter et al. 2020), and about 10−100 times higher than that of massive star-forming galaxies at intermediate redshifts (Förster Schreiber et al. 2019; Swinbank et al. 2019; Freeman et al. 2019; Newman et al. 2012), and low-z luminous and ultraluminous infrared galaxies (U/LIRGs, Arribas et al. 2014). Moreover, BOSS-EUVLG1 is deficient in terms of dust and is underluminous in the infrared (log(LIR[L⊙]) < 10.91), yielding a IR/UV luminosity ratio of log(LIR/LUV) < − 1.2. The low extinction and metallicity values measured in BOSS-EUVLG1 suggest that this galaxy is likely in an early and short-lived evolutionary stage as it continues to build its stellar mass, providing metal and dust enrichment, before evolving into the infrared luminous dust-enshrouded starburst phase (MC20).
The observed properties of the outflow in BOSS-EUVLG1 are presented in Fig. 3 together with those of low- and intermediate-z star-forming galaxies using the same tracer (Hα). The ratio of the Hα broad- to narrow-component flux (FB/FN) in BOSS-EUVLG1 is 0.34 ± 0.12. This value is lower than the ratios measured in low- and intermediate-z star-forming galaxies where the broad line component carries a larger fraction of the flux with typical FB/FN values in the 0.4 to 0.8 range, independent of the SFR albeit with a large scatter (upper left panel of Fig. 3). However, the FB/FN ratio as a function of the stellar mass follows the expected trend for star-forming galaxies (upper right panel of Fig. 3). These results clearly indicate that the ionized outflow component in BOSS-EUVLG1 is low for its specific star formation rate (sSFR = 106 ± 15 Gyr−1, MC20). BOSS-EUVLG1 is an extreme starburst galaxy with sSFR about 30−50 times higher than that of Main Sequence star-forming galaxies at z ∼ 2 (Rodighiero et al. 2014). There is no clear explanation for the low contribution of the outflowing component to the total Hα emission. It appears as if the strength of the outflow is dominated by the stellar mass of the system – and not by the strength of the star formation as given by the SFR. This is somehow unexpected given that in a starburst-dominated outflow, the mass, momentum, and energy rates of the outflow are directly related to the total star formation rate (Heckman et al. 1990; Veilleux et al. 2005). In addition, outflows in z ∼ 2 star-forming galaxies are characterized by high FB/FN values close to 1 for all galaxies with SFR surface brightness above 1 M⊙ yr−1 kpc−2 and low values of about 0.2 for the galaxies below the SFR brightness threshold, independent of the stellar mass (Newman et al. 2012). According to its sSFR, SFR, and upper limit on the size of the Hα emission, BOSS-EUVLG1 would be classified in the category of galaxies with strong outflows characterized by a FB/FN ratio of about 1. However, the measured value of 0.34 is closer to that of weak outflows in starbursts with a low-SFR surface brightness. Explanations of the low FB/FN value in BOSS-EUVLG1 could invoke a qualitatively different multi-phase structure and composition of the interstellar medium and outflows in IR-luminous, namely, dustier, more molecular, and metal-rich; and in UV-luminous starbursts such as BOSS-EUVLG1, namely, dust-, molecular-, and metal-poor (see discussion in Sects. 5.2 and 5.3). The Spitzer UltraVISTA ultra-deep stripes survey (SMUVS survey; Caputi et al. 2017) has identified a number IR-luminous, Hα-excess galaxies at redshifts of ∼4−5 with extreme sSFR. An important fraction (15%) of the SMUVS galaxies in the 9.4 ≤ log(Mstellar[M⊙]) ≤ 11.0 mass range appear to be in a starburst phase with SFR of hundreds to thousands M⊙ yr−1 and sSFR well above 24 Gyr−1, and up to 120 Gyr−1. Future Hα spectroscopy of EUVLGs and SMUVS galaxies in the starburst phase from the ground and with JWST (for z > 3) are required to confirm the strength of outflows in extreme starbursts (i.e., very large sSFR) and investigate it in depth to identify the potential differences between the classes of IR- and UV-luminous.
![]() |
Fig. 3. Observational properties of the ionized outflow in BOSS-EUVLG1 (yellow star) together with samples of star-forming galaxies (SFGs) at intermediate redshifts including 0.6 < z < 2.7 SFGs (red dots) from the KMOS3D survey (Förster Schreiber et al. 2019, values of Vout are a priv. comm.), SFGs at z ∼ 1 (magenta dots) from KMOS observations (Swinbank et al. 2019), SFGs at z ∼ 2 (green dots) from the MOSDEF survey Freeman et al. (2019), and SFGs at z∼2 (blue dots) from the SINS and zC-SINF surveys (Newman et al. 2012). Low-z (ultra-)luminous galaxies (cyan and orange dots) are also represented (Arribas et al. 2014). Upper panels: flux ratio of the broad to narrow Hα line components (FB/FN) as a function of the SFR (left) and Mstellar (right). Bottom panels: outflow velocity as a function of the SFR (left) and Mstellar (right). The continuum line in the Vout to SFR relation defines the best fit (Vout = 148 × SFR0.23) to the BOSS-EUVLG1, low- and intermediate-z data points. The dashed line represents the Vout to Hα-derived SFR relation derived for the low-z sample of U/LIRGs without AGNs (Arribas et al. 2014). |
The outflow velocity in BOSS-EUVlG1 is in agreement with the expected value as extrapolated from previous Vout ∝ SFRα relations derived in low-z U/LIRGs (α = 0.21−0.24, Arribas et al. 2014 and references therein). This ratio is also consistent with the dependence observed in the ionized outflows of intermediate-z star-forming galaxies (lower-left panel of Fig. 3). The best fit, which includes BOSS-EUVLG1 as well as low- and intermediate-z galaxies (i.e., covering a factor of 400 in SFR) gives a slope (α) of 0.23, which is consistent with the predicted dependence with the bolometric luminosity of the starburst (Vout ∝ , Heckman et al. 2000). However, the outflow velocity in BOSS-EUVLG1, 573 km s−1, stands out against the measured weak dependence of the outflowing velocities as a function of the stellar mass. Velocities, Vout, lower than 400 km s−1 are measured in galaxies with log(Mstellar[M⊙]) of 10.0 (lower right panel of Fig. 3). In summary, the observed properties of the ionized outflowing gas in BOSS-EUVLG1 appear to be a consequence of its extremely high sSFR, that is, a very massive, young, and compact starburst in a galaxy with relatively low stellar mass. On the one hand, the Vout follows the functional dependence observed in star-forming galaxies at low and intermediate redshifts. On the other hand, the amount of outflowing material, as identified by its low FB/FN, appears to be consistent with the ratio expected by its low stellar mass.
The blueshifts identified in the UV ISM absorption lines of BOSS-EUVLG1 are also indicative of large-scale galactic outflows, similar to that observed in many other star-forming galaxies at low (e.g., Heckman et al. 2000; Rubin et al. 2014) and high redshift (e.g., Shapley et al. 2003; Steidel et al. 2010). The right panel of Fig. 2 shows a comparison between the blueshift of the low-ionization ISM absorption lines (relative to the systemic redshift) and the UV absolute magnitude MUV of BOSS-EUVLG1 and other UV-luminous star-forming galaxies at 2 < z < 3 (most of them strongly magnified and where Δv can be measured on an individual basis; Pettini et al. 2000; Quider et al. 2009, 2010; Dessauges-Zavadsky et al. 2010; Patrício et al. 2016; Marques-Chaves et al. 2018, 2020b). This figure clearly shows a trend among Δvlow-ion and MUV, indicating that UV-luminous star-forming galaxies, such as BOSS-EUVLG1, experience outflows caused by supernova explosions and stellar winds, with the peak velocity of these outflows increasing with the absolute UV luminosity of the galaxy.
A more direct comparison of other relevant physical quantities of the ionized outflow, such as the rates of outflowing mass, kinetic energy, and momentum, as well as the mass-loading factors, are difficult to establish. These quantities have a strong dependence on the value adopted for the electron density and size of the Hα emitting nebulae. In some samples, such as the MOSDEF survey of z ∼ 2 of star-forming galaxies (Freeman et al. 2019), a fixed radius of 3 kpc and an electron density of 50 cm−3 is adopted for the outflowing material. Other studies of outflows in mean sequence star-forming galaxies at z ∼ 1 (Swinbank et al. 2019) adopt a radius of 10 kpc and a electron density of 75 cm−3 for the composite spectra of different galaxies of the sample. Finally, the recently completed KMOS-3D survey of z ∼ 0.6 − 2.7 (Förster Schreiber et al. 2019) derived an electron density of 380 cm−3 from the stack spectra of the star-forming driven outflows. The electron densities measured or adopted in all these intermediate redshift samples are factors of ∼4−30 lower than the value (1500 cm−3) measured in BOSS-EUVLG1, while the sizes are consistently larger by factors of 3−10, making any comparison between the ionized outflow in BOSS-EUVLG1 with other galaxy samples rather uncertain in absolute terms. However, if the differences in the size and electron density were real, these would indicate that the outflow in BOSS-EUVLG1 is much more compact and involves a much denser gas than in mean sequence star-forming galaxies at intermediate redshifts. Accurate measurements in both IR-luminous and UV-luminous star-forming are needed to establish whether or not these two population of galaxies have intrinsically different physical properties in their outflows.
5.2. The multi-phase nature of outflows
The outflows predicted in cosmological simulations of high-z galaxies are extended over regions of several kpc in size (effective radius of 1−7 kpc) but with maximal velocities of less than 200 km s−1 (Ceverino et al. 2016), that is, a factor 2−3 smaller than the Vout measured in BOSS-EUVLG1. However, according to most recent high angular resolution (100−200 pc) hydrodynamical simulations (TNG50; Nelson et al. 2019), outflowing maximal velocities tend to increase with the stellar mass of the galaxy, such that Vmax of 500−600 km s−1 are expected in z ∼ 2 − 3 galaxies with stellar masses of about 1010 M⊙. According to these simulations, the mass-loading factor would be about 4 for a typical 1010 M⊙ galaxy at z ∼ 2, namely, a factor of about 100 larger than the value derived for BOSS-EUVLG1 (η = 0.05 ± 0.03). This apparent discrepancy could be explained by the multi-phase nature of the outflows in galaxies. The Hα line used in this study traces only the warm 104 K optical emitting gas, while it is known that some outflows carry also relevant amounts of mass, momentum, and energy in their molecular phase (Veilleux et al. 2020), and in the diffuse hot X-ray emitting gas (e.g., M 82; Strickland & Heckman 2009; Lopez et al. 2020; also Heckman & Thompson 2017 for a review). In fact, according to the TNG50 simulations, the ouflowing mass rate for galaxies with masses of 3 × 1010 M⊙, or less, is dominated by the diffuse X-ray emitting gas with the peak of the outflowing mass at temperatures above 106 K (Nelson et al. 2019). The X-ray emitting mass ouflowing rate in these galaxies is predicted to be at least a factor of 10 higher than that of the warm emitting gas (Nelson et al. 2019). It could well be that the Hα outflow identified in BOSS-EUVLG1 could be tracing only the warm, lower velocity phase of a much powerful outflow with a total mass-loading factor (i.e., all velocities and gas phases) close to one or even higher (Nelson et al. 2019) and, therefore, relevant for the quenching or temporal delay of the star formation in these galaxies.
5.3. Prospects for detecting X-ray emitting outflowing gas with ATHENA
ATHENA, the future European Space Agency X-ray satellite, could be able to detect the X-ray emission associated with the massive starburst at the redshift of BOSS-EUVLG1. One of its instruments, the Wide Field Imager (WFI, Meidinger et al. 2018), will be able to detect sources with a flux of 1 × 10−16 erg s−1 cm−2 in the SOFT X-ray (0.5−2 keV) with a signal-to-noise ratio of 5σ in 10 h (Athena Community Office, priv. comm.). At the redshift of BOSS-EUVLG1, this limit corresponds to a luminosity LX-SOFT = 5 × 1042 erg s−1. Well-known nearby starbursts such as M 82 and NGC 3079 have extended X-ray emitting halos of several kpc (Strickland et al. 2004) with LX-SOFT of 4.1 × 1039 erg s−1 and 8.3 × 1039 erg s−1, and total integrated (i.e., halo plus disk and nucleus) LX-SOFT of 4.3 × 1040 erg s−1 and 3.2 × 1040 erg s−1, respectively. These galaxies are forming stars at a rate of 8−9 M⊙ yr−1 (Strickland et al. 2004), about hundred times less than in BOSS-EUVLG1. Therefore, if the efficiency in the conversion of the radiative and mechanical energy of the starburst into heating of the interstellar medium in BOSS-EUVLG1 is similar to that of M 82 and NGC 3079, a LX-SOFT of 3−4 × 1042 erg s−1 for the integrated emission would be expected, and therefore could be detected with ATHENA. This luminosity has to be considered as a lower limit to the total X-ray emission if the ratio of the Hα to SOFT (0.3−2 keV) X-ray luminosities (LHα/LX-SOFT) for the extended emission in starbursts (Colina et al. 2004) is taken into account. The LHα/LX-SOFT ratio covers a wide range of values from about 0.25 (NGC 253), to 0.5 (NGC 3079, Arp 220), and 6 (M 82). Applying similar ratios to BOSS-EUVLG1, a total SOFT X-ray luminosity in the 8−180 × 1042 erg s−1 would be expected for the outflowing ionized gas with an Hα luminosity of 3.1 × 1043 erg s−1. However, the imaging capabilities of the WFI instrument precludes any direct, spatially resolved, detection of the extended emission associated with the outflowing material. WFI will be providing X-ray images over a 40 × 40 arcmin2 field-of-view with an angular resolution (FWHM) of about 5 arcsec (Meidinger et al. 2018), that is, about 40 kpc at the redshift of BOSS-EUVLG1. Our Hα spectroscopic observations place an upper limit of a few kpc to the size of the extended emission. In addition, X-ray emitting halos is nearby starbursts have sizes of a few kpc (M 82, NGC 253, NGC 3079) to about 25 kpc (Arp 220). In summary, the detection with a high significance of the X-ray emission in starburst galaxies similar to BOSS-EUVLG1 at redshifts of 2−3 will be possible in the future with the ATHENA satellite, but the X-ray emitting halos associated with the starburst-generated outflowing material will not be spatially resolved.
6. Summary and outlook
We report the detection of a massive and fast ionized gas outflow in BOSS-EUVLG1, the most UV luminous star-forming galaxy detected so far in the Universe. BOSS-EUVLG1 is the brightest member of a new class of galaxies, the Extremely UV Luminous Galaxies (EUVLGs), with (unobscured) star formation rates above 450 M⊙ yr−1, recently detected in the BOSS survey. The Hα line in emission has been used to characterized the ionized outflow of BOSS-EUVLG1. Hα presents a broad component blueshifted by −139 ± 87 km s−1 and with a width (FWHM) of 511 ± 145 km s−1. The presence of the outflow in BOSS-EUVLG1 is also supported by blueshifted UV ISM absorption lines. The Hα outflow is described by a broad-to-narrow component flux ratio (FB/FN) of 0.34 ± 0.12, an outflowing velocity (Vout) of 573 ± 151 km s−1, a total mass, log(Mout[M⊙]), of 7.94 ± 0.15, and an outflowing mass rate (Ṁout) of 44 ± 20 M⊙ yr−1. The velocity of the outflow in BOSS-EUVLG1 follows the expected value in the Vout-SFR diagram predicted from the relation derived for low- and intermediate-z star-forming galaxies. The fraction of the flux carried out by the outflow is within the value expected for its stellar mass in the FB/FN ratio to stellar mass diagram.
The momentum and energy involved in the ionized gas outflow do not require the presence of an additional energy source other than the starburst itself, such as an AGN. The velocity of the outflow (Vout) is such (573 km s−1) that only a small fraction (≤15%) of the outflowing mass would be able to escape and enrich the surrounding circumgalactic medium. The mass-loading factor is very low (η = 0.05 ± 0.03), indicating that, although it is massive, the ionized phase of the outflow does not play a relevant role in the quenching of the star formation. Other phases, such as the X-ray emitting gas, might be responsible for most of the outflowing mass and energy, as predicted by some recent high angular-resolution cosmological simulations. This will require confirmation in the future thanks to sensitive X-ray observations with Athena. The properties of the ionized outflow in the class of extreme UV-luminous galaxies (EUVLGs) recently identified in the BOSS survey will contribute to characterizing it further, using the Hα emission line as a tracer.
MPFIT Python version used: https://github.com/segasai/astrolibpy/tree/master/mpfit
Acknowledgments
The authors thank to the anonymous referee for useful comments, to Natascha Forster-Schreiber for kindly sharing the data from the KMOS-3D survey used in the plots, and to Francisco Carrera and the ATHENA Community Office for providing us with the expected sensitivity for the future WFI/ATHENA instrument. This work is based on observations made with the Gran Telescopio Canarias (GTC) installed in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, in the island of La Palma. This work was supported by the Spanish State Research Agency (AEI) under grants ESP2015-65597-C4-4-R, ESP2017-83197, ESP2017-86852-C4-2-R, PID2019-106280GB-I00, and MDM-2017-0737 Unidad de Excelencia “María de Maeztu”- Centro de Astrobiología (CSIC-INTA).
References
- Abolfathi, B., Aguado, D. S., Aguilar, G., et al. 2018, ApJS, 235, 42 [NASA ADS] [CrossRef] [Google Scholar]
- Arribas, S., Colina, L., Bellocchi, E., Maiolino, R., & Villar-Martín, M. 2014, A&A, 568, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bellocchi, E., Arribas, S., Colina, L., & Miralles-Caballero, D. 2013, A&A, 557, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Caputi, K. I., Deshmukh, S., Ashby, M. L. N., et al. 2017, ApJ, 849, 45 [NASA ADS] [CrossRef] [Google Scholar]
- Cazzoli, S., Arribas, S., Maiolino, R., & Colina, L. 2016, A&A, 590, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ceverino, D., Arribas, S., Colina, L., et al. 2016, MNRAS, 460, 2731 [Google Scholar]
- Chabrier, G. 2003, ApJ, 586, L133 [NASA ADS] [CrossRef] [Google Scholar]
- Cicone, C., Maiolino, R., Sturm, E., et al. 2014, A&A, 562, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cicone, C., Maiolino, R., & Marconi, A. 2016, A&A, 588, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Colina, L., Sparks, W. B., & Macchetto, F. 1991, ApJ, 370, 102 [NASA ADS] [CrossRef] [Google Scholar]
- Colina, L., Arribas, S., & Clements, D. 2004, ApJ, 602, 181 [NASA ADS] [CrossRef] [Google Scholar]
- Concas, A., Popesso, P., Brusa, M., Mainieri, V., & Thomas, D. 2019, A&A, 622, A188 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de la Cruz, L. M., Schneider, E. E., Ostriker, E. C., et al. 2020, ArXiv e-prints [arXiv:2012.02212] [Google Scholar]
- Dessauges-Zavadsky, M., D’Odorico, S., Schaerer, D., et al. 2010, A&A, 510, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Emonts, B. H. C., Colina, L., Piqueras-López, J., et al. 2017, A&A, 607, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Erb, D. K. 2015, Nature, 523, 169 [NASA ADS] [CrossRef] [Google Scholar]
- Finnerty, L., Larson, K., Soifer, B. T., et al. 2020, ApJ, 905, 16 [Google Scholar]
- Fiore, F., Feruglio, C., Shankar, F., et al. 2017, A&A, 601, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Förster Schreiber, N. M., Übler, H., Davies, R. L., et al. 2019, ApJ, 875, 21 [NASA ADS] [CrossRef] [Google Scholar]
- Freeman, W. R., Siana, B., Kriek, M., et al. 2019, ApJ, 873, 102 [Google Scholar]
- Harrison, C. M., Alexander, D. M., Mullaney, J. R., & Swinbank, A. M. 2014, MNRAS, 441, 3306 [Google Scholar]
- Harrison, C. M., Costa, T., Tadhunter, C. N., et al. 2018, Nat. Astron., 2, 198 [NASA ADS] [CrossRef] [Google Scholar]
- Heckman, T. M., & Thompson, T. A. 2017, ArXiv e-prints [arXiv:1701.09062] [Google Scholar]
- Heckman, T. M., Armus, L., & Miley, G. K. 1990, ApJS, 74, 833 [NASA ADS] [CrossRef] [Google Scholar]
- Heckman, T. M., Lehnert, M. D., Strickland, D. K., & Armus, L. 2000, ApJS, 129, 493 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Hill, M. J., & Zakamska, N. L. 2014, MNRAS, 439, 2701 [NASA ADS] [CrossRef] [Google Scholar]
- Horne, K. 1986, PASP, 98, 609 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kennicutt, R. C., Jr. 1998, ApJ, 498, 541 [NASA ADS] [CrossRef] [Google Scholar]
- Kewley, L. J., Nicholls, D. C., Sutherland, R., et al. 2019, ApJ, 880, 16 [NASA ADS] [CrossRef] [Google Scholar]
- Kim, C.-G., Ostriker, E. C., Somerville, R. S., et al. 2020, ApJ, 900, 61 [Google Scholar]
- King, A., & Pounds, K. 2015, ARA&A, 53, 115 [NASA ADS] [CrossRef] [Google Scholar]
- Lopez, L. A., Mathur, S., Nguyen, D. D., Thompson, T. A., & Olivier, G. M. 2020, ApJ, 904, 152 [Google Scholar]
- Lutz, D., Sturm, E., Janssen, A., et al. 2020, A&A, 633, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Markwardt, C. B. 2009, in Astronomical Data Analysis Software and Systems XVIII, eds. D. A. Bohlender, D. Durand, & P. Dowler, ASP Conf. Ser., 411, 251 [Google Scholar]
- Marques-Chaves, R., Pérez-Fournon, I., Gavazzi, R., et al. 2018, ApJ, 854, 151 [NASA ADS] [CrossRef] [Google Scholar]
- Marques-Chaves, R., Álvarez-Márquez, J., Colina, L., et al. 2020a, MNRAS, 499, L105 [Google Scholar]
- Marques-Chaves, R., Pérez-Fournon, I., Shu, Y., et al. 2020b, MNRAS, 492, 1257 [Google Scholar]
- Meidinger, N., Nandra, K., & Plattner, M. 2018, in Space Telescopes and Instrumentation 2018: Ultraviolet to Gamma Ray, eds. J. W. A. den Herder, S. Nikzad, & K. Nakazawa, SPIE Conf. Ser., 10699, 106991F [Google Scholar]
- Mitchell, P. D., Schaye, J., Bower, R. G., & Crain, R. A. 2020, MNRAS, 494, 3971 [CrossRef] [Google Scholar]
- Naab, T., & Ostriker, J. P. 2017, ARA&A, 55, 59 [Google Scholar]
- Nelson, D., Pillepich, A., Springel, V., et al. 2019, MNRAS, 490, 3234 [NASA ADS] [CrossRef] [Google Scholar]
- Newman, S. F., Genzel, R., Förster-Schreiber, N. M., et al. 2012, ApJ, 761, 43 [NASA ADS] [CrossRef] [Google Scholar]
- Patrício, V., Richard, J., Verhamme, A., et al. 2016, MNRAS, 456, 4191 [NASA ADS] [CrossRef] [Google Scholar]
- Pereira-Santaella, M., Colina, L., García-Burillo, S., et al. 2018, A&A, 616, A171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pettini, M., Steidel, C. C., Adelberger, K. L., Dickinson, M., & Giavalisco, M. 2000, ApJ, 528, 96 [NASA ADS] [CrossRef] [Google Scholar]
- Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Quider, A. M., Pettini, M., Shapley, A. E., & Steidel, C. C. 2009, MNRAS, 398, 1263 [NASA ADS] [CrossRef] [Google Scholar]
- Quider, A. M., Shapley, A. E., Pettini, M., Steidel, C. C., & Stark, D. P. 2010, MNRAS, 402, 1467 [NASA ADS] [CrossRef] [Google Scholar]
- Reuter, C., Vieira, J. D., Spilker, J. S., et al. 2020, ApJ, 902, 78 [NASA ADS] [CrossRef] [Google Scholar]
- Rix, S. A., Pettini, M., Leitherer, C., et al. 2004, ApJ, 615, 98 [NASA ADS] [CrossRef] [Google Scholar]
- Rodighiero, G., Renzini, A., Daddi, E., et al. 2014, MNRAS, 443, 19 [NASA ADS] [CrossRef] [Google Scholar]
- Rubin, K. H. R., Prochaska, J. X., Koo, D. C., et al. 2014, ApJ, 794, 156 [NASA ADS] [CrossRef] [Google Scholar]
- Rupke, D. 2018, Galaxies, 6, 138 [NASA ADS] [CrossRef] [Google Scholar]
- Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
- Shapley, A. E., Steidel, C. C., Pettini, M., & Adelberger, K. L. 2003, ApJ, 588, 65 [NASA ADS] [CrossRef] [Google Scholar]
- Somerville, R. S., & Davé, R. 2015, ARA&A, 53, 51 [Google Scholar]
- Spilker, J. S., Phadke, K. A., Aravena, M., et al. 2020a, ApJ, 905, 85 [Google Scholar]
- Spilker, J. S., Aravena, M., Phadke, K. A., et al. 2020b, ApJ, 905, 86 [Google Scholar]
- Steidel, C. C., Erb, D. K., Shapley, A. E., et al. 2010, ApJ, 717, 289 [NASA ADS] [CrossRef] [Google Scholar]
- Strickland, D. K., & Heckman, T. M. 2009, ApJ, 697, 2030 [NASA ADS] [CrossRef] [Google Scholar]
- Strickland, D. K., Heckman, T. M., Colbert, E. J. M., Hoopes, C. G., & Weaver, K. A. 2004, ApJS, 151, 193 [NASA ADS] [CrossRef] [Google Scholar]
- Sugahara, Y., Ouchi, M., Harikane, Y., et al. 2019, ApJ, 886, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Swinbank, A. M., Harrison, C. M., Tiley, A. L., et al. 2019, MNRAS, 487, 381 [CrossRef] [Google Scholar]
- Valentini, M., Murante, G., Borgani, S., et al. 2020, MNRAS, 491, 2779 [Google Scholar]
- Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 43, 769 [NASA ADS] [CrossRef] [Google Scholar]
- Veilleux, S., Meléndez, M., Sturm, E., et al. 2013, ApJ, 776, 27 [NASA ADS] [CrossRef] [Google Scholar]
- Veilleux, S., Maiolino, R., Bolatto, A. D., & Aalto, S. 2020, A&ARv, 28, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Whitaker, K. E., Franx, M., Bezanson, R., et al. 2015, ApJ, 811, L12 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
![]() |
Fig. 1. Observed Hα emission of BOSS-EUVLG1. Upper panel: near-IR 2D calibrated spectrum centered on the Hα emission. Bottom panel: observed Hα emission line profile of BOSS-EUVLG1, together with the two-components Gaussian fit and its residuals. Gray line: 1D observed Hα spectrum. Black line: best two-component Gaussian Hα fitted model. Red and blue lines: narrow and broad Hα components. Green line: [N II]6548,6583 Å lines. Gray dotted line: illustration of observed spectral line-spread function, derived from a sky line at 2.25 μm, normalized to the peak of the Hα emission, and centered within a velocity of −1200 km s−1. The χ2, calculated from −1000 to 750 km s−1, is also included. |
In the text |
![]() |
Fig. 2. View of the UV absorption lines and the UV peak velocity of BOSS-EUVLG1. Left panel: comparison between the Hα emission and the UV ISM absorption lines. Upper left panel: normalized Hα spectrum (black line), together with the narrow (red line) and broad (blue line) Hα components. Bottom left panel: UV ISM absorption lines for two different ionization states. Orange line: combination of UV ISM absorption lines of partially neutral elements (Si II 1260 Å, C II 1334 Å, and Si II 1526 Å). Green line: combination of UV ISM absorption lines of high ionization elements (Si IV 1393,1402 Å). Both lines are normalized in a spectral windows clean of absorption and emission features. The peak velocities of both UV absorption ionization states and the broad Hα component are shown as color-coded vertical lines. Right panel: peak velocity of the UV low-ionization ISM absorption lines as a function of the UV absolute magnitude. Yellow star: BOSS-EUVLG1. Blue dots: Lyman alpha emitters (Patrício et al. 2016; Marques-Chaves et al. 2020b). Red squares: Lyman break galaxies (Pettini et al. 2000; Quider et al. 2009, 2010; Dessauges-Zavadsky et al. 2010; Marques-Chaves et al. 2018). |
In the text |
![]() |
Fig. 3. Observational properties of the ionized outflow in BOSS-EUVLG1 (yellow star) together with samples of star-forming galaxies (SFGs) at intermediate redshifts including 0.6 < z < 2.7 SFGs (red dots) from the KMOS3D survey (Förster Schreiber et al. 2019, values of Vout are a priv. comm.), SFGs at z ∼ 1 (magenta dots) from KMOS observations (Swinbank et al. 2019), SFGs at z ∼ 2 (green dots) from the MOSDEF survey Freeman et al. (2019), and SFGs at z∼2 (blue dots) from the SINS and zC-SINF surveys (Newman et al. 2012). Low-z (ultra-)luminous galaxies (cyan and orange dots) are also represented (Arribas et al. 2014). Upper panels: flux ratio of the broad to narrow Hα line components (FB/FN) as a function of the SFR (left) and Mstellar (right). Bottom panels: outflow velocity as a function of the SFR (left) and Mstellar (right). The continuum line in the Vout to SFR relation defines the best fit (Vout = 148 × SFR0.23) to the BOSS-EUVLG1, low- and intermediate-z data points. The dashed line represents the Vout to Hα-derived SFR relation derived for the low-z sample of U/LIRGs without AGNs (Arribas et al. 2014). |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.