ABSTRACT
Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed ...quasars with a median redshift of 1.8 in a 4.6 deg2 overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-yr long multicolour (griz) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to r ∼23.5 mag) powered by less massive black holes (with masses ≳ 108.5M⊙) at high redshift for the first time. We find five candidates with significant (at >99.74 per cent single-frequency significance in at least two bands with a global p-value of ∼7 × 10−4–3 × 10−3 accounting for the look-elsewhere effect) periodicity with observed periods of ∼3–5 yr (i.e. 1–2 yr in rest frame) having ∼4–6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of ${\sim }0.8^{+0.5}_{-0.3}$ per cent or ${\sim }1.1^{+0.7}_{-0.5}$ quasar per deg2. Our detection rate is 4–80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. We discuss implications on the future direct detection of low-frequency gravitational waves. Continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins.
Abstract
We quantify evolution in the cluster-scale stellar mass–halo mass (SMHM) relation’s parameters using 2323 clusters and brightest central galaxies (BCGs) over the redshift range 0.03 ≤
z
≤ ...0.60. The precision on the inferred SMHM parameters is improved by including the magnitude gap (
m
gap
) between the BCG and fourth-brightest cluster member (M14) as a third parameter in the SMHM relation. At fixed halo mass, accounting for
m
gap
, through a stretch parameter, reduces the SMHM relation’s intrinsic scatter. To explore this redshift range, we use clusters, BCGs, and cluster members identified using the Sloan Digital Sky Survey C4 and redMaPPer cluster catalogs and the Dark Energy Survey redMaPPer catalog. Through this joint analysis, we detect no systematic differences in BCG stellar mass,
m
gap
, and cluster mass (inferred from richness) between the data sets. We utilize the Pareto function to quantify each parameter’s evolution. We confirm prior findings of negative evolution in the SMHM relation’s slope (3.5
σ
), and detect negative evolution in the stretch parameter (4.0
σ
) and positive evolution in the offset parameter (5.8
σ
). This observed evolution, combined with the absence of BCG growth, when stellar mass is measured within 50 kpc, suggests that this evolution results from changes in the cluster’s
m
gap
. For this to occur, late-term growth must be in the intracluster light surrounding the BCG. We also compare the observed results to IllustrisTNG 300-1 cosmological hydrodynamic simulations and find modest qualitative agreement. However, the simulations lack the evolutionary features detected in the real data.
ABSTRACT
Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational ...wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here, we report discovery of a periodicity (p = 1607 ± 7 d) at 99.95 per cent significance (with a global p value of ∼10−3 accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67−002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ∼20-yr time baseline spans ∼4.6 cycles of the observed 4.4-yr (rest frame 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion discs. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disc. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302−102. Furthermore, the frequency dependence of the variability amplitudes disfavours Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna.
The size of the dust torus in active galactic nuclei (AGNs) and their high-luminosity counterparts, quasars, can be inferred from the time delay between UV/optical accretion disk continuum ...variability and the response in the mid-infrared (MIR) torus emission. This dust reverberation mapping (RM) technique has been successfully applied to ∼70 z 0.3 AGNs and quasars. Here we present first results of our dust RM program for distant quasars covered in the Sloan Digital Sky Survey Stripe 82 region combining ∼20 yr ground-based optical light curves with 10 yr MIR light curves from the WISE satellite. We measure a high-fidelity lag between W1 band (3.4 m) and g band for 587 quasars over 0.3 z 2 ( ) and two orders of magnitude in quasar luminosity. They tightly follow (intrinsic scatter ∼0.17 dex in lag) the IR lag-luminosity relation observed for z < 0.3 AGNs, revealing a remarkable size-luminosity relation for the dust torus over more than four decades in AGN luminosity, with little dependence on additional quasar properties such as Eddington ratio and variability amplitude. This study motivates further investigations in the utility of dust RM for cosmology and strongly endorses a compelling science case for the combined 10 yr Vera C. Rubin Observatory Legacy Survey of Space and Time (optical) and 5 yr Nancy Grace Roman Space Telescope 2 m light curves in a deep survey for low-redshift AGN dust RM with much lower luminosities and shorter, measurable IR lags. The compiled optical and MIR light curves for 7384 quasars in our parent sample are made public with this work.
We detect the kinematic Sunyaev-Zel'dovich (kSZ) effect with a statistical significance of 4.2σ by combining a cluster catalogue derived from the first year data of the Dark Energy Survey with cosmic ...microwave background temperature maps from the South Pole Telescope Sunyaev-Zel'dovich Survey. This measurement is performed with a differential statistic that isolates the pairwise kSZ signal, providing the first detection of the large-scale, pairwise motion of clusters using redshifts derived from photometric data. By fitting the pairwise kSZ signal to a theoretical template, we measure the average central optical depth of the cluster sample,
$\bar{\tau }_e = (3.75 \pm 0.89)\times 10^{-3}$
. We compare the extracted signal to realistic simulations and find good agreement with respect to the signal to noise, the constraint on
$\bar{\tau }_e$
, and the corresponding gas fraction. High-precision measurements of the pairwise kSZ signal with future data will be able to place constraints on the baryonic physics of galaxy clusters, and could be used to probe gravity on scales ≳100 Mpc.
Exploiting the sensitivity of the IRAM NOrthern Extended Millimeter Array (NOEMA) and its ability to process large instantaneous bandwidths, we have studied the morphology and other properties of the ...molecular gas and dust in the star forming galaxy, H-ATLAS J131611.5+281219 (HerBS-89a), at
z
= 2.95. High angular resolution (0
.
″3) images reveal a partial 1
.
″0 diameter Einstein ring in the dust continuum emission and the molecular emission lines of
12
CO(9−8) and H
2
O(2
02
− 1
11
). Together with lower angular resolution (0
.
″6) images, we report the detection of a series of molecular lines including the three fundamental transitions of the molecular ion OH
+
, namely (1
1
− 0
1
), (1
2
− 0
1
), and (1
0
− 0
1
), seen in absorption; the molecular ion CH
+
(1 − 0) seen in absorption, and tentatively in emission; two transitions of amidogen (NH
2
), namely (2
02
− 1
11
) and (2
20
− 2
11
) seen in emission; and HCN(11 − 10) and/or NH(1
2
− 0
1
) seen in absorption. The NOEMA data are complemented with Very Large Array data tracing the
12
CO(1 − 0) emission line, which provides a measurement of the total mass of molecular gas and an anchor for a CO excitation analysis. In addition, we present
Hubble
Space Telescope imaging that reveals the foreground lensing galaxy in the near-infrared (1.15
μ
m). Together with photometric data from the Gran Telescopio Canarias, we derive a photometric redshift of
z
phot
= 0.9
−0.5
+0.3
for the foreground lensing galaxy. Modeling the lensing of HerBS-89a, we reconstruct the dust continuum (magnified by a factor
μ
≃ 5.0) and molecular emission lines (magnified by
μ
∼ 4 − 5) in the source plane, which probe scales of ∼0
.
″1 (or 800 pc). The
12
CO(9 − 8) and H
2
O(2
02
− 1
11
) emission lines have comparable spatial and kinematic distributions; the source-plane reconstructions do not clearly distinguish between a one-component and a two-component scenario, but the latter, which reveals two compact rotating components with sizes of ≈1 kpc that are likely merging, more naturally accounts for the broad line widths observed in HerBS-89a. In the core of HerBS-89a, very dense gas with
n
H
2
∼ 10
7 − 9
cm
−3
is revealed by the NH
2
emission lines and the possible HCN(11 − 10) absorption line. HerBS-89a is a powerful star forming galaxy with a molecular gas mass of
M
mol
= (2.1 ± 0.4) × 10
11
M
⊙
, an infrared luminosity of
L
IR
= (4.6 ± 0.4) × 10
12
L
⊙
, and a dust mass of
M
dust
= (2.6 ± 0.2) × 10
9
M
⊙
, yielding a dust-to-gas ratio
δ
GDR
≈ 80. We derive a star formation rate SFR = 614 ± 59
M
⊙
yr
−1
and a depletion timescale
τ
depl
= (3.4 ± 1.0) × 10
8
years. The OH
+
and CH
+
absorption lines, which trace low (∼100 cm
−3
) density molecular gas, all have their main velocity component red-shifted by Δ
V
∼ 100 km s
−1
relative to the global CO reservoir. We argue that these absorption lines trace a rare example of gas inflow toward the center of a galaxy, indicating that HerBS-89a is accreting gas from its surroundings.
ABSTRACT
Local primordial non-Gaussianity (PNG) is a promising observable of the underlying physics of inflation, characterized by $f_{\rm NL}^{\rm loc}$. We present the methodology to measure ...$f_{\rm NL}^{\rm loc}$ from the Dark Energy Survey (DES) data using the two-point angular correlation function (ACF) with scale-dependent bias. One of the focuses of the work is the integral constraint (IC). This condition appears when estimating the mean number density of galaxies from the data and is key in obtaining unbiased $f_{\rm NL}^{\rm loc}$ constraints. The methods are analysed for two types of simulations: ∼246 goliat-png N-body small area simulations with fNL equal to −100 and 100, and 1952 Gaussian ice-cola mocks with fNL = 0 that follow the DES angular and redshift distribution. We use the ensemble of goliat-png mocks to show the importance of the IC when measuring PNG, where we recover the fiducial values of fNL within the 1σ when including the IC. In contrast, we found a bias of ΔfNL ∼ 100 when not including it. For a DES-like scenario, we forecast a bias of ΔfNL ∼ 23, equivalent to 1.8σ, when not using the IC for a fiducial value of fNL = 100. We use the ice-cola mocks to validate our analysis in a realistic DES-like set-up finding it robust to different analysis choices: best-fitting estimator, the effect of IC, BAO damping, covariance, and scale choices. We forecast a measurement of fNL within σ(fNL) = 31 when using the DES-Y3 BAO sample, with the ACF in the 1 deg < θ < 20 deg range.
ABSTRACT
We present a sample of 706, z < 1.5 active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an ...area of 4.64 deg2. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 d, a 6-yr baseline, and a single-epoch imaging depth of up to g ∼ 24.5. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast }\lt 10^{9.5}\, {\rm M}_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability time-scales (∼ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole (BH) masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalogue, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower BH masses in deep fields, which may be more ‘pristine’ analogues of supermassive BH seeds.
Abstract
Current and future cosmological analyses with Type Ia supernovae (SNe Ia) face three critical challenges: (i) measuring the redshifts from the SNe or their host galaxies; (ii) classifying ...the SNe without spectra; and (iii) accounting for correlations between the properties of SNe Ia and their host galaxies. We present here a novel approach that addresses each of these challenges. In the context of the Dark Energy Survey (DES), we analyze an SN Ia sample with host galaxies in the redMaGiC galaxy catalog, a selection of luminous red galaxies. redMaGiC photo-
z
estimates are expected to be accurate to
σ
Δ
z
/(1+
z
)
∼ 0.02. The DES-5YR photometrically classified SN Ia sample contains approximately 1600 SNe, and 125 of these SNe are in redMaGiC galaxies. We demonstrate that redMaGiC galaxies almost exclusively host SNe Ia, reducing concerns relating to classification uncertainties. With this subsample, we find similar Hubble scatter (to within ∼0.01 mag) using photometric redshifts in place of spectroscopic redshifts. With detailed simulations, we show that the bias due to using redMaGiC photo-
z
s on the measurement of the dark energy equation of state
w
is up to Δ
w
∼ 0.01–0.02. With real data, we measure a difference in
w
when using the redMaGiC photo-
z
s versus the spec-
z
s of Δ
w
= 0.005. Finally, we discuss how SNe in redMaGiC galaxies appear to comprise a more standardizable population, due to a weaker relation between color and luminosity (
β
) compared to the DES-3YR population by ∼5
σ
. These results establish the feasibility of performing redMaGiC SN cosmology with photometric survey data in the absence of spectroscopic data.
ABSTRACT
Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback ...strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 × 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 × 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2σ.