Showing 51–100 of 596 results for author: Lewis, G

The Instantaneous Redshift Difference of Gravitationally Lensed Images: Theory and Observational Prospects

Authors: Chengyi Wang, Krzysztof Bolejko, Geraint F. Lewis

Abstract: Due to the expansion of our Universe, the redshift of distant objects changes with time. Although the amplitude of this redshift drift is small, it will be measurable with a decade-long campaigns on the next generation of telescopes. Here we present an alternative view of the redshift drift which captures the expansion of the universe in single epoch observations of the multiple images of gravitat… ▽ More Due to the expansion of our Universe, the redshift of distant objects changes with time. Although the amplitude of this redshift drift is small, it will be measurable with a decade-long campaigns on the next generation of telescopes. Here we present an alternative view of the redshift drift which captures the expansion of the universe in single epoch observations of the multiple images of gravitationally lensed sources. Considering a sufficiently massive lens, with an associated time delay of order decades, simultaneous photons arriving at a detector would have been emitted decades earlier in one image compared to another, leading to an instantaneous redshift difference between the images. We also investigate the effect of peculiar velocities on the redshift difference in the observed images. Whilst still requiring the observational power of the next generation of telescopes and instruments, the advantage of such a single epoch detection over other redshift drift measurements is that it will be less susceptible to systematic effects that result from requiring instrument stability over decade-long campaigns. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 6 pages, 5 figures

Predicting Swarm Equatorial Plasma Bubbles via Machine Learning and Shapley Values

Authors: S. A. Reddy, C. Forsyth, A. Aruliah, A. Smith, J. Bortnik, E. Aa, D. O. Kataria, G. Lewis

Abstract: In this study we present AI Prediction of Equatorial Plasma Bubbles (APE), a machine learning model that can accurately predict the Ionospheric Bubble Index (IBI) on the Swarm spacecraft. IBI is a correlation ($R^2$) between perturbations in plasma density and the magnetic field, whose source can be Equatorial Plasma Bubbles (EPBs). EPBs have been studied for a number of years, but their day-to-da… ▽ More In this study we present AI Prediction of Equatorial Plasma Bubbles (APE), a machine learning model that can accurately predict the Ionospheric Bubble Index (IBI) on the Swarm spacecraft. IBI is a correlation ($R^2$) between perturbations in plasma density and the magnetic field, whose source can be Equatorial Plasma Bubbles (EPBs). EPBs have been studied for a number of years, but their day-to-day variability has made predicting them a considerable challenge. We build an ensemble machine learning model to predict IBI. We use data from 2014-22 at a resolution of 1sec, and transform it from a time-series into a 6-dimensional space with a corresponding EPB $R^2$ (0-1) acting as the label. APE performs well across all metrics, exhibiting a skill, association and root mean squared error score of 0.96, 0.98 and 0.08 respectively. The model performs best post-sunset, in the American/Atlantic sector, around the equinoxes, and when solar activity is high. This is promising because EPBs are most likely to occur during these periods. Shapley values reveal that F10.7 is the most important feature in driving the predictions, whereas latitude is the least. The analysis also examines the relationship between the features, which reveals new insights into EPB climatology. Finally, the selection of the features means that APE could be expanded to forecasting EPBs following additional investigations into their onset. △ Less

Submitted 30 September, 2023; v1 submitted 27 September, 2022; originally announced September 2022.

Comments: 13 Pages, 9 Figures

Journal ref: Journal of Geophysical Research: Space Physics (2023): e2022JA031183

Multi-level Adversarial Spatio-temporal Learning for Footstep Pressure based FoG Detection

Authors: Kun Hu, Shaohui Mei, Wei Wang, Kaylena A. Ehgoetz Martens, Liang Wang, Simon J. G. Lewis, David D. Feng, Zhiyong Wang

Abstract: Freezing of gait (FoG) is one of the most common symptoms of Parkinson's disease, which is a neurodegenerative disorder of the central nervous system impacting millions of people around the world. To address the pressing need to improve the quality of treatment for FoG, devising a computer-aided detection and quantification tool for FoG has been increasingly important. As a non-invasive technique… ▽ More Freezing of gait (FoG) is one of the most common symptoms of Parkinson's disease, which is a neurodegenerative disorder of the central nervous system impacting millions of people around the world. To address the pressing need to improve the quality of treatment for FoG, devising a computer-aided detection and quantification tool for FoG has been increasingly important. As a non-invasive technique for collecting motion patterns, the footstep pressure sequences obtained from pressure sensitive gait mats provide a great opportunity for evaluating FoG in the clinic and potentially in the home environment. In this study, FoG detection is formulated as a sequential modelling task and a novel deep learning architecture, namely Adversarial Spatio-temporal Network (ASTN), is proposed to learn FoG patterns across multiple levels. A novel adversarial training scheme is introduced with a multi-level subject discriminator to obtain subject-independent FoG representations, which helps to reduce the over-fitting risk due to the high inter-subject variance. As a result, robust FoG detection can be achieved for unseen subjects. The proposed scheme also sheds light on improving subject-level clinical studies from other scenarios as it can be integrated with many existing deep architectures. To the best of our knowledge, this is one of the first studies of footstep pressure-based FoG detection and the approach of utilizing ASTN is the first deep neural network architecture in pursuit of subject-independent representations. Experimental results on 393 trials collected from 21 subjects demonstrate encouraging performance of the proposed ASTN for FoG detection with an AUC 0.85. △ Less

Submitted 22 September, 2022; originally announced September 2022.

Data Leakage in Notebooks: Static Detection and Better Processes

Authors: Chenyang Yang, Rachel A Brower-Sinning, Grace A. Lewis, Christian Kästner

Abstract: Data science pipelines to train and evaluate models with machine learning may contain bugs just like any other code. Leakage between training and test data can lead to overestimating the model's accuracy during offline evaluations, possibly leading to deployment of low-quality models in production. Such leakage can happen easily by mistake or by following poor practices, but may be tedious and cha… ▽ More Data science pipelines to train and evaluate models with machine learning may contain bugs just like any other code. Leakage between training and test data can lead to overestimating the model's accuracy during offline evaluations, possibly leading to deployment of low-quality models in production. Such leakage can happen easily by mistake or by following poor practices, but may be tedious and challenging to detect manually. We develop a static analysis approach to detect common forms of data leakage in data science code. Our evaluation shows that our analysis accurately detects data leakage and that such leakage is pervasive among over 100,000 analyzed public notebooks. We discuss how our static analysis approach can help both practitioners and educators, and how leakage prevention can be designed into the development process. △ Less

Submitted 7 September, 2022; originally announced September 2022.

OzDES Reverberation Mapping Program: Mg II Lags and R-L relation

Authors: Zhefu Yu, Paul Martini, A. Penton, T. M. Davis, C. S. Kochanek, G. F. Lewis, C. Lidman, U. Malik, R. Sharp, B. E. Tucker, M. Aguena, J. Annis, E. Bertin, S. Bocquet, D. Brooks, A. Carnero Rosell, D. Carollo, M. Carrasco Kind, J. Carretero, M. Costanzi, L. N. da Costa, M. E. S. Pereira, J. De Vicente, H. T. Diehl, P. Doel , et al. (33 additional authors not shown)

Abstract: The correlation between the broad line region radius and continuum luminosity ($R-L$ relation) of active galactic nuclei (AGN) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At $z \sim 1-2$, where AGN activity peaks, the $R-L$ relation is constrained by the reverberation mapping (RM) lags of the Mg II line. We present 25 Mg II lags from the Australian Dark Energy… ▽ More The correlation between the broad line region radius and continuum luminosity ($R-L$ relation) of active galactic nuclei (AGN) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At $z \sim 1-2$, where AGN activity peaks, the $R-L$ relation is constrained by the reverberation mapping (RM) lags of the Mg II line. We present 25 Mg II lags from the Australian Dark Energy Survey (OzDES) RM project based on six years of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the re-scaled, re-sampled versions of the observed lightcurves of local, well-measured AGN. Our sample significantly increases the number of Mg II lags and extends the $R-L$ relation to higher redshifts and luminosities. The relative iron line strength $\mathcal{R}_{\rm Fe}$ has little impact on the $R-L$ relation. The best-fit Mg II $R-L$ relation has a slope $α= 0.39 \pm 0.08$ with an intrinsic scatter $σ_{\rm rl} = 0.15^{+0.03}_{-0.02}$. The slope is consistent with previous measurements and shallower than the H$β$ $R-L$ relation. The intrinsic scatter of the new $R-L$ relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H$β$ $R-L$ relation. Our new $R-L$ relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon. △ Less

Submitted 1 August, 2023; v1 submitted 10 August, 2022; originally announced August 2022.

Comments: 16 pages, 8 figures; MNRAS, Volume 522, pp.4132

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

Streams on FIRE: Populations of Detectable Stellar Streams in the Milky Way and FIRE

Authors: Nora Shipp, Nondh Panithanpaisal, Lina Necib, Robyn Sanderson, Denis Erkal, Ting S. Li, Isaiah B. Santistevan, Andrew Wetzel, Lara R. Cullinane, Alexander P. Ji, Sergey E. Koposov, Kyler Kuehn, Geraint F. Lewis, Andrew B. Pace, Daniel B. Zucker, Joss Bland-Hawthorn, Emily C. Cunningham, Stacy Y. Kim, Sophia Lilleengen, Jorge Moreno, Sanjib Sharma

Abstract: We present the first detailed study comparing the populations of stellar streams in cosmological simulations to observed Milky Way dwarf galaxy streams. In particular, we compare streams identified around Milky Way analogs in the FIRE-2 simulations to stellar streams observed by the Southern Stellar Stream Spectroscopic Survey (S5). For an accurate comparison between the stream populations, we pro… ▽ More We present the first detailed study comparing the populations of stellar streams in cosmological simulations to observed Milky Way dwarf galaxy streams. In particular, we compare streams identified around Milky Way analogs in the FIRE-2 simulations to stellar streams observed by the Southern Stellar Stream Spectroscopic Survey (S5). For an accurate comparison between the stream populations, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to > 110 kpc) and surviving only at larger apocenters (> 40 kpc) than those observed in the Milky Way. We find that the population of high-stellar mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as satellites in DES-like observations, since their tidal tails are too low-surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low surface brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: 24 pages, 13 figures, 3 tables, submitted to ApJ

Concerning Colour: The Effect of Environment on Type Ia Supernova Colour in the Dark Energy Survey

Authors: L. Kelsey, M. Sullivan, P. Wiseman, P. Armstrong, R. Chen, D. Brout, T. M. Davis, M. Dixon, C. Frohmaier, L. Galbany, O. Graur, R. Kessler, C. Lidman, A. Möller, B. Popovic, B. Rose, D. Scolnic, M. Smith, M. Vincenzi, T. M. C. Abbott, M. Aguena, S. Allam, O. Alves, J. Annis, D. Bacon , et al. (45 additional authors not shown)

Abstract: Recent analyses have found intriguing correlations between the colour ($c$) of type Ia supernovae (SNe Ia) and the size of their 'mass-step', the relationship between SN Ia host galaxy stellar mass ($M_\mathrm{stellar}$) and SN Ia Hubble residual, and suggest that the cause of this relationship is dust. Using 675 photometrically-classified SNe Ia from the Dark Energy Survey 5-year sample, we study… ▽ More Recent analyses have found intriguing correlations between the colour ($c$) of type Ia supernovae (SNe Ia) and the size of their 'mass-step', the relationship between SN Ia host galaxy stellar mass ($M_\mathrm{stellar}$) and SN Ia Hubble residual, and suggest that the cause of this relationship is dust. Using 675 photometrically-classified SNe Ia from the Dark Energy Survey 5-year sample, we study the differences in Hubble residual for a variety of global host galaxy and local environmental properties for SN Ia subsamples split by their colour. We find a $3σ$ difference in the mass-step when comparing blue ($c<0$) and red ($c>0$) SNe. We observe the lowest r.m.s. scatter ($\sim0.14$ mag) in the Hubble residual for blue SNe in low mass/blue environments, suggesting that this is the most homogeneous sample for cosmological analyses. By fitting for $c$-dependent relationships between Hubble residuals and $M_\mathrm{stellar}$, approximating existing dust models, we remove the mass-step from the data and find tentative $\sim 2σ$ residual steps in rest-frame galaxy $U-R$ colour. This indicates that dust modelling based on $M_\mathrm{stellar}$ may not fully explain the remaining dispersion in SN Ia luminosity. Instead, accounting for a $c$-dependent relationship between Hubble residuals and global $U-R$, results in $\leq1σ$ residual steps in $M_\mathrm{stellar}$ and local $U-R$, suggesting that $U-R$ provides different information about the environment of SNe Ia compared to $M_\mathrm{stellar}$, and motivating the inclusion of galaxy $U-R$ colour in SN Ia distance bias correction. △ Less

