Showing 1–50 of 200 results for author: Baxter, J

A detailed study of the very-high-energy Crab pulsar emission with the LST-1

Authors: CTA-LST Project, :, K. Abe, S. Abe, A. Abhishek, F. Acero, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González , et al. (272 additional authors not shown)

Abstract: Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov… ▽ More Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims: To study the Crab pulsar emission with the LST-1, improving and complementing the results from other telescopes. These observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods: We analyzed a total of $\sim$103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles less than 50 degrees. A new analysis of the Fermi-LAT data was also performed, including $\sim$14 years of observations. Results: The Crab pulsar phaseogram, long-term light-curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for P1 and up to 700 GeV for P2. The pulsed emission is detected with a significance of 15.2$σ$. The two characteristic emission peaks of the Crab pulsar are clearly detected (>10$σ$), as well as the so-called bridge emission (5.7$σ$). We find that both peaks are well described by power laws, with spectral indices of $\sim$3.44 and $\sim$3.03 respectively. The joint analysis of Fermi-LAT and LST-1 data shows a good agreement between both instruments in the overlapping energy range. The detailed results obtained in the first observations of the Crab pulsar with LST-1 show the potential that CTAO will have to study this type of sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

Constraints on Lorentz invariance violation from the extraordinary Mrk 421 flare of 2014 using a novel analysis method

Authors: MAGIC Collaboration, S. Abe, J. Abhir, A. Abhishek, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete , et al. (192 additional authors not shown)

Abstract: The Lorentz Invariance Violation (LIV), a proposed consequence of certain quantum gravity (QG) scenarios, could instigate an energy-dependent group velocity for ultra-relativistic particles. This energy dependence, although suppressed by the massive QG energy scale $E_\mathrm{QG}$, expected to be on the level of the Planck energy $1.22 \times 10^{19}$ GeV, is potentially detectable in astrophysica… ▽ More The Lorentz Invariance Violation (LIV), a proposed consequence of certain quantum gravity (QG) scenarios, could instigate an energy-dependent group velocity for ultra-relativistic particles. This energy dependence, although suppressed by the massive QG energy scale $E_\mathrm{QG}$, expected to be on the level of the Planck energy $1.22 \times 10^{19}$ GeV, is potentially detectable in astrophysical observations. In this scenario, the cosmological distances traversed by photons act as an amplifier for this effect. By leveraging the observation of a remarkable flare from the blazar Mrk\,421, recorded at energies above 100 GeV by the MAGIC telescopes on the night of April 25 to 26, 2014, we look for time delays scaling linearly and quadratically with the photon energies. Using for the first time in LIV studies a binned-likelihood approach we set constraints on the QG energy scale. For the linear scenario, we set $95\%$ lower limits $E_\mathrm{QG}>2.7\times10^{17}$ GeV for the subluminal case and $E_\mathrm{QG}> 3.6 \times10^{17}$ GeV for the superluminal case. For the quadratic scenario, the $95\%$ lower limits for the subluminal and superluminal cases are $E_\mathrm{QG}>2.6 \times10^{10}$ GeV and $E_\mathrm{QG}>2.5\times10^{10}$ GeV, respectively. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Inference of Causal Networks using a Topological Threshold

Authors: Filipe Barroso, Diogo Gomes, Gareth J. Baxter

Abstract: We propose a constraint-based algorithm, which automatically determines causal relevance thresholds, to infer causal networks from data. We call these topological thresholds. We present two methods for determining the threshold: the first seeks a set of edges that leaves no disconnected nodes in the network; the second seeks a causal large connected component in the data. We tested these methods… ▽ More We propose a constraint-based algorithm, which automatically determines causal relevance thresholds, to infer causal networks from data. We call these topological thresholds. We present two methods for determining the threshold: the first seeks a set of edges that leaves no disconnected nodes in the network; the second seeks a causal large connected component in the data. We tested these methods both for discrete synthetic and real data, and compared the results with those obtained for the PC algorithm, which we took as the benchmark. We show that this novel algorithm is generally faster and more accurate than the PC algorithm. The algorithm for determining the thresholds requires choosing a measure of causality. We tested our methods for Fisher Correlations, commonly used in PC algorithm (for instance in \cite{kalisch2005}), and further proposed a discrete and asymmetric measure of causality, that we called Net Influence, which provided very good results when inferring causal networks from discrete data. This metric allows for inferring directionality of the edges in the process of applying the thresholds, speeding up the inference of causal DAGs. △ Less

Submitted 21 April, 2024; originally announced April 2024.

Comments: 17 pages, 12 figures

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

The variability patterns of the TeV blazar PG 1553+113 from a decade of MAGIC and multi-band observations

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (242 additional authors not shown)

Abstract: PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m a… ▽ More PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

A Systematic Review of Low-Rank and Local Low-Rank Matrix Approximation in Big Data Medical Imaging

Authors: Sisipho Hamlomo, Marcellin Atemkeng, Yusuf Brima, Chuneeta Nunhokee, Jeremy Baxter

Abstract: The large volume and complexity of medical imaging datasets are bottlenecks for storage, transmission, and processing. To tackle these challenges, the application of low-rank matrix approximation (LRMA) and its derivative, local LRMA (LLRMA) has demonstrated potential. A detailed analysis of the literature identifies LRMA and LLRMA methods applied to various imaging modalities, and the challenges… ▽ More The large volume and complexity of medical imaging datasets are bottlenecks for storage, transmission, and processing. To tackle these challenges, the application of low-rank matrix approximation (LRMA) and its derivative, local LRMA (LLRMA) has demonstrated potential. A detailed analysis of the literature identifies LRMA and LLRMA methods applied to various imaging modalities, and the challenges and limitations associated with existing LRMA and LLRMA methods are addressed. We note a significant shift towards a preference for LLRMA in the medical imaging field since 2015, demonstrating its potential and effectiveness in capturing complex structures in medical data compared to LRMA. Acknowledging the limitations of shallow similarity methods used with LLRMA, we suggest advanced semantic image segmentation for similarity measure, explaining in detail how it can be used to measure similar patches and its feasibility. We note that LRMA and LLRMA are mainly applied to unstructured medical data, and we propose extending their application to different medical data types, including structured and semi-structured. This paper also discusses how LRMA and LLRMA can be applied to regular data with missing entries and the impact of inaccuracies in predicting missing values and their effects. We discuss the impact of patch size and propose the use of random search (RS) to determine the optimal patch size. To enhance feasibility, a hybrid approach using Bayesian optimization and RS is proposed, which could improve the application of LRMA and LLRMA in medical imaging. △ Less

Submitted 27 May, 2024; v1 submitted 21 February, 2024; originally announced February 2024.

SuperBench: Improving Cloud AI Infrastructure Reliability with Proactive Validation

Authors: Yifan Xiong, Yuting Jiang, Ziyue Yang, Lei Qu, Guoshuai Zhao, Shuguang Liu, Dong Zhong, Boris Pinzur, Jie Zhang, Yang Wang, Jithin Jose, Hossein Pourreza, Jeff Baxter, Kushal Datta, Prabhat Ram, Luke Melton, Joe Chau, Peng Cheng, Yongqiang Xiong, Lidong Zhou

Abstract: Reliability in cloud AI infrastructure is crucial for cloud service providers, prompting the widespread use of hardware redundancies. However, these redundancies can inadvertently lead to hidden degradation, so called "gray failure", for AI workloads, significantly affecting end-to-end performance and concealing performance issues, which complicates root cause analysis for failures and regressions… ▽ More Reliability in cloud AI infrastructure is crucial for cloud service providers, prompting the widespread use of hardware redundancies. However, these redundancies can inadvertently lead to hidden degradation, so called "gray failure", for AI workloads, significantly affecting end-to-end performance and concealing performance issues, which complicates root cause analysis for failures and regressions. We introduce SuperBench, a proactive validation system for AI infrastructure that mitigates hidden degradation caused by hardware redundancies and enhances overall reliability. SuperBench features a comprehensive benchmark suite, capable of evaluating individual hardware components and representing most real AI workloads. It comprises a Validator which learns benchmark criteria to clearly pinpoint defective components. Additionally, SuperBench incorporates a Selector to balance validation time and issue-related penalties, enabling optimal timing for validation execution with a tailored subset of benchmarks. Through testbed evaluation and simulation, we demonstrate that SuperBench can increase the mean time between incidents by up to 22.61x. SuperBench has been successfully deployed in Azure production, validating hundreds of thousands of GPUs over the last two years. △ Less

Submitted 7 June, 2024; v1 submitted 9 February, 2024; originally announced February 2024.

Comments: USENIX ATC '24

Performance and first measurements of the MAGIC Stellar Intensity Interferometer

Authors: MAGIC Collaboration, S. Abe, J. Abhir, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, E. Bernardini, M. Bernardos, J. Bernete, A. Berti , et al. (195 additional authors not shown)

Abstract: In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the sys… ▽ More In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the system was upgraded and now features a real-time, dead-time-free, 4-channel, GPU-based correlator. These hardware modifications allow seamless transitions between MAGIC's standard very-high-energy gamma-ray observations and optical interferometry measurements within seconds. We establish the feasibility and potential of employing IACTs as competitive optical Intensity Interferometers with minimal hardware adjustments. The measurement of a total of 22 stellar diameters are reported, 9 corresponding to reference stars with previous comparable measurements, and 13 with no prior measurements. A prospective implementation involving telescopes from the forthcoming Cherenkov Telescope Array Observatory's northern hemisphere array, such as the first prototype of its Large-Sized Telescopes, LST-1, is technically viable. This integration would significantly enhance the sensitivity of the current system and broaden the UV-plane coverage. This advancement would enable the system to achieve competitive sensitivity with the current generation of long-baseline optical interferometers over blue wavelengths. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 18 pages, 13 figures, submitted to MNRAS

Applying Simulation-Based Inference to Spectral and Spatial Information from the Galactic Center Gamma-Ray Excess

Authors: Katharena Christy, Eric J. Baxter, Jason Kumar

Abstract: The two most favored explanations of the Fermi Galactic Center gamma-ray excess (GCE) are millisecond pulsars and self annihilation of the smooth dark matter halo of the galaxy. In order to distinguish between these possibilities, we would like to optimally use all information in the available data, including photon direction and energy information. To date, analyses of the GCE have generally trea… ▽ More The two most favored explanations of the Fermi Galactic Center gamma-ray excess (GCE) are millisecond pulsars and self annihilation of the smooth dark matter halo of the galaxy. In order to distinguish between these possibilities, we would like to optimally use all information in the available data, including photon direction and energy information. To date, analyses of the GCE have generally treated directional and energy information separately, or have ignored one or the other completely. Here, we develop a method for analyzing the GCE that relies on simulation-based inference with neural posterior models to jointly analyze photon directional and spectral information while correctly accounting for the spatial and energy resolution of the telescope, here assumed to be the Fermi Large Area Telescope (LAT). Our results also have implications for analyses of the diffuse gamma-ray background, which we discuss. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 17 pages, 5 figures

GODMAX: Modeling gas thermodynamics and matter distribution using JAX

Authors: Shivam Pandey, Jaime Salcido, Chun-Hao To, J. Colin Hill, Dhayaa Anbajagane, Eric J. Baxter, Ian G. McCarthy

Abstract: We introduce GODMAX (Gas thermODynamics and Matter distribution using jAX), a novel code designed to calculate correlations between the cosmological matter distribution and various gas thermodynamic quantities. Utilizing the extensive ANTILLES suite of 200 hydrodynamical simulations with a diverse range of baryonic feedback strengths, we jointly fit the 3D profiles of total matter distribution, el… ▽ More We introduce GODMAX (Gas thermODynamics and Matter distribution using jAX), a novel code designed to calculate correlations between the cosmological matter distribution and various gas thermodynamic quantities. Utilizing the extensive ANTILLES suite of 200 hydrodynamical simulations with a diverse range of baryonic feedback strengths, we jointly fit the 3D profiles of total matter distribution, electron density, and electron pressure across various halo masses and redshifts. By accommodating significant variations in gas profiles expected due to baryonic feedback, solving exact hydrostatic equilibrium equation and offering flexible modeling of non-thermal pressure support, GODMAX has the capability to jointly fit all these profiles within the measurement uncertainties. This advancement enables, for the first time, robust joint analyses of multiple cosmic probes, including the kinetic and thermal Sunyaev-Zel'dovich effect, weak lensing, and X-ray observations. Furthermore, the model accurately captures correlations between the total matter power suppression due to baryonic feedback and local average thermodynamic quantities, such as the baryon fraction and integrated tSZ effect, in high-mass halos, aligning with observations from hydrodynamical simulations. Looking ahead, we forecast the expected constraints on cosmological and baryonic parameters from upcoming weak lensing catalogs from the LSST and tSZ maps from the Simons Observatory. This analysis underscores the importance of cross-correlations between weak lensing and tSZ in enhancing parameter constraints by resolving major systematic uncertainties due to baryonic physics. The GODMAX code leverages the JAX library, resulting in a fully differentiable halo model with native GPU compilation support. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: 17 pages, 9 figures

Inferring galaxy cluster masses from cosmic microwave background lensing with neural simulation based inference

Authors: Eric J. Baxter, Shivam Pandey

Abstract: Gravitational lensing by massive galaxy clusters distorts the observed cosmic microwave background (CMB) on arcminute scales, and these distortions carry information about cluster masses. Standard approaches to extracting cluster mass constraints from the CMB cluster lensing signal are either sub-optimal, ignore important physical or observational effects, are computationally intractable, or requi… ▽ More Gravitational lensing by massive galaxy clusters distorts the observed cosmic microwave background (CMB) on arcminute scales, and these distortions carry information about cluster masses. Standard approaches to extracting cluster mass constraints from the CMB cluster lensing signal are either sub-optimal, ignore important physical or observational effects, are computationally intractable, or require additional work to turn the lensing measurements into constraints on cluster masses. We apply simulation based inference (SBI) using neural likelihood models to the problem. We show that in circumstances where the exact likelihood can be computed, the SBI constraints on cluster masses are in agreement with the exact likelihood, demonstrating that the SBI constraints are close to optimal. In scenarios where the exact likelihood cannot be feasibly computed, SBI still recovers unbiased estimates of individual cluster masses and combined constraints from multiple clusters. SBI will be a powerful tool for constraining the masses of galaxy clusters detected by future cosmic surveys. Code to run the analyses presented here will be made publicly available. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 11 pages, 2 figures

Insights into the broad-band emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, J. Bernete, A. Berti, J. Besenrieder , et al. (239 additional authors not shown)

Abstract: We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and… ▽ More We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several instruments covering the optical and radio bands. During the IXPE pointings, the VHE state is close to the average behavior with a 0.2-1 TeV flux of 20%-50% the emission of the Crab Nebula. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings in March 2022 with a synchrotron peak frequency >1 keV. For the third IXPE pointing in July 2022, the synchrotron peak shifts towards lower energies and the optical/X-ray polarization degrees drop. The X-ray polarization is systematically higher than at lower energies, suggesting an energy-stratification of the jet. While during the IXPE epochs the polarization angle in the X-ray, optical and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. We model the broad-band spectra simultaneous to the IXPE pointings assuming a compact zone dominating in the X-rays and VHE, and an extended zone stretching further downstream the jet dominating the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change of magnetization and/or emission region size, which directly connects the shift of the synchrotron peak to lower energies with the drop in polarization degree. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. 19 pages, 9 figures. Corresponding authors: Lea Heckmann, Axel Arbet Engels, David Paneque

Constraints on axion-like particles with the Perseus Galaxy Cluster with MAGIC

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete, A. Berti , et al. (189 additional authors not shown)

Abstract: Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the r… ▽ More Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the radio galaxy in the center of the cluster: NGC 1275. By modeling the magnetic field surrounding this target, we searched for spectral indications of ALP presence. Despite finding no statistical evidence of ALP signatures, we were able to exclude ALP models in the sub-micro electronvolt range. Our analysis improved upon previous work by calculating the full likelihood and statistical coverage for all considered models across the parameter space. Consequently, we achieved the most stringent limits to date for ALP masses around 50 neV, with cross sections down to $g_{aγ} = 3 \times 10^{-12}$ GeV$^{-1}$. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 25 pages, 10 figures, accepted for publication in Physics of the Dark Universe

First characterization of the emission behavior of Mrk421 from radio to VHE gamma rays with simultaneous X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland , et al. (229 additional authors not shown)

Abstract: We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured… ▽ More We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured 0.2-1 TeV flux is close to the quiescent state and ranges from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the VHE and X-ray emission are positively correlated at a $4σ$ significance level. The IXPE measurements unveil a X-ray polarization degree that is a factor of 2-5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift-XRT monitoring unveils an X-ray flux increase with a clear spectral hardening. It suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counter-clockwise), implying important changes in the particle acceleration efficiency on $\sim$hour timescales. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 29 pages, 22 figures. Corresponding authors: Axel Arbet Engels, Felix Schmuckermaier, David Paneque

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

MAGIC detection of GRB 201216C at $z=1.1$

Authors: H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (195 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus it is the farthest one detected at very high energies. The emission above \SI{70}{\GeV} of GRB~201216C… ▽ More Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus it is the farthest one detected at very high energies. The emission above \SI{70}{\GeV} of GRB~201216C is modelled together with multi-wavelength data within a synchrotron and synchrotron-self Compton (SSC) scenario. We find that SSC can explain the broadband data well from the optical to the very-high-energy band. For the late-time radio data, a different component is needed to account for the observed emission. Differently from previous GRBs detected in the very-high-energy range, the model for GRB~201216C strongly favors a wind-like medium. The model parameters have values similar to those found in past studies of the afterglows of GRBs detected up to GeV energies. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 13 pages, 6 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 3, January 2024, Pages 5856-5867

Multi-year characterisation of the broad-band emission from the intermittent extreme BL Lac 1ES~2344+514

Authors: H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland, O. Blanch , et al. (210 additional authors not shown)

Abstract: The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, f… ▽ More The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $ν_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $ν_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 682, A114 (2024)

Performance of the joint LST-1 and MAGIC observations evaluated with Crab Nebula data

Authors: H. Abe, K. Abe, S. Abe, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, T. Aniello, S. Ansoldi, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, D. Baack, A. Babić, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batković , et al. (344 additional authors not shown)

Abstract: Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both syste… ▽ More Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both systems. Methods. We describe the joint LST-1+MAGIC analysis pipeline and use simultaneous Crab Nebula observations and Monte Carlo simulations to assess the performance of the three-telescope system. The addition of the LST-1 telescope allows the recovery of events in which one of the MAGIC images is too dim to survive analysis quality cuts. Results. Thanks to the resulting increase in the collection area and stronger background rejection, we find a significant improvement in sensitivity, allowing the detection of 30% weaker fluxes in the energy range between 200 GeV and 3 TeV. The spectrum of the Crab Nebula, reconstructed in the energy range ~60 GeV to ~10 TeV, is in agreement with previous measurements. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 680, A66 (2023)

A Strategy for Identifying Periodic Sources Contributing to the Galactic Center Excess

Authors: Eric J. Baxter, Jason Kumar

Abstract: The origin of the Galactic Center gamma-ray excess has not been conclusively determined after over a decade of careful study. The two most widely discussed possibilities are a population of millisecond pulsars (MSPs), and annihilation of dark matter particles. In contrast with annihilating dark matter, MSPs are expected to produce periodic emission. We show that even though the number of photons c… ▽ More The origin of the Galactic Center gamma-ray excess has not been conclusively determined after over a decade of careful study. The two most widely discussed possibilities are a population of millisecond pulsars (MSPs), and annihilation of dark matter particles. In contrast with annihilating dark matter, MSPs are expected to produce periodic emission. We show that even though the number of photons contributing to the excess is small, there is potentially sufficient information in the data from Fermi to detect a periodic MSP signal. Such a detection would definitively prove that at least some fraction of the excess is due to MSPs. We argue that this conclusion is robust to potential timing perturbations of the gamma-ray photons, such as those due to Earth's orbit, even if the number of parameters that must be used to model the perturbations is $\sim 7$. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 7 pages, 2 figures

Cosmological shocks around galaxy clusters: A coherent investigation with DES, SPT & ACT

Authors: D. Anbajagane, C. Chang, E. J. Baxter, S. Charney, M. Lokken, M. Aguena, S. Allam, O. Alves, A. Amon, R. An, F. Andrade-Oliveira, D. Bacon, N. Battaglia, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, L. Bleem, S. Bocquet, J. R. Bond, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, R. Chen, A. Choi , et al. (89 additional authors not shown)

Abstract: We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around… ▽ More We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around $10^5$ clusters with mass and redshift ranges $10^{13.7} < M_{\rm 200m}/M_\odot < 10^{15.5}$ and $0.1 < z < 2$, and the total sky coverage of the maps is $\approx 15,000 \,\,{\rm deg}^2$. We find a clear pressure deficit at $R/R_{\rm 200m}\approx 1.1$ in SZ profiles around both ACT and SPT clusters, estimated at $6σ$ significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2-5 compared to predictions, and we discuss possible interpretations of this behavior. An oriented stacking of clusters -- where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogs -- shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius. △ Less

Submitted 12 December, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: [v2]: Version accepted to MNRAS

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

Late Time Modification of Structure Growth and the S8 Tension

Authors: Meng-Xiang Lin, Bhuvnesh Jain, Marco Raveri, Eric J. Baxter, Chihway Chang, Marco Gatti, Sujeong Lee, Jessica Muir

Abstract: The $S_8$ tension between low-redshift galaxy surveys and the primary CMB signals a possible breakdown of the $Λ$CDM model. Recently differing results have been obtained using low-redshift galaxy surveys and the higher redshifts probed by CMB lensing, motivating a possible time-dependent modification to the growth of structure. We investigate a simple phenomenological model in which the growth of… ▽ More The $S_8$ tension between low-redshift galaxy surveys and the primary CMB signals a possible breakdown of the $Λ$CDM model. Recently differing results have been obtained using low-redshift galaxy surveys and the higher redshifts probed by CMB lensing, motivating a possible time-dependent modification to the growth of structure. We investigate a simple phenomenological model in which the growth of structure deviates from the $Λ$CDM prediction at late times, in particular as a simple function of the dark energy density. Fitting to galaxy lensing, CMB lensing, BAO, and Supernovae datasets, we find significant evidence - 2.5 - 3$σ$, depending on analysis choices - for a non-zero value of the parameter quantifying a deviation from $Λ$CDM. The preferred model, which has a slower growth of structure below $z\sim 1$, improves the joint fit to the data over $Λ$CDM. While the overall fit is improved, there is weak evidence for galaxy and CMB lensing favoring different changes in the growth of structure. △ Less

Submitted 15 March, 2024; v1 submitted 30 August, 2023; originally announced August 2023.

Comments: Match the published version. Add more references. 6 pages, 5 figures, 1 table

Journal ref: Phys.Rev.D 109 (2024) 6, 063523

Beyond the 3rd moment: A practical study of using lensing convergence CDFs for cosmology with DES Y3

Authors: D. Anbajagane, C. Chang, A. Banerjee, T. Abel, M. Gatti, V. Ajani, A. Alarcon, A. Amon, E. J. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Chen, A. Choi, C. Davis, J. DeRose, H. T. Diehl, S. Dodelson, C. Doux, A. Drlica-Wagner, K. Eckert, J. Elvin-Poole , et al. (73 additional authors not shown)

Abstract: Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at… ▽ More Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at multiple scales as a summary of the galaxy lensing convergence field. Using a suite of N-body lightcone simulations, we show the CDFs' constraining power is modestly better than that of the 2nd and 3rd moments of the field, as they approximately capture the information from all moments of the field in a concise data vector. We then study the practical aspects of applying the CDFs to observational data, using the first three years of the Dark Energy Survey (DES Y3) data as an example, and compute the impact of different systematics on the CDFs. The contributions from the point spread function are 2-3 orders of magnitude below the cosmological signal, while those from reduced shear approximation contribute $\lesssim 1\%$ to the signal. Source clustering effects and baryon imprints contribute $1-10\%$. Enforcing scale cuts to limit systematics-driven biases in parameter constraints degrades these constraints a noticeable amount, and this degradation is similar for the CDFs and the moments. We also detect correlations between the observed convergence field and the shape noise field at $13σ$. We find that the non-Gaussian correlations in the noise field must be modeled accurately to use the CDFs, or other statistics sensitive to all moments, as a rigorous cosmology tool. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 21 pages, 12 figures

Cosmological constraints from the tomography of DES-Y3 galaxies with CMB lensing from ACT DR4

Authors: G. A. Marques, M. S. Madhavacheril, O. Darwish, S. Shaikh, M. Aguena, O. Alves, S. Avila, D. Bacon, E. J. Baxter, K. Bechtol, M. R. Becker, E. Bertin, J. Blazek, J. Richard Bond, D. Brooks, H. Cai, E. Calabrese, A. Carnero Rosell, M. Carrasco Kind J. Carretero, R. Cawthon, M. Crocce, L. N. da Costa, M. E. S. Pereira, J. De Vicente, S. Desai , et al. (70 additional authors not shown)

Abstract: We present a measurement of the cross-correlation between the MagLim galaxies selected from the Dark Energy Survey (DES) first three years of observations (Y3) and cosmic microwave background (CMB) lensing from the Atacama Cosmology Telescope (ACT) Data Release 4 (DR4), reconstructed over $\sim 436$ sq.deg. of the sky. Our galaxy sample, which covers $\sim 4143$ sq.deg., is divided into six redshi… ▽ More We present a measurement of the cross-correlation between the MagLim galaxies selected from the Dark Energy Survey (DES) first three years of observations (Y3) and cosmic microwave background (CMB) lensing from the Atacama Cosmology Telescope (ACT) Data Release 4 (DR4), reconstructed over $\sim 436$ sq.deg. of the sky. Our galaxy sample, which covers $\sim 4143$ sq.deg., is divided into six redshift bins spanning the redshift range of $0.20<z<1.05$. We adopt a blinding procedure until passing all consistency and systematics tests. After imposing scale cuts for the cross-power spectrum measurement, we reject the null hypothesis of no correlation at 9.1σ. We constrain cosmological parameters from a joint analysis of galaxy and CMB lensing-galaxy power spectra considering a flat \LCDM model, marginalized over 23 astrophysical and systematic nuisance parameters. We find the clustering amplitude $S_8\equiv σ_8 (Ω_m/0.3)^{0.5} = 0.75^{+0.04}_{-0.05}$. In addition, we constrain the linear growth of cosmic structure as a function of redshift. Our results are consistent with recent DES Y3 analyses and suggest a preference for a lower $S_8$ compared to results from measurements of CMB anisotropies by the Planck satellite, although at a mild level ($< 2 σ$) of statistical significance. △ Less

Submitted 11 October, 2023; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 42 pages, 17 figures

Observations of the Crab Nebula and Pulsar with the Large-Sized Telescope Prototype of the Cherenkov Telescope Array

Authors: CTA-LST Project, :, H. Abe, K. Abe, S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini , et al. (467 additional authors not shown)

Abstract: CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing a… ▽ More CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing astronomical observations in November 2019, during its commissioning phase, and it has been taking data since then. We present the first LST-1 observations of the Crab Nebula, the standard candle of very-high energy gamma-ray astronomy, and use them, together with simulations, to assess the basic performance parameters of the telescope. The data sample consists of around 36 hours of observations at low zenith angles collected between November 2020 and March 2022. LST-1 has reached the expected performance during its commissioning period - only a minor adjustment of the preexisting simulations was needed to match the telescope behavior. The energy threshold at trigger level is estimated to be around 20 GeV, rising to $\simeq 30$ GeV after data analysis. Performance parameters depend strongly on energy, and on the strength of the gamma-ray selection cuts in the analysis: angular resolution ranges from 0.12 to 0.40 degrees, and energy resolution from 15 to 50%. Flux sensitivity is around 1.1% of the Crab Nebula flux above 250 GeV for a 50-h observation (12% for 30 minutes). The spectral energy distribution (in the 0.03 - 30 TeV range) and the light curve obtained for the Crab Nebula agree with previous measurements, considering statistical and systematic uncertainties. A clear periodic signal is also detected from the pulsar at the center of the Nebula. △ Less

Submitted 19 July, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted in ApJ. v3: updated author list and acknowledgements, fixed typos and other minor issues

The Kinematic Sunyaev-Zel'dovich Effect with ACT, DES, and BOSS: a Novel Hybrid Estimator

Authors: M. Mallaby-Kay, S. Amodeo, J. C. Hill, M. Aguena, S. Allam, O. Alves, J. Annis, N. Battaglia, E. S. Battistelli, E. J. Baxter, K. Bechtol, M. R. Becker, E. Bertin, J. R. Bond, D. Brooks, E. Calabrese, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, A. Choi, M. Crocce, L. N. da Costa, M. E. S. Pereira, J. De Vicente, S. Desai , et al. (58 additional authors not shown)

Abstract: The kinematic and thermal Sunyaev-Zel'dovich (kSZ and tSZ) effects probe the abundance and thermodynamics of ionized gas in galaxies and clusters. We present a new hybrid estimator to measure the kSZ effect by combining cosmic microwave background temperature anisotropy maps with photometric and spectroscopic optical survey data. The method interpolates a velocity reconstruction from a spectroscop… ▽ More The kinematic and thermal Sunyaev-Zel'dovich (kSZ and tSZ) effects probe the abundance and thermodynamics of ionized gas in galaxies and clusters. We present a new hybrid estimator to measure the kSZ effect by combining cosmic microwave background temperature anisotropy maps with photometric and spectroscopic optical survey data. The method interpolates a velocity reconstruction from a spectroscopic catalog at the positions of objects in a photometric catalog, which makes it possible to leverage the high number density of the photometric catalog and the precision of the spectroscopic survey. Combining this hybrid kSZ estimator with a measurement of the tSZ effect simultaneously constrains the density and temperature of free electrons in the photometrically selected galaxies. Using the 1000 deg2 of overlap between the Atacama Cosmology Telescope (ACT) Data Release 5, the first three years of data from the Dark Energy Survey (DES), and the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12, we detect the kSZ signal at 4.8$σ$ and reject the null (no-kSZ) hypothesis at 5.1$σ$. This corresponds to 2.0$σ$ per 100,000 photometric objects with a velocity field based on a spectroscopic survey with 1/5th the density of the photometric catalog. For comparison, a recent ACT analysis using exclusively spectroscopic data from BOSS measured the kSZ signal at 2.1$σ$ per 100,000 objects. Our derived constraints on the thermodynamic properties of the galaxy halos are consistent with previous measurements. With future surveys, such as the Dark Energy Spectroscopic Instrument and the Rubin Observatory Legacy Survey of Space and Time, we expect that this hybrid estimator could result in measurements with significantly better signal-to-noise than those that rely on spectroscopic data alone. △ Less

Submitted 14 August, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Comments: 19 pages, 15 figures - matches published version

Journal ref: Phys. Rev. D 108, 023516 - Published 18 July 2023

A Finite-Difference Time-Domain approach for dispersive magnetic media

Authors: Jasmin Graf, Joshua Baxter, Sanchar Sharma, Silvia Viola Kusminskiy, Lora Ramunno

Abstract: We extend the Finite-Difference Time-Domain method to treat dispersive magnetic media by incorporating magneto-optical effects through a frequency-dependent permittivity tensor. For benchmarking our method, we consider the light scattering on a magnetic sphere in the Mie regime. We first derive the analytical scattering expressions which predict a peak broadening in the scattering efficiency due t… ▽ More We extend the Finite-Difference Time-Domain method to treat dispersive magnetic media by incorporating magneto-optical effects through a frequency-dependent permittivity tensor. For benchmarking our method, we consider the light scattering on a magnetic sphere in the Mie regime. We first derive the analytical scattering expressions which predict a peak broadening in the scattering efficiency due to the atomic energy level splitting in the presence of a magnetic field, together with an additional rotated part in the scattered field profile due to the Faraday rotation. We show that our numerical method is able to capture the main scattering features and discuss its limitations and possible improvements in accuracy. △ Less

Submitted 5 May, 2023; originally announced May 2023.

The impact of halo concentration on the Sunyaev Zel'dovich effect signal from massive galaxy clusters

Authors: Eric J. Baxter, Shivam Pandey, Susmita Adhikari, Weiguang Cui, Tae-hyeon Shin, Qingyang Li, Elena Rasia

Abstract: The Sunyaev Zel'dovich (SZ) effect is sensitive to the pressure of ionized gas inside galaxy clusters, which is in turn controlled largely by the gravitational potential of the cluster. Changing the concentration parameter describing the cluster mass distribution impacts the gravitational potential and thus the cluster SZ signal, with implications for cosmological and other analyses of SZ-selected… ▽ More The Sunyaev Zel'dovich (SZ) effect is sensitive to the pressure of ionized gas inside galaxy clusters, which is in turn controlled largely by the gravitational potential of the cluster. Changing the concentration parameter describing the cluster mass distribution impacts the gravitational potential and thus the cluster SZ signal, with implications for cosmological and other analyses of SZ-selected clusters. We investigate the concentration-SZ relation in theory and simulations. We find that the impact of concentration on the inner SZ profile ($R \lesssim 0.75 R_{200c}$) can be captured with standard polytropic gas models. However, we find that such models do a poor job of reproducing the outer SZ profiles ($R \gtrsim 0.75 R_{200c}$) and the relation between the integrated SZ signal, $Y$, and concentration. This disagreement results from a sharp truncation of the gas pressure profile near the splashback radius, likely caused by virial shocks. We develop a simple description of the truncation that leads to a good match with the simulated SZ profiles out to several $R_{200c}$ for clusters of varying mass and concentration, and that also accurately predicts the concentration-$Y$ relationship. Finally, we determine how inference of the linear bias parameter and splashback radius for SZ-selected clusters can be biased by ignoring the concentration dependence of the SZ signal, finding that bias to the former is essentially negligible, while bias to the latter can be as much as 2\%. △ Less

Submitted 27 November, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 16 pages, 11 figures; replaced to match version accepted by MNRAS. Revised version uses an updated model for nonthermal pressure support, but conclusions are not impacted

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

Understanding the leading indicators of hospital admissions from COVID-19 across successive waves in the UK

Authors: Jonathon Mellor, Christopher E Overton, Martyn Fyles, Liam Chawner, James Baxter, Tarrion Baird, Thomas Ward

Abstract: Following the UK Government's Living with COVID-19 Strategy and the end of universal testing, hospital admissions are an increasingly important measure of COVID-19 pandemic pressure. Understanding leading indicators of admissions at National Health Service (NHS) Trust, regional and national geographies help health services plan capacity needs and prepare for ongoing pressures. We explored the spat… ▽ More Following the UK Government's Living with COVID-19 Strategy and the end of universal testing, hospital admissions are an increasingly important measure of COVID-19 pandemic pressure. Understanding leading indicators of admissions at National Health Service (NHS) Trust, regional and national geographies help health services plan capacity needs and prepare for ongoing pressures. We explored the spatio-temporal relationships of leading indicators of hospital pressure across successive waves of SARS-CoV-2 incidence in England. This includes an analysis of internet search volume values from Google Trends, NHS triage calls and online queries, the NHS COVID-19 App, lateral flow devices and the ZOE App. Data sources were analysed for their feasibility as leading indicators using linear and non-linear methods; granger causality, cross correlations and dynamic time warping at fine spatial scales. Consistent temporal and spatial relationships were found for some of the leading indicators assessed across resurgent waves of COVID-19. Google Trends and NHS queries consistently led admissions in over 70% of Trusts, with lead times ranging from 5-20 days, whereas an inconsistent relationship was found for the ZOE app, NHS COVID-19 App, and rapid testing, that diminished with granularity, showing limited autocorrelation of leads between -7 to 7 days. This work shows that novel syndromic surveillance data has utility for understanding the expected hospital burden at fine spatial scales. The analysis shows at low level geographies that some surveillance sources can predict hospital admissions, though care must be taken in relying on the lead times and consistency between waves. △ Less

Submitted 16 August, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Inverse design of optical pulse shapes for time-varying photonics

Authors: Joshua Baxter, Lora Ramunno

Abstract: There has been an explosion of interest in time-varying photonics due to the recent discovery and design of materials and metamaterials with strong, time-varying, nonlinear optical responses. This opens the door to novel optical phenomena including reciprocity breaking, frequency translation, and amplification that can be enhanced by optimizing the light-matter interaction. Although there has been… ▽ More There has been an explosion of interest in time-varying photonics due to the recent discovery and design of materials and metamaterials with strong, time-varying, nonlinear optical responses. This opens the door to novel optical phenomena including reciprocity breaking, frequency translation, and amplification that can be enhanced by optimizing the light-matter interaction. Although there has been recent interest in applying topology-based inverse design to this problem, we have decided to take a more novel approach. In this article, we will introduce a method for the inverse design of optical pulse shapes to enhance their interaction with time-varying media. We test our objective-first approach by maximizing the transmittance of optical pulses of equal intensity through time-varying media. Indeed, without requiring a change in pulse energy, we demonstrate large, broadband enhancements in the pulse energy transmission through the thin-films, including gain. Our final test includes maximizing pulse transmission through indium tin oxide, a time-varying medium when strongly pumped in its ENZ band. Through this work, we hope to inspire exploration of this new degree of freedom. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Inferring the impact of feedback on the matter distribution using the Sunyaev Zel'dovich effect: Insights from CAMELS simulations and ACT+DES data

Authors: Shivam Pandey, Kai Lehman, Eric J. Baxter, Yueying Ni, Daniel Anglés-Alcázar, Shy Genel, Francisco Villaescusa-Navarro, Ana Maria Delgado, Tiziana di Matteo

Abstract: Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Ex… ▽ More Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Extending the results of previous work to smaller halo masses and higher wavenumber, $k$, we find that the baryon fraction in halos contains significant information about the impact of feedback on the matter power spectrum. We explore how the thermal Sunyaev Zel'dovich (tSZ) signal from group-scale halos contains similar information. Using recent Dark Energy Survey (DES) weak lensing and Atacama Cosmology Telescope (ACT) tSZ cross-correlation measurements and models trained on CAMELS, we obtain $10\%$ constraints on feedback effects on the power spectrum at $k \sim 5\, h/{\rm Mpc}$. We show that with future surveys, it will be possible to constrain baryonic effects on the power spectrum to $\mathcal{O}(<1\%)$ at $k = 1\, h/{\rm Mpc}$ and $\mathcal{O}(3\%)$ at $k = 5\, h/{\rm Mpc}$ using the methods that we introduce here. Finally, we investigate the impact of feedback on the matter bispectrum, finding that tSZ observables are highly informative in this case. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 18 pages, 15 figures. Comments are welcome

Search for Gamma-ray Spectral Lines from Dark Matter Annihilation up to 100 TeV towards the Galactic Center with MAGIC

Authors: MAGIC Collaboration, H. Abe, S. Abe, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya , et al. (188 additional authors not shown)

Abstract: Line-like features in TeV $γ$-rays constitute a ''smoking gun'' for TeV-scale particle dark matter and new physics. Probing the Galactic Center region with ground-based Cherenkov telescopes enables the search for TeV spectral features in immediate association with a dense dark matter reservoir at a sensitivity out of reach for satellite $γ$-ray detectors, and direct detection and collider experime… ▽ More Line-like features in TeV $γ$-rays constitute a ''smoking gun'' for TeV-scale particle dark matter and new physics. Probing the Galactic Center region with ground-based Cherenkov telescopes enables the search for TeV spectral features in immediate association with a dense dark matter reservoir at a sensitivity out of reach for satellite $γ$-ray detectors, and direct detection and collider experiments. We report on 223 hours of observations of the Galactic Center region with the MAGIC stereoscopic telescope system reaching $γ$-ray energies up to 100 TeV. We improved the sensitivity to spectral lines at high energies using large-zenith-angle observations and a novel background modeling method within a maximum-likelihood analysis in the energy domain. No line-like spectral feature is found in our analysis. Therefore, we constrain the cross section for dark matter annihilation into two photons to $\langle σv \rangle \lesssim 5 \times 10^{-28}\,\mathrm{cm^3\,s^{-1}}$ at 1 TeV and $\langle σv \rangle \lesssim 1 \times 10^{-25}\,\mathrm{cm^3\,s^{-1}}$ at 100 TeV, achieving the best limits to date for a dark matter mass above 20 TeV and a cuspy dark matter profile at the Galactic Center. Finally, we use the derived limits for both cuspy and cored dark matter profiles to constrain supersymmetric wino models. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: Accepted in Phys. Rev. Lett., includes Supplemental materials. 9+4 pages, 3+3 figures, 3 tables. Corresponding authors: T. Inada, D. Kerszberg, M. Hütten

Report number: KEK-TH-2487, KEK-Cosmo-0307

MAGIC observations provide compelling evidence of the hadronic multi-TeV emission from the putative PeVatron SNR G106.3+2.7

Authors: MAGIC Collaboration, H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder , et al. (192 additional authors not shown)

Abstract: The SNR G106.3+2.7, detected at 1--100 TeV energies by different $γ$-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape which can be divided into a head and a tail region with different physical conditions. However, it is not identified in which region the 100 TeV emission is produced due to the limited position accuracy and/or angular resolution of exi… ▽ More The SNR G106.3+2.7, detected at 1--100 TeV energies by different $γ$-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape which can be divided into a head and a tail region with different physical conditions. However, it is not identified in which region the 100 TeV emission is produced due to the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear whether the origin of the $γ$-ray emission is leptonic or hadronic. With the better angular resolution provided by these new MAGIC data compared to earlier $γ$-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1$^\circ$ resolution at TeV energies. We detected extended $γ$-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detected a significant $γ$-ray emission with energies above 6.0 TeV from the tail region only suggests that the emissions above 10 TeV, detected with air shower experiments (Milagro, HAWC, Tibet AS$γ$ and LHAASO), are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of $\sim 1$ PeV for the tail region. Such a high energy emission in this middle-aged SNR (4--10 kyr) can be explained by considering the scenario that protons escaping from the SNR in the past interact with surrounding dense gases at present. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 13 pages, 7 figures, Accepted for publication in A&A

Journal ref: A&A 671, A12 (2023)

Birdwatch: Crowd Wisdom and Bridging Algorithms can Inform Understanding and Reduce the Spread of Misinformation

Authors: Stefan Wojcik, Sophie Hilgard, Nick Judd, Delia Mocanu, Stephen Ragain, M. B. Fallin Hunzaker, Keith Coleman, Jay Baxter

Abstract: We present an approach for selecting objectively informative and subjectively helpful annotations to social media posts. We draw on data from on an online environment where contributors annotate misinformation and simultaneously rate the contributions of others. Our algorithm uses a matrix-factorization (MF) based approach to identify annotations that appeal broadly across heterogeneous user group… ▽ More We present an approach for selecting objectively informative and subjectively helpful annotations to social media posts. We draw on data from on an online environment where contributors annotate misinformation and simultaneously rate the contributions of others. Our algorithm uses a matrix-factorization (MF) based approach to identify annotations that appeal broadly across heterogeneous user groups - sometimes referred to as "bridging-based ranking." We pair these data with a survey experiment in which individuals are randomly assigned to see annotations to posts. We find that annotations selected by the algorithm improve key indicators compared with overall average and crowd-generated baselines. Further, when deployed on Twitter, people who saw annotations selected through this bridging-based approach were significantly less likely to reshare social media posts than those who did not see the annotations. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Multi-messenger characterization of Mrk 501 during historically low X-ray and $γ$-ray activity

Authors: MAGIC collaboration, H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder , et al. (300 additional authors not shown)

Abstract: We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the… ▽ More We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the highest occurring at X-rays and very-high-energy (VHE) $γ$-rays. A significant correlation ($>$3$σ$) between X-rays and VHE $γ$-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between Swift-XRT and Fermi-LAT. We additionally find correlations between high-energy $γ$-rays and radio, with the radio lagging by more than 100 days, placing the $γ$-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE $γ$-rays from mid-2017 to mid-2019 with a stable VHE flux ($>$0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2-year-long low-state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock. △ Less

Submitted 5 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

Journal ref: ApJS 266 37 (2023)

Multi-wavelength study of the galactic PeVatron candidate LHAASO J2108+5157

Authors: S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini, M. I. Bernardos, J. Bernete Medrano, A. Berti, P. Bhattacharjee, N. Biederbeck , et al. (245 additional authors not shown)

Abstract: LHAASO J2108+5157 is one of the few known unidentified Ultra-High-Energy (UHE) gamma-ray sources with no Very-High-Energy (VHE) counterpart, recently discovered by the LHAASO collaboration. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good quality data. In… ▽ More LHAASO J2108+5157 is one of the few known unidentified Ultra-High-Energy (UHE) gamma-ray sources with no Very-High-Energy (VHE) counterpart, recently discovered by the LHAASO collaboration. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good quality data. In addition, we analyzed 12 years of Fermi-LAT data, to better constrain emission of its High-Energy (HE) counterpart 4FGL J2108.0+5155. We found an excess (3.7 sigma) in the LST-1 data at energies E > 3 TeV. Further analysis in the whole LST-1 energy range assuming a point-like source, resulted in a hint (2.2 sigma) of hard emission which can be described with a single power law with photon index Gamma = 1.6 +- 0.2 between 0.3 - 100 TeV. We did not find any significant extended emission which could be related to a Supernova Remnant (SNR) or Pulsar Wind Nebula (PWN) in the XMM-Newton data, which puts strong constraints on possible synchrotron emission of relativistic electrons. The LST-1 and LHAASO observations can be explained as inverse Compton-dominated leptonic emission of relativistic electrons with a cutoff energy of $100^{+70}_{-30}$ TeV. The low magnetic field in the source imposed by the X-ray upper limits on synchrotron emission is compatible with a hypothesis of a PWN or a TeV halo. The lack of a pulsar in the neighborhood of the UHE source is a challenge to the PWN/TeV-halo scenario. The UHE gamma rays can also be explained as $π^0$ decay-dominated hadronic emission due to interaction of relativistic protons with one of the two known molecular clouds in the direction of the source. The hard spectrum in the LST-1 band is compatible with protons escaping a shock around a middle-aged SNR because of their high low-energy cut-off. △ Less

Submitted 16 March, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 26 pages, 9 figures, 5 tables; Astronomy and Astrophysics (accepted on 6 March 2023)

Journal ref: A&A 673, A75 (2023)

Gamma-ray observations of MAXI J1820+070 during the 2018 outburst

Authors: H. Abe, S. Abe, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari , et al. (418 additional authors not shown)

Abstract: MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS expe… ▽ More MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA. △ Less

Submitted 6 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

Report of the Topical Group on Dark Energy and Cosmic Acceleration: Complementarity of Probes and New Facilities for Snowmass 2021

Authors: Brenna Flaugher, Vivian Miranda, David J. Schlegel, Adam J. Anderson, Felipe Andrade-Oliveira, Eric J. Baxter, Amy N. Bender, Lindsey E. Bleem, Chihway Chang, Clarence C. Chang, Thomas Y. Chen, Kyle S. Dawson, Seth W. Digel, Alex Drlica-Wagner, Simone Ferraro, Alyssa Garcia, Katrin Heitmann, Alex G. Kim, Eric V. Linder, Sayan Mandal, Rachel Mandelbaum, Phil Marshall, Joel Meyers, Laura Newburgh, Peter E. Nugent , et al. (5 additional authors not shown)

Abstract: The mechanism(s) driving the early- and late-time accelerated expansion of the Universe represent one of the most compelling mysteries in fundamental physics today. The path to understanding the causes of early- and late-time acceleration depends on fully leveraging ongoing surveys, developing and demonstrating new technologies, and constructing and operating new instruments. This report presents… ▽ More The mechanism(s) driving the early- and late-time accelerated expansion of the Universe represent one of the most compelling mysteries in fundamental physics today. The path to understanding the causes of early- and late-time acceleration depends on fully leveraging ongoing surveys, developing and demonstrating new technologies, and constructing and operating new instruments. This report presents a multi-faceted vision for the cosmic survey program in the 2030s and beyond that derives from these considerations. Cosmic surveys address a wide range of fundamental physics questions, and are thus a unique and powerful component of the HEP experimental portfolio. △ Less

Submitted 18 September, 2022; originally announced September 2022.

Comments: Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021); Topical Group Report for CF06 (Cosmic Frontier Topical Group on Dark Energy and Cosmic Acceleration: Complementarity of Probes and New Facilities

Enabling scalable clinical interpretation of ML-based phenotypes using real world data

Authors: Owen Parsons, Nathan E Barlow, Janie Baxter, Karen Paraschin, Andrea Derix, Peter Hein, Robert Dürichen

Abstract: The availability of large and deep electronic healthcare records (EHR) datasets has the potential to enable a better understanding of real-world patient journeys, and to identify novel subgroups of patients. ML-based aggregation of EHR data is mostly tool-driven, i.e., building on available or newly developed methods. However, these methods, their input requirements, and, importantly, resulting ou… ▽ More The availability of large and deep electronic healthcare records (EHR) datasets has the potential to enable a better understanding of real-world patient journeys, and to identify novel subgroups of patients. ML-based aggregation of EHR data is mostly tool-driven, i.e., building on available or newly developed methods. However, these methods, their input requirements, and, importantly, resulting output are frequently difficult to interpret, especially without in-depth data science or statistical training. This endangers the final step of analysis where an actionable and clinically meaningful interpretation is needed.This study investigates approaches to perform patient stratification analysis at scale using large EHR datasets and multiple clustering methods for clinical research. We have developed several tools to facilitate the clinical evaluation and interpretation of unsupervised patient stratification results, namely pattern screening, meta clustering, surrogate modeling, and curation. These tools can be used at different stages within the analysis. As compared to a standard analysis approach, we demonstrate the ability to condense results and optimize analysis time. In the case of meta clustering, we demonstrate that the number of patient clusters can be reduced from 72 to 3 in one example. In another stratification result, by using surrogate models, we could quickly identify that heart failure patients were stratified if blood sodium measurements were available. As this is a routine measurement performed for all patients with heart failure, this indicated a data bias. By using further cohort and feature curation, these patients and other irrelevant features could be removed to increase the clinical meaningfulness. These examples show the effectiveness of the proposed methods and we hope to encourage further research in this field. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: 27 pages, 14 figures

MSC Class: 92C50; 92C60; 68T09; 62H30 ACM Class: J.3; I.5; I.2

Weak percolation on multiplex networks with overlapping edges

Authors: G. J. Baxter, R. A. da Costa, S. N. Dorogovtsev, J. F. F. Mendes

Abstract: We solve the weak percolation problem for multiplex networks with overlapping edges. In weak percolation, a vertex belongs to a connected component if at least one of its neighbors in each of the layers is in this component. This is a weaker condition than for a mutually connected component in interdependent networks, in which any two vertices must be connected by a path within each of the layers.… ▽ More We solve the weak percolation problem for multiplex networks with overlapping edges. In weak percolation, a vertex belongs to a connected component if at least one of its neighbors in each of the layers is in this component. This is a weaker condition than for a mutually connected component in interdependent networks, in which any two vertices must be connected by a path within each of the layers. The effect of the overlaps on weak percolation turns out to be opposite to that on the giant mutually connected component. While for the giant mutually connected component, overlaps do not change the critical phenomena, our theory shows that in two layers any (nonzero) concentration of overlaps drives the weak percolation transition to the ordinary percolation universality class. In three layers, the phase diagram of the problem contains two lines -- of a continuous phase transition and of a discontinuous one -- connected in various ways depending on how the layers overlap. In the case of only doubled overlapped edges, two of the end points of these lines coincide, resulting in a tricritical point like that seen in heterogeneous $k$-core percolation. △ Less

Submitted 26 June, 2022; originally announced June 2022.

Comments: 21 pages, 8 figures

Joint analysis of DES Year 3 data and CMB lensing from SPT and Planck III: Combined cosmological constraints

Authors: T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, J. Annis, B. Ansarinejad, S. Avila, D. Bacon, E. J. Baxter, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, E. Bertin, J. Blazek, L. E. Bleem, S. Bocquet, D. Brooks, E. Buckley-Geer, D. L. Burke, H. Camacho, A. Campos, J. E. Carlstrom , et al. (146 additional authors not shown)

Abstract: We present cosmological constraints from the analysis of two-point correlation functions between galaxy positions and galaxy lensing measured in Dark Energy Survey (DES) Year 3 data and measurements of cosmic microwave background (CMB) lensing from the South Pole Telescope (SPT) and Planck. When jointly analyzing the DES-only two-point functions and the DES cross-correlations with SPT+Planck CMB l… ▽ More We present cosmological constraints from the analysis of two-point correlation functions between galaxy positions and galaxy lensing measured in Dark Energy Survey (DES) Year 3 data and measurements of cosmic microwave background (CMB) lensing from the South Pole Telescope (SPT) and Planck. When jointly analyzing the DES-only two-point functions and the DES cross-correlations with SPT+Planck CMB lensing, we find $Ω_{\rm m} = 0.344\pm 0.030$ and $S_8 \equiv σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.773\pm 0.016$, assuming $Λ$CDM. When additionally combining with measurements of the CMB lensing autospectrum, we find $Ω_{\rm m} = 0.306^{+0.018}_{-0.021}$ and $S_8 = 0.792\pm 0.012$. The high signal-to-noise of the CMB lensing cross-correlations enables several powerful consistency tests of these results, including comparisons with constraints derived from cross-correlations only, and comparisons designed to test the robustness of the galaxy lensing and clustering measurements from DES. Applying these tests to our measurements, we find no evidence of significant biases in the baseline cosmological constraints from the DES-only analyses or from the joint analyses with CMB lensing cross-correlations. However, the CMB lensing cross-correlations suggest possible problems with the correlation function measurements using alternative lens galaxy samples, in particular the redMaGiC galaxies and high-redshift MagLim galaxies, consistent with the findings of previous studies. We use the CMB lensing cross-correlations to identify directions for further investigating these problems. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 20 pages, 15 figures

Report number: FERMILAB-PUB-22-475-PPD

Searching for velocity-dependent dark matter annihilation signals from extragalactic halos

Authors: Eric J. Baxter, Jason Kumar, Aleczander D. Paul, Jack Runburg

Abstract: We consider gamma-ray signals of dark matter annihilation in extragalactic halos in the case where dark matter annihilates from a $p$-wave or $d$-wave state. In these scenarios, signals from extragalactic halos are enhanced relative to other targets, such as the Galactic Center or dwarf spheroidal galaxies, because the typical relative speed of the dark matter is larger in extragalactic halos. We… ▽ More We consider gamma-ray signals of dark matter annihilation in extragalactic halos in the case where dark matter annihilates from a $p$-wave or $d$-wave state. In these scenarios, signals from extragalactic halos are enhanced relative to other targets, such as the Galactic Center or dwarf spheroidal galaxies, because the typical relative speed of the dark matter is larger in extragalactic halos. We perform a mock data analysis of gamma rays produced by dark matter annihilation in halos detected by the Sloan Digital Sky Survey. We include a model for uncorrelated galactic and extragalactic gamma ray backgrounds, as well as a simple model for backgrounds due to astrophysical processes in the extragalactic halos detected by the survey. We find that, for models which are still allowed by other gamma ray searches, searches of extragalactic halos with the current Fermi exposure can produce evidence for dark matter annihilation, though it is difficult to distinguish the $p$-wave and $d$-wave scenarios. With a factor $10\times$ larger exposure, though, discrimination of the velocity-dependence is possible. △ Less

Submitted 15 September, 2022; v1 submitted 4 May, 2022; originally announced May 2022.

Comments: 20 pages, 5 figures; replaced to match version accepted by JCAP

Journal ref: Journal of Cosmology and Astroparticle Physics, Volume 2022, Issue 09, id.026, 21 pp

Ultrafast imaging of polariton propagation and interactions

Authors: Ding Xu, Arkajit Mandal, James M. Baxter, Shan-Wen Cheng, Inki Lee, Haowen Su, Song Liu, David R. Reichman, Milan Delor

Abstract: Semiconductor excitations can hybridize with cavity photons to form exciton-polaritons (EPs) with remarkable properties, including light-like energy flow combined with matter-like interactions. To fully harness these properties, EPs must retain ballistic, coherent transport despite matter-mediated interactions with lattice phonons. Here we develop a nonlinear momentum-resolved optical approach tha… ▽ More Semiconductor excitations can hybridize with cavity photons to form exciton-polaritons (EPs) with remarkable properties, including light-like energy flow combined with matter-like interactions. To fully harness these properties, EPs must retain ballistic, coherent transport despite matter-mediated interactions with lattice phonons. Here we develop a nonlinear momentum-resolved optical approach that directly images EPs in real space on femtosecond scales in a range of polaritonic architectures. We focus our analysis on EP propagation in layered halide perovskite microcavities. We reveal that EP-phonon interactions lead to a large renormalization of EP velocities at high excitonic fractions at room temperature. Despite these strong EP-phonon interactions, ballistic transport is maintained for up to half-exciton EPs, in agreement with quantum simulations of dynamic disorder shielding through light-matter hybridization. Above 50% excitonic character, rapid decoherence leads to diffusive transport. Our work provides a general framework to precisely balance EP coherence, velocity, and nonlinear interactions. △ Less

Submitted 15 June, 2023; v1 submitted 2 May, 2022; originally announced May 2022.

Determining the Hubble Constant without the Sound Horizon: A $3.6\%$ Constraint on $H_0$ from Galaxy Surveys, CMB Lensing and Supernovae

Authors: Oliver H. E. Philcox, Gerrit S. Farren, Blake D. Sherwin, Eric J. Baxter, Dillon J. Brout

Abstract: Many theoretical resolutions to the so-called "Hubble tension" rely on modifying the sound horizon at recombination, $r_s$, and thus the acoustic scale used as a standard ruler in the cosmic microwave background (CMB) and large scale structure (LSS) datasets. As shown in a number of recent works, these observables can also be used to compute $r_s$-independent constraints on $H_0$ by making use of… ▽ More Many theoretical resolutions to the so-called "Hubble tension" rely on modifying the sound horizon at recombination, $r_s$, and thus the acoustic scale used as a standard ruler in the cosmic microwave background (CMB) and large scale structure (LSS) datasets. As shown in a number of recent works, these observables can also be used to compute $r_s$-independent constraints on $H_0$ by making use of the horizon scale at matter-radiation equality, $k_{\rm eq}$, which has different sensitivity to high redshift physics than $r_s$. As such, $r_s$- and $k_{\rm eq}$-based measurements of $H_0$ (within a $Λ$CDM framework) may differ if there is new physics present pre-recombination. In this work, we present the tightest constraints on the latter from current data, finding $H_0=64.8^{+2.2}_{-2.5}$ at 68% CL (in $\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ units) from a combination of BOSS galaxy power spectra, Planck CMB lensing, and the newly released Pantheon+ supernova constraints, as well as physical priors on the baryon density, neutrino mass, and spectral index. The BOSS and Planck measurements have different degeneracy directions, leading to the improved combined constraints, with a bound of $H_0 = 67.1^{+2.5}_{-2.9}$ ($63.6^{+2.9}_{-3.6}$) from BOSS (Planck) alone. The results show some dependence on the neutrino mass bounds, with the constraint broadening to $H_0 = 68.0^{+2.9}_{-3.2}$ if we instead impose a weak prior on $\sum m_ν$ from terrestrial experiments, or shifting to $H_0 = 64.6\pm2.4$ if the neutrino mass is fixed to its minimal value. Even without dependence on the sound horizon, our results are in $\approx 3σ$ tension with those obtained from the Cepheid-calibrated distance ladder, which begins to cause problems for new physics models that vary $H_0$ by changing acoustic physics or the expansion history immediately prior to recombination. △ Less

Submitted 13 September, 2022; v1 submitted 6 April, 2022; originally announced April 2022.

Comments: 11 pages, 3 figures, accepted by Phys. Rev. D

Joint analysis of DES Year 3 data and CMB lensing from SPT and Planck II: Cross-correlation measurements and cosmological constraints

Authors: C. Chang, Y. Omori, E. J. Baxter, C. Doux, A. Choi, S. Pandey, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, K. Bechtol, M. R. Becker, G. M. Bernstein, F. Bianchini, J. Blazek, L. E. Bleem, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, R. Chen, J. Cordero, T. M. Crawford, M. Crocce , et al. (141 additional authors not shown)

Abstract: Cross-correlations of galaxy positions and galaxy shears with maps of gravitational lensing of the cosmic microwave background (CMB) are sensitive to the distribution of large-scale structure in the Universe. Such cross-correlations are also expected to be immune to some of the systematic effects that complicate correlation measurements internal to galaxy surveys. We present measurements and model… ▽ More Cross-correlations of galaxy positions and galaxy shears with maps of gravitational lensing of the cosmic microwave background (CMB) are sensitive to the distribution of large-scale structure in the Universe. Such cross-correlations are also expected to be immune to some of the systematic effects that complicate correlation measurements internal to galaxy surveys. We present measurements and modeling of the cross-correlations between galaxy positions and galaxy lensing measured in the first three years of data from the Dark Energy Survey with CMB lensing maps derived from a combination of data from the 2500 deg$^2$ SPT-SZ survey conducted with the South Pole Telescope and full-sky data from the Planck satellite. The CMB lensing maps used in this analysis have been constructed in a way that minimizes biases from the thermal Sunyaev Zel'dovich effect, making them well suited for cross-correlation studies. The total signal-to-noise of the cross-correlation measurements is 23.9 (25.7) when using a choice of angular scales optimized for a linear (nonlinear) galaxy bias model. We use the cross-correlation measurements to obtain constraints on cosmological parameters. For our fiducial galaxy sample, which consist of four bins of magnitude-selected galaxies, we find constraints of $Ω_{m} = 0.272^{+0.032}_{-0.052}$ and $S_{8} \equiv σ_8 \sqrt{Ω_{m}/0.3}= 0.736^{+0.032}_{-0.028}$ ($Ω_{m} = 0.245^{+0.026}_{-0.044}$ and $S_{8} = 0.734^{+0.035}_{-0.028}$) when assuming linear (nonlinear) galaxy bias in our modeling. Considering only the cross-correlation of galaxy shear with CMB lensing, we find $Ω_{m} = 0.270^{+0.043}_{-0.061}$ and $S_{8} = 0.740^{+0.034}_{-0.029}$. Our constraints on $S_8$ are consistent with recent cosmic shear measurements, but lower than the values preferred by primary CMB measurements from Planck. △ Less

Submitted 31 March, 2022; v1 submitted 23 March, 2022; originally announced March 2022.

Comments: 25 pages, 19 figures, submitted to PRD

Report number: FERMILAB-PUB-22-098-PPD

Joint analysis of DES Year 3 data and CMB lensing from SPT and Planck I: Construction of CMB Lensing Maps and Modeling Choices

Authors: Y. Omori, E. J. Baxter, C. Chang, O. Friedrich, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, K. Bechtol, M. R. Becker, G. M. Bernstein, J. Blazek, L. E. Bleem, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, R. Chen, A. Choi, J. Cordero, T. M. Crawford, M. Crocce, C. Davis, J. DeRose , et al. (138 additional authors not shown)

Abstract: Joint analyses of cross-correlations between measurements of galaxy positions, galaxy lensing, and lensing of the cosmic microwave background (CMB) offer powerful constraints on the large-scale structure of the Universe. In a forthcoming analysis, we will present cosmological constraints from the analysis of such cross-correlations measured using Year 3 data from the Dark Energy Survey (DES), and… ▽ More Joint analyses of cross-correlations between measurements of galaxy positions, galaxy lensing, and lensing of the cosmic microwave background (CMB) offer powerful constraints on the large-scale structure of the Universe. In a forthcoming analysis, we will present cosmological constraints from the analysis of such cross-correlations measured using Year 3 data from the Dark Energy Survey (DES), and CMB data from the South Pole Telescope (SPT) and Planck. Here we present two key ingredients of this analysis: (1) an improved CMB lensing map in the SPT-SZ survey footprint, and (2) the analysis methodology that will be used to extract cosmological information from the cross-correlation measurements. Relative to previous lensing maps made from the same CMB observations, we have implemented techniques to remove contamination from the thermal Sunyaev Zel'dovich effect, enabling the extraction of cosmological information from smaller angular scales of the cross-correlation measurements than in previous analyses with DES Year 1 data. We describe our model for the cross-correlations between these maps and DES data, and validate our modeling choices to demonstrate the robustness of our analysis. We then forecast the expected cosmological constraints from the galaxy survey-CMB lensing auto and cross-correlations. We find that the galaxy-CMB lensing and galaxy shear-CMB lensing correlations will on their own provide a constraint on $S_8=σ_8 \sqrt{Ω_{\rm m}/0.3}$ at the few percent level, providing a powerful consistency check for the DES-only constraints. We explore scenarios where external priors on shear calibration are removed, finding that the joint analysis of CMB lensing cross-correlations can provide constraints on the shear calibration amplitude at the 5 to 10% level. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 30 pages, 20 figures, To be submitted to PRD

Report number: FERMILAB-PUB-22-194-PPD

Snowmass2021: Opportunities from Cross-survey Analyses of Static Probes

Authors: Eric J. Baxter, Chihway Chang, Andrew Hearin, Jonathan Blazek, Lindsey E. Bleem, Simone Ferraro, Mustapha Ishak, Kirit S. Karkare, Alexie Leauthaud, Jia Liu, Rachel Mandelbaum, Joel Meyers, Azadeh Moradinezhad Dizgah, Daisuke Nagai, Jeffrey A. Newman, Yuuki Omori, Neelima Sehgal, Martin White, Joe Zuntz, Marcelo A. Alvarez, Camille Avestruz, Federico Bianchini, Sebastian Bocquet, Boris Bolliet, John E. Carlstrom , et al. (15 additional authors not shown)

Abstract: Cosmological data in the next decade will be characterized by high-precision, multi-wavelength measurements of thousands of square degrees of the same patches of sky. By performing multi-survey analyses that harness the correlated nature of these datasets, we will gain access to new science, and increase the precision and robustness of science being pursued by each individual survey. However, effe… ▽ More Cosmological data in the next decade will be characterized by high-precision, multi-wavelength measurements of thousands of square degrees of the same patches of sky. By performing multi-survey analyses that harness the correlated nature of these datasets, we will gain access to new science, and increase the precision and robustness of science being pursued by each individual survey. However, effective application of such analyses requires a qualitatively new level of investment in cross-survey infrastructure, including simulations, associated modeling, coordination of data sharing, and survey strategy. The scientific gains from this new level of investment are multiplicative, as the benefits can be reaped by even present-day instruments, and can be applied to new instruments as they come online. △ Less

Submitted 16 May, 2022; v1 submitted 13 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Effect of initial infection size on network SIR model

Authors: G. Machado, G. J. Baxter

Abstract: We consider the effect of a nonvanishing fraction of initially infected nodes (seeds) on the SIR epidemic model on random networks. This is relevant when, for example, the number of arriving infected individuals is large, but also to the modeling of a large number of infected individuals, but also to more general situations such as the spread of ideas in the presence of publicity campaigns. This m… ▽ More We consider the effect of a nonvanishing fraction of initially infected nodes (seeds) on the SIR epidemic model on random networks. This is relevant when, for example, the number of arriving infected individuals is large, but also to the modeling of a large number of infected individuals, but also to more general situations such as the spread of ideas in the presence of publicity campaigns. This model is frequently studied by mapping to a bond percolation problem, in which edges in the network are occupied with the probability, $p$, of eventual infection along an edge connecting an infected individual to a susceptible neighbor. This approach allows one to calculate the total final size of the infection and epidemic threshold in the limit of a vanishingly small seed fraction. We show, however, that when the initial infection occupies a nonvanishing fraction $f$ of the network, this method yields ambiguous results, as the correspondence between edge occupation and contagion transmission no longer holds. We propose instead to measure the giant component of recovered individuals within the original contact network. This has an unambiguous interpretation and correctly captures the dependence of the epidemic size on $f$. We give exact equations for the size of the epidemic and the epidemic threshold in the infinite size limit. We observe a second order phase transition as in the original formulation, however with an epidemic threshold which decreases with increasing $f$. When the seed fraction $f$ tends to zero we recover the standard results. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 10 pages, 6 figures

Compensating trajectory bias for unsupervised patient stratification using adversarial recurrent neural networks

Authors: Avelino Javer, Owen Parsons, Oliver Carr, Janie Baxter, Christian Diedrich, Eren Elçi, Steffen Schaper, Katrin Coboeken, Robert Dürichen

Abstract: Electronic healthcare records are an important source of information which can be used in patient stratification to discover novel disease phenotypes. However, they can be challenging to work with as data is often sparse and irregularly sampled. One approach to solve these limitations is learning dense embeddings that represent individual patient trajectories using a recurrent neural network autoe… ▽ More Electronic healthcare records are an important source of information which can be used in patient stratification to discover novel disease phenotypes. However, they can be challenging to work with as data is often sparse and irregularly sampled. One approach to solve these limitations is learning dense embeddings that represent individual patient trajectories using a recurrent neural network autoencoder (RNN-AE). This process can be susceptible to unwanted data biases. We show that patient embeddings and clusters using previously proposed RNN-AE models might be impacted by a trajectory bias, meaning that results are dominated by the amount of data contained in each patients trajectory, instead of clinically relevant details. We investigate this bias on 2 datasets (from different hospitals) and 2 disease areas as well as using different parts of the patient trajectory. Our results using 2 previously published baseline methods indicate a particularly strong bias in case of an event-to-end trajectory. We present a method that can overcome this issue using an adversarial training scheme on top of a RNN-AE. Our results show that our approach can reduce the trajectory bias in all cases. △ Less

Submitted 14 December, 2021; originally announced December 2021.