Showing 1–50 of 535 results for author: Benson, A

Design and Validation of a Cold Load for Characterization of CMB-S4 Detectors

Authors: Cesiley L. King, Ian Gullet, Adam J. Anderson, Bradford A. Benson, Rick Bihary, Haichen Fan, Johanna M. Nagy, Hogan Nguyen, John E. Ruhl, Sara M. Simon

Abstract: We present the design and validation of a variable temperature cryogenic blackbody source, hereinafter called a cold load, that will be used to characterize detectors to be deployed by CMB-S4, the next-generation ground-based cosmic microwave background (CMB) experiment. Although cold loads have been used for detector characterization by previous CMB experiments, this cold load has three novel des… ▽ More We present the design and validation of a variable temperature cryogenic blackbody source, hereinafter called a cold load, that will be used to characterize detectors to be deployed by CMB-S4, the next-generation ground-based cosmic microwave background (CMB) experiment. Although cold loads have been used for detector characterization by previous CMB experiments, this cold load has three novel design features: (1) the ability to operate from the 1 K stage of a dilution refrigerator (DR), (2) a 3He gas-gap heat switch to reduce cooling time, and (3) the ability to couple small external optical signals to measure detector optical time constants under low optical loading. The efficacy of this design was validated using a 150 GHz detector array previously deployed by the Spider experiment. Thermal tests showed that the cold load can be heated to temperatures required for characterizing CMB-S4's detectors without significantly impacting the temperatures of other cryogenic stages when mounted to the DR's 1 K stage. Additionally, optical tests demonstrated that external signals can be coupled to a detector array through the cold load without imparting a significant optical load on the detectors, which will enable measurements of the CMB-S4 detectors' optical time constants. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Advancing Stellar Streams as a Dark Matter Probe -- I: Evolution of the CDM subhalo population

Authors: Paul Menker, Andrew Benson

Abstract: Stellar streams, long thin streams of stars, have been used as sensitive probes of dark matter substructure for over two decades. Gravitational interactions between dark matter substructures and streams lead to the formation of low density "gaps" in streams, with any given stream typically containing no more than a few such gaps. Prior models for the statistics of such gaps have relied on several… ▽ More Stellar streams, long thin streams of stars, have been used as sensitive probes of dark matter substructure for over two decades. Gravitational interactions between dark matter substructures and streams lead to the formation of low density "gaps" in streams, with any given stream typically containing no more than a few such gaps. Prior models for the statistics of such gaps have relied on several simplifying assumptions for the properties of the subhalo population in the cold dark matter scenario. With the expected forthcoming increase in the number of streams, and gaps, observed, in this work we develop a more detailed model for the statistics of subhalos interacting with streams, and test some of the assumptions made in prior works. Instead of using simple fits to N-body estimates of subhalo population statistics at z = 0 as in previous work, we make use of realizations of time-dependent subhalo populations generated from a fully physical model, incorporating structure formation, and subhalo orbital evolution, including tidal heating and stripping physics, which has been carefully calibrated to match results of cosmological N-body simulations. We find that this model predicts up to 60% more gaps on average in Pal-5-like streams than prior works. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 11 pages, 7 figures, 1 table

A central limit theorem with explicit Lyapunov exponent and variance for products of $2\times2$ random non-invertible matrices

Authors: Audrey Benson, Hunter Gould, Phanuel Mariano, Grace Newcombe, Joshua Vaidman

Abstract: The theory of products of random matrices and Lyapunov exponents have been widely studied and applied in the fields of biology, dynamical systems, economics, engineering and statistical physics. We consider the product of an i.i.d. sequence of $2\times 2$ random non-invertible matrices with real entries. Given some mild moment assumptions we prove an explicit formula for the Lyapunov exponent and… ▽ More The theory of products of random matrices and Lyapunov exponents have been widely studied and applied in the fields of biology, dynamical systems, economics, engineering and statistical physics. We consider the product of an i.i.d. sequence of $2\times 2$ random non-invertible matrices with real entries. Given some mild moment assumptions we prove an explicit formula for the Lyapunov exponent and prove a central limit theorem with an explicit formula for the variance in terms of the entries of the matrices. We also give examples where exact values for the Lyapunov exponent and variance are computed. An important example where non-invertible matrices are essential is the random Hill's equation, which has numerous physical applications, including the astrophysical orbit problem. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 11 pages

JWST Lensed quasar dark matter survey II: Strongest gravitational lensing limit on the dark matter free streaming length to date

Authors: Ryan E. Keeley, Anna M. Nierenberg, Daniel Gilman, Charles Gannon, Simon Birrer, Tommaso Treu, Andrew J. Benson, Xiaolong Du, K. N. Abazajian, T. Anguita, V. N. Bennert, S. G. Djorgovski, K. K. Gupta, S. F. Hoenig, A. Kusenko, C. Lemon, M. Malkan, V. Motta, L. A. Moustakas, M. S. H. Oh, D. Sluse, D. Stern, R. H. Wechsler

Abstract: This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to pertu… ▽ More This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to perturbations by populations of halos down to masses $\sim 10^6$ M$_{\odot}$. Using these warm dust flux-ratio measurements in combination with 5 previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhalos. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6} M_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly$α$ forest and Milky Way satellite galaxies. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Surface Brightness Fluctuations in Two SPT clusters: a Pilot Study

Authors: Charles E. Romero, Massimo Gaspari, Gerrit Schellenberger, Bradford A. Benson, Lindsey E. Bleem, Esra Bulbul, Matthias Klein, Ralph Kraft, Paul Nulsen, Christian L. Reichardt, Laura Salvati, Taweewat Somboonpanyakul, Yuanyuan Su

Abstract: Studies of surface brightness fluctuations in the intracluster medium (ICM) present an indirect probe of turbulent properties such as the turbulent velocities, injection scales, and the slope of the power spectrum of fluctuations towards smaller scales. With the advancement of Sunyaev-Zel'dovich (SZ) studies and surveys relative to X-ray observations, we seek to investigate surface brightness fluc… ▽ More Studies of surface brightness fluctuations in the intracluster medium (ICM) present an indirect probe of turbulent properties such as the turbulent velocities, injection scales, and the slope of the power spectrum of fluctuations towards smaller scales. With the advancement of Sunyaev-Zel'dovich (SZ) studies and surveys relative to X-ray observations, we seek to investigate surface brightness fluctuations in a sample of SPT-SZ clusters which also have archival \textit{XMM-Newton} data. Here we present a pilot study of two typical clusters in that sample: SPT-CLJ0232-4421 and SPT-CLJ0638-5358. We infer injection scales larger than 500 kpc in both clusters and Mach numbers $\approx 0.5$ in SPT-CLJ0232-4421 and Mach numbers $\approx 0.6 - 1.6$ in SPT-CLJ0638-5358, which has a known shock. We find hydrostatic bias values for $M_{500}$ less than 0.2 for SPT-CLJ0232-4421 and less than 0.1 for SPT-CLJ0638-5358. These results show the importance to assess its quantitative values via a detailed multiwavelength approach and suggest that the drivers of turbulence may occur at quite larger scales. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: Accepted to ApJ; 34 pages, 23 figures, and 14 tables

Mass calibration of DES Year-3 clusters via SPT-3G CMB cluster lensing

Authors: B. Ansarinejad, S. Raghunathan, T. M. C. Abbott, P. A. R. Ade, M. Aguena, O. Alves, A. J. Anderson, F. Andrade-Oliveira, M. Archipley, L. Balkenhol, K. Benabed, A. N. Bender, B. A. Benson, E. Bertin, F. Bianchini, L. E. Bleem, S. Bocquet, F. R. Bouchet, D. Brooks, L. Bryant, D. L. Burke, E. Camphuis, J. E. Carlstrom, A. Carnero Rosell, J. Carretero , et al. (120 additional authors not shown)

Abstract: We measure the stacked lensing signal in the direction of galaxy clusters in the Dark Energy Survey Year 3 (DES Y3) redMaPPer sample, using cosmic microwave background (CMB) temperature data from SPT-3G, the third-generation CMB camera on the South Pole Telescope (SPT). We estimate the lensing signal using temperature maps constructed from the initial 2 years of data from the SPT-3G 'Main' survey,… ▽ More We measure the stacked lensing signal in the direction of galaxy clusters in the Dark Energy Survey Year 3 (DES Y3) redMaPPer sample, using cosmic microwave background (CMB) temperature data from SPT-3G, the third-generation CMB camera on the South Pole Telescope (SPT). We estimate the lensing signal using temperature maps constructed from the initial 2 years of data from the SPT-3G 'Main' survey, covering 1500 deg$^2$ of the Southern sky. We then use this signal as a proxy for the mean cluster mass of the DES sample. In this work, we employ three versions of the redMaPPer catalogue: a Flux-Limited sample containing 8865 clusters, a Volume-Limited sample with 5391 clusters, and a Volume&Redshift-Limited sample with 4450 clusters. For the three samples, we find the mean cluster masses to be ${M}_{200{\rm{m}}}=1.66\pm0.13$ [stat.]$\pm0.03$ [sys.], $1.97\pm0.18$ [stat.]$\pm0.05$ [sys.], and $2.11\pm0.20$ [stat.]$\pm0.05$ [sys.]$\times{10}^{14}\ {\rm{M}}_{\odot }$, respectively. This is a factor of $\sim2$ improvement relative to the precision of measurements with previous generations of SPT surveys and the most constraining cluster mass measurements using CMB cluster lensing to date. Overall, we find no significant tensions between our results and masses given by redMaPPer mass-richness scaling relations of previous works, which were calibrated using CMB cluster lensing, optical weak lensing, and velocity dispersion measurements from various combinations of DES, SDSS and Planck data. We then divide our sample into 3 redshift and 3 richness bins, finding no significant tensions with optical weak-lensing calibrated masses in these bins. We forecast a $5.7\%$ constraint on the mean cluster mass of the DES Y3 sample with the complete SPT-3G surveys when using both temperature and polarization data and including an additional $\sim1400$ deg$^2$ of observations from the 'Extended' SPT-3G survey. △ Less

Submitted 12 June, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: 23 pages, 9 figures, accepted for publication in JCAP. Minor changes and corrections have been made relative to v1

Testing the $\mathbfΛ$CDM Cosmological Model with Forthcoming Measurements of the Cosmic Microwave Background with SPT-3G

Authors: K. Prabhu, S. Raghunathan, M. Millea, G. Lynch, P. A. R. Ade, E. Anderes, A. J. Anderson, B. Ansarinejad, M. Archipley, L. Balkenhol, K. Benabed, A. N. Bender, B. A. Benson, F. Bianchini, L. E. Bleem, F. R. Bouchet, L. Bryant, E. Camphuis, J. E. Carlstrom, T. W. Cecil, C. L. Chang, P. Chaubal, P. M. Chichura, T. -L. Chou, A. Coerver , et al. (76 additional authors not shown)