Submitted 28 February, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 19 pages, 8 figures. Published in MNRAS

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

Journal ref: MNRAS, February 2023, Volume 519, Issue 2, Pages 3046-3063

Joint HST, VLT/MUSE and XMM-Newton observations to constrain the mass distribution of the two strong lensing galaxy clusters: MACS J0242.5-2132 & MACS J0949.8+1708

Authors: Joseph F. V. Allingham, Mathilde Jauzac, David J. Lagattuta, Guillaume Mahler, Céline Bœhm, Geraint F. Lewis, Dominique Eckert, Alastair Edge, Stefano Ettori

Abstract: We present the strong lensing analysis of two galaxy clusters: MACS J0242.5-2132 (MACS J0242, $z=0.313$) and MACS J0949.8+1708 (MACS J0949, $z=0.383$). Their total matter distributions are constrained thanks to the powerful combination of observations with the Hubble Space Telescope and the MUSE instrument. Using these observations, we precisely measure the redshift of six multiple image systems i… ▽ More We present the strong lensing analysis of two galaxy clusters: MACS J0242.5-2132 (MACS J0242, $z=0.313$) and MACS J0949.8+1708 (MACS J0949, $z=0.383$). Their total matter distributions are constrained thanks to the powerful combination of observations with the Hubble Space Telescope and the MUSE instrument. Using these observations, we precisely measure the redshift of six multiple image systems in MACS J0242, and two in MACS J0949. We also include four multiple image systems in the latter cluster identified in HST imaging without MUSE redshift measurements. For each cluster, our best-fit mass model consists of a single cluster-scale halo, and 57 (170) galaxy-scale halos for MACS J0242 (MACS J0949). Multiple images positions are predicted with a $rms$ 0.39 arcsec and 0.15 arcsec for MACS J0242 and MACS J0949 models respectively. From these mass models, we derive aperture masses of $M(R<$200 kpc$) = 1.67_{-0.05}^{+0.03}\times 10^{14} M_{\odot}$, and $M(R<$200 kpc$) = 2.00_{-0.20}^{+0.05}\times 10^{14} M_{\odot}$. Combining our analysis with X-ray observations from the XMM-Newton Observatory, we show that MACS J0242 appears to be a relatively relaxed cluster, while conversely, MACS J0949 shows a relaxing post-merger state. At 200 kpc, X-ray observations suggest the hot gas fraction to be respectively $f_g = 0.115^{+0.003}_{-0.004}$ and $0.053^{+0.007}_{-0.006}$ for MACS J0242 and MACS J0949. MACS J0242 being relaxed, its density profile is very well fitted by a NFW distribution, in agreement with X-ray observations. Finally, the strong lensing analysis of MACS J0949 suggests a flat dark matter density distribution in the core, between 10 and 100 kpc. This appears consistent with X-ray observations. △ Less

Submitted 19 April, 2023; v1 submitted 21 July, 2022; originally announced July 2022.

Comments: 20 pages, 11 figures, published in MNRAS

Journal ref: MNRAS, 522, 1118-1137 (2023)

Core-collapse Supernovae in the Dark Energy Survey: Luminosity Functions and Host Galaxy Demographics

Authors: M. Grayling, C. P. Gutiérrez, M. Sullivan, P. Wiseman, M. Vincenzi, L. Galbany, A. Möller, D. Brout, T. M. Davis, C. Frohmaier, O. Graur, L. Kelsey, C. Lidman, B. Popovic, M. Smith, M. Toy, B. E. Tucker, Z. Zontou, T. M. C. Abbott, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Annis, J. Asorey, D. Bacon , et al. (51 additional authors not shown)

Abstract: We present the luminosity functions and host galaxy properties of the Dark Energy Survey (DES) core-collapse supernova (CCSN) sample, consisting of 69 Type II and 50 Type Ibc spectroscopically and photometrically-confirmed supernovae over a redshift range $0.045<z<0.25$. We fit the observed DES $griz$ CCSN light-curves and K-correct to produce rest-frame $R$-band light curves. We compare the sampl… ▽ More We present the luminosity functions and host galaxy properties of the Dark Energy Survey (DES) core-collapse supernova (CCSN) sample, consisting of 69 Type II and 50 Type Ibc spectroscopically and photometrically-confirmed supernovae over a redshift range $0.045<z<0.25$. We fit the observed DES $griz$ CCSN light-curves and K-correct to produce rest-frame $R$-band light curves. We compare the sample with lower-redshift CCSN samples from Zwicky Transient Facility (ZTF) and Lick Observatory Supernova Search (LOSS). Comparing luminosity functions, the DES and ZTF samples of SNe II are brighter than that of LOSS with significances of 3.0$σ$ and 2.5$σ$ respectively. While this difference could be caused by redshift evolution in the luminosity function, simpler explanations such as differing levels of host extinction remain a possibility. We find that the host galaxies of SNe II in DES are on average bluer than in ZTF, despite having consistent stellar mass distributions. We consider a number of possibilities to explain this -- including galaxy evolution with redshift, selection biases in either the DES or ZTF samples, and systematic differences due to the different photometric bands available -- but find that none can easily reconcile the differences in host colour between the two samples and thus its cause remains uncertain. △ Less

Submitted 22 March, 2023; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: 28 pages, 13 figures, 5 tables. Accepted by MNRAS Dec 2022

Journal ref: MNRAS 520 (2023) 684-701

Using Host Galaxy Spectroscopy to Explore Systematics in the Standardisation of Type Ia Supernovae

Authors: M. Dixon, C. Lidman, J. Mould, L. Kelsey, D. Brout, A. Möller, P. Wiseman, M. Sullivan, L. Galbany, T. M. Davis, M. Vincenzi, D. Scolnic, G. F. Lewis, M. Smith, R. Kessler, A. Duffy, E. Taylor, C. Flynn, T. M. C. Abbott, M. Aguena, S. Allam, F. Andrade-Oliveir, J. Annis, J. Asorey, E. Bertin , et al. (53 additional authors not shown)

Abstract: We use stacked spectra of the host galaxies of photometrically identified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilising full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and prope… ▽ More We use stacked spectra of the host galaxies of photometrically identified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilising full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and properties of the host galaxies that rely on spectral absorption features ($< 1.3σ$), such as stellar population age, metallicity, and mass-to-light ratio. However, we find significant trends between the Hubble residuals and the strengths of [OII] ($4.4σ$) and the Balmer emission lines ($3σ$). These trends are weaker than the well known trend between Hubble residuals and host galaxy stellar mass ($7.2σ$) that is derived from broad band photometry. After light curve corrections, we see fainter SNe Ia residing in galaxies with larger line strengths. We also find a trend (3$σ$) between Hubble residual and the Balmer decrement (a measure of reddening by dust) using H$β$ and H$γ$. The trend, quantified by correlation coefficients, is slightly more significant in the redder SNe Ia, suggesting that bluer SNe Ia are relatively unaffected by dust in the interstellar medium of the host and that dust contributes to current Hubble diagram scatter impacting the measurement of cosmological parameters. △ Less

Submitted 24 October, 2022; v1 submitted 24 June, 2022; originally announced June 2022.

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

Report number: DES-2018-0379 FERMILAB-PUB-22-428-PPD

The Dark Energy Survey Supernova Program results: Type Ia Supernova brightness correlates with host galaxy dust

Authors: Cole Meldorf, Antonella Palmese, Dillon Brout, Rebecca Chen, Daniel Scolnic, Lisa Kelsey, Lluís Galbany, Will Hartley, Tamara Davis, Alex Drlica-Wagner, Maria Vincenzi, James Annis, Mitchell Dixon, Or Graur, Alex Kim, Christopher Lidman, Anais Möller, Peter Nugent, Benjamin Rose, Mathew Smith, Sahar Allam, H. Thomas Diehl, Douglas Tucker, Jacobo Asorey, Josh Calcino , et al. (46 additional authors not shown)

Abstract: Cosmological analyses with type Ia supernovae (SNe Ia) often assume a single empirical relation between color and luminosity ($β$) and do not account for varying host-galaxy dust properties. However, from studies of dust in large samples of galaxies, it is known that dust attenuation can vary significantly. Here we take advantage of state-of-the-art modeling of galaxy properties to characterize du… ▽ More Cosmological analyses with type Ia supernovae (SNe Ia) often assume a single empirical relation between color and luminosity ($β$) and do not account for varying host-galaxy dust properties. However, from studies of dust in large samples of galaxies, it is known that dust attenuation can vary significantly. Here we take advantage of state-of-the-art modeling of galaxy properties to characterize dust parameters (dust attenuation $A_V$, and a parameter describing the dust law slope $R_V$) for the Dark Energy Survey (DES) SN Ia host galaxies using the publicly available \texttt{BAGPIPES} code. Utilizing optical and infrared data of the hosts alone, we find three key aspects of host dust that impact SN Ia cosmology: 1) there exists a large range ($\sim1-6$) of host $R_V$ 2) high stellar mass hosts have $R_V$ on average $\sim0.7$ lower than that of low-mass hosts 3) there is a significant ($>3σ$) correlation between the Hubble diagram residuals of red SNe Ia that when corrected for reduces scatter by $\sim13\%$ and the significance of the ``mass step'' to $\sim1σ$. These represent independent confirmations of recent predictions based on dust that attempted to explain the puzzling ``mass step'' and intrinsic scatter ($σ_{\rm int}$) in SN Ia analyses. We also find that red-sequence galaxies have both lower and more peaked dust law slope distributions on average in comparison to non red-sequence galaxies. We find that the SN Ia $β$ and $σ_{\rm int}$ both differ by $>3σ$ when determined separately for red-sequence galaxy and all other galaxy hosts. The agreement between fitted host-$R_V$ and SN Ia $β$ \& $σ_{\rm int}$ suggests that host dust properties play a major role in SN Ia color-luminosity standardization and supports the claim that SN Ia intrinsic scatter is driven by $R_V$ variation. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 22 pages. Submitted to MNRAS

Report number: DES-2021-0641 FERMILAB-PUB-21-051-AE DES-2021-0641 FERMILAB-PUB-21-051-AE FERMILAB-PUB-21-051-AE

The Hubble Space Telescope Survey of M31 Satellite Galaxies I. RR Lyrae-based Distances and Refined 3D Geometric Structure

Authors: Alessandro Savino, Daniel R. Weisz, Evan D. Skillman, Andrew Dolphin, Nitya Kallivayalil, Andrew Wetzel, Jay Anderson, Gurtina Besla, Michael Boylan-Kolchin, James S. Bullock, Andrew A. Cole, Michelle L. M. Collins, M. C. Cooper, Alis J. Deason, Aaron L. Dotter, Mark Fardal, Annette M. N. Ferguson, Tobias K. Fritz, Marla C. Geha, Karoline M. Gilbert, Puragra Guhathakurta, Rodrigo Ibata, Michael J. Irwin, Myoungwon Jeon, Evan Kirby , et al. (11 additional authors not shown)

Abstract: We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days… ▽ More We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system. △ Less

Submitted 12 September, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: 32 pages, 15 figures, accepted for publication on ApJ

The PAndAS View of the Andromeda Satellite System. III. Dwarf galaxy detection limits

Authors: Amandine Doliva-Dolinsky, Nicolas F. Martin, Guillaume F. Thomas, Annette M. N. Ferguson, Rodrigo A. Ibata, Geraint F. Lewis, Dougal Mackey, Alan W. McConnachie, Zhen Yuan

