Showing 1–50 of 107 results for author: Efstathiou, A

SMART: spectral energy distributions Markov chain analysis with radiative transfer models

Authors: Charalambia Varnava, Andreas Efstathiou

Abstract: In this paper we present the publicly available open-source spectral energy distribution (SED) fitting code SMART (Spectral energy distributions Markov chain Analysis with Radiative Transfer models). Implementing a Bayesian Markov chain Monte Carlo (MCMC) method, SMART fits the ultraviolet to millimetre SEDs of galaxies exclusively with radiative transfer models that currently constitute four type… ▽ More In this paper we present the publicly available open-source spectral energy distribution (SED) fitting code SMART (Spectral energy distributions Markov chain Analysis with Radiative Transfer models). Implementing a Bayesian Markov chain Monte Carlo (MCMC) method, SMART fits the ultraviolet to millimetre SEDs of galaxies exclusively with radiative transfer models that currently constitute four types of pre-computed libraries, which describe the starburst, active galactic nucleus (AGN) torus, host galaxy and polar dust components. An important novelty of SMART is that, although it fits SEDs exclusively with radiative transfer models, it takes comparable time to popular energy balance methods to run. Here we describe the key features of SMART and test it by fitting the multi-wavelength SEDs of the 42 local ultraluminous infrared galaxies (ULIRGs) that constitute the HERschel Ultraluminous Infrared Galaxy Survey (HERUS) sample. The Spitzer spectroscopy data of the HERUS ULIRGs are included in the fitting at a spectral resolution, which is matched to that of the radiative transfer models. We also present other results that highlight the performance and versatility of SMART. SMART promises to be a useful tool for studying galaxy evolution in the JWST era. SMART is developed in PYTHON and is available at https://github.com/ch-var/SMART.git. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Reconstructing the mid-infrared spectra of galaxies using ultraviolet to submillimeter photometry and Deep Generative Networks

Authors: Agapi Rissaki, Orestis Pavlou, Dimitris Fotakis, Vicky Papadopoulou Lesta, Andreas Efstathiou

Abstract: The mid-infrared spectra of galaxies are rich in features such as the Polycyclic Aromatic Hydrocarbon (PAH) and silicate dust features which give valuable information about the physics of galaxies and their evolution. For example they can provide information about the relative contribution of star formation and accretion from a supermassive black hole to the power output of galaxies. However, the… ▽ More The mid-infrared spectra of galaxies are rich in features such as the Polycyclic Aromatic Hydrocarbon (PAH) and silicate dust features which give valuable information about the physics of galaxies and their evolution. For example they can provide information about the relative contribution of star formation and accretion from a supermassive black hole to the power output of galaxies. However, the mid-infrared spectra are currently available for a very small fraction of galaxies that have been detected in deep multi-wavelength surveys of the sky. In this paper we explore whether Deep Generative Network methods can be used to reconstruct mid-infrared spectra in the 5-35μm range using the limited multi-wavelength photometry in ~20 bands from the ultraviolet to the submillimeter which is typically available in extragalactic surveys. For this purpose we use simulated spectra computed with a combination of radiative transfer models for starbursts, active galactic nucleus (AGN) tori and host galaxies. We find that our method using Deep Generative Networks, namely Generative Adversarial Networks and Generative Latent Optimization models, can efficiently produce high quality reconstructions of mid-infrared spectra in ~70% of the cases. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: Published in Astronomy and Computing (Volume 47, April 2024, 100823)

Journal ref: Astronomy and Computing 47 (2024) 100823

Structures Of Dust and gAs (SODA): Constraining the innermost dust properties of II Zw96 with JWST observations of H$_2$O and CO

Authors: I. García-Bernete, M. Pereira-Santaella, E. González-Alfonso, D. Rigopoulou, A. Efstathiou, F. R. Donnan, N. Thatte

Abstract: We analyze JWST NIRSpec$+$MIRI/MRS observations of the infrared (IR) gas-phase molecular bands of the most enshrouded source (D1) within the interacting system and luminous IR galaxy II Zw 096. We report the detection of rovibrational lines of H$_2$O $ν_2$=1-0 ($\sim$5.3-7.2 $μ$m) and $^{12}$CO $ν$=1-0 ($\sim$4.45-4.95 $μ$m) in D1. The CO band shows the R- and P-branches in emission and the spectr… ▽ More We analyze JWST NIRSpec$+$MIRI/MRS observations of the infrared (IR) gas-phase molecular bands of the most enshrouded source (D1) within the interacting system and luminous IR galaxy II Zw 096. We report the detection of rovibrational lines of H$_2$O $ν_2$=1-0 ($\sim$5.3-7.2 $μ$m) and $^{12}$CO $ν$=1-0 ($\sim$4.45-4.95 $μ$m) in D1. The CO band shows the R- and P-branches in emission and the spectrum of the H$_2$O band shows the P-branch in emission and the R-branch in absorption. The H$_2$O R-branch in absorption unveils an IR-bright embedded compact source in D1 and the CO broad component features a highly turbulent environment. From both bands, we also identified extended intense star-forming (SF) activity associated with circumnuclear photodissociation regions (PDRs), consistent with the strong emission of the ionised 7.7 $μ$m polycyclic aromatic hydrocarbon band in this source. By including the 4.5-7.0 $μ$m continuum information derived from the H$_2$O and CO analysis, we modelled the IR emission of D1 with a dusty torus and SF component. The torus is very compact (diameter of $\sim$3 pc at 5 $μ$m) and characterised by warm dust ($\sim$ 370 K), giving an IR surface brightness of $\sim$3.6$\times$10$^{8}$ L$_{\rm sun}$/pc$^2$. This result suggests the presence of a dust-obscured active galactic nucleus (AGN) in D1, which has an exceptionally high covering factor that prevents the direct detection of AGN emission. Our results open a new way to investigate the physical conditions of inner dusty tori via modelling the observed IR molecular bands. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A (official date of acceptance: 27/12/2023; submitted 27/11/2023). 11 pages, 10 figures

The Galaxy Activity, Torus, and Outflow Survey (GATOS). III: Revealing the inner icy structure in local AGN

Authors: I. García-Bernete, A. Alonso-Herrero, D. Rigopoulou, M. Pereira-Santaella, T. Shimizu, R. Davies, F. R. Donnan, P. F. Roche, O. González-Martín, C. Ramos Almeida, E. Bellocchi, P. Boorman, F. Combes, A. Efstathiou, D. Esparza-Arredondo, S. García-Burillo, E. González-Alfonso, E. K. S. Hicks, S. Hönig, A. Labiano, N. A. Levenson, E. López-Rodríguez, C. Ricci, C. Packham, D. Rouan , et al. (2 additional authors not shown)

Abstract: We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1… ▽ More We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1 $μ$m). All the nuclei show the 9.7 $μ$m silicate band in absorption. We compare the strength of the 9.7 and 18 $μ$m silicate features with torus model predictions. The observed silicate features are generally well explained by clumpy and smooth torus models. We report the detection of the 6 $μ$m dirty water ice band (i.e., a mix of water and other molecules such as CO and CO$_2$) at sub-arcsecond scales ($\sim$0.26 arcsec at 6 $μ$m; inner $\sim$50 pc) in a sample of local AGN with different levels of nuclear obscuration in the range log N$_{\rm H}^{\rm X-Ray}$(cm$^{-2}$)$\sim22-25$. We find a good correlation between the 6 $μ$m water ice optical depths and N$_{\rm H}^{\rm X-Ray}$. This result indicates that the water ice absorption might be a reliable tracer of the nuclear intrinsic obscuration in AGN. The weak water ice absorption in less obscured AGN (log N$_H^{X-ray}$ (cm$^{-2}$)$\lesssim$23.0 cm$^{-2}$) might be related to the hotter dust temperature ($>$T$_{sub}^{H_2O}\sim$110 K) expected to be reached in the outer layers of the torus due to their more inhomogeneous medium. Our results suggest it might be necessary to include the molecular content, such as, H$_2$O, aliphatic hydrocarbons (CH-) and more complex PAH molecules in torus models to better constrain key parameters such as the torus covering factor (i.e. nuclear obscuration). △ Less

Submitted 4 December, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A. 13 pages, 12 Figures

Graph Theoretical Analysis of local ultraluminous infrared galaxies and quasars

Authors: Orestis Pavlou, Ioannis Michos, Vicky Papadopoulou Lesta, Michalis Papadopoulos, Evangelos S. Papaefthymiou, Andreas Efstathiou

Abstract: We present a methodological framework for studying galaxy evolution by utilizing Graph Theory and network analysis tools. We study the evolutionary processes of local ultraluminous infrared galaxies (ULIRGs) and quasars and the underlying physical processes, such as star formation and active galactic nucleus (AGN) activity, through the application of Graph Theoretical analysis tools. We extract, p… ▽ More We present a methodological framework for studying galaxy evolution by utilizing Graph Theory and network analysis tools. We study the evolutionary processes of local ultraluminous infrared galaxies (ULIRGs) and quasars and the underlying physical processes, such as star formation and active galactic nucleus (AGN) activity, through the application of Graph Theoretical analysis tools. We extract, process and analyse mid-infrared spectra of local (z < 0.4) ULIRGs and quasars between 5-38 microns through internally developed Python routines, in order to generate similarity graphs, with the nodes representing ULIRGs being grouped together based on the similarity of their spectra. Additionally, we extract and compare physical features from the mid-IR spectra, such as the polycyclic aromatic hydrocarbons (PAHs) emission and silicate depth absorption features, as indicators of the presence of star-forming regions and obscuring dust, in order to understand the underlying physical mechanisms of each evolutionary stage of ULIRGs. Our analysis identifies five groups of local ULIRGs based on their mid-IR spectra, which is quite consistent with the well established fork classification diagram by providing a higher level classification. We demonstrate how graph clustering algorithms and network analysis tools can be utilized as unsupervised learning techniques for revealing direct or indirect relations between various galaxy properties and evolutionary stages, which provides an alternative methodology to previous works for classification in galaxy evolution. Additionally, our methodology compares the output of several graph clustering algorithms in order to demonstrate the best-performing Graph Theoretical tools for studying galaxy evolution. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: Accepted for publication in Astronomy and Computing

Journal ref: Astronomy and Computing 45C (2023) 100742

Observational evidence for cosmological coupling of black holes and its implications for an astrophysical source of dark energy

Authors: Duncan Farrah, Kevin S. Croker, Gregory Tarlé, Valerio Faraoni, Sara Petty, Jose Afonso, Nicolas Fernandez, Kurtis A. Nishimura, Chris Pearson, Lingyu Wang, Michael Zevin, David L Clements, Andreas Efstathiou, Evanthia Hatziminaoglou, Mark Lacy, Conor McPartland, Lura K Pitchford, Nobuyuki Sakai, Joel Weiner

Abstract: Observations have found black holes spanning ten orders of magnitude in mass across most of cosmic history. The Kerr black hole solution is however provisional as its behavior at infinity is incompatible with an expanding universe. Black hole models with realistic behavior at infinity predict that the gravitating mass of a black hole can increase with the expansion of the universe independently of… ▽ More Observations have found black holes spanning ten orders of magnitude in mass across most of cosmic history. The Kerr black hole solution is however provisional as its behavior at infinity is incompatible with an expanding universe. Black hole models with realistic behavior at infinity predict that the gravitating mass of a black hole can increase with the expansion of the universe independently of accretion or mergers, in a manner that depends on the black hole's interior solution. We test this prediction by considering the growth of supermassive black holes in elliptical galaxies over $0<z\lesssim2.5$. We find evidence for cosmologically coupled mass growth among these black holes, with zero cosmological coupling excluded at 99.98% confidence. The redshift dependence of the mass growth implies that, at $z\lesssim7$, black holes contribute an effectively constant cosmological energy density to Friedmann's equations. The continuity equation then requires that black holes contribute cosmologically as vacuum energy. We further show that black hole production from the cosmic star formation history gives the value of $Ω_Λ$ measured by Planck while being consistent with constraints from massive compact halo objects. We thus propose that stellar remnant black holes are the astrophysical origin of dark energy, explaining the onset of accelerating expansion at $z \sim 0.7$. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 10 pages, 3 figures, published in ApJ Letters

