Showing 1–50 of 94 results for author: Rapetti, D

TOI 762 A b and TIC 46432937 b: Two Giant Planets Transiting M Dwarf Stars

Authors: Joel D. Hartman, Daniel Bayliss, Rafael Brahm, Edward M. Bryant, Andrés Jordán, Gáspár Á. Bakos, Melissa J. Hobson, Elyar Sedaghati, Xavier Bonfils, Marion Cointepas, Jose Manuel Almenara, Khalid Barkaoui, Mathilde Timmermans, George Dransfield, Elsa Ducrot, Sebastián Zúñiga-Fernández, Matthew J. Hooton, Peter Pihlmann Pedersen, Francisco J. Pozuelos, Amaury H. M. J. Triaud, Michaël Gillon, Emmanuel Jehin, William C. Waalkes, Zachory K. Berta-Thompson, Steve B. Howell , et al. (11 additional authors not shown)

Abstract: We present the discovery of TOI 762 A b and TIC 46432937 b, two giant planets transiting M dwarf stars. Transits of both systems were first detected from observations by the NASA TESS mission, and the transiting objects are confirmed as planets through high-precision radial velocity (RV) observations carried out with VLT/ESPRESSO. TOI 762 A b is a warm sub-Saturn with a mass of 0.251 +- 0.042 M_J,… ▽ More We present the discovery of TOI 762 A b and TIC 46432937 b, two giant planets transiting M dwarf stars. Transits of both systems were first detected from observations by the NASA TESS mission, and the transiting objects are confirmed as planets through high-precision radial velocity (RV) observations carried out with VLT/ESPRESSO. TOI 762 A b is a warm sub-Saturn with a mass of 0.251 +- 0.042 M_J, a radius of 0.744 +- 0.017 R_J, and an orbital period of 3.4717 d. It transits a mid-M dwarf star with a mass of 0.442 +- 0.025 M_S and a radius of 0.4250 +- 0.0091 R_S. The star TOI 762 A has a resolved binary star companion TOI 762 B that is separated from TOI 762 A by 3.2" (~ 319 AU) and has an estimated mass of 0.227 +- 0.010 M_S. The planet TIC 46432937 b is a warm Super-Jupiter with a mass of 3.20 +- 0.11 M_J and radius of 1.188 +- 0.030 R_J. The planet's orbital period is P = 1.4404 d, and it undergoes grazing transits of its early M dwarf host star, which has a mass of 0.563 +- 0.029 M_S and a radius of 0.5299 +- 0.0091 R_S. TIC 46432937 b is one of the highest mass planets found to date transiting an M dwarf star. TIC 46432937 b is also a promising target for atmospheric observations, having the highest Transmission Spectroscopy Metric or Emission Spectroscopy Metric value of any known warm Super-Jupiter (mass greater than 3.0 M_J, equilibrium temperature below 1000 K). △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 29 pages, 7 figures, 8 tables, accepted for publication in AAS Journals

Gliese 12 b: A temperate Earth-sized planet at 12 pc ideal for atmospheric transmission spectroscopy

Authors: M. Kuzuhara, A. Fukui, J. H. Livingston, J. A. Caballero, J. P. de Leon, T. Hirano, Y. Kasagi, F. Murgas, N. Narita, M. Omiya, Jaume Orell-Miquel, E. Palle, Q. Changeat, E. Esparza-Borges, H. Harakawa, C. Hellier, Yasunori Hori, Kai Ikuta, H. T. Ishikawa, T. Kodama, T. Kotani, T. Kudo, J. C. Morales, M. Mori, E. Nagel , et al. (81 additional authors not shown)

Abstract: Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We repor… ▽ More Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We report the discovery of an Earth-sized planet transiting the nearby (12 pc) inactive M3.0 dwarf Gliese 12 (TOI-6251) with an orbital period ($P_{\rm{orb}}$) of 12.76 days. The planet, Gliese 12b, was initially identified as a candidate with an ambiguous $P_{\rm{orb}}$ from TESS data. We confirmed the transit signal and $P_{\rm{orb}}$ using ground-based photometry with MuSCAT2 and MuSCAT3, and validated the planetary nature of the signal using high-resolution images from Gemini/NIRI and Keck/NIRC2 as well as radial velocity (RV) measurements from the InfraRed Doppler instrument on the Subaru 8.2 m telescope and from CARMENES on the CAHA 3.5 m telescope. X-ray observations with XMM-Newton showed the host star is inactive, with an X-ray-to-bolometric luminosity ratio of $\log L_{\rm X}/L_{\rm bol} \approx -5.7$. Joint analysis of the light curves and RV measurements revealed that Gliese 12b has a radius of 0.96 $\pm$ 0.05 $R_\oplus$, a 3$σ$ mass upper limit of 3.9 $M_\oplus$, and an equilibrium temperature of 315 $\pm$ 6 K assuming zero albedo. The transmission spectroscopy metric (TSM) value of Gliese 12b is close to the TSM values of the TRAPPIST-1 planets, adding Gliese 12b to the small list of potentially terrestrial, temperate planets amenable to atmospheric characterization with JWST. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 29 pages (20 pages in main body), 13 figures (10 figures in main body). Equal contributions from M. K. and A. F.. Accepted for Publication in ApJL at 2024 March 21

Journal ref: Published on 2024 May 23 by Astrophysical Journal Letters (ApJL) 967 L21

Discovery of two warm mini-Neptunes with contrasting densities orbiting the young K3V star TOI-815

Authors: Angelica Psaridi, Hugh Osborn, François Bouchy, Monika Lendl, Léna Parc, Nicolas Billot, Christopher Broeg, Sérgio G. Sousa, Vardan Adibekyan, Omar Attia, Andrea Bonfanti, Hritam Chakraborty, Karen A. Collins, Jeanne Davoult, Elisa Delgado-Mena, Nolan Grieves, Tristan Guillot, Alexis Heitzmann, Ravit Helled, Coel Hellier, Jon M. Jenkins, Henrik Knierim, Andreas Krenn, JackJ. Lissauer, Rafael Luque , et al. (108 additional authors not shown)

Abstract: We present the discovery and characterization of two warm mini-Neptunes transiting the K3V star TOI-815 in a K-M binary system. Analysis of the spectra and rotation period reveal it to be a young star with an age of $200^{+400}_{-200}$Myr. TOI-815b has a 11.2-day period and a radius of 2.94$\pm$0.05$\it{R_{\rm\mathrm{\oplus}}}$ with transits observed by TESS, CHEOPS, ASTEP, and LCOGT. The outer pl… ▽ More We present the discovery and characterization of two warm mini-Neptunes transiting the K3V star TOI-815 in a K-M binary system. Analysis of the spectra and rotation period reveal it to be a young star with an age of $200^{+400}_{-200}$Myr. TOI-815b has a 11.2-day period and a radius of 2.94$\pm$0.05$\it{R_{\rm\mathrm{\oplus}}}$ with transits observed by TESS, CHEOPS, ASTEP, and LCOGT. The outer planet, TOI-815c, has a radius of 2.62$\pm$0.10$\it{R_{\rm\mathrm{\oplus}}}$, based on observations of three non-consecutive transits with TESS, while targeted CHEOPS photometry and radial velocity follow-up with ESPRESSO were required to confirm the 35-day period. ESPRESSO confirmed the planetary nature of both planets and measured masses of 7.6$\pm$1.5 $\it{M_{\rm \mathrm{\oplus}}}$ ($ρ_\mathrm{P}$=1.64$^{+0.33}_{-0.31}$gcm$^{-3}$) and 23.5$\pm$2.4$\it{M_{\rm\mathrm{\oplus}}}$ ($ρ_\mathrm{P}$=7.2$^{+1.1}_{-1.0}$gcm$^{-3}$) respectively. Thus, the planets have very different masses, unlike the usual similarity of masses in compact multi-planet systems. Moreover, our statistical analysis of mini-Neptunes orbiting FGK stars suggests that weakly irradiated planets tend to have higher bulk densities compared to those suffering strong irradiation. This could be ascribed to their cooler atmospheres, which are more compressed and denser. Internal structure modeling of TOI-815b suggests it likely has a H-He atmosphere constituting a few percent of the total planet mass, or higher if the planet is assumed to have no water. In contrast, the measured mass and radius of TOI-815c can be explained without invoking any atmosphere, challenging planetary formation theories. Finally, we infer from our measurements that the star is viewed close to pole-on, which implies a spin-orbit misalignment at the 3$σ$ level. △ Less

Submitted 30 January, 2024; v1 submitted 28 January, 2024; originally announced January 2024.

Comments: 24 pages, 27 figures, 6 tables

A resonant sextuplet of sub-Neptunes transiting the bright star HD 110067

Authors: R. Luque, H. P. Osborn, A. Leleu, E. Pallé, A. Bonfanti, O. Barragán, T. G. Wilson, C. Broeg, A. Collier Cameron, M. Lendl, P. F. L. Maxted, Y. Alibert, D. Gandolfi, J. -B. Delisle, M. J. Hooton, J. A. Egger, G. Nowak, M. Lafarga, D. Rapetti, J. D. Twicken, J. C. Morales, I. Carleo, J. Orell-Miquel, V. Adibekyan, R. Alonso , et al. (127 additional authors not shown)

Abstract: Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial con… ▽ More Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial conditions and environment. Those in resonance (with their orbital periods related by a ratio of small integers) are particularly valuable because they imply a system architecture practically unchanged since its birth. Here, we present the observations of six transiting planets around the bright nearby star HD 110067. We find that the planets follow a chain of resonant orbits. A dynamical study of the innermost planet triplet allowed the prediction and later confirmation of the orbits of the rest of the planets in the system. The six planets are found to be sub-Neptunes with radii ranging from 1.94 to 2.85 Re. Three of the planets have measured masses, yielding low bulk densities that suggest the presence of large hydrogen-dominated atmospheres. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Published in Nature on November 30, 2023. Supplementary Information can be found in the online version of the paper in the journal

Journal ref: Nature 623, 932-937 (2023)

The SOPHIE search for northern extrasolar planets-XIX. A system including a cold sub-Neptune potentially transiting a V = 6.5 star HD88986

Authors: N. Heidari, I. Boisse, N. C. Hara, T. G. Wilson, F. Kiefer, G. Hébrard, F. Philipot, S. Hoyer, K. G. Stassun, G. W. Henry, N. C. Santos, L. Acuña, D. Almasian, L. Arnold, N. Astudillo-Defru, O. Attia, X. Bonfils, F. Bouchy, V. Bourrier, B. Collet, P. Cortés-Zuleta, A. Carmona, X. Delfosse, S. Dalal, M. Deleuil , et al. (29 additional authors not shown)

Abstract: Transiting planets with orbital periods longer than 40 d are extremely rare among the 5000+ planets discovered so far. The lack of discoveries of this population poses a challenge to research into planetary demographics, formation, and evolution. Here, we present the detection and characterization of HD88986b, a potentially transiting sub-Neptune, possessing the longest orbital period among known… ▽ More Transiting planets with orbital periods longer than 40 d are extremely rare among the 5000+ planets discovered so far. The lack of discoveries of this population poses a challenge to research into planetary demographics, formation, and evolution. Here, we present the detection and characterization of HD88986b, a potentially transiting sub-Neptune, possessing the longest orbital period among known transiting small planets (< 4 R$_{\oplus}$) with a precise mass measurement ($σ_M/M$ > 25%). Additionally, we identified the presence of a massive companion in a wider orbit around HD88986. Our analysis reveals that HD88986b, based on two potential single transits on sector 21 and sector 48 which are both consistent with the predicted transit time from the RV model, is potentially transiting. The joint analysis of RV and photometric data show that HD88986b has a radius of 2.49$\pm$0.18 R$_{\oplus}$, a mass of 17.2$^{+4.0}_{-3.8}$ M$_{\oplus}$, and it orbits every 146.05$^{+0.43}_{-0.40}$ d around a subgiant HD88986 which is one of the closest and brightest exoplanet host stars (G2V type, R=1.543 $\pm$0.065 R$_{\odot}$, V=$6.47\pm 0.01$ mag, distance=33.37$\pm$0.04 pc). The nature of the outer, massive companion is still to be confirmed; a joint analysis of RVs, Hipparcos, and Gaia astrometric data shows that with a 3$σ$ confidence interval, its semi-major axis is between 16.7 and 38.8 au and its mass is between 68 and 284 M$_{Jup}$. HD88986b's wide orbit suggests the planet did not undergo significant mass loss due to extreme-ultraviolet radiation from its host star. Therefore, it probably maintained its primordial composition, allowing us to probe its formation scenario. Furthermore, the cold nature of HD88986b (460$\pm$8 K), thanks to its long orbital period, will open up exciting opportunities for future studies of cold atmosphere composition characterization. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 37 pages, accepted to be published in A&A