Abstract: We determine the detection limits of the search for dwarf galaxies in the Pan-Andromeda Archaeological Survey (PAndAS) using the algorithm developed by the PAndAS team. The recovery fractions of artificial dwarf galaxies are, as expected, a strong function of physical size and luminosity and, to a lesser extent, distance. We show that these recovery fractions vary strongly with location in the sur… ▽ More We determine the detection limits of the search for dwarf galaxies in the Pan-Andromeda Archaeological Survey (PAndAS) using the algorithm developed by the PAndAS team. The recovery fractions of artificial dwarf galaxies are, as expected, a strong function of physical size and luminosity and, to a lesser extent, distance. We show that these recovery fractions vary strongly with location in the surveyed area because of varying levels of contamination from both the Milky Way foreground stars and the stellar halo of Andromeda. We therefore provide recovery fractions that are a function of size, luminosity, and location within the survey on a scale of 1 square degree. Overall, the effective surface brightness for a 50-percent detection rate range between 28 and 30 mag per square arcsecond. This is in line with expectations for a search that relies on photometric data that are as deep as the PAndAS survey. The derived detection limits are an essential ingredient on the path to constraining the global properties of Andromeda's system of satellite dwarf galaxies and, more broadly, to provide constraints on dwarf galaxy formation and evolution in a cosmological context. △ Less

Submitted 5 May, 2022; originally announced May 2022.

Comments: 17 pages, 10 figures, accepted in ApJ

The effect of the deforming dark matter haloes of the Milky Way and the Large Magellanic Cloud on the Orphan-Chenab stream

Authors: Sophia Lilleengen, Michael S. Petersen, Denis Erkal, Jorge Peñarrubia, Sergey E. Koposov, Ting S. Li, Lara R. Cullinane, Alexander P. Ji, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Andrew B. Pace, Nora Shipp, Daniel B. Zucker, Joss Bland-Hawthorn, Tariq Hilmi

Abstract: It has recently been shown that the Large Magellanic Cloud (LMC) has a substantial effect on the Milky Way's stellar halo and stellar streams. Here, we explore how deformations of the Milky Way and LMC's dark matter haloes affect stellar streams, and whether these effects are observable. In particular, we focus on the Orphan-Chenab (OC) stream which passes particularly close to the LMC, and spans… ▽ More It has recently been shown that the Large Magellanic Cloud (LMC) has a substantial effect on the Milky Way's stellar halo and stellar streams. Here, we explore how deformations of the Milky Way and LMC's dark matter haloes affect stellar streams, and whether these effects are observable. In particular, we focus on the Orphan-Chenab (OC) stream which passes particularly close to the LMC, and spans a large portion of the Milky Way's halo. We represent the Milky Way--LMC system using basis function expansions that capture their evolution in an $N$-body simulation. We present the properties of this system, such as the evolution of the densities and force fields of each galaxy. The OC stream is evolved in this time-dependent, deforming potential, and we investigate the effects of the various moments of the Milky Way and the LMC. We find that the simulated OC stream is strongly influenced by the deformations of both the Milky Way and the LMC, and that this effect is much larger than current observational errors. In particular, the Milky Way dipole has the biggest impact on the stream, followed by the evolution of the LMC's monopole, and the LMC's quadrupole. Detecting these effects would confirm a key prediction of collisionless, cold dark matter, and would be a powerful test of alternative dark matter and alternative gravity models. △ Less

Submitted 10 November, 2022; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: 17 pages, 10 figures. Accepted for publication in MNRAS

Ringing the universe with cosmic emptiness: void properties through a combined analysis of stacked weak gravitational and Doppler lensing

Authors: Md Rasel Hossen, Sonia Akter Ema, Krzysztof Bolejko, Geraint F. Lewis

Abstract: An essential aspect of cosmic voids is that these underdense regions provide complementary information about the properties of our Universe. Unlike dense regions, voids are avoided by matter and are less contaminated by baryonic processes. The first step to understanding the properties of cosmic voids is to correctly infer their mass profiles. In the literature, various techniques have been implem… ▽ More An essential aspect of cosmic voids is that these underdense regions provide complementary information about the properties of our Universe. Unlike dense regions, voids are avoided by matter and are less contaminated by baryonic processes. The first step to understanding the properties of cosmic voids is to correctly infer their mass profiles. In the literature, various techniques have been implemented. In this paper, we review them and implement a new technique that is based on Doppler lensing. We use a relativistic $N$-body code, \textsc{Gevolution}, to generate cosmological mass perturbations and implement a three-dimensional ray-tracing technique, which follows the evolution of a ray-bundles. We focus on the various properties of cosmic voids (e.g. void size function, 2-point correlation function, and the density profile of voids), and compare the results with their universal trends. We show that when weak-lensing is combined with the Doppler lensing we obtain even tighter constraints than weak-lensing alone. We also obtain better agreement between density profiles within central parts of voids inferred from lensing and density profiles inferred from halo tracers. The implication of the result relevant to the ongoing and prospective low-redshift spectroscopic surveys is briefly discussed. △ Less

Submitted 1 May, 2022; originally announced May 2022.

Comments: 14 pages, 10 figures, Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 513, 5575, 2022

The Pristine Inner Galaxy Survey (PIGS) IV: A photometric metallicity analysis of the Sagittarius dwarf spheroidal galaxy

Authors: Sara Vitali, Anke Arentsen, Else Starkenburg, Paula Jofré, Nicolas F. Martin, David S. Aguado, Raymond Carlberg, Jonay I. González Hernández, Rodrigo Ibata, Georges Kordopatis, Khyati Malhan, Pau Ramos, Federico Sestito, Zhen Yuan, Sven Buder, Geraint F. Lewis, Zhen Wan, Daniel B. Zucker

Abstract: We present a comprehensive metallicity analysis of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using $Pristine\,CaHK$ photometry. We base our member selection on $Gaia$ EDR3 astrometry applying a magnitude limit at $G_{0} = 17.3$, and our population study on the metallicity-sensitive photometry from the $Pristine$ Inner Galaxy Survey (PIGS). Working with photometric metallicities instead of… ▽ More We present a comprehensive metallicity analysis of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using $Pristine\,CaHK$ photometry. We base our member selection on $Gaia$ EDR3 astrometry applying a magnitude limit at $G_{0} = 17.3$, and our population study on the metallicity-sensitive photometry from the $Pristine$ Inner Galaxy Survey (PIGS). Working with photometric metallicities instead of spectroscopic metallicities allows us to cover an unprecedented large area ($\sim 100$ square degrees) of the dwarf galaxy, and to study the spatial distribution of its members as function of metallicity with little selection effects. Our study compares the spatial distributions of a metal-poor population of 9719 stars with [Fe/H] $< -1.3$ and a metal rich one of 30115 stars with [Fe/H] $> -1.0$. The photometric Sgr sample also allows us to assemble the largest sample of 1150 very metal-poor Sgr candidates ([Fe/H] $< -2.0$). By investigating and fitting the spatial properties of the metal-rich and metal-poor population, we find a negative metallicity gradient which extends up to 12 degrees from the Sgr center (or $\sim 5.5$ kpc at the distance of Sgr), the limit of our footprint. We conclude that the relative number of metal-poor stars increases in the outer areas of the galaxy, while the central region is dominated by metal-rich stars. These finding suggest an outside-in formation process and are an indication of the extended formation history of Sgr, which has been affected by the tidal interaction between Sgr and the Milky Way. △ Less

Submitted 4 October, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: The paper was accepted on 20 September by MNRAS

Density Matrix Renormalization Group with Tensor Processing Units

Authors: Martin Ganahl, Jackson Beall, Markus Hauru, Adam G. M. Lewis, Jae Hyeon Yoo, Yijian Zou, Guifre Vidal

Abstract: Google's Tensor Processing Units (TPUs) are integrated circuits specifically built to accelerate and scale up machine learning workloads. They can perform fast distributed matrix multiplications and therefore be repurposed for other computationally intensive tasks. In this work we demonstrate the use of TPUs for accelerating and scaling up the density matrix renormalization group (DMRG), a powerfu… ▽ More Google's Tensor Processing Units (TPUs) are integrated circuits specifically built to accelerate and scale up machine learning workloads. They can perform fast distributed matrix multiplications and therefore be repurposed for other computationally intensive tasks. In this work we demonstrate the use of TPUs for accelerating and scaling up the density matrix renormalization group (DMRG), a powerful numerical approach to compute the ground state of a local quantum many-body Hamiltonian. The cost of DMRG scales with system size $N$ as $O(ND^3)$, where the so-called bond dimension $D$ regulates how expressive the underlying matrix product state (MPS) variational ansatz is. We consider lattice models in two spatial dimensions, with square lattices of size $10\times 10$ (free fermions) and $20\times 20$ (transverse field Ising model), for which the required MPS bond dimension is known to scale at least as $\exp(\sqrt{N})$. Using half of a TPU v3 pod (namely $1,\!024$ TPU v3 cores) we reached an unprecedentedly large bond dimension $D = 2^{16} = 65,\!536$, for which optimizing a single MPS tensor took about 2 minutes. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Dynamics of Transmon Ionization

Authors: Ross Shillito, Alexandru Petrescu, Joachim Cohen, Jackson Beall, Markus Hauru, Martin Ganahl, Adam G. M. Lewis, Guifre Vidal, Alexandre Blais

Abstract: Qubit measurement and control in circuit QED rely on microwave drives, with higher drive amplitudes ideally leading to faster processes. However, degradation in qubit coherence time and readout fidelity has been observed even under moderate drive amplitudes corresponding to few photons populating the measurement resonator. Here, we numerically explore the dynamics of a driven transmon-resonator sy… ▽ More Qubit measurement and control in circuit QED rely on microwave drives, with higher drive amplitudes ideally leading to faster processes. However, degradation in qubit coherence time and readout fidelity has been observed even under moderate drive amplitudes corresponding to few photons populating the measurement resonator. Here, we numerically explore the dynamics of a driven transmon-resonator system under strong and nearly resonant measurement drives, and find clear signatures of transmon ionization where the qubit escapes out of its cosine potential. Using a semiclassical model, we interpret this ionization as resulting from resonances occurring at specific resonator photon populations. We find that the photon populations at which these spurious transitions occur are strongly parameter dependent and that they can occur at low resonator photon population, something which may explain the experimentally observed degradation in measurement fidelity. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 10 pages, 8 figures

The GALAH Survey: A New Sample of Extremely Metal-Poor Stars Using A Machine Learning Classification Algorithm

Authors: Arvind C. N. Hughes, Lee R. Spitler, Daniel B. Zucker, Thomas Nordlander, Jeffrey Simpson, Gary S. Da Costa, Yuan-Sen Ting, Chengyuan Li, Joss Bland-Hawthorn, Sven Buder, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Michael R. Hayden, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Katharine J. Schlesinger, Sanjib Sharma, Tomaz Zwitter, The GALAH Collaboration

Abstract: Extremely Metal-Poor (EMP) stars provide a valuable probe of early chemical enrichment in the Milky Way. Here we leverage a large sample of $\sim600,000$ high-resolution stellar spectra from the GALAH survey plus a machine learning algorithm to find 54 candidates with estimated [Fe/H]~$\leq$~-3.0, 6 of which have [Fe/H]~$\leq$~-3.5. Our sample includes $\sim 20 \%$ main sequence EMP candidates, un… ▽ More Extremely Metal-Poor (EMP) stars provide a valuable probe of early chemical enrichment in the Milky Way. Here we leverage a large sample of $\sim600,000$ high-resolution stellar spectra from the GALAH survey plus a machine learning algorithm to find 54 candidates with estimated [Fe/H]~$\leq$~-3.0, 6 of which have [Fe/H]~$\leq$~-3.5. Our sample includes $\sim 20 \%$ main sequence EMP candidates, unusually high for \emp surveys. We find the magnitude-limited metallicity distribution function of our sample is consistent with previous work that used more complex selection criteria. The method we present has significant potential for application to the next generation of massive stellar spectroscopic surveys, which will expand the available spectroscopic data well into the millions of stars. △ Less

Submitted 8 August, 2022; v1 submitted 21 March, 2022; originally announced March 2022.

Comments: 27 pages, 20 figures, accepted for publication in ApJ, candidate table available at this https://github.com/arvhug/GALAH---TSNE_EMP

LensingETC: a tool to optimize multi-filter imaging campaigns of galaxy-scale strong lensing systems

Authors: Anowar J. Shajib, Karl Glazebrook, Tania Barone, Geraint F. Lewis, Tucker Jones, Kim-Vy H. Tran, Elizabeth Buckley-Geer, Thomas E. Collett, Joshua Frieman, Colin Jacobs