Journal ref: ApJL 944 L31 (2023)

A Preferential Growth Channel for Supermassive Black Holes in Elliptical Galaxies at z<2

Authors: Duncan Farrah, Sara Petty, Kevin Croker, Gregory Tarle, Michael Zevin, Evanthia Hatziminaoglou, Francesco Shankar, Lingyu Wang, David L Clements, Andreas Efstathiou, Mark Lacy, Kurtis A. Nishimura, Jose Afonso, Chris Pearson, Lura K Pitchford

Abstract: The assembly of stellar and supermassive black hole (SMBH) mass in elliptical galaxies since $z\sim1$ can help to diagnose the origins of locally-observed correlations between SMBH mass and stellar mass. We therefore construct three samples of elliptical galaxies, one at $z\sim0$ and two at $0.7\lesssim z \lesssim2.5$, and quantify their relative positions in the $M_{BH}-M_*$ plane. Using a Bayesi… ▽ More The assembly of stellar and supermassive black hole (SMBH) mass in elliptical galaxies since $z\sim1$ can help to diagnose the origins of locally-observed correlations between SMBH mass and stellar mass. We therefore construct three samples of elliptical galaxies, one at $z\sim0$ and two at $0.7\lesssim z \lesssim2.5$, and quantify their relative positions in the $M_{BH}-M_*$ plane. Using a Bayesian analysis framework, we find evidence for translational offsets in both stellar mass and SMBH mass between the local sample and both higher redshift samples. The offsets in stellar mass are small, and consistent with measurement bias, but the offsets in SMBH mass are much larger, reaching a factor of seven between $z\sim1$ and $z\sim0$. The magnitude of the SMBH offset may also depend on redshift, reaching a factor of $\sim20$ at $z\sim 2$. The result is robust against variation in the high and low redshift samples and changes in the analysis approach. The magnitude and redshift evolution of the offset are challenging to explain in terms of selection and measurement biases. We conclude that either there is a physical mechanism that preferentially grows SMBHs in elliptical galaxies at $z\lesssim 2$, or that selection and measurement biases are both underestimated, and depend on redshift. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: ApJ, accepted

Classification of local ultraluminous infrared galaxies and quasars with kernel principal component analysis

Authors: Evangelos S. Papaefthymiou, Ioannis Michos, Orestis Pavlou, Vicky Papadopoulou Lesta, Andreas Efstathiou

Abstract: We present a new diagnostic diagram for local ultraluminous infrared galaxies (ULIRGs) and quasars, analysing particularly the Spitzer Space Telescope's Infrared Spectrograph (IRS) spectra of 102 local ULIRGs and 37 Palomar Green quasars. Our diagram is based on a special non-linear mapping of these data, employing the Kernel Principal Component Analysis method. The novelty of this map lies in the… ▽ More We present a new diagnostic diagram for local ultraluminous infrared galaxies (ULIRGs) and quasars, analysing particularly the Spitzer Space Telescope's Infrared Spectrograph (IRS) spectra of 102 local ULIRGs and 37 Palomar Green quasars. Our diagram is based on a special non-linear mapping of these data, employing the Kernel Principal Component Analysis method. The novelty of this map lies in the fact that it distributes the galaxies under study on the surface of a well-defined ellipsoid, which, in turn, links basic concepts from geometry to physical properties of the galaxies. Particularly, we have found that the equatorial direction of the ellipsoid corresponds to the evolution of the power source of ULIRGs, starting from the pre-merger phase, moving through the starburst-dominated coalescing stage towards the active galactic nucleus (AGN)-dominated phase, and finally terminating with the post-merger quasar phase. On the other hand, the meridian directions distinguish deeply obscured power sources of the galaxies from unobscured ones. These observations have also been verified by comparison with simulated ULIRGs and quasars using radiative transfer models. The diagram correctly identifies unique galaxies with extreme features that lie distinctly away from the main distribution of the galaxies. Furthermore, special two-dimensional projections of the ellipsoid recover almost monotonic variations of the two main physical properties of the galaxies, the silicate and PAH features. This suggests that our diagram naturally extends the well-known Spoon diagram and it can serve as a diagnostic tool for existing and future infrared spectroscopic data, such as those provided by the James Webb Space Telescope. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Published in Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society 517 3 (2022) 4162-4174

Detecting and analysing the topology of the cosmic web with spatial clustering algorithms I: Methods

Authors: Dimitrios Kelesis, Spyros Basilakos, Vicky Papadopoulou Lesta, Dimitris Fotakis, Andreas Efstathiou

Abstract: In this paper we explore the use of spatial clustering algorithms as a new computational approach for modeling the cosmic web. We demonstrate that such algorithms are efficient in terms of computing time needed. We explore three distinct spatial methods which we suitably adjust for (i) detecting the topology of the cosmic web and (ii) categorizing various cosmic structures as voids, walls, cluster… ▽ More In this paper we explore the use of spatial clustering algorithms as a new computational approach for modeling the cosmic web. We demonstrate that such algorithms are efficient in terms of computing time needed. We explore three distinct spatial methods which we suitably adjust for (i) detecting the topology of the cosmic web and (ii) categorizing various cosmic structures as voids, walls, clusters and superclusters based on a variety of topological and physical criteria such as the physical distance between objects, their masses and local densities. The methods explored are (1) a new spatial method called Gravity Lattice ; (2) a modified version of another spatial clustering algorithm, the ABACUS; and (3) the well known spatial clustering algorithm HDBSCAN. We utilize HDBSCAN in order to detect cosmic structures and categorize them using their overdensity. We demonstrate that the ABACUS method can be combined with the classic DTFE method to obtain similar results in terms of the achieved accuracy with about an order of magnitude less computation time. To further solidify our claims, we draw insights from the computer science domain and compare the quality of the results with and without the application of our method. Finally, we further extend our experiments and verify their effectiveness by showing their ability to scale well with different cosmic web structures that formed at different redshifts. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Stellar and black hole assembly in z<0.3 infrared-luminous mergers: intermittent starbursts vs. super-Eddington accretion

Authors: Duncan Farrah, Andreas Efstathiou, Jose Afonso, Jeronimo Bernard-Salas, Joe Cairns, David L Clements, Kevin Croker, Evanthia Hatziminaoglou, Maya Joyce, Mark Lacy, Vianney Lebouteiller, Alix Lieblich, Carol Lonsdale, Seb Oliver, Chris Pearson, Sara Petty, Lura K Pitchford, Dimitra Rigopoulou, Michael Rowan-Robinson, Jack Runburg, Henrik Spoon, Aprajita Verma, Lingyu Wang

Abstract: We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show i… ▽ More We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show increased black hole mass to stellar mass ratios, consistent with the transition from AGN to ellipticals and implying substantial black hole mass growth over the course of the merger. Star formation rates are enhanced relative to the local main sequence, by factors of ~100 in the starburst and ~1.8 in the host. The starburst star formation rates appear distinct to star formation in the main sequence at all redshifts up to at least z~5. Starbursts may prefer late-stage mergers, but are observed at any merger stage. We do not find evidence that the starbursts in these low-redshift systems substantially increase the total stellar mass, with a soft upper limit on the stellar mass increase from starburst activity of about a factor of two. In contrast, 12 objects show evidence for super-Eddington accretion, associated with late-stage mergers, suggesting that many AGN in infrared-luminous mergers go through a super-Eddington phase. The super-Eddington phase may increase black hole mass by up to an order of magnitude at an accretion efficiency of 42+/-33% over a period of 44+/-22Myr. Our results imply that super-Eddington accretion is an important black hole growth channel in infrared-luminous galaxies at all redshifts. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: MNRAS Accepted

A new look at local ultraluminous infrared galaxies: the atlas and radiative transfer models of their complex physics

Authors: A. Efstathiou, D. Farrah, J. Afonso, D. L. Clements, E. González-Alfonso, M. Lacy, S. Oliver, V. Papadopoulou Lesta, C. Pearson, D. Rigopoulou, M. Rowan-Robinson, H. W. W. Spoon, A. Verma, L. Wang

