-
The LSST era of supermassive black holes accretion-disk reverberation mapping
Authors:
Andjelka B. Kovacevic,
Viktor Radovic,
Dragana Ilic,
Luka C. Popovic,
Roberto J. Assef,
Paula Sanchez-Saez,
Robert Nikutta,
Claudia M. Raiteri,
Ilsang Yoon,
Yasaman Homayouni,
Yan-Rong Li,
Neven Caplar,
Bozena Czerny,
Swayamtrupta Panda,
Claudio Ricci,
Isidora Jankov,
Hermine Landt,
Christian Wolf,
Jelena Kovacevic-Dojcinovic,
Masa Lakicevic,
Djorđe Savic,
Oliver Vince,
Sasa Simic,
Iva Cvorovic-Hajdinjak,
Sladjana Marceta-Mandic
Abstract:
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will detect an unprecedentedly large sample of actively accreting supermassive black holes with typical accretion disk (AD) sizes of a few light days. This brings us to face challenges in the reverberation mapping (RM) measurement of AD sizes in active galactic nuclei (AGNs) using interband continuum delays. We examine the effect…
▽ More
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will detect an unprecedentedly large sample of actively accreting supermassive black holes with typical accretion disk (AD) sizes of a few light days. This brings us to face challenges in the reverberation mapping (RM) measurement of AD sizes in active galactic nuclei (AGNs) using interband continuum delays. We examine the effect of LSST cadence strategies on AD RM using our metric AGNTimeLagMetric. It accounts for redshift, cadence, the magnitude limit, and magnitude corrections for dust extinction. Running our metric on different LSST cadence strategies, we produce an atlas of the performance estimations for LSST photometric RM measurements. We provide an upper limit on the estimated number of quasars for which the AD time lag can be computed within 0<z<7 using the features of our metric. We forecast that the total counts of such objects will increase as the mean sampling rate of the survey decreases. The AD time lag measurements are expected for >1000 sources in each Deep Drilling field (DDF, 10 sq. deg) in any filter, with the redshift distribution of these sources peaking at z~1. We find the LSST observation strategies with a good cadence (~ 5 days) and a long cumulative season (~9 yr), as proposed for LSST DDF, are favored for the AD size measurement. We create synthetic LSST light curves for the most suitable DDF cadences and determine RM time lags to demonstrate the impact of the best cadences based on the proposed metric.
△ Less
Submitted 12 August, 2022;
originally announced August 2022.
-
Stochastic Modelling of Star Formation Histories III. Constraints from Physically-Motivated Gaussian Processes
Authors:
Kartheik G. Iyer,
Joshua S. Speagle,
Neven Caplar,
John C. Forbes,
Eric Gawiser,
Joel Leja,
Sandro Tacchella
Abstract:
Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow,…
▽ More
Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow, equilibrium, and dynamical processes influencing GMC creation and destruction using the power spectral density (PSD) formalism. We use the connection between the PSD and auto-covariance function (ACF) for general stochastic processes to reformulate this model as an auto-covariance function, which we use to model variability in galaxy star formation histories (SFHs) using physically-motivated Gaussian Processes in log SFR space. Using stellar population synthesis models, we then explore how changes in model stochasticity can affect spectral signatures across galaxy populations with properties similar to the Milky Way and present-day dwarfs as well as at higher redshifts. We find that, even at fixed scatter, perturbations to the stochasticity model (changing timescales vs overall variability) leave unique spectral signatures across both idealized and more realistic galaxy populations. Distributions of spectral features including H$α$ and UV-based SFR indicators, H$δ$ and Ca-H,K absorption line strengths, D$_n$(4000) and broadband colors provide testable predictions for galaxy populations from present and upcoming surveys with Hubble, Webb \& Roman. The Gaussian process SFH framework provides a fast, flexible implementation of physical covariance models for the next generation of SED modeling tools. Code to reproduce our results can be found at https://github.com/kartheikiyer/GP-SFH
△ Less
Submitted 11 August, 2022;
originally announced August 2022.
-
From Data to Software to Science with the Rubin Observatory LSST
Authors:
Katelyn Breivik,
Andrew J. Connolly,
K. E. Saavik Ford,
Mario Jurić,
Rachel Mandelbaum,
Adam A. Miller,
Dara Norman,
Knut Olsen,
William O'Mullane,
Adrian Price-Whelan,
Timothy Sacco,
J. L. Sokoloski,
Ashley Villar,
Viviana Acquaviva,
Tomas Ahumada,
Yusra AlSayyad,
Catarina S. Alves,
Igor Andreoni,
Timo Anguita,
Henry J. Best,
Federica B. Bianco,
Rosaria Bonito,
Andrew Bradshaw,
Colin J. Burke,
Andresa Rodrigues de Campos
, et al. (75 additional authors not shown)
Abstract:
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the po…
▽ More
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science.
To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems.
This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science.
△ Less
Submitted 4 August, 2022;
originally announced August 2022.
-
Dwarf AGNs from Variability for the Origins of Seeds (DAVOS): Intermediate-mass black hole demographics from optical synoptic surveys
Authors:
Colin J. Burke,
Yue Shen,
Xin Liu,
Priyamvada Natarajan,
Neven Caplar,
Jillian M. Bellovary,
Z. Franklin Wang
Abstract:
We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black h…
▽ More
We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black hole (BH) mass to generate mock light curves. Adopting several different models for the BH occupation function, including one for off-nuclear IMBHs, we quantify differences in the predicted local AGN mass and luminosity functions in dwarf galaxies. As a result, we are able to model the variable fraction of AGNs as a function of physical host properties, such as host galaxy stellar mass, in the presence of complex selection effects. We find that our adopted occupation fractions for the "heavy" and "light" initial BH seeding scenarios can be distinguished with variability data at the $2-3 σ$ level for galaxy host stellar masses below $\sim 10^8 M_\odot$ with the Vera C. Rubin Observatory. We demonstrate the prevalence of a selection bias whereby recovered IMBH masses fall, on average, above the predicted value from the local host galaxy - BH mass scaling relation with the strength of the bias dependent on the survey sensitivity. The methodology developed in this work can be used more broadly to forecast and correct for selection effects for AGN demographic studies in synoptic surveys. Finally, we show that a targeted $\sim$ hourly cadence program over a few nights with the Rubin Observatory can provide strong constraints on IMBH masses given their expected rapid variability timescales.
△ Less
Submitted 10 August, 2023; v1 submitted 8 July, 2022;
originally announced July 2022.
-
On possible proxies of AGN light curves cadence selection in future time domain surveys
Authors:
Andjelka B. Kovacevic,
Dragana Ilic,
Luka C. Popovic,
Viktor Radovic,
Isidora Jankov,
Ilsang Yoon,
Neven Caplar,
Iva Cvorovic-Hajdinjak,
Sasa Simic
Abstract:
Motivated by upcoming photometric and spectroscopic surveys (Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), Manuakea Spectroscopic Explorer), we design the statistical proxies to measure the cadence effects on active galactic nuclei (AGN) variability-observables (time-lags, periodicity, and structure-function (SF)). We constructed a multiple-regression model to statistically ide…
▽ More
Motivated by upcoming photometric and spectroscopic surveys (Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), Manuakea Spectroscopic Explorer), we design the statistical proxies to measure the cadence effects on active galactic nuclei (AGN) variability-observables (time-lags, periodicity, and structure-function (SF)). We constructed a multiple-regression model to statistically identify the cadence-formal error pattern knowing AGN time-lags and periodicity from different surveys. We defined the simple metric for the SF's properties, accounting for the 'observed' SF's deviation relative to those obtained from the homogenously-sampled light curves. We tested the regression models on different observing strategies: the optical dataset of long light-curves of eight AGN with peculiarities and the artificial datasets based on several idealized and LSST-like cadences. The SFs metric is assessed on synthetic datasets. The regression models (for both data types) predict similar cadences for time-lags and oscillation detection, whereas for light curves with low variability ($\sim 10\%$), cadences for oscillation detection differ. For higher variability ($\sim20\%$), predicted cadences are larger than for $F_{var}\sim 10\%$. The predicted cadences are decreasing with redshift. SFs with dense and homogenous cadences are more likely to behave similarly. SFs with oscillatory signals are sensitive to the cadences, possibly impacting LSST-like operation strategy.
The proposed proxies can help to select spectroscopic and photometric-surveys cadence strategies, and they will be tested further in larger samples of objects.
△ Less
Submitted 31 May, 2021;
originally announced May 2021.
-
LSST AGN SC Cadence Note: Two metrics on AGN variability observable
Authors:
Andjelka Kovacevic,
Dragana Ilic,
Isidora Jankov,
Luka C. Popovic,
Ilsang Yoon,
Viktor Radovic,
Neven Caplar,
Iva Cvorovic-Hajdinjak
Abstract:
We have developed two metrics related to AGN variability observables (time-lags, periodicity, and Structure Function (SF)) to evaluate LSST OpSim FBS 1.5, 1.6, 1.7 performance in AGN time-domain analysis. For this purpose, we generate an ensemble of AGN light curves based on AGN empirical relations and LSST OpSim cadences. Although our metrics show that denser LSST cadences produce more reliable t…
▽ More
We have developed two metrics related to AGN variability observables (time-lags, periodicity, and Structure Function (SF)) to evaluate LSST OpSim FBS 1.5, 1.6, 1.7 performance in AGN time-domain analysis. For this purpose, we generate an ensemble of AGN light curves based on AGN empirical relations and LSST OpSim cadences. Although our metrics show that denser LSST cadences produce more reliable time-lag, periodicity, and SF measurements, the discrepancies in the performance between different LSST OpSim cadences are not drastic based on Kullback-Leibler divergence. This is complementary to Yu and Richards results on DCR and SF metrics, extending them to include the point of view of AGN variability.
△ Less
Submitted 26 May, 2021;
originally announced May 2021.
-
The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution
Authors:
Kartheik G. Iyer,
Sandro Tacchella,
Shy Genel,
Christopher C. Hayward,
Lars Hernquist,
Alyson M. Brooks,
Neven Caplar,
Romeel Davé,
Benedikt Diemer,
John C. Forbes,
Eric Gawiser,
Rachel S. Somerville,
Tjitske K. Starkenburg
Abstract:
Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIR…
▽ More
Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
Stochastic modelling of star-formation histories II: star-formation variability from molecular clouds and gas inflow
Authors:
Sandro Tacchella,
John C. Forbes,
Neven Caplar
Abstract:
A key uncertainty in galaxy evolution is the physics regulating star formation, ranging from small-scale processes related to the life-cycle of molecular clouds within galaxies to large-scale processes such as gas accretion onto galaxies. We study the imprint of such processes on the time-variability of star formation with an analytical approach tracking the gas mass of galaxies ("regulator model"…
▽ More
A key uncertainty in galaxy evolution is the physics regulating star formation, ranging from small-scale processes related to the life-cycle of molecular clouds within galaxies to large-scale processes such as gas accretion onto galaxies. We study the imprint of such processes on the time-variability of star formation with an analytical approach tracking the gas mass of galaxies ("regulator model"). Specifically, we quantify the strength of the fluctuation in the star-formation rate (SFR) on different timescales, i.e. the power spectral density (PSD) of the star-formation history, and connect it to gas inflow and the life-cycle of molecular clouds. We show that in the general case the PSD of the SFR has three breaks, corresponding to the correlation time of the inflow rate, the equilibrium timescale of the gas reservoir of the galaxy, and the average lifetime of individual molecular clouds. On long and intermediate timescales (relative to the dynamical timescale of the galaxy), the PSD is typically set by the variability of the inflow rate and the interplay between outflows and gas depletion. On short timescales, the PSD shows an additional component related to the life-cycle of molecular clouds, which can be described by a damped random walk with a power-law slope of $β\approx2$ at high frequencies with a break near the average cloud lifetime. We discuss star-formation "burstiness" in a wide range of galaxy regimes, study the evolution of galaxies about the main sequence ridgeline, and explore the applicability of our method for understanding the star-formation process on cloud-scale from galaxy-integrated measurements.
△ Less
Submitted 6 July, 2020; v1 submitted 16 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.
-
Observational nonstationarity of AGN variability: The only way to go is down!
Authors:
Neven Caplar,
Theodore Pena,
Sean D. Johnson,
Jenny E. Greene
Abstract:
To gain insights into long-term Active Galactic Nuclei (AGN) variability, we analyze an AGN sample from the Sloan Digital Sky Survey (SDSS) and compare their photometry with observations from the Hyper Suprime-Cam survey (HSC) observed $\langle 14.85 \rangle$ years after SDSS. On average, the AGN are fainter in HSC than SDSS. We demonstrate that the difference is not due to subtle differences in t…
▽ More
To gain insights into long-term Active Galactic Nuclei (AGN) variability, we analyze an AGN sample from the Sloan Digital Sky Survey (SDSS) and compare their photometry with observations from the Hyper Suprime-Cam survey (HSC) observed $\langle 14.85 \rangle$ years after SDSS. On average, the AGN are fainter in HSC than SDSS. We demonstrate that the difference is not due to subtle differences in the SDSS versus HSC filters or photometry. The decrease in mean brightness is redshift dependent, consistent with expectations for a change that is a function of the rest-frame time separation between observations. At a given redshift, the mean decrease in brightness is stronger for more luminous AGN and for objects with longer time separation between measurements. We demonstrate that the dependence on redshift and luminosity of measured mean brightness decrease is consistent with simple models of Eddington ratio variability in AGN on long (Myr, Gyr) timescales. We show how our results can be used to constrain the variability and demographic properties of AGN populations.
△ Less
Submitted 22 January, 2020; v1 submitted 13 January, 2020;
originally announced January 2020.
-
A forward modelling approach to AGN variability -- Method description and early applications
Authors:
Lia F. Sartori,
Benny Trakhtenbrot,
Kevin Schawinski,
Neven Caplar,
Ezequiel Treister,
Ce Zhang
Abstract:
We present a numerical framework for the variability of active galactic nuclei (AGN), which links the variability of AGN over a broad range of timescales and luminosities to the observed properties of the AGN population as a whole, and particularly the Eddington ratio distribution function (ERDF). We have implemented our framework on GPU architecture, relying on previously published time series ge…
▽ More
We present a numerical framework for the variability of active galactic nuclei (AGN), which links the variability of AGN over a broad range of timescales and luminosities to the observed properties of the AGN population as a whole, and particularly the Eddington ratio distribution function (ERDF). We have implemented our framework on GPU architecture, relying on previously published time series generating algorithms. After extensive tests that characterise several intrinsic and numerical aspects of the simulations, we describe some applications used for current and future time domain surveys and for the study of extremely variable sources (e.g., "changing look" or flaring AGN). Specifically, we define a simulation setup which reproduces the AGN variability observed in the PTF/iPTF survey, and use it to forward model longer light curves of the kind that may be observed within the LSST main survey. Thanks to our effcient implementations, these simulations are able to cover for example over 1 Myr with a roughly weekly cadence. We envision that this framework will become highly valuable to prepare for, and best exploit, data from upcoming time domain surveys, such as for example LSST.
△ Less
Submitted 13 September, 2019;
originally announced September 2019.
-
Stochastic modeling of star-formation histories I: the scatter of the star-forming main sequence
Authors:
Neven Caplar,
Sandro Tacchella
Abstract:
We present a framework for modelling the star-formation histories of galaxies as a stochastic process. We define this stochastic process through a power spectrum density with a functional form of a broken power-law. Star-formation histories are correlated on short timescales, the strength of this correlation described by a power-law slope, $α$, and they decorrelate to resemble white noise over a t…
▽ More
We present a framework for modelling the star-formation histories of galaxies as a stochastic process. We define this stochastic process through a power spectrum density with a functional form of a broken power-law. Star-formation histories are correlated on short timescales, the strength of this correlation described by a power-law slope, $α$, and they decorrelate to resemble white noise over a timescale that is proportional to the timescale of the break in the power spectrum density, $τ_{\rm break}$. We use this framework to explore the properties of the stochastic process that, we assume, gives rise to the log-normal scatter about the relationship between star-formation rate and stellar mass, the so-called galaxy star-forming main sequence. Specifically, we show how the measurements of the normalisation and width ($σ_{\rm MS}$) of the main sequence, measured in several passbands that probe different timescales, give a constraint on the parameters of the underlying power spectrum density. We first derive these results analytically for a simplified case where we model observations by averaging over the recent star-formation history. We then run numerical simulations to find results for more realistic observational cases. As a proof of concept, we use observational estimates of the main sequence scatter at $z\sim0$ and $M_{\star}\approx10^{10}~M_{\odot}$ measured in H$α$, UV+IR and the u-band, and show that combination of these point to $τ_{\rm break}=178^{+104}_{-66}$ Myr, when assuming $α=2$. This implies that star-formation histories of galaxies lose "memory" of their previous activity on a timescale of $\sim200$ Myr, highlighting the importance of baryonic effects that act over the dynamical timescales of galaxies.
△ Less
Submitted 22 January, 2019;
originally announced January 2019.
-
AGN evolution from galaxy evolution viewpoint - II
Authors:
Neven Caplar,
Simon J. Lilly,
Benny Trakhtenbrot
Abstract:
In order to relate the observed evolution of the galaxy stellar mass function and the luminosity function of active galactic nuclei (AGN), we explore a co-evolution scenario in which AGN are associated only with the very last phases of the star-forming life of a galaxy. We derive analytically the connections between the parameters of the observed quasar luminosity functions and galaxy mass functio…
▽ More
In order to relate the observed evolution of the galaxy stellar mass function and the luminosity function of active galactic nuclei (AGN), we explore a co-evolution scenario in which AGN are associated only with the very last phases of the star-forming life of a galaxy. We derive analytically the connections between the parameters of the observed quasar luminosity functions and galaxy mass functions. The $(m_{\rm bh}/m_{*})_{Qing}$ associated with quenching is given by the ratio of the global black hole accretion rate density (BHARD) and star-formation rate density (SFRD) at the epoch in question. Observational data on the SFRD and BHARD suggests $(m_{\rm bh}/m_{*})_{Qing} \propto (1+z)^{1.5}$ below redshift 2. This evolution reproduces the observed mass-luminosity plane of SDSS quasars, and also reproduces the local $m_{\rm bh}/m_{*}$ relation in passive galaxies. The characteristic Eddington ratio, $λ^*$, is derived from both the BHARD/SFRD ratio and the evolving $L^*$ of the AGN population. This increases up to $z \sim 2$ as $λ^* \propto (1+z)^{2.5}$ but at higher redshifts, $λ^*$ stabilizes at the physically interesting Eddington limit, $λ^* \sim 1$. The new model may be thought of as an opposite extreme to our earlier co-evolution scenario in Caplar et al. 2015. The main observable difference between the two co-evolution scenarios, presented here and in Caplar et al. 2015, is in the active fraction of low mass star-forming galaxies. We compare the predictions with the data from deep multi-wavelength surveys and find that the "quenching" scenario developed in the current paper is much to be preferred.
△ Less
Submitted 9 October, 2018;
originally announced October 2018.
-
A model for AGN variability on multiple timescales
Authors:
Lia F. Sartori,
Kevin Schawinski,
Benny Trakhtenbrot,
Neven Caplar,
Ezequiel Treister,
Michael J. Koss,
C. Megan Urry,
Ce Zhang
Abstract:
We present a framework to link and describe AGN variability on a wide range of timescales, from days to billions of years. In particular, we concentrate on the AGN variability features related to changes in black hole fuelling and accretion rate. In our framework, the variability features observed in different AGN at different timescales may be explained as realisations of the same underlying stat…
▽ More
We present a framework to link and describe AGN variability on a wide range of timescales, from days to billions of years. In particular, we concentrate on the AGN variability features related to changes in black hole fuelling and accretion rate. In our framework, the variability features observed in different AGN at different timescales may be explained as realisations of the same underlying statistical properties. In this context, we propose a model to simulate the evolution of AGN light curves with time based on the probability density function (PDF) and power spectral density (PSD) of the Eddington ratio ($L/L_{\rm Edd}$) distribution. Motivated by general galaxy population properties, we propose that the PDF may be inspired by the $L/L_{\rm Edd}$ distribution function (ERDF), and that a single (or limited number of) ERDF+PSD set may explain all observed variability features. After outlining the framework and the model, we compile a set of variability measurements in terms of structure function (SF) and magnitude difference. We then combine the variability measurements on a SF plot ranging from days to Gyr. The proposed framework enables constraints on the underlying PSD and the ability to link AGN variability on different timescales, therefore providing new insights into AGN variability and black hole growth phenomena.
△ Less
Submitted 13 March, 2018; v1 submitted 15 February, 2018;
originally announced February 2018.
-
Galaxy Zoo: major galaxy mergers are not a significant quenching pathway
Authors:
Anna K. Weigel,
Kevin Schawinski,
Neven Caplar,
Alfredo Carpineti,
Ross E. Hart,
Sugata Kaviraj,
William C. Keel,
Sandor J. Kruk,
Chris J. Lintott,
Robert C. Nichol,
Brooke D. Simmons,
Rebecca J. Smethurst
Abstract:
We use stellar mass functions to study the properties and the significance of quenching through major galaxy mergers. In addition to SDSS DR7 and Galaxy Zoo 1 data, we use samples of visually selected major galaxy mergers and post merger galaxies. We determine the stellar mass functions of the stages that we would expect major merger quenched galaxies to pass through on their way from the blue clo…
▽ More
We use stellar mass functions to study the properties and the significance of quenching through major galaxy mergers. In addition to SDSS DR7 and Galaxy Zoo 1 data, we use samples of visually selected major galaxy mergers and post merger galaxies. We determine the stellar mass functions of the stages that we would expect major merger quenched galaxies to pass through on their way from the blue cloud to the red sequence: 1: major merger, 2: post merger, 3: blue early type, 4: green early type and 5: red early type. Based on the similar mass function shapes we conclude that major mergers are likely to form an evolutionary sequence from star formation to quiescence via quenching. Relative to all blue galaxies, the major merger fraction increases as a function of stellar mass. Major merger quenching is inconsistent with the mass and environment quenching model. At z~0 major merger quenched galaxies are unlikely to constitute the majority of galaxies that transition the green valley. Furthermore, between z~0-0.5 major merger quenched galaxies account for 1-5% of all quenched galaxies at a given stellar mass. Major galaxy mergers are therefore not a significant quenching pathway, neither at z~0 nor within the last 5 Gyr. The majority of red galaxies must have been quenched through an alternative quenching mechanism which causes a slow blue to red evolution.
△ Less
Submitted 2 August, 2017;
originally announced August 2017.
-
AGN and their host galaxies in the local Universe: two mass independent Eddington ratio distribution functions characterize black hole growth
Authors:
Anna K. Weigel,
Kevin Schawinski,
Neven Caplar,
O. Ivy Wong,
Ezequiel Treister,
Benny Trakhtenbrot
Abstract:
We use a phenomenological model to show that black hole growth in the local Universe (z < 0.1) can be described by two separate, mass independent Eddington ratio distribution functions (ERDFs). We assume that black holes can be divided into two independent groups: those with radiatively efficient accretion, primarily hosted by optically blue and green galaxies, and those with radiatively inefficie…
▽ More
We use a phenomenological model to show that black hole growth in the local Universe (z < 0.1) can be described by two separate, mass independent Eddington ratio distribution functions (ERDFs). We assume that black holes can be divided into two independent groups: those with radiatively efficient accretion, primarily hosted by optically blue and green galaxies, and those with radiatively inefficient accretion, which are mainly found in red galaxies. With observed galaxy stellar mass functions as input, we show that the observed AGN luminosity functions can be reproduced by using mass independent, broken power law shaped ERDFs. We use the observed hard X-ray and 1.4 GHz radio luminosity functions to constrain the ERDF for radiatively efficient and inefficient AGN, respectively. We also test alternative ERDF shapes and mass dependent models. Our results are consistent with a mass independent AGN fraction and AGN hosts being randomly drawn from the galaxy population. We argue that the ERDF is not shaped by galaxy-scale effects, but by how efficiently material can be transported from the inner few parsecs to the accretion disc. Our results are incompatible with the simplest form of mass quenching where massive galaxies host higher accretion rate AGN. Furthermore, if reaching a certain Eddington ratio is a sufficient condition for maintenance mode, it can occur in all red galaxies, not just the most massive ones.
△ Less
Submitted 17 July, 2017;
originally announced July 2017.
-
Optical variability of AGN in the PTF/iPTF survey
Authors:
Neven Caplar,
Simon J. Lilly,
Benny Trakhtenbrot
Abstract:
We characterize the optical variability of quasars in the intermediate Palomar Transient Factory (iPTF) and Palomar Transient Factory (PTF) surveys. We re-calibrate the $r$-band light curves for $\sim$28,000 luminous, broad-line AGNs from the SDSS, producing a total of $\sim$2.4 million photometric data points. We utilize both the structure function (SF) and power spectrum density (PSD) formalisms…
▽ More
We characterize the optical variability of quasars in the intermediate Palomar Transient Factory (iPTF) and Palomar Transient Factory (PTF) surveys. We re-calibrate the $r$-band light curves for $\sim$28,000 luminous, broad-line AGNs from the SDSS, producing a total of $\sim$2.4 million photometric data points. We utilize both the structure function (SF) and power spectrum density (PSD) formalisms to search for links between the optical variability and the physical parameters of the accreting supermassive black holes that power the quasars. The excess variance (SF$^{2}$) of the quasar sample tends to zero at very short time separations, validating our re-calibration of the time-series data. We find that the the amplitude of variability at a given time-interval, or equivalently the time-scale of variability to reach a certain amplitude, is most strongly correlated with luminosity with weak or no dependence on black hole mass and redshift. For a variability level of SF($τ$)=0.07 mag, the time-scale has a dependency of $τ\propto L^{0.4}$. This is broadly consistent with the expectation from a simple Keplerian accretion disk model, which provides $τ\propto L^{0.5}$. The PSD analysis also reveals that many quasar light curves are steeper than a damped random walk. We find a correlation between the steepness of the PSD slopes, specifically the fraction of slopes steeper than 2.5, and black hole mass, although we cannot exclude the possibility that luminosity or Eddington ratio are the drivers of this effect. This effect is also seen in the SF analysis of the (i)PTF data, and in a PSD analysis of quasars in the SDSS Stripe 82.
△ Less
Submitted 9 November, 2016;
originally announced November 2016.
-
Quantitative Evaluation of Gender Bias in Astronomical Publications from Citation Counts
Authors:
Neven Caplar,
Sandro Tacchella,
Simon Birrer
Abstract:
We analyze the role of first (leading) author gender on the number of citations that a paper receives, on the publishing frequency and on the self-citing tendency. We consider a complete sample of over 200,000 publications from 1950 to 2015 from five major astronomy journals. We determine the gender of the first author for over 70% of all publications. The fraction of papers which have a female fi…
▽ More
We analyze the role of first (leading) author gender on the number of citations that a paper receives, on the publishing frequency and on the self-citing tendency. We consider a complete sample of over 200,000 publications from 1950 to 2015 from five major astronomy journals. We determine the gender of the first author for over 70% of all publications. The fraction of papers which have a female first author has increased from less than 5% in the 1960s to about 25% today. We find that the increase of the fraction of papers authored by females is slowest in the most prestigious journals such as Science and Nature. Furthermore, female authors write 19$\pm$7% fewer papers in seven years following their first paper than their male colleagues. At all times papers with male first authors receive more citations than papers with female first authors. This difference has been decreasing with time and amounts to $\sim$6% measured over the last 30 years. To account for the fact that the properties of female and male first author papers differ intrinsically, we use a random forest algorithm to control for the non-gender specific properties of these papers which include seniority of the first author, number of references, total number of authors, year of publication, publication journal, field of study and region of the first author's institution. We show that papers authored by females receive 10.4$\pm$0.9% fewer citations than what would be expected if the papers with the same non-gender specific properties were written by the male authors. Finally, we also find that female authors in our sample tend to self-cite more, but that this effect disappears when controlled for non-gender specific variables.
△ Less
Submitted 27 October, 2016;
originally announced October 2016.
-
AGN evolution from a galaxy evolution viewpoint
Authors:
Neven Caplar,
Simon J. Lilly,
Benny Trakhtenbrot
Abstract:
We explore the connections between the evolving galaxy and AGN populations. We present a simple phenomenological model that links the evolving galaxy mass function and the evolving quasar luminosity function, which makes specific and testable predictions for the distribution of host galaxy masses for AGN of different luminosities. We show that the $φ^{*}$ normalisations of the galaxy mass function…
▽ More
We explore the connections between the evolving galaxy and AGN populations. We present a simple phenomenological model that links the evolving galaxy mass function and the evolving quasar luminosity function, which makes specific and testable predictions for the distribution of host galaxy masses for AGN of different luminosities. We show that the $φ^{*}$ normalisations of the galaxy mass function and of the AGN luminosity function closely track each other over a wide range of redshifts, implying a constant "duty cycle" of AGN activity. The strong redshift evolution in the AGN $L^*$ can be produced by either an evolution in the distribution of Eddington ratios, or in the $m_{bh}/m_{*}$ mass ratio, or both. To try to break this degeneracy we look at the distribution of AGN in the SDSS ($m_{bh},L$) plane, showing that an evolving ratio $m_{bh}/m_{*} \propto (1+z)^2$ reproduces the observed data and also reproduces the local relations which connect the black hole population with the host galaxies for both quenched and star-forming populations. We stress that observational studies that compare the masses of black holes in active galaxies at high redshift with those in quiescent galaxies locally will always see much weaker evolution. Evolution of this form would produce, or could be produced by, a redshift-independent $m_{bh} - σ$ relation and could explain why the local $m_{bh} - σ$ relation is tighter than $m_{bh} - m_{*}$ even if $σ$ is not directly linked to black hole growth. Irrespective of the evolution of $m_{bh}/m_{*}$, the model reproduces both the appearance of "downsizing" and the so-called "sub-Eddington boundary" without any mass-dependence in the evolution of black hole growth rates.
△ Less
Submitted 31 August, 2015; v1 submitted 13 November, 2014;
originally announced November 2014.
-
Generalized models of unification of dark matter and dark energy
Authors:
Neven Caplar,
Hrvoje Stefancic
Abstract:
A model of unification of dark matter and dark energy based on the modeling of the speed of sound as a function of the parameter of the equation of state is introduced. It is found that the model in which the speed of sound depends on the power of the parameter of the equation of state, $c_s^2=α(-w)^γ$, contains the generalized Chaplygin gas models as its subclass. An effective scalar field descri…
▽ More
A model of unification of dark matter and dark energy based on the modeling of the speed of sound as a function of the parameter of the equation of state is introduced. It is found that the model in which the speed of sound depends on the power of the parameter of the equation of state, $c_s^2=α(-w)^γ$, contains the generalized Chaplygin gas models as its subclass. An effective scalar field description of the model is obtained in a parametric form which in some cases can be translated into a closed form solution for the scalar field potential. A constraint on model parameters is obtained using the observational data on the Hubble parameter at different redshifts.
△ Less
Submitted 3 April, 2013; v1 submitted 2 August, 2012;
originally announced August 2012.