Abstract: Imaging data is the principal observable required to use galaxy-scale strong lensing in a multitude of applications in extragalactic astrophysics and cosmology. In this paper, we develop Lensing Exposure Time Calculator (LensingETC) to optimize the efficiency of telescope time usage when planning multi-filter imaging campaigns for galaxy-scale strong lenses. This tool simulates realistic data tail… ▽ More Imaging data is the principal observable required to use galaxy-scale strong lensing in a multitude of applications in extragalactic astrophysics and cosmology. In this paper, we develop Lensing Exposure Time Calculator (LensingETC) to optimize the efficiency of telescope time usage when planning multi-filter imaging campaigns for galaxy-scale strong lenses. This tool simulates realistic data tailored to specified instrument characteristics and then automatically models them to assess the power of the data in constraining lens model parameters. We demonstrate a use case of this tool by optimizing a two-filter observing strategy (in IR and UVIS) within the limited exposure time per system allowed by a Hubble Space Telescope (HST) Snapshot program. We find that higher resolution is more advantageous to gain constraining power on the lensing observables, when there is a trade-off between signal-to-noise ratio and resolution; e.g., between the UVIS and IR filters of the HST. We also find that, whereas a point spread function (PSF) with sub-Nyquist sampling allows the sample mean for a model parameter to be robustly recovered for both galaxy-galaxy and point-source lensing systems, a sub-Nyquist sampled PSF introduces a larger scatter than a Nyquist sampled one in the deviation from the ground truth for point-source lens systems. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: 10 pages, 7 figures, 3 tables

Observational window effects on multi-object Reverberation Mapping

Authors: Umang Malik, Rob Sharp, Paul Martini, Tamara M. Davis, Brad E. Tucker, Zhefu Yu, Andrew Penton, Geraint F. Lewis, Josh Calcino

Abstract: Contemporary reverberation mapping campaigns are employing wide-area photometric data and high-multiplex spectroscopy to efficiently monitor hundreds of active galactic nuclei (AGN). However, the interaction of the window function(s) imposed by the observation cadence with the reverberation lag and AGN variability time scales (intrinsic to each source over a range of luminosities) impact our abili… ▽ More Contemporary reverberation mapping campaigns are employing wide-area photometric data and high-multiplex spectroscopy to efficiently monitor hundreds of active galactic nuclei (AGN). However, the interaction of the window function(s) imposed by the observation cadence with the reverberation lag and AGN variability time scales (intrinsic to each source over a range of luminosities) impact our ability to recover these fundamental physical properties. Time dilation effects due to the sample source redshift distribution introduces added complexity. We present comprehensive analysis of the implications of observational cadence, seasonal gaps and campaign baseline duration (i.e., the survey window function) for reverberation lag recovery. We find the presence of a significant seasonal gap dominates the efficacy of any given campaign strategy for lag recovery across the parameter space, particularly for those sources with observed-frame lags above 100 days. Using the OzDES survey as a baseline, we consider the implications of this analysis for the 4MOST/TiDES campaign providing concurrent follow-up of the LSST deep-drilling fields, as well as upcoming programs. We conclude that the success of such surveys will be critically limited by the seasonal visibility of some potential field choices, but show significant improvement from extending the baseline. Optimising the sample selection to fit the window function will improve survey efficacy. △ Less

Submitted 15 August, 2022; v1 submitted 8 March, 2022; originally announced March 2022.

Comments: 16 pages, 16 figures. Accepted by MNRAS

Measuring Cosmological Parameters with Type Ia Supernovae in redMaGiC galaxies

Authors: R. Chen, D. Scolnic, E. Rozo, E. S. Rykoff, B. Popovic, R. Kessler, M. Vincenzi, T. M. Davis, P. Armstrong, D. Brout, L. Galbany, L. Kelsey, C. Lidman, A. Möller, B. Rose, M. Sako, M. Sullivan, G. Taylor, P. Wiseman, J. Asorey, A. Carr, C. Conselice, K. Kuehn, G. F. Lewis, E. Macaulay , et al. (60 additional authors not shown)

Abstract: Current and future cosmological analyses with Type Ia Supernovae (SNe Ia) face three critical challenges: i) measuring redshifts from the supernova or its host galaxy; ii) classifying SNe without spectra; and iii) accounting for correlations between the properties of SNe Ia and their host galaxies. We present here a novel approach that addresses each challenge. In the context of the Dark Energy Su… ▽ More Current and future cosmological analyses with Type Ia Supernovae (SNe Ia) face three critical challenges: i) measuring redshifts from the supernova or its host galaxy; ii) classifying SNe without spectra; and iii) accounting for correlations between the properties of SNe Ia and their host galaxies. We present here a novel approach that addresses each challenge. In the context of the Dark Energy Survey (DES), we analyze a SNIa sample with host galaxies in the redMaGiC galaxy catalog, a selection of Luminous Red Galaxies. Photo-$z$ estimates for these galaxies are expected to be accurate to $σ_{Δz/(1+z)}\sim0.02$. The DES-5YR photometrically classified SNIa sample contains approximately 1600 SNe and 125 of these SNe are in redMaGiC galaxies. We demonstrate that redMaGiC galaxies almost exclusively host SNe Ia, reducing concerns with classification uncertainties. With this subsample, we find similar Hubble scatter (to within $\sim0.01$ mag) using photometric redshifts in place of spectroscopic redshifts. With detailed simulations, we show the bias due to using photo-$z$s from redMaGiC host galaxies on the measurement of the dark energy equation-of-state $w$ is up to $Δw \sim 0.01-0.02$. With real data, we measure a difference in $w$ when using redMaGiC photometric redshifts versus spectroscopic redshifts of $Δw = 0.005$. Finally, we discuss how SNe in redMaGiC galaxies appear to be a more standardizable population due to a weaker relation between color and luminosity ($β$) compared to the DES-3YR population by $\sim5σ$; this finding is consistent with predictions that redMaGiC galaxies exhibit lower reddening ratios ($\textrm{R}_\textrm{V}$) than the general population of SN host galaxies. These results establish the feasibility of performing redMaGiC SN cosmology with photometric survey data in the absence of spectroscopic data. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: 15 pages, 8 figures. Submitting to ApJ, comments welcome

Journal ref: ApJ 938 62 (2022)

Large scale quantum chemistry with Tensor Processing Units

Authors: Ryan Pederson, John Kozlowski, Ruyi Song, Jackson Beall, Martin Ganahl, Markus Hauru, Adam G. M. Lewis, Yi Yao, Shrestha Basu Mallick, Volker Blum, Guifre Vidal

Abstract: We demonstrate the use of Google's cloud-based Tensor Processing Units (TPUs) to accelerate and scale up conventional (cubic-scaling) density functional theory (DFT) calculations. Utilizing 512 TPU cores, we accomplish the largest such DFT computation to date, with 247848 orbitals, corresponding to a cluster of 10327 water molecules with 103270 electrons, all treated explicitly. Our work thus pave… ▽ More We demonstrate the use of Google's cloud-based Tensor Processing Units (TPUs) to accelerate and scale up conventional (cubic-scaling) density functional theory (DFT) calculations. Utilizing 512 TPU cores, we accomplish the largest such DFT computation to date, with 247848 orbitals, corresponding to a cluster of 10327 water molecules with 103270 electrons, all treated explicitly. Our work thus paves the way towards accessible and systematic use of conventional DFT, free of any system-specific constraints, at unprecedented scales. △ Less

Submitted 12 December, 2022; v1 submitted 2 February, 2022; originally announced February 2022.

The Dark Energy Survey 5-year photometrically identified Type Ia Supernovae

Authors: A. Möller, M. Smith, M. Sako, M. Sullivan, M. Vincenzi, P. Wiseman, P. Armstrong, J. Asorey, D. Brout, D. Carollo, T. M. Davis, C. Frohmaier, L. Galbany, K. Glazebrook, L. Kelsey, R. Kessler, G. F. Lewis, C. Lidman, U. Malik, R. C. Nichol, D. Scolnic, B. E. Tucker, T. M. C. Abbott, M. Aguena, S. Allam , et al. (58 additional authors not shown)

Abstract: As part of the cosmology analysis using Type Ia Supernovae (SN Ia) in the Dark Energy Survey (DES), we present photometrically identified SN Ia samples using multi-band light-curves and host galaxy redshifts. For this analysis, we use the photometric classification framework SuperNNova (SNN; Möller et al. 2019) trained on realistic DES-like simulations. For reliable classification, we process the… ▽ More As part of the cosmology analysis using Type Ia Supernovae (SN Ia) in the Dark Energy Survey (DES), we present photometrically identified SN Ia samples using multi-band light-curves and host galaxy redshifts. For this analysis, we use the photometric classification framework SuperNNova (SNN; Möller et al. 2019) trained on realistic DES-like simulations. For reliable classification, we process the DES SN programme (DES-SN) data and introduce improvements to the classifier architecture, obtaining classification accuracies of more than 98 per cent on simulations. This is the first SN classification to make use of ensemble methods, resulting in more robust samples. Using photometry, host galaxy redshifts, and a classification probability requirement, we identify 1,863 SNe Ia from which we select 1,484 cosmology-grade SNe Ia spanning the redshift range of 0.07 < z < 1.14. We find good agreement between the light-curve properties of the photometrically-selected sample and simulations. Additionally, we create similar SN Ia samples using two types of Bayesian Neural Network classifiers that provide uncertainties on the classification probabilities. We test the feasibility of using these uncertainties as indicators for out-of-distribution candidates and model confidence. Finally, we discuss the implications of photometric samples and classification methods for future surveys such as Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). △ Less

Submitted 19 July, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 18 pages, 8 figures. Accepted in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 4, August 2022, Pages 5159 5177

Metallicity distribution of the progenitor of the Giant Stellar Stream in the Andromeda Galaxy

Authors: Stanislav Milošević, Miroslav Mićić, Geraint F. Lewis

Abstract: The Giant Stellar Stream (GSS) is a prominent tidal feature in the halo of the Andromeda Galaxy (M31), representing the ongoing destruction of a satellite galaxy. In this paper, we investigate the formation of the GSS through detailed numerical simulations of the tidal disruption of a progenitor system. Assuming that the stream was created in a single merger event between M31 and a dwarf spheroida… ▽ More The Giant Stellar Stream (GSS) is a prominent tidal feature in the halo of the Andromeda Galaxy (M31), representing the ongoing destruction of a satellite galaxy. In this paper, we investigate the formation of the GSS through detailed numerical simulations of the tidal disruption of a progenitor system. Assuming that the stream was created in a single merger event between M31 and a dwarf spheroidal galaxy with stellar mass of $10^{9}M_{\odot}$, we successfully reproduce the dynamical properties of the GSS. As the metallicity distribution along the stream has been well determined from the observations (PAndAS and AMIGA data sets), we use Monte Carlo simulations to reconstruct the original metallicity distribution of the dwarf progenitor. We find that a progenitor dwarf galaxy with a negative radial metallicity gradient, $Δ$FeH = -$0.3 \pm 0.2$, reproduces the observed GSS properties at a time between 2.4 and 2.9 Gyrs into the merger. We also show that the observed double peak metallicity distribution along the stream is a transitory structure caused by unique merger circumstances where two groups of streaming stars are moving in opposite directions, intersecting to produce the peaks. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Comments: 13 pages, 14 figures, MNRAS, accepted

The Hierarchical Structure of Galactic Haloes: Generalised N-Dimensional Clustering with CluSTAR-ND

Authors: William H. Oliver, Pascal J. Elahi, Geraint F. Lewis

Abstract: We present CluSTAR-ND, a fast hierarchical galaxy/(sub)halo finder that produces {\bf Clu}stering {\bf S}tructure via {\bf T}ransformative {\bf A}ggregation and {\bf R}ejection in {\bf N}-{\bf D}imensions. It is designed to improve upon Halo-OPTICS -- an algorithm that automatically detects and extracts significant astrophysical clusters from the 3D spatial positions of simulation particles -- by… ▽ More We present CluSTAR-ND, a fast hierarchical galaxy/(sub)halo finder that produces {\bf Clu}stering {\bf S}tructure via {\bf T}ransformative {\bf A}ggregation and {\bf R}ejection in {\bf N}-{\bf D}imensions. It is designed to improve upon Halo-OPTICS -- an algorithm that automatically detects and extracts significant astrophysical clusters from the 3D spatial positions of simulation particles -- by decreasing run-times, possessing the capability for metric adaptivity, and being readily applicable to data with any number of features. We directly compare these algorithms and find that not only does CluSTAR-ND produce a similarly robust clustering structure, it does so in a run-time that is at least $3$ orders of magnitude faster. In optimising CluSTAR-ND's clustering performance, we have also carefully calibrated $4$ of the $7$ CluSTAR-ND parameters which -- unless specified by the user -- will be automatically and optimally chosen based on the input data. We conclude that CluSTAR-ND is a robust astrophysical clustering algorithm that can be leveraged to find stellar satellite groups on large synthetic or observational data sets. △ Less