Journal ref: A&A 681, A55 (2024)

Validating posteriors obtained by an emulator when jointly-fitting mock data of the global 21-cm signal and high-z galaxy UV luminosity function

Authors: J. Dorigo Jones, D. Rapetti, J. Mirocha, J. J. Hibbard, J. O. Burns, N. Bassett

Abstract: Although neural-network-based emulators enable efficient parameter estimation in 21-cm cosmology, the accuracy of such constraints is poorly understood. We employ nested sampling to fit mock data of the global 21-cm signal and high-$z$ galaxy ultraviolet luminosity function (UVLF) and compare for the first time the emulated posteriors obtained using the global signal emulator ${\tt globalemu}$ to… ▽ More Although neural-network-based emulators enable efficient parameter estimation in 21-cm cosmology, the accuracy of such constraints is poorly understood. We employ nested sampling to fit mock data of the global 21-cm signal and high-$z$ galaxy ultraviolet luminosity function (UVLF) and compare for the first time the emulated posteriors obtained using the global signal emulator ${\tt globalemu}$ to the `true' posteriors obtained using the full model on which the emulator is trained using ${\tt ARES}$. Of the eight model parameters we employ, four control the star formation efficiency (SFE), and thus can be constrained by UVLF data, while the remaining four control UV and X-ray photon production, and the minimum virial temperature of star-forming halos ($T_{\rm min}$), and thus are uniquely probed by reionization and 21-cm measurements. For noise levels of 50 and 250 mK in the 21-cm data being jointly-fit, the emulated and `true' posteriors are consistent to within $1σ$. However, at lower noise levels of 10 and 25 mK, ${\tt globalemu}$ overpredicts $T_{\rm min}$ and underpredicts $γ_{\rm lo}$, an SFE parameter, by $\approx3-4σ$, while the `true' ${\tt ARES}$ posteriors capture their fiducial values within $1σ$. We find that jointly-fitting the mock UVLF and 21-cm data significantly improves constraints on the SFE parameters by breaking degeneracies in the ${\tt ARES}$ parameter space. Our results demonstrate the astrophysical constraints that can be expected for global 21-cm experiments for a range of noise levels from pessimistic to optimistic, and also the potential for probing redshift evolution of SFE parameters by including UVLF data. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 19 pages, 8 figures, 2 tables. Accepted by ApJ

TOI-4860 b, a short-period giant planet transiting an M3.5 dwarf

Authors: J. M. Almenara, X. Bonfils, E. M. Bryant, A. Jordán, G. Hébrard, E. Martioli, A. C. M. Correia, N. Astudillo-Defru, C. Cadieux, L. Arnold, É. Artigau, G. Á. Bakos, S. C. C. Barros, D. Bayliss, F. Bouchy, G. Boué, R. Brahm, A. Carmona, D. Charbonneau, D. R. Ciardi, R. Cloutier, M. Cointepas, N. J. Cook, N. B. Cowan, X. Delfosse , et al. (25 additional authors not shown)

Abstract: We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03… ▽ More We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03 R$_J$ and an orbital period of 1.52 days. With high-resolution spectroscopy taken by the CFHT/SPIRou and ESO/ESPRESSO spectrographs, we refined the host star parameters ([Fe/H] = 0.27 $\pm$ 0.12) and measured the mass of the planet (0.273 $\pm$ 0.006 M$_J$). Based on these measurements, TOI-4860 b joins the small set of massive planets ($>$80 M$_E$) found around mid to late M dwarfs ($<$0.4 R$_\odot$), providing both an interesting challenge to planet formation theory and a favourable target for further atmospheric studies with transmission spectroscopy. We identified an additional signal in the radial velocity data that we attribute to an eccentric planet candidate ($e=0.66\pm0.09$) with an orbital period of $427\pm7$~days and a minimum mass of $1.66\pm 0.26$ M$_J$, but additional data would be needed to confirm this. △ Less

Submitted 12 January, 2024; v1 submitted 2 August, 2023; originally announced August 2023.

Comments: 16 pages, 14 figures, accepted for publication in A&A

TOI-1859b: A 64-Day Warm Jupiter on an Eccentric and Misaligned Orbit

Authors: Jiayin Dong, Songhu Wang, Malena Rice, George Zhou, Chelsea X. Huang, Rebekah I. Dawson, Gudmundur K. Stefánsson, Samuel Halverson, Shubham Kanodia, Suvrath Mahadevan, Michael W. McElwain, Jaime A. Alvarado-Montes, Joe P. Ninan, Paul Robertson, Arpita Roy, Christian Schwab, Sarah E. Logsdon, Ryan C. Terrien, Karen A. Collins, Gregor Srdoc, Ramotholo Sefako, Didier Laloum, David W. Latham, Allyson Bieryla, Paul A. Dalba , et al. (9 additional authors not shown)

Abstract: Warm Jupiters are close-in giant planets with relatively large planet-star separations (i.e., $10< a/R_\star <100$). Given their weak tidal interactions with their host stars, measurements of stellar obliquity may be used to probe the initial obliquity distribution and dynamical history for close-in gas giants. Using spectroscopic observations, we confirm the planetary nature of TOI-1859b and dete… ▽ More Warm Jupiters are close-in giant planets with relatively large planet-star separations (i.e., $10< a/R_\star <100$). Given their weak tidal interactions with their host stars, measurements of stellar obliquity may be used to probe the initial obliquity distribution and dynamical history for close-in gas giants. Using spectroscopic observations, we confirm the planetary nature of TOI-1859b and determine the stellar obliquity of TOI-1859 to be $λ= 38.9^{+2.8}_{-2.7}°$ relative to its planetary companion using the Rossiter-McLaughlin effect. TOI-1859b is a 64-day warm Jupiter orbiting around a late-F dwarf and has an orbital eccentricity of $0.57^{+0.12}_{-0.16}$, inferred purely from transit light curves. The eccentric and misaligned orbit of TOI-1859b is likely an outcome of dynamical interactions, such as planet-planet scattering and planet-disk resonance crossing. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 15 pages, 7 figures, 1 table; accepted to ApJL

Fitting and Comparing Galactic Foreground Models for Unbiased 21-cm Cosmology

Authors: Joshua J. Hibbard, David Rapetti, Jack O. Burns, Nivedita Mahesh, Neil Bassett

Abstract: Accurate detection of the cosmological 21-cm global signal requires galactic foreground models which can remove power over ~$10^6$. Although foreground and global signal models unavoidably exhibit overlap in their vector-spaces inducing bias error in the extracted signal, a second source of bias and error arises from inadequate foreground models, i.e. models which cannot fit spectra down to the no… ▽ More Accurate detection of the cosmological 21-cm global signal requires galactic foreground models which can remove power over ~$10^6$. Although foreground and global signal models unavoidably exhibit overlap in their vector-spaces inducing bias error in the extracted signal, a second source of bias and error arises from inadequate foreground models, i.e. models which cannot fit spectra down to the noise level of the signal. We therefore test the level to which seven commonly employed foreground models -- including nonlinear and linear forward-models, polynomials, and maximally-smooth polynomials -- fit realistic simulated mock foreground spectra, as well as their dependence upon model inputs. The mock spectra are synthesized for an EDGES-like experiment and we compare all models' goodness-of-fit and preference using a Kolomogorov-Smirnov test of the noise-normalized residuals in order to compare models with differing, and sometimes indeterminable, degrees of freedom. For a single LST bin spectrum and p-value threshold of $p=0.05$, the nonlinear-forward model with 4 parameters is preferred ($p=0.99$), while the linear forward-model fits well with 6-7 parameters ($p=0.94,0.97$ respectively). The polynomials and maximally-smooth polynomials, like those employed by the EDGES and SARAS3 experiments, cannot produce good fits with 5 parameters for the experimental simulations in this work ($p<10^{-6}$). However, we find that polynomials with 6 parameters pass the KS-test ($p=0.4$), although a 9 parameter fit produces the highest p-value ($p\sim0.67$). When fitting multiple LST bins simultaneously, we find that the linear forward-model outperforms (a higher p-value) the nonlinear for 2, 5 and 10 LST bins. Importantly, the KS-test consistently identifies best-fit \textit{and} preferred models. △ Less

Submitted 30 October, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: 26 pages, 12 figures, accepted for publication by ApJ, minor revisions and edits from previous version

TimeSOAP: Tracking high-dimensional fluctuations in complex molecular systems via time-variations of SOAP spectra

Authors: Cristina Caruso, Annalisa Cardellini, Martina Crippa, Daniele Rapetti, Giovanni M. Pavan