Abstract: We forecast constraints on cosmological parameters enabled by three surveys conducted with SPT-3G, the third-generation camera on the South Pole Telescope. The surveys cover separate regions of 1500, 2650, and 6000 ${\rm deg}^{2}$ to different depths, in total observing 25% of the sky. These regions will be measured to white noise levels of roughly 2.5, 9, and 12 $μ{\rm K-arcmin}$, respectively, i… ▽ More We forecast constraints on cosmological parameters enabled by three surveys conducted with SPT-3G, the third-generation camera on the South Pole Telescope. The surveys cover separate regions of 1500, 2650, and 6000 ${\rm deg}^{2}$ to different depths, in total observing 25% of the sky. These regions will be measured to white noise levels of roughly 2.5, 9, and 12 $μ{\rm K-arcmin}$, respectively, in CMB temperature units at 150 GHz by the end of 2024. The survey also includes measurements at 95 and 220 GHz, which have noise levels a factor of ~1.2 and 3.5 times higher than 150 GHz, respectively, with each band having a polarization noise level ~$\sqrt{\text{2}}$ times higher than the temperature noise. We use a novel approach to obtain the covariance matrices for jointly and optimally estimated gravitational lensing potential bandpowers and unlensed CMB temperature and polarization bandpowers. We demonstrate the ability to test the $Λ{\rm CDM}$ model via the consistency of cosmological parameters constrained independently from SPT-3G and Planck data, and consider the improvement in constraints on $Λ{\rm CDM}$ extension parameters from a joint analysis of SPT-3G and Planck data. The $Λ{\rm CDM}$ cosmological parameters are typically constrained with uncertainties up to ~2 times smaller with SPT-3G data, compared to Planck, with the two data sets measuring significantly different angular scales and polarization levels, providing additional tests of the standard cosmological model. △ Less

Submitted 5 July, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: 26 pages; 13 figures; Accepted for publication in ApJ; Minor edits have been made

Connection between galaxy morphology and dark-matter halo structure I: a running threshold for thin discs and size predictors from the dark sector

Authors: Jinning Liang, Fangzhou Jiang, Houjun Mo, Andrew Benson, Avishai Dekel, Noa Tavron, Philip F. Hopkins, Luis C. Ho

Abstract: We present a series of studies on the connection between galaxy morphology and the structure of host dark-matter (DM) haloes using cosmological simulations. In this work, we introduce a new kinematic decomposition scheme that features physical identification of morphological components, enabling robust separation of thin and thick discs; and measure a wide range of halo properties, including their… ▽ More We present a series of studies on the connection between galaxy morphology and the structure of host dark-matter (DM) haloes using cosmological simulations. In this work, we introduce a new kinematic decomposition scheme that features physical identification of morphological components, enabling robust separation of thin and thick discs; and measure a wide range of halo properties, including their locations in the cosmic web, internal structures, and assembly histories. Our analysis of the TNG50 simulation reveals that the orbital-circularity threshold for disc differentiation varies across galaxies, with systematic trends in mass and redshift, so the widely used decomposition method with constant circularity cuts is oversimplified and underestimates thin disc at JWST redshifts. The energy threshold between the stellar halo and the inner galaxy is also a function of mass and redshift, minimizing at the sub-Galactic halo mass, where the circularity threshold peaks. Revisiting the issue of galaxy size predictor, we show that disc sizes in TNG50 exhibit correlations with three structural parameters besides virial mass and redshift: 1) a positive correlation with halo spin $λ$ across redshifts -- stronger than previously reported for zoom-in simulations but still weaker than the simple $r_{1/2}/R_{\rm vir} \propto λ$ scaling; 2) an anti-correlation with DM concentration $c$ that is well described by $r_{1/2}/R_{\rm vir} \propto c^{-0.7}$ even when $c$ is measured in the DM only run; 3) more actively accreting haloes having slightly larger discs, as well as more significant stellar haloes and lower thin-to-thick ratio. Disc mass fraction is higher in rounder haloes and in cosmic knots and filaments, implying that disc development needs both stable halo conditions and continuous material supply. Our methodology is public and adaptable to other simulations. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 20 pages, 17 figures, submitted to MNRAS

Tidal evolution of cored and cuspy dark matter halos

Authors: Xiaolong Du, Andrew Benson, Zhichao Carton Zeng, Tommaso Treu, Annika H. G. Peter, Charlie Mace, Fangzhou Jiang, Shengqi Yang, Charles Gannon, Daniel Gilman, Anna. M. Nierenberg, Ethan O. Nadler

Abstract: The internal structure and abundance of dark matter halos and subhalos are powerful probes of the nature of dark matter. In order to compare observations with dark matter models, accurate theoretical predictions of these quantities are needed. We present a fast and accurate method to describe the tidal evolution of subhalos within their parent halo, based on a semi-analytic approach. We first cons… ▽ More The internal structure and abundance of dark matter halos and subhalos are powerful probes of the nature of dark matter. In order to compare observations with dark matter models, accurate theoretical predictions of these quantities are needed. We present a fast and accurate method to describe the tidal evolution of subhalos within their parent halo, based on a semi-analytic approach. We first consider idealized N-body simulations of subhalos within their host halo, using a generalized mass density profile that describes their properties in a variety of dark matter models at infall, including popular warm, cold, and self-interacting ones. Using these simulations we construct tidal "tracks" for the evolution of subhalos based on their conditions at infall. Second, we use the results of these simulations to build semi-analytic models (SAMs) for tidal effects, including stripping and heating and implement them within the code GALACTICUS. Our SAMs can accurately predict the tidal evolution of both cored and cuspy subhalos, including the bound mass and density profiles, providing a powerful and efficient tool for studying the post-infall properties of subhalos in different dark matter models. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 25 pages, 50 figures

Star formation beyond galaxies: widespread in-situ formation of intra-cluster stars in TNG50

Authors: Niusha Ahvazi, Laura V. Sales, Julio F. Navarro, Andrew Benson, Alessandro Boselli, Richard D'Souza

Abstract: We study the fraction of the intra-cluster light (ICL) formed in-situ in the three most massive clusters of the TNG50 simulation, with virial masses $\sim 10^{14}$ M$_{\odot}$. We find that a significant fraction of ICL stars ($8\%$-$28\%$) are born in-situ. This amounts to a total stellar mass comparable to the central galaxy itself. Contrary to simple expectations, only a sub-dominant fraction o… ▽ More We study the fraction of the intra-cluster light (ICL) formed in-situ in the three most massive clusters of the TNG50 simulation, with virial masses $\sim 10^{14}$ M$_{\odot}$. We find that a significant fraction of ICL stars ($8\%$-$28\%$) are born in-situ. This amounts to a total stellar mass comparable to the central galaxy itself. Contrary to simple expectations, only a sub-dominant fraction of these in-situ ICL stars are born in the central regions and later re-distributed to more energetic orbits during mergers. Instead, many in-situ ICL stars form directly hundreds of kiloparsecs away from the central galaxy, in clouds condensing out of the circum-cluster medium. The simulations predict a present-date diffuse star formation rate of $\sim$1 $\mathrm{M}_{\odot}$/yr, with higher rates at higher redshifts. The diffuse star forming component of the ICL is filamentary in nature, extends for hundreds of kiloparsecs and traces the distribution of neutral gas in the cluster host halo. We discuss briefly how numerical details of the baryonic treatment in the simulation may play a role in this result and conclude that a sensitivity of $1.6 \times 10^{-19} - 2.6 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ in H$_α$ flux (beyond current observational capabilities) would be necessary to detect this diffuse star-forming component in galaxy clusters. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 10 pages, 8 figures - comments are welcome

First Constraints on the Epoch of Reionization Using the non-Gaussianity of the Kinematic Sunyaev-Zel{'}dovich Effect from the South Pole Telescope and {\it Herschel}-SPIRE Observations

Authors: S. Raghunathan, P. A. R. Ade, A. J. Anderson, B. Ansarinejad, M. Archipley, J. E. Austermann, L. Balkenhol, J. A. Beall, K. Benabed, A. N. Bender, B. A. Benson, F. Bianchini, L. E. Bleem, J. Bock, F. R. Bouchet, L. Bryant, E. Camphuis, J. E. Carlstrom, T. W. Cecil, C. L. Chang, P. Chaubal, H. C. Chiang, P. M. Chichura, T. -L. Chou, R. Citron , et al. (97 additional authors not shown)

Abstract: We report results from an analysis aimed at detecting the trispectrum of the kinematic Sunyaev-Zel{'}dovich (kSZ) effect by combining data from the South Pole Telescope (SPT) and {\it Herschel}-SPIRE experiments over a 100 ${\rm deg}^{2}$ field. The SPT observations combine data from the previous and current surveys, namely SPTpol and SPT-3G, to achieve depths of 4.5, 3, and 16 $μ{\rm K-arcmin}$ i… ▽ More We report results from an analysis aimed at detecting the trispectrum of the kinematic Sunyaev-Zel{'}dovich (kSZ) effect by combining data from the South Pole Telescope (SPT) and {\it Herschel}-SPIRE experiments over a 100 ${\rm deg}^{2}$ field. The SPT observations combine data from the previous and current surveys, namely SPTpol and SPT-3G, to achieve depths of 4.5, 3, and 16 $μ{\rm K-arcmin}$ in bands centered at 95, 150, and 220 GHz. For SPIRE, we include data from the 600 and 857 GHz bands. We reconstruct the velocity-induced large-scale correlation of the small-scale kSZ signal with a quadratic estimator that uses two cosmic microwave background (CMB) temperature maps, constructed by optimally combining data from all the frequency bands. We reject the null hypothesis of a zero trispectrum at $10.3σ$ level. However, the measured trispectrum contains contributions from both the kSZ and other undesired components, such as CMB lensing and astrophysical foregrounds, with kSZ being sub-dominant. We use the \textsc{Agora} simulations to estimate the expected signal from CMB lensing and astrophysical foregrounds. After accounting for the contributions from CMB lensing and foreground signals, we do not detect an excess kSZ-only trispectrum and use this non-detection to set constraints on reionization. By applying a prior based on observations of the Gunn-Peterson trough, we obtain an upper limit on the duration of reionization of $Δz_{\rm re, 50} < 4.5$ (95\% C.L). We find these constraints are fairly robust to foregrounds assumptions. This trispectrum measurement is independent of, but consistent with, {\it Planck}'s optical depth measurement. This result is the first constraint on the epoch of reionization using the non-Gaussian nature of the kSZ signal. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 15 pages, 5 figures (3 in main text and 2 in Appendix); To be submitted to PRL; Comments welcome; Data products and plotting scripts can be downloaded from https://github.com/sriniraghunathan/kSZ_4pt_SPT_SPIRE

Convergence Tests of Self-Interacting Dark Matter Simulations

Authors: Charlie Mace, Zhichao Carton Zeng, Annika H. G. Peter, Xiaolong Du, Shengqi Yang, Andrew Benson, Mark Vogelsberger

Abstract: Self-interacting dark matter (SIDM) theory predicts that dark matter halos experience core-collapse, a process where the halo's inner region rapidly increases in density and decreases in size. The N-body simulations used to study this process can suffer from numerical errors when simulation parameters are selected incorrectly. Optimal choices for simulation parameters are well studied for cold dar… ▽ More Self-interacting dark matter (SIDM) theory predicts that dark matter halos experience core-collapse, a process where the halo's inner region rapidly increases in density and decreases in size. The N-body simulations used to study this process can suffer from numerical errors when simulation parameters are selected incorrectly. Optimal choices for simulation parameters are well studied for cold dark matter (CDM), but are not deeply understood when self-interactions are included. In order to perform reliable N-body simulations and model core-collapse accurately we must understand the potential numerical errors, how to diagnose them, and what parameter selections must be made to reduce them. We use the \texttt{Arepo} N-body code to perform convergence tests of core-collapsing SIDM halos across a range of halo concentrations and SIDM cross-sections, and quantify potential numerical issues related to mass resolution, timestep size, and gravitational softening length. Our tests discover that halos with fewer than $10^5$ simulation particles, a resolution typically not met by subhalos in N-body simulations, suffer from significant discreteness noise that leads to variation and extreme outliers in the collapse rate. At our lowest resolution of $N=10^4$ particles, this collapse time variation can reach as high as 20\%. At this low resolution we also find a bias in collapse times and a small number of extreme outliers. Additionally, we find that simulations which run far beyond the age of the Universe, which have been used to calibrate SIDM gravothermal fluid models in previous work, have a sensitivity to the timestep size that is not present in shorter simulations or simulations using only CDM. Our work shows that choices of simulation parameters that yield converged results for some halo masses and SIDM models do not necessarily yield convergence for others. △ Less

Submitted 12 April, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: 16 pages, 11 figures; added new references, additional discussion, and new figure (fig 1)

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Flaring Stars in a Non-targeted mm-wave Survey with SPT-3G

Authors: C. Tandoi, S. Guns, A. Foster, P. A. R. Ade, A. J. Anderson, B. Ansarinejad, M. Archipley, L. Balkenhol, K. Benabed, A. N. Bender, B. A. Benson, F. Bianchini, L. E. Bleem, F. R. Bouchet, L. Bryant, E. Camphuis, J. E. Carlstrom, T. W. Cecil, C. L. Chang, P. Chaubal, P. M. Chichura, T. -L. Chou, A. Coerver, T. M. Crawford, A. Cukierman , et al. (74 additional authors not shown)

Abstract: We present a flare star catalog from four years of non-targeted millimeter-wave survey data from the South Pole Telescope (SPT). The data were taken with the SPT-3G camera and cover a 1500-square-degree region of the sky from $20^{h}40^{m}0^{s}$ to $3^{h}20^{m}0^{s}$ in right ascension and $-42^{\circ}$ to $-70^{\circ}$ in declination. This region was observed on a nearly daily cadence from 2019-2… ▽ More We present a flare star catalog from four years of non-targeted millimeter-wave survey data from the South Pole Telescope (SPT). The data were taken with the SPT-3G camera and cover a 1500-square-degree region of the sky from $20^{h}40^{m}0^{s}$ to $3^{h}20^{m}0^{s}$ in right ascension and $-42^{\circ}$ to $-70^{\circ}$ in declination. This region was observed on a nearly daily cadence from 2019-2022 and chosen to avoid the plane of the galaxy. A short-duration transient search of this survey yields 111 flaring events from 66 stars, increasing the number of both flaring events and detected flare stars by an order of magnitude from the previous SPT-3G data release. We provide cross-matching to Gaia DR3, as well as matches to X-ray point sources found in the second ROSAT all-sky survey. We have detected flaring stars across the main sequence, from early-type A stars to M dwarfs, as well as a large population of evolved stars. These stars are mostly nearby, spanning 10 to 1000 parsecs in distance. Most of the flare spectral indices are constant or gently rising as a function of frequency at 95/150/220 GHz. The timescale of these events can range from minutes to hours, and the peak $νL_ν$ luminosities range from $10^{27}$ to $10^{31}$ erg s$^{-1}$ in the SPT-3G frequency bands. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Forecasts for Galaxy Formation and Dark Matter Constraints from Dwarf Galaxy Surveys

Authors: Ethan O. Nadler, Vera Gluscevic, Trey Driskell, Risa H. Wechsler, Leonidas A. Moustakas, Andrew Benson, Yao-Yuan Mao

Abstract: The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)--mass cosmological zoom-in simulations and an empiric… ▽ More The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)--mass cosmological zoom-in simulations and an empirical galaxy--halo connection model, and assess sensitivity to galaxy formation and DM signals when marginalizing over galaxy--halo connection uncertainties. We find that a survey of all satellites around one MW-mass host can constrain a galaxy formation cutoff at peak virial masses of $M_{50}=10^8~M_{\mathrm{\odot}}$ at the $1σ$ level; however, a tail toward low $M_{50}$ prevents a $2σ$ measurement. In this scenario, combining hosts with differing bright satellite abundances significantly reduces uncertainties on $M_{50}$ at the $1σ$ level, but the $2σ$ tail toward low $M_{50}$ persists. We project that observations of one (two) complete satellite populations can constrain warm DM models with $m_{\mathrm{WDM}}\approx 10~\mathrm{keV}$ ($20~\mathrm{keV}$). Subhalo mass function (SHMF) suppression can be constrained to $\approx 70\%$, $60\%$, and $50\%$ that in cold dark matter (CDM) at peak virial masses of $10^8$, $10^9$, and $10^{10}~M_{\mathrm{\odot}}$, respectively; SHMF enhancement constraints are weaker ($\approx 20$, $4$, and $2$ times that in CDM, respectively) due to galaxy--halo connection degeneracies. These results motivate searches for faint dwarf galaxies beyond the MW and indicate that ongoing missions like Euclid and upcoming facilities including the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope will probe new galaxy formation and DM physics. △ Less