Submitted 15 June, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 19 pages, 11 figures, 1 table. Accepted to MNRAS

Large Scale Distributed Linear Algebra With Tensor Processing Units

Authors: Adam G. M. Lewis, Jackson Beall, Martin Ganahl, Markus Hauru, Shrestha Basu Mallick, Guifre Vidal

Abstract: We have repurposed Google Tensor Processing Units (TPUs), application-specific chips developed for machine learning, into large-scale dense linear algebra supercomputers. The TPUs' fast inter-core interconnects (ICI)s, physically two-dimensional network topology, and high-bandwidth memory (HBM) permit distributed matrix multiplication algorithms to rapidly become computationally bound. In this reg… ▽ More We have repurposed Google Tensor Processing Units (TPUs), application-specific chips developed for machine learning, into large-scale dense linear algebra supercomputers. The TPUs' fast inter-core interconnects (ICI)s, physically two-dimensional network topology, and high-bandwidth memory (HBM) permit distributed matrix multiplication algorithms to rapidly become computationally bound. In this regime, the matrix-multiply units (MXU)s dominate the runtime, yielding impressive scaling, performance, and raw size: operating in float32 precision, a full 2048-core pod of third generation TPUs can multiply two matrices with linear size $N= 220= 1 048 576$ in about 2 minutes. Via curated algorithms emphasizing large, single-core matrix multiplications, other tasks in dense linear algebra can similarly scale. As examples, we present (i) QR decomposition; (ii) resolution of linear systems; and (iii) the computation of matrix functions by polynomial iteration, demonstrated by the matrix polar factorization. △ Less

Submitted 16 December, 2021; originally announced December 2021.

Comments: 12 pages, 8 figures

The GALAH Survey: Improving our understanding of confirmed and candidate planetary systems with large stellar surveys

Authors: Jake T. Clark, Duncan J. Wright, Robert A. Wittenmyer, Jonathan Horner, Natalie R. Hinkel, Mathieu Clerté, Brad D. Carter, Sven Buder, Michael R. Hayden, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Katharine J. Schlesinger, Sanjib Sharma, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Tomaž Zwitter , et al. (2 additional authors not shown)

Abstract: Pioneering photometric, astrometric, and spectroscopic surveys are helping exoplanetary scientists better constrain the fundamental properties of stars within our galaxy, and the planets these stars host. In this study, we use the third data release from the stellar spectroscopic GALAH Survey, coupled with astrometric data of eDR3 from the \textit{Gaia} satellite, and other data from NASA's Exopla… ▽ More Pioneering photometric, astrometric, and spectroscopic surveys are helping exoplanetary scientists better constrain the fundamental properties of stars within our galaxy, and the planets these stars host. In this study, we use the third data release from the stellar spectroscopic GALAH Survey, coupled with astrometric data of eDR3 from the \textit{Gaia} satellite, and other data from NASA's Exoplanet Archive, to refine our understanding of 279 confirmed and candidate exoplanet host stars and their exoplanets. This homogenously analysed data set comprises 105 confirmed exoplanets, along with 146 K2 candidates, 95 TESS Objects of Interest (TOIs) and 52 Community TOIs (CTOIs). Our analysis significantly shifts several previously (unknown) planet parameters while decreasing the uncertainties for others; Our radius estimates suggest that 35 planet candidates are more likely brown dwarfs or stellar companions due to their new radius values. We are able to refine the radii and masses of WASP-47 e, K2-106 b, and CoRoT-7 b to their most precise values yet, to less than 2.3\% and 8.5\% respectively. We also use stellar rotational values from GALAH to show that most planet candidates will have mass measurements that will be tough to obtain with current ground-based spectrographs. With GALAH's chemical abundances, we show through chemo-kinematics that there are five planet-hosts that are associated with the galaxy's thick disc, including NGTS-4, K2-183 and K2-337. Finally, we show there is no statistical difference between the chemical properties of hot Neptune and hot rocky exoplanet hosts, with the possibility that short-period rocky worlds might be the remnant cores of hotter, gaseous worlds. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 21 pages with 11 figures and 5 tables. Accepted into MNRAS with appendix tables available within this arXiv submission

Maximum Likelihood Estimation of Differentiated Products Demand Systems

Authors: Greg Lewis, Bora Ozaltun, Georgios Zervas

Abstract: We discuss estimation of the differentiated products demand system of Berry et al (1995) (BLP) by maximum likelihood estimation (MLE). We derive the maximum likelihood estimator in the case where prices are endogenously generated by firms that set prices in Bertrand-Nash equilibrium. In Monte Carlo simulations the MLE estimator outperforms the best-practice GMM estimator on both bias and mean squa… ▽ More We discuss estimation of the differentiated products demand system of Berry et al (1995) (BLP) by maximum likelihood estimation (MLE). We derive the maximum likelihood estimator in the case where prices are endogenously generated by firms that set prices in Bertrand-Nash equilibrium. In Monte Carlo simulations the MLE estimator outperforms the best-practice GMM estimator on both bias and mean squared error when the model is correctly specified. This remains true under some forms of misspecification. In our simulations, the coverage of the ML estimator is close to its nominal level, whereas the GMM estimator tends to under-cover. We conclude the paper by estimating BLP on the car data used in the original Berry et al (1995) paper, obtaining similar estimates with considerably tighter standard errors. △ Less

Submitted 24 November, 2021; originally announced November 2021.

The Dark Energy Survey Supernova Program: Cosmological biases from supernova photometric classification

Authors: M. Vincenzi, M. Sullivan, A. Möller, P. Armstrong, B. A. Bassett, D. Brout, D. Carollo, A. Carr, T. M. Davis, C. Frohmaier, L. Galbany, K. Glazebrook, O. Graur, L. Kelsey, R. Kessler, E. Kovacs, G. F. Lewis, C. Lidman, U. Malik, R. C. Nichol, B. Popovic, M. Sako, D. Scolnic, M. Smith, G. Taylor , et al. (59 additional authors not shown)

Abstract: Cosmological analyses of samples of photometrically-identified Type Ia supernovae (SNe Ia) depend on understanding the effects of 'contamination' from core-collapse and peculiar SN Ia events. We employ a rigorous analysis on state-of-the-art simulations of photometrically identified SN Ia samples and determine cosmological biases due to such 'non-Ia' contamination in the Dark Energy Survey (DES) 5… ▽ More Cosmological analyses of samples of photometrically-identified Type Ia supernovae (SNe Ia) depend on understanding the effects of 'contamination' from core-collapse and peculiar SN Ia events. We employ a rigorous analysis on state-of-the-art simulations of photometrically identified SN Ia samples and determine cosmological biases due to such 'non-Ia' contamination in the Dark Energy Survey (DES) 5-year SN sample. As part of the analysis, we test on our DES simulations the performance of SuperNNova, a photometric SN classifier based on recurrent neural networks. Depending on the choice of non-Ia SN models in both the simulated data sample and training sample, contamination ranges from 0.8-3.5 %, with the efficiency of the classification from 97.7-99.5 %. Using the Bayesian Estimation Applied to Multiple Species (BEAMS) framework and its extension 'BEAMS with Bias Correction' (BBC), we produce a redshift-binned Hubble diagram marginalised over contamination and corrected for selection effects and we use it to constrain the dark energy equation-of-state, $w$. Assuming a flat universe with Gaussian $Ω_M$ prior of $0.311\pm0.010$, we show that biases on $w$ are $<0.008$ when using SuperNNova and accounting for a wide range of non-Ia SN models in the simulations. Systematic uncertainties associated with contamination are estimated to be at most $σ_{w, \mathrm{syst}}=0.004$. This compares to an expected statistical uncertainty of $σ_{w,\mathrm{stat}}=0.039$ for the DES-SN sample, thus showing that contamination is not a limiting uncertainty in our analysis. We also measure biases due to contamination on $w_0$ and $w_a$ (assuming a flat universe), and find these to be $<$0.009 in $w_0$ and $<$0.108 in $w_a$, hence 5 to 10 times smaller than the statistical uncertainties expected from the DES-SN sample. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Simulation of quantum many-body dynamics with Tensor Processing Units: Floquet prethermalization

Authors: Alan Morningstar, Markus Hauru, Jackson Beall, Martin Ganahl, Adam G. M. Lewis, Vedika Khemani, Guifre Vidal

Abstract: Tensor Processing Units (TPUs) are specialized hardware accelerators developed by Google to support large-scale machine-learning tasks, but they can also be leveraged to accelerate and scale other linear-algebra-intensive computations. In this paper we demonstrate the usage of TPUs for massively parallel, classical simulations of quantum many-body dynamics on long timescales. We apply our methods… ▽ More Tensor Processing Units (TPUs) are specialized hardware accelerators developed by Google to support large-scale machine-learning tasks, but they can also be leveraged to accelerate and scale other linear-algebra-intensive computations. In this paper we demonstrate the usage of TPUs for massively parallel, classical simulations of quantum many-body dynamics on long timescales. We apply our methods to study the phenomenon of Floquet prethermalization, i.e., exponentially slow heating in quantum spin chains subject to high-frequency periodic driving. We simulate the dynamics of L=34 qubits for over $10^5$ Floquet periods, corresponding to circuits with millions of two-qubit gates. The circuits simulated have no additional symmetries and represent a pure-state evolution in the full $2^L$-dimensional Hilbert space. This is achieved by distributing the computation over 128 TPU cores. On that size TPU cluster, we find speedups in wall-clock runtime of 230x and 15x when compared to reference CPU and single-GPU simulations, respectively, for shorter 30-qubit simulations that can be handled by all three platforms. We study the computational cost of the simulations, as a function of both the number of qubits and the number of TPU cores used, up to our maximum capacity of L=40 qubits, which requires a ``full pod" of 2048 TPU cores with tens of terabytes of memory in total. For these simulations, an 8-TPU-core machine is comparable to a single A100 GPU, and thus the full TPU pod is comparable to a machine with hundreds of GPUs. However, the TPU pod is more energy and cost efficient, and readily accessible (via Google Cloud), unlike such large many-GPU configurations. We also study the accumulation of numerical error as a function of circuit depth in very deep circuits. Our work demonstrates that TPUs can offer significant advantages for state-of-the-art simulations of quantum many-body dynamics. △ Less

Submitted 11 February, 2022; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: 12 pages, 7 figures; v2 contains substantial improvements including GPU simulations

Journal ref: PRX Quantum 3, 020331 (2022)

A barred Milky Way surrogate from an N-body simulation

Authors: T. Tepper-Garcia, J. Bland-Hawthorn, E. Vasiliev, E. Athanassoula, O. Gerhard, A. Quillen, P. McMillan, K. Freeman, G. F. Lewis, R. Teyssier, S. Sharma, M. R. Hayden, S. Buder

Abstract: We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a se… ▽ More We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a self-consistent dynamical model of the MW defined by the best available parameters for the dark matter halo, the stellar disc and the bulge. For the known (stellar and gas) disc mass (4.5 x 10^10 Msun) and disc mass fraction at R ~ 2.2 R_d (f_d ~ 0.3 - 0.6), the low mass limit does not yield a bar in a Hubble time. The high mass limit adopted here produces a box/peanut bar within about 2 Gyr with the correct mass (~10^10 Msun), size (~5 kpc) and peak pattern speed (~ 40-45 km/s/kpc). In agreement with earlier work, the bar formation timescale scales inversely with f_d (i.e. log [T/Gyr] ~ 0.60/f_d - 0.83 for 1 < f_d < 0.3). The disc radial heating is strong, but, in contrast to earlier claims, we find that disc vertical heating outside of the box/peanut bulge structure is negligible. The synthetic barred MW exhibits long-term stability, except for the slow decline (roughly -2 km/s/kpc/Gyr) of the bar pattern speed, consistent with recent estimates. If our model is indicative of the Milky Way, we estimate that the bar first emerged 3-4 Gyr ago. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: We welcome comments and missing references