Abstract: Many molecular systems and physical phenomena are controlled by local fluctuations and microscopic dynamical rearrangements of the constitutive interacting units that are often difficult to detect. This is the case, for example, of phase transitions, phase equilibria, nucleation events, and defect propagation, to mention a few. A detailed comprehension of local atomic environments and of their dyn… ▽ More Many molecular systems and physical phenomena are controlled by local fluctuations and microscopic dynamical rearrangements of the constitutive interacting units that are often difficult to detect. This is the case, for example, of phase transitions, phase equilibria, nucleation events, and defect propagation, to mention a few. A detailed comprehension of local atomic environments and of their dynamic rearrangements is essential to understand such phenomena, and also to draw structure-property relationships useful to unveil how to control complex molecular systems. Considerable progresses in the development of advanced high-dimensional structural descriptors (e.g., Smooth Overlap of Atomic Position (SOAP), etc.) have certainly enhanced the representation of atomic-scale simulations data. However, despite such efforts, local dynamic environment rearrangements remain still difficult to elucidate. Here, exploiting the structural-rich description of atomic environments of SOAP and building on the concept of time-dependent local variations, we developed a time-dependent SOAP-based descriptor, TimeSOAP (tSOAP), which essentially tracks the time variations in the local SOAP environments surrounding each molecule (i.e., each SOAP center) in complex molecular systems along ensemble trajectories. We demonstrate how analysis of the time-series {tSOAP data and of their time-derivatives allows detecting dynamics domains and tracking instantaneous changes of local atomic arrangements (i.e., local fluctuations) in a variety of molecular systems. The approach is simple and general, and we expect will help to shed light on a variety of complex dynamical phenomena. △ Less

Submitted 21 February, 2023; v1 submitted 19 February, 2023; originally announced February 2023.

Comments: 20 pages, 9 figures

LuSEE 'Night': The Lunar Surface Electromagnetics Experiment

Authors: Stuart D. Bale, Neil Bassett, Jack O. Burns, Johnny Dorigo Jones, Keith Goetz, Christian Hellum-Bye, Sven Hermann, Joshua Hibbard, Milan Maksimovic, Ryan McLean, Raul Monsalve, Paul O'Connor, Aaron Parsons, Marc Pulupa, Rugved Pund, David Rapetti, Kaja M. Rotermund, Ben Saliwanchik, Anze Slosar, David Sundkvist, Aritoki Suzuki

Abstract: The Lunar Surface Electromagnetics Explorer 'LuSEE Night' is a low frequency radio astronomy experiment that will be delivered to the farside of the Moon by the NASA Commercial Lunar Payload Services (CLPS) program in late 2025 or early 2026. The payload system is being developed jointly by NASA and the US Department of Energy (DOE) and consists of a 4 channel, 50 MHz Nyquist baseband receiver sys… ▽ More The Lunar Surface Electromagnetics Explorer 'LuSEE Night' is a low frequency radio astronomy experiment that will be delivered to the farside of the Moon by the NASA Commercial Lunar Payload Services (CLPS) program in late 2025 or early 2026. The payload system is being developed jointly by NASA and the US Department of Energy (DOE) and consists of a 4 channel, 50 MHz Nyquist baseband receiver system and 2 orthogonal $\sim$6m tip-to-tip electric dipole antennas. LuSEE Night will enjoy standalone operations through the lunar night, without the electromagnetic interference (EMI) of an operating lander system and antipodal to our noisy home planet. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: summary paper submitted to URSI GASS 2023

HIP 33609 b: An Eccentric Brown Dwarf Transiting a V=7.3 Rapidly Rotating B-Star

Authors: Noah Vowell, Joseph E. Rodriguez, Samuel N. Quinn, George Zhou, Andrew Vanderburg, Andrew W. Mann, Matthew J. Hooton, Keivan G. Stassun, Saburo Howard, Allyson Bieryla, David W. Latham, Steve B. Howell, Tristan Guillot, Carl Ziegler, Karen A. Collins, Theron W. Carmichael, Jon M. Jenkins, Avi Shporer, Lyu ABE, Philippe Bendjoya, Jonathan L. Bush, Marco Buttu, Kevin I. Collins, Jason D. Eastman, Matthew J. Fields , et al. (19 additional authors not shown)

Abstract: We present the discovery and characterization of HIP 33609 b, a transiting warm brown dwarf orbiting a late B star, discovered by NASA's Transiting Exoplanet Survey Satellite TESS as TOI-588 b. HIP 33609 b is a large (R$_{b}$ = 1.580$_{-0.070}^{+0.074}$ R$_{J}$) brown dwarf on a highly eccentric (e = 0.560$_{-0.031}^{+0.029}$) orbit with a 39-day period. The host star is a bright (V = 7.3 mag), T… ▽ More We present the discovery and characterization of HIP 33609 b, a transiting warm brown dwarf orbiting a late B star, discovered by NASA's Transiting Exoplanet Survey Satellite TESS as TOI-588 b. HIP 33609 b is a large (R$_{b}$ = 1.580$_{-0.070}^{+0.074}$ R$_{J}$) brown dwarf on a highly eccentric (e = 0.560$_{-0.031}^{+0.029}$) orbit with a 39-day period. The host star is a bright (V = 7.3 mag), T$_{eff}$ = 10,400$_{-660}^{+800}$ K star with a mass of M$_{*}$ = 2.383$_{-0.095}^{+0.10}$ M$_{\odot}$ and radius of R$_{*}$ = 1.863$_{-0.082}^{+0.087}$ R$_{\odot}$, making it the hottest transiting brown dwarf host star discovered to date. We obtained radial velocity measurements from the CHIRON spectrograph confirming the companion's mass of M$_{b}$ = 68.0$_{-7.1}^{+7.4}$ M$_{J}$ as well as the host star's rotation rate ($vsini_{*} = 55.6 \pm 1.8$ km/s). We also present the discovery of a new comoving group of stars, designated as MELANGE-6, and determine that HIP 33609 is a member. We use a combination of rotation periods and isochrone models fit to the cluster members to estimate an age of 150 $\pm$ 25 Myr. With a measured mass, radius, and age, HIP 33609 b becomes a benchmark for substellar evolutionary models. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 16 pages, 5 figures, 4 tables, Submitted to AAS Journals

Constraining a Model of the Radio Sky Below 6 MHz Using the Parker Solar Probe/FIELDS Instrument in Preparation for Upcoming Lunar-based Experiments

Authors: Neil Bassett, David Rapetti, Bang D. Nhan, Brent Page, Jack O. Burns, Marc Pulupa, Stuart D. Bale

Abstract: We present a Bayesian analysis of data from the FIELDS instrument on board the Parker Solar Probe (PSP) spacecraft with the aim of constraining low frequency ($\lesssim$ 6 MHz) sky in preparation for several upcoming lunar-based experiments. We utilize data recorded during PSP's ``coning roll'' maneuvers, in which the axis of the spacecraft is pointed 45$^{\circ}$ off of the Sun. The spacecraft th… ▽ More We present a Bayesian analysis of data from the FIELDS instrument on board the Parker Solar Probe (PSP) spacecraft with the aim of constraining low frequency ($\lesssim$ 6 MHz) sky in preparation for several upcoming lunar-based experiments. We utilize data recorded during PSP's ``coning roll'' maneuvers, in which the axis of the spacecraft is pointed 45$^{\circ}$ off of the Sun. The spacecraft then rotates about a line between the Sun and the spacecraft with a period of 24 minutes. We reduce the data into two formats: roll-averaged, in which the spectra are averaged over the roll, and phase-binned, in which the spectra are binned according to the phase of the roll. We construct a forward model of the FIELDS observations that includes numerical simulations of the antenna beam, an analytic emissivity function of the galaxy, and estimates of the absorption due to free electrons. Fitting 5 parameters, we find that the roll-averaged data can be fit well by this model and we obtain posterior parameter constraints that are in general agreement with previous estimates. The model is not, however, able to fit the phase-binned data well, likely due to limitations such as the lack of non-smooth emission structure at both small and large scales, enforced symmetry between the northern and southern galactic hemispheres, and large uncertainties in the free electron density. This suggests that significant improvement in the low frequency sky model is needed in order to fully and accurately represent the sky at frequencies below 6 MHz. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 18 pages, 10 figures, 5 tables. Under review in the Astrophysical Journal

Innate Dynamics and Identity Crisis of a Metal Surface Unveiled by Machine Learning of Atomic Environments

Authors: Matteo Cioni, Daniela Polino, Daniele Rapetti, Luca Pesce, Massimo Delle Piane, Giovanni M. Pavan

Abstract: Metals are traditionally considered hard matter. However, it is well known that their atomic lattices may become dynamic and undergo reconfigurations even well-below the melting temperature. The innate atomic dynamics of metals is directly related to their bulk and surface properties. Understanding their complex structural dynamics is thus important for many applications but is not easy. Here we r… ▽ More Metals are traditionally considered hard matter. However, it is well known that their atomic lattices may become dynamic and undergo reconfigurations even well-below the melting temperature. The innate atomic dynamics of metals is directly related to their bulk and surface properties. Understanding their complex structural dynamics is thus important for many applications but is not easy. Here we report deep-potential molecular dynamics simulations allowing to resolve at atomic-resolution the complex dynamics of various types of copper (Cu) surfaces, used as an example, near the Hüttig ($\sim1/3$ of melting) temperature. The development of a deep neural network potential trained on DFT calculations provides a dynamically-accurate force field that we use to simulate large atomistic models of different Cu surface types. A combination of high-dimensional structural descriptors and unsupervised machine learning allows identifying and tracking all the atomic environments (AEs) emerging in the surfaces at finite temperatures. We can directly observe how AEs that are non-native in a specific (ideal) surface, but that are instead typical of other surface types, continuously emerge/disappear in that surface in relevant regimes in dynamic equilibrium with the native ones. Our analyses allow estimating the lifetime of all the AEs populating these Cu surfaces and to reconstruct their dynamic interconversions networks. This reveals the elusive identity of these metal surfaces, which preserve their identity only in part and in part transform into something else in relevant conditions. This also proposes a concept of "statistical identity" for metal surfaces, which is key for understanding their behaviors and properties. △ Less

Submitted 21 February, 2023; v1 submitted 29 July, 2022; originally announced July 2022.

On the Origin of Diffuse Radio Emission in Abell 85 -- Insights from new GMRT Observations

Authors: Majidul Rahaman, Ramij Raja, Abhirup Datta, Jack O Burns, David Rapetti

Abstract: Extended, steep, and ultra-steep spectrum radio emission in a galaxy cluster is usually associated with recent mergers. Simulations show that radio phoenixes are aged radio galaxy lobes whose emission reactivates when a low Mach shock compresses it. A85 hosts a textbook example of a radio phoenix at about 320 kpc southwest of the cluster center. We present a new high resolution 325 MHz GMRT radio… ▽ More Extended, steep, and ultra-steep spectrum radio emission in a galaxy cluster is usually associated with recent mergers. Simulations show that radio phoenixes are aged radio galaxy lobes whose emission reactivates when a low Mach shock compresses it. A85 hosts a textbook example of a radio phoenix at about 320 kpc southwest of the cluster center. We present a new high resolution 325 MHz GMRT radio map illustrating this radio phoenix's complex and filamentary structure. The full extent of the radio structure is revealed for the first time from these radio images of A85. Using archival \textit{Chandra} X-ray observations, we applied an automated 2-D shock finder to the X-ray surface brightness and Adaptive Circular Binning (ACB) temperature maps which confirmed a bow shock at the location of the radio phoenix. We also compared the Mach number from the X-ray data with the radio-derived Mach number in the same region using multi-frequency radio observations and find that they are consistent within the 1$σ$ error level. △ Less

Submitted 30 June, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS (11 Pages, 6 Figures, 4 Tables)

X-ray morphology of cluster-mass haloes in self-interacting dark matter

Authors: Xuejian Shen, Thejs Brinckmann, David Rapetti, Mark Vogelsberger, Adam Mantz, Jesús Zavala, Steven W. Allen

Abstract: We perform cosmological zoom-in simulations of $19$ relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq 1\operatorname{-}3\times 10^{15}\,{\rm M}_{\odot}$. Both the dark matter and the intracluster ga… ▽ More We perform cosmological zoom-in simulations of $19$ relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq 1\operatorname{-}3\times 10^{15}\,{\rm M}_{\odot}$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of $84$ relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\operatorname{-}0.2\,R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(σ/m)= 0.5\operatorname{-}1~{\rm cm}^2/{\rm g}$ than CDM and SIDM with $(σ/m)=0.1\,{\rm cm}^2/{\rm g}$. Our statistical analysis indicates that the latter two models are disfavored at the $68\%$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations. △ Less

Submitted 1 November, 2022; v1 submitted 31 January, 2022; originally announced February 2022.

Comments: Accepted by MNRAS

Constraining Warm Dark Matter and Pop III stars with the Global 21-cm Signal

Authors: Joshua J. Hibbard, Jordan Mirocha, David Rapetti, Neil Bassett, Jack O. Burns, Keith Tauscher

Abstract: Upcoming ground and space-based experiments may have sufficient accuracy to place significant constraints upon high-redshift star formation, Reionization, and dark matter (DM) using the global 21-cm signal of the intergalactic medium. In the early universe, when the relative abundance of low-mass DM halos is important, measuring the global signal would place constraints on the damping of structure… ▽ More Upcoming ground and space-based experiments may have sufficient accuracy to place significant constraints upon high-redshift star formation, Reionization, and dark matter (DM) using the global 21-cm signal of the intergalactic medium. In the early universe, when the relative abundance of low-mass DM halos is important, measuring the global signal would place constraints on the damping of structure formation caused by DM having a higher relic velocity (warm dark matter, or WDM) than in cold dark matter (CDM). Such damping, however, can be mimicked by altering the star formation efficiency (SFE) and difficult to detect because of the presence of Pop III stars with unknown properties. We study these various cases and their degeneracies with the WDM mass parameter $m_X$ using a Fisher matrix analysis. We study the $m_X = 7$ keV case and a star-formation model that parametrizes the SFE as a strong function of halo mass and include several variations of this model along with three different input noise levels for the likelihood; we also use a minimum halo virial temperature for collapse near the molecular cooling threshold. We find that when the likelihood includes only Pop II stars, $m_X$ is constrained to an uncertainty of $\sim 0.4$ keV for all models and noise levels at 68$\%$ CI. When the likelihood includes weak Pop III stars, $m_X \sim 0.3$ keV, and if Pop III star formation is relatively efficient, $m_X \sim 0.1$ keV uncertainty, with tight Pop III star-formation parameter constraints. Our results show that the global 21-cm signal is a promising test-bed for WDM models, even in the presence of strong degeneracies with astrophysical parameters. △ Less

Submitted 31 March, 2022; v1 submitted 7 January, 2022; originally announced January 2022.

Comments: 22 pages, 9 figures, updated to version accepted by ApJ

Cosmological Constraints from Gas Mass Fractions of Massive, Relaxed Galaxy Clusters

Authors: Adam B. Mantz, Steven W. Allen, Rebecca E. A. Canning, Lucie Baumont, Bradford Benson, Lindsey E. Bleem, Steven R. Ehlert, Benjamin Floyd, Ricardo Herbonnet, Patrick L. Kelly, Shuang Liang, Anja von der Linden, Michael McDonald, David A. Rapetti, Robert W. Schmidt, Norbert Werner, Adam Wright

Abstract: We present updated cosmological constraints from measurements of the gas mass fractions ($f_{gas}$) of massive, dynamically relaxed galaxy clusters. Our new data set has greater leverage on models of dark energy, thanks to the addition of the Perseus Cluster at low redshifts, two new clusters at redshifts $z>0.97$, and significantly longer observations of four clusters at $0.6<z<0.9$. Our low-reds… ▽ More We present updated cosmological constraints from measurements of the gas mass fractions ($f_{gas}$) of massive, dynamically relaxed galaxy clusters. Our new data set has greater leverage on models of dark energy, thanks to the addition of the Perseus Cluster at low redshifts, two new clusters at redshifts $z>0.97$, and significantly longer observations of four clusters at $0.6<z<0.9$. Our low-redshift ($z<0.16$) $f_{gas}$ data, combined with the cosmic baryon fraction measured from the cosmic microwave background (CMB), imply a Hubble constant of $h = 0.722 \pm 0.067$. Combining the full $f_{gas}$ data set with priors on the cosmic baryon density and the Hubble constant, we constrain the dark energy density to be $Ω_Λ= 0.865 \pm 0.119$ in non-flat $Λ$CDM (cosmological constant) models, and its equation of state to be $w = -1.13_{-0.20}^{+0.17}$ in flat, constant-w models, respectively 41 and 29 per cent tighter than our previous work, and comparable to the best constraints available from other probes. Combining $f_{gas}$, CMB, supernova, and baryon acoustic oscillation data, we also constrain models with global curvature and evolving dark energy. For the massive, relaxed clusters employed here, we find the scaling of $f_{gas}$ with mass to be consistent with a constant, with an intrinsic scatter that corresponds to just 3 per cent in distance. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: To be published in MNRAS. See https://github.com/abmantz/fgas-cosmo for code and https://github.com/abmantz/fgas-2021-paper for data presented in figures/tables

Lost Horizon: Quantifying the Effect of Local Topography on Global 21-cm Cosmology Data Analysis

Authors: Neil Bassett, David Rapetti, Keith Tauscher, Bang D. Nhan, David D. Bordenave, Joshua J. Hibbard, Jack O. Burns

Abstract: We present an investigation of the horizon and its effect on global 21-cm observations and analysis. We find that the horizon cannot be ignored when modeling low frequency observations. Even if the sky and antenna beam are known exactly, forward models cannot fully describe the beam-weighted foreground component without accurate knowledge of the horizon. When fitting data to extract the 21-cm sign… ▽ More We present an investigation of the horizon and its effect on global 21-cm observations and analysis. We find that the horizon cannot be ignored when modeling low frequency observations. Even if the sky and antenna beam are known exactly, forward models cannot fully describe the beam-weighted foreground component without accurate knowledge of the horizon. When fitting data to extract the 21-cm signal, a single time-averaged spectrum or independent multi-spectrum fits may be able to compensate for the bias imposed by the horizon. However, these types of fits lack constraining power on the 21-cm signal, leading to large uncertainties on the signal extraction, in some cases larger in magnitude than the 21-cm signal itself. A significant decrease in signal uncertainty can be achieved by performing multi-spectrum fits in which the spectra are modeled simultaneously with common parameters. The cost of this greatly increased constraining power, however, is that the time dependence of the horizon's effect, which is more complex than its spectral dependence, must be precisely modeled to achieve a good fit. To aid in modeling the horizon, we present an algorithm and Python package for calculating the horizon profile from a given observation site using elevation data. We also address several practical concerns such as pixelization error, uncertainty in the horizon profile, and foreground obstructions such as surrounding buildings and vegetation. We demonstrate that our training set-based analysis pipeline can account for all of these factors to model the horizon well enough to precisely extract the 21-cm signal from simulated observations. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: 19 pages, 11 figures, 1 table

Global 21-cm signal extraction from foreground and instrumental effects IV: Accounting for realistic instrument uncertainties and their overlap with foreground and signal models

Authors: Keith Tauscher, David Rapetti, Bang D. Nhan, Alec Handy, Neil Bassett, Joshua Hibbard, David Bordenave, Richard F. Bradley, Jack O. Burns

Abstract: All 21-cm signal experiments rely on electronic receivers that affect the data via both multiplicative and additive biases through the receiver's gain and noise temperature. While experiments attempt to remove these biases, the residuals of their imperfect calibration techniques can still confuse signal extraction algorithms. In this paper, the fourth and final installment of our pipeline series,… ▽ More All 21-cm signal experiments rely on electronic receivers that affect the data via both multiplicative and additive biases through the receiver's gain and noise temperature. While experiments attempt to remove these biases, the residuals of their imperfect calibration techniques can still confuse signal extraction algorithms. In this paper, the fourth and final installment of our pipeline series, we present a technique for fitting out receiver effects as efficiently as possible. The fact that the gain and global signal, which are multiplied in the observation equation, must both be modeled implies that the model of the data is nonlinear in its parameters, making numerical sampling the only way to explore the parameter distribution rigorously. However, multi-spectra fits, which are necessary to extract the signal confidently as demonstrated in the third paper of the series, often require large numbers of foreground parameters, increasing the dimension of the posterior distribution that must be explored and therefore causing numerical sampling inefficiencies. Building upon techniques in the second paper of the series, we outline a method to explore the full parameter distribution by numerically sampling a small subset of the parameters and analytically marginalizing over the others. We test this method in simulation using a type-I Chebyshev band-pass filter gain model and a fast signal model based on a spline between local extrema. The method works efficiently, converging quickly to the posterior signal parameter distribution. The final signal uncertainties are of the same order as the noise in the data. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 15 pages, 7 figures, submitted to ApJ

An X-ray and Radio Study of the Hubble Frontier Field Cluster Abell S1063

Authors: Majidul Rahaman, Ramij Raja, Abhirup Datta, Jack O Burns, Brian Alden, David Rapetti

Abstract: We present results from \textit{Chandra} X-ray observations and 325 MHz Giant Metrewave Radio Telescope (GMRT) observations of the massive and X-ray luminous cluster of galaxies Abell S1063. We report the detection of large-scale \lq\lq excess brightness\rq\rq\ in the residual \textit{Chandra} X-ray surface brightness map, which extends at least 2.7 Mpc towards the north-east from the center of th… ▽ More We present results from \textit{Chandra} X-ray observations and 325 MHz Giant Metrewave Radio Telescope (GMRT) observations of the massive and X-ray luminous cluster of galaxies Abell S1063. We report the detection of large-scale \lq\lq excess brightness\rq\rq\ in the residual \textit{Chandra} X-ray surface brightness map, which extends at least 2.7 Mpc towards the north-east from the center of the cluster. We also present a high fidelity X-ray flux and temperature map using \textit{Chandra} archival data of 122 ksec, which shows the disturbed morphology in the cluster. The residual flux map shows the first observational confirmation of the merging axis proposed by earlier simulation by \citet{Gomez2012AJ....144...79G}. The average temperature within $R_{500}$ is $11.7 \pm 0.56$ keV, which makes AS1063 one of the hottest clusters in the nearby Universe. The integrated radio flux density at 325 MHz is found to be $62.0\pm6.3$ mJy. The integrated spectrum of the radio halo follows a power-law with a spectral index $α=-1.43\pm 0.13$. The radio halo is found to be significantly under-luminous, which favored for both the hadronic as well as the turbulent re-acceleration mechanism for its origin. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: Accepted for publication in MNRAS (12 pages, 8 figures, 2 tables)

A Lunar Farside Low Radio Frequency Array for Dark Ages 21-cm Cosmology

Authors: Jack Burns, Gregg Hallinan, Tzu-Ching Chang, Marin Anderson, Judd Bowman, Richard Bradley, Steven Furlanetto, Alex Hegedus, Justin Kasper, Jonathan Kocz, Joseph Lazio, Jim Lux, Robert MacDowall, Jordan Mirocha, Issa Nesnas, Jonathan Pober, Ronald Polidan, David Rapetti, Andres Romero-Wolf, Anže Slosar, Albert Stebbins, Lawrence Teitelbaum, Martin White

Abstract: An array of low-frequency dipole antennas on the lunar farside surface will probe a unique, unexplored epoch in the early Universe called the Dark Ages. It begins at Recombination when neutral hydrogen atoms formed, first revealed by the cosmic microwave background. This epoch is free of stars and astrophysics, so it is ideal to investigate high energy particle processes including dark matter, ear… ▽ More An array of low-frequency dipole antennas on the lunar farside surface will probe a unique, unexplored epoch in the early Universe called the Dark Ages. It begins at Recombination when neutral hydrogen atoms formed, first revealed by the cosmic microwave background. This epoch is free of stars and astrophysics, so it is ideal to investigate high energy particle processes including dark matter, early Dark Energy, neutrinos, and cosmic strings. A NASA-funded study investigated the design of the instrument and the deployment strategy from a lander of 128 pairs of antenna dipoles across a 10 kmx10 km area on the lunar surface. The antenna nodes are tethered to the lander for central data processing, power, and data transmission to a relay satellite. The array, named FARSIDE, would provide the capability to image the entire sky in 1400 channels spanning frequencies from 100 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, the Earth's auroral kilometric radiation, and plasma noise from the solar wind. It is thus the only location within the inner solar system from which sky noise limited observations can be carried out at sub-MHz frequencies. Through precision calibration via an orbiting beacon and exquisite foreground characterization, the farside array would measure the Dark Ages global 21-cm signal at redshifts z~35-200. It will also be a pathfinder for a larger 21-cm power spectrum instrument by carefully measuring the foreground with high dynamic range. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 14 pages, 3 figures, 1 table, response to DOE request for information on lunar farside radio telescope to explore the early universe. arXiv admin note: substantial text overlap with arXiv:1911.08649

Global 21-cm Cosmology from the Farside of the Moon

Authors: Jack Burns, Stuart Bale, Richard Bradley, Z. Ahmed, S. W. Allen, J. Bowman, S. Furlanetto, R. MacDowall, J. Mirocha, B. Nhan, M. Pivovaroff, M. Pulupa, D. Rapetti, A. Slosar, K. Tauscher

Abstract: One of the last unexplored windows to the cosmos, the Dark Ages and Cosmic Dawn, can be opened using a simple low frequency radio telescope from the stable, quiet lunar farside to measure the Global 21-cm spectrum. This frontier remains an enormous gap in our knowledge of the Universe. Standard models of physics and cosmology are untested during this critical epoch. The messenger of information ab… ▽ More One of the last unexplored windows to the cosmos, the Dark Ages and Cosmic Dawn, can be opened using a simple low frequency radio telescope from the stable, quiet lunar farside to measure the Global 21-cm spectrum. This frontier remains an enormous gap in our knowledge of the Universe. Standard models of physics and cosmology are untested during this critical epoch. The messenger of information about this period is the 1420 MHz (21-cm) radiation from the hyperfine transition of neutral hydrogen, Doppler-shifted to low radio astronomy frequencies by the expansion of the Universe. The Global 21-cm spectrum uniquely probes the cosmological model during the Dark Ages plus the evolving astrophysics during Cosmic Dawn, yielding constraints on the first stars, on accreting black holes, and on exotic physics such as dark matter-baryon interactions. A single low frequency radio telescope can measure the Global spectrum between ~10-110 MHz because of the ubiquity of neutral hydrogen. Precise characterizations of the telescope and its surroundings are required to detect this weak, isotropic emission of hydrogen amidst the bright foreground Galactic radiation. We describe how two antennas will permit observations over the full frequency band: a pair of orthogonal wire antennas and a 0.3-m$^3$ patch antenna. A four-channel correlation spectropolarimeter forms the core of the detector electronics. Technology challenges include advanced calibration techniques to disentangle covariances between a bright foreground and a weak 21-cm signal, using techniques similar to those for the CMB, thermal management for temperature swings of >250C, and efficient power to allow operations through a two-week lunar night. This simple telescope sets the stage for a lunar farside interferometric array to measure the Dark Ages power spectrum. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Ensuring Robustness in Training Set Based Global 21-cm Cosmology Analysis

Authors: Neil Bassett, David Rapetti, Keith Tauscher, Jack Burns, Joshua Hibbard

Abstract: We present a methodology for ensuring the robustness of our analysis pipeline in separating the global 21-cm hydrogen cosmology signal from large systematics based on singular value decomposition (SVD) of training sets. We show how traditional goodness-of-fit metrics such as the $χ^2$ statistic that assess the fit to the full data may not be able to detect a suboptimal extraction of the 21-cm sign… ▽ More We present a methodology for ensuring the robustness of our analysis pipeline in separating the global 21-cm hydrogen cosmology signal from large systematics based on singular value decomposition (SVD) of training sets. We show how traditional goodness-of-fit metrics such as the $χ^2$ statistic that assess the fit to the full data may not be able to detect a suboptimal extraction of the 21-cm signal when it is fit alongside one or more additional components due to significant covariance between them. However, we find that comparing the number of SVD eigenmodes for each component chosen by the pipeline for a given fit to the distribution of eigenmodes chosen for synthetic data realizations created from training set curves can detect when one or more of the training sets is insufficient to optimally extract the signal. Furthermore, this test can distinguish which training set (e.g. foreground, 21-cm signal) needs to be modified in order to better describe the data and improve the quality of the 21-cm signal extraction. We also extend this goodness-of-fit testing to cases where a prior distribution derived from the training sets is applied and find that, in this case, the $χ^2$ statistic as well as the recently introduced $ψ^2$ statistic are able to detect inadequacies in the training sets due to the increased restrictions imposed by the prior. Crucially, the tests described in this paper can be performed when analyzing any type of observations with our pipeline. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: 16 pages, 6 figures, 2 tables. Submitted to the Astrophysical Journal

Modelling the Galactic Foreground and Beam Chromaticity for Global 21-cm Cosmology

Authors: Joshua J. Hibbard, Keith Tauscher, David Rapetti, Jack O. Burns

Abstract: In order to characterize and model the beam-weighted foreground for global 21-cm signal experiments, we present a methodology for generating basis eigenvectors that combines analytical and observational models of both the galactic spectral index and sky brightness temperature with simulations of beams having various angular and spectral dependencies and pointings. Each combination creates a unique… ▽ More In order to characterize and model the beam-weighted foreground for global 21-cm signal experiments, we present a methodology for generating basis eigenvectors that combines analytical and observational models of both the galactic spectral index and sky brightness temperature with simulations of beams having various angular and spectral dependencies and pointings. Each combination creates a unique beam-weighted foreground. By generating eigenvectors to fit each foreground model using Singular Value Decomposition (SVD), we examine the effects of varying the components of the beam-weighted foreground. We find that the eigenvectors for modelling an achromatic, isotropic beam -- the ideal case -- are nearly identical regardless of the unweighted foreground model used, and are practicably indistinguishable from polynomial-based models. When anisotropic, chromatic beams weight the foreground, however, a coupling is introduced between the spatial and spectral structure of the foreground which distorts the eigenvectors away from the polynomial models and induces a dependence of the basis upon the exact features of the beam (chromaticity, pattern, pointing) and foreground (spectral index, sky brightness temperature map). We find that the beam has a greater impact upon the eigenvectors than foreground models. Any model which does not account for its distortion may produce RMS uncertainties on the order of $\sim 10$ - $10^3$ Kelvin for six-parameter, single spectrum fits. If the beam is incorporated directly using SVD and training sets, however, the resultant eigenvectors yield milli-Kelvin level uncertainties. Given a sufficiently detailed description of the sky, our methodology can be applied to any particular experiment with a suitably characterized beam for the purpose of generating accurate beam-weighted foreground models. △ Less

Submitted 1 November, 2020; originally announced November 2020.

Comments: 16 pages, 10 figures

Cosmological Constraints from DES Y1 Cluster Abundances and SPT Multi-wavelength data

Authors: M. Costanzi, A. Saro, S. Bocquet, T. M. C. Abbott, M. Aguena, S. Allam, A. Amara, J. Annis, S. Avila, D. Bacon, B. A. Benson, S. Bhargava, D. Brooks, E. Buckley-Geer, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, A. Choi, L. N. da Costa, M. E. S. Pereira, J. De Vicente, S. Desai, H. T. Diehl, J. P. Dietrich , et al. (65 additional authors not shown)

Abstract: We perform a joint analysis of the counts of redMaPPer clusters selected from the Dark Energy Survey (DES) Y1 data and multi-wavelength follow-up data collected within the 2500 deg$^2$ South Pole Telescope (SPT) SZ survey. The SPT follow-up data, calibrating the richness--mass relation of the optically selected redMaPPer catalog, enable the cosmological exploitation of the DES cluster abundance da… ▽ More We perform a joint analysis of the counts of redMaPPer clusters selected from the Dark Energy Survey (DES) Y1 data and multi-wavelength follow-up data collected within the 2500 deg$^2$ South Pole Telescope (SPT) SZ survey. The SPT follow-up data, calibrating the richness--mass relation of the optically selected redMaPPer catalog, enable the cosmological exploitation of the DES cluster abundance data. To explore possible systematics related to the modeling of projection effects, we consider two calibrations of the observational scatter on richness estimates: a simple Gaussian model which account only for the background contamination (BKG), and a model which further includes contamination and incompleteness due to projection effects (PRJ). Assuming either a $Λ$CDM+$\sum m_ν$ or $w$CDM+$\sum m_ν$ cosmology, and for both scatter models, we derive cosmological constraints consistent with multiple cosmological probes of the low and high redshift Universe, and in particular with the SPT cluster abundance data. This result demonstrates that the DES Y1 and SPT cluster counts provide consistent cosmological constraints, if the same mass calibration data set is adopted. It thus supports the conclusion of the DES Y1 cluster cosmology analysis which interprets the tension observed with other cosmological probes in terms of systematics affecting the stacked weak lensing analysis of optically--selected low--richness clusters. Finally, we analyse the first combined optically-SZ selected cluster catalogue obtained by including the SPT sample above the maximum redshift probed by the DES Y1 redMaPPer sample. Besides providing a mild improvement of the cosmological constraints, this data combination serves as a stricter test of our scatter models: the PRJ model, providing scaling relations consistent between the two abundance and multi-wavelength follow-up data, is favored over the BKG model. △ Less

Submitted 15 February, 2022; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: matches published version

Journal ref: Phys. Rev. D 103, 043522 (2021)

Formulating and critically examining the assumptions of global 21-cm signal analyses: How to avoid the false troughs that can appear in single spectrum fits

Authors: Keith Tauscher, David Rapetti, Jack O. Burns

Abstract: The assumptions inherent to global 21-cm signal analyses are rarely delineated. In this paper, we formulate a general list of suppositions underlying a given claimed detection of the global 21-cm signal. Then, we specify the form of these assumptions for two different analyses: 1) the one performed by the EDGES team showing an absorption trough in brightness temperature that they modeled separatel… ▽ More The assumptions inherent to global 21-cm signal analyses are rarely delineated. In this paper, we formulate a general list of suppositions underlying a given claimed detection of the global 21-cm signal. Then, we specify the form of these assumptions for two different analyses: 1) the one performed by the EDGES team showing an absorption trough in brightness temperature that they modeled separately from the sky foreground and 2) a new, so-called Minimum Assumption Analysis (MAA), that makes the most conservative assumptions possible for the signal. We show fits using the EDGES analysis on various beam-weighted foreground simulations from the EDGES latitude with no signal added. Depending on the beam used, these simulations produce large false troughs due to the invalidity of the foreground model to describe the combination of beam chromaticity and the shape of the Galactic plane in the sky, the residuals of which are captured by the ad hoc flattened Gaussian signal model. On the other hand, the MAA provides robust fits by including many spectra at different time bins and allowing any possible 21-cm spectrum to be modeled exactly. We present uncertainty levels and example signal reconstructions found with the MAA for different numbers of time bins. With enough time bins, one can determine the true 21-cm signal with the MAA to $<10$ times the noise level. △ Less