Abstract: We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo (MCMC) code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (A… ▽ More We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo (MCMC) code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (AGN) tori and host galaxy. The Spitzer IRS spectroscopy data are included in the fitting. Our bayesian SED fitting method takes comparable time to popular energy balance methods but it is more physically motivated and versatile. All HERUS galaxies harbor high rates of star formation but we also find bolometrically significant AGN in all of the galaxies of the sample. We estimate the correction of the luminosities of the AGN in the ULIRGs due to the anisotropic emission of the torus and find that it could be up to a factor of $\sim10$ for nearly edge-on tori. We present a comparison of our results with the smooth torus model of Fritz et al. and the two-phase models of Siebenmorgen et al. and SKIRTOR. We find that the CYGNUS AGN torus models fit significantly better the SEDs of our sample compared to all other models. We find no evidence that strong AGN appear either at the beginning or end of a starburst episode or that starbursts and AGN affect each other. IRAS 01003-2238 and Mrk 1014 show evidence for dual AGN in their SED fits suggesting a minimum dual AGN fraction in the sample of 5%. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: Published in Monthly Notices of the Royal Astronomical Society. The CYGNUS models used in this paper are available at https://arc.euc.ac.cy/cygnus/

Journal ref: 2022MNRAS.512.5183E

A technique to select the most obscured galaxy nuclei

Authors: I. García-Bernete, D. Rigopoulou, S. Aalto, H. W. W. Spoon, A. Hernán-Caballero, A. Efstathiou, P. F. Roche, S. König

Abstract: Compact obscured nuclei (CONs) are mainly found in local U/LIRGs. In the local Universe, these sources are generally selected through the detection of the HCN-vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and continuum ratios. Using Spitzer/IRS… ▽ More Compact obscured nuclei (CONs) are mainly found in local U/LIRGs. In the local Universe, these sources are generally selected through the detection of the HCN-vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and continuum ratios. Using Spitzer/IRS spectra of a representative sample of local ULIRGs (z<0.27), we examine their PAH and underlying continuum emission ratios. For deeply embedded sources, we find that the 9.7 micron silicate absorption band has a particularly pronounced effect on the 11.3 micron PAH feature. The low flux level in the nuclear silicate absorption band enhances the 11.3 micron PAH feature contrast (high PAH equivalent width) compared to that of the other PAH features. The technique has been extended to include the use of the continuum ratios. However, the latter are affected both by the extinction coming from the host galaxy as well as the nuclear region, whereas the foreground extinction is cancelled out when using the PAH equivalent width ratios. We apply our method to the HERUS and GOALS samples and classify as CON candidates 14 ULIRGs and 10 LIRGs, corresponding to 30% of ULIRGs and 7% of LIRGs from these samples. We find that the observed continuum ratios of CON-dominated sources can be explained by assuming torus models with a tapered disk geometry and a smooth dust distribution. This suggests that the nuclear dusty structure of CONs has an extremely high dust coverage. We also demonstrate that the use of mid-IR color-color diagrams is an effective way to select CON-dominated sources at different redshifts. In particular, the combination of filters of the JWST/MIRI will enable the selection of CONs out to z~1.5. This will allow extending the selection of CONs to high redshifts where U/LIRGs are more numerous. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 13 pages, 11 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 663, A46 (2022)

Energetic nuclear transients in luminous and ultraluminous infrared galaxies

Authors: T. M. Reynolds, S. Mattila, E. Kankare, A. Efstathiou, E. Kool, S. Ryder, L. Peña-Moñino, M. A. Pérez-Torres

Abstract: Energetic nuclear outbursts have been discovered in luminous and ultraluminous infrared galaxies (U/LIRGs) at unexpectedly high rates. To investigate this population of transients, we performed a search in mid-IR data from the Wide-field Infrared Survey Explorer (WISE) satellite and its NEOWISE survey to detect and characterise luminous and smoothly evolving transients in a sample of 215 U/LIRGs.… ▽ More Energetic nuclear outbursts have been discovered in luminous and ultraluminous infrared galaxies (U/LIRGs) at unexpectedly high rates. To investigate this population of transients, we performed a search in mid-IR data from the Wide-field Infrared Survey Explorer (WISE) satellite and its NEOWISE survey to detect and characterise luminous and smoothly evolving transients in a sample of 215 U/LIRGs. We report three new transients, all with $ΔL > 10^{43}$ erg s$^{-1}$, in addition to two previously known cases. Their host galaxies are all part of major galaxy mergers, and through radiative transfer model fitting we find that all have a significant contribution from an active galactic nucleus (AGN). We characterised the transients through measurements of their luminosities and resulting energetics, all of which are between 10$^{50.9}$ erg and 10$^{52.2}$ erg. The IR emission of the five transients was found to be consistent with re-radiation by the hot dust of emission at shorter wavelengths, presumably originating from an accretion event, onto the supermassive black hole. The corresponding transient rate of (1.6-4.6)$\times$10$^{-3}$ / yr / galaxy is over an order of magnitude higher than the rate of large amplitude flares shown by AGN in the optical. We suggest that the observed transients are part of a dust-obscured population of tidal disruption events (TDEs) that have remained out of the reach of optical surveys due to the obscuring dust. In one case, this is supported by our radio observations. We also discuss other plausible explanations. The observed rate of events is significantly higher than optical TDE rates, which can be expected in U/LIRG hosts undergoing a major galaxy merger with increased stellar densities in the nuclear regions. Continued searches for such transients and their multi-wavelength follow-up is required to constrain their rate and nature. △ Less

Submitted 11 October, 2022; v1 submitted 8 February, 2022; originally announced February 2022.

Comments: 15 pages; 8 figures, published by A&A. Re-uploaded to reflect published version. Main changes from previous version are expanded radio analysis and slightly altered rate due to additional year of NEOWISE data

Journal ref: A&A 664, A158 (2022)

Radio Spectra of Luminous, Heavily Obscured WISE-NVSS Selected Quasars

Authors: Pallavi Patil, Mark Whittle, Kristina Nyland, Carol Lonsdale, Mark Lacy, Amy E. Kimball, Colin Lonsdale, Wendy Peters, Tracy E. Clarke, Andreas Efstathiou, Simona Giacintucci, Minjin Kim, Lauranne Lanz, Dipanjan Mukherjee, Emil Polisensky

Abstract: We present radio spectra spanning $0.1 - 10$ GHz for the sample of heavily obscured luminous quasars with extremely red mid-infrared-optical colors and compact radio emission. The spectra are constructed from targeted 10 GHz observations and archival radio survey data, which together yield $6-11$ flux density measurements for each object. Our suite of Python tools for modeling the radio spectra is… ▽ More We present radio spectra spanning $0.1 - 10$ GHz for the sample of heavily obscured luminous quasars with extremely red mid-infrared-optical colors and compact radio emission. The spectra are constructed from targeted 10 GHz observations and archival radio survey data, which together yield $6-11$ flux density measurements for each object. Our suite of Python tools for modeling the radio spectra is publicly available on Github. Our primary result is that most (61%) of the sample have peaked or curved radio spectra and many (36%) could be classified as Gigahertz Peaked Spectrum (GPS) sources. This indicates compact emission regions likely arising from recently triggered radio jets. Assuming synchrotron self-absorption (SSA) generates the peaks, we infer compact source sizes ($3 - 100$ pc) with strong magnetic fields ($6 - 100$ mG) and young ages ($30 - 10^4$ years). Conversely, free-free absorption (FFA) could also create peaks due to the high column densities associated with the deeply embedded nature of the sample. However, we find no correlations between the existence or frequency of the peaks and any parameters of the MIR emission. The high-frequency spectral indices are steep ($α\approx -1$) and correlate, weakly, with the ratio of MIR photon energy density to magnetic energy density, suggesting that the spectral steepening could arise from inverse Compton scattering off the intense MIR photon field. This study provides a foundation for combining multi-frequency and mixed-resolution radio survey data for understanding the impact of young radio jets on the ISM and star formation rates of their host galaxies. △ Less

Submitted 21 July, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 48 pages, 17 figures, published in Astrophysical Journal

The viewing angle in AGN SED models, a data-driven analysis

Authors: Andrés Felipe Ramos Padilla, Lingyu Wang, Katarzyna Małek, Andreas Efstathiou, Guang Yang

Abstract: The validity of the unified active galactic nuclei (AGN) model has been challenged in the last decade, especially when different types of AGNs are considered to only differ in the viewing angle to the torus. We aim to assess the importance of the viewing angle in classifying different types of Seyfert galaxies in spectral energy distribution (SED) modelling. We retrieve photometric data from publi… ▽ More The validity of the unified active galactic nuclei (AGN) model has been challenged in the last decade, especially when different types of AGNs are considered to only differ in the viewing angle to the torus. We aim to assess the importance of the viewing angle in classifying different types of Seyfert galaxies in spectral energy distribution (SED) modelling. We retrieve photometric data from publicly available astronomical databases: CDS and NED, to model SEDs with X-CIGALE in a sample of 13 173 Seyfert galaxies located at redshift range from $z=0$ to $z=3.5$, with a median redshift of $z\approx0.2$. We assess whether the estimated viewing angle from the SED models reflects different Seyfert classifications. Two AGN models with either a smooth or clumpy torus structure are adopted in this paper. We find that the viewing angle in Type-1 AGNs is better constrained than in Type-2 AGNs. Limiting the viewing angles representing these two types of AGNs do not affect the physical parameter estimates such as star-formation rate (SFR) or AGN fractional contribution ($f_{\rm{AGN}}$). In addition, the viewing angle is not the most discriminating physical parameter to differentiate Seyfert types. We suggest that the observed and intrinsic AGN disc luminosity can: i) be used in $z<0.5$ studies to distinguish between Type-1 and Type-2 AGNs, and ii) explain the probable evolutionary path between these AGN types. Finally, we propose the use of X-CIGALE for AGN galaxy classification tasks. All data from the 13 173 SED fits are available at https://doi.org/10.5281/zenodo.5221764 △ Less

Submitted 29 November, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: 19 pages, 10 figures, 6 tables (plus appendix). Accepted for publication in MNRAS

The Nature of Hyperluminous Infrared Galaxies

Authors: F. Gao, L. Wang, A. Efstathiou, K. Małek, P. N. Best, M. Bonato, D. Farrah, R. Kondapally, I. McCheyne, H. J. A. Röttgering

Abstract: We make use of multi-wavelength data of a large hyperluminous infrared (HLIRG) sample to derive their main physical properties, e.g., stellar mass, star-formation rate (SFR), volume density, contribution to the cosmic stellar mass density and to the cosmic SFR density. We also study the black hole (BH) growth rate and its relationship with the SFR of the host galaxy. We select 526 HLIRGs in three… ▽ More We make use of multi-wavelength data of a large hyperluminous infrared (HLIRG) sample to derive their main physical properties, e.g., stellar mass, star-formation rate (SFR), volume density, contribution to the cosmic stellar mass density and to the cosmic SFR density. We also study the black hole (BH) growth rate and its relationship with the SFR of the host galaxy. We select 526 HLIRGs in three deep fields (Bo$ö$tes, Lockman-Hole, ELAIS-N1) and adopt two spectral energy distribution (SED) fitting codes, CIGALE, which assumes energy balance, and CYGNUS, which is based on radiative transfer models and does not adopt energy balance principle. We use two different active galactic nucleus (AGN) models in CIGALE and three AGN models in CYGNUS to compare the results estimated using different SED fitting codes and different AGN models. The stellar mass, total IR luminosity and AGN luminosity agree well between different models with a typical median offset of 0.1 dex. The SFR estimates show the largest dispersions (up to 0.5 dex). This dispersion has an impact on the subsequent analysis, which may suggest that previous contradictory results could partly be due to different choices of methods. HLIRGs are ultra-massive galaxies with 99% of them having stellar masses larger than $10^{11} M_{\odot}$. Our results reveal a higher space density of ultra-massive galaxies than found in previous surveys or predicted by simulations. We find that HLIRGs contribute more to the cosmic SFR density as redshift increases. In terms of BH growth, the two SED fitting methods provide different results. We can see a clear trend in which SFR decreases as AGN luminosity increases when using CYGNUS estimates, possibly implying quenching by AGN, while this trend is much weaker when using CIGALE estimates. This difference is also influenced by the dispersion between SlFR estimates obtained by the two codes. △ Less

Submitted 2 August, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 21 pages, 15 figures. Accepted in A&A

Journal ref: A&A 654, A117 (2021)

HELP: The Herschel Extragalactic Legacy Project

Authors: R. Shirley, K. Duncan, M. C. Campos Varillas, P. D. Hurley, K. Malek, Y. Roehlly, M. W. L. Smith, H. Aussel, T. Bakx, V. Buat, D. Burgarella, N. Christopher, S. Duivenvoorden, S. Eales, A. Efstathiou, E. A. Gonzalez Solares, M. Griffin, M. Jarvis, B. Lo Faro, L. Marchetti, I. McCheyne, A. Papadopoulos, K. Penner, E. Pons, M. Prescott , et al. (7 additional authors not shown)

Abstract: We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the H… ▽ More We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or "open" methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenisation of calibration, meta data and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg2 and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP database. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions. △ Less

Submitted 12 May, 2021; originally announced May 2021.

A hyperluminous obscured quasar at a redshift of z ~ 4.3

Authors: Andreas Efstathiou, Katarzyna Malek, Denis Burgarella, Peter Hurley, Seb Oliver, Veronique Buat, Raphael Shirley, Steven Duivenvoorden, Vicky Papadopoulou Lesta, Duncan Farrah, Kenneth J. Duncan, Maria del Carmen Campos Varillas

Abstract: In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies a… ▽ More In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) multi-component radiative transfer models. We find that its emission is dominated by an obscured quasar with a predicted total 1-1000um luminosity of $3.91^{+1.69}_{-0.55} \times 10^{13} L_\odot$ and an active galactic nucleus (AGN) fraction of ~89%. We also fit HELP_J100156.75+022344.7 with the Code Investigating GALaxy Emission (CIGALE) code and find a similar result. This is only the second z > 4 hyperluminous obscured quasar discovered to date. The discovery of HELP_J100156.75+022344.7 in the ~ 2deg^2 COSMOS field implies that a large number of obscured hyperluminous quasars may lie in the HELP fields which cover ~ 1300deg^2. If this is confirmed, tension between supermassive black hole evolution models and observations will be alleviated. We estimate the space density of objects like HELP_J100156.75+022344.7 at z ~ 4.5 to be $\sim 1.8 \times 10^{-8}$Mpc$^{-3}$. This is slightly higher than the space density of coeval hyperluminous optically selected quasars suggesting that the obscuring torus in z > 4 quasars may have a covering factor $\gtrsim 50\%$. △ Less