Mapping the cosmic mass distribution with stacked weak gravitational lensing and Doppler lensing

Authors: Md Rasel Hossen, Sonia Akter Ema, Krzysztof Bolejko, Geraint F. Lewis

Abstract: Dark matter halos represent the highest density peaks in the matter distribution. Conversely, cosmic voids are under-dense patches of the universe. Probing the mass distribution of the universe requires various approaches, including weak gravitational lensing that subtly modifies the shape of distant sources, and Doppler lensing that changes the apparent size and magnitude of objects due to peculi… ▽ More Dark matter halos represent the highest density peaks in the matter distribution. Conversely, cosmic voids are under-dense patches of the universe. Probing the mass distribution of the universe requires various approaches, including weak gravitational lensing that subtly modifies the shape of distant sources, and Doppler lensing that changes the apparent size and magnitude of objects due to peculiar velocities. In this work, we adopt both gravitational and Doppler lensing effects to study the underlying matter distribution in and around cosmic voids/halos. We use the relativistic $N$-body code \texttt{gevolution}, to generate the mass perturbations and develop a new ray-tracing code that relies on the design of the ray bundle method. We consider three categories of halo masses and void radii, and extract the cosmological information by stacking weak-lensing and Doppler lensing signals around voids/halos. The results of this paper show that the most optimal strategy that combines both gravitational and Doppler lensing effects to map the mass distribution should focus on the redshift range $z\approx 0.3-0.4$. The recommendation of this paper is that future spectroscopic surveys should focus on these redshifts and utilise the gravitational and Doppler lensing techniques to extract information about underlying matter distribution across the cosmic web, especially inside cosmic voids. This could provide a complimentary cosmological analysis for ongoing or future low-redshift spectroscopic surveys. △ Less

Submitted 8 December, 2021; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: 13 pages, 10 figures, corresponds to the published version

The density distributions of cosmic structures: impact of the local environment on weak-lensing convergence

Authors: Sonia Akter Ema, Md Rasel Hossen, Krzysztof Bolejko, Geraint F. Lewis

Abstract: Whilst the underlying assumption of the Friedman-Lemaître-Robertson-Walker (FLRW) cosmological model is that matter is homogeneously distributed throughout the universe, gravitational influences over the life of the universe have resulted in mass clustered on a range of scales. Hence we expect that, in our inhomogeneous universe, the view of an observer will be influenced by the location and local… ▽ More Whilst the underlying assumption of the Friedman-Lemaître-Robertson-Walker (FLRW) cosmological model is that matter is homogeneously distributed throughout the universe, gravitational influences over the life of the universe have resulted in mass clustered on a range of scales. Hence we expect that, in our inhomogeneous universe, the view of an observer will be influenced by the location and local environment. Here we analyse the one-point probability distribution functions and angular power spectra of weak-lensing (WL) convergence and magnification numerically to investigate the influence of our local environment on WL statistics in relativistic $N$-body simulations. To achieve this, we numerically solve the null geodesic equations which describe the propagation of light bundles backwards in time from today, and develop a ray-tracing algorithm, and from these calculate various WL properties. Our findings demonstrate how cosmological observations of large-scale structure through WL can be impacted by the locality of the observer. We also calculate the constraints on the cosmological parameters as a function of redshift from the theoretical and numerical study of the angular power spectrum of WL convergence. This study concludes the minimal redshift for the constraint on the parameter $Ω_m$ ($H_0$) is $z \sim 0.2$ $(z \sim 0.6 )$ beyond which the local environment's effect is negligible and the data from WL surveys are more meaningful above that redshift. The outcomes of this study will have direct consequences for future surveys, where percent-level-precision is necessary. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 11 pages, 10 figures, Accepted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 509, Issue 2, January 2022, Pages 3004

Collaboration Challenges in Building ML-Enabled Systems: Communication, Documentation, Engineering, and Process

Authors: Nadia Nahar, Shurui Zhou, Grace Lewis, Christian Kästner

Abstract: The introduction of machine learning (ML) components in software projects has created the need for software engineers to collaborate with data scientists and other specialists. While collaboration can always be challenging, ML introduces additional challenges with its exploratory model development process, additional skills and knowledge needed, difficulties testing ML systems, need for continuous… ▽ More The introduction of machine learning (ML) components in software projects has created the need for software engineers to collaborate with data scientists and other specialists. While collaboration can always be challenging, ML introduces additional challenges with its exploratory model development process, additional skills and knowledge needed, difficulties testing ML systems, need for continuous evolution and monitoring, and non-traditional quality requirements such as fairness and explainability. Through interviews with 45 practitioners from 28 organizations, we identified key collaboration challenges that teams face when building and deploying ML systems into production. We report on common collaboration points in the development of production ML systems for requirements, data, and integration, as well as corresponding team patterns and challenges. We find that most of these challenges center around communication, documentation, engineering, and process and collect recommendations to address these challenges. △ Less

Submitted 10 February, 2022; v1 submitted 19 October, 2021; originally announced October 2021.

Comments: 22 pages, 10 figures, 5 tables

$S^5$: The Orbital and Chemical Properties of One Dozen Stellar Streams

Authors: Ting S. Li, Alexander P. Ji, Andrew B. Pace, Denis Erkal, Sergey E. Koposov, Nora Shipp, Gary S. Da Costa, Lara R. Cullinane, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Jeffrey D. Simpson, Daniel B. Zucker, Peter S. Ferguson, Sarah L. Martell, Joss Bland-Hawthorn, Eduardo Balbinot, Kiyan Tavangar, Alex Drlica-Wagner, Gayandhi M. De Silva, Joshua D. Simon, S5 Collaboration

Abstract: We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors, using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ($S^5$), proper motions from $Gaia$ EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full… ▽ More We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors, using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ($S^5$), proper motions from $Gaia$ EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full 6D kinematics and metallicities. All streams have heliocentric distances between ${\sim}10-50$ kpc. The velocity and metallicity dispersions show that half of the stream progenitors were disrupted dwarf galaxies (DGs), while the other half originated from disrupted globular clusters (GCs), hereafter referred to as DG and GC streams. Based on the mean metallicities of the streams and the mass-metallicity relation, the luminosities of the progenitors of the DG streams range between Carina and Ursa Major I ($-9.5\lesssim M_V\lesssim-5.5$). Four of the six GC streams have mean metallicities of [Fe/H]$< -2$, more metal-poor than typical Milky Way (MW) GCs at similar distances. Interestingly, the 300S and Jet GC streams are the only streams on retrograde orbits in our dozen stream sample. Finally, we compare the orbital properties of the streams with known DGs and GCs in the MW, finding several possible associations. Some streams appear to have been accreted with the recently discovered Gaia-Enceladus-Sausage system, and others suggest that GCs were formed in and accreted together with the progenitors of DG streams whose stellar masses are similar to Draco to Carina ($\sim10^5-10^6M_\odot$). △ Less

Submitted 2 January, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: 24 pages, 8 Figures, 2 Tables. Accepted for publication in ApJ

From the Fire: A Deeper Look at the Phoenix Stream

Authors: K. Tavangar, P. Ferguson, N. Shipp, A. Drlica-Wagner, S. Koposov, D. Erkal, E. Balbinot, J. García-Bellido, K. Kuehn, G. F. Lewis, T. S. Li, S. Mau, A. B. Pace, A. H. Riley, T. M. C. Abbott, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Annis, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero , et al. (42 additional authors not shown)

Abstract: We use six years of data from the Dark Energy Survey to perform a detailed photometric characterization of the Phoenix stellar stream, a 15-degree long, thin, dynamically cold, low-metallicity stellar system in the southern hemisphere. We use natural splines, a non-parametric modeling technique, to simultaneously fit the stream track, width, and linear density. This updated stream model allows us… ▽ More We use six years of data from the Dark Energy Survey to perform a detailed photometric characterization of the Phoenix stellar stream, a 15-degree long, thin, dynamically cold, low-metallicity stellar system in the southern hemisphere. We use natural splines, a non-parametric modeling technique, to simultaneously fit the stream track, width, and linear density. This updated stream model allows us to improve measurements of the heliocentric distance ($17.4 \pm 0.1\,{\rm (stat.)} \pm 0.8\,{\rm (sys.)}$ kpc) and distance gradient ($-0.009 \pm 0.006$ kpc deg$^{-1}$) of Phoenix, which corresponds to a small change of $0.13 \pm 0.09$ kpc in heliocentric distance along the length of the stream. We measure linear intensity variations on degree scales, as well as deviations in the stream track on $\sim 2$-degree scales, suggesting that the stream may have been disturbed during its formation and/or evolution. We recover three peaks and one gap in linear intensity along with fluctuations in the stream track. Compared to other thin streams, the Phoenix stream shows more fluctuations and, consequently, the study of Phoenix offers a unique perspective on gravitational perturbations of stellar streams. We discuss possible sources of perturbations to Phoenix including baryonic structures in the Galaxy and dark matter subhalos. △ Less

Submitted 15 February, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 18 Pages, 7 Figures, Matches published version

Journal ref: The Astrophysical Journal, Volume 925, Number 2, (2022) 118-132

Deriving the bulk properties of solar wind electrons observed by Solar Orbiter: A preliminary study of electron plasma thermodynamics

Authors: Georgios Nicolaou, Robert T. Wicks, Christopher J. Owen, Dhiren O. Kataria, Anekallu Chandrasekhar, Gethyn R. Lewis, Daniel Verscharen, Vito Fortunato, Gennaro Mele, Rossana DeMarco, Roberto Bruno

Abstract: We demonstrate the calculation of solar wind electron bulk parameters from recent observations by Solar Wind Analyser Electron Analyser System on board Solar Orbiter. We use our methods to derive the electron bulk parameters in a time interval of a few hours. We attempt a preliminary examination of the polytropic behavior of the electrons by analyzing the derived electron density and temperature.… ▽ More We demonstrate the calculation of solar wind electron bulk parameters from recent observations by Solar Wind Analyser Electron Analyser System on board Solar Orbiter. We use our methods to derive the electron bulk parameters in a time interval of a few hours. We attempt a preliminary examination of the polytropic behavior of the electrons by analyzing the derived electron density and temperature. Moreover, we discuss the challenges in analyzing the observations due to the spacecraft charging and photo-electron contamination in the energy range < 10 eV. Aims: We derive bulk parameters of thermal solar wind electrons by analyzing Solar Orbiter observations and we investigate if there is any typical polytropic model that applies to the electron density and temperature fluctuations. Methods: We use the appropriate transformations to convert the observations to velocity distribution functions in the instrument frame. We then derive the electron bulk parameters by a) calculating the statistical moments of the constructed velocity distribution functions and b) by fitting the constructed distributions with analytical expressions. We firstly test our methods by applying them to an artificial data-set, which we produce by using the forward modeling technique. Results: The forward model validates the analysis techniques which we use to derive the electron bulk parameters. The calculation of the statistical moments and the fitting method determines bulk parameters that are identical within uncertainty to the input parameters we use to simulate the plasma electrons in the first place. An application of our analysis technique to the data reveals a nearly isothermal electron "core". The results are affected by the spacecraft potential and the photo-electron contamination, which we need to characterize in detail in future analyses. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Journal ref: A&A 656, A10 (2021)

The GALAH Survey: Chemical tagging and chrono-chemodynamics of accreted halo stars with GALAH+ DR3 and $Gaia$ eDR3

Authors: Sven Buder, Karin Lind, Melissa K. Ness, Diane K. Feuillet, Danny Horta, Stephanie Monty, Tobias Buck, Thomas Nordlander, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Michael R. Hayden, Janez Kos, Sarah L. Martell, Geraint F. Lewis, Jane Lin, Katharine. J. Schlesinger, Sanjib Sharma, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Tomaz Zwitter, Ioana Ciuca , et al. (5 additional authors not shown)