Submitted 5 June, 2020; v1 submitted 30 April, 2020; originally announced May 2020.

Comments: 19 pages, 4 figures, accepted to ApJ. Since previous version, added frequency correlation structure under MAA analysis and laid out extension of MAA to motion-induced dipole of 21-cm signal

Global 21-cm signal extraction from foreground and instrumental effects III: Utilizing drift-scan time dependence and full Stokes measurements

Authors: Keith Tauscher, David Rapetti, Jack O. Burns

Abstract: When using valid foreground and signal models, the uncertainties on extracted signals in global 21-cm signal experiments depend principally on the overlap between signal and foreground models. In this paper, we investigate two strategies for decreasing this overlap: (i) utilizing time dependence by fitting multiple drift-scan spectra simultaneously and (ii) measuring all four Stokes parameters ins… ▽ More When using valid foreground and signal models, the uncertainties on extracted signals in global 21-cm signal experiments depend principally on the overlap between signal and foreground models. In this paper, we investigate two strategies for decreasing this overlap: (i) utilizing time dependence by fitting multiple drift-scan spectra simultaneously and (ii) measuring all four Stokes parameters instead of only the total power, Stokes I. Although measuring polarization requires different instruments than are used in most existing experiments, all existing experiments can utilize drift-scan measurements merely by averaging their data differently. In order to evaluate the increase in constraining power from using these two techniques, we define a method for connecting Root-Mean-Square (RMS) uncertainties to probabilistic confidence levels. Employing simulations, we find that fitting only one total power spectrum leads to RMS uncertainties at the few K level, while fitting multiple time-binned, drift-scan spectra yields uncertainties at the $\lesssim 10$ mK level. This significant improvement only appears if the spectra are modeled with one set of basis vectors, instead of using multiple sets of basis vectors that independently model each spectrum. Assuming that they are simulated accurately, measuring all four Stokes parameters also leads to lower uncertainties. These two strategies can be employed simultaneously and fitting multiple time bins of all four Stokes parameters yields the best precision measurements of the 21-cm signal, approaching the noise level in the data. △ Less

