There is broad consensus that Type Ia supernovae (SNe Ia) are the thermonuclear explosions of C/O white dwarfs (WDs) in binary systems, but their progenitors and explosion mechanisms remain uncertain. Their earliest light curves probe the outermost ejecta and provide important constraints on the explosion. We analyse the ZTF DR2 sample of SNe Ia to search for objects exhibiting early flux excess (EEx SNe Ia). We employ two complementary approaches: (i) identifying deviations from models without early excess using power-law and SALT2 fits, and (ii) comparing observations with double-detonation (DD) and companion-interaction (CI) models. The latter identify 145 and 199 candidates, respectively, although the methods disagree substantially on both the objects selected and the total number of candidates. Combining all methods, we identify 17 robust EEx SN Ia candidates, including six over-luminous events (five 91T-like and one 03fg-like). The strongest excesses reach ∼10–17% of the peak flux. The best candidates have, on average, higher stretch and preferentially occur in lower-mass, bluer host galaxies, suggesting younger stellar populations. We estimate that EEx SNe Ia comprise ≲25% of all SNe Ia, consistent with previous work, with over-luminous subtypes showing a higher relative incidence than normal SNe Ia. Among the DD models, the preferred solutions favour intermediate WD core masses (1 M_⊙), low shell-burning fractions (20%), and heavier He-burning products (^56Ni), whereas the CI models show no clear parameter trends. Overall, we find no compelling evidence that either scenario is preferred as the origin of the observed early flux excesses.

With the 4-metre Multi-Object Spectroscopic Telescope (4MOST) expected to provide an influx of transient spectra when it begins observations in early 2026 we consider the potential for real-time classification of these spectra. We investigate three extant spectroscopic transient classifiers: the Deep Automated Supernova and Host classifier (dash), Next Generation SuperFit (ngsf), and SuperNova IDentification (snid), with a focus on comparing the completeness and purity of the transient samples they produce. We manually simulate fibre losses critical for accurately determining host contamination and use the 4MOST Exposure Time Calculator to produce realistic, 4MOST-like, host-galaxy contaminated spectra. We investigate the three classifiers individually and in all possible combinations. We find that a combination of dash and ngsf can produce a supernova (SN) Ia sample with a purity of 99.9 per cent, while successfully classifying 70 per cent of SNe Ia. However, it struggles to classify non-SN Ia transients. We investigate photometric cuts to transient magnitude and the transient's fraction of total fibre flux, finding that both can be used to improve non-SN Ia transient classification completeness by 8-44 per cent with SNe Ibc benefitting the most and superluminous (SL) SNe the least. Finally, we present an example classification plan for live classification and the predicted purities and completeness across five transient classes: Ia, Ibc, II, SL, and non-SN transients. We find that it is possible to classify 75 per cent of input spectra with >70 per cent purity in all classes except non-SN transients. Precise values can be varied using different classifiers and photometric cuts to suit the needs of a given study.

We present the La Silla Schmidt Southern Survey (LS4), a new wide-field, time-domain survey to be conducted with the 1 m ESO Schmidt telescope. The 268 megapixel LS4 camera mosaics 32 2k x 4k fully depleted CCDs, providing a similar to 20 deg2 field of view with 1 '' pixel-1 resolution. The LS4 camera will have excellent performance at longer wavelengths: in a standard 45 s exposure the expected 5 sigma limiting magnitudes in g, i, z are similar to 21.5, similar to 20.9, and similar to 20.3 mag (AB), respectively. The telescope design requires a novel filter holder that fixes different bandpasses over each quadrant of the detector. Two quadrants will have i band, while the other two will be g and z band with color information obtained by dithering targets across the different quadrants. The majority (90%) of the observing time will be used to conduct a public survey that monitors the extragalactic sky at both moderate (3 days) and high (1 day) cadence, as well as focused observations within the Galactic plane and bulge. Alerts from the public survey will be broadcast to the community via established alert brokers. LS4 will run concurrently with the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). The combination of LS4+LSST will enable detailed holistic monitoring of many nearby transients: high-cadence LS4 observations will resolve the initial rise and peak of the light curve while less-frequent but deeper observations by LSST will characterize the years before and after explosion. Here, we summarize the primary science objectives of LS4 including microlensing events in the Galaxy, extragalactic transients powered by massive black holes or stellar explosions, the search for electromagnetic counterparts to multi-messenger events, and supernova cosmology.

Dust associated with various stellar sources in galaxies at all cosmic epochs remains a controversial topic, particularly whether supernovae (SNe) play an important role in dust production. We report evidence of dust formation in the cold, dense shell behind the ejecta-circumstellar medium (CSM) interaction in the Type Ia-CSM SN 2018evt three years after the explosion, characterized by a rise in the mid-infrared (MIR) emission accompanied by an accelerated decline in the optical radiation of the SN. Such a dust-formation picture is also corroborated by the concurrent evolution of the profiles of the Ha emission line. Our model suggests enhanced CSM dust concentration at increasing distances from the SN as compared to what can be expected from the density profile of the mass loss from a steady stellar wind. By the time of the last MIR observations at day +1041, a total amount of 1.2+-0.2x10^{-2} Msun of new dust has been formed by SN 2018evt, making SN 2018evt one of the most prolific dust factories among SNe with evidence of dust formation. The unprecedented witness of the intense production procedure of dust may shed light on the perceptions of dust formation in cosmic history.

Low-density cosmic voids gravitationally lens the cosmic microwave background (CMB), leaving a negative imprint on the CMB convergence kappa. This effect provides insight into the distribution of matter within voids, and can also be used to study the growth of structure. We measure this lensing imprint by cross-correlating the Planck CMB lensing convergence map with voids identified in the Dark Energy Survey Year 3 (DES Y3) data set, covering approximately 4200 deg(2) of the sky. We use two distinct void-finding algorithms: a 2D void-finder that operates on the projected galaxy density field in thin redshift shells, and a new code, Voxel, which operates on the full 3D map of galaxy positions. We employ an optimal matched filtering method for cross-correlation, using the Marenostrum Institut de Ciencies ` de l'Espai N-body simulation both to establish the template for the matched filter and to calibrate detection significances. Using the DES Y3 photometric luminous red galaxy sample, we measure A(kappa) , the amplitude of the observed lensing signal relative to the simulation template, obtaining A(kappa) = 1.03 +/- 0.22 (4.6 sigma significance) for Voxel and A(kappa) = 1.02 +/- 0.17 (5.9 sigma significance) for 2D voids, both consistent with Lambda cold dark matter expectations. We additionally invert the 2D void-finding process to identify superclusters in the projected density field, for which we measure A(kappa) = 0.87 +/- 0.15 (5.9 sigma significance). The leading source of noise in our measurements is Planck noise, implying that data from the Atacama Cosmology Telescope, South Pole Telescope and CMB-S4 will increase sensitivity and allow for more precise measurements.

We present the first homogeneous release of several thousand Type Ia supernovae (SNe Ia), all having spectroscopic classification, and spectroscopic redshifts for half the sample. This release, named the "DR2", contains 3628 nearby (z < 0.3) SNe Ia discovered, followed and classified by the Zwicky Transient Facility survey between March 2018 and December 2020. Of these, 3000 have good-to-excellent sampling and 2667 pass standard cosmology light-curve quality cuts. This release is thus the largest SN Ia release to date, increasing by an order of magnitude the number of well characterized low-redshift objects. With the "DR2", we also provide a volume-limited (z < 0.06) sample of nearly a thousand SNe Ia. With such a large, homogeneous and well controlled dataset, we are studying key current questions on SN cosmology, such as the linearity SNe Ia standardization, the SN and host dependencies, the diversity of the SN Ia population, and the accuracy of the current light-curve modeling. These, and more, are studied in detail in a series of articles associated with this release. Alongside the SN Ia parameters, we publish our force-photometry gri-band light curves, 5138 spectra, local and global host properties, observing logs, and a python tool to ease use and access of these data. The photometric accuracy of the "DR2" is not yet suited for cosmological parameter inference, which will follow as "DR2.5" release. We nonetheless demonstrate that the multi-thousand SN Ia Hubble Diagram has a typical 0.15 mag scatter.

Measurements of the growth rate of structures at $z < 0.1$ with peculiar velocity surveys have the potential of testing the validity of general relativity on cosmic scales. In this work, we present growth-rate measurements from realistic simulated sets of type-Ia supernovae (SNe Ia) from the Zwicky Transient Facility (ZTF). We describe our simulation methodology, the light-curve fitting and peculiar velocity estimation. Using the maximum likelihood method, we derive constraints on $f\sigma_8$ using only ZTF SN Ia peculiar velocities. We carefully tested the method and we quantified biases due to selection effects (photometric detection, spectroscopic follow-up for typing) on several independent realizations. We simulated the equivalent of 6 years of ZTF data, and considering an unbiased spectroscopically typed sample at $z < 0.06$, we obtained unbiased estimates of $f\sigma_8$ with an average uncertainty of 19% precision. We also investigated the information gain in applying bias correction methods. Our results validate our framework which can be used on real ZTF data.

BACKGROUND:We aimed to develop a machine learning algorithm to predict the presence of a culprit lesion in patients with out-of-hospital cardiac arrest (OHCA).METHODS:We used the King's Out-of-Hospital Cardiac Arrest Registry, a retrospective cohort of 398 patients admitted to King's College Hospital between May 2012 and December 2017. The primary outcome was the presence of a culprit coronary artery lesion, for which a gradient boosting model was optimized to predict. The algorithm was then validated in two independent European cohorts comprising 568 patients.RESULTS:A culprit lesion was observed in 209/309 (67.4%) patients receiving early coronary angiography in the development, and 199/293 (67.9%) in the Ljubljana and 102/132 (61.1%) in the Bristol validation cohorts, respectively. The algorithm, which is presented as a web application, incorporates nine variables including age, a localizing feature on electrocardiogram (ECG) (≥2 mm of ST change in contiguous leads), regional wall motion abnormality, history of vascular disease and initial shockable rhythm. This model had an area under the curve (AUC) of 0.89 in the development and 0.83/0.81 in the validation cohorts with good calibration and outperforms the current gold standard-ECG alone (AUC: 0.69/0.67/0/67).CONCLUSIONS:A novel simple machine learning-derived algorithm can be applied to patients with OHCA, to predict a culprit coronary artery disease lesion with high accuracy.

Type Ia supernovae (SNe Ia) are thermonuclear explosions of degenerate white dwarf stars destabilized by mass accretion from a companion star 1 , but the nature of their progenitors remains poorly understood. A way to discriminate between progenitor systems is through radio observations; a non-degenerate companion star is expected to lose material through winds 2 or binary interaction 3 before explosion, and the supernova ejecta crashing into this nearby circumstellar material should result in radio synchrotron emission. However, despite extensive efforts, no type Ia supernova (SN Ia) has ever been detected at radio wavelengths, which suggests a clean environment and a companion star that is itself a degenerate white dwarf star 4 , 5 . Here we report on the study of SN 2020eyj, a SN Ia showing helium-rich circumstellar material, as demonstrated by its spectral features, infrared emission and, for the first time in a SN Ia to our knowledge, a radio counterpart. On the basis of our modelling, we conclude that the circumstellar material probably originates from a single-degenerate binary system in which a white dwarf accretes material from a helium donor star, an often proposed formation channel for SNe Ia (refs. 6 , 7 ). We describe how comprehensive radio follow-up of SN 2020eyj-like SNe Ia can improve the constraints on their progenitor systems.

In the local universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (H i ), radio surveys to probe the cosmic evolution of H i in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep H i survey, we report the first untargeted detection of an OHM at z > 0.5, LADUMA J033046.20−275518.1 (nicknamed “Nkalakatha”). The host system, WISEA J033046.26−275518.3, is an infrared-luminous radio galaxy whose optical redshift z ≈ 0.52 confirms the MeerKAT emission-line detection as OH at a redshift z OH = 0.5225 ± 0.0001 rather than H i at lower redshift. The detected spectral line has 18.4 σ peak significance, a width of 459 ± 59 km s −1 , and an integrated luminosity of (6.31 ± 0.18 [statistical] ± 0.31 [systematic]) × 10 3 L ⊙ , placing it among the most luminous OHMs known. The galaxy’s far-infrared luminosity L FIR = (1.576 ±0.013) × 10 12 L ⊙ marks it as an ultraluminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step toward a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts.

Erik C. Kool, Joel Johansson, Jesper Sollerman, Javier Moldón, Takashi J. Moriya, Steve Schulze, Laura Chomiuk, Miguel Pérez-Torres, Chelsea Harris, Seppo Mattila, Peter Lundqvist, Matthew Graham, Sheng Yang, Daniel A. Perley, Nora Linn Strotjohann, Christoffer Fremling, Avishay Gal-Yam, Jeremy Lezmy, Kate Maguire, Conor Omand, Mathew Smith, Igor Andreoni, Eric C. Bellm, Joshua S. Bloom, Kishalay De, Steven L. Groom, Mansi M. Kasliwal, Frank J. Masci, Michael S. Medford, Sungmin Park, Josiah Purdum, Thomas M. Reynolds, Reed Riddle, Estelle Robert, Stuart D. Ryder, Yashvi Sharma, Daniel Stern [affiliations can be found after the references] Type Ia supernovae are thought to be carbon-oxygen white dwarf stars that explode after accreting material from a companion star, but despite extensive studies the nature of the companion star is still poorly understood, as is the explosion mechanism. In the single degenerate scenario, the companion is a non-degenerate star that loses material through winds and/or binary interaction, and a few Type Ia supernovae have shown evidence for hydrogen-rich circumstellar material. We present here the study of SN 2020eyj, a unique Type Ia supernova showing delayed interaction with helium-rich, but hydrogen-poor, circumstellar material. This material surrounding SN 2020eyj is revealed by its unusual light curve and infrared emission, narrow helium emission lines and, for the first time ever in a Type Ia supernova, also a radio counterpart. The circumstellar material likely originates from the companion star, providing the first direct evidence for a, so far hypothesized, single degenerate progenitor system composed of a white dwarf and a helium donor star.

The discovery of the electromagnetic counterpart to the binary neutron star (NS) merger GW170817 has opened the era of gravitational-wave multimessenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multiwavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of NS mergers and other gravitational-wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory’s Legacy Survey of Space and Time can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving NSs (∼tens per year) out to distances of several hundred megaparsecs. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of NS and other compact-object mergers, and yet unknown classes of gravitational-wave events.

The current Cepheid-calibrated distance ladder measurement of H 0 is reported to be in tension with the values inferred from the cosmic microwave background (CMB), assuming standard cosmology. However, some tip of the red giant branch (TRGB) estimates report H 0 in better agreement with the CMB. Hence, it is critical to reduce systematic uncertainties in local measurements to understand the Hubble tension. In this paper, we propose a uniform distance ladder between the second and third rungs, combining Type Ia supernovae (SNe Ia) observed by the Zwicky Transient Facility (ZTF) with a TRGB calibration of their absolute luminosity. A large, volume-limited sample of both calibrator and Hubble flow SNe Ia from the same survey minimizes two of the largest sources of systematics: host-galaxy bias and nonuniform photometric calibration. We present results from a pilot study using the existing TRGB distance to the host galaxy of ZTF SN Ia SN 2021rhu (aka ZTF21abiuvdk) in NGC7814. Combining the ZTF calibrator with a volume-limited sample from the first data release of ZTF Hubble flow SNe Ia, we infer H 0 = 76.94 ± 6.4 km s −1 Mpc −1 , an 8.3% measurement. The error budget is dominated by the single object calibrating the SN Ia luminosity in this pilot study. However, the ZTF sample includes already five other SNe Ia within ∼20 Mpc for which TRGB distances can be obtained with the Hubble Space Telescope. Finally, we present the prospects of building this distance ladder out to 80 Mpc with James Webb Space Telescope observations of more than 100 ZTF SNe Ia.

ABSTRACT Forthcoming time-domain surveys, such as the Rubin Observatory Legacy Survey of Space and Time, will vastly increase samples of supernovae (SNe) and other optical transients, requiring new data-driven techniques to analyse their photometric light curves. Here, we present the ‘Python for Intelligent Supernova-COsmology Light-curve Analysis’ (PISCOLA ), an open source data-driven light-curve fitter using Gaussian Processes that can estimate rest-frame light curves of transients without the need for an underlying light-curve template. We test PISCOLA on large-scale simulations of type Ia SNe (SNe Ia) to validate its performance, and show it successfully retrieves rest-frame peak magnitudes for average survey cadences of up to 7 d. We also compare to the existing SN Ia light-curve fitter SALT2 on real data, and find only small (but significant) disagreements for different light-curve parameters. As a proof-of-concept of an application of PISCOLA , we decomposed and analysed the PISCOLA rest-frame light curves of SNe Ia from the Pantheon SN Ia sample with Non-Negative Matrix Factorization. Our new parametrization provides a similar performance to existing light-curve fitters such as SALT2. We further derived a SN Ia colour law from PISCOLA fits over ∼3500–7000 Å, and find agreement with the SALT2 colour law and with reddening laws with total-to-selective extinction ratio RV ≲ 3.1.

Studying galaxies at different cosmic epochs entails several observational effects that need to be taken into account to compare populations across a large time span in a consistent manner. We use a sample of 166 nearby galaxies that hosted type Ia supernovae (SNe Ia) and have been observed with the integral field spectrograph MUSE through the AMUSING survey. Here, we present a study of the systematic errors and bias in the host stellar mass with increasing redshifts that are generally overlooked in SNe Ia cosmological analyses. We simulate observations at different redshifts (0.1<z<2.0) using four photometric bands (griz, similar to the Dark Energy Survey-SN program) to then estimate the host galaxy properties across cosmic time. We find that stellar masses are systematically underestimated as we move towards higher redshifts, due mostly to different rest-frame wavelength coverage, with differences reaching 0.3 dex at z~1. We have used the newly derived corrections as a function of redshift to correct the stellar masses of a known sample of SN Ia hosts and derive cosmological parameters. We show that these corrections have a small impact on the derived cosmological parameters. The most affected is the value of the mass step $\Delta_M$, which is reduced by $\sim$0.004 (6% lower). The dark energy equation of state parameter $w$ changes by $\Delta w \sim $0.006 (0.6% higher) and the value of $\Omega_m$ increases at most by 0.001 ($\sim$0.3%), all within the derived uncertainties of the model. While the systematic error found in the estimate of the host stellar mass does not significantly affect the derived cosmological parameters, it is an important source of a systematic error that one should correct for as we enter a new era of precision cosmology.

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 similar to 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 similar to 20-yr time baseline spans similar to 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.

We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev–Zel’dovich effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius is consistent with CDM-only simulations and is located at 2.4 − 0.4 + 0.3 Mpc h − 1 . We split the galaxies on color and find significant differences in their profile shapes. Red and green-valley galaxies show a splashback-like minimum in their slope profile consistent with theory, while the bluest galaxies show a weak feature at a smaller radius. We develop a mapping of galaxies to subhalos in simulations and assign colors based on infall time onto their hosts. We find that the shift in location of the steepest slope and different profile shapes can be mapped to the average time of infall of galaxies of different colors. The steepest slope traces a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in clusters over 3.2 Gyr, green galaxies about 2.2 Gyr, while blue galaxies have been accreted most recently and have not reached apocenter. Using the full radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching occurs after a delay of about a gigayear and that galaxies quench rapidly thereafter with an exponential timescale of 0.6 Gyr.

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σ.

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 deg(2) 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 similar to 23.5 mag) powered by less massive black holes (with masses greater than or similar to 10(8.5)M(circle dot)) 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 similar to 7 x 10(-4)-3 x 10(-3) accounting for the look-elsewhere effect) periodicity with observed periods of similar to 3-5 yr (i.e. 1-2 yr in rest frame) having similar to 4-6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of similar to 0.8(-0.3)(+0.5) per cent or similar to 1.1(-0.5)(+0.7) quasar per deg(2). 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.