Abstract: Since the advent of $Gaia$ astrometry, it is possible to identify massive accreted systems within the Galaxy through their unique dynamical signatures. One such system, $Gaia$-Sausage-Enceladus (GSE), appears to be an early "building block" given its virial mass $> 10^{10}\,\mathrm{M_\odot}$ at infall ($z\sim1-3$). In order to separate the progenitor population from the background stars, we invest… ▽ More Since the advent of $Gaia$ astrometry, it is possible to identify massive accreted systems within the Galaxy through their unique dynamical signatures. One such system, $Gaia$-Sausage-Enceladus (GSE), appears to be an early "building block" given its virial mass $> 10^{10}\,\mathrm{M_\odot}$ at infall ($z\sim1-3$). In order to separate the progenitor population from the background stars, we investigate its chemical properties with up to 30 element abundances from the GALAH+ Survey Data Release 3 (DR3). To inform our choice of elements for purely chemically selecting accreted stars, we analyse 4164 stars with low-$α$ abundances and halo kinematics. These are most different to the Milky Way stars for abundances of Mg, Si, Na, Al, Mn, Fe, Ni, and Cu. Based on the significance of abundance differences and detection rates, we apply Gaussian mixture models to various element abundance combinations. We find the most populated and least contaminated component, which we confirm to represent GSE, contains 1049 stars selected via [Na/Fe] vs. [Mg/Mn] in GALAH+ DR3. We provide tables of our selections and report the chrono-chemodynamical properties (age, chemistry, and dynamics). Through a previously reported clean dynamical selection of GSE stars, including $30 < \sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} < 55$, we can characterise an unprecedented 24 abundances of this structure with GALAH+ DR3. Our chemical selection allows us to prevent circular reasoning and characterise the dynamical properties of the GSE, for example mean $\sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} = 26_{-14}^{+9}$. We find only $(29\pm1)\%$ of the GSE stars within the clean dynamical selection region. Our methodology will improve future studies of accreted structures and their importance for the formation of the Milky Way. △ Less

Submitted 5 January, 2022; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: 29 pages, 18 figures, 6 tables; publicshed in MNRAS; all code to recreate results (based on public data described in arXiv:2011.02505) and figures available at https://github.com/svenbuder/Accreted-stars-in-GALAH-DR3

Journal ref: MNRAS, 510, 2407 (2022)

Signature of a massive rotating metal-poor star imprinted in the Phoenix stellar stream

Authors: Andrew R. Casey, Alexander P. Ji, Terese T. Hansen, Ting S. Li, Sergey E. Koposov, Gary S. Da Costa, Joss Bland-Hawthorn, Lara Cullinane, Denis Erkal, Geraint F. Lewis, Kyler Kuehn, Dougal Mackey, Sarah L. Martell, Andrew B. Pace, Jeffrey D. Simpson, Daniel B. Zucker

Abstract: The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low mass globular cluster. In this work we use Magellan/MIKE high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ($σ(\mathrm{[Fe~II/H]}) = 0.04^{+0.11}_{-0.03}$) and… ▽ More The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low mass globular cluster. In this work we use Magellan/MIKE high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ($σ(\mathrm{[Fe~II/H]}) = 0.04^{+0.11}_{-0.03}$) and a large peak-to-peak range in [Na/Fe] and [Al/Fe] abundance ratios, consistent with the light element abundance patterns seen in the most metal-poor globular clusters. However, unlike any other globular cluster, we also find an intrinsic spread in [Sr II/Fe] spanning $\sim$1 dex, while [Ba II/Fe] shows nearly no intrinsic spread ($σ(\mathrm{[Ba~II/H]}) = {0.03}^{+0.10}_{-0.02}$). This abundance signature is best interpreted as slow neutron capture element production from a massive fast-rotating metal-poor star ($15-20 \mathrm{M}_\odot$, $v_\mathrm{ini}/v_\mathrm{crit} = 0.4$, $[\mathrm{Fe/H}] = -3.8$). The low inferred cluster mass suggests the system would have been unable to retain supernovae ejecta, implying that any massive fast-rotating metal-poor star that enriched the interstellar medium must have formed and evolved before the globular cluster formed. Neutron capture element production from asymptotic giant branch stars or magneto-rotational instabilities in core-collapse supernovae provide poor fits to the observations. We also report one Phoenix stream star to be a lithium-rich giant ($A(\mathrm{Li}) = 3.1 \pm 0.1$). At $[\mathrm{Fe/H}] = -2.93$ it is among the most metal-poor lithium-rich giants known. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: Accepted to ApJ 2021-07-09

Multi-wavelength Optical and NIR Variability Analysis of the Blazar PKS 0027-426

Authors: E. Guise, S. F. Hönig, T. Almeyda, K. Horne, M. Kishimoto, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Asorey, M. Banerji, E. Bertin, B. Boulderstone, D. Brooks, D. L. Burke, A. Carnero Rosell, D. Carollo, M. Carrasco Kind, J. Carretero, M. Costanzi, L. N. da Costa, T. M. Davis, J. De Vicente, P. Doel, S. Everett, I. Ferrero , et al. (40 additional authors not shown)

Abstract: We present multi-wavelength spectral and temporal variability analysis of PKS 0027-426 using optical griz observations from DES (Dark Energy Survey) between 2013-2018 and VOILETTE (VEILS Optical Light curves of Extragalactic TransienT Events) between 2018-2019 and near infrared (NIR) JKs observations from VEILS (VISTAExtragalactic Infrared Legacy Survey) between 2017-2019. Multiple methods of cros… ▽ More We present multi-wavelength spectral and temporal variability analysis of PKS 0027-426 using optical griz observations from DES (Dark Energy Survey) between 2013-2018 and VOILETTE (VEILS Optical Light curves of Extragalactic TransienT Events) between 2018-2019 and near infrared (NIR) JKs observations from VEILS (VISTAExtragalactic Infrared Legacy Survey) between 2017-2019. Multiple methods of cross-correlation of each combination of light curve provides measurements of possible lags between optical-optical, optical-NIR, and NIR-NIR emission, for each observation season and for the entire observational period. Inter-band time lag measurements consistently suggest either simultaneous emission or delays between emission regions on timescales smaller than the cadences of observations. The colour-magnitude relation between each combination of filters was also studied to determine the spectral behaviour of PKS 0027-426. Our results demonstrate complex colour behaviour that changes between bluer when brighter (BWB), stable when brighter (SWB) and redder when brighter (RWB) trends over different timescales and using different combinations of optical filters. Additional analysis of the optical spectra is performed to provide further understanding of this complex spectral behaviour. △ Less

Submitted 25 November, 2021; v1 submitted 30 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS. 33 pages, 28 figures, 11 tables

Clock Transitions Guard Against Spin Decoherence in Singlet Fission

Authors: Sina G. Lewis, Kori E. Smyser, Joel D. Eaves

Abstract: Short coherence times present a primary obstacle in quantum computing and sensing applications. In atomic systems, clock transitions (CTs), formed from avoided crossings in an applied Zeeman field, can substantially increase coherence times. We show how CTs can dampen intrinsic and extrinsic sources of quantum noise in molecules. Conical intersections between two periodic potentials form CTs in el… ▽ More Short coherence times present a primary obstacle in quantum computing and sensing applications. In atomic systems, clock transitions (CTs), formed from avoided crossings in an applied Zeeman field, can substantially increase coherence times. We show how CTs can dampen intrinsic and extrinsic sources of quantum noise in molecules. Conical intersections between two periodic potentials form CTs in electron paramagnetic resonance experiments of the spin-polarized singlet fission photoproduct. We report on a pair of CTs for a two-chromophore molecule in terms of the Zeeman field strength, molecular orientation relative to the field, and molecular geometry. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 18 pages, 3 figures --- Submitted to Journal of Chemical Physics

The K2 Galactic Archaeology Program Data Release 3: Age-abundance patterns in C1-C8 and C10-C18

Authors: Joel C. Zinn, Dennis Stello, Yvonne Elsworth, Rafael A. García, Thomas Kallinger, Savita Mathur, Benoît Mosser, Marc Hon, Lisa Bugnet, Caitlin Jones, Claudia Reyes, Sanjib Sharma, Ralph Schönrich, Jack T. Warfield, Rodrigo Luger, Andrew Vanderburg, Chiaki Kobayashi, Marc H. Pinsonneault, Jennifer A. Johnson, Daniel Huber, Sven Buder, Meridith Joyce, Joss Bland-Hawthorn, Luca Casagrande, Geraint F. Lewis , et al. (6 additional authors not shown)

Abstract: We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1-C8 and C10-C18. We provide asteroseismic radius and mass coefficients, $κ_R$ and $κ_M$, for $\sim 19,000$ red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 s… ▽ More We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1-C8 and C10-C18. We provide asteroseismic radius and mass coefficients, $κ_R$ and $κ_M$, for $\sim 19,000$ red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 stellar parameters are calibrated to be on an absolute parallactic scale based on Gaia DR2, with red giant branch and red clump evolutionary state classifications provided via a machine-learning approach. Combining these stellar parameters with GALAH DR3 spectroscopy, we determine asteroseismic ages with precisions of $\sim 20-30\%$ and compare age-abundance relations to Galactic chemical evolution models among both low- and high-$α$ populations for $α$, light, iron-peak, and neutron-capture elements. We confirm recent indications in the literature of both increased Ba production at late Galactic times, as well as significant contribution to r-process enrichment from prompt sources associated with, e.g., core-collapse supernovae. With an eye toward other Galactic archaeology applications, we characterize K2 GAP DR3 uncertainties and completeness using injection tests, suggesting K2 GAP DR3 is largely unbiased in mass/age and with uncertainties of $2.9\%\,(\rm{stat.})\,\pm0.1\%\,(\rm{syst.})$ & $6.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ in $κ_R$ & $κ_M$ for red giant branch stars and $4.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ & $11\%\,(\rm{stat.})\,\pm0.9\%\,(\rm{syst.})$ for red clump stars. We also identify percent-level asteroseismic systematics, which are likely related to the time baseline of the underlying data, and which therefore should be considered in TESS asteroseismic analysis. △ Less

Submitted 10 March, 2022; v1 submitted 11 August, 2021; originally announced August 2021.

Comments: Published in ApJ; machine-readable tables available through the open access online journal article

Journal ref: ApJ 926 191 (2022)

Measuring the Mass of the Large Magellanic Cloud with Stellar Streams Observed by ${S}^5$

Authors: Nora Shipp, Denis Erkal, Alex Drlica-Wagner, Ting S. Li, Andrew B. Pace, Sergey E. Koposov, Lara R. Cullinane, Gary S. Da Costa, Alexander P. Ji, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Jeffrey D. Simpson, Zhen Wan, Daniel B. Zucker, Joss Bland-Hawthorn, Peter S. Ferguson, Sophia Lilleengen, S5 Collaboration

Abstract: Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey (${S}^5$), Gaia EDR3, and the Dark Energy Survey (DES), to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble o… ▽ More Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey (${S}^5$), Gaia EDR3, and the Dark Energy Survey (DES), to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble of streams, we find a mass of the LMC ranging from 14 to $19 \times 10^{10}\ \mathrm{M}_{\odot}$, probed over a range of closest approach times and distances. With the most constraining stream (Orphan-Chenab), we measure an LMC mass of $18.8^{+ 3.5}_{- 4.0} \times 10^{10}\ \mathrm{M}_{\odot}$, probed at a closest approach time of 310 Myr and a closest approach distance of 25.4 kpc. This mass is compatible with previous measurements, showing that a consistent picture is emerging of the LMC's influence on structures in the Milky Way. Using this sample of streams, we find that the LMC's effect depends on the relative orientation of the stream and LMC at their point of closest approach. To better understand this, we present a simple model based on the impulse approximation and we show that the LMC's effect depends both on the magnitude of the velocity kick imparted to the stream and the direction of this kick. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: 23 pages, 10 figures, submitted to ApJ

Whistler instability driven by the sunward electron deficit in the solar wind

Authors: Laura Berčič, Daniel Verscharen, Christopher J. Owen, Lucas Colomban, Matthieu Kretzschmar, Thomas Chust, Milan Maksimović, Dhiren Kataria, Etienne Behar, Matthieu Berthomier, Roberto Bruno, Vito Fortunato, Christopher W. Kelly, Yuri. V. Khotyaintsev, Gethyn R. Lewis, Stefano Livi, Philippe Louarn, Gennaro Mele, Georgios Nicolaou, Gillian Watson, Robert T. Wicks