Submitted 23 April, 2020; v1 submitted 11 March, 2020; originally announced March 2020.

Comments: 16 pages, 6 figures, accepted for publication in ApJ

Characterizing the Radio Quiet Region Behind the Lunar Farside for Low Radio Frequency Experiments

Authors: Neil Bassett, David Rapetti, Jack O. Burns, Keith Tauscher, Robert MacDowall

Abstract: Low radio frequency experiments performed on Earth are contaminated by both ionospheric effects and radio frequency interference (RFI) from Earth-based sources. The lunar farside provides a unique environment above the ionosphere where RFI is heavily attenuated by the presence of the Moon. We present electrodynamics simulations of the propagation of radio waves around and through the Moon in order… ▽ More Low radio frequency experiments performed on Earth are contaminated by both ionospheric effects and radio frequency interference (RFI) from Earth-based sources. The lunar farside provides a unique environment above the ionosphere where RFI is heavily attenuated by the presence of the Moon. We present electrodynamics simulations of the propagation of radio waves around and through the Moon in order to characterize the level of attenuation on the farside. The simulations are performed for a range of frequencies up to 100 kHz, assuming a spherical lunar shape with an average, constant density. Additionally, we investigate the role of the topography and density profile of the Moon in the propagation of radio waves and find only small effects on the intensity of RFI. Due to the computational demands of performing simulations at higher frequencies, we propose a model for extrapolating the width of the quiet region above 100 kHz that also takes into account height above the lunar surface as well as the intensity threshold chosen to define the quiet region. This model, which we make publicly available through a Python package, allows the size of the radio quiet region to be easily calculated both in orbit or on the surface, making it directly applicable for lunar satellites as well as surface missions. △ Less

Submitted 11 June, 2020; v1 submitted 6 March, 2020; originally announced March 2020.

Comments: 14 pages, 8 figures, accepted to Advances in Space Research. Additional physical insight such as the skin depth of the lunar material added since previous version. Inset of figure 5 revised to be more legible. Results unchanged

Diffuse radio emission in the galaxy cluster SPT-CL J2031-4037: a steep spectrum intermediate radio halo?

Authors: Ramij Raja, Majidul Rahaman, Abhirup Datta, Jack O. Burns, Huib T. Intema, R. J. van Weeren, Eric J. Hallman, David Rapetti, Surajit Paul