Submitted 17 May, 2024; v1 submitted 18 January, 2024; originally announced January 2024.

Comments: 38 pages, 23 figures, 2 tables. Updated to published version

Journal ref: ApJ 967, 61 (2024)

SPT Clusters with DES and HST Weak Lensing. II. Cosmological Constraints from the Abundance of Massive Halos

Authors: S. Bocquet, S. Grandis, L. E. Bleem, M. Klein, J. J. Mohr, T. Schrabback, T. M. C. Abbott, P. A. R. Ade, M. Aguena, A. Alarcon, S. Allam, S. W. Allen, O. Alves, A. Amon, A. J. Anderson, J. Annis, B. Ansarinejad, J. E. Austermann, S. Avila, D. Bacon, M. Bayliss, J. A. Beall, K. Bechtol, M. R. Becker, A. N. Bender , et al. (171 additional authors not shown)

Abstract: We present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel'dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d… ▽ More We present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel'dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d surveys, and comprises 1,005 confirmed clusters in the redshift range $0.25-1.78$ over a total sky area of 5,200 deg$^2$. We use DES Year 3 weak-lensing data for 688 clusters with redshifts $z<0.95$ and HST weak-lensing data for 39 clusters with $0.6<z<1.7$. The weak-lensing measurements enable robust mass measurements of sample clusters and allow us to empirically constrain the SZ observable--mass relation. For a flat $Λ$CDM cosmology, and marginalizing over the sum of massive neutrinos, we measure $Ω_\mathrm{m}=0.286\pm0.032$, $σ_8=0.817\pm0.026$, and the parameter combination $σ_8\,(Ω_\mathrm{m}/0.3)^{0.25}=0.805\pm0.016$. Our measurement of $S_8\equivσ_8\,\sqrt{Ω_\mathrm{m}/0.3}=0.795\pm0.029$ and the constraint from Planck CMB anisotropies (2018 TT,TE,EE+lowE) differ by $1.1σ$. In combination with that Planck dataset, we place a 95% upper limit on the sum of neutrino masses $\sum m_ν<0.18$ eV. When additionally allowing the dark energy equation of state parameter $w$ to vary, we obtain $w=-1.45\pm0.31$ from our cluster-based analysis. In combination with Planck data, we measure $w=-1.34^{+0.22}_{-0.15}$, or a $2.2σ$ difference with a cosmological constant. We use the cluster abundance to measure $σ_8$ in five redshift bins between 0.25 and 1.8, and we find the results to be consistent with structure growth as predicted by the $Λ$CDM model fit to Planck primary CMB data. △ Less

Submitted 21 June, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

Analytical strong line diagnostics and their redshift evolution

Authors: Shengqi Yang, Adam Lidz, Andrew Benson, Swathya Singh Chauhan, Aaron Smith, Hui Li

Abstract: The JWST is allowing new measurements of gas-phase metallicities in galaxies between cosmic noon and cosmic dawn through observations of multiple rest-frame optical/ultraviolet [OIII], [OII], and hydrogen Balmer series lines. The most robust approach to such measurements uses luminosity ratios between the excited auroral transition, [OIII] 4364A, and the lower [OIII] 5008A/4960A lines to determine… ▽ More The JWST is allowing new measurements of gas-phase metallicities in galaxies between cosmic noon and cosmic dawn through observations of multiple rest-frame optical/ultraviolet [OIII], [OII], and hydrogen Balmer series lines. The most robust approach to such measurements uses luminosity ratios between the excited auroral transition, [OIII] 4364A, and the lower [OIII] 5008A/4960A lines to determine the gas temperature. The ratio of the luminosities in the latter transitions to those in hydrogen Balmer series lines then yield relatively clean metallicity estimates. In the absence of detection of the, often faint, [OIII] auroral line, however, the ratios of various [OIII], [OII], [NII], and Balmer lines are instead used to determine metallicities. Here we present a refined approach for extracting metallicities from these ``strong line diagnostics''. Our method exploits empirical correlations between the temperature of OIII/OII regions and gas-phase metallicity, which lie close to theoretical expectations from thermal equilibrium calculations. We then show, from first principles, how to extract metallicities from traditional strong line diagnostics, R2, R3, R23, O3O2, and N2O2. We show that these ratios depend also on ionization correction factors, but that these can be determined self-consistently along with the metallicities. We quantify the success of our method using metallicities derived from galaxies with auroral line determinations and show that it generally works better than previous empirical approaches. The scatter in the observed line ratios and redshift evolution are largely explained by O3O2 variations. We provide publicly available routines for extracting metallicities from strong line diagnostics using our methodology. △ Less

Submitted 6 June, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 11 pages, 6 figures

Galaxy Clusters Discovered via the Thermal Sunyaev-Zel'dovich Effect in the 500-square-degree SPTpol Survey

Authors: L. E. Bleem, M. Klein, T. M. C. Abbott, P. A. R. Ade, M. Aguena, O. Alves, A. J. Anderson, F. Andrade-Oliveira, B. Ansarinejad, M. Archipley, M. L. N. Ashby, J. E. Austermann, D. Bacon, J. A. Beall, A. N. Bender, B. A. Benson, F. Bianchini, S. Bocquet, D. Brooks, D. L. Burke, M. Calzadilla, J. E. Carlstrom, A. Carnero Rosell, J. Carretero, C. L. Chang , et al. (103 additional authors not shown)

Abstract: We present a catalog of 689 galaxy cluster candidates detected at significance $ξ>4$ via their thermal Sunyaev-Zel'dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and \spitzer \ satellites, to confirm 544 of these candidates as clusters with… ▽ More We present a catalog of 689 galaxy cluster candidates detected at significance $ξ>4$ via their thermal Sunyaev-Zel'dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and \spitzer \ satellites, to confirm 544 of these candidates as clusters with $\sim94\%$ purity. The sample has an approximately redshift-independent mass threshold at redshift $z>0.25$ and spans $1.5 \times 10^{14} < M_{500c} < 9.1 \times 10^{14}$ $M_\odot/h_{70}$ \ and $0.03<z\lesssim1.6$ in mass and redshift, respectively; 21\% of the confirmed clusters are at $z>1$. We use external radio data from the Sydney University Molonglo Sky Survey (SUMSS) to estimate contamination to the SZ signal from synchrotron sources. The contamination reduces the recovered $ξ$ by a median value of 0.032, or $\sim0.8\%$ of the $ξ=4$ threshold value, and $\sim7\%$ of candidates have a predicted contamination greater than $Δξ= 1$. With the exception of a small number of systems $(<1\%)$, an analysis of clusters detected in single-frequency 95 and 150 GHz data shows no significant contamination of the SZ signal by emission from dusty or synchrotron sources. This cluster sample will be a key component in upcoming astrophysical and cosmological analyses of clusters. The SPTpol millimeter-wave maps and associated data products used to produce this sample are available at https://pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html, and the NASA LAMBDA website. An interactive sky server with the SPTpol maps and Dark Energy Survey data release 2 images is also available at NCSA https://skyviewer.ncsa.illinois.edu. △ Less

Submitted 8 February, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: Matches version accepted by OJA. 19 pages + references, 14 figures, cluster candidate table provided in Appendix. Data products available at https://pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html and an interactive sky server at https://skyviewer.ncsa.illinois.edu

Journal ref: Open Journal of Astrophysics, Volume 7, 2024