Abstract: Solar wind electrons play an important role in the energy balance of the solar wind acceleration by carrying energy into interplanetary space in the form of electron heat flux. The heat flux is stored in the complex electron velocity distribution functions (VDFs) shaped by expansion, Coulomb collisions, and field-particle interactions. We investigate how the suprathermal electron deficit in the an… ▽ More Solar wind electrons play an important role in the energy balance of the solar wind acceleration by carrying energy into interplanetary space in the form of electron heat flux. The heat flux is stored in the complex electron velocity distribution functions (VDFs) shaped by expansion, Coulomb collisions, and field-particle interactions. We investigate how the suprathermal electron deficit in the anti-strahl direction, which was recently discovered in the near-Sun solar wind, drives a kinetic instability and creates whistler waves with wave vectors that are quasi-parallel to the direction of the background magnetic field. We combine high-cadence measurements of electron pitch-angle distribution functions and electromagnetic waves provided by Solar Orbiter during its first orbit. Our case study is based on a burst-mode data interval from the Electrostatic Analyser System (SWA-EAS) at a distance of 112 $R_S$ (0.52 au) from the Sun, during which several whistler wave packets were detected by Solar Orbiter's Radio and Plasma Waves (RPW) instrument. The sunward deficit creates kinetic conditions under which the quasi-parallel whistler wave is driven unstable. We directly test our predictions for the existence of these waves through solar wind observations. We find whistler waves that are quasi-parallel and almost circularly polarised, propagating away from the Sun, coinciding with a pronounced sunward deficit in the electron VDF. The cyclotron-resonance condition is fulfilled for electrons moving in the direction opposite to the direction of wave propagation, with energies corresponding to those associated with the sunward deficit. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Journal ref: A&A 656, A31 (2021)

Dim but not entirely dark: Extracting the Galactic Center Excess' source-count distribution with neural nets

Authors: Florian List, Nicholas L. Rodd, Geraint F. Lewis

Abstract: The two leading hypotheses for the Galactic Center Excess (GCE) in the $\textit{Fermi}$ data are an unresolved population of faint millisecond pulsars (MSPs) and dark-matter (DM) annihilation. The dichotomy between these explanations is typically reflected by modeling them as two separate emission components. However, point-sources (PSs) such as MSPs become statistically degenerate with smooth Poi… ▽ More The two leading hypotheses for the Galactic Center Excess (GCE) in the $\textit{Fermi}$ data are an unresolved population of faint millisecond pulsars (MSPs) and dark-matter (DM) annihilation. The dichotomy between these explanations is typically reflected by modeling them as two separate emission components. However, point-sources (PSs) such as MSPs become statistically degenerate with smooth Poisson emission in the ultra-faint limit (formally where each source is expected to contribute much less than one photon on average), leading to an ambiguity that can render questions such as whether the emission is PS-like or Poissonian in nature ill-defined. We present a conceptually new approach that describes the PS and Poisson emission in a unified manner and only afterwards derives constraints on the Poissonian component from the so obtained results. For the implementation of this approach, we leverage deep learning techniques, centered around a neural network-based method for histogram regression that expresses uncertainties in terms of quantiles. We demonstrate that our method is robust against a number of systematics that have plagued previous approaches, in particular DM / PS misattribution. In the $\textit{Fermi}$ data, we find a faint GCE described by a median source-count distribution (SCD) peaked at a flux of $\sim4 \times 10^{-11} \ \text{counts} \ \text{cm}^{-2} \ \text{s}^{-1}$ (corresponding to $\sim3 - 4$ expected counts per PS), which would require $N \sim \mathcal{O}(10^4)$ sources to explain the entire excess (median value $N = \text{29,300}$ across the sky). Although faint, this SCD allows us to derive the constraint $η_P \leq 66\%$ for the Poissonian fraction of the GCE flux $η_P$ at 95% confidence, suggesting that a substantial amount of the GCE flux is due to PSs. △ Less

Submitted 15 December, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 36+9 pages, 15+7 figures, main results in Figs. 8 and 12. v2 matches published version in Phys. Rev. D

Journal ref: Phys. Rev. D 104 (2021), 123022

Velocity Dispersions of Clusters in the Dark Energy Survey Y3 redMaPPer Catalog

Authors: V. Wetzell, T. E. Jeltema, B. Hegland, S. Everett, P. A. Giles, R. Wilkinson, A. Farahi, M. Costanzi, D. L. Hollowood, E. Upsdell, A. Saro, J. Myles, A. Bermeo, S. Bhargava, C. A. Collins, D. Cross, O. Eiger, G. Gardner, M. Hilton, J. Jobel, P. Kelly, D. Laubner, A. R. Liddle, R. G. Mann, V. Martinez , et al. (74 additional authors not shown)

Abstract: We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $λ$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlation… ▽ More We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $λ$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlations between cluster velocity dispersion, richness, X-ray temperature and luminosity as well as central galaxy velocity offsets. The velocity dispersion-richness relation exhibits a bimodal distribution. The majority of clusters follow scaling relations between velocity dispersion, richness, and X-ray properties similar to those found for previous samples; however, there is a significant population of clusters with velocity dispersions which are high for their richness. These clusters account for roughly 22\% of the $λ< 70$ systems in our sample, but more than half (55\%) of $λ< 70$ clusters at $z>0.5$. A couple of these systems are hot and X-ray bright as expected for massive clusters with richnesses that appear to have been underestimated, but most appear to have high velocity dispersions for their X-ray properties likely due to line-of-sight structure. These results suggest that projection effects contribute significantly to redMaPPer selection, particularly at higher redshifts and lower richnesses. The redMaPPer determined richnesses for the velocity dispersion outliers are consistent with their X-ray properties, but several are X-ray undetected and deeper data is needed to understand their nature. △ Less

Submitted 9 June, 2022; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 22 pages, accepted to MNRAS

Dark Energy Survey Year 3 Results: Galaxy Sample for BAO Measurement

Authors: A. Carnero Rosell, M. Rodriguez-Monroy, M. Crocce, J. Elvin-Poole, A. Porredon, I. Ferrero, J. Mena-Fernandez, R. Cawthon, J. De Vicente, E. Gaztanaga, A. J. Ross, E. Sanchez, I. Sevilla-Noarbe, O. Alves, F. Andrade-Oliveira, J. Asorey, S. Avila, A. Brandao-Souza, H. Camacho, K. C. Chan, A. Ferte, J. Muir, W. Riquelme, R. Rosenfeld, D. Sanchez Cid , et al. (84 additional authors not shown)

Abstract: In this paper we present and validate the galaxy sample used for the analysis of the Baryon Acoustic Oscillation signal (BAO) in the Dark Energy Survey (DES) Y3 data. The definition is based on a colour and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.5, while ensuring a high quality photometric redshift determination. The sample covers $\approx 4100$ sq… ▽ More In this paper we present and validate the galaxy sample used for the analysis of the Baryon Acoustic Oscillation signal (BAO) in the Dark Energy Survey (DES) Y3 data. The definition is based on a colour and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.5, while ensuring a high quality photometric redshift determination. The sample covers $\approx 4100$ square degrees to a depth of $i = 22.3 \ (AB)$ at $10σ$. It contains 7,031,993 galaxies in the redshift range from $z$= 0.6 to 1.1, with a mean effective redshift of 0.835. Photometric redshifts are estimated with the machine learning algorithm DNF, and are validated using the VIPERS PDR2 sample. We find a mean redshift bias of $z_{\mathrm{bias}} \approx 0.01$ and a mean uncertainty, in units of $1+z$, of $σ_{68} \approx 0.03$. We evaluate the galaxy population of the sample, showing it is mostly built upon Elliptical to Sbc types. Furthermore, we find a low level of stellar contamination of $\lesssim 4\%$. We present the method used to mitigate the effect of spurious clustering coming from observing conditions and other large-scale systematics. We apply it to the DES Y3 BAO sample and calculate sample weights that are used to get a robust estimate of the galaxy clustering signal. This paper is one of a series dedicated to the analysis of the BAO signal in the DES Y3 data. In the companion papers, Ferrero et al. (2021) and DES Collaboration (2021), we present the galaxy mock catalogues used to calibrate the analysis and the angular diameter distance constraints obtained through the fitting to the BAO scale, respectively. The galaxy sample, masks and additional material will be released in the public DES data repository upon acceptance. △ Less

Submitted 21 October, 2021; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: Accepted for publication in MNRAS. 24 pages, 24 figures

Report number: FERMILAB-PUB-21-286-AE

Dark Energy Survey Year 3 Results: A 2.7% measurement of Baryon Acoustic Oscillation distance scale at redshift 0.835

Authors: DES Collaboration, T. M. C. Abbott, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Asorey, S. Avila, G. M. Bernstein, E. Bertin, A. Brandao-Souza, D. Brooks, D. L. Burke, J. Calcino, H. Camacho, A. Carnero Rosell, D. Carollo, M. Carrasco Kind, J. Carretero, F. J. Castander, R. Cawthon, K. C. Chan, A. Choi, C. Conselice, M. Costanzi, M. Crocce , et al. (86 additional authors not shown)

Abstract: We present angular diameter measurements obtained by measuring the position of Baryon Acoustic Oscillations (BAO) in an optimised sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg$^2$ with $0.6 < z_{\rm photo} < 1.1$ and a typical redshift uncertainty of $0.03(1+z)$. The sample selec… ▽ More We present angular diameter measurements obtained by measuring the position of Baryon Acoustic Oscillations (BAO) in an optimised sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg$^2$ with $0.6 < z_{\rm photo} < 1.1$ and a typical redshift uncertainty of $0.03(1+z)$. The sample selection is the same as in the BAO measurement with the first year of DES data, but the analysis presented here uses three times the area, extends to higher redshift and makes a number of improvements, including a fully analytical BAO template, the use of covariances from both theory and simulations, and an extensive pre-unblinding protocol. We used two different statistics: angular correlation function and power spectrum, and validate our pipeline with an ensemble of over 1500 realistic simulations. Both statistics yield compatible results. We combine the likelihoods derived from angular correlations and spherical harmonics to constrain the ratio of comoving angular diameter distance $D_M$ at the effective redshift of our sample to the sound horizon scale at the drag epoch. We obtain $D_M(z_{\rm eff}=0.835)/r_{\rm d} = 18.92 \pm 0.51$, which is consistent with, but smaller than, the Planck prediction assuming flat \lcdm, at the level of $2.3 σ$. The analysis was performed blind and is robust to changes in a number of analysis choices. It represents the most precise BAO distance measurement from imaging data to date, and is competitive with the latest transverse ones from spectroscopic samples at $z>0.75$. When combined with DES 3x2pt + SNIa, they lead to improvements in $H_0$ and $Ω_m$ constraints by $\sim 20\%$ △ Less

Submitted 18 March, 2022; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 24 pages, 10 figures, matches version accepted by PRD

Journal ref: Physical Review D, 2022, Volume 105, Issue 4, article id.043512

Kinematics of Antlia 2 and Crater 2 from The Southern Stellar Stream Spectroscopic Survey (S5)

Authors: Alexander P. Ji, Sergey E. Koposov, Ting S. Li, Denis Erkal, Andrew B. Pace, Joshua D. Simon, Vasily Belokurov, Lara R. Cullinane, Gary S. Da Costa, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Nora Shipp, Jeffrey D. Simpson, Daniel B. Zucker, Terese T. Hansen, Joss Bland-Hawthorn, S5 Collaboration

Abstract: We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia… ▽ More We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia EDR3 proper motions, detects a clear velocity gradient in Antlia 2 and a tentative velocity gradient in Crater 2. The velocity gradient magnitudes and directions are consistent with particle stream simulations of tidal disruption. Furthermore, the orbit and kinematics of Antlia 2 require a model that includes the reflex motion of the Milky Way induced by the Large Magellanic Cloud. We also find that Antlia 2's metallicity was previously overestimated, so it lies on the empirical luminosity-metallicity relation and is likely only now experiencing substantial stellar mass loss. Current dynamical models of Antlia 2 require it to have lost over 90% of its stars to tides, in tension with the low stellar mass loss implied by the updated metallicity. Overall, the new kinematic measurements support a tidal disruption scenario for the origin of these large and extended dwarf spheroidal galaxies. △ Less

Submitted 21 September, 2021; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: 23 pages (+appendix), 11 figures, accepted to ApJ. Updated to accepted version (1 new figure, minor changes to text) Video of Antlia 2 tidal disruption here: https://www.youtube.com/watch?v=wvt-Q4kRq90