-
Galaxy interactions are the dominant trigger for local type 2 quasars
Authors:
J. C. S. Pierce,
C. Tadhunter,
C. Ramos Almeida,
P. Bessiere,
J. V. Heaton,
S. L. Ellison,
G. Speranza,
Y. Gordon,
C. O'Dea,
L. Grimmett,
L. Makrygianni
Abstract:
The triggering mechanism for the most luminous, quasar-like active galactic nuclei (AGN) remains a source of debate, with some studies favouring triggering via galaxy mergers, but others finding little evidence to support this mechanism. Here, we present deep Isaac Newton Telescope/Wide Field Camera imaging observations of a complete sample of 48 optically-selected type 2 quasars $-$ the QSOFEED s…
▽ More
The triggering mechanism for the most luminous, quasar-like active galactic nuclei (AGN) remains a source of debate, with some studies favouring triggering via galaxy mergers, but others finding little evidence to support this mechanism. Here, we present deep Isaac Newton Telescope/Wide Field Camera imaging observations of a complete sample of 48 optically-selected type 2 quasars $-$ the QSOFEED sample (L$_{\rm [OIII]}>$10$^{8.5}$$L_{\odot}$; $z < 0.14$). Based on visual inspection by eight classifiers, we find clear evidence that galaxy interactions are the dominant triggering mechanism for quasar activity in the local universe, with 65$^{+6}_{-7}$ per cent of the type 2 quasar hosts showing morphological features consistent with galaxy mergers or encounters, compared with only 22$^{+5}_{-4}$ per cent of a stellar-mass- and redshift-matched comparison sample of non-AGN galaxies $-$ a 5$σ$ difference. The type 2 quasar hosts are a factor 3.0$^{+0.5}_{-0.8}$ more likely to be morphologically disturbed than their matched non-AGN counterparts, similar to our previous results for powerful 3CR radio AGN of comparable [OIII] emission-line luminosity and redshift. In contrast to the idea that quasars are triggered at the peaks of galaxy mergers as the two nuclei coalesce, and only become visible post-coalescence, the majority of morphologically-disturbed type 2 quasar sources in our sample are observed in the pre-coalescence phase (61$^{+8}_{-9}$ per cent). We argue that much of the apparent ambiguity that surrounds observational results in this field is a result of differences in the surface brightness depths of the observations, combined with the effects of cosmological surface brightness dimming.
△ Less
Submitted 27 March, 2023;
originally announced March 2023.
-
Do AGN triggering mechanisms vary with radio power? II. The importance of mergers as a function of radio power and optical luminosity
Authors:
J. C. S. Pierce,
C. N. Tadhunter,
Y. Gordon,
C. Ramos Almeida,
S. L. Ellison,
C. O'Dea,
L. Grimmett,
L. Makrygianni,
P. S. Bessiere,
P. Doña Girón
Abstract:
Investigation of the triggering mechanisms of radio AGN is important for improving our general understanding of galaxy evolution. In the first paper in this series, detailed morphological analysis of high-excitation radio galaxies (HERGs) with intermediate radio powers suggested that the importance of triggering via galaxy mergers and interactions increases strongly with AGN radio power and weakly…
▽ More
Investigation of the triggering mechanisms of radio AGN is important for improving our general understanding of galaxy evolution. In the first paper in this series, detailed morphological analysis of high-excitation radio galaxies (HERGs) with intermediate radio powers suggested that the importance of triggering via galaxy mergers and interactions increases strongly with AGN radio power and weakly with optical emission-line luminosity. Here, we use an online classification interface to expand our morphological analysis to a much larger sample of 155 active galaxies (3CR radio galaxies, radio-intermediate HERGs and Type 2 quasars) that covers a broad range in both 1.4 GHz radio power and [OIII]$λ$5007 emission-line luminosity. All active galaxy samples are found to exhibit excesses in their rates of morphological disturbance relative to 378 stellar-mass- and redshift-matched non-active control galaxies classified randomly and blindly alongside them. These excesses are highest for the 3CR HERGs (4.7$σ$) and Type 2 quasar hosts (3.7$σ$), supporting the idea that galaxy mergers provide the dominant triggering mechanism for these subgroups. When the full active galaxy sample is considered, there is clear evidence to suggest that the enhancement in the rate of disturbance relative to the controls increases strongly with [OIII]$λ$5007 emission-line luminosity but not with 1.4 GHz radio power. Evidence that the dominant AGN host types change from early-type galaxies at high radio powers to late-type galaxies at low radio powers is also found, suggesting that triggering by secular, disk-based processes holds more importance for lower-power radio AGN.
△ Less
Submitted 4 November, 2021;
originally announced November 2021.
-
The post-Herschel view of intrinsic AGN emission: constructing templates for galaxy and AGN emission at IR wavelengths
Authors:
E. Bernhard,
C. Tadhunter,
J. R. Mullaney,
L. P. Grimmett,
D. J. Rosario,
D. M. Alexander
Abstract:
Measuring the star-forming properties of AGN hosts is key to our understanding of galaxy formation and evolution. However, this topic remains debated, partly due to the difficulties in separating the infrared (i.e. 1--1000 $μ$m) emission into AGN and star-forming components. Taking advantage of archival far-infrared data from Herschel, we present a new set of AGN and galaxy infrared templates, and…
▽ More
Measuring the star-forming properties of AGN hosts is key to our understanding of galaxy formation and evolution. However, this topic remains debated, partly due to the difficulties in separating the infrared (i.e. 1--1000 $μ$m) emission into AGN and star-forming components. Taking advantage of archival far-infrared data from Herschel, we present a new set of AGN and galaxy infrared templates, and introduce the spectral energy distribution fitting code IRAGNSEP. Both can be used to measure infrared host galaxy properties, free of AGN contamination. To build these, we used a sample of 100 local ($z$ < 0.3), low-to-high luminosity AGNs (i.e. $L_{\rm bol}~\sim~10^{42--46}~\rm erg~s^{-1}$), selected from the 105-month Swift - BAT X-ray survey, which have archival Spitzer - IRS spectra and Herschel photometry. We first built a set of seven galaxy templates using a sample of 55 star-forming galaxies selected via infrared diagnostics. Using these templates, combined with a flexible model for the AGN contribution, we extracted the intrinsic infrared emission of our AGN sample. We further demonstrate that we can reduce the diversity in the intrinsic shapes of AGN spectral energy distributions down to a set of three AGN templates, of which two represent AGN continuum, and one represents silicate emission. Our results indicate that, on average, the contribution of AGNs to the far-infrared ($λ~\gtrsim$ 50 $μ$m) is not as high as suggested by some recent work. We further show that the need for two infrared AGN continuum templates could be related to nuclear obscuration, where one of our templates appears dominated by the emission of the extended polar dust.
△ Less
Submitted 10 February, 2021;
originally announced February 2021.
-
Exoplanet detection and its dependence on stochastic sampling of the stellar Initial Mass Function
Authors:
Amy L. Bottrill,
Molly E. Haigh,
Madeleine R. A. Hole,
Sarah C. M. Theakston,
Rosa B. Allen,
Liam P. Grimmett,
Richard J. Parker
Abstract:
Young Moving Groups (YMGs) are close (<100pc), coherent collections of young (<100Myr) stars that appear to have formed in the same star-forming molecular cloud. As such we would expect their individual initial mass functions (IMFs) to be similar to other star-forming regions, and by extension the Galactic field. Their close proximity to the Sun and their young ages means that YMGs are promising l…
▽ More
Young Moving Groups (YMGs) are close (<100pc), coherent collections of young (<100Myr) stars that appear to have formed in the same star-forming molecular cloud. As such we would expect their individual initial mass functions (IMFs) to be similar to other star-forming regions, and by extension the Galactic field. Their close proximity to the Sun and their young ages means that YMGs are promising locations to search for young forming exoplanets. However, due to their low numbers of stars, stochastic sampling of the IMF means their stellar populations could vary significantly. We determine the range of planet-hosting stars (spectral types A, G and M) possible from sampling the IMF multiple times, and find that some YMGs appear deficient in M-dwarfs. We then use these data to show that the expected probability of detecting terrestrial magma ocean planets is highly dependent on the exact numbers of stars produced through stochastic sampling of the IMF.
△ Less
Submitted 4 June, 2020;
originally announced June 2020.
-
The evolving AGN duty cycle in galaxies since z$\sim$3 as encoded in the X-ray luminosity function
Authors:
I. Delvecchio,
E. Daddi,
J. Aird,
J. R. Mullaney,
E. Bernhard,
L. P. Grimmett,
R. Carraro,
A. Cimatti,
G. Zamorani,
N. Caplar,
F. Vito,
D. Elbaz,
G. Rodighiero
Abstract:
We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($λ_{\rm EDD}$) distributions, normalised via empirical…
▽ More
We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($λ_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $λ_{\rm EDD}$ distribution, and with a positive shift of the break $λ^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($λ_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $λ_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly.
△ Less
Submitted 20 February, 2020;
originally announced February 2020.
-
A binning-free method reveals a continuous relationship between galaxies' AGN power and offset from main sequence
Authors:
L. P. Grimmett,
J. R. Mullaney,
E. P. Bernhard,
C. M. Harrison,
D. M. Alexander,
F. Stanley,
V. A. Masoura,
K. Walters
Abstract:
Studies investigating the relationship between AGN power and the star formation rates (SFRs) of their host galaxies often rely on averaging techniques -- such as stacking -- to incorporate information from non-detections. However, averages, and especially means, can be strongly affected by outliers and can therefore give a misleading indication of the "typical" case. Recently, a number of studies…
▽ More
Studies investigating the relationship between AGN power and the star formation rates (SFRs) of their host galaxies often rely on averaging techniques -- such as stacking -- to incorporate information from non-detections. However, averages, and especially means, can be strongly affected by outliers and can therefore give a misleading indication of the "typical" case. Recently, a number of studies have taken a step further by binning their sample in terms of AGN power (approximated by the 2-10keV luminosity of the AGN), and investigating how the SFR distribution differs between these bins. These bin thresholds are often weakly motivated, and binning implicitly assumes that sources within the same bin have similar (or even identical) properties. In this paper, we investigate whether the distribution of host SFRs -- relative to the locus of the star-forming main sequence (i.e., $R_{\rm MS}$) -- changes continuously as a function of AGN power. We achieve this by using a hierarchical Bayesian model that completely removes the need to bin in AGN power. In doing so, we find strong evidence that the $R_{\rm MS}$ distribution changes with 2-10keV X-ray luminosity. The results suggest that higher X_ray luminosity AGNs have a tighter physical connection to the star-forming process than lower X-ray luminosity AGNs, at least within the $0.8 < z < 1.2$ redshift range considered here.
△ Less
Submitted 4 May, 2020; v1 submitted 30 January, 2020;
originally announced January 2020.
-
The galaxy's gas content regulated by the dark matter halo mass results in a super-linear M$_{\rm BH}$-M$_{\star}$ relation
Authors:
I. Delvecchio,
E. Daddi,
F. Shankar,
J. R. Mullaney,
G. Zamorani,
J. Aird,
E. Bernhard,
A. Cimatti,
D. Elbaz,
M. Giavalisco,
L. P. Grimmett
Abstract:
Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS…
▽ More
Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $ε\sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$).
△ Less
Submitted 17 October, 2019;
originally announced October 2019.
-
Dynamical evolution of star-forming regions: III. Unbound stars and predictions for Gaia
Authors:
Christina Schoettler,
Richard J. Parker,
Becky Arnold,
Liam P. Grimmett,
Jos de Bruijne,
Nicholas J. Wright
Abstract:
We use $N$-body simulations to probe the early phases of the dynamical evolution of star-forming regions and focus on mass and velocity distributions of unbound stars. In this parameter space study, we vary the initial virial ratio and degree of spatial and kinematic substructure and analyse the fraction of stars that become unbound in two different mass classes (above and below 8 M$_{\odot}$). We…
▽ More
We use $N$-body simulations to probe the early phases of the dynamical evolution of star-forming regions and focus on mass and velocity distributions of unbound stars. In this parameter space study, we vary the initial virial ratio and degree of spatial and kinematic substructure and analyse the fraction of stars that become unbound in two different mass classes (above and below 8 M$_{\odot}$). We find that the fraction of unbound stars differs depending on the initial conditions. After 10 Myr, in initially highly subvirial, substructured simulations, the high-mass and lower-mass unbound fractions are similar at $\sim$23 per cent. In initially virialised, substructured simulations, we find only $\sim$16 per cent of all high-mass stars are unbound, whereas $\sim$37 per cent of all lower-mass stars are. The velocity distributions of unbound stars only show differences for extremely different initial conditions. The distributions are dominated by large numbers of lower-mass stars becoming unbound just above the escape velocity of $\sim$3 km s$^{-1}$ with unbound high-mass stars moving faster on average than lower-mass unbound stars. We see no high-mass runaway stars (velocity > 30 km s$^{-1}$) from any of our initial conditions and only an occasional lower-mass runaway star from initially subvirial/substructured simulations. In our simulations, we find a small number of lower-mass walkaway stars (with velocity 5-30 km s$^{-1}$) from all of our initial conditions. These walkaway stars should be observable around many nearby star-forming regions with Gaia.
△ Less
Submitted 24 May, 2019;
originally announced May 2019.
-
Revealing the differences in the SMBH accretion rate distributions of starburst and non-starburst galaxies
Authors:
L. P. Grimmett,
J. R. Mullaney,
S. Jin,
E. Bernhard,
E. Daddi,
K. Walters
Abstract:
We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from no…
▽ More
We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from non-detections in the form of upper limits. When we use our inferred distributions to calculate the ratios of high to low sLx AGN (corresponding to above and below $0.1λ_{\text{Edd}}$, respectively) we find that starbursts have significantly higher proportions of high sLx AGN compared to their non-starburst counterparts. These findings help explain the increase in average X-ray luminosity in bins of increasing SFR reported by previous studies.
△ Less
Submitted 16 April, 2019;
originally announced April 2019.
-
Inferring a difference in the star-forming properties of lower versus higher X-ray luminosity AGNs
Authors:
E. Bernhard,
L. P. Grimmett,
J. R. Mullaney,
E. Daddi,
C. Tadhunter,
S. Jin
Abstract:
We explore the distribution of RMS=SFR/SFR_MS (where SFR_MS is the star formation rate of "Main Sequence" star-forming galaxies) for AGN hosts at z=1. We split our sample into two bins of X-ray luminosity divided at Lx=2x10^43erg s-1 to investigate whether the RMS distribution changes as a function of AGN power. Our main results suggest that, when the RMS distribution of AGN hosts is modelled as a…
▽ More
We explore the distribution of RMS=SFR/SFR_MS (where SFR_MS is the star formation rate of "Main Sequence" star-forming galaxies) for AGN hosts at z=1. We split our sample into two bins of X-ray luminosity divided at Lx=2x10^43erg s-1 to investigate whether the RMS distribution changes as a function of AGN power. Our main results suggest that, when the RMS distribution of AGN hosts is modelled as a log-normal distribution (i.e. the same shape as that of MS galaxies), galaxies hosting more powerful X-ray AGNs (i.e. Lx>2x10^43erg s-1) display a narrower RMS distribution that is shifted to higher values compared to their lower Lx counterparts. In addition, we find that more powerful X-ray AGNs have SFRs that are more consistent with that of MS galaxies compared to lower Lx AGNs. Despite this, the mean SFRs (as opposed to RMS) measured from these distributions are consistent with the previously observed flat relationship between SFR and Lx. Our results suggest that the typical star-forming properties of AGN hosts change with Lx , and that more powerful AGNs typically reside in more MS-like star-forming galaxies compared to lower Lx AGNs.
△ Less
Submitted 14 November, 2018;
originally announced November 2018.
-
Evidence for a mass-dependent AGN Eddington ratio distribution via the flat relationship between SFR and AGN luminosity
Authors:
E. Bernhard,
J. R. Mullaney,
J. Aird,
R. C. Hickox,
M. L. Jones,
F. Stanley,
L. P. Grimmett,
E. Daddi
Abstract:
The lack of a strong correlation between AGN X-ray luminosity ($L_X$; a proxy for AGN power) and the star formation rate (SFR) of their host galaxies has recently been attributed to stochastic AGN variability. Studies using population synthesis models have incorporated this by assuming a broad, universal (i.e. does not depend on the host galaxy properties) probability distribution for AGN specific…
▽ More
The lack of a strong correlation between AGN X-ray luminosity ($L_X$; a proxy for AGN power) and the star formation rate (SFR) of their host galaxies has recently been attributed to stochastic AGN variability. Studies using population synthesis models have incorporated this by assuming a broad, universal (i.e. does not depend on the host galaxy properties) probability distribution for AGN specific X-ray luminosities (i.e. the ratio of $L_X$ to host stellar mass; a common proxy for Eddington ratio). However, recent studies have demonstrated that this universal Eddington ratio distribution fails to reproduce the observed X-ray luminosity functions beyond z$\sim$1.2. Furthermore, empirical studies have recently shown that the Eddington ratio distribution may instead depend upon host galaxy properties, such as SFR and/or stellar mass. To investigate this further we develop a population synthesis model in which the Eddington ratio distribution is different for star-forming and quiescent host galaxies. We show that, although this model is able to reproduce the observed X-ray luminosity functions out to z$\sim$2, it fails to simultaneously reproduce the observed flat relationship between SFR and X-ray luminosity. We can solve this, however, by incorporating a mass dependency in the AGN Eddington ratio distribution for star-forming host galaxies. Overall, our models indicate that a relative suppression of low Eddington ratios ($λ_{\rm Edd}\lesssim$0.1) in lower mass galaxies (M<$10^{10-11}$Msun) is required to reproduce both the observed X-ray luminosity functions and the observed flat SFR/X-ray relationship.
△ Less
Submitted 19 February, 2018; v1 submitted 23 January, 2018;
originally announced January 2018.