Abstract: The advent of sensitive low frequency radio observations has revealed a number of diffuse radio objects with peculiar properties that are challenging our understanding about the physics of the intracluster medium. Here, we report the discovery of a steep spectrum radio halo surrounding the central Brightest Cluster Galaxy (BCG) in the galaxy cluster SPT-CL J2031-4037. This cluster is morphological… ▽ More The advent of sensitive low frequency radio observations has revealed a number of diffuse radio objects with peculiar properties that are challenging our understanding about the physics of the intracluster medium. Here, we report the discovery of a steep spectrum radio halo surrounding the central Brightest Cluster Galaxy (BCG) in the galaxy cluster SPT-CL J2031-4037. This cluster is morphologically disturbed yet has a weak cool core, an example of cool core/non-cool core transition system, which harbours a radio halo of $\sim 0.7$ Mpc in size. The halo emission detected at 1.7 GHz is less extended compared to that in the 325 MHz observation, and the spectral index of the part of the halo visible at 325 MHz to 1.7 GHz frequencies was found to be $-1.35 \pm 0.07$. Also, $P_{1.4\ \mathrm{GHz}}$ was found to be $0.77 \times 10^{24}$ W Hz$^{-1}$ which falls in the region where radio mini-halos, halo upper limits and ultra-steep spectrum (USS) halos are found in the $P_{1.4\ \mathrm{GHz}} - L_\mathrm{X}$ plane. Additionally, simulations presented in the paper provide support to the scenario of the steep spectrum. The diffuse radio emission found in this cluster may be a steep spectrum "intermediate" or "hybrid" radio halo which is transitioning into a mini-halo. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 6 pages, 3 figures; Accepted for publication in MNRAS Letter

Probing the origin of diffuse radio emission in the cool-core of the Phoenix galaxy cluster

Authors: Ramij Raja, Majidul Rahaman, Abhirup Datta, Jack O. Burns, Brian Alden, H. T. Intema, R. J. van Weeren, Eric J. Hallman, David Rapetti, Surajit Paul

Abstract: Cool core galaxy clusters are considered to be dynamically relaxed clusters with regular morphology and highly X-ray luminous central region. However, cool core clusters can also be sites for merging events that exhibit cold fronts in X-ray and mini-halos in radio. We present recent radio/X-ray observations of the Phoenix Cluster or SPT-CL J2344-4243 at the redshift of $z=0.596$. Using archival {\… ▽ More Cool core galaxy clusters are considered to be dynamically relaxed clusters with regular morphology and highly X-ray luminous central region. However, cool core clusters can also be sites for merging events that exhibit cold fronts in X-ray and mini-halos in radio. We present recent radio/X-ray observations of the Phoenix Cluster or SPT-CL J2344-4243 at the redshift of $z=0.596$. Using archival {\it Chandra} X-ray observations, we detect spiraling cool gas around the cluster core as well as discover two cold fronts near the core. It is perhaps the most distant galaxy cluster to date known to host cold fronts. Also, we present JVLA\footnote{Jansky Very Large Array\\ \url{https://science.nrao.edu/facilities/vla}} 1.52 GHz observations of the minihalo, previously discovered at 610 MHz with GMRT\footnote{Giant Metrewave Radio Telescope\\ \url{http://www.gmrt.ncra.tifr.res.in}} observations in the center of the Phoenix galaxy cluster. The minihalo flux density at 1.52 GHz is $9.65 \pm 0.97$ mJy with the spectral index between 610 MHz and 1.52 GHz being $-0.98 \pm 0.16$\footnote{$S_ν = ν^α$ where $S_ν$}. A possible origin of these radio sources is turbulence induced by sloshing of the gas in the cluster core. △ Less

Submitted 17 December, 2019; originally announced December 2019.

Comments: 13 pages, 6 figures, accepted by ApJ

XXL Survey groups and clusters in the Hyper Suprime-Cam Survey. Scaling relations between X-ray properties and weak lensing mass

Authors: Mauro Sereno, Keiichi Umetsu, Stefano Ettori, Dominique Eckert, Fabio Gastaldello, Paul Giles, Maggie Lieu, Ben Maughan, Nobuhiro Okabe, Mark Birkinshaw, I-Non Chiu, Yutaka Fujita, Satoshi Miyazaki, David Rapetti, Elias Koulouridis, Marguerite Pierre

Abstract: Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak l… ▽ More Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak lensing mass show imprints of radiative cooling and AGN feedback in groups. From the multi-variate analysis, we found some evidence for steeper than self-similar slopes for gas mass ($β_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($β_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($β_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $ρ_{X_1X_2|m}\sim0.34$) due to dynamical state and merger history of the halos. Positive correlations with the weak lensing mass (median correlation factor $ρ_{m_\text{wl}X|m}\sim0.35$) can be connected to triaxiality and orientation. Comparison of weak lensing and hydrostatic masses suggests a small role played by non-thermal pressure support ($9\pm17\%$). △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: 19 pages; in press on MNRAS. One of two companion papers on HSC-XXL surveys (Keiichi Umetsu et al., arXiv:1909.10524)

Global 21-cm Signal Extraction from Foreground and Instrumental Effects II: Efficient and Self-Consistent Technique for Constraining Nonlinear Signal Models

Authors: David Rapetti, Keith Tauscher, Jordan Mirocha, Jack O. Burns

Abstract: We present the completion of a data analysis pipeline that self-consistently separates global 21-cm signals from large systematics using a pattern recognition technique. In the first paper of this series, we obtain optimal basis vectors from signal and foreground training sets to linearly fit both components with the minimal number of terms that best extracts the signal given its overlap with the… ▽ More We present the completion of a data analysis pipeline that self-consistently separates global 21-cm signals from large systematics using a pattern recognition technique. In the first paper of this series, we obtain optimal basis vectors from signal and foreground training sets to linearly fit both components with the minimal number of terms that best extracts the signal given its overlap with the foreground. In this second paper, we utilize the spectral constraints derived in the first paper to calculate the full posterior probability distribution of any signal parameter space of choice. The spectral fit provides the starting point for a Markov Chain Monte Carlo (MCMC) engine that samples the signal without traversing the foreground parameter space. At each MCMC step, we marginalize over the weights of all linear foreground modes and suppress those with unimportant variations by applying priors gleaned from the training set. This method drastically reduces the number of MCMC parameters, augmenting the efficiency of exploration, circumvents the need for selecting a minimal number of foreground modes, and allows the complexity of the foreground model to be greatly increased to simultaneously describe many observed spectra without requiring extra MCMC parameters. Using two nonlinear signal models, one based on EDGES observations and the other on phenomenological frequencies and temperatures of theoretically expected extrema, we demonstrate the success of this methodology by recovering the input parameters from multiple randomly simulated signals at low radio frequencies (10-200 MHz), while rigorously accounting for realistically modeled beam-weighted foregrounds. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 19 pages, 11 figures, 5 tables. Submitted to ApJ. Code available at https://bitbucket.org/ktausch/pylinex

NASA Probe Study Report: Farside Array for Radio Science Investigations of the Dark ages and Exoplanets (FARSIDE)

Authors: Jack O. Burns, Gregg Hallinan, Jim Lux, Lawrence Teitelbaum, Jonathon Kocz, Robert MacDowall, Richard Bradley, David Rapetti, Wenbo Wu, Steven Furlanetto, Alex Austin, Andres Romero-Wolf, Tzu-Ching Chang, Judd Bowman, Justin Kasper, Marin Anderson, Zhongwen Zhen, Jonathan Pober, Jordan Mirocha

Abstract: This is the final report submitted to NASA for a Probe-class concept study of the "Farside Array for Radio Science Investigations of the Dark ages and Exoplanets" (FARSIDE), a low radio frequency interferometric array on the farside of the Moon. The design study focused on the instrument, a deployment rover, the lander and base station, and delivered an architecture broadly consistent with the req… ▽ More This is the final report submitted to NASA for a Probe-class concept study of the "Farside Array for Radio Science Investigations of the Dark ages and Exoplanets" (FARSIDE), a low radio frequency interferometric array on the farside of the Moon. The design study focused on the instrument, a deployment rover, the lander and base station, and delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dipole antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway, or an alternative relay satellite. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 150 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts from 50-100. It will also be a pathfinder for a larger 21-cm power spectrum instrument by carefully measuring the foreground with high dynamic range. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: 50 pages, NASA Probe final study report. arXiv admin note: text overlap with arXiv:1907.05407

Algorithms and Statistical Models for Scientific Discovery in the Petabyte Era

Authors: Brian Nord, Andrew J. Connolly, Jamie Kinney, Jeremy Kubica, Gautaum Narayan, Joshua E. G. Peek, Chad Schafer, Erik J. Tollerud, Camille Avestruz, G. Jogesh Babu, Simon Birrer, Douglas Burke, João Caldeira, Douglas A. Caldwell, Joleen K. Carlberg, Yen-Chi Chen, Chuanfei Dong, Eric D. Feigelson, V. Zach Golkhou, Vinay Kashyap, T. S. Li, Thomas Loredo, Luisa Lucie-Smith, Kaisey S. Mandel, J. R. Martínez-Galarza , et al. (13 additional authors not shown)

Abstract: The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our… ▽ More The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/). △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1905.05116

Report number: FERMILAB-FN-1093-A-AE-SCD

Weak lensing Analysis of X-Ray-selected XXL Galaxy Groups and Clusters with Subaru HSC Data

Authors: Keiichi Umetsu, Mauro Sereno, Maggie Lieu, Hironao Miyatake, Elinor Medezinski, Atsushi J. Nishizawa, Paul Giles, Fabio Gastaldello, Ian G. McCarthy, Martin Kilbinger, Mark Birkinshaw, Stefano Ettori, Nobuhiro Okabe, I-Non Chiu, Jean Coupon, Dominique Eckert, Yutaka Fujita, Yuichi Higuchi, Elias Koulouridis, Ben Maughan, Satoshi Miyazaki, Masamune Oguri, Florian Pacaud, Marguerite Pierre, David Rapetti , et al. (1 additional authors not shown)

Abstract: We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of in… ▽ More We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be β= -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as σ(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5σ$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures. △ Less

Submitted 4 March, 2020; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: Version matching the one published in ApJ. We recommend to use statistically corrected mass estimates (M200MT, M500MT) of Table 2 for a given individual cluster. One of two companion papers presenting initial HSC-XXL results (Mauro Sereno et al., arXiv:1912.02827)

Journal ref: ApJ, 890, 148 (2020)

Research and Development for HI Intensity Mapping

Authors: Zeeshan Ahmed, David Alonso, Mustafa A. Amin, Réza Ansari, Evan J. Arena, Kevin Bandura, Adam Beardsley, Philip Bull, Emanuele Castorina, Tzu-Ching Chang, Romeel Davé, Joshua S. Dillon, Alexander van Engelen, Aaron Ewall-Wice, Simone Ferraro, Simon Foreman, Josef Frisch, Daniel Green, Gilbert Holder, Daniel Jacobs, Dionysios Karagiannis, Alexander A. Kaurov, Lloyd Knox, Emily Kuhn, Adrian Liu , et al. (29 additional authors not shown)

Abstract: Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and an… ▽ More Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift. △ Less

Submitted 29 July, 2019; originally announced July 2019.

Comments: 10 pages + references, 2 figures, 1 table; APC white paper submitted to the Astro2020 decadal survey. arXiv admin note: text overlap with arXiv:1810.09572

A Detection of CMB-Cluster Lensing using Polarization Data from SPTpol

Authors: S. Raghunathan, S. Patil, E. Baxter, B. A. Benson, L. E. Bleem, T. M. Crawford, G. P. Holder, T. McClintock, C. L. Reichardt, T. N. Varga, N. Whitehorn, P. A. R. Ade, S. Allam, A. J. Anderson, J. E. Austermann, S. Avila, J. S. Avva, D. Bacon, J. A. Beall, A. N. Bender, F. Bianchini, S. Bocquet, D. Brooks, D. L. Burke, J. E. Carlstrom , et al. (104 additional authors not shown)

Abstract: We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization… ▽ More We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization gradient. Using data from the SPTpol 500 deg$^{2}$ survey at the locations of roughly 18,000 clusters with richness $λ\ge 10$ from the Dark Energy Survey (DES) Year-3 full galaxy cluster catalog, we detect lensing at $4.8σ$. The mean stacked mass of the selected sample is found to be $(1.43 \pm 0.4)\ \times 10^{14}\ {\rm M_{\odot}}$ which is in good agreement with optical weak lensing based estimates using DES data and CMB-lensing based estimates using SPTpol temperature data. This measurement is a key first step for cluster cosmology with future low-noise CMB surveys, like CMB-S4, for which CMB polarization will be the primary channel for cluster lensing measurements. △ Less

Submitted 24 September, 2019; v1 submitted 19 July, 2019; originally announced July 2019.

Comments: 10 pages, 3 figures, 1 table; typos fixed; accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 123, 181301 (2019)

FARSIDE: A Low Radio Frequency Interferometric Array on the Lunar Farside

Authors: Jack Burns, Gregg Hallinan, Jim Lux, Andres Romero-Wolf, Tzu-Ching Chang, Jonathan Kocz, Judd Bowman, Robert MacDowall, Justin Kasper, Richard Bradley, Marin Anderson, David Rapetti

Abstract: FARSIDE (Farside Array for Radio Science Investigations of the Dark ages and Exoplanets) is a Probe-class concept to place a low radio frequency interferometric array on the farside of the Moon. A NASA-funded design study, focused on the instrument, a deployment rover, the lander and base station, delivered an architecture broadly consistent with the requirements for a Probe mission. This notional… ▽ More FARSIDE (Farside Array for Radio Science Investigations of the Dark ages and Exoplanets) is a Probe-class concept to place a low radio frequency interferometric array on the farside of the Moon. A NASA-funded design study, focused on the instrument, a deployment rover, the lander and base station, delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dual polarization antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 100 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets, including the hypothetical Planet Nine. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts z=50-100. The unique observational window offered by FARSIDE would enable an abundance of additional science ranging from sounding of the lunar subsurface to characterization of the interstellar medium in the solar system neighborhood. △ Less

Submitted 11 July, 2019; originally announced July 2019.

Comments: 13 pages, 9 figures, 2 tables, Astro2020 Decadal Survey whitepaper

The Galaxy Cluster 'Pypeline' for X-ray Temperature Maps: ClusterPyXT

Authors: Brian Alden, Eric J. Hallman, David Rapetti, Jack O. Burns, Abhirup Datta

Abstract: ClusterPyXT is a new software pipeline to generate spectral temperature, X-ray surface brightness, pressure, and density maps from X-ray observations of galaxy clusters. These data products help elucidate the physics of processes occurring within clusters of galaxies, including turbulence, shock fronts, nonthermal phenomena, and the overall dynamics of cluster mergers. ClusterPyXT automates the cr… ▽ More ClusterPyXT is a new software pipeline to generate spectral temperature, X-ray surface brightness, pressure, and density maps from X-ray observations of galaxy clusters. These data products help elucidate the physics of processes occurring within clusters of galaxies, including turbulence, shock fronts, nonthermal phenomena, and the overall dynamics of cluster mergers. ClusterPyXT automates the creation of these data products with minimal user interaction, and allows for rapid analyses of archival data with user defined parameters and the ability to straightforwardly incorporate additional observations. In this paper, we describe in detail the use of this code and release it as an open source Python project on GitHub. △ Less

Submitted 3 June, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: 11 pages, 7 figures, Published in Astronomy and Computing, Volume 27, April 2019

Journal ref: Astronomy and Computing, Volume 27, April 2019, Pages 147-155

Tests of Gravity with Galaxy Clusters

Authors: Matteo Cataneo, David Rapetti

Abstract: Changes in the law of gravity have far-reaching implications for the formation and evolution of galaxy clusters, and appear as peculiar signatures in their mass-observable relations, structural properties, internal dynamics, and abundance. We review the outstanding progress made in recent years towards constraining deviations from General Relativity with galaxy clusters, and give an overview of th… ▽ More Changes in the law of gravity have far-reaching implications for the formation and evolution of galaxy clusters, and appear as peculiar signatures in their mass-observable relations, structural properties, internal dynamics, and abundance. We review the outstanding progress made in recent years towards constraining deviations from General Relativity with galaxy clusters, and give an overview of the yet untapped information becoming accessible with forthcoming surveys that will map large portions of the sky in great detail and unprecedented depth. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 45 pages, 11 figures. This is a pre-referee version of the article published as part of the special issue "Modified Gravity: Progresses and Outlook of Theories, Numerical Techniques and Observational Tests". To read its accepted version and to access the other review articles of the special issue visit https://www.worldscientific.com/toc/ijmpd/27/15

Journal ref: International Journal of Modern Physics D, Vol. 27, No. 15, 1848006 (2018)

Dark Cosmology: Investigating Dark Matter & Exotic Physics in the Dark Ages using the Redshifted 21-cm Global Spectrum

Authors: Jack O. Burns, S. Bale, N. Bassett, J. Bowman, R. Bradley, A. Fialkov, S. Furlanetto, M. Hecht, M. Klein-Wolt, C. Lonsdale, R. MacDowall, J. Mirocha, Julian B. Munoz, B. Nhan, J. Pober, D. Rapetti, A. Rogers, K. Tauscher

Abstract: The Dark Ages, probed by the redshifted 21-cm signal, is the ideal epoch for a new rigorous test of the standard LCDM cosmological model. Divergences from that model would indicate new physics, such as dark matter decay (heating) or baryonic cooling beyond that expected from adiabatic expansion of the Universe. In the early Universe, most of the baryonic matter was in the form of neutral hydrogen… ▽ More The Dark Ages, probed by the redshifted 21-cm signal, is the ideal epoch for a new rigorous test of the standard LCDM cosmological model. Divergences from that model would indicate new physics, such as dark matter decay (heating) or baryonic cooling beyond that expected from adiabatic expansion of the Universe. In the early Universe, most of the baryonic matter was in the form of neutral hydrogen (HI), detectable via its ground state's spin-flip transition. A measurement of the redshifted 21-cm spectrum maps the history of the HI gas through the Dark Ages and Cosmic Dawn and up to the Epoch of Reionization (EoR). The Experiment to Detect the Global EoR Signature (EDGES) recently reported an absorption trough at 78 MHz (redshift z of 17), similar in frequency to expectations for Cosmic Dawn, but about 3 times deeper than was thought possible from standard cosmology and adiabatic cooling of HI. Interactions between baryons and slightly-charged dark matter particles with electron-like mass provide a potential explanation of this difference but other cooling mechanisms are also being investigated to explain these results. The Cosmic Dawn trough is affected by cosmology and the complex astrophysical history of the first luminous objects. Another trough is expected during the Dark Ages, prior to the formation of the first stars and thus determined entirely by cosmological phenomena (including dark matter). Observations on or in orbit above the Moon's farside can investigate this pristine epoch (15-40 MHz; z=100-35), which is inaccessible from Earth. A single cross-dipole antenna or compact array can measure the amplitude of the 21-cm spectrum to the level required to distinguish (at >5$σ$}) the standard cosmological model from that of additional cooling derived from current EDGES results. This observation constitutes a powerful, clean probe of exotic physics in the Dark Ages. △ Less

Submitted 16 February, 2019; originally announced February 2019.

Comments: 7 pages, 6 figures, Science whitepaper submitted to the Astro2020 Decadal Survey

Cluster Cosmology Constraints from the 2500 deg$^2$ SPT-SZ Survey: Inclusion of Weak Gravitational Lensing Data from Magellan and the Hubble Space Telescope

Authors: S. Bocquet, J. P. Dietrich, T. Schrabback, L. E. Bleem, M. Klein, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, M. Bayliss, B. A. Benson, M. Brodwin, E. Bulbul, R. E. A. Canning, R. Capasso, J. E. Carlstrom, C. L. Chang, I. Chiu, H-M. Cho, A. Clocchiatti, T. M. Crawford, A. T. Crites, T. de Haan, S. Desai, M. A. Dobbs , et al. (55 additional authors not shown)

Abstract: We derive cosmological constraints using a galaxy cluster sample selected from the 2500~deg$^2$ SPT-SZ survey. The sample spans the redshift range $0.25< z<1.75$ and contains 343 clusters with SZ detection significance $ξ>5$. The sample is supplemented with optical weak gravitational lensing measurements of 32 clusters with $0.29<z<1.13$ (from Magellan and HST) and X-ray measurements of 89 cluster… ▽ More We derive cosmological constraints using a galaxy cluster sample selected from the 2500~deg$^2$ SPT-SZ survey. The sample spans the redshift range $0.25< z<1.75$ and contains 343 clusters with SZ detection significance $ξ>5$. The sample is supplemented with optical weak gravitational lensing measurements of 32 clusters with $0.29<z<1.13$ (from Magellan and HST) and X-ray measurements of 89 clusters with $0.25<z<1.75$ (from Chandra). We rely on minimal modeling assumptions: i) weak lensing provides an accurate means of measuring halo masses, ii) the mean SZ and X-ray observables are related to the true halo mass through power-law relations in mass and dimensionless Hubble parameter $E(z)$ with a-priori unknown parameters, iii) there is (correlated, lognormal) intrinsic scatter and measurement noise relating these observables to their mean relations. We simultaneously fit for these astrophysical modeling parameters and for cosmology. Assuming a flat $νΛ$CDM model, in which the sum of neutrino masses is a free parameter, we measure $Ω_\mathrm{m}=0.276\pm0.047$, $σ_8=0.781\pm0.037$, and $σ_8(Ω_\mathrm{m}/0.3)^{0.2}=0.766\pm0.025$. The redshift evolution of the X-ray $Y_\mathrm{X}$-mass and $M_\mathrm{gas}$-mass relations are both consistent with self-similar evolution to within $1σ$. The mass-slope of the $Y_\mathrm{X}$-mass relation shows a $2.3σ$ deviation from self-similarity. Similarly, the mass-slope of the $M_\mathrm{gas}$-mass relation is steeper than self-similarity at the $2.5σ$ level. In a $νw$CDM cosmology, we measure the dark energy equation of state parameter $w=-1.55\pm0.41$ from the cluster data. We perform a measurement of the growth of structure since redshift $z\sim1.7$ and find no evidence for tension with the prediction from General Relativity. We provide updated redshift and mass estimates for the SPT sample. (abridged) △ Less

Submitted 20 May, 2019; v1 submitted 4 December, 2018; originally announced December 2018.

Comments: Accepted for publication in ApJ (v2 is accepted version), the catalog can be found at https://pole.uchicago.edu/public/data/sptsz-clusters/

Measurement of the Splashback Feature around SZ-selected Galaxy Clusters with DES, SPT and ACT

Authors: T. Shin, S. Adhikari, E. J. Baxter, C. Chang, B. Jain, N. Battaglia, L. Bleem, S. Bocquet, J. DeRose, D. Gruen, M. Hilton, A. Kravtsov, T. McClintock, E. Rozo, E. S. Rykoff, T. N. Varga, R. H. Wechsler, H. Wu, S. Aiola, S. Allam, K. Bechtol, B. A. Benson, E. Bertin, J. R. Bond, M. Brodwin , et al. (85 additional authors not shown)

Abstract: We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the ob… ▽ More We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these clusters with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases than optically selected clusters. We find that the measured $r_{\rm sp}$ for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters, $r_{\rm sp}$ is $\sim$ $2σ$ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogs and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy color, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster. △ Less

Submitted 24 May, 2019; v1 submitted 14 November, 2018; originally announced November 2018.

Comments: 17 pages, 13 figures, published in MNRAS

Assessment of the Projection-induced Polarimetry Technique for Constraining the Foreground Spectrum in Global 21 cm Cosmology

Authors: Bang D. Nhan, David D. Bordenave, Richard F. Bradley, Jack O. Burns, Keith Tauscher, David Rapetti, Patricia J. Klima

Abstract: Detecting the cosmological sky-averaged (global) 21 cm signal as a function of observed frequency will provide a powerful tool to study the ionization and thermal history of the intergalactic medium (IGM) in the early Universe ($\sim$ 400 million years after the Big Bang). The greatest challenge in conventional total-power global 21 cm experiments is the removal of the foreground synchrotron emiss… ▽ More Detecting the cosmological sky-averaged (global) 21 cm signal as a function of observed frequency will provide a powerful tool to study the ionization and thermal history of the intergalactic medium (IGM) in the early Universe ($\sim$ 400 million years after the Big Bang). The greatest challenge in conventional total-power global 21 cm experiments is the removal of the foreground synchrotron emission ($\sim 10^3$-$10^4$ K) to uncover the weak cosmological signal (tens to hundreds of mK), especially since the intrinsic smoothness of the foreground spectrum is corrupted by instrumental effects. Although the EDGES team has recently reported an absorption profile at 78 MHz in the sky-averaged spectrum, it is necessary to confirm this detection with an independent approach. The projection effect from observing anisotropic foreground source emission with a wide-view antenna pointing at the North Celestial Pole (NCP) can induce a net polarization, referred as the Projection-Induced Polarization Effect (PIPE). Due to Earth's rotation, observation centered at the circumpolar region will impose a dynamic sky modulation on the net polarization's waveforms which is unique to the foreground component. In this study, we review the implementation practicality and underlying instrumental effects of this new polarimetry-based technique with detailed numerical simulation and a testbed instrument, the Cosmic Twilight Polarimeter (CTP). In addition, we explore an SVD-based analysis approach for separating the foreground and instrumental effects from the background global 21 cm signal using the sky-modulated PIPE. △ Less

Submitted 5 August, 2019; v1 submitted 12 November, 2018; originally announced November 2018.

Comments: 19 pages, 24 figures, 1 table, accepted by ApJ

Mass Calibration of Optically Selected DES clusters using a Measurement of CMB-Cluster Lensing with SPTpol Data

Authors: S. Raghunathan, S. Patil, E. Baxter, B. A. Benson, L. E. Bleem, T. L. Chou, T. M. Crawford, G. P. Holder, T. McClintock, C. L. Reichardt, E. Rozo, T. N. Varga, T. M. C. Abbott, P. A. R. Ade, S. Allam, A. J. Anderson, J. Annis, J. E. Austermann, S. Avila, J. A. Beall, K. Bechtol, A. N. Bender, G. Bernstein, E. Bertin, F. Bianchini , et al. (107 additional authors not shown)

Abstract: We use cosmic microwave background (CMB) temperature maps from the 500 deg$^{2}$ SPTpol survey to measure the stacked lensing convergence of galaxy clusters from the Dark Energy Survey (DES) Year-3 redMaPPer (RM) cluster catalog. The lensing signal is extracted through a modified quadratic estimator designed to be unbiased by the thermal Sunyaev-Zel{'}dovich (tSZ) effect. The modified estimator us… ▽ More We use cosmic microwave background (CMB) temperature maps from the 500 deg$^{2}$ SPTpol survey to measure the stacked lensing convergence of galaxy clusters from the Dark Energy Survey (DES) Year-3 redMaPPer (RM) cluster catalog. The lensing signal is extracted through a modified quadratic estimator designed to be unbiased by the thermal Sunyaev-Zel{'}dovich (tSZ) effect. The modified estimator uses a tSZ-free map, constructed from the SPTpol 95 and 150 GHz datasets, to estimate the background CMB gradient. For lensing reconstruction, we employ two versions of the RM catalog: a flux-limited sample containing 4003 clusters and a volume-limited sample with 1741 clusters. We detect lensing at a significance of 8.7$σ$(6.7$σ$) with the flux(volume)-limited sample. By modeling the reconstructed convergence using the Navarro-Frenk-White profile, we find the average lensing masses to be $M_{200m}$ = ($1.62^{+0.32}_{-0.25}$ [stat.] $\pm$ 0.04 [sys.]) and ($1.28^{+0.14}_{-0.18}$ [stat.] $\pm$ 0.03 [sys.]) $\times\ 10^{14}\ M_{\odot}$ for the volume- and flux-limited samples respectively. The systematic error budget is much smaller than the statistical uncertainty and is dominated by the uncertainties in the RM cluster centroids. We use the volume-limited sample to calibrate the normalization of the mass-richness scaling relation, and find a result consistent with the galaxy weak-lensing measurements from DES (Mcclintock et al. 2018). △ Less

Submitted 20 February, 2019; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: 19 pages, 6 figures, published in ApJ

A Ground Plane Artifact that Induces an Absorption Profile in Averaged Spectra from Global 21-cm Measurements - with Possible Application to EDGES

Authors: Richard F. Bradley, Keith Tauscher, David Rapetti, Jack O. Burns

Abstract: Most of the current Global 21-cm experiments include ground screens that help moderate effects from the Earth. In this paper, we report on a possible systematic artifact within the ground plane that may produce broad absorption features in the spectra observed by these experiments. Using analytical approximations and numerical modeling, the origin of the artifact and its impact on the sky-averaged… ▽ More Most of the current Global 21-cm experiments include ground screens that help moderate effects from the Earth. In this paper, we report on a possible systematic artifact within the ground plane that may produce broad absorption features in the spectra observed by these experiments. Using analytical approximations and numerical modeling, the origin of the artifact and its impact on the sky-averaged spectrum are described. The publicly released EDGES dataset, from which a 78 MHz absorption feature was recently suggested, is used to probe for the potential presence of ground plane resonances. While the lack of a noise level for the EDGES spectrum makes traditional goodness-of-fit statistics unattainable, the rms residual can be used to assess the relative goodness of fits performed under similar circumstances. The fit to the EDGES spectrum using a model with a simple 2-term foreground and three cavity-mode resonances is compared to a fit to the same spectrum with a model used by the EDGES team consisting of a 5-term foreground and a flattened Gaussian signal. The fits with the physically motivated resonance and empirical flattened Gaussian models have rms residuals of 20.8 mK (11 parameters) and 24.5 mK (9 parameters), respectively, allowing us to conclude that ground plane resonances constitute another plausible explanation for the EDGES data. △ Less

Submitted 5 March, 2019; v1 submitted 21 October, 2018; originally announced October 2018.

Comments: 12 pages, 9 figures, accepted to ApJ

The XXL Survey XX: The 365 cluster catalogue

Authors: C. Adami, P. Giles, E. Koulouridis, F. Pacaud, C. A. Caretta, M. Pierre, D. Eckert, M. E. Ramos-Ceja, F. Gastaldello, S. Fotopoulou, V. Guglielmo, C. Lidman, T. Sadibekova, A. Iovino, B. Maughan, L. Chiappetti, S. Alis, B. Altieri, I. Baldry, D. Bottini, M. Birkinshaw, M. Bremer, M. Brown, O. Cucciati, S. Driver , et al. (32 additional authors not shown)

Abstract: In the currently debated context of using clusters of galaxies as cosmological probes, the need for well-defined cluster samples is critical. The XXL Survey has been specifically designed to provide a well characterised sample of some 500 X-ray detected clusters suitable for cosmological studies. The main goal of present article is to make public and describe the properties of the cluster catalogu… ▽ More In the currently debated context of using clusters of galaxies as cosmological probes, the need for well-defined cluster samples is critical. The XXL Survey has been specifically designed to provide a well characterised sample of some 500 X-ray detected clusters suitable for cosmological studies. The main goal of present article is to make public and describe the properties of the cluster catalogue in its present state, as well as of associated catalogues as super-clusters and fossil groups. We release a sample containing 365 clusters in total. We give the details of the follow-up observations and explain the procedure adopted to validate the cluster spectroscopic redshifts. Considering the whole XXL cluster sample, we have provided two types of selection, both complete in a particular sense: one based on flux-morphology criteria, and an alternative based on the [0.5-2] keV flux within one arcmin of the cluster centre. We have also provided X-ray temperature measurements for 80$\%$ of the clusters having a flux larger than 9$\times$10$^{-15}$$\rm \thinspace erg \, s^{-1} \, cm^{-2}$. Our cluster sample extends from z$\sim$0 to z$\sim$1.2, with one cluster at z$\sim$2. Clusters were identified through a mean number of six spectroscopically confirmed cluster members. Our updated luminosity function and luminosity-temperature relation are compatible with our previous determinations based on the 100 brightest clusters, but show smaller uncertainties. We also present an enlarged list of super-clusters and a sample of 18 possible fossil groups. This intermediate publication is the last before the final release of the complete XXL cluster catalogue when the ongoing C2 cluster spectroscopic follow-up is complete. It provides a unique inventory of medium-mass clusters over a 50~\dd\ area out to z$\sim$1. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: A&A, in press Shortened abstract Full catalogues available online, see paper

A new goodness-of-fit statistic and its application to 21-cm cosmology

Authors: Keith Tauscher, David Rapetti, Jack O. Burns

Abstract: The reduced chi-squared statistic is a commonly used goodness-of-fit measure, but it cannot easily detect features near the noise level, even when a large amount of data is available. In this paper, we introduce a new goodness-of-fit measure that we name the reduced psi-squared statistic. It probes the two-point correlations in the residuals of a fit, whereas chi-squared accounts for only the abso… ▽ More The reduced chi-squared statistic is a commonly used goodness-of-fit measure, but it cannot easily detect features near the noise level, even when a large amount of data is available. In this paper, we introduce a new goodness-of-fit measure that we name the reduced psi-squared statistic. It probes the two-point correlations in the residuals of a fit, whereas chi-squared accounts for only the absolute values of each residual point, not considering the relationship between these points. The new statistic maintains sensitivity to individual outliers, but is superior to chi-squared in detecting wide, low level features in the presence of a large number of noisy data points. After presenting this new statistic, we show an instance of its use in the context of analyzing radio spectroscopic data for 21-cm cosmology experiments. We perform fits to simulated data with four components: foreground emission, the global 21-cm signal, an idealized instrument systematic, and noise. This example is particularly timely given the ongoing efforts to confirm the first observational result for this signal, where this work found its original motivation. In addition, we release a Python script dubbed $\texttt{psipy}$ which allows for quick, efficient calculation of the reduced psi-squared statistic on arbitrary data arrays, to be applied in any field of study. △ Less

Submitted 20 November, 2018; v1 submitted 28 September, 2018; originally announced October 2018.

Comments: 22 pages, 10 figures, submitted to JCAP, psipy code to calculate new statistic available at https://bitbucket.org/ktausch/psipy

Galaxy populations in the most distant SPT-SZ clusters - I. Environmental quenching in massive clusters at $1.4\lesssim z\lesssim1.7$

Authors: V. Strazzullo, M. Pannella, J. J. Mohr, A. Saro, M. L. N. Ashby, M. B. Bayliss, S. Bocquet, E. Bulbul, G. Khullar, A. B. Mantz, S. A. Stanford, B. A. Benson, L. E. Bleem, M. Brodwin, R. E. A. Canning, R. Capasso, I. Chiu, A. H. Gonzalez, N. Gupta, J. Hlavacek-Larrondo, M. Klein, M. McDonald, E. Noordeh, D. Rapetti, C. Reichardt , et al. (3 additional authors not shown)

Abstract: We present first results from a galaxy population study in the highest redshift galaxy clusters identified in the 2500 deg$^2$ South Pole Telescope Sunyaev Zel'dovich effect (SPT-SZ) survey. The cluster selection is to first order independent of galaxy properties, making the SPT-SZ sample particularly well suited for cluster galaxy population studies. We carry out a 4-band imaging campaign with th… ▽ More We present first results from a galaxy population study in the highest redshift galaxy clusters identified in the 2500 deg$^2$ South Pole Telescope Sunyaev Zel'dovich effect (SPT-SZ) survey. The cluster selection is to first order independent of galaxy properties, making the SPT-SZ sample particularly well suited for cluster galaxy population studies. We carry out a 4-band imaging campaign with the {\it Hubble} and {\it Spitzer} Space Telescopes of the five $z\gtrsim 1.4$, S/N$_{SZE}>$5 clusters, that are among the rarest most massive clusters known at this redshift. All five show clear overdensities of red galaxies whose colors agree with the initial cluster redshift estimates. The highest redshift cluster in this sample, SPT-CLJ0459-4947 at $z\sim1.72$, is the most distant $M_{500}>10^{14}~M_{\odot}$ ICM-selected cluster discovered thus far, and is one of only three known clusters in this mass range at $z\gtrsim 1.7$, regardless of selection. Based on UVJ-like photometric classification of quiescent and star-forming galaxies, the passive fraction in the cluster central regions ($r/r_{500}<0.7$) is higher than in the field at the same redshift, with corresponding environmental quenching efficiencies typically in the range $\sim0.5-0.8$ for stellar masses $\log(M/M_{\odot})>10.85$. We have explored the impact of emission from star formation on the selection of this sample, concluding that all five clusters studied here would still have been detected with S/N$_{SZE}>$5, even if they had the same passive fraction as measured in the field. Our results thus point towards an efficient suppression of star formation in massive galaxies in the central regions of the most massive clusters, occurring already earlier than $z\sim1.5$. [Abridged] △ Less

Submitted 6 February, 2019; v1 submitted 25 July, 2018; originally announced July 2018.

Comments: 29 pages, 16 figures, A&A in press. Replaced to match published version, results unchanged

Journal ref: A&A 622, A117 (2019)