The SPT-Chandra BCG Spectroscopic Survey I: Evolution of the Entropy Threshold for Cooling and Feedback in Galaxy Clusters Over the Last 10 Gyr

Authors: Michael S. Calzadilla, Michael McDonald, Bradford A. Benson, Lindsey E. Bleem, Judith H. Croston, Megan Donahue, Alastair C. Edge, Benjamin Floyd, Gordon P. Garmire, Julie Hlavacek-Larrondo, Minh T. Huynh, Gourav Khullar, Ralph P. Kraft, Brian R. McNamara, Allison G. Noble, Charles E. Romero, Florian Ruppin, Taweewat Somboonpanyakul, G. Mark Voit

Abstract: We present a multi-wavelength study of the brightest cluster galaxies (BCGs) in a sample of the 95 most massive galaxy clusters selected from South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) survey. Our sample spans a redshift range of 0.3 < z < 1.7, and is complete with optical spectroscopy from various ground-based observatories, as well as ground and space-based imaging from optical, X-ray and… ▽ More We present a multi-wavelength study of the brightest cluster galaxies (BCGs) in a sample of the 95 most massive galaxy clusters selected from South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) survey. Our sample spans a redshift range of 0.3 < z < 1.7, and is complete with optical spectroscopy from various ground-based observatories, as well as ground and space-based imaging from optical, X-ray and radio wavebands. At z~0, previous studies have shown a strong correlation between the presence of a low-entropy cool core and the presence of star-formation and a radio-loud AGN in the central BCG. We show for the first time that a central entropy threshold for star formation persists out to z~1. The central entropy (measured in this work at a radius of 10 kpc) below which clusters harbor star-forming BCGs is found to be as low as $K_\mathrm{10 ~ kpc} = 35 \pm 4$ keV cm$^2$ at z < 0.15 and as high as $K_\mathrm{10 ~ kpc} = 52 \pm 11$ keV cm$^2$ at z~1. We find only marginal (~1$σ$) evidence for evolution in this threshold. In contrast, we do not find a similar high-z analog for an entropy threshold for feedback, but instead measure a strong evolution in the fraction of radio-loud BCGs in high-entropy cores as a function of redshift. This could imply that the cooling-feedback loop was not as tight in the past, or that some other fuel source like mergers are fueling the radio sources more often with increasing redshift, making the radio luminosity an increasingly unreliable proxy for radio jet power. We also find that our SZ-based sample is missing a small (~4%) population of the most luminous radio sources ($νL_ν > 10^{42}$ erg/s), likely due to radio contamination suppressing the SZ signal with which these clusters are detected. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 22 pages. 10 figures, 2 tables. Submitted to ApJ

LATIS: Constraints on the Galaxy-halo Connection at $z \sim 2.5$ from Galaxy-galaxy and Galaxy-Ly$α$ Clustering

Authors: Andrew B. Newman, Mahdi Qezlou, Nima Chartab, Gwen C. Rudie, Guillermo A. Blanc, Simeon Bird, Andrew J. Benson, Daniel D. Kelson, Brian C. Lemaux

Abstract: The connection between galaxies and dark matter halos is often quantified using the stellar mass-halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at $z \gtrsim 2$ can also constrain the SMHM relation via the g… ▽ More The connection between galaxies and dark matter halos is often quantified using the stellar mass-halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at $z \gtrsim 2$ can also constrain the SMHM relation via the galaxy autocorrelation function and through the cross-correlation between galaxies and Ly$α$ absorption measured in transverse sightlines; however, such studies are very few and have produced some unexpected or inconclusive results. We use $\sim$3000 spectra of $z\sim2.5$ galaxies from the Lyman-alpha Tomography IMACS Survey (LATIS) to measure the galaxy-galaxy and galaxy-Ly$α$ correlation functions in four bins of stellar mass spanning $10^{9.2} \lesssim M_* / M_{\odot} \lesssim 10^{10.5}$. Parallel analyses of the MultiDark N-body and ASTRID hydrodynamic cosmological simulations allow us to model the correlation functions, estimate covariance matrices, and infer halo masses. We find that results of the two methods are mutually consistent and are broadly in accord with standard SMHM relations. This consistency demonstrates that we are able to accurately measure and model Ly$α$ transmission fluctuations $δ_F$ in LATIS. We also show that the galaxy-Ly$α$ cross-correlation, a free byproduct of optical spectroscopic galaxy surveys at these redshifts, can constrain halo masses with similar precision to galaxy-galaxy clustering. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: Accepted to ApJ

SPT Clusters with DES and HST Weak Lensing. I. Cluster Lensing and Bayesian Population Modeling of Multi-Wavelength Cluster Datasets

Authors: S. Bocquet, S. Grandis, L. E. Bleem, M. Klein, J. J. Mohr, M. Aguena, A. Alarcon, S. Allam, S. W. Allen, O. Alves, A. Amon, B. Ansarinejad, D. Bacon, M. Bayliss, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, M. Brodwin, D. Brooks, A. Campos, R. E. A. Canning, J. E. Carlstrom, A. Carnero Rosell, M. Carrasco Kind , et al. (108 additional authors not shown)

Abstract: We present a Bayesian population modeling method to analyze the abundance of galaxy clusters identified by the South Pole Telescope (SPT) with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). We discuss and validate the modeling choices with a particular focus on a robust, weak-lensing-based mass calibrati… ▽ More We present a Bayesian population modeling method to analyze the abundance of galaxy clusters identified by the South Pole Telescope (SPT) with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). We discuss and validate the modeling choices with a particular focus on a robust, weak-lensing-based mass calibration using DES data. For the DES Year 3 data, we report a systematic uncertainty in weak-lensing mass calibration that increases from 1% at $z=0.25$ to 10% at $z=0.95$, to which we add 2% in quadrature to account for uncertainties in the impact of baryonic effects. We implement an analysis pipeline that joins the cluster abundance likelihood with a multi-observable likelihood for the Sunyaev-Zel'dovich effect, optical richness, and weak-lensing measurements for each individual cluster. We validate that our analysis pipeline can recover unbiased cosmological constraints by analyzing mocks that closely resemble the cluster sample extracted from the SPT-SZ, SPTpol ECS, and SPTpol 500d surveys and the DES Year 3 and HST-39 weak-lensing datasets. This work represents a crucial prerequisite for the subsequent cosmological analysis of the real dataset. △ Less

Submitted 21 June, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

Till the core collapses: the evolution and properties of self-interacting dark matter subhalos

Authors: Zhichao Carton Zeng, Annika H. G. Peter, Xiaolong Du, Shengqi Yang, Andrew Benson, Francis-Yan Cyr-Racine, Fangzhou Jiang, Charlie Mace, R. Benton Metcalf

Abstract: One of the hottest questions in the cosmology of self-interacting dark matter (SIDM) is whether scatterings can induce detectable core-collapse in halos by the present day. Because gravitational tides can accelerate core-collapse, the most promising targets to observe core-collapse are satellite galaxies and subhalo systems. However, simulating small subhalos is computationally intensive, especial… ▽ More One of the hottest questions in the cosmology of self-interacting dark matter (SIDM) is whether scatterings can induce detectable core-collapse in halos by the present day. Because gravitational tides can accelerate core-collapse, the most promising targets to observe core-collapse are satellite galaxies and subhalo systems. However, simulating small subhalos is computationally intensive, especially when subhalos start to core-collapse. In this work, we present a hierarchical framework for simulating a population of SIDM subhalos, which reduces the computation time to linear order in the total number of subhalos. With this method, we simulate substructure lensing systems with multiple velocity-dependent SIDM models, and show how subhalo evolution depends on the SIDM model, subhalo mass and orbits. We find that an SIDM cross section of $\gtrsim 200$ cm$^2$/g at velocity scales relevant for subhalos' internal heat transfer is needed for a significant fraction of subhalos to core-collapse in a typical lens system at redshift $z=0.5$, and that core-collapse has unique observable features in lensing. We show quantitatively that core-collapse in subhalos is typically accelerated compared to field halos, except when the SIDM cross section is non-negligible ($\gtrsim \mathcal{O}(1)$ cm$^2$/g) at subhalos' orbital velocities, in which case evaporation by the host can delay core-collapse. This suggests that substructure lensing can be used to probe velocity-dependent SIDM models, especially if line-of-sight structures (field halos) can be distinguished from lens-plane subhalos. Intriguingly, we find that core-collapse in subhalos can explain the recently reported ultra-steep density profiles of substructures found by lensing with the \emph{Hubble Space Telescope} △ Less

Submitted 4 November, 2023; v1 submitted 15 October, 2023; originally announced October 2023.

Comments: 35 pages, 19 figures, comments are welcome; added subsection Results F for v2 submission; submitted to PRD

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

JWST lensed quasar dark matter survey I: Description and First Results

Authors: A. M. Nierenberg, R. E. Keeley, D. Sluse, D. Gilman, S. Birrer, T. Treu, K. N. Abazajian, T. Anguita, A. J. Benson, V. N. Bennert, S. G. Djorgovski, X. Du, C. D. Fassnacht, S. F. Hoenig, A. Kusenko, C. Lemon, M. Malkan, V. Motta, L. A. Moustakas, D. Stern, R. H. Wechsler

Abstract: The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light… ▽ More The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light is emitted from a physical region of the quasar with an angular scale of milli-arcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10$μ$m. The JWST Mid-Infrared Instrument (MIRI) is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3% precision. Our simulations suggest that this precision makes it feasible to detect the presence of 10$^7$ M$_\odot$ dark matter halos at cosmological distances. Such halos are expected to be completely dark in Cold Dark Matter models. △ Less

Submitted 18 September, 2023; originally announced September 2023.

SPT-SZ MCMF: An extension of the SPT-SZ catalog over the DES region

Authors: M. Klein, J. J. Mohr, S. Bocquet, M. Aguena, S. W. Allen, O. Alves, B. Ansarinejad, M. L. N. Ashby, D. Bacon, M. Bayliss, B. A. Benson, L. E. Bleem, M. Brodwin, D. Brooks, E. Bulbul, D. L. Burke, R. E. A. Canning, J. E. Carlstrom, A. Carnero Rosell, J. Carretero, C. L. Chang, C. Conselice, M. Costanzi, A. T. Crites, L. N. da Costa , et al. (82 additional authors not shown)

Abstract: We present an extension to a Sunyaev-Zel'dovich Effect (SZE) selected cluster catalog based on observations from the South Pole Telescope (SPT); this catalog extends to lower signal-to-noise than the previous SPT-SZ catalog and therefore includes lower mass clusters. Optically derived redshifts, centers, richnesses and morphological parameters together with catalog contamination and completeness s… ▽ More We present an extension to a Sunyaev-Zel'dovich Effect (SZE) selected cluster catalog based on observations from the South Pole Telescope (SPT); this catalog extends to lower signal-to-noise than the previous SPT-SZ catalog and therefore includes lower mass clusters. Optically derived redshifts, centers, richnesses and morphological parameters together with catalog contamination and completeness statistics are extracted using the multi-component matched filter algorithm (MCMF) applied to the S/N>4 SPT-SZ candidate list and the Dark Energy Survey (DES) photometric galaxy catalog. The main catalog contains 811 sources above S/N=4, has 91% purity and is 95% complete with respect to the original SZE selection. It contains 50% more total clusters and twice as many clusters above z=0.8 in comparison to the original SPT-SZ sample. The MCMF algorithm allows us to define subsamples of the desired purity with traceable impact on catalog completeness. As an example, we provide two subsamples with S/N>4.25 and S/N>4.5 for which the sample contamination and cleaning-induced incompleteness are both as low as the expected Poisson noise for samples of their size. The subsample with S/N>4.5 has 98% purity and 96% completeness, and will be included in a combined SPT cluster and DES weak-lensing cosmological analysis. We measure the number of false detections in the SPT-SZ candidate list as function of S/N, finding that it follows that expected from assuming Gaussian noise, but with a lower amplitude compared to previous estimates from simulations. △ Less

Submitted 4 October, 2023; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 16 pages, 17 figures, submitted to MNRAS

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

A comprehensive model for the formation and evolution of the faintest Milky Way dwarf satellites

Authors: Niusha Ahvazi, Andrew Benson, Laura V. Sales, Ethan O. Nadler, Sachi Weerasooriya, Xiaolong Du, Mia Sauda Bovill

Abstract: In this study, we modify the semi-analytic model Galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass-metallicity relation for dwarf galaxies requires us to include H$_2$ cooling, an updated UV background radiation model, and to introduce a model for the… ▽ More In this study, we modify the semi-analytic model Galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass-metallicity relation for dwarf galaxies requires us to include H$_2$ cooling, an updated UV background radiation model, and to introduce a model for the metal content of the intergalactic medium. By fine-tuning various model parameters and incorporating empirical constraints, we have tailored the model to match the statistical properties of Milky Way dwarf galaxies, such as their luminosity function and size$-$mass relation. We have validated our modified semi-analytic framework by undertaking a comparative analysis of the resulting galaxy-halo connection. We predict a total of $300 ^{+75} _{-99}$ satellites with an absolute $V$-band magnitude (M$_{V}$) less than $0$ within $300$ kpc from our Milky Way-analogs. The fraction of subhalos that host a galaxy at least this bright drops to $50\%$ by a halo peak mass of $\sim 8.9 \times 10^{7}$ M$_{\odot}$, consistent with the occupation fraction inferred from the latest observations of Milky Way satellite population. △ Less

Submitted 18 March, 2024; v1 submitted 25 August, 2023; originally announced August 2023.

Comments: 20 pages, 11 figures. Accepted for publication in MNRAS

Devouring The Centaurus A Satellites: Modeling Dwarf Galaxies with Galacticus

Authors: Sachi Weerasooriya, Mia Sauda Bovill, Matthew A. Taylor, Andrew J. Benson, Cameron Leahy

Abstract: For the first time, systematic studies of dwarf galaxies are being conducted throughout the Local Volume, including Centaurus A (NGC 5128), which is the nearest elliptical galaxy. Given Centaurus As mass (roughly ten times that of the Milky Way), AGN activity, and recent major mergers, investigating these dwarfs and their star formation physics is imperative. However, simulating the faintest dwarf… ▽ More For the first time, systematic studies of dwarf galaxies are being conducted throughout the Local Volume, including Centaurus A (NGC 5128), which is the nearest elliptical galaxy. Given Centaurus As mass (roughly ten times that of the Milky Way), AGN activity, and recent major mergers, investigating these dwarfs and their star formation physics is imperative. However, simulating the faintest dwarfs in a massive galaxy like Centaurus A with sufficient resolution in a hydrodynamic simulation is computationally expensive and currently unfeasible. In this study, we seek to reproduce Centaurus A dwarfs using the same star formation physics as the Milky Way. We employ the semi-analytic model Galacticus to model dwarfs within a 600 kpc region. Utilizing astrophysical prescriptions and parameters matching the Milky Way satellites, we explore predictions for various properties and star formation histories (SFHs) to investigate environmental effects. We also reproduce cumulative luminosity and luminosity metallicity relations consistent with observations for the overall Centaurus A satellite population, while predicting half-light radii, velocity dispersion, and SFHs for the dwarf galaxies in Centaurus A. The agreement between our predicted SFHs for Centaurus A dwarfs and those of the Milky Way implies the presence of universal processes governing star formation in these galaxies. Overall, our findings shed light on the star formation physics of dwarf galaxies in the Centaurus A system, revealing insights into their properties and dependence on the host environment. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 12 pages, 8 figures

The progenitors of the intra-cluster light and intra-cluster globular clusters in galaxy groups and clusters

Authors: Niusha Ahvazi, Laura V. Sales, Jessica E. Doppel, Andrew Benson, Richard D'Souza, Vicente Rodriguez-Gomez

Abstract: We use the IllustrisTNG50 cosmological hydrodynamical simulation, complemented by a catalog of tagged globular clusters, to investigate the properties and build up of two extended luminous components: the intra-cluster light (ICL) and the intra-cluster globular clusters (ICGC). We select the 39 most massive groups and clusters in the box, spanning the range of virial masses… ▽ More We use the IllustrisTNG50 cosmological hydrodynamical simulation, complemented by a catalog of tagged globular clusters, to investigate the properties and build up of two extended luminous components: the intra-cluster light (ICL) and the intra-cluster globular clusters (ICGC). We select the 39 most massive groups and clusters in the box, spanning the range of virial masses $5 \times 10^{12} < \rm M_{200}/\rm M_{\odot} < 2 \times 10^{14}$. We find good agreement between predictions from the simulations and current observational estimates of the fraction of mass in the ICL and its radial extension. The stellar mass of the ICL is only $\sim10\%-20\%$ of the stellar mass in the central galaxy but encodes useful information on the assembly history of the group or cluster. About half the ICL in all our systems is brought in by galaxies in a narrow stellar mass range, $M_*=10^{10}-10^{11}$ $\rm M_{\odot}$. However, the contribution of low-mass galaxies ($M_*<10^{10}$ $\rm M_{\odot}$) to the build-up of the ICL varies broadly from system to system, $\sim 5\%-45\%$, a feature that might be recovered from the observable properties of the ICL at $z=0$. At fixed virial mass, systems where the accretion of dwarf galaxies plays an important role have shallower metallicity profiles, less metal content and a lower stellar mass in the ICL than systems where the main contributors are more massive galaxies. We show that intra-cluster GCs are also good tracers of this history, representing a valuable alternative when diffuse light is not detectable. △ Less

Submitted 23 August, 2023; originally announced August 2023.

A Measurement of Gravitational Lensing of the Cosmic Microwave Background Using SPT-3G 2018 Data

Authors: Z. Pan, F. Bianchini, W. L. K. Wu, P. A. R. Ade, Z. Ahmed, E. Anderes, A. J. Anderson, B. Ansarinejad, M. Archipley, K. Aylor, L. Balkenhol, P. S. Barry, R. Basu Thakur, K. Benabed, A. N. Bender, B. A. Benson, L. E. Bleem, F. R. Bouchet, L. Bryant, K. Byrum, E. Camphuis, J. E. Carlstrom, F. W. Carter, T. W. Cecil, C. L. Chang , et al. (111 additional authors not shown)

Abstract: We present a measurement of gravitational lensing over 1500 deg$^2$ of the Southern sky using SPT-3G temperature data at 95 and 150 GHz taken in 2018. The lensing amplitude relative to a fiducial Planck 2018 $Λ$CDM cosmology is found to be $1.020\pm0.060$, excluding instrumental and astrophysical systematic uncertainties. We conduct extensive systematic and null tests to check the robustness of th… ▽ More We present a measurement of gravitational lensing over 1500 deg$^2$ of the Southern sky using SPT-3G temperature data at 95 and 150 GHz taken in 2018. The lensing amplitude relative to a fiducial Planck 2018 $Λ$CDM cosmology is found to be $1.020\pm0.060$, excluding instrumental and astrophysical systematic uncertainties. We conduct extensive systematic and null tests to check the robustness of the lensing measurements, and report a minimum-variance combined lensing power spectrum over angular multipoles of $50<L<2000$, which we use to constrain cosmological models. When analyzed alone and jointly with primary cosmic microwave background (CMB) spectra within the $Λ$CDM model, our lensing amplitude measurements are consistent with measurements from SPT-SZ, SPTpol, ACT, and Planck. Incorporating loose priors on the baryon density and other parameters including uncertainties on a foreground bias template, we obtain a $1σ$ constraint on $σ_8 Ω_{\rm m}^{0.25}=0.595 \pm 0.026$ using the SPT-3G 2018 lensing data alone, where $σ_8$ is a common measure of the amplitude of structure today and $Ω_{\rm m}$ is the matter density parameter. Combining SPT-3G 2018 lensing measurements with baryon acoustic oscillation (BAO) data, we derive parameter constraints of $σ_8 = 0.810 \pm 0.033$, $S_8 \equiv σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.836 \pm 0.039$, and Hubble constant $H_0 =68.8^{+1.3}_{-1.6}$ km s$^{-1}$ Mpc$^{-1}$. Using CMB anisotropy and lensing measurements from SPT-3G only, we provide independent constraints on the spatial curvature of $Ω_{K} = 0.014^{+0.023}_{-0.026}$ (95% C.L.) and the dark energy density of $Ω_Λ= 0.722^{+0.031}_{-0.026}$ (68% C.L.). When combining SPT-3G lensing data with SPT-3G CMB anisotropy and BAO data, we find an upper limit on the sum of the neutrino masses of $\sum m_ν< 0.30$ eV (95% C.L.). △ Less

Submitted 29 January, 2024; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: Bandpower and likelihood data available at https://pole.uchicago.edu/public/data/spt3g_2018_lensing/

Journal ref: Physical Review D 108.12 (2023): 122005

Do assumptions about the central density of subhaloes affect dark matter annihilation and lensing calculations?

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: A growing body of evidence suggests that the central density of cuspy dark matter subhaloes is conserved in minor mergers. However, empirical models of subhalo evolution, calibrated from simulations, often assume a drop in the central density. Since empirical models of subhaloes are used in galaxy-galaxy lensing studies and dark matter annihilation calculations, we explore the consequences of assu… ▽ More A growing body of evidence suggests that the central density of cuspy dark matter subhaloes is conserved in minor mergers. However, empirical models of subhalo evolution, calibrated from simulations, often assume a drop in the central density. Since empirical models of subhaloes are used in galaxy-galaxy lensing studies and dark matter annihilation calculations, we explore the consequences of assuming different subhalo models. We find that dark matter annihilation calculations are very sensitive to the assumed subhalo mass profile, and different models can give more than a magnitude difference in the J-factor and boost factor in individual haloes. On the other hand, the shear and convergence profiles used in galaxy-galaxy lensing are sensitive to the initial profile assumed (e.g., NFW versus Einato) but are otherwise well-approximated by a simple model in which the original profile is sharply truncated. We conclude that since the innermost parts of haloes are difficult to resolve in simulations, it is important to have a theoretical understanding of how subhaloes evolve to make accurate predictions of the dark matter annihilation signal. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: 15 pages, 13 figures. Submitted to MNRAS

Anisotropic strong lensing as a probe of dark matter self-interactions

Authors: Birendra Dhanasingham, Francis-Yan Cyr-Racine, Charlie Mace, Annika H. G. Peter, Andrew Benson

Abstract: Galaxy-scale strongly lensed systems have been shown to provide a unique technique for exploring the underlying physics of dark matter at sub-galactic scales. In the past, much attention was given to detecting and studying individual haloes in a strong lens system. In addition to the subhaloes, line-of-sight haloes contribute significantly to the small perturbations in lensed images. In prior work… ▽ More Galaxy-scale strongly lensed systems have been shown to provide a unique technique for exploring the underlying physics of dark matter at sub-galactic scales. In the past, much attention was given to detecting and studying individual haloes in a strong lens system. In addition to the subhaloes, line-of-sight haloes contribute significantly to the small perturbations in lensed images. In prior work, we demonstrated that these line-of-sight haloes imprint a distinctive anisotropic signature and hence give rise to a detectable non-zero parity-even quadrupole moment in the effective convergence field's two-point correlation function. In this study, we show that these line-of-sight haloes also produce a non-zero curl component of the effective deflection field with a parity-odd quadrupole moment of the two-point function. These multipole moments have the ability to statistically separate line-of-sight haloes from dark matter substructure. In this paper, we examine how these multipole moments evolve in the presence of warm dark matter and self-interacting dark matter in terms of central density evolution and dark matter halo abundance. Importantly, we show that these different multipole moments display exquisite sensitivity to both the amplitude and the velocity dependence of the dark matter self-interaction cross-section. Our approach opens the door for strong lensing observations to probe dark matter self-interaction over a broad range of relative velocities. △ Less

Submitted 2 November, 2023; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 13 pages, 9 figures + appendices. Matches the published version in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2023, Volume 526, Issue 4, pp.5455-5470

A quantitative comparison between velocity dependent SIDM cross sections constrained by the gravothermal and isothermal models

Authors: Shengqi Yang, Fangzhou Jiang, Andrew Benson, Yi-Ming Zhong, Charlie Mace, Xiaolong Du, Zhichao Carton Zeng, Annika H. G. Peter, Moritz S. Fischer

Abstract: One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest -- the effects of these proc… ▽ More One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest -- the effects of these processes can be highly degenerate with those of dark matter self-interactions at small scales. In this work we select 9 isolated galaxies and brightest cluster galaxies (BCGs) with baryonic components small enough such that the baryonic gravitational potentials do not significantly influence the halo gravothermal evolution processes. We then constrain the parameters of Rutherford and Moller scattering cross section models with the measured rotation curves and stellar kinematics through the gravothermal fluid formalism and isothermal method. Cross sections constrained by the two methods are consistent at $1σ$ confidence level, but the isothermal method prefers cross sections greater than the gravothermal approach constraints by a factor of $\sim3$. △ Less

Submitted 26 April, 2024; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 15 pages, 11 figures

Constrain the Dark Matter Distribution of Ultra-diffuse Galaxies with Globular-Cluster Mass Segregation: A Case Study with NGC5846-UDG1

Authors: Jinning Liang, Fangzhou Jiang, Shany Danieli, Andrew Benson, Phil Hopkins

Abstract: The properties of globular clusters (GCs) contain valuable information of their host galaxies and dark-matter halos. In the remarkable example of ultra-diffuse galaxy, NGC5846-UDG1, the GC population exhibits strong radial mass segregation, indicative of dynamical-friction-driven orbital decay, which opens the possibility of using imaging data alone to constrain the dark-matter content of the gala… ▽ More The properties of globular clusters (GCs) contain valuable information of their host galaxies and dark-matter halos. In the remarkable example of ultra-diffuse galaxy, NGC5846-UDG1, the GC population exhibits strong radial mass segregation, indicative of dynamical-friction-driven orbital decay, which opens the possibility of using imaging data alone to constrain the dark-matter content of the galaxy. To explore this possibility, we develop a semi-analytical model of GC evolution, which starts from the initial mass function, the initial structure-mass relation, and the initial spatial distribution of the GC progenitors, and follows the effects of dynamical friction, tidal evolution, and two-body relaxation. Using Markov Chain Monte Carlo, we forward-model the GCs in a NGC5846-UDG1-like potential to match the observed GC mass, size, and spatial distributions, and to constrain the profile of the host halo and the origin of the GCs. We find that, with the assumptions of zero mass segregation when the star clusters were born, NGC5846-UDG1 is relatively dark-matter poor compared to what is expected from stellar-to-halo-mass relations, and its halo concentration is lower than the cosmological average, irrespective of having a cuspy or a cored profile. Its GC population has an initial spatial distribution more extended than the smooth stellar distribution. We discuss the results in the context of scaling laws of galaxy-halo connections, and warn against naively using the GC-abundance-halo-mass relation to infer the halo mass of UDGs. Our model is generally applicable to GC-rich dwarf galaxies, and is publicly available at https://github.com/JiangFangzhou/GCevo. △ Less

Submitted 9 January, 2024; v1 submitted 27 April, 2023; originally announced April 2023.

Comments: 27 pages, 15 figures, ApJ accepted

Efficient simulations of ionized ISM emission lines: A detailed comparison between the FIRE high-redshift suite and observations

Authors: Shengqi Yang, Adam Lidz, Aaron Smith, Andrew Benson, Hui Li

Abstract: The Atacama Large Millimeter/Submillimeter Array (ALMA) in the sub-millimeter and the James Webb Space Telescope (JWST) in the infrared have achieved robust spectroscopic detections of emission lines from the interstellar medium (ISM) in some of the first galaxies. These unprecedented measurements provide valuable information regarding the ISM properties, stellar populations, galaxy morphologies,… ▽ More The Atacama Large Millimeter/Submillimeter Array (ALMA) in the sub-millimeter and the James Webb Space Telescope (JWST) in the infrared have achieved robust spectroscopic detections of emission lines from the interstellar medium (ISM) in some of the first galaxies. These unprecedented measurements provide valuable information regarding the ISM properties, stellar populations, galaxy morphologies, and kinematics in these high-redshift galaxies and, in principle, offer powerful tests of state-of-the-art galaxy formation models, as implemented in hydrodynamical simulations. To facilitate direct comparisons between simulations and observations, we develop a fast post-processing pipeline for predicting the line emission from the HII regions around simulated star particles, accounting for spatial variations in the surrounding gas density, metallicity, and incident radiation spectrum. Our ISM line emission model currently captures H$α$, H$β$, and all of the [OIII] and [OII] lines targeted by ALMA and the JWST at $z>6$. We illustrate the power of this approach by applying our line emission model to the publicly available Feedback In Realistic Environment (FIRE) high-$z$ simulation suite and perform a detailed comparison with current observations. We show that the FIRE mass--metallicity relation is in $1σ$ agreement with ALMA/JWST measurements after accounting for the inhomogeneities in ISM properties. We also quantitatively validate the one-zone model description, which is widely used for interpreting [OIII] and H$β$ line luminosity measurements. This model is publicly available and can be implemented on top of a broad range of galaxy formation simulations for comparison with JWST and ALMA measurements. △ Less

Submitted 4 August, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 16 pages, 15 figures. (Update gas temperature model in Section 3. Test one-zone model in more realistic JWST and ALMA mock measurements in Section 4.)

Detecting Anomalous Microflows in IoT Volumetric Attacks via Dynamic Monitoring of MUD Activity

Authors: Ayyoob Hamza, Hassan Habibi Gharakheili, Theophilus A. Benson, Gustavo Batista, Vijay Sivaraman

Abstract: IoT networks are increasingly becoming target of sophisticated new cyber-attacks. Anomaly-based detection methods are promising in finding new attacks, but there are certain practical challenges like false-positive alarms, hard to explain, and difficult to scale cost-effectively. The IETF recent standard called Manufacturer Usage Description (MUD) seems promising to limit the attack surface on IoT… ▽ More IoT networks are increasingly becoming target of sophisticated new cyber-attacks. Anomaly-based detection methods are promising in finding new attacks, but there are certain practical challenges like false-positive alarms, hard to explain, and difficult to scale cost-effectively. The IETF recent standard called Manufacturer Usage Description (MUD) seems promising to limit the attack surface on IoT devices by formally specifying their intended network behavior. In this paper, we use SDN to enforce and monitor the expected behaviors of each IoT device, and train one-class classifier models to detect volumetric attacks. Our specific contributions are fourfold. (1) We develop a multi-level inferencing model to dynamically detect anomalous patterns in network activity of MUD-compliant traffic flows via SDN telemetry, followed by packet inspection of anomalous flows. This provides enhanced fine-grained visibility into distributed and direct attacks, allowing us to precisely isolate volumetric attacks with microflow (5-tuple) resolution. (2) We collect traffic traces (benign and a variety of volumetric attacks) from network behavior of IoT devices in our lab, generate labeled datasets, and make them available to the public. (3) We prototype a full working system (modules are released as open-source), demonstrates its efficacy in detecting volumetric attacks on several consumer IoT devices with high accuracy while maintaining low false positives, and provides insights into cost and performance of our system. (4) We demonstrate how our models scale in environments with a large number of connected IoTs (with datasets collected from a network of IP cameras in our university campus) by considering various training strategies (per device unit versus per device type), and balancing the accuracy of prediction against the cost of models in terms of size and training time. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: 18 pages, 13 figures

ACM Class: C.2

Gravitational Bose-Einstein Condensation of Vector/Hidden Photon Dark Matter

Authors: Jiajun Chen, Xiaolong Du, Mingzhen Zhou, Andrew Benson, David J. E. Marsh

Abstract: We study the gravitational Bose-Einstein condensation of a massive vector field in the kinetic regime and the non-relativistic limit using non-linear dynamical numerical methods. Gravitational condensation leads to the spontaneous formation of solitons. We measure the condensation time and growth rate, and compare to analytical models in analogy to the scalar case. We find that the condensation ti… ▽ More We study the gravitational Bose-Einstein condensation of a massive vector field in the kinetic regime and the non-relativistic limit using non-linear dynamical numerical methods. Gravitational condensation leads to the spontaneous formation of solitons. We measure the condensation time and growth rate, and compare to analytical models in analogy to the scalar case. We find that the condensation time of the vector field depends on the correlation between its different components. For fully correlated configurations, the condensation time is the same as that for a scalar field. On the other hand, uncorrelated or partially correlated configurations condense slower than the scalar case. As the vector soliton grows, it can acquire a net spin angular momentum even if the total spin angular momentum of the initial conditions is zero. △ Less

Submitted 25 May, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 10 pages, 11 figures

Characterization of MKIDs for CMB observation at 220 GHz with the South Pole Telescope

Authors: Karia R. Dibert, Peter S. Barry, Adam J. Anderson, Bradford A. Benson, Thomas Cecil, Clarence L. Chang, Kyra N. Fichman, Kirit Karkare, Juliang Li, Tyler Natoli, Zhaodi Pan, Maclean Rouble, Erik Shirokoff, Matthew Young

Abstract: We present an updated design of the 220 GHz microwave kinetic inductance detector (MKID) pixel for SPT-3G+, the next-generation camera for the South Pole Telescope. We show results of the dark testing of a 63-pixel array with mean inductor quality factor $Q_i = 4.8 \times 10^5$, aluminum inductor transition temperature $T_c = 1.19$ K, and kinetic inductance fraction $α_k = 0.32$. We optically char… ▽ More We present an updated design of the 220 GHz microwave kinetic inductance detector (MKID) pixel for SPT-3G+, the next-generation camera for the South Pole Telescope. We show results of the dark testing of a 63-pixel array with mean inductor quality factor $Q_i = 4.8 \times 10^5$, aluminum inductor transition temperature $T_c = 1.19$ K, and kinetic inductance fraction $α_k = 0.32$. We optically characterize both the microstrip-coupled and CPW-coupled resonators, and find both have a spectral response close to prediction with an optical efficiency of $η\sim 70\%$. However, we find slightly lower optical response on the lower edge of the band than predicted, with neighboring dark detectors showing more response in this region, though at level consistent with less than 5\% frequency shift relative to the optical detectors. The detectors show polarized response consistent with expectations, with a cross-polar response of $\sim 10\%$ for both detector orientations. △ Less

Submitted 3 April, 2023; originally announced April 2023.

Comments: 6 pages, 5 figures, ASC 2022 proceedings

Noise Optimization for MKIDs with Different Design Geometries and Material Selections

Authors: Z. Pan, K. R. Dibert, J. Zhang, P. S. Barry, A. J. Anderson, A. N. Bender, B. A. Benson, T. Cecil, C. L. Chang, R. Gualtieri, J. Li, M. Lisovenko, V. Novosad, M. Rouble, G. Wang, V. Yefremenko

Abstract: The separation and optimization of noise components is critical to microwave-kinetic inductance detector (MKID) development. We analyze the effect of several changes to the lumped-element inductor and interdigitated capacitor geometry on the noise performance of a series of MKIDs intended for millimeter-wavelength experiments. We extract the contributions from two-level system noise in the dielect… ▽ More The separation and optimization of noise components is critical to microwave-kinetic inductance detector (MKID) development. We analyze the effect of several changes to the lumped-element inductor and interdigitated capacitor geometry on the noise performance of a series of MKIDs intended for millimeter-wavelength experiments. We extract the contributions from two-level system noise in the dielectric layer, the generation-recombination noise intrinsic to the superconducting thin-film, and system white noise from each detector noise power spectrum and characterize how these noise components depend on detector geometry, material, and measurement conditions such as driving power and temperature. We observe a reduction in the amplitude of two-level system noise with both an elevated sample temperature and an increased gap between the fingers within the interdigitated capacitors for both aluminum and niobium detectors. We also verify the expected reduction of the generation-recombination noise and associated quasiparticle lifetime with reduced inductor volume. This study also iterates over different materials, including aluminum, niobium, and aluminum manganese, and compares the results with an underlying physical model. △ Less

Submitted 3 April, 2023; originally announced April 2023.

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Characterizing the ELG luminosity functions in the nearby Universe

Authors: G. Favole, V. Gonzalez-Perez, Y. Ascasibar, P. Corcho-Caballero, A. D. Montero-Dorta, A. J. Benson, J. Comparat, S. A. Cora, D. Croton, H. Guo, D. Izquierdo-Villalba, A. Knebe, Á. Orsi, D. Stoppacher, C. A. Vega-Martínez

Abstract: Nebular emission lines are powerful diagnostics for the physical processes at play in galaxy formation and evolution. Moreover, emission-line galaxies (ELGs) are one of the main targets of current and forthcoming spectroscopic cosmological surveys. We investigate the contributions to the line luminosity functions (LFs) of different galaxy populations in the local Universe, providing a benchmark fo… ▽ More Nebular emission lines are powerful diagnostics for the physical processes at play in galaxy formation and evolution. Moreover, emission-line galaxies (ELGs) are one of the main targets of current and forthcoming spectroscopic cosmological surveys. We investigate the contributions to the line luminosity functions (LFs) of different galaxy populations in the local Universe, providing a benchmark for future surveys of earlier cosmic epochs. The large statistics of the observations from the SDSS DR7 Main galaxy sample and the MPA-JHU spectral catalogue enabled us to precisely measure the H$α$, H$β$, [OII], [OIII], and, for the first time, the [NII], and [SII] emission-line LFs over ~2.4 Gyrs in the low-z Universe, 0.02<z<0.22. We present a generalized 1/Vmax LF estimator capable of simultaneously correcting for spectroscopic, r-band magnitude, and emission-line incompleteness. We studied the contribution to the LF of different types of ELGs classified using two methods: (i) the value of the specific star formation rate (sSFR), and (ii) the line ratios on the Baldwin-Phillips-Terlevich (BPT) and the WHAN (i.e., H$α$ equivalent width versus the [NII]/H$α$ line ratio) diagrams. The ELGs in our sample are mostly star forming, with 84 per cent having sSFR>10$^{-11}$/yr. When classifying ELGs using the BPT+WHAN diagrams, we find that 63.3 per cent are star forming, only 0.03 are passively evolving, and 1.3 have nuclear activity (Seyfert). The rest are low-ionization narrow emission-line regions (LINERs) and composite ELGs. We found that a Saunders function is the most appropriate to describe all of the emission-line LFs, both observed and dust-extinction-corrected. They are dominated by star-forming regions, except for the bright end of the [OIII] and [NII] LFs (i.e. L[NII]>10$^{42}$ erg/s, L[OIII]>10$^{43}$ erg/s), where the contribution of Seyfert galaxies is not negligible. △ Less

Submitted 14 December, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Non-Gaussian assembly bias from a semi-analytic galaxy formation model

Authors: M. Marinucci, V. Desjacques, A. Benson

Abstract: We use $z=1$ mock galaxy catalogues produced with the semi-analytic code GALACTICUS to study the dependence of the non-Gaussian bias parameter $b_φ$ on the mass assembly history of the host halos. We generate large sets of merger trees and measure the non-Gaussian assembly bias $Δb_φ$ for galaxies selected by color magnitude and emission line luminosities. For galaxies selected by $g-r$ color, we… ▽ More We use $z=1$ mock galaxy catalogues produced with the semi-analytic code GALACTICUS to study the dependence of the non-Gaussian bias parameter $b_φ$ on the mass assembly history of the host halos. We generate large sets of merger trees and measure the non-Gaussian assembly bias $Δb_φ$ for galaxies selected by color magnitude and emission line luminosities. For galaxies selected by $g-r$ color, we find a large assembly bias consistent with the analysis of Barreira et al. (2020) based on hydro-dynamical simulations of galaxy formation. This effect arises from the fact that a larger value of the normalization amplitude $σ_8$ implies a faster mass assembly (at fixed halo mass) and, therefore, older and redder galaxies. On the contrary, for galaxies selected by their H$α$ luminosity, we do not detect a significant assembly bias, at least at $z=1$ and in the halo mass range $3\times10^{10} < M < 10^{12}\ M_\odot$ considered here. This is presumably due to the fact that emission line strengths are mainly sensitive to the instantaneous star formation rate, which appears to depend weakly on $σ_8$ at $z=1$. This indicates that the non-Gaussian assembly bias should be less of a concern for future emission line galaxy surveys. We also investigate, for the first time, the sensitivity of the non-Gaussian assembly bias to a change in the parameters of the galaxy formation model that control the AGN and stellar feedback as well as the star formation rate. When these parameters change within a factor of two from their fiducial value, they induce variations up to order unity in the measured $Δb_φ$, but the overall trends with color or luminosity remain the same. However, since these results may be sensitive to the choice of galaxy formation model, it will be prudent to extend this analysis to other semi-analytic models in addition to halo mass and redshift. △ Less

Submitted 18 March, 2023; originally announced March 2023.

Comments: 11 pages. Comments are welcome

Simultaneous Millimeter-wave, Gamma-ray, and Optical Monitoring of the Blazar PKS 2326-502 During a Flaring State

Authors: J. C. Hood II, A. Simpson, A. McDaniel, A. Foster, P. A. R. Ade, M. Ajello, A. J. Anderson, J. E. Austermann, J. A. Beall, A. N. Bender, B. A. Benson, F. Bianchini, L. E. Bleem, J. E. Carlstrom, C. L. Chang, P. Chaubal, H. C. Chiang, T-L. Chou, R. Citron, C. Corbett Moran, T. M. Crawford, A. T. Crites, T. de Haan, M. A. Dobbs, W. Everett , et al. (44 additional authors not shown)

Abstract: Including millimeter-wave (mm-wave) data in multi-wavelength studies of the variability of active galactic nuclei (AGN) can provide insights into AGN physics that are not easily accessible at other wavelengths. We demonstrate in this work the potential of cosmic microwave background (CMB) telescopes to provide long-term, high-cadence mm-wave AGN monitoring over large fractions of sky. We report on… ▽ More Including millimeter-wave (mm-wave) data in multi-wavelength studies of the variability of active galactic nuclei (AGN) can provide insights into AGN physics that are not easily accessible at other wavelengths. We demonstrate in this work the potential of cosmic microwave background (CMB) telescopes to provide long-term, high-cadence mm-wave AGN monitoring over large fractions of sky. We report on a pilot study using data from the SPTpol instrument on the South Pole Telescope (SPT), which was designed to observe the CMB at arcminute and larger angular scales. Between 2013 and 2016, SPTpol was used primarily to observe a single 500 deg^2 field, covering the entire field several times per day with detectors sensitive to radiation in bands centered at 95 and 150 GHz. We use SPT 150 GHz observations to create AGN light curves, and we compare these mm-wave light curves to those at other wavelengths, in particular gamma-ray and optical. In this Letter, we focus on a single source, PKS 2326-502, which has extensive, day-timescale monitoring data in gamma-ray, optical, and now mm-wave between 2013 and 2016. We find PKS 2326-502 to be in a flaring state in the first two years of this monitoring, and we present a search for evidence of correlated variability between mm-wave, optical R band, and gamma-ray observations. This pilot study is paving the way for AGN monitoring with current and upcoming CMB experiments such as SPT-3G, Simons Observatory, and CMB-S4, including multi-wavelength studies with facilities such as VRO-LSST. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 9 pages, 3 figures, accepted to Astrophysical Journal Letters

Soliton Merger Rates and Enhanced Axion Dark Matter Decay

Authors: Xiaolong Du, David J. E. Marsh, Miguel Escudero, Andrew Benson, Diego Blas, Charis Kaur Pooni, Malcolm Fairbairn

Abstract: Solitons are observed to form in simulations of dark matter (DM) halos consisting of bosonic fields. We use the extended Press-Schechter formalism to compute the mass function of solitons, assuming various forms for the relationship between halo mass and soliton mass. We further provide a new calculation of the rate of soliton major mergers. Solitons composed of axion DM are unstable above a criti… ▽ More Solitons are observed to form in simulations of dark matter (DM) halos consisting of bosonic fields. We use the extended Press-Schechter formalism to compute the mass function of solitons, assuming various forms for the relationship between halo mass and soliton mass. We further provide a new calculation of the rate of soliton major mergers. Solitons composed of axion DM are unstable above a critical mass, and decay to either relativistic axions or photons, depending on the values of the coupling constants. We use the computed soliton major merger rate to predict the enhanced DM decay rate due to soliton instability. For certain values of currently allowed axion parameters, the energy injection into the intergalactic medium from soliton decays to photons is comparable to or larger than the energy injection due to core collapse supernovae at $z>10$. A companion paper explores the phenomenology of such an energy injection. △ Less

Submitted 12 February, 2024; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: 18 pages, 18 figures. Version 2: matches published version. Add an appendix discussing the timescale of axion star merging process due to dynamical friction

Report number: KCL-PH-TH-2023-03, CERN-TH-2023-009

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

Pushing the Limits of Detectability: Mixed Dark Matter from Strong Gravitational Lenses

Authors: Ryan E. Keeley, Anna M. Nierenberg, Daniel Gilman, Simon Birrer, Andrew Benson, Tommaso Treu

Abstract: One of the frontiers for advancing what is known about dark matter lies in using strong gravitational lenses to characterize the population of the smallest dark matter halos. There is a large volume of information in strong gravitational lens images -- the question we seek to answer is to what extent we can refine this information. To this end, we forecast the detectability of a mixed warm and col… ▽ More One of the frontiers for advancing what is known about dark matter lies in using strong gravitational lenses to characterize the population of the smallest dark matter halos. There is a large volume of information in strong gravitational lens images -- the question we seek to answer is to what extent we can refine this information. To this end, we forecast the detectability of a mixed warm and cold dark matter scenario using the anomalous flux ratio method from strong gravitational lensed images. The halo mass function of the mixed dark matter scenario is suppressed relative to cold dark matter but still predicts numerous low-mass dark matter halos relative to warm dark matter. Since the strong lens signal is a convolution over a range of dark matter halo masses and since the signal is sensitive to the specific configuration of dark matter halos, not just the halo mass function, degeneracies between different forms of suppression in the halo mass function, relative to cold dark matter, can arise. We find that, with a set of lenses with different configurations of the main deflector and hence different sensitivities to different mass ranges of the halo mass function, the different forms of suppression of the halo mass function between the warm dark matter model and the mixed dark matter model can be distinguished with $40$ lenses with Bayesian odds of 29.4:1. △ Less

Submitted 30 August, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 8 pages, 7 figures

Evaluating the sustainability of a de facto harvest strategy for British Columbia's Spot Prawn (Pandalus platyceros) fishery in the presence of environmental drivers of recruitment and hyperstable catch rates

Authors: Steven P. Rossi, Sean P. Cox, Samuel D. N. Johnson, Ashleen J. Benson

Abstract: The Spot Prawn trap fishery off the west coast of British Columbia (BC) is managed using a fixed escapement strategy that aims to prevent recruitment overfishing while maximizing expected long-term yield by closing the fishery when the catch rate of spawners, projected to the following spring, drops below 1.7 spawners per trap (the de jure rule). We develop a management strategy evaluation framewo… ▽ More The Spot Prawn trap fishery off the west coast of British Columbia (BC) is managed using a fixed escapement strategy that aims to prevent recruitment overfishing while maximizing expected long-term yield by closing the fishery when the catch rate of spawners, projected to the following spring, drops below 1.7 spawners per trap (the de jure rule). We develop a management strategy evaluation framework for BC's Spot Prawn fishery that examines the expected performance of the management procedure implemented in practice (the de facto rule), which was significantly more conservative than the de jure rule, usually closing the fishery when spawner catch rates were at least twice as high as specified by the de jure rule. Simulations indicate that the de facto spawner index rule using average empirical March 31st targets from 2000 to 2019 maintains most stocks near or above 0.8 BMSY with or without accounting for environmental effects and/or increasing future SST on recruitment. Abundance indices were found to be strongly hyperstable, with fishing efficiency 1.5 to 3.0 times higher under low biomass than high biomass. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 86 pages, 12 figures, 33 appendix figures

Evidence for AGN-Regulated Cooling in Clusters at $z \sim 1.4$: A Multi-Wavelength View of SPT-CL J0607-4448

Authors: Megan Masterson, Michael McDonald, Behzad Ansarinejad, Matthew Bayliss, Bradford A. Benson, Lindsey E. Bleem, Michael S. Calzadilla, Alastair C. Edge, Benjamin Floyd, Keunho J. Kim, Gourav Khullar, Taweewat Somboonpanyakul

Abstract: We present a multi-wavelength analysis of the galaxy cluster SPT-CL J0607-4448 (SPT0607), which is one of the most distant clusters discovered by the South Pole Telescope (SPT) at $z=1.4010\pm0.0028$. The high-redshift cluster shows clear signs of being relaxed with well-regulated feedback from the active galactic nucleus (AGN) in the brightest cluster galaxy (BCG). Using Chandra X-ray data, we co… ▽ More We present a multi-wavelength analysis of the galaxy cluster SPT-CL J0607-4448 (SPT0607), which is one of the most distant clusters discovered by the South Pole Telescope (SPT) at $z=1.4010\pm0.0028$. The high-redshift cluster shows clear signs of being relaxed with well-regulated feedback from the active galactic nucleus (AGN) in the brightest cluster galaxy (BCG). Using Chandra X-ray data, we construct thermodynamic profiles and determine the properties of the intracluster medium. The cool core nature of the cluster is supported by a centrally-peaked density profile and low central entropy ($K_0=18_{-9}^{+11}$ keV cm$^2$), which we estimate assuming an isothermal temperature profile due to the limited spectral information given the distance to the cluster. Using the density profile and gas cooling time inferred from the X-ray data, we find a mass cooling rate of $\dot{M}_\mathrm{cool}=100_{-60}^{+90}~M_\odot$ yr$^{-1}$. From optical spectroscopy and photometry around the [O II] emission line, we estimate that the BCG star formation rate is SFR$_\mathrm{[O~II]}=1.7_{-0.6}^{+1.0}~M_\odot$ yr$^{-1}$, roughly two orders of magnitude lower than the predicted mass cooling rate. In addition, using ATCA radio data at 2.1 GHz, we measure a radio jet power of $P_\mathrm{cav}=3.2_{-1.3}^{+2.1}\times10^{44}$ erg s$^{-1}$, which is consistent with the X-ray cooling luminosity ($L_\mathrm{cool}=1.9_{-0.5}^{+0.2}\times10^{44}$ erg s$^{-1}$ within $r_\mathrm{cool}=43$ kpc). These findings suggest that SPT0607 is a relaxed, cool core cluster with AGN-regulated cooling at an epoch shortly after cluster formation, implying that the balance between cooling and feedback can be reached quickly. We discuss implications for these findings on the evolution of AGN feedback in galaxy clusters. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: 13 pages, 5 figures, accepted for publication in ApJ

Growing the First Galaxies' Merger Trees

Authors: Ethan O. Nadler, Andrew Benson, Trey Driskell, Xiaolong Du, Vera Gluscevic

Abstract: Modelling the growth histories of specific galaxies often involves generating the entire population of objects that arise in a given cosmology and selecting systems with appropriate properties. This approach is highly inefficient when targeting rare systems such as the extremely luminous high-redshift galaxy candidates detected by JWST. Here, we present a novel framework for generating merger tree… ▽ More Modelling the growth histories of specific galaxies often involves generating the entire population of objects that arise in a given cosmology and selecting systems with appropriate properties. This approach is highly inefficient when targeting rare systems such as the extremely luminous high-redshift galaxy candidates detected by JWST. Here, we present a novel framework for generating merger trees with branches that are guaranteed to achieve a desired halo mass at a chosen redshift. This method augments extended Press Schechter theory solutions with constrained random processes known as Brownian bridges and is implemented in the open-source semi-analytic model $\texttt{Galacticus}$. We generate ensembles of constrained merger trees to predict the growth histories of seven high-redshift JWST galaxy candidates, finding that these systems most likely merge $\approx 2~\mathrm{Gyr}$ after the observation epoch and occupy haloes of mass $\gtrsim 10^{14}~M_{\mathrm{\odot}}$ today. These calculations are thousands of times more efficient than existing methods, are analytically controlled, and provide physical insights into the evolution of haloes with rapid early growth. Our constrained merger tree implementation is publicly available at http://github.com/galacticusorg/galacticus. △ Less

Submitted 23 March, 2023; v1 submitted 16 December, 2022; originally announced December 2022.

Comments: 20 pages, 12 figures. Updated to published version

Journal ref: MNRAS, 521, 3201 (2023)

A Measurement of the CMB Temperature Power Spectrum and Constraints on Cosmology from the SPT-3G 2018 TT/TE/EE Data Set

Authors: L. Balkenhol, D. Dutcher, A. Spurio Mancini, A. Doussot, K. Benabed, S. Galli, P. A. R. Ade, A. J. Anderson, B. Ansarinejad, M. Archipley, A. N. Bender, B. A. Benson, F. Bianchini, L. E. Bleem, F. R. Bouchet, L. Bryant, E. Camphuis, J. E. Carlstrom, T. W. Cecil, C. L. Chang, P. Chaubal, P. M. Chichura, T. -L. Chou, A. Coerver, T. M. Crawford , et al. (62 additional authors not shown)

Abstract: We present a sample-variance-limited measurement of the temperature power spectrum ($TT$) of the cosmic microwave background (CMB) using observations of a $\sim\! 1500 \,\mathrm{deg}^2$ field made by SPT-3G in 2018. We report multifrequency power spectrum measurements at 95, 150, and 220GHz covering the angular multipole range $750 \leq \ell < 3000$. We combine this $TT$ measurement with the publi… ▽ More We present a sample-variance-limited measurement of the temperature power spectrum ($TT$) of the cosmic microwave background (CMB) using observations of a $\sim\! 1500 \,\mathrm{deg}^2$ field made by SPT-3G in 2018. We report multifrequency power spectrum measurements at 95, 150, and 220GHz covering the angular multipole range $750 \leq \ell < 3000$. We combine this $TT$ measurement with the published polarization power spectrum measurements from the 2018 observing season and update their associated covariance matrix to complete the SPT-3G 2018 $TT/TE/EE$ data set. This is the first analysis to present cosmological constraints from SPT $TT$, $TE$, and $EE$ power spectrum measurements jointly. We blind the cosmological results and subject the data set to a series of consistency tests at the power spectrum and parameter level. We find excellent agreement between frequencies and spectrum types and our results are robust to the modeling of astrophysical foregrounds. We report results for $Λ$CDM and a series of extensions, drawing on the following parameters: the amplitude of the gravitational lensing effect on primary power spectra $A_\mathrm{L}$, the effective number of neutrino species $N_{\mathrm{eff}}$, the primordial helium abundance $Y_{\mathrm{P}}$, and the baryon clumping factor due to primordial magnetic fields $b$. We find that the SPT-3G 2018 $T/TE/EE$ data are well fit by $Λ$CDM with a probability-to-exceed of $15\%$. For $Λ$CDM, we constrain the expansion rate today to $H_0 = 68.3 \pm 1.5\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ and the combined structure growth parameter to $S_8 = 0.797 \pm 0.042$. The SPT-based results are effectively independent of Planck, and the cosmological parameter constraints from either data set are within $<1\,σ$ of each other. (abridged) △ Less

Submitted 27 July, 2023; v1 submitted 11 December, 2022; originally announced December 2022.

Comments: 35 Pages, 17 Figures, 11 Tables

Mapping gas around massive galaxies: cross-correlation of DES Y3 galaxies and Compton-$y$-maps from SPT and Planck

Authors: J. Sánchez, Y. Omori, C. Chang, L. E. Bleem, T. Crawford, A. Drlica-Wagner, S. Raghunathan, G. Zacharegkas, T. M. C. Abbott, M. Aguena, A. Alarcon, S. Allam, O. Alves, A. Amon, S. Avila, E. Baxter, K. Bechtol, B. A. Benson, G. M. Bernstein, E. Bertin, S. Bocquet, D. Brooks, D. L. Burke, A. Campos, J. E. Carlstrom , et al. (102 additional authors not shown)

Abstract: We cross-correlate positions of galaxies measured in data from the first three years of the Dark Energy Survey with Compton-$y$-maps generated using data from the South Pole Telescope (SPT) and the {\it Planck} mission. We model this cross-correlation measurement together with the galaxy auto-correlation to constrain the distribution of gas in the Universe. We measure the hydrostatic mass bias or,… ▽ More We cross-correlate positions of galaxies measured in data from the first three years of the Dark Energy Survey with Compton-$y$-maps generated using data from the South Pole Telescope (SPT) and the {\it Planck} mission. We model this cross-correlation measurement together with the galaxy auto-correlation to constrain the distribution of gas in the Universe. We measure the hydrostatic mass bias or, equivalently, the mean halo bias-weighted electron pressure $\langle b_{h}P_{e}\rangle$, using large-scale information. We find $\langle b_{h}P_{e}\rangle$ to be $[0.16^{+0.03}_{-0.04},0.28^{+0.04}_{-0.05},0.45^{+0.06}_{-0.10},0.54^{+0.08}_{-0.07},0.61^{+0.08}_{-0.06},0.63^{+0.07}_{-0.08}]$ meV cm$^{-3}$ at redshifts $z \sim [0.30, 0.46, 0.62,0.77, 0.89, 0.97]$. These values are consistent with previous work where measurements exist in the redshift range. We also constrain the mean gas profile using small-scale information, enabled by the high-resolution of the SPT data. We compare our measurements to different parametrized profiles based on the cosmo-OWLS hydrodynamical simulations. We find that our data are consistent with the simulation that assumes an AGN heating temperature of $10^{8.5}$K but are incompatible with the model that assumes an AGN heating temperature of $10^{8.0}$K. These comparisons indicate that the data prefer a higher value of electron pressure than the simulations within $r_{500c}$ of the galaxies' halos. △ Less

Submitted 18 October, 2022; v1 submitted 16 October, 2022; originally announced October 2022.

Comments: 20 pages, 9 figures. Submitted to MNRAS

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