Submitted 8 April, 2022; v1 submitted 14 March, 2021; originally announced March 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters. The CYGNUS models used in this work are available at https://arc.euc.ac.cy/cygnus-project-arc/

Journal ref: 2021MNRAS503L..11E

Core-collapse supernova subtypes in luminous infrared galaxies

Authors: E. Kankare, A. Efstathiou, R. Kotak, E. C. Kool, T. Kangas, D. O'Neill, S. Mattila, P. Vaisanen, R. Ramphul, M. Mogotsi, S. D. Ryder, S. Parker, T. Reynolds, M. Fraser, A. Pastorello, E. Cappellaro, P. A. Mazzali, P. Ochner, L. Tomasella, M. Turatto, J. Kotilainen, H. Kuncarayakti, M. A. Perez-Torres, Z. Randriamanakoto, C. Romero-Canizales , et al. (12 additional authors not shown)

Abstract: We present the discovery and follow-up observations of two CCSNe that occurred in the luminous infrared galaxy (LIRG), NGC3256. The first, SN2018ec, was discovered using the ESO HAWK-I/GRAAL adaptive optics seeing enhancer, and was classified as a Type Ic with a host galaxy extinction of $A_V=2.1^{+0.3}_{-0.1}$ mag. The second, AT2018cux, was discovered during the course of follow-up observations… ▽ More We present the discovery and follow-up observations of two CCSNe that occurred in the luminous infrared galaxy (LIRG), NGC3256. The first, SN2018ec, was discovered using the ESO HAWK-I/GRAAL adaptive optics seeing enhancer, and was classified as a Type Ic with a host galaxy extinction of $A_V=2.1^{+0.3}_{-0.1}$ mag. The second, AT2018cux, was discovered during the course of follow-up observations of SN2018ec, and is consistent with a sub-luminous Type IIP classification with an $A_V=2.1 \pm 0.4$ mag of host extinction. A third CCSN, PSNJ10275082-4354034 in NGC3256, has previously been reported in 2014, and we recovered the source in late time archival HST imaging. Based on template light-curve fitting, we favour a Type IIn classification for it with modest host galaxy extinction of $A_V=0.3^{+0.4}_{-0.3}$ mag. We also extend our study with follow-up data of the recent Type IIb SN2019lqo and Type Ib SN2020fkb that occurred in the LIRG system Arp299 with host extinctions of $A_V=2.1^{+0.1}_{-0.3}$ and $A_V=0.4^{+0.1}_{-0.2}$ mag, respectively. Motivated by the above, we inspected, for the first time, a sample of 29 CCSNe located within a projected distance of 2.5 kpc from the host galaxy nuclei in a sample of 16 LIRGs. We find that, if star formation within these galaxies is modelled assuming a global starburst episode and normal IMF, there is evidence of a correlation between the starburst age and the CCSN subtype. We infer that the two subgroups of 14 H-poor (Type IIb/Ib/Ic/Ibn) and 15 H-rich (Type II/IIn) CCSNe have different underlying progenitor age distributions, with the H-poor progenitors being younger at 3$σ$ significance. However, we do note that the available sample sizes of CCSNe and host LIRGs are so far small, and the statistical comparisons between subgroups do not take into account possible systematic or model errors related to the estimated starburst ages. (abridged) △ Less

Submitted 26 February, 2021; originally announced February 2021.

Comments: 25 pages, 21 figures, accepted for publication in A&A

Solving Inverse Problems for Spectral Energy Distributions with Deep Generative Networks

Authors: Agapi Rissaki, Orestis Pavlou, Dimitris Fotakis, Vicky Papadopoulou, Andreas Efstathiou

Abstract: We propose an end-to-end approach for solving inverse problems for a class of complex astronomical signals, namely Spectral Energy Distributions (SEDs). Our goal is to reconstruct such signals from scarce and/or unreliable measurements. We achieve that by leveraging a learned structural prior in the form of a Deep Generative Network. Similar methods have been tested almost exclusively for images w… ▽ More We propose an end-to-end approach for solving inverse problems for a class of complex astronomical signals, namely Spectral Energy Distributions (SEDs). Our goal is to reconstruct such signals from scarce and/or unreliable measurements. We achieve that by leveraging a learned structural prior in the form of a Deep Generative Network. Similar methods have been tested almost exclusively for images which display useful properties (e.g., locality, periodicity) that are implicitly exploited. However, SEDs lack such properties which make the problem more challenging. We manage to successfully extend the methods to SEDs using a Generative Latent Optimization model trained with significantly fewer and corrupted data. △ Less

Submitted 9 December, 2020; originally announced December 2020.

Comments: Accepted to NeurIPS 2020 Workshop on Machine Learning and the Physical Sciences

ACM Class: J.2; I.5.4; I.2.6

Star formation and nuclear activity in luminous infrared galaxies: An infrared through radio review

Authors: Miguel Perez-Torres, Seppo Mattila, Almudena Alonso-Herrero, Susanne Aalto, Andreas Efstathiou

Abstract: Nearby galaxies offer unique laboratories allowing multi-wavelength spatially resolved studies of the interstellar medium, star formation and nuclear activity across a broad range of physical conditions. In particular, detailed studies of individual local luminous infrared galaxies (LIRGs) are crucial for gaining a better understanding of these processes and for developing and testing models that… ▽ More Nearby galaxies offer unique laboratories allowing multi-wavelength spatially resolved studies of the interstellar medium, star formation and nuclear activity across a broad range of physical conditions. In particular, detailed studies of individual local luminous infrared galaxies (LIRGs) are crucial for gaining a better understanding of these processes and for developing and testing models that are used to explain statistical studies of large populations of such galaxies at high redshift for which it is currently impossible to reach a sufficient physical resolution. Here, we provide an overview of the impact of spatially resolved infrared, sub-millimetre and radio observations in the study of the interstellar medium, star formation and active galactic nuclei as well as their interplay in local LIRGs. We also present an overview of the modelling of their spectral energy distributions using state-of-the-art radiative transfer codes. These contribute necessary and powerful 'workhorse' tools for the study of LIRGs (and their more luminous counterparts) at higher redshifts which are unresolved in observations. We describe how spatially-resolved time domain observations have recently opened a new window to study the nuclear activity in LIRGs. We describe in detail the observational characteristics of Arp 299 which is one of the best studied local LIRGs and exemplifies the power of the combination of high-resolution observations at infrared to radio wavelengths together with radiative transfer modelling used to explain the spectral energy distributions of its different components. We summarise the previous achievements obtained using high-spatial resolution observations and provide an outlook into what we can expect to achieve with future facilities. △ Less

Submitted 10 October, 2020; originally announced October 2020.

Comments: 104 pages, 25 figures, Invited review article for The Astronomy and Astrophysics Review, accepted for publication

AT 2017gbl: a dust obscured TDE candidate in a luminous infrared galaxy

Authors: E. C. Kool, T. M. Reynolds, S. Mattila, E. Kankare, M. A. Perez-Torres, A. Efstathiou, S. Ryder, C. Romero-Canizales, W. Lu, T. Heikkila, G. E. Anderson, M. Berton, J. Bright, G. Cannizzaro, D. Eappachen, M. Fraser, M. Gromadzki, P. G. Jonker, H. Kuncarayakti, P. Lundqvist, K. Maeda, R. M. McDermid, A. M. Medling, S. Moran, A. Reguitti , et al. (4 additional authors not shown)

Abstract: We present the discovery with Keck of the extremely infrared (IR) luminous transient AT 2017gbl, coincident with the Northern nucleus of the luminous infrared galaxy (LIRG) IRAS 23436+5257. Our extensive multi-wavelength follow-up spans ~900 days, including photometry and spectroscopy in the optical and IR, and (very long baseline interferometry) radio and X-ray observations. Radiative transfer mo… ▽ More We present the discovery with Keck of the extremely infrared (IR) luminous transient AT 2017gbl, coincident with the Northern nucleus of the luminous infrared galaxy (LIRG) IRAS 23436+5257. Our extensive multi-wavelength follow-up spans ~900 days, including photometry and spectroscopy in the optical and IR, and (very long baseline interferometry) radio and X-ray observations. Radiative transfer modelling of the host galaxy spectral energy distribution and long-term pre-outburst variability in the mid-IR indicate the presence of a hitherto undetected dust obscured active galactic nucleus (AGN). The optical and near-IR spectra show broad 2000 km/s hydrogen, He I and O I emission features that decrease in flux over time. Radio imaging shows a fast evolving compact source of synchrotron emission spatially coincident with AT 2017gbl. We infer a lower limit for the radiated energy of 7.3 x 10^50 erg from the IR photometry. An extremely energetic supernova would satisfy this budget, but is ruled out by the radio counterpart evolution. Instead, we propose AT 2017gbl is related to an accretion event by the central supermassive black hole, where the spectral signatures originate in the AGN broad line region and the IR photometry is consistent with re-radiation by polar dust. Given the fast evolution of AT 2017gbl, we deem a tidal disruption event (TDE) of a star a more plausible scenario than a dramatic change in the AGN accretion rate. This makes AT 2017gbl the third TDE candidate to be hosted by a LIRG, in contrast to the so far considered TDE population discovered at optical wavelengths and hosted preferably by post-starburst galaxies. △ Less

Submitted 4 August, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: Accepted for publication by MNRAS. 32 pages, 22 figures

High-resolution VLA Imaging of Obscured Quasars: Young Radio Jets Caught in a Dense ISM

Authors: Pallavi Patil, Kristina Nyland, Mark Whittle, Carol Lonsdale, Mark Lacy, Colin Lonsdale, Dipanjan Mukherjee, A. C. Trapp, Amy E. Kimball, Lauranne Lanz, Belinda J. Wilkes, Andrew Blain, Jeremy J. Harwood, Andreas Efstathiou, Catherine Vlahakis

Abstract: We present new sub-arcsecond-resolution Karl G. Jansky Very Large Array (VLA) imaging at 10 GHz of 155 ultra-luminous ($L_{\rm bol}\sim10^{11.7-14.2} L_\odot$) and heavily obscured quasars with redshifts $z \sim0.4-3$. The sample was selected to have extremely red mid-infrared (MIR)-optical color ratios based on data from Wide-Field Infrared Survey Explorer (WISE) along with a detection of bright,… ▽ More We present new sub-arcsecond-resolution Karl G. Jansky Very Large Array (VLA) imaging at 10 GHz of 155 ultra-luminous ($L_{\rm bol}\sim10^{11.7-14.2} L_\odot$) and heavily obscured quasars with redshifts $z \sim0.4-3$. The sample was selected to have extremely red mid-infrared (MIR)-optical color ratios based on data from Wide-Field Infrared Survey Explorer (WISE) along with a detection of bright, unresolved radio emission from the NRAO VLA Sky Survey (NVSS) or Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) Survey. Our high-resolution VLA observations have revealed that the majority of the sources in our sample (93 out of 155) are compact on angular scales $<0.2^{\prime \prime}$ ($\leq 1.7$ kpc at $z \sim2$). The radio luminosities, linear extents, and lobe pressures of our sources are similar to young radio active galactic nuclei (AGN; e.g., Gigahertz Peaked Spectrum, GPS, and Compact Steep Spectrum, CSS, sources), but their space density is considerably lower. Application of a simple adiabatic lobe expansion model suggests relatively young dynamical ages ($\sim10^{4-7}$ years), relatively high ambient ISM densities ($\sim1-10^4$ cm$^{-3}$), and modest lobe expansion speeds ($\sim30-10,000$ km s$^{-1}$). Thus, we find our sources to be consistent with a population of newly triggered, young jets caught in a unique evolutionary stage in which they still reside within the dense gas reservoirs of their hosts. Based on their radio luminosity function and dynamical ages, we estimate only $\sim20\%$ of classical large scale FRI/II radio galaxies could have evolved directly from these objects. We speculate that the WISE-NVSS sources might first become GPS or CSS sources, of which some might ultimately evolve into larger radio galaxies. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: Accepted for publication in ApJ, 38 pages, 14 figures, 6 tables

Have we seen all the galaxies that comprise the cosmic infrared background at 250\,$μ$m $\le λ\le$ 500\,$μ$m?

Authors: S. Duivenvoorden, S. Oliver, M. Bethermin, D. L. Clements, G. De Zotti, A. Efstathiou, D. Farrah, P. D. Hurley, R. J. Ivison, G. Lagache, D. Scott, R. Shirley, L. Wang, M. Zemcov

Abstract: The cosmic infrared background (CIB) provides a fundamental observational constraint on the star-formation history of galaxies over cosmic history. We estimate the contribution to the CIB from catalogued galaxies in the COSMOS field by using a novel map fitting technique on the \textit{Herschel} SPIRE maps. Prior galaxy positions are obtained using detections over a large range in wavelengths in t… ▽ More The cosmic infrared background (CIB) provides a fundamental observational constraint on the star-formation history of galaxies over cosmic history. We estimate the contribution to the CIB from catalogued galaxies in the COSMOS field by using a novel map fitting technique on the \textit{Herschel} SPIRE maps. Prior galaxy positions are obtained using detections over a large range in wavelengths in the $K_{\rm s}$--3\,GHz range. Our method simultaneously fits the galaxies, the system foreground, and the leakage of flux from galaxies located in masked areas and corrects for an "over-fitting" effect not previously accounted for in stacking methods. We explore the contribution to the CIB as a function of galaxy survey wavelength and depth. We find high contributions to the CIB with the deep $r$ ($m_{\rm AB} \le 26.5$), $K_{\rm s}$ ($m_{\rm AB} \le 24.0$) and 3.6\,$μ$m ($m_{\rm AB} \le 25.5$) catalogues. We combine these three deep catalogues and find a total CIB contributions of 10.5 $\pm$ 1.6, 6.7 $\pm$ 1.5 and 3.1 $\pm$ 0.7\,nWm$^{-2}$sr$^{-1}$ at 250, 350 and 500\,$μ$m, respectively. Our CIB estimates are consistent with recent phenomenological models, prior based SPIRE number counts and with (though more precise than) the diffuse total measured by FIRAS. Our results raise the interesting prospect that the CIB contribution at $λ\le 500\,μ$m from known galaxies has converged. Future large-area surveys like those with the Large Synoptic Survey Telescope are therefore likely to resolve a substantial fraction of the population responsible for the CIB at 250\,$μ$m $\leq λ\leq$ 500\,$μ$m. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: 16 pages, 14 figures

HELP: A catalogue of 170 million objects, selected at 0.36-4.5 $μ$m, from 1270 deg.$^{2}$ of prime extragalactic fields

Authors: Raphael Shirley, Yannick Roehlly, Peter D Hurley, Veronique Buat, María del Carmen Campos Varillas, Steven Duivenvoorden, Kenneth J Duncan, Andreas Efstathiou, Duncan Farrah, Eduardo González Solares, Katarzyna Małek, Lucia Marchetti, Ian McCheyne, Andreas Papadopoulos, Estelle Pons, Roberto Scipioni, Mattia Vaccari, Seb Oliver

Abstract: We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170… ▽ More We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170 million objects is constructed by a positional cross match of 51 public surveys. This high resolution optical, near-infrared, and mid-infrared catalogue is designed for photometric redshift estimation, extraction of fluxes in lower resolution far-infrared maps, and spectral energy distribution modelling. It collates, standardises, and provides value added derived quantities including corrected aperture magnitudes and astrometry correction over the Herschel extragalactic wide fields for the first time. $grizy$ fluxes are available on all fields with $g$ band data reaching $5σ$ point-source depths in a 2 arcsec aperture of 23.5, 24.4, and 24.6 (AB) mag at the 25th, 50th, and 75th percentiles, by area covered, across all HELP fields. It has $K$ or $K_s$ coverage over 1146 deg.$^2$ with depth percentiles of 20.2, 20.4, and 21.0 mag respectively. The IRAC Ch 1 band is available over 273 deg.$^2$ with depth percentiles of 17.7, 21.4, and 22.2 mag respectively. This paper defines the "masterlist" objects for the first data release (DR1) of HELP. This large sample of standardised total and corrected aperture fluxes, uniform quality flags, and completeness measures provides large well understood statistical samples over the full Herschel extragalactic sky. △ Less

Submitted 22 October, 2019; v1 submitted 9 September, 2019; originally announced September 2019.

Journal ref: MNRAS 490.1 (2019): 634-656

The mid-infrared and CO gas properties of an extreme star-forming FeLoBAL quasar

Authors: Lura K. Pitchford, Duncan Farrah, Katherine Alatalo, José Afonso, Andreas Efstathiou, Evanthia Hatziminaoglou, Mark Lacy, Tanya Urrutia, Giulio Violino

Abstract: We present a detailed study of a high-redshift iron low-ionization broad absorption line (FeLoBAL) quasar (SDSS1214 at $z = 1.046$), including new interferometric $^{12}$CO $J$=2-1 observations, optical through far-infrared photometry, and mid-infrared spectroscopy. The CO line is well-fit by a single Gaussian centered 40 kms$^{-1}$ away from the systemic velocity and implies a total molecular gas… ▽ More We present a detailed study of a high-redshift iron low-ionization broad absorption line (FeLoBAL) quasar (SDSS1214 at $z = 1.046$), including new interferometric $^{12}$CO $J$=2-1 observations, optical through far-infrared photometry, and mid-infrared spectroscopy. The CO line is well-fit by a single Gaussian centered 40 kms$^{-1}$ away from the systemic velocity and implies a total molecular gas mass of $M_\textrm{gas} = 7.3 \times 10^{10} \textrm{M}_\odot$. The infrared SED requires three components: an active galactic nucleus (AGN) torus, an AGN polar dust component, and a starburst. The starburst dominates the infrared emission with a luminosity of log($L_\textrm{SB}[\textrm{L}_\odot]) = 12.91^{+0.02}_{-0.02}$, implying a star formation rate of about 2000 $\textrm{M}_{\odot}$yr$^{-1}$, the highest known among FeLoBAL quasars. The AGN torus and polar dust components are less luminous, at log($L_\textrm{AGN}[\textrm{L}_\odot]) = 12.36^{+0.14}_{-0.15}$ and log($L_\textrm{dust}[\textrm{L}_\odot]) = 11.75^{+0.26}_{-0.46}$, respectively. If all of the molecular gas is used to fuel the ongoing star formation, then the lower limit on the subsequent duration of the starburst is 40 Myr. We do not find conclusive evidence that the AGN is affecting the CO gas reservoir. The properties of SDSS1214 are consistent with it representing the endpoint of an obscured starburst transitioning through a LoBAL phase to that of a classical quasar. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: 10 pages, 7 figures, 3 tables, accepted for publication in MNRAS

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

Authors: Ph. André, A. Hughes, V. Guillet, F. Boulanger, A. Bracco, E. Ntormousi, D. Arzoumanian, A. J. Maury, J. -Ph. Bernard, S. Bontemps, I. Ristorcelli, J. M. Girart, F. Motte, K. Tassis, E. Pantin, T. Montmerle, D. Johnstone, S. Gabici, A. Efstathiou, Shantanu Basu, M. Béthermin, H. Beuther, J. Braine, J. Di Francesco, E. Falgarone , et al. (31 additional authors not shown)

Abstract: SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This… ▽ More SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: Accepted for publication in PASA; 37 pages, 16 figures

Journal ref: PASA, 36, e029 (2019)

HELP: modelling the spectral energy distributions of Herschel detected galaxies in the ELAIS N1 field

Authors: K. Malek, V. Buat, Y. Roehlly, D. Burgarella, P. D. Hurley, R. Shirley, K. Duncan, A. Efstathiou, A. Papadopoulos, M. Vaccari, D. Farrah, L. Marchetti, S. Oliver

Abstract: The HELP project focuses on the data from ESA's Herschel mission, which covered over 1300$deg^2$ and is preparing to publish a multi-wavelength catalogue of millions of objects. Our main goal is to find the best approach to simultaneously fitting SEDs of millions of galaxies across a wide redshift range to obtain homogeneous estimates of the main physical parameters of detected IR galaxies. We per… ▽ More The HELP project focuses on the data from ESA's Herschel mission, which covered over 1300$deg^2$ and is preparing to publish a multi-wavelength catalogue of millions of objects. Our main goal is to find the best approach to simultaneously fitting SEDs of millions of galaxies across a wide redshift range to obtain homogeneous estimates of the main physical parameters of detected IR galaxies. We perform SED fitting on the UV/NIR to FIR emission of 42 047 galaxies from the pilot HELP field: ELAIS N1. We use the latest release of CIGALE, a galaxy SED fitting code relying on energy balance, to deliver the main physical parameters such as $M_{star}$, SFR, and $L_{dust}$. We implement additional quality criteria to the fits by calculating $χ^2$ values for the stellar and dust part of the spectra independently. These criteria allow us to identify the best fits and to identify peculiar galaxies. We perform the SED fitting by assuming three different dust attenuation laws separately allowing us to test the impact of the assumed law on estimated physical parameters. We implemented two additional quality value checks for the SED fitting method based on $M_{star}$ estimation and energy budget. This method allows us to identify possible objects with incorrect matching in the catalogue and peculiar galaxies; we found 351 possible candidates of lensed galaxies using two complementary $χ^2$s criteria (stellar and IR) and $z_{phot}$ calculated for the IR part of the spectrum only. We find that the attenuation law has an important impact on the $M_{star}$ estimate (on average leading to disparities of a factor of two). We derive the relation between $M_{star}$ estimates obtained by different attenuation laws and we find the best recipe for our sample. We also make independent estimates of the total $L_{dust}$ parameter from stellar emission by fitting the galaxies with and without IR data separately. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: 17 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 620, A50 (2018)

Dust attenuation and Halpha emission in a sample of galaxies observed with Herschel at $0.6 < z <1.6$

Authors: V. Buat, M. Boquien, K. Malek, D. Corre, H. salas, Y. Roehlly, R. Shirley, A. Efstathiou

Abstract: Dust attenuation shapes the spectral energy distribution of galaxies. It is particularly true for dusty galaxies in which stars experience a heavy attenuation. The combination of UV-to-IR photometry with the spectroscopic measurement of the H$α$ recombination line helps to quantify dust attenuation of the whole stellar population and its wavelength dependence. We selected an IR complete sample of… ▽ More Dust attenuation shapes the spectral energy distribution of galaxies. It is particularly true for dusty galaxies in which stars experience a heavy attenuation. The combination of UV-to-IR photometry with the spectroscopic measurement of the H$α$ recombination line helps to quantify dust attenuation of the whole stellar population and its wavelength dependence. We selected an IR complete sample of galaxies in the COSMOS 3D-HST CANDELS field detected with the Herschel satellite with a signal to noise ratio larger than five. Optical to NIR photometry is available as well as NIR spectroscopy for each source. We reduced the sample to the redshift range $0.6 < z < 1.6$ to include the H$α$ line in the G141 grism spectra. We have used a new version of the CIGALE code to fit simultaneously the continuum and H$α$ line emission of the 34 selected galaxies. Using flexible attenuation laws with free parameters, we are able to measure the shape of the attenuation curve for each galaxy as well as the amount of attenuation of each stellar population, the former being in general steeper than the starburst law in the UV-optical with a large variation of the slope among galaxies. The attenuation of young stars or nebular continuum is found on average about twice the attenuation affecting older stars, again with a large variation. Our model with power-laws, based on a modification of the Charlot and Fall recipe, gives results in better agreement with the radiative transfer models than the global modification of the slope of the Calzetti law. △ Less

Submitted 1 September, 2018; originally announced September 2018.

Comments: 12 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 619, A135 (2018)

AKARI and IRAS: From Beam Corrections to SEDs

Authors: David L. Clements, Michael Rowan-Robinson, Chris Pearson, Jose Afonso, Vianney Labouteiller, Duncan Farrah, Andreas Efstathiou, Josh Greenslade, Lingyu Wang

Abstract: There is significant scientific value to be gained from combining AKARI fluxes with data at other far-IR wavelengths from the IRAS and Herschel missions. To be able to do this we must ensure that there are no systematic differences between the datasets that need to be corrected before the fluxes are compatible with each other. One such systematic effect identified in the BSCv1 data is the issue of… ▽ More There is significant scientific value to be gained from combining AKARI fluxes with data at other far-IR wavelengths from the IRAS and Herschel missions. To be able to do this we must ensure that there are no systematic differences between the datasets that need to be corrected before the fluxes are compatible with each other. One such systematic effect identified in the BSCv1 data is the issue of beam corrections. We determine these for the BSCv2 data by correlating ratios of appropriate IRAS and AKARI bands with the difference in 2MASS J band extended and point source magnitudes for sources cross matched between the IRAS FSC, AKARI BSCv2 and 2MASS catalogs. We find significant correlations (p<< 10^-13) indicating that beam corrections are necessary in the 65 and 90 micron bands. We then use these corrected fluxes to supplement existing data in spectral energy distribution (SED) fits for ultraluminous infrared galaxies in the HERUS survey. The addition of AKARI fluxes makes little difference to the results of simple (T,beta) fits to the SEDs of these sources, though there is a general decrease in reduced chi^2 values. The utility of the extra AKARI data, however, is in allowing physically more realistic SED models with more parameters to be fit to the data. We also extend our analysis of beam correction issues in the AKARI data by examining the Herschel Reference Sample, which have Herschel photometry from 100 to 500 microns and which are more spatially extended than the HERUS ULIRGs. Thirty-four of the HRS sources have good Herschel SEDs and matching data from AKARI. This investigation finds that our simple 2MASS-based beam correction scheme is inadequate for these larger and more complex sources. There are also indications that additional beam corrections at 140 and 160 microns are needed for these sources, extended on scales >1 arcminute. △ Less

Submitted 10 August, 2018; originally announced August 2018.

Comments: Accepted for publication in the PASJ special issue on AKARI. 3 Figures

A dust-enshrouded tidal disruption event with a resolved radio jet in a galaxy merger

Authors: S. Mattila, M. Pérez-Torres, A. Efstathiou, P. Mimica, M. Fraser, E. Kankare, A. Alberdi, M. Á. Aloy, T. Heikkilä, P. G. Jonker, P. Lundqvist, I. Martí-Vidal, W. P. S. Meikle, C. Romero-Cañizales, S. J. Smartt, S. Tsygankov, E. Varenius, A. Alonso-Herrero, M. Bondi, C. Fransson, R. Herrero-Illana, T. Kangas, R. Kotak, N. Ramírez-Olivencia, P. Väisänen , et al. (11 additional authors not shown)

Abstract: Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its… ▽ More Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: Published online by Science on 14 Jun 2018

The emission and distribution of dust of the torus of NGC 1068

Authors: Enrique Lopez-Rodriguez, Lindsay Fuller, Almudena Alonso-Herrero, Andreas Efstathiou, Kohei Ichikawa, Nancy A. Levenson, Chris Packham, James Radomski, Cristina Ramos Almeida, Dominic J. Benford, Marc Berthoud, Ryan Hamilton, Doyal Harper, Attila Kovavcs, Fabio P. Santos, J. Staguhn, Terry Herter

Abstract: We present observations of NGC 1068 covering the $19.7-53.0$ $μ$m wavelength range using FORCAST and HAWC+ onboard SOFIA. Using these observations, high-angular resolution infrared (IR) and sub-mm observations, we find an observational turn-over of the torus emission in the $30-40$ $μ$m wavelength range with a characteristic temperature of $70-100$ K. This component is clearly different from the d… ▽ More We present observations of NGC 1068 covering the $19.7-53.0$ $μ$m wavelength range using FORCAST and HAWC+ onboard SOFIA. Using these observations, high-angular resolution infrared (IR) and sub-mm observations, we find an observational turn-over of the torus emission in the $30-40$ $μ$m wavelength range with a characteristic temperature of $70-100$ K. This component is clearly different from the diffuse extended emission in the narrow line and star formation regions at 10-100 $μ$m within the central 700 pc. We compute $2.2-432$ $μ$m 2D images using the best inferred \textsc{clumpy} torus model based on several nuclear spectral energy distribution (SED) coverages. We find that when $1-20$ $μ$m SED is used, the inferred result gives a small torus size ($<4$ pc radius) and a steep radial dust distribution. The computed torus using the $1-432$ $μ$m SED provides comparable torus sizes, $5.1^{+0.4}_{-0.4}$ pc radius, and morphology to the recently resolved 432 $μ$m ALMA observations. This result indicates that the $1-20$ $μ$m wavelength range is not able to probe the full extent of the torus. The characterization of the turn-over emission of the torus using the $30-60$ $μ$m wavelength range is sensitive to the detection of cold dust in the torus. The morphology of the dust emission in our 2D image at 432 $μ$m is spatially coincident with the cloud distribution, while the morphology of the emission in the $1-20$ $μ$m wavelength range shows an elongated morphology perpendicular to the cloud distribution. We find that our 2D \textsc{clumpy} torus image at 12 $μ$m can produce comparable results to those observed using IR interferometry. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Comments: 12 pages, 7 Figures, Accepted for publication in ApJ

Red, redder, reddest: SCUBA-2 imaging of colour-selected \textit{Herschel} sources

Authors: S. Duivenvoorden, S. Oliver, J. M. Scudder, J. Greenslade, D. A. Riechers, S. M. Wilkins, V. Buat, S. C. Chapman, D. L. Clements, A. Cooray, K. E. K. Coppin, H. Dannerbauer, G. De Zotti, J. S. Dunlop, S. A. Eales, A. Efstathiou, D. Farrah, J. E. Geach, W. S. Holland, P. D. Hurley, R. J. Ivison, L. Marchetti, G. Petitpas, M. T. Sargent, D. Scott , et al. (6 additional authors not shown)

Abstract: High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $μ$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We d… ▽ More High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $μ$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$μ$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$μ$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $σ_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$). △ Less

Submitted 22 January, 2018; originally announced January 2018.

Comments: 22 pages, 14 figures and 3 tables

HERUS: The Far-IR/Submm Spectral Energy Distributions of Local ULIRGs & Photometric Atlas

Authors: D. L. Clements, C. Pearson, D. Farrah, J. Greenslade, Jeronimo Bernard-Salas, E. Gonzalez-Alfonso, J. Afonso, A. Efstathiou, D. Rigopoulou, V. Lebouteiller, P. D. Hurley, H. Spoon

Abstract: We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60$μ$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500$μ$m. We describe these observations, present the results, and combine the new observations with data fro… ▽ More We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60$μ$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500$μ$m. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-IR spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parameterised modified black body model where temperature and emissivity $β$ are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9+/-4.7 K compared to 21.3+/-3.4 K) but a similar $β$ distribution (median $β$ = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al., 2014) galaxies. Dust masses are found to be in the range of 10^7.5 to 10^9 Msun significantly higher than that of Herschel Reference Sample (HRS) sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850$μ$m. These latter sources generally have cooler dust and/or redder 100-to-250 $μ$m colours than our 60$μ$m-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: Accepted for publication in MNRAS

Unbiased large spectroscopic surveys of galaxies selected by SPICA using dust bands

Authors: H. Kaneda, D. Ishihara, S. Oyabu, M. Yamagishi, T. Wada, L. Armus, M. Baes, V. Charmandaris, B. Czerny, A. Efstathiou, J. A. Fern'andez-Ontiveros, A. Ferrara, E. Gonz'alez-Alfonso, M. Griffin, C. Gruppioni, E. Hatziminaoglou, M. Imanishi, K. Kohno, J. Kwon, T. Nakagawa, T. Onaka, F. Pozzi, D. Scott, J. -D. T. Smith, L. Spinoglio , et al. (5 additional authors not shown)

Abstract: The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic… ▽ More The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic surveys of galaxies that could be performed by SPICA. The surveys are designed to address the question how star formation and black hole accretion activities evolved over cosmic time through spectral diagnostics of the physical conditions of the interstellar/circumnuclear media in galaxies. On the basis of results obtained with Herschel far-IR photometric surveys of distant galaxies and Spitzer and AKARI near- to mid-IR spectroscopic observations of nearby galaxies, we estimate the numbers of the galaxies at redshift z > 0.5, which are expected to be detected in the PAH features or dust continuum by a wide (10 deg^2) or deep (1 deg^2) blind survey, both for a given observation time of 600 hours. As by-products of the wide blind survey, we also expect to detect debris disks, through the mid-IR excess above the photospheric emission of nearby main-sequence stars, and we estimate their number. We demonstrate that the SPICA mid-IR surveys will efficiently provide us with unprecedentedly large spectral samples, which can be studied further in the far-IR with SPICA. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: Paper accepted for publication on PASA on 18th October 2017, as part of the SPICA Special Issue

Feedback and feeding in the context of galaxy evolution with SPICA: direct characterization of molecular outflows and inflows

Authors: E. González-Alfonso, L. Armus, F. J. Carrera, V. Charmandaris, A. Efstathiou, E. Egami, J. A. Fernández-Ontiveros, J. Fischer, G. L. Granato, C. Gruppioni, E. Hatziminaoglou, M. Imanishi, N. Isobe, H. Kaneda, D. Koziel-Wierzbowska, M. A. Malkan, J. Martin-Pintado, S. Mateos, H. Matsuhara, G. Miniutti, T. Nakagawa, F. Pozzi, F. Rico-Villas, G. Rodighiero, P. Roelfsema , et al. (6 additional authors not shown)

Abstract: A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas.… ▽ More A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: This paper belongs to the SPICA Special Issue on PASA. Accepted for publication in PASA

Tracing the evolution of dust obscured star-formation and accretion back to the reionisation epoch with SPICA

Authors: C. Gruppioni, L. Ciesla, E. Hatziminaoglou, F. Pozzi, G. Rodighiero, P. Santini, L. Armus, M. Baes, J. Braine, V. Charmandaris, D. L. Clements, N. Christopher, H. Dannerbauer, A. Efstathiou, E. Egami, J. A. Fernandez-Ontiveros, F. Fontanot, A. Franceschini, E. Gonzalez-Alfonso, M. Griffin, H. Kaneda, L. Marchetti, P. Monaco, T. Nakagawa, T. Onaka , et al. (13 additional authors not shown)

Abstract: Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the earl… ▽ More Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: This paper belongs to the SPICA Special Issue on PASA

Galaxy evolution studies with the SPace IR telescope for Cosmology and Astrophysics (SPICA): the power of IR spectroscopy

Authors: L. Spinoglio, A. Alonso-Herrero, L. Armus, M. Baes, J. Bernard-Salas, S. Bianchi, M. Bocchio, A. Bolatto, C. M. Bradford, J. Braine, F. J. Carrera, L. Ciesla, D. L. Clements, H. Dannerbauer, Y. Doi, A. Efstathiou, E. Egami, J. A. Fernandez-Ontiveros, A. Ferrara, J. Fischer, A. Franceschini, S. Gallerani, M. Giard, E. Gonzalez-Alfonso, C. Gruppioni , et al. (43 additional authors not shown)

Abstract: IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Larg… ▽ More IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: Paper accepted for publication on PASA on 21st September 2017, as part of the SPICA Special Issue

First results from GeMS/GSAOI for project SUNBIRD: Supernovae UNmasked By Infra-Red Detection

Authors: E. C. Kool, S. Ryder, E. Kankare, S. Mattila, T. Reynolds, R. M. McDermid, M. A. Pérez-Torres, R. Herrero-Illana, M. Schirmer, A. Efstathiou, F. E. Bauer, J. Kotilainen, P. Väisänen, C. Baldwin, C. Romero-Cañizales, A. Alberdi

Abstract: Core collapse supernova (CCSN) rates suffer from large uncertainties as many CCSNe exploding in regions of bright background emission and significant dust extinction remain unobserved. Such a shortfall is particularly prominent in luminous infrared galaxies (LIRGs), which have high star formation (and thus CCSN) rates and host bright and crowded nuclear regions, where large extinctions and reduced… ▽ More Core collapse supernova (CCSN) rates suffer from large uncertainties as many CCSNe exploding in regions of bright background emission and significant dust extinction remain unobserved. Such a shortfall is particularly prominent in luminous infrared galaxies (LIRGs), which have high star formation (and thus CCSN) rates and host bright and crowded nuclear regions, where large extinctions and reduced search detection efficiency likely lead to a significant fraction of CCSNe remaining undiscovered. We present the first results of project SUNBIRD (Supernovae UNmasked By Infra-Red Detection), where we aim to uncover CCSNe that otherwise would remain hidden in the complex nuclear regions of LIRGs, and in this way improve the constraints on the fraction that is missed by optical seeing-limited surveys. We observe in the near-infrared 2.15 μm $K_s$-band, which is less affected by dust extinction compared to the optical, using the multi-conjugate adaptive optics imager GeMS/GSAOI on Gemini South, allowing us to achieve a spatial resolution that lets us probe close in to the nuclear regions. During our pilot program and subsequent first full year we have discovered three CCSNe and one candidate with projected nuclear offsets as small as 200 pc. When compared to the total sample of LIRG CCSNe discovered in the near-IR and optical, we show that our method is singularly effective in uncovering CCSNe in nuclear regions and we conclude that the majority of CCSNe exploding in LIRGs are not detected as a result of dust obscuration and poor spatial resolution. △ Less

Submitted 24 September, 2017; originally announced September 2017.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society, 19 pages, 11 figures

Characterizing the UV-to-NIR shape of the dust attenuation curve of IR luminous galaxies up to z$\sim$2

Authors: B. Lo Faro, V. Buat, Y. Roehlly, J. Álvarez-Márquez, D. Burgarella, L. Silva, A. Efstathiou

Abstract: In this work we investigate the far-UV to NIR shape of the dust attenuation curve of a sample of IR selected dust obscured (U)LIRGs at z$\sim$2. The spectral energy distributions (SEDs) are fitted with CIGALE, a physically-motivated spectral synthesis model based on energy balance. Its flexibility allows us to test a wide range of different analytical prescriptions for the dust attenuation curve,… ▽ More In this work we investigate the far-UV to NIR shape of the dust attenuation curve of a sample of IR selected dust obscured (U)LIRGs at z$\sim$2. The spectral energy distributions (SEDs) are fitted with CIGALE, a physically-motivated spectral synthesis model based on energy balance. Its flexibility allows us to test a wide range of different analytical prescriptions for the dust attenuation curve, including the well-known Calzetti and Charlot & Fall curves, and modified versions of them. The attenuation curves computed under the assumption of our reference double power-law model are in very good agreement with those derived, in previous works, with radiative transfer (RT) SED fitting. We investigate the position of our galaxies in the IRX-$β$ diagram and find this to be consistent with grayer slopes, on average, in the UV. We also find evidence for a flattening of the attenuation curve in the NIR with respect to more classical Calzetti-like recipes. This larger NIR attenuation yields larger derived stellar masses from SED fitting, by a median factor of $\sim$ 1.4 and up to a factor $\sim$10 for the most extreme cases. The star formation rate appears instead to be more dependent on the total amount of attenuation in the galaxy. Our analysis highlights the need for a flexible attenuation curve when reproducing the physical properties of a large variety of objects. △ Less

Submitted 31 July, 2017; originally announced July 2017.

Comments: 20 pages, 14 figures

Star formation and AGN activity in a sample of local Luminous Infrared Galaxies through multi-wavelength characterization

Authors: Rubén Herrero-Illana, Miguel Á. Pérez-Torres, Zara Randriamanakoto, Antxon Alberdi, Andreas Efstathiou, Petri Väisänen, Erkki Kankare, Erik Kool, Seppo Mattila, Rajin Ramphul, Stuart Ryder

Abstract: Nuclear starbursts and AGN activity are the main heating processes in luminous infrared galaxies (LIRGs) and their relationship is fundamental to understand galaxy evolution. In this paper, we study the star-formation and AGN activity of a sample of 11 local LIRGs imaged with subarcsecond angular resolution at radio (8.4GHz) and near-infrared ($2.2μ$m) wavelengths. This allows us to characterize t… ▽ More Nuclear starbursts and AGN activity are the main heating processes in luminous infrared galaxies (LIRGs) and their relationship is fundamental to understand galaxy evolution. In this paper, we study the star-formation and AGN activity of a sample of 11 local LIRGs imaged with subarcsecond angular resolution at radio (8.4GHz) and near-infrared ($2.2μ$m) wavelengths. This allows us to characterize the central kpc of these galaxies with a spatial resolution of $\simeq100$pc. In general, we find a good spatial correlation between the radio and the near-IR emission, although radio emission tends to be more concentrated in the nuclear regions. Additionally, we use an MCMC code to model their multi-wavelength spectral energy distribution (SED) using template libraries of starburst, AGN and spheroidal/cirrus models, determining the luminosity contribution of each component, and finding that all sources in our sample are starburst-dominated, except for NGC6926 with an AGN contribution of $\simeq64$\%. Our sources show high star formation rates ($40$ to $167M_\odot\mathrm{yr}^{-1}$), supernova rates (0.4 to $2.0\mathrm{SN}\mathrm{yr}^{-1}$), and similar starburst ages (13 to $29\mathrm{Myr}$), except for the young starburst (9Myr) in NGC6926. A comparison of our derived star-forming parameters with estimates obtained from different IR and radio tracers shows an overall consistency among the different star formation tracers. AGN tracers based on mid-IR, high-ionization line ratios also show an overall agreement with our SED model fit estimates for the AGN. Finally, we use our wide-band VLA observations to determine pixel-by-pixel radio spectral indices for all galaxies in our sample, finding a typical median value ($α\simeq-0.8$) for synchrotron-powered LIRGs. △ Less

Submitted 3 July, 2017; v1 submitted 26 May, 2017; originally announced May 2017.

Comments: Accepted for publication in MNRAS. 20 pages, 12 figures

The role of the most luminous, obscured AGN in galaxy assembly at z~2

Authors: Duncan Farrah, Sara Petty, Brian Connolly, Andrew Blain, Andreas Efstathiou, Mark Lacy, Daniel Stern, Sean Lake, Tom Jarrett, Carrie Bridge, Peter Eisenhardt, Dominic Benford, Suzy Jones, Chao-Wei Tsai, Roberto Assef, Jingwen Wu, Leonidas Moustakas

Abstract: We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least… ▽ More We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least in most part from AGN activity. The AGN fractional luminosities are higher than those in either sub-millimeter galaxies, or AGN selected via other mid-infrared criteria. Adopting the $G$, M$_{20}$ and $A$ morphological parameters, together with traditional classification boundaries, infers that three quarters of the sample as mergers. Our sample do not, however, show any correlation between the considered morphological parameters and either infrared luminosity or AGN fractional luminosity. Moreover, their asymmetries and effective radii are distributed identically to those of massive galaxies at $z\sim2$. We conclude that our sample is not preferentially associated with mergers, though a significant merger fraction is still plausible. Instead, we propose that our sample are examples of the massive galaxy population at $z\sim2$ that harbor a briefly luminous, "flickering" AGN, and in which the $G$ and M$_{20}$ values have been perturbed, due to either the AGN, and/or the earliest formation stages of a bulge in an inside-out manner. Furthermore, we find that the mass assembly of the central black holes in our sample leads the mass assembly of any bulge component. Finally, we speculate that our sample represent a small fraction of the immediate antecedents of compact star-forming galaxies at $z\sim2$. △ Less

Submitted 6 June, 2017; v1 submitted 7 May, 2017; originally announced May 2017.

Comments: ApJ, accepted. Updated to reflect the accepted version

Unveiling the AGN in IC 883: discovery of a parsec-scale radio jet

Authors: Cristina Romero-Cañizales, Antxon Alberdi, Claudio Ricci, Patricia Arévalo, Miguel A. Pérez-Torres, John E. Conway, Rob J. Beswick, Marco Bondi, Tom W. B. Muxlow, Megan K. Argo, Franz E. Bauer, Andreas Efstathiou, Rubén Herrero-Illana, Seppo Mattila, Stuart D. Ryder

Abstract: IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transi… ▽ More IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transient radio source with luminosity typical of bright supernovae, give further evidence of the ongoing star formation activity, which dominates the energetics of the system. At sub-parsec scales, the radio nucleus has a core-jet morphology with the jet being a newly ejected component showing a subluminal proper motion of 0.6c-1c. The AGN contributes less than two per cent of the total IR luminosity of the system. The corresponding Eddington factor is ~1E-3, suggesting this is a low-accretion rate engine, as often found in LLAGNs. However, its high bolometric luminosity (~10E44erg/s) agrees better with a normal AGN. This apparent discrepancy may just be an indication of the transition nature of the nucleus from a system dominated by star-formation, to an AGN-dominated system. The nucleus has a strongly inverted spectrum and a turnover at ~4.4GHz, thus qualifying as a candidate for the least luminous (L_5.0GHz ~ 6.3E28erg/s/Hz) and one of the youngest (~3000yr) gigahertz-peaked spectrum (GPS) sources. If the GPS origin for the IC883 nucleus is confirmed, then advanced mergers in the LIRG category are potentially key environments to unveil the evolution of GPS sources into more powerful radio galaxies. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: 13 pages, 8 figures. Accepted for publication in MNRAS

HERUS: A CO Atlas from SPIRE Spectroscopy of local ULIRGs

Authors: Chris Pearson, Dimitra Rigopoulou, Peter Hurley, Duncan Farrah, Jose Afonso, Jeronimo Bernard-Salas, Colin Borys, David L. Clements, Diane Cormier, Andreas Efstathiou, Eduardo Gonzalez-Alfonso, Vianney Lebouteiller, Henrik Spoon

Abstract: We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectr… ▽ More We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671μ$m with photometry in the SPIRE bands at 250$μ$m, 350$μ$m and 500$μ$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $μ$m, [CI] 609 $μ$m, and [NII] 205 $μ$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.} △ Less

Submitted 19 October, 2016; originally announced October 2016.

HELP: star formation as function of galaxy environmentwith Herschel

Authors: S. Duivenvoorden, S. Oliver, V. Buat, B. Darvish, A. Efstathiou, D. Farrah, M. Griffin, P. D. Hurley, E. Ibar, M. Jarvis, A. Papadopoulos, M. T. Sargent, D. Scott, J. M. Scudder, M. Symeonidis, M. Vaccari, M. P. Viero, L. Wang

Abstract: The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust obscured star formation over 1300 deg$^2$. With this unprecedented combination of data sets, it is possible to investigate how the star formation vs stellar mass relation (main-sequence) of star-forming galaxies depends… ▽ More The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust obscured star formation over 1300 deg$^2$. With this unprecedented combination of data sets, it is possible to investigate how the star formation vs stellar mass relation (main-sequence) of star-forming galaxies depends on environment. In this pilot study we explore this question between 0.1 < z < 3.2 using data in the COSMOS field. We estimate the local environment from a smoothed galaxy density field using the full photometric redshift probability distribution. We estimate star formation rates by stacking the SPIRE data from the Herschel Multi-tiered Extragalactic Survey (HerMES). Our analysis rules out the hypothesis that the main-sequence for star-forming systems is independent of environment at 1.5 < z < 2, while a simple model in which the mean specific star formation rate declines with increasing environmental density gives a better description. However, we cannot exclude a simple hypothesis in which the main-sequence for star-forming systems is independent of environment at z < 1.5 and z > 2. We also estimate the evolution of the star formation rate density in the COSMOS field and our results are consistent with previous measurements at z < 1.5 and z > 2 but we find a $1.4^{+0.3}_{-0.2}$ times higher peak value of the star formation rate density at $z \sim 1.9$. △ Less

Submitted 16 June, 2016; originally announced June 2016.

Comments: 14 pages, 7 figures. MNRAS accepted

Star Formation in Luminous Quasars at 2<z<3

Authors: Kathryn Harris, Duncan Farrah, Bernhard Schulz, Evanthia Hatziminaoglou, Marco Viero, Nick Anderson, Matthieu Bethermin, Scott Chapman, David L. Clements, Asantha Cooray, Andreas Efstathiou, Anne Feltre, Peter Hurley, Eduardo Ibar, Mark Lacy, Sebastian Oliver, Mathew J. Page, Ismael Perez-Fournon, Sara M. Petty, Lura K. Pitchford, Dimitra Rigopoulou, Douglas Scott, Myrto Symeonidis, Joaquin Vieira, Lingyu Wang

Abstract: We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr… ▽ More We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: MNRAS, accepted

Radio Jet Feedback and Star Formation in Heavily Obscured Quasars at Redshifts ~0.3-3, I: ALMA Observations

Authors: Carol J. Lonsdale, Mark Lacy, Amy E. Kimball, Andrew Blain, Mark Whittle, Belinda Wilkes, Dan Stern, Jim Condon, Minjin Kim, Roberto J. Assef, Chao-Wei Tsai, Andreas Efstathiou, Suzy Jones, Peter Eisenhardt, Carrie Bridge, Jinwen Wu, Colin J. Lonsdale, Kristen Jones, Tom Jarrett, Robyn Smith

Abstract: We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47<z<2.85), luminous (11.7 < log L(bol) (Lsun) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid-infrared (M… ▽ More We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47<z<2.85), luminous (11.7 < log L(bol) (Lsun) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid-infrared (MIR) colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are detected at 870 microns, and we find that the sample has large mid- to far-infrared luminosity ratios consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are 24.7 < log P3.0 GHz (W/Hz) < 27.3, and all sources are radio-intermediate or radio-loud. BH mass estimates are 7.7 < log M(BH) (Msun) < 10.2. The rest frame 1-5 um SEDs are very similar to the "Hot DOGs" (Hot Dust Obscured Galaxies), and steeper (redder) than almost any other known extragalactic sources. ISM masses estimated for the ALMA detected sources are 9.9 < log M(ISM) (Msun) < 11.75 assuming a dust temperature of 30K. The cool dust emission is consistent with star formation rates (SFRs) reaching several thousand Msun/yr, depending on the assumed dust temperature, however we cannot rule out the alternative that the AGN powers all the emission in some cases. Our best constrained source has radiative transfer solutions with ~ equal contributions from an obscured AGN and a young (10-15 Myr) compact starburst. △ Less

Submitted 14 September, 2015; v1 submitted 1 September, 2015; originally announced September 2015.

Comments: 29 pages, 8 figures. To appear in Astrophysical Journal. Update on Sept 14 to correct the ALMA proposal id. to ADS/JAO.ALMA#2011.0.00397.S and to add a missing acknowledgement

Self-consistent 2-phase AGN torus models: SED library for observers

Authors: Ralf Siebenmorgen, Frank Heymann, Andreas Efstathiou

Abstract: We assume that dust near active galactic nuclei (AGN) is distributed in a torus-like geometry, which may be described by a clumpy medium or a homogeneous disk or as a combination of the two (i.e. a 2-phase medium). The dust particles considered are fluffy and have higher submillimeter emissivities than grains in the diffuse ISM. The dust-photon interaction is treated in a fully self-consistent thr… ▽ More We assume that dust near active galactic nuclei (AGN) is distributed in a torus-like geometry, which may be described by a clumpy medium or a homogeneous disk or as a combination of the two (i.e. a 2-phase medium). The dust particles considered are fluffy and have higher submillimeter emissivities than grains in the diffuse ISM. The dust-photon interaction is treated in a fully self-consistent three dimensional radiative transfer code. We provide an AGN library of spectral energy distributions (SEDs). Its purpose is to quickly obtain estimates of the basic parameters of the AGN, such as the intrinsic luminosity of the central source, the viewing angle, the inner radius, the volume filling factor and optical depth of the clouds, and the optical depth of the disk midplane, and to predict the flux at yet unobserved wavelengths. The procedure is simple and consists of finding an element in the library that matches the observations. We discuss the general properties of the models and in particular the 10mic. silicate band. The AGN library accounts well for the observed scatter of the feature strengths and wavelengths of the peak emission. AGN extinction curves are discussed and we find that there is no direct one-to-one link between the observed extinction and the wavelength dependence of the dust cross sections. We show that objects of the library cover the observed range of mid IR colors of known AGN. The validity of the approach is demonstrated by matching the SEDs of a number of representative objects: Four Seyferts and two quasars for which we present new Herschel photometry, two radio galaxies, and one hyperluminous infrared galaxy. Strikingly, for the five luminous objects we find pure AGN models fit the SED without a need to postulate starburst activity. △ Less

Submitted 18 August, 2015; originally announced August 2015.

Comments: A&A accepted by referee, AGN library available at http://www.eso.org/~rsiebenm/agn_models/index.html

Combining Physical galaxy models with radio observations to constrain the SFRs of high-z dusty star forming galaxies

Authors: B. Lo Faro, L. Silva, A. Franceschini, N. Miller, A. Efstathiou

Abstract: We complement our previous analysis of a sample of z~1-2 luminous and ultra-luminous infrared galaxies ((U)LIRGs), by adding deep VLA radio observations at 1.4 GHz to a large data-set from the far-UV to the sub-mm, including Spitzer and Herschel data. Given the relatively small number of (U)LIRGs in our sample with high S/N radio data, and to extend our study to a different family of galaxies, we… ▽ More We complement our previous analysis of a sample of z~1-2 luminous and ultra-luminous infrared galaxies ((U)LIRGs), by adding deep VLA radio observations at 1.4 GHz to a large data-set from the far-UV to the sub-mm, including Spitzer and Herschel data. Given the relatively small number of (U)LIRGs in our sample with high S/N radio data, and to extend our study to a different family of galaxies, we also include 6 well sampled near IR-selected BzK galaxies at z~1.5. From our analysis based on the radiative transfer spectral synthesis code GRASIL, we find that, while the IR luminosity may be a biased tracer of the star formation rate (SFR) depending on the age of stars dominating the dust heating, the inclusion of the radio flux offers significantly tighter constraints on SFR. Our predicted SFRs are in good agreement with the estimates based on rest-frame radio luminosity and the Bell (2003) calibration. The extensive spectro-photometric coverage of our sample allows us to set important constraints on the SF history of individual objects. For essentially all galaxies we find evidence for a rather continuous SFR and a peak epoch of SF preceding that of the observation by a few Gyrs. This seems to correspond to a formation redshift of z~5-6. We finally show that our physical analysis may affect the interpretation of the SFR-M* diagram, by possibly shifting, with respect to previous works, the position of the most dust obscured objects to higher M* and lower SFRs. △ Less

Submitted 10 December, 2014; originally announced December 2014.

Comments: 26 pages, 15 figures, 3 tables, accepted for publication in MNRAS on Dec. 4th, 2014

A multi-wavelength view of the central kiloparsec region in the Luminous Infrared Galaxy NGC1614

Authors: Rubén Herrero-Illana, Miguel Á. Pérez-Torres, Almudena Alonso-Herrero, Antxon Alberdi, Luis Colina, Andreas Efstathiou, Lorena Hernández-García, Daniel Miralles-Caballero, Petri Väisänen, Christopher C. Packham, Vinesh Rajpaul, Albert A. Zijlstra

Abstract: The Luminous Infrared Galaxy NGC1614 hosts a prominent circumnuclear ring of star formation. However, the nature of the dominant emitting mechanism in its central ~100 pc is still under debate. We present sub-arcsecond angular resolution radio, mid-infrared, Pa-alpha, optical, and X-ray observations of NGC1614, aimed at studying in detail both the circumnuclear ring and the nuclear region. The 8.4… ▽ More The Luminous Infrared Galaxy NGC1614 hosts a prominent circumnuclear ring of star formation. However, the nature of the dominant emitting mechanism in its central ~100 pc is still under debate. We present sub-arcsecond angular resolution radio, mid-infrared, Pa-alpha, optical, and X-ray observations of NGC1614, aimed at studying in detail both the circumnuclear ring and the nuclear region. The 8.4 GHz continuum emission traced by the Very Large Array (VLA) and the Gemini/T-ReCS 8.7 micron emission, as well as the Pa-alpha line emission, show remarkable morphological similarities within the star-forming ring, suggesting that the underlying emission mechanisms are tightly related. We used an HST/NICMOS Pa-alpha map of similar resolution to our radio maps to disentangle the thermal free-free and non-thermal synchrotron radio emission, from which we obtained the intrinsic synchrotron power-law for each individual region within the central kpc of NGC1614. The radio ring surrounds a relatively faint, steep-spectrum source at the very center of the galaxy, suggesting that the central source is not powered by an AGN, but rather by a compact (r < 90 pc) starburst. Chandra X-ray data also show that the central kpc region is dominated by starburst activity, without requiring the existence of an AGN. We also used publicly available infrared data to model-fit the spectral energy distribution of both the starburst ring and a putative AGN in NGC1614. In summary, we conclude that there is no need to invoke an AGN to explain the observed bolometric properties of the galaxy. △ Less

Submitted 27 March, 2014; originally announced March 2014.

Comments: 13 pages, 7 figures, 5 tables. Accepted for publication in ApJ