Researcher profile

Luke Zoltan Kelley

Luke Zoltan Kelley contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

preprint2022arXiv

Probing the $z\gtrsim6$ quasars in a universe with IllustrisTNG physics: Impact of gas-based black hole seeding models

We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z\gtrsim6$ quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in IllustrisTNG. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim10^{12}M_{\odot}/h$ halos by $z=7$) using a $(9~\mathrm{Mpc}/h)^3$ simulated volume. BH growth is maximal within halos that are compact and have a low tidal field. For these halos, we consider an array of gas-based seeding prescriptions wherein $M_{\mathrm{seed}}=10^4-10^6~M_{\odot}/h$ seeds are inserted in halos above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $\tilde{M}_{\mathrm{h}}$ and $\tilde{M}_{\mathrm{sf,mp}}$, respectively, in units of $M_{\mathrm{seed}}$). We find that a seed model with $\tilde{M}_{\mathrm{sf,mp}}=5$ and $\tilde{M}_{\mathrm{h}}=3000$ successfully produces a $z\sim6$ quasar with $\sim10^9~M_{\odot}$ mass and $\sim10^{47}~\mathrm{ergs~s^ {-1}}$ luminosity. BH mergers play a crucial role at $z\gtrsim9$, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. When more stringent seeding conditions are applied (for e.g., $\tilde{M}_{\mathrm{sf,mp}}=1000$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $z\gtrsim6$ quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors $\gtrsim10$) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results show that progenitors of $z\sim6$ quasars have distinct BH merger histories for different seeding models, which will be distinguishable with LISA observations.

preprint2022arXiv

Running Late: Testing Delayed Supermassive Black Hole Growth Models Against the Quasar Luminosity Function

Observations of massive galaxies at low redshift have revealed approximately linear scaling relations between the mass of a supermassive black hole (SMBH) and properties of its host galaxy. How these scaling relations evolve with redshift and whether they extend to lower-mass galaxies however remain open questions. Recent galaxy formation simulations predict a delayed, or "two-phase", growth of SMBHs: slow, highly intermittent BH growth due to repeated gas ejection by stellar feedback in low-mass galaxies, followed by more sustained gas accretion that eventually brings BHs onto the local scaling relations. The predicted two-phase growth implies a steep increase, or "kink", in BH-galaxy scaling relations at a stellar mass $M_{*}\sim5\times10^{10} M_{\odot}$. We develop a parametric, semi-analytic model to compare different SMBH growth models against observations of the quasar luminosity function (QLF) at $z\sim0.5-4$. We compare models in which the relation between SMBH mass and galaxy mass is purely linear versus two-phase models. The models are anchored to the observed galaxy stellar mass function, and the BH mass functions at different redshifts are consistently connected by the accretion rates contributing to the QLF. The best fits suggest that two-phase evolution is significantly preferred by the QLF data over a purely linear scaling relation. Moreover, when the model parameters are left free, the two-phase model fits imply a transition mass consistent with that predicted by simulations. Our analysis motivates further observational tests, including measurements of BH masses and AGN activity at the low-mass end, which could more directly test two-phase SMBH growth.

preprint2022arXiv

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.

preprint2021arXiv

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature.

preprint2020arXiv

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

preprint2020arXiv

The effect of differential accretion on the Gravitational Wave Background and the present day MBH Binary population

Massive black hole binaries (MBHBs) form as a consequence of galaxy mergers. However, it is still unclear whether they typically merge within a Hubble time, and how accretion may affect their evolution. These questions will be addressed by pulsar timing arrays (PTAs), which aim to detect the GW background (GWB) emitted by MBHBs during the last Myrs of inspiral. Here we investigate the influence of differential accretion on MBHB merger rates, chirp masses and the resulting GWB spectrum. We evolve a MBHB sample from the Illustris hydrodynamic cosmological simulation using semi-analytic models and for the first time self-consistently evolve their masses with binary accretion models. In all models, MBHBs coalesce with median total masses up to $1.5 \times 10^8 M_{\odot}$, up to $3-4$ times larger than in models neglecting accretion. In our model with the largest plausible impact, the median mass ratio of coalescing MBHBs increases by a factor $3.6$, the coalescence rate by $52.3\%$, and the GWB amplitude by a factor $4.0$, yielding a dimensionless GWB strain $A_{yr^{-1}} = 1 \times 10^{-15}$. Our model that favours accretion onto the primary MBH reduces the median mass ratio of coalescing MBHBs by a factor of $2.9$, and yields a GWB amplitude $A_{yr^{-1}} = 3.1 \times 10^{-16}$. This is nearly indistinguishable from our model neglecting accretion, despite higher MBHB masses at coalescence. \textbf{We further predict binary separation and mass ratio distributions of stalled MBHBs in the low-redshift universe, and find that these depend sensitively on binary accretion models. This presents the potential for combined EM and GW observational constraints on merger rates and accretion models of MBHB populations.}

preprint2019arXiv

Probing Massive Black Hole Binary Populations with LISA

ESA and NASA are moving forward with plans to launch LISA around 2034. With data from the Illustris cosmological simulation, we provide analysis of LISA detection rates accompanied by characterization of the merging massive black hole population. Massive black holes of total mass $\sim10^5-10^{10} M_\odot$ are the focus of this study. We evolve Illustris massive black hole mergers, which form at separations on the order of the simulation resolution ($\sim$kpc scales), through coalescence with two different treatments for the binary massive black hole evolutionary process. The coalescence times of the population, as well as physical properties of the black holes, form a statistical basis for each evolutionary treatment. From these bases, we Monte Carlo synthesize many realizations of the merging massive black hole population to build mock LISA detection catalogs. We analyze how our massive black hole binary evolutionary models affect detection rates and the associated parameter distributions measured by LISA. With our models, we find massive black hole binary detection rates with LISA of $\sim0.5-1$ yr$^{-1}$ for massive black holes with masses greater than $10^5M_\odot$. This should be treated as a lower limit primarily because our massive black hole sample does not include masses below $10^5M_\odot$, which may significantly add to the observed rate. We suggest reasons why we predict lower detection rates compared to much of the literature.