Showing 1–50 of 93 results for author: Dawson, W

Astrometric Microlensing by Primordial Black Holes with The Roman Space Telescope

Authors: James Fardeen, Peter McGill, Scott E. Perkins, William A. Dawson, Natasha S. Abrams, Jessica R. Lu, Ming-Feng Ho, Simeon Bird

Abstract: Primordial Black Holes (PBHs) could explain some fraction of dark matter and shed light on many areas of early-universe physics. Despite over half a century of research interest, a PBH population has so far eluded detection. The most competitive constraints on the fraction of dark matter comprised of PBHs ($f_{\rm DM}$) in the $(10^{-9}-10)M_{\odot}$ mass-ranges come from photometric microlensing… ▽ More Primordial Black Holes (PBHs) could explain some fraction of dark matter and shed light on many areas of early-universe physics. Despite over half a century of research interest, a PBH population has so far eluded detection. The most competitive constraints on the fraction of dark matter comprised of PBHs ($f_{\rm DM}$) in the $(10^{-9}-10)M_{\odot}$ mass-ranges come from photometric microlensing and bound $f_{\rm DM}\lesssim10^{-2}-10^{-1}$. With the advent of the Roman Space Telescope with its sub-milliarcsecond (mas) astrometric capabilities and its planned Galactic Bulge Time Domain Survey (GBTDS), detecting astrometric microlensing signatures will become routine. Compared with photometric microlensing, astrometric microlensing signals are sensitive to different lens masses-distance configurations and contains different information, making it a complimentary lensing probe. At sub-mas astrometric precision, astrometric microlensing signals are typically detectable at larger lens-source separations than photometric signals, suggesting a microlensing detection channel of pure astrometric events. We use a Galactic simulation to predict the number of detectable microlensing events during the GBTDS via this pure astrometric microlensing channel. Assuming an absolute astrometric precision floor for bright stars of 0.1 mas for the GBTDS, we find that the number of detectable events peaks at $\approx 10^{3} f_{\rm DM}$ for a population of $ 1 M_{\odot}$ PBHs and tapers to $\approx 10f_{\rm DM}$ and $\approx 100f_{\rm DM}$ at $10^{-4}M_{\odot}$ and $10^{3}M_{\odot}$, respectively. Accounting for the distinguishability of PBHs from Stellar lenses, we conclude the GBTDS will be sensitive to a PBH population at $f_{\rm DM}$ down to $\approx10^{-1}-10^{-3}$ for $(10^{-1}-10^{2})M_{\odot}$ likely yielding novel PBH constraints. △ Less

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

Comments: 21 pages, 11 figures, accepted to AAS Journals

Spatially resolved microlensing timescale distributions across the Galactic bulge with the VVV survey

Authors: Zofia Kaczmarek, Peter McGill, N. Wyn Evans, Leigh C. Smith, Nathan Golovich, Eamonn Kerins, David Specht, William A. Dawson

Abstract: We analyze 1602 microlensing events found in the VISTA Variables in the Via Lactea (VVV) near-infrared (NIR) survey data. We obtain spatially-resolved, efficiency-corrected timescale distributions across the Galactic bulge ($|\ell|<10^\circ,$ $|b|<5^\circ$), using a Bayesian hierarchical model. Spatially-resolved peaks and means of the timescale distributions, along with their marginal distributio… ▽ More We analyze 1602 microlensing events found in the VISTA Variables in the Via Lactea (VVV) near-infrared (NIR) survey data. We obtain spatially-resolved, efficiency-corrected timescale distributions across the Galactic bulge ($|\ell|<10^\circ,$ $|b|<5^\circ$), using a Bayesian hierarchical model. Spatially-resolved peaks and means of the timescale distributions, along with their marginal distributions in strips of longitude and latitude, are in agreement at a 1$σ$ level with predictions based on the Besançon model of the Galaxy. We find that the event timescales in the central bulge fields ($|\ell| < 5^\circ$) are on average shorter than the non-central ($|\ell| > 5^\circ$) fields, with the average peak of the lognormal timescale distribution at 23.6 $\pm$ 1.9 days for the central fields and 29.0 $\pm$ 3.0 days for the non-central fields. Our ability to probe the structure of the Bulge with this sample of NIR microlensing events is limited by the VVV survey's sparse cadence and relatively small number of detected microlensing events compared to dedicated optical surveys. Looking forward to future surveys, we investigate the capability of the Roman telescope to detect spatially-resolved asymmetries in the timescale distributions. We propose two pairs of Roman fields, centred on ($\ell = \pm 9,5^\circ$, $b=-0.125^\circ$) and ($\ell = -5^\circ$, $b=\pm 1.375^\circ$) as good targets to measure the asymmetry in longitude and latitude, respectively. △ Less

Submitted 11 February, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

Comments: 13 pages, 11 figures, accepted for publication in MNRAS

Exploratory Data Science on Supercomputers for Quantum Mechanical Calculations

Authors: William Dawson, Louis Beal, Laura E. Ratcliff, Martina Stella, Takahito Nakajima, Luigi Genovese

Abstract: Literate programming - the bringing together of program code and natural language narratives - has become a ubiquitous approach in the realm of data science. This methodology is appealing as well for the domain of Density Functional Theory (DFT) calculations, particularly for interactively developing new methodologies and workflows. However, effective use of literate programming is hampered by old… ▽ More Literate programming - the bringing together of program code and natural language narratives - has become a ubiquitous approach in the realm of data science. This methodology is appealing as well for the domain of Density Functional Theory (DFT) calculations, particularly for interactively developing new methodologies and workflows. However, effective use of literate programming is hampered by old programming paradigms and the difficulties associated with using High Performance Computing (HPC) resources. Here we present two Python libraries that aim to remove these hurdles. First, we describe the PyBigDFT library, which can be used to setup materials or molecular systems and provides high-level access to the wavelet based BigDFT code. We then present the related remotemanager library, which is able to serialize and execute arbitrary Python functions on remote supercomputers. We show how together these libraries enable transparent access to HPC based DFT calculations and can serve as building blocks for rapid prototyping and data exploration. △ Less

Submitted 14 October, 2023; originally announced October 2023.

Disentangling the Black Hole Mass Spectrum with Photometric Microlensing Surveys

Authors: Scott Ellis Perkins, Peter McGill, William Dawson, Natasha S. Abrams, Casey Y. Lam, Ming-Feng Ho, Jessica R. Lu, Simeon Bird, Kerianne Pruett, Nathan Golovich, George Chapline

Abstract: From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However,… ▽ More From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However, photometric microlensing observables are degenerate for different types of lenses, and one can seldom classify an event as involving either a compact object or stellar lens on its own. To address this difficulty, we apply a Bayesian framework that treats lens type probabilistically and jointly with a lens population model. This method allows lens population characteristics to be inferred despite intrinsic uncertainty in the lens-class of any single event. We investigate this method's effectiveness on a simulated ground-based photometric survey in the context of characterizing a hypothetical population of primordial black holes (PBHs) with an average mass of $30 M_{\odot}$. On simulated data, our method outperforms current black hole (BH) lens identification pipelines and characterizes different subpopulations of lenses while jointly constraining the PBH contribution to dark matter to ${\approx}25$\%. Key to robust inference, our method can marginalize over population model uncertainty. We find the lower mass cutoff for stellar origin BHs, a key observable in understanding the BH mass gap, particularly difficult to infer in our simulations. This work lays the foundation for cutting-edge PBH abundance constraints to be extracted from current photometric microlensing surveys. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 31 pages, 18 figures, submitted to AAS

Roman CCS White Paper: Characterizing the Galactic population of isolated black holes

Authors: Casey Y. Lam, Natasha Abrams, Jeff Andrews, Etienne Bachelet, Arash Bahramian, David Bennett, Valerio Bozza, Floor Broekgaarden, Sukanya Chakrabarti, William Dawson, Kareem El-Badry, Maya Fishbach, Giacomo Fragione, Scott Gaudi, Abhimat Gautam, Ryosuke Hirai, Daniel Holz, Matthew Hosek Jr., Macy Huston, Tharindu Jayasinghe, Samson Johnson, Daisuke Kawata, Naoki Koshimoto, Jessica R. Lu, Ilya Mandel , et al. (12 additional authors not shown)

Abstract: Although there are estimated to be 100 million isolated black holes (BHs) in the Milky Way, only one has been found so far, resulting in significant uncertainty about their properties. The Galactic Bulge Time Domain Survey provides the only opportunity in the coming decades to grow this catalog by order(s) of magnitude. This can be achieved if 1) Roman's astrometric potential is fully realized in… ▽ More Although there are estimated to be 100 million isolated black holes (BHs) in the Milky Way, only one has been found so far, resulting in significant uncertainty about their properties. The Galactic Bulge Time Domain Survey provides the only opportunity in the coming decades to grow this catalog by order(s) of magnitude. This can be achieved if 1) Roman's astrometric potential is fully realized in the observation strategy and software pipelines, 2) Roman's observational gaps of the Bulge are minimized, and 3) observations with ground-based facilities are taken of the Bulge to fill in gaps during non-Bulge seasons. A large sample of isolated BHs will enable a broad range of astrophysical questions to be answered, such as massive stellar evolution, origin of gravitational wave sources, supernova physics, and the growth of supermassive BHs, maximizing Roman's scientific return. △ Less

Submitted 29 June, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 20 pages. Submitted in response to Nancy Grace Roman Space Telescope white paper call: https://roman.gsfc.nasa.gov/science/ccs_white_papers.html. v2 fixes a typo in Figure 5 axis label (days --> year)

Complexity Reduction in Density Functional Theory: Locality in Space and Energy

Authors: William Dawson, Eisuke Kawashima, Laura E. Ratcliff, Muneaki Kamiya, Luigi Genovese, Takahito Nakajima

Abstract: We present recent developments of the NTChem program for performing large scale hybrid Density Functional Theory calculations on the supercomputer Fugaku. We combine these developments with our recently proposed Complexity Reduction Framework to assess the impact of basis set and functional choice on its measures of fragment quality and interaction. We further exploit the all electron representati… ▽ More We present recent developments of the NTChem program for performing large scale hybrid Density Functional Theory calculations on the supercomputer Fugaku. We combine these developments with our recently proposed Complexity Reduction Framework to assess the impact of basis set and functional choice on its measures of fragment quality and interaction. We further exploit the all electron representation to study system fragmentation in various energy envelopes. Building off this analysis, we propose two algorithms for computing the orbital energies of the Kohn-Sham Hamiltonian. We demonstrate these algorithms can efficiently be applied to systems composed of thousands of atoms and as an analysis tool that reveals the origin of spectral properties. △ Less

Submitted 28 April, 2023; v1 submitted 16 January, 2023; originally announced January 2023.

Comments: Accepted Manuscript

Journal ref: J. Chem. Phys. 158, 164114 (2023)

Primordial Black Hole Dark Matter Simulations Using PopSyCLE

Authors: K. Pruett, W. Dawson, M. S. Medford, C. Lam, J. R. Lu, S. Perkins, P. McGill, N. Golovich, G. Chapline

Abstract: Primordial black holes (PBHs), theorized to have originated in the early universe, are speculated to be a viable form of dark matter. If they exist, they should be detectable through photometric and astrometric signals resulting from gravitational microlensing of stars in the Milky Way. Population Synthesis for Compact-object Lensing Events, or PopSyCLE, is a simulation code that enables users to… ▽ More Primordial black holes (PBHs), theorized to have originated in the early universe, are speculated to be a viable form of dark matter. If they exist, they should be detectable through photometric and astrometric signals resulting from gravitational microlensing of stars in the Milky Way. Population Synthesis for Compact-object Lensing Events, or PopSyCLE, is a simulation code that enables users to simulate microlensing surveys, and is the first of its kind to include both photometric and astrometric microlensing effects, which are important for potential PBH detection and characterization. To estimate the number of observable PBH microlensing events we modify PopSyCLE to include a dark matter halo consisting of PBHs. We detail our PBH population model, and demonstrate our PopSyCLE + PBH results through simulations of the OGLE-IV and Roman microlensing surveys. We provide a proof-of-concept analysis for adding PBHs into PopSyCLE, and thus include many simplifying assumptions, such as $f_{\text{DM}}$, the fraction of dark matter composed of PBHs, and $\bar{m}_{\text{PBH}}$, mean PBH mass. Assuming $\bar{m}_{\text{PBH}}=30$ $M_{\odot}$, we find $\sim$ 3.6$f_{\text{DM}}$ times as many PBH microlensing events than stellar evolved black hole events, a PBH average peak Einstein crossing time of $\sim$ 91.5 days, estimate on order of $10^2f_{\text{DM}}$ PBH events within the 8 year OGLE-IV results, and estimate Roman to detect $\sim$ 1,000$f_{\text{DM}}$ PBH microlensing events throughout its planned microlensing survey. △ Less

Submitted 8 April, 2024; v1 submitted 12 November, 2022; originally announced November 2022.

Roadmap on Electronic Structure Codes in the Exascale Era

Authors: Vikram Gavini, Stefano Baroni, Volker Blum, David R. Bowler, Alexander Buccheri, James R. Chelikowsky, Sambit Das, William Dawson, Pietro Delugas, Mehmet Dogan, Claudia Draxl, Giulia Galli, Luigi Genovese, Paolo Giannozzi, Matteo Giantomassi, Xavier Gonze, Marco Govoni, Andris Gulans, François Gygi, John M. Herbert, Sebastian Kokott, Thomas D. Kühne, Kai-Hsin Liou, Tsuyoshi Miyazaki, Phani Motamarri , et al. (16 additional authors not shown)

Abstract: Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources d… ▽ More Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Submitted as a roadmap article to Modelling and Simulation in Materials Science and Engineering; Address any correspondence to Vikram Gavini (vikramg@umich.edu) and Danny Perez (danny_perez@lanl.gov)

Report of the Topical Group on Cosmic Probes of Dark Matter for Snowmass 2021

Authors: Alex Drlica-Wagner, Chanda Prescod-Weinstein, Hai-Bo Yu, Andrea Albert, Mustafa Amin, Arka Banerjee, Masha Baryakhtar, Keith Bechtol, Simeon Bird, Simon Birrer, Torsten Bringmann, Regina Caputo, Sukanya Chakrabarti, Thomas Y. Chen, Djuna Croon, Francis-Yan Cyr-Racine, William A. Dawson, Cora Dvorkin, Vera Gluscevic, Daniel Gilman, Daniel Grin, Renée Hložek, Rebecca K. Leane, Ting S. Li, Yao-Yuan Mao , et al. (15 additional authors not shown)

Abstract: Cosmological and astrophysical observations currently provide the only robust, positive evidence for dark matter. Cosmic probes of dark matter, which seek to determine the fundamental properties of dark matter through observations of the cosmos, have emerged as a promising means to reveal the nature of dark matter. This report summarizes the current status and future potential of cosmic probes to… ▽ More Cosmological and astrophysical observations currently provide the only robust, positive evidence for dark matter. Cosmic probes of dark matter, which seek to determine the fundamental properties of dark matter through observations of the cosmos, have emerged as a promising means to reveal the nature of dark matter. This report summarizes the current status and future potential of cosmic probes to inform our understanding of the fundamental nature of dark matter in the coming decade. △ Less

Submitted 13 December, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: Report of the CF3 Topical Group for Snowmass 2021; 35 pages, 10 figures, many references. V3 updates Fig 3-2 and the author list

Probing Disorder in 2CzPN using Core and Valence States

Authors: Nathalie K. Fernando, Martina Stella, William Dawson, Takahito Nakajima, Luigi Genovese, Anna Regoutz, Laura E. Ratcliff

Abstract: Molecules which exhibit thermally activated delayed fluorescence (TADF) show great promise for use in efficient, environmentally-friendly OLEDs, and thus the design of new TADF emitters is an active area of research. However, when used in devices, they are typically in the form of disordered thin films, where both the external molecular environment and thermally-induced internal variations in para… ▽ More Molecules which exhibit thermally activated delayed fluorescence (TADF) show great promise for use in efficient, environmentally-friendly OLEDs, and thus the design of new TADF emitters is an active area of research. However, when used in devices, they are typically in the form of disordered thin films, where both the external molecular environment and thermally-induced internal variations in parameters such as the torsion angle can strongly influence their electronic structure. In this work, we use density functional theory and X-ray photoelectron spectroscopy to investigate the impact of disorder on both core and valence states in the TADF emitter 2CzPN. By simulating gas phase molecules displaying varying levels of disorder, we assess the relative sensitivity of the different states to factors such as varying torsion angle. The theoretical results for both core and valence states show good agreement with experiment, thereby also highlighting the advantages of our approach for interpreting experimental spectra of large aromatic molecules, which are too complex to interpret based solely on experimental data. △ Less

Submitted 10 June, 2022; originally announced June 2022.

Snowmass2021 Cosmic Frontier White Paper:Primordial Black Hole Dark Matter

Authors: Simeon Bird, Andrea Albert, Will Dawson, Yacine Ali-Haimoud, Adam Coogan, Alex Drlica-Wagner, Qi Feng, Derek Inman, Keisuke Inomata, Ely Kovetz, Alexander Kusenko, Benjamin V. Lehmann, Julian B. Munoz, Rajeev Singh, Volodymyr Takhistov, Yu-Dai Tsai

Abstract: Primordial Black Holes (PBHs) are a viable candidate to comprise some or all of the dark matter and provide a unique window into the high-energy physics of the early universe. This white paper discusses the scientific motivation, current status, and future reach of observational searches for PBHs. Future observational facilities supported by DOE, NSF, and NASA will provide unprecedented sensitivit… ▽ More Primordial Black Holes (PBHs) are a viable candidate to comprise some or all of the dark matter and provide a unique window into the high-energy physics of the early universe. This white paper discusses the scientific motivation, current status, and future reach of observational searches for PBHs. Future observational facilities supported by DOE, NSF, and NASA will provide unprecedented sensitivity to PBHs. However, devoted analysis pipelines and theoretical modeling are required to fully leverage these novel data. The search for PBHs constitutes a low-cost, high-reward science case with significant impact on the high energy physics community. △ Less

Submitted 1 July, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 22 pages, 4 figures, submission to the Snowmass 2021 process. Update to add references

Positivity Preserving Density Matrix Minimization at Finite Temperatures via Square Root

Authors: Jacob M. Leamer, William Dawson, Denys I. Bondar

Abstract: We present a Wave Operator Minimization (WOM) method for calculating the Fermi-Dirac density matrix for electronic structure problems at finite temperature while preserving physicality by construction using the wave operator, i.e., the square root of the density matrix. WOM models cooling a state initially at infinite temperature down to the desired finite temperature. We consider both the grand c… ▽ More We present a Wave Operator Minimization (WOM) method for calculating the Fermi-Dirac density matrix for electronic structure problems at finite temperature while preserving physicality by construction using the wave operator, i.e., the square root of the density matrix. WOM models cooling a state initially at infinite temperature down to the desired finite temperature. We consider both the grand canonical (constant chemical potential) and canonical (constant number of electrons) ensembles. Additionally, we show that the number of steps required for convergence is independent of the number of atoms in the system. We hope that the discussion and results presented in this article reinvigorates interest in density matrix minimization methods. △ Less

Submitted 6 December, 2023; v1 submitted 11 March, 2021; originally announced March 2021.

Comments: 10 pages, 5 figures, 1 algorithm, supplementary information 5 pages, 3 figures

Journal ref: J Chem Phys 160 (2024)

Hybrid Information-driven Multi-agent Reinforcement Learning

Authors: William A. Dawson, Ruben Glatt, Edward Rusu, Braden C. Soper, Ryan A. Goldhahn

Abstract: Information theoretic sensor management approaches are an ideal solution to state estimation problems when considering the optimal control of multi-agent systems, however they are too computationally intensive for large state spaces, especially when considering the limited computational resources typical of large-scale distributed multi-agent systems. Reinforcement learning (RL) is a promising alt… ▽ More Information theoretic sensor management approaches are an ideal solution to state estimation problems when considering the optimal control of multi-agent systems, however they are too computationally intensive for large state spaces, especially when considering the limited computational resources typical of large-scale distributed multi-agent systems. Reinforcement learning (RL) is a promising alternative which can find approximate solutions to distributed optimal control problems that take into account the resource constraints inherent in many systems of distributed agents. However, the RL training can be prohibitively inefficient, especially in low-information environments where agents receive little to no feedback in large portions of the state space. We propose a hybrid information-driven multi-agent reinforcement learning (MARL) approach that utilizes information theoretic models as heuristics to help the agents navigate large sparse state spaces, coupled with information based rewards in an RL framework to learn higher-level policies. This paper presents our ongoing work towards this objective. Our preliminary findings show that such an approach can result in a system of agents that are approximately three orders of magnitude more efficient at exploring a sparse state space than naive baseline metrics. While the work is still in its early stages, it provides a promising direction for future research. △ Less

Submitted 1 February, 2021; originally announced February 2021.

Comments: Published at Workshop on Challenges and Opportunities for Multi-Agent Reinforcement Learning (COMARL AAAI 2021). This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under contract DE-AC52-07NA27344. Lawrence Livermore National Security, LLC through the support of LDRD 20-SI-005. LLNL-CONF-816423

Exemplary Merging Clusters: Weak-lensing and X-ray Analysis of the Double Radio Relic Merging Galaxy Clusters MACS 1752.0+4440 and ZWCL 1856.8+6616

Authors: Kyle Finner, Kim HyeongHan, M. James Jee, David Wittman, William R. Forman, Reinout J. van Weeren, Nathan R. Golovich, William A. Dawson, Alexander Jones, Francesco de Gasperin, Christine Jones

Abstract: The investigation of merging galaxy clusters that exhibit radio relics is strengthening our understanding of the formation and evolution of galaxy clusters, the nature of dark matter, the intracluster medium, and astrophysical particle acceleration. Each merging cluster provides only a single view of the cluster formation process and the variety of merging clusters is vast. Clusters hosting double… ▽ More The investigation of merging galaxy clusters that exhibit radio relics is strengthening our understanding of the formation and evolution of galaxy clusters, the nature of dark matter, the intracluster medium, and astrophysical particle acceleration. Each merging cluster provides only a single view of the cluster formation process and the variety of merging clusters is vast. Clusters hosting double radio relics are rare and extremely important because they allow tight constraints on the merger scenario. We present a weak-lensing and X-ray analysis of MACSJ1752.0+4440 ($z$=0.365) and ZWCL1856.8+6616 ($z$=0.304), two double radio relic clusters. Our weak-lensing mass estimates show that each cluster is a major merger with approximately 1:1 mass ratio. The total mass of MACSJ1752.0+4440 (ZWCL1856.8+6616) is $M_{200}=14.7^{+3.8}_{-3.3}\times10^{14}\ $M$_\odot$ ($M_{200}=2.4^{+0.9}_{-0.7}\times10^{14}\ $M$_\odot$). We find that these two clusters have comparable features in their weak-lensing and gas distributions, even though the systems have vastly different total masses. From the likeness of the X-ray morphologies and the remarkable symmetry of the radio relics, we propose that both systems underwent nearly head-on collisions. However, revelations from the hot-gas features and our multiwavelength data analysis suggest that ZWCL1856.8+6618 is likely at a later merger phase than MACSJ1752.0+4440. We postulate that the SW radio relic in MACSJ1752.0+4440 is a result of particle re-acceleration. △ Less

Submitted 11 October, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: 21 pages, 12 figures, Submitted to ApJ

A Reanalysis of Public Galactic Bulge Gravitational Microlensing Events from OGLE-III and IV

Authors: Nathan Golovich, William A. Dawson, Fran Bartolić, Casey Y. Lam, Jessica R. Lu, Michael S. Medford, Michael D. Schneider, George Chapline, Edward F. Schlafly, Alex Drlica-Wagner, Kerianne Pruett

Abstract: Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. Surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light… ▽ More Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. Surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light curves as well as over-simplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian Process model. We used light curves for nearly 10,000 OGLE-III and IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the timescale distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and previous analyses over-estimated the number of long timescale ($t_E>100$ days) events due to their over simplistic models ignoring parallax effects and stellar variability. We use our fits to identify hundreds of events that are likely black holes. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: Submitted version, in review, 33 pages, 18 figures, MCMC posterior samples available by publisher after acceptance

RELICS: Properties of z>5.5 Galaxies Inferred from Spitzer and Hubble Imaging Including A Candidate z~6.8 Strong [OIII] Emitter

Authors: Victoria Strait, Marusa Bradac, Dan Coe, Brian C. Lemaux, Adam Carnall, Larry Bradley, Debora Pelliccia, Keren Sharon, Adi Zitrin, Ana Acebron, Chloe Neufeld, Felipe Andrade-Santos, Roberto J. Avila, Brenda L. Frye, Guillaume Mahler, Mario Nonino, Sara Ogaz, Masamune Oguri, Masami Ouchi, Rachel Paterno-Mahler, Daniel P. Stark, Ramesh Mainali, Pascal A. Oesch, Michele Trenti, Daniela Carrasco , et al. (4 additional authors not shown)

Abstract: We present constraints on the physical properties (including stellar mass, age, and star formation rate) of 207 $6\lesssim z \lesssim8$ galaxy candidates from the Reionization Lensing Cluster Survey (RELICS) and companion Spitzer-RELICS surveys. We measure photometry using T-PHOT and perform spectral energy distribution fitting using EA$z$Y and BAGPIPES. Of the 207 candidates for which we could su… ▽ More We present constraints on the physical properties (including stellar mass, age, and star formation rate) of 207 $6\lesssim z \lesssim8$ galaxy candidates from the Reionization Lensing Cluster Survey (RELICS) and companion Spitzer-RELICS surveys. We measure photometry using T-PHOT and perform spectral energy distribution fitting using EA$z$Y and BAGPIPES. Of the 207 candidates for which we could successfully measure Spitzer fluxes, 23 were demoted to likely low redshift ($z<4$). Among the remaining high redshift candidates, we find intrinsic stellar masses between $1\times10^6\rm{M_{\odot}}$ and $4\times10^9\rm{M_\odot}$, and rest-frame UV absolute magnitudes between $-22.6$ and $-14.5$ mag. While our sample is mostly comprised of $L_{UV}/L^*_{UV}<1$ galaxies, there are a number of brighter objects in the sample, extending to $L_{UV}/L^*_{UV}\sim2$. The galaxies in our sample span approximately four orders of magnitude in stellar mass and star-formation rates, and exhibit ages ranging from maximally young to maximally old. We highlight 11 galaxies which have detections in Spitzer/IRAC imaging and redshift estimates $z\geq6.5$, several of which show evidence for some combination of evolved stellar populations, large contributions of nebular emission lines, and/or dust. Among these is PLCKG287+32-2013, one of the brightest $z\sim7$ candidates known (AB mag 24.9) with a Spitzer 3.6$μ$m flux excess suggesting strong [OIII] + H-$β$ emission ($\sim$1000Å rest-frame equivalent width). We discuss the possible uses and limits of our sample and present a public catalog of Hubble 0.4--1.6$μ$m + Spitzer 3.6$μ$m and 4.5$μ$m photometry along with physical property estimates for all 207 objects in the sample. Because of their apparent brightnesses, high redshifts, and variety of stellar populations, these objects are excellent targets for follow-up with James Webb Space Telescope. △ Less

Submitted 3 February, 2021; v1 submitted 31 August, 2020; originally announced September 2020.

Comments: 20 pages, 9 figures

Gravitational Microlensing Event Statistics for the Zwicky Transient Facility

Authors: Michael S. Medford, Jessica R. Lu, William A. Dawson, Casey Y. Lam, Nathan R. Golovich, Edward F. Schlafly, Peter Nugent

Abstract: Microlensing surveys have discovered thousands of events with almost all events discovered within the Galactic bulge or toward the Magellanic clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky every few nights including the Galactic plane. ZTF observes $\sim10^9$ stars in g-band and r… ▽ More Microlensing surveys have discovered thousands of events with almost all events discovered within the Galactic bulge or toward the Magellanic clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky every few nights including the Galactic plane. ZTF observes $\sim10^9$ stars in g-band and r-band and can significantly contribute to the observed microlensing population. We predict that ZTF will observe $\sim$1100 microlensing events in three years of observing within $10^\circ$ degrees latitude of the Galactic plane, with $\sim$500 events in the outer Galaxy ($\ell \geq 10^\circ$). This yield increases to $\sim$1400 ($\sim$800) events by combining every three ZTF exposures, $\sim$1800 ($\sim$900) events if ZTF observes for a total of five years, and $\sim$2400 ($\sim$1300) events for a five year survey with post-processing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: 19 pages, 13 figures, 5 tables, accepted to ApJ (6/4/2020), microlensing simulation catalogs available at https://portal.nersc.gov/project/uLens/Galactic_Microlensing_Distributions

Journal ref: The Astrophysical Journal, 2020, Volume 897, Number 2, 144

ELSI -- An Open Infrastructure for Electronic Structure Solvers

Authors: Victor Wen-zhe Yu, Carmen Campos, William Dawson, Alberto García, Ville Havu, Ben Hourahine, William P Huhn, Mathias Jacquelin, Weile Jia, Murat Keçeli, Raul Laasner, Yingzhou Li, Lin Lin, Jianfeng Lu, Jonathan Moussa, Jose E Roman, Álvaro Vázquez-Mayagoitia, Chao Yang, Volker Blum

Abstract: Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infi… ▽ More Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn-Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. Finally, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets. △ Less

Submitted 4 July, 2020; v1 submitted 31 December, 2019; originally announced December 2019.

Journal ref: Computer Physics Communications 256 (2020) 107459

PopSyCLE: A New Population Synthesis Code for Compact Object Microlensing Events

Authors: Casey Y. Lam, Jessica R. Lu, Matthew W. Hosek, Jr., William A. Dawson, Nathan R. Golovich

Abstract: We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating th… ▽ More We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating the optimal way to find black holes (BHs) with microlensing. Candidate BHs have typically been selected from wide-field photometric microlensing surveys, such as OGLE, by selecting events with long Einstein crossing times ($t_E>120$ days). These events can be selected at closest approach and monitored astrometrically in order to constrain the mass of each lens; PopSyCLE predicts a BH detection rate of ~40% for such a program. We find that the detection rate can be enhanced to ~85% by selecting events with both $t_E>120$ days and a microlensing parallax of $π_E<0.08$. Unfortunately, such a selection criterion cannot be applied during the event as $π_E$ requires both pre- and post-peak photometry. However, historical microlensing events from photometric surveys can be revisited using this new selection criteria in order to statistically constrain the abundance of BHs in the Milky Way. The future WFIRST microlensing survey provides both precise photometry and astrometry and will yield individual masses of $\mathcal{O}(100-1000) black holes, which is at least an order of magnitude more than is possible with individual candidate follow-up with current facilities. The resulting sample of BH masses from WFIRST will begin to constrain the shape of the black hole present-day mass function, BH multiplicity, and BH kick velocity distributions. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: 33 pages, 22 figures, accepted to ApJ

RELICS: A Very Large ($θ_{E}\sim40"$) Cluster Lens -- RXC J0032.1+1808

Authors: Ana Acebron, Adi Zitrin, Dan Coe, Guillaume Mahler, Keren Sharon, Masamune Oguri, Maruša Bradač, Larry Bradley, Brenda Frye, Christine J. Forman, Victoria Strait, Yuanyuan Su, Keiichi Umetsu, Felipe Andrade-Santos, Roberto J. Avila, Daniela Carrasco, Catherine Cerny, Nicole G. Czakon, William A. Dawson, Carter Fox, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Shotaro Kikuchihara, Daniel Lam , et al. (16 additional authors not shown)

Abstract: Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $θ_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $θ_{E}\sim40\arcsec$ or above (for… ▽ More Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $θ_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $θ_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $θ_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties. △ Less

Submitted 11 May, 2020; v1 submitted 5 December, 2019; originally announced December 2019.

Comments: 23 pages, accepted for publication in ApJ

Complexity Reduction in Density Functional Theory Calculations of Large Systems: System Partitioning and Fragment Embedding

Authors: William Dawson, Stephan Mohr, Laura E. Ratcliff, Takahito Nakajima, Luigi Genovese

Abstract: With the development of low order scaling methods for performing Kohn-Sham Density Functional Theory, it is now possible to perform fully quantum mechanical calculations of systems containing tens of thousands of atoms. However, with an increase in the size of system treated comes an increase in complexity, making it challenging to analyze such large systems and determine the cause of emergent pro… ▽ More With the development of low order scaling methods for performing Kohn-Sham Density Functional Theory, it is now possible to perform fully quantum mechanical calculations of systems containing tens of thousands of atoms. However, with an increase in the size of system treated comes an increase in complexity, making it challenging to analyze such large systems and determine the cause of emergent properties. To address this issue, in this paper we present a systematic complexity reduction methodology which can break down large systems into their constituent fragments, and quantify inter-fragment interactions. The methodology proposed here requires no a priori information or user interaction, allowing a single workflow to be automatically applied to any system of interest. We apply this approach to a variety of different systems, and show how it allows for the derivation of new system descriptors, the design of QM/MM partitioning schemes, and the novel application of graph metrics to molecules and materials. △ Less

Submitted 30 March, 2020; v1 submitted 19 November, 2019; originally announced November 2019.

Chandra Observations of the Spectacular A3411-12 Merger Event

Authors: Felipe Andrade-Santos, Reinout J. van Weeren, Gabriella Di Gennaro, David Wittman, Dongsu Ryu, Dharam Vir Lal, Vinicius M. Placco, Kevin Fogarty, M. James Jee, Andra Stroe, David Sobral, William R. Forman, Christine Jones, Ralph P. Kraft, Stephen S. Murray, Marcus Brüggen, Hyesung Kang, Rafael Santucci, Nathan Golovich, William Dawson

Abstract: We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The… ▽ More We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.14 \pm 0.01$. We compute the shock strength using density jumps to conclude that the Mach number of the merging subcluster is small ($M \leq 1.15_{-0.09}^{+0.14}$). We also present pseudo-density, projected temperature, pseudo-pressure, and pseudo-entropy maps. Based on the pseudo-entropy map we conclude that the cluster is undergoing a mild merger, consistent with the small Mach number. On the other hand, radio relics extend over Mpc scale in the A3411-12 system, which strongly suggests that a population of energetic electrons already existed over extended regions of the cluster. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 23 pages, 17 figures, 4 tables. Accepted for publication in ApJ

Joint Survey Processing of LSST, Euclid and WFIRST: Enabling a broad array of astrophysics and cosmology through pixel level combinations of datasets

Authors: R. Chary, G. Brammer, P. Capak, W. Dawson, A. Faisst, S. Fajardo-Acosta, H. C. Ferguson, C. J. Grillmair, S. Hemmati, A. Koekemoer, B. Lee, R. Lupton, S. Malhotra, P. Melchior, I. Momcheva, J. Newman, J. Masiero, R. Paladini, A. Prakash, J. Rhodes, B. Rusholme, M. Schneider, N. Stickley, A. Smith, W. M. Wood-Vasey , et al. (1 additional authors not shown)

Abstract: Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper… ▽ More Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: Astro2020 APC White Paper, 11 pages

Astro2020: From Stars to Compact Objects: The Initial-Final Mass Relation

Authors: Jessica R. Lu, Casey Lam, Will Dawson, B. Scott Gaudi, Nathan Golovich, Michael Medford, Fatima Abdurrahman, Rachael L. Beaton

Abstract: One of the key phases of stellar evolution that remains poorly understood is stellar death. We lack a predictive model for how a star of a given mass explodes and what kind of remnant it leaves behind (i.e. the initial-final mass relation, IFMR). Progress has been limited due to the difficulty in finding and weighing black holes and neutron stars in large numbers. Technological advances that allow… ▽ More One of the key phases of stellar evolution that remains poorly understood is stellar death. We lack a predictive model for how a star of a given mass explodes and what kind of remnant it leaves behind (i.e. the initial-final mass relation, IFMR). Progress has been limited due to the difficulty in finding and weighing black holes and neutron stars in large numbers. Technological advances that allow for sub-milliarcsecond astrometry in crowded fields have opened a new window for finding black holes and neutron stars: astrometric gravitational lensing. Finding and weighing a sample of compact objects with astrometric microlensing will allow us to place some of the first constraints on the present-day mass function of isolated black holes and neutron stars, their multiplicity, and their kick velocities. All of these are fundamental inputs into understanding the death phase of stellar evolution, improving supernovae models, and interpreting LIGO detections in an astrophysical context. To achieve these goals, we require large area surveys, such as the WFIRST exoplanet microlensing survey, to photometrically identify long-duration (>120 day), un-blended microlensing events as candidate compact objects. We also require high-precision astrometric follow-up monitoring using extremely large telescopes, equipped with adaptive optics, such as TMT and GMT. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: Astro2020 Science White Paper

Primordial Black Hole Microlensing: The Einstein Crossing Time Distribution

Authors: Jessica R. Lu, Casey Y. Lam, Michael Medford, William Dawson, Nathan Golovich

Abstract: Gravitational microlensing is one of the few means of finding primordial black holes (PBHs), if they exist. Recent LIGO detections of 30 Msun black holes have re-invigorated the search for PBHs in the 10-100 Msun mass regime. Unfortunately, individual PBH microlensing events cannot easily be distinguished from stellar lensing events from photometry alone. However, the distribution of microlensing… ▽ More Gravitational microlensing is one of the few means of finding primordial black holes (PBHs), if they exist. Recent LIGO detections of 30 Msun black holes have re-invigorated the search for PBHs in the 10-100 Msun mass regime. Unfortunately, individual PBH microlensing events cannot easily be distinguished from stellar lensing events from photometry alone. However, the distribution of microlensing timescales (tE, the Einstein radius crossing time) can be analyzed in a statistical sense using models of the Milky Way with and without PBHs. While previous works have presented both theoretical models and observational constrains for PBHs (e.g. Calcino et al. 2018; Niikura et al. 2019), surprisingly, they rarely show the observed quantity -- the tE distribution -- for different abundances of PBHs relative to the total dark matter mass (fPBH). We present a simple calculation of how the tE distribution changes between models with and without PBHs. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: RNAAS

Journal ref: 2019, RNAAS, Vol 3, 58

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

RELICS: Reionization Lensing Cluster Survey

Authors: Dan Coe, Brett Salmon, Marusa Bradac, Larry D. Bradley, Keren Sharon, Adi Zitrin, Ana Acebron, Catherine Cerny, Nathalia Cibirka, Victoria Strait, Rachel Paterno-Mahler, Guillaume Mahler, Roberto J. Avila, Sara Ogaz, Kuang-Han Huang, Debora Pelliccia, Daniel P. Stark, Ramesh Mainali, Pascal A. Oesch, Michele Trenti, Daniela Carrasco, William A. Dawson, Steven A. Rodney, Louis-Gregory Strolger, Adam G. Riess , et al. (32 additional authors not shown)

Abstract: Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed p… ▽ More Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals. △ Less

Submitted 5 March, 2019; originally announced March 2019.

Comments: 29 pages, 6 figures, submitted to ApJ. For reduced images, catalogs, lens models, and more, see relics.stsci.edu

Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope

Authors: Alex Drlica-Wagner, Yao-Yuan Mao, Susmita Adhikari, Robert Armstrong, Arka Banerjee, Nilanjan Banik, Keith Bechtol, Simeon Bird, Kimberly K. Boddy, Ana Bonaca, Jo Bovy, Matthew R. Buckley, Esra Bulbul, Chihway Chang, George Chapline, Johann Cohen-Tanugi, Alessandro Cuoco, Francis-Yan Cyr-Racine, William A. Dawson, Ana Díaz Rivero, Cora Dvorkin, Denis Erkal, Christopher D. Fassnacht, Juan García-Bellido, Maurizio Giannotti , et al. (75 additional authors not shown)

Abstract: Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We d… ▽ More Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ . △ Less

Submitted 24 April, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 96 pages, 22 figures, 1 table

Report number: FERMILAB-PUB-19-048-A-AE

Evidence for a merger induced shock wave in ZwCl\,0008.8+5215 with {\it Chandra} and {\it Suzaku}

Authors: G. Di Gennaro, R. J. van Weeren, F. Andrade-Santos, H. Akamatsu, S. W. Randall, W. Forman, R. P. Kraft, G. Brunetti, W. A. Dawson, N. Golovich, C. Jones

Abstract: We present the results from new deep {\it Chandra} ($\sim410$~ks) and {\it Suzaku} ($\sim180$ ks) observations of the merging galaxy cluster ZwCl\,0008.8+5215 ($z=0.104$). Previous radio observations revealed the presence of a double radio relic located diametrically west and east of the cluster center. Using our new {\it Chandra} data, we find evidence for the presence of a shock at the location… ▽ More We present the results from new deep {\it Chandra} ($\sim410$~ks) and {\it Suzaku} ($\sim180$ ks) observations of the merging galaxy cluster ZwCl\,0008.8+5215 ($z=0.104$). Previous radio observations revealed the presence of a double radio relic located diametrically west and east of the cluster center. Using our new {\it Chandra} data, we find evidence for the presence of a shock at the location of the western relic, RW, with a Mach number $\mathcal{M}_{S_X}=1.48^{+0.50}_{-0.32}$ from the density jump. We also measure $\mathcal{M}_{T_X}=2.35^{+0.74}_{-0.55}$ and $\mathcal{M}_{T_X}=2.02^{+0.74}_{-0.47}$ from the temperature jump, with {\it Chandra} and {\it Suzaku} respectively. These values are consistent with the Mach number estimate from a previous study of the radio spectral index, under the assumption of diffusive shock acceleration ($\mathcal{M}_{\rm RW}=2.4^{+0.4}_{-0.2}$). Interestingly, the western radio relic does not entirely trace the X-ray shock. A possible explanation is that the relic traces fossil plasma from nearby radio galaxies which is re-accelerated at the shock. For the eastern relic we do not detect an X-ray surface brightness discontinuity, despite the fact that radio observations suggest a shock with $\mathcal{M}_{\rm RE}=2.2^{+0.2}_{-0.1}$. The low surface brightness and reduced integration time for this region might have prevented the detection. {\it Chandra} surface brightness profile suggests $\mathcal{M}\lesssim1.5$, while {\it Suzaku} temperature measurements found $\mathcal{M}_{T_X}=1.54^{+0.65}_{-0.47}$. Finally, we also detect a merger induced cold front on the western side of the cluster, behind the shock that traces the western relic. △ Less

Submitted 1 February, 2019; originally announced February 2019.

Comments: 21 pages, 15 figures, 4 tables; accepted for publication in ApJ

Multi-wavelength Analysis of the Merging Galaxy Cluster A115

Authors: Mincheol Kim, M. James Jee, Kyle Finner, Nathan Golovich, David M. Wittman, R. J. van Weeren, W. A. Dawson

Abstract: A115 is a merging galaxy cluster at $z\sim0.2$ with a number of remarkable features including a giant ($\sim2.5$ Mpc) radio relic, two asymmetric X-ray peaks with trailing tails, and a peculiar line-of-sight velocity structure. We present a multi-wavelength study of A115 using optical imaging data from Subaru, X-ray data from $Chandra$, and spectroscopic data from the Keck/DEIMOS and MMT/Hectospec… ▽ More A115 is a merging galaxy cluster at $z\sim0.2$ with a number of remarkable features including a giant ($\sim2.5$ Mpc) radio relic, two asymmetric X-ray peaks with trailing tails, and a peculiar line-of-sight velocity structure. We present a multi-wavelength study of A115 using optical imaging data from Subaru, X-ray data from $Chandra$, and spectroscopic data from the Keck/DEIMOS and MMT/Hectospec instruments. Our weak-lensing analysis shows that the cluster is comprised of two subclusters whose mass centroids are in excellent agreement with the two BCG positions ($\lesssim10$"). By modeling A115 with a superposition of two Navarro-Frenk-White halos, we determine the masses of the northern and southern subclusters to be $M_{200}=1.58_{-0.49}^{+0.56}\times 10^{14} \text{M}_{\odot}$ and $3.15_{-0.71}^{+0.79}\times 10^{14} \text{M}_{\odot}$, respectively. Combining the two halos, we estimate the total cluster mass to be $M_{200}=6.41_{-1.04}^{+1.08}\times10^{14} \text{M}_{\odot}$ at $R_{200}=1.67_{-0.09}^{+0.10}$ Mpc. These weak-lensing masses are significantly (a factor of 3-10) lower than what is implied by the X-ray and optical spectroscopic data. We attribute the difference to the gravitational and hydrodynamic disruption caused by the collision between the two subclusters. △ Less

Submitted 5 March, 2019; v1 submitted 20 December, 2018; originally announced December 2018.

Comments: 20 pages, 17 figures, Accepted for publication in ApJ

Unique Science from a Coordinated LSST-WFIRST Survey of the Galactic Bulge

Authors: R. A. Street, M. B. Lund, M. Donachie, S. Khakpash, N. Golovich, M. Penny, D. Bennett, W. A. Dawson, J. Pepper, M. Rabus, P. Szkody, W. I. Clarkson, R. Di Stefano, N. Rattenbury, M. P. G. Hundertmark, Y. Tsapras, S. Ridgway, K. Stassun, V. Bozza, A. Bhattacharya, S. Calchi Novati, Y. Shvartzvald

Abstract: NASA's WFIRST mission will perform a wide-field, NIR survey of the Galactic Bulge to search for exoplanets via the microlensing techniques. As the mission is due to launch in the mid-2020s, around half-way through the LSST Main Survey, we have a unique opportunity to explore synergistic science from two landmark programs. LSST can survey the entire footprint of the WFIRST microlensing survey in a… ▽ More NASA's WFIRST mission will perform a wide-field, NIR survey of the Galactic Bulge to search for exoplanets via the microlensing techniques. As the mission is due to launch in the mid-2020s, around half-way through the LSST Main Survey, we have a unique opportunity to explore synergistic science from two landmark programs. LSST can survey the entire footprint of the WFIRST microlensing survey in a single Deep Drilling Field. Here we explore the great scientific potential of this proposal and recommend the most effective observing strategies. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: White Paper submitted to the 2018 call for LSST cadence proposals. arXiv admin note: text overlap with arXiv:1812.03137

The Diverse Science Return from a Wide-Area Survey of the Galactic Plane

Authors: R. A. Street, M. B. Lund, S. Khakpash, M. Donachie, W. A. Dawson, N. Golovich, L. Wyrzykowski, P. Szkody, T. Naylor, M. Penny, N. Rattenbury, M. Dall'Ora, W. I. Clarkson, D. Bennett, J. Pepper, M. Rabus, M. P. G. Hundertmark, Y. Tsapras, R. Di Stefano, S. Ridgway, M. Liu, E. Lin

Abstract: The overwhelming majority of objects visible to LSST lie within the Galactic Plane. Though many previous surveys have avoided this region for fear of stellar crowding, LSST's spatial resolution combined with its state-of-the-art Difference Image Analysis mean that it can conduct a high cadence survey of most of the Galaxy for the first time. Here we outline the many areas of science that would gre… ▽ More The overwhelming majority of objects visible to LSST lie within the Galactic Plane. Though many previous surveys have avoided this region for fear of stellar crowding, LSST's spatial resolution combined with its state-of-the-art Difference Image Analysis mean that it can conduct a high cadence survey of most of the Galaxy for the first time. Here we outline the many areas of science that would greatly benefit from an LSST survey that included the Galactic Plane, Magellanic Clouds and Bulge at a cadence of 2-3 d. Particular highlights include measuring the mass spectrum of black holes, and mapping the population of exoplanets in the Galaxy in relation to variations in star forming environments. But the same survey data will provide a goldmine for a wide range of science, and we explore possible survey strategies which maximize the scientific return for a number of fields including young stellar objects, cataclysmic variables and Neptune Trojans. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: White Paper submitted in response to the 2018 call for LSST cadence proposals

RELICS: Strong Lensing Analysis of MACS J0417.5-1154 and Predictions for Observing the Magnified High-Redshift Universe with JWST

Authors: Guillaume Mahler, Keren Sharon, Carter Fox, Dan Coe, Mathilde Jauzac, Victoria Strait, Alastair Edge, Ana Acebron, Felipe Andrade-Santos, Roberto J. Avila, Maruša Bradač, Larry D. Bradley, Daniela Carrasco, Catherine Cerny, Nathália Cibirka, Nicole G. Czakon, William A. Dawson, Brenda L. Frye, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Christine Jones, Shotaro Kikuchihara, 15 Daniel Lam, Rachael Livermore , et al. (15 additional authors not shown)

Abstract: Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in p… ▽ More Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST). △ Less

Submitted 10 March, 2019; v1 submitted 31 October, 2018; originally announced October 2018.

Comments: 20 pages, 8 figures. Accepted to ApJ

RELICS: High-Resolution Constraints on the Inner Mass Distribution of the z=0.83 Merging Cluster RXJ0152.7-1357 from strong lensing

Authors: Ana Acebron, May Alon, Adi Zitrin, Guillaume Mahler, Dan Coe, Keren Sharon, Nathália Cibirka, Maruša Bradač, Michele Trenti, Keiichi Umetsu, Felipe Andrade-Santos, Roberto J. Avila, Larry Bradley, Daniela Carrasco, Catherine Cerny, Nicole G. Czakon, William A. Dawson, Brenda Frye, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Christine Jones, Shotaro Kikuchihara, Daniel Lam, Rachael C. Livermore , et al. (15 additional authors not shown)

Abstract: Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by… ▽ More Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive. △ Less

Submitted 26 February, 2019; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: 15 Pages, 7 Figures, 4 Tables Accepted for publication in ApJ

Deep Very Large Array observations of the merging cluster CIZA J2242.8+5301: continuum and spectral imaging

Authors: G. Di Gennaro, R. J. van Weeren, M. Hoeft, H. Kang, D. Ryu, L. Rudnick, W. Forman, H. J. A. Rottgering, M. Bruggen, W. A. Dawson, N. Golovich, D. N. Hoang, H. T. Intema, C. Jones, R. P. Kraft, T. W. Shimwell, A. Stroe

Abstract: Despite progress in understanding radio relics, there are still open questions regarding the underlying particle acceleration mechanisms. In this paper we present deep 1--4 GHz VLA observations of CIZA\,J2242.8+5301 ($z=0.1921$), a double radio relic cluster characterized by small projection on the plane of the sky. Our VLA observations reveal, for the first time, the complex morphology of the dif… ▽ More Despite progress in understanding radio relics, there are still open questions regarding the underlying particle acceleration mechanisms. In this paper we present deep 1--4 GHz VLA observations of CIZA\,J2242.8+5301 ($z=0.1921$), a double radio relic cluster characterized by small projection on the plane of the sky. Our VLA observations reveal, for the first time, the complex morphology of the diffuse sources and the filamentary structure of the northern relic. We discover new faint diffuse radio emission extending north of the main northern relic. Our Mach number estimates for the northern and southern relics, based on the radio spectral index map obtained using the VLA observations and existing LOFAR and GMRT data, are consistent with previous radio and X-ray studies ($\mathcal{M}_{\rm RN}=2.58\pm0.17$ and $\mathcal{M}_{\rm RS}=2.10\pm0.08$). However, color-color diagrams and modelings suggest a flatter injection spectral index than the one obtained from the spectral index map, indicating that projection effects might be not entirely negligible. The southern relic consists of five "arms". Embedded in it, we find a tailed radio galaxy which seems to be connected to the relic. A spectral index flattening, where the radio tail connects to the relic, is also measured. We propose that the southern relic may trace AGN fossil electrons that are re-accelerated at a shock, with an estimated strength of $\mathcal{M}=2.4$. High-resolution mapping of other tailed radio galaxies also supports a scenario where AGN fossil electrons are revived by the merger event and could be related to the formation of some diffuse cluster radio emission. △ Less

Submitted 30 July, 2018; originally announced August 2018.

Comments: 29 pages, 19 figures, 5 tables, accepted for publication in the Astrophysical Journal

Merging Cluster Collaboration: A Panchromatic Atlas of Radio Relic Mergers

Authors: Nathan Golovich, William A. Dawson, David M. Wittman, Reinout J. van Weeren, Felipe Andrade-Santos, M. James Jee, Bryant Benson, Francesco de Gasperin, Tiziana Venturi, Annalisa Bonafede, David Sobral, Georgiana A. Ogrean, Brian C. Lemaux, Maruša Bradač, Marcus Brüggen, Annika H. G. Peter

Abstract: Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of… ▽ More Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of sight velocity difference between likely merging subclusters. We find that the line-of-sight velocity components of the mergers are typically small compared with the maximum three dimensional relative velocity (usually $<1000$ km s$^{-1}$ and often consistent with zero). This suggests that the merger axes of these systems are generally in or near the plane of the sky matching findings in magneto-hydrodynamical simulations. In 28 of the 29 systems we identify substructures in the galaxy population aligned with the radio relic(s) and presumed associated merger induced shock. From this ensemble, we identify eight systems to include in a `gold' sample that is prime for further observation, modeling, and simulation study. Additional papers will present weak lensing mass maps and dynamical modeling for each merging system, ultimately leading to new insight into a wide range of astrophysical phenomena at some of the largest scales in the universe. △ Less

Submitted 27 June, 2018; originally announced June 2018.

Comments: 71 pages, 42 figures, Submitted to ApJ

RELICS: A Strong Lens Model for SPT-CLJ0615-5746, a z=0.972 Cluster

Authors: Rachel Paterno-Mahler, Keren Sharon, Dan Coe, Guillaume Mahler, Catherine Cerny, Traci Johnson, Tim Schrabback, Felipe Andrade-Santos, Roberto J. Avila, Marusa Bradac, Larry D. Bradley, Daniela Carrasco, Nicole G. Czakon, William A. Dawson, Brenda L. Frye, Austin T. Hoag, Kuang-Han Huang, Christine Jones, Daniel Lam, Rachael Livermore, Lorenzo Lovisari, Ramesh Mainali, Pascal A. Oesch, Sara Ogaz, Matthew Past , et al. (8 additional authors not shown)

Abstract: We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use a… ▽ More We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known. △ Less

Submitted 9 July, 2018; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: 15 pages, 8 figures 4 tables. ApJ Accepted

Radio observations of the double-relic galaxy cluster Abell 1240

Authors: D. N. Hoang, T. W. Shimwell, R. J. van Weeren, H. T. Intema, H. J. A. Röttgering, F. Andrade-Santos, H. Akamatsu, A. Bonafede, G. Brunetti, W. A. Dawson, N. Golovich, P. N. Best, A. Botteon, M. Brüggen, R. Cassano, F. de Gasperin, M. Hoeft, A. Stroe, G. J. White

Abstract: We present LOFAR $120-168$ MHz images of the merging galaxy cluster Abell 1240 that hosts double radio relics. In combination with the GMRT $595-629$ MHz and VLA $2-4$ GHz data, we characterised the spectral and polarimetric properties of the radio emission. The spectral indices for the relics steepen from their outer edges towards the cluster centre and the electric field vectors are approximatel… ▽ More We present LOFAR $120-168$ MHz images of the merging galaxy cluster Abell 1240 that hosts double radio relics. In combination with the GMRT $595-629$ MHz and VLA $2-4$ GHz data, we characterised the spectral and polarimetric properties of the radio emission. The spectral indices for the relics steepen from their outer edges towards the cluster centre and the electric field vectors are approximately perpendicular to the major axes of the relics. The results are consistent with the picture that these relics trace large-scale shocks propagating outwards during the merger. Assuming diffusive shock acceleration (DSA), we obtain shock Mach numbers of $\mathcal{M}=2.4$ and $2.3$ for the northern and southern shocks, respectively. For $\mathcal{M}\lesssim3$ shocks, a pre-existing population of mildly relativistic electrons is required to explain the brightness of the relics due to the high ($>10$ per cent) particle acceleration efficiency required. However, for $\mathcal{M}\gtrsim4$ shocks the required efficiency is $\gtrsim1\%$ and $\gtrsim0.5\%$, respectively, which is low enough for shock acceleration directly from the thermal pool. We used the fractional polarization to constrain the viewing angle to $\geqslant(53\pm3)^\circ$ and $\geqslant(39\pm5)^\circ$ for the northern and southern shocks, respectively. We found no evidence for diffuse emission in the cluster central region. If the halo spans the entire region between the relics ($\sim1.8\,\text{Mpc}$) our upper limit on the power is $P_\text{1.4 GHz}=(1.4\pm0.6)\times10^{23}\,\text{W}\text{Hz}^{-1}$ which is approximately equal to the anticipated flux from a cluster of this mass. However, if the halo is smaller than this, our constraints on the power imply that the halo is underluminous. △ Less

Submitted 30 April, 2018; originally announced April 2018.

RELICS: Strong Lensing analysis of the galaxy clusters Abell S295, Abell 697, MACS J0025.4-1222, and MACS J0159.8-0849

Authors: Nathália Cibirka, Ana Acebron, Adi Zitrin, Dan Coe, Irene Agulli, Felipe Andrade-Santos, Maruša Bradač, Brenda Frye, Rachael C. Livermore, Guillaume Mahler, Brett Salmon, Keren Sharon, Michele Trenti, Keiichi Umetsu, Roberto Avila, Larry Bradley, Daniela Carrasco, Catherine Cerny, Nicole G. Czakon, William A. Dawson, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Christine Jones, Shotaro Kikuchihara , et al. (15 additional authors not shown)

Abstract: We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models f… ▽ More We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($μ>5$) $\sim 1-3$ arcmin$^2$, A($μ>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$α$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$α$ rest-frame equivalent width of $52\pm22$ Å, and an observed Gaussian width of $117\pm 15$ km/s. △ Less

Submitted 11 July, 2018; v1 submitted 26 March, 2018; originally announced March 2018.

Comments: 23 pages, 16 figures; V2, accepted for publication in ApJ

RELICS: Strong-lensing analysis of the massive clusters MACS J0308.9+2645 and PLCK G171.9-40.7

Authors: Ana Acebron, Nathália Cibirka, Adi Zitrin, Dan Coe, Irene Agulli, Keren Sharon, Maruša Bradač, Brenda Frye, Rachael C. Livermore, Guillaume Mahler, Brett Salmon, Keiichi Umetsu, Larry Bradley, Felipe Andrade-Santos, Roberto Avila, Daniela Carrasco, Catherine Cerny, Nicole G. Czakon, William A. Dawson, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Christine Jones, Shotaro Kikuchihara, Daniel Lam , et al. (15 additional authors not shown)

Abstract: Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK… ▽ More Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($θ_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $μ>5$ to $\sim1$ arcmin$^{2}$ with $μ>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes. △ Less

Submitted 2 May, 2018; v1 submitted 1 March, 2018; originally announced March 2018.

Comments: 14 pages, 6 figures, 4 tables, accepted for publication in ApJ

A Candidate $z\sim10$ Galaxy Strongly Lensed into a Spatially Resolved Arc

Authors: Brett Salmon, Dan Coe, Larry Bradley, Marusa Bradač, Kuang-Han Huang, Victoria Strait, Pascal Oesch, Rachel Paterno-Mahler, Adi Zitrin, Ana Acebron, Nathália Cibirka, Shotaro Kikuchihara, Masamune Oguri, Gabriel B. Brammer, Keren Sharon, Michele Trenti, Roberto J. Avila, Sara Ogaz, Felipe Andrade-Santos, Daniela Carrasco, Catherine Cerny, William Dawson, Brenda L. Frye, Austin Hoag, Christine Jones , et al. (5 additional authors not shown)

Abstract: The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lens… ▽ More The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of μ~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: 7 pages, 4 figures. Submitted to ApJ Letters

Deep VLA observations of the cluster 1RXS J0603.3+4214 in the frequency range 1-2 GHz

Authors: K. Rajpurohit, M. Hoeft, R. J. van Weeren, L. Rudnick, H. J. A. Röttgering, W. R. Forman, M. Brüggen, J. H. Croston, F. Andrade-Santos, W. A. Dawson, H. T. Intema, R. P. Kraft, C. Jones, M. James Jee

Abstract: We report L-band VLA observations of 1RXS J0603.3+4214, a cluster that hosts a bright radio relic, known as the Toothbrush, and an elongated giant radio halo. These new observations allow us to study the surface brightness distribution down to one arcsec resolution with very high sensitivity. Our images provide an unprecedented detailed view of the Toothbrush, revealing enigmatic filamentary struc… ▽ More We report L-band VLA observations of 1RXS J0603.3+4214, a cluster that hosts a bright radio relic, known as the Toothbrush, and an elongated giant radio halo. These new observations allow us to study the surface brightness distribution down to one arcsec resolution with very high sensitivity. Our images provide an unprecedented detailed view of the Toothbrush, revealing enigmatic filamentary structures. To study the spectral index distribution, we complement our analysis with published LOFAR and GMRT observations. The bright `brush' of the Toothbrush shows a prominent narrow ridge to its north with a sharp outer edge. The spectral index at the ridge is in the range $-0.70\leqα\leq-0.80$. We suggest that the ridge is caused by projection along the line of sight. With a simple toy model for the smallest region of the ridge, we conclude that the magnetic field is below $5\,\rmμG$ and varies significantly across the shock front. Our model indicates that the actual Mach number is higher than that obtained from the injection index and agrees well with the one derived from the overall spectrum, namely ${\cal M}=3.78^{+0.3}_{-0.2}$. The radio halo shows an average spectral index of $α=-1.16\pm0.05$ and a slight gradient from north to south. The southernmost part of the halo is steeper and possibly related to a shock front. Excluding the southernmost part, the halo morphology agrees very well with the X-ray morphology. A power-law correlation is found between the radio and X-ray surface brightness △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: 23 pages, 21 figures, accepted for publication in ApJ

Merging Cluster Collaboration: Optical and Spectroscopic Survey of a Radio-Selected Sample of Twenty Nine Merging Galaxy Clusters

Authors: Nathan Golovich, William A. Dawson, David M. Wittman, Myungkook J. Jee, Bryant Benson, Brian Lemaux, Reinout J. van Weeren, Felipe Andrade-Santos, David Sobral, Francesco de Gasperin, Marcus Bruggen, Marusa Bradac, Kyle Finner, Annika Peter

Abstract: Multi-band photometric and multi-object spectroscopic surveys of merging galaxy clusters allow for the characterization of the distributions of constituent dark matter and galaxy populations, constraints on the dynamics of the merging subclusters, and an understanding of galaxy evolution of member galaxies. We present deep photometric observations from Subaru/SuprimeCam and a catalog of $\sim$5400… ▽ More Multi-band photometric and multi-object spectroscopic surveys of merging galaxy clusters allow for the characterization of the distributions of constituent dark matter and galaxy populations, constraints on the dynamics of the merging subclusters, and an understanding of galaxy evolution of member galaxies. We present deep photometric observations from Subaru/SuprimeCam and a catalog of $\sim$5400 spectroscopic cluster members from Keck/DEIMOS across 29 merging galaxy clusters ranging in redshift from $z=0.07$ to $0.55$. The ensemble is compiled based on the presence of radio relics, which highlight cluster scale collisionless shocks in the intra-cluster medium. Together with the spectroscopic and photometric information, the velocities, timescales, and geometries of the respective merging events may be tightly constrained. In this preliminary analysis, the velocity distributions of 28 of the 29 clusters are shown to be well fit by single Gaussians. This indicates that radio relic mergers largely occur transverse to the line of sight and/or near apocenter. In this paper, we present our optical and spectroscopic surveys, preliminary results, and a discussion of the value of radio relic mergers for developing accurate dynamical models of each system. △ Less

Submitted 3 November, 2017; originally announced November 2017.

Comments: 20 pages, 7 figures, 7 tables, Submitted to the Astrophysical Journal Supplements on 3 November 2017

RELICS: Strong Lens Models for Five Galaxy Clusters From the Reionization Lensing Cluster Survey

Authors: Catherine Cerny, Keren Sharon, Felipe Andrade-Santos, Roberto J. Avila, Marusa Bradac, Larry D. Bradley, Daniela Carrasco, Dan Coe, Nicole G. Czakon, William A. Dawson, Brenda L. Frye, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Christine Jones, Daniel Lam, Lorenzo Lovisari, Ramesh Mainali, Pascal A. Oesch, Sara Ogaz, Matthew Past, Rachel Paterno-Mahler, Avery Peterson, Adam G. Riess, Steven A. Rodney , et al. (9 additional authors not shown)

Abstract: Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imagin… ▽ More Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imaging of the clusters RXC J0142.9+4438, Abell 2537, Abell 2163, RXC J2211.7-0349, and ACT-CLJ0102-49151. The derived lensing magnification is essential for estimating the intrinsic properties of high-redshift galaxy candidates, and properly accounting for the survey volume. We report on new spectroscopic redshifts of multiply imaged lensed galaxies behind these clusters, which are used as constraints, and detail our strategy to reduce systematic uncertainties due to lack of spectroscopic information. In addition, we quantify the uncertainty on the lensing magnification due to statistical and systematic errors related to the lens modeling process, and find that in all but one cluster, the magnification is constrained to better than 20% in at least 80% of the field of view, including statistical and systematic uncertainties. The five clusters presented in this paper span the range of masses and redshifts of the clusters in the RELICS program. We find that they exhibit similar strong lensing efficiencies to the clusters targeted by the Hubble Frontier Fields within the WFC3/IR field of view. Outputs of the lens models are made available to the community through the Mikulski Archive for Space Telescopes △ Less

Submitted 15 April, 2018; v1 submitted 25 October, 2017; originally announced October 2017.

Comments: Accepted to ApJ

The Reionization Lensing Cluster Survey (RELICS) and the Brightest High-z Galaxies

Authors: Brett Salmon, Dan Coe, Larry Bradley, Rychard Bouwens, Marusa Bradac, Kuang-Han Huang, Pascal Oesch, Daniel Stark, Keren Sharon, Michele Trenti, Roberto J. Avila, Sara Ogaz, Felipe Andrade-Santos, Daniela Carrasco, Catherine Cerny, William Dawson, Brenda L. Frye, Austin Hoag, Traci Lin Johnson, Christine Jones, Daniel Lam, Lorenzo Lovisari, Ramesh Mainali, Matt Past, Rachel Paterno-Mahler , et al. (9 additional authors not shown)

Abstract: Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200… ▽ More Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope. △ Less

Submitted 20 April, 2018; v1 submitted 24 October, 2017; originally announced October 2017.

Comments: 20 pages, 15 figures, submitted to ApJ. v2 edit to Table 1

MC$^2$: Subaru and Hubble Space Telescope Weak-Lensing Analysis of the Double Radio Relic Galaxy Cluster PLCK G287.0+32.9

Authors: Kyle Finner, M. James Jee, Nathan Golovich, David Wittman, William Dawson, Daniel Gruen, Anton M. Koekemoer, Brian C. Lemaux, Stella Seitz

Abstract: The second most significant detection of the Planck Sunyaev Zel'dovich survey, PLCK~G287.0+32.9 ($z=0.385$) boasts two similarly bright radio relics and a radio halo. One radio relic is located $\sim 400$ kpc northwest of the X-ray peak and the other $\sim 2.8$ Mpc to the southeast. This large difference suggests that a complex merging scenario is required. A key missing puzzle for the merging sce… ▽ More The second most significant detection of the Planck Sunyaev Zel'dovich survey, PLCK~G287.0+32.9 ($z=0.385$) boasts two similarly bright radio relics and a radio halo. One radio relic is located $\sim 400$ kpc northwest of the X-ray peak and the other $\sim 2.8$ Mpc to the southeast. This large difference suggests that a complex merging scenario is required. A key missing puzzle for the merging scenario reconstruction is the underlying dark matter distribution in high resolution. We present a joint Subaru Telescope and {\it Hubble Space Telescope} weak-lensing analysis of the cluster. Our analysis shows that the mass distribution features four significant substructures. Of the substructures, a primary cluster of mass $M_{200\text{c}}=1.59^{+0.25}_{-0.22}\times 10^{15} \ h^{-1}_{70} \ \text{M}_{\odot}$ dominates the weak-lensing signal. This cluster is likely to be undergoing a merger with one (or more) subcluster whose mass is approximately a factor of 10 lower. One candidate is the subcluster of mass $M_{200\text{c}}=1.16^{+0.15}_{-0.13}\times 10^{14} \ h^{-1}_{70} \ \text{M}_{\odot}$ located $\sim 400$ kpc to the southeast. The location of this subcluster suggests that its interaction with the primary cluster could be the source of the NW radio relic. Another subcluster is detected $\sim 2$ Mpc to the SE of the X-ray peak with mass $M_{200\text{c}}=1.68^{+0.22}_{-0.20}\times 10^{14} \ h^{-1}_{70} \ \text{M}_{\odot}$. This SE subcluster is in the vicinity of the SE radio relic and may have created the SE radio relic during a past merger with the primary cluster. The fourth subcluster, $M_{200\text{c}}=1.87^{+0.24}_{-0.22}\times 10^{14} \ h^{-1}_{70} \ \text{M}_{\odot}$, is northwest of the X-ray peak and beyond the NW radio relic. △ Less

Submitted 17 November, 2017; v1 submitted 6 October, 2017; originally announced October 2017.

Comments: 19 pages, 14 figures; Accepted to ApJ

US Cosmic Visions: New Ideas in Dark Matter 2017: Community Report

Authors: Marco Battaglieri, Alberto Belloni, Aaron Chou, Priscilla Cushman, Bertrand Echenard, Rouven Essig, Juan Estrada, Jonathan L. Feng, Brenna Flaugher, Patrick J. Fox, Peter Graham, Carter Hall, Roni Harnik, JoAnne Hewett, Joseph Incandela, Eder Izaguirre, Daniel McKinsey, Matthew Pyle, Natalie Roe, Gray Rybka, Pierre Sikivie, Tim M. P. Tait, Natalia Toro, Richard Van De Water, Neal Weiner , et al. (226 additional authors not shown)

Abstract: This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017. This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017. △ Less

Submitted 14 July, 2017; originally announced July 2017.

Comments: 102 pages + references

MC^2: A deeper look at ZwCl 2341.1+0000 with Bayesian galaxy clustering and weak lensing analyses

Authors: Bryant Benson, David M. Wittman, Nathan Golovich, M. James Jee, Reinout J. van Weeren, William A. Dawson

Abstract: ZwCl 2341.1+0000, a merging galaxy cluster with disturbed X-ray morphology and widely separated ($\sim$3 Mpc) double radio relics, was thought to be an extremely massive ($10-30 \times 10^{14} M_\odot$) and complex system with little known about its merger history. We present JVLA 2-4 GHz observations of the cluster, along with new spectroscopy from our Keck/DEIMOS survey, and apply Gaussian Mixtu… ▽ More ZwCl 2341.1+0000, a merging galaxy cluster with disturbed X-ray morphology and widely separated ($\sim$3 Mpc) double radio relics, was thought to be an extremely massive ($10-30 \times 10^{14} M_\odot$) and complex system with little known about its merger history. We present JVLA 2-4 GHz observations of the cluster, along with new spectroscopy from our Keck/DEIMOS survey, and apply Gaussian Mixture Modeling to the three-dimensional distribution of 227 confirmed cluster galaxies. After adopting the Bayesian Information Criterion to avoid overfitting, which we discover can bias total dynamical mass estimates high, we find that a three-substructure model with a total dynamical mass estimate of $9.39 \pm 0.81 \times 10^{14} M_\odot$ is favored. We also present deep Subaru imaging and perform the first weak lensing analysis on this system, obtaining a weak lensing mass estimate of $5.57 \pm 2.47 \times 10^{14} M_\odot$. This is a more robust estimate because it does not depend on the dynamical state of the system, which is disturbed due to the merger. Our results indicate that ZwCl 2341.1+0000 is a multiple merger system comprised of at least three substructures, with the main merger that produced the radio relics occurring near to the plane of the sky, and a younger merger in the North occurring closer to the line of sight. Dynamical modeling of the main merger reproduces observed quantities (relic positions and polarizations, subcluster separation and radial velocity difference), if the merger axis angle of $\sim$10$^{+34}_{-6}$ degrees and the collision speed at pericenter is $\sim$1900$^{+300}_{-200}$ km/s. △ Less

Submitted 30 June, 2017; originally announced July 2017.

Journal ref: The Astrophysical Journal, Volume 841, Number 1, 2017

Suppressed Star Formation by a Merging Cluster System

Authors: Alison S. Mansheim, Brian C. Lemaux, Adam R. Tomczak, Lori M. Lubin, Nicholas Rumbaugh, Po-Feng Wu, Roy R. Gal, Lu Shen, William A. Dawson, Gordon K. Squires

Abstract: We examine the effects of an impending cluster merger on galaxies in the large scale structure (LSS) RX J0910 at $z =1.105$. Using multi-wavelength data, including 102 spectral members drawn from the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey and precise photometric redshifts, we calculate star formation rates and map the specific star formation rate density of… ▽ More We examine the effects of an impending cluster merger on galaxies in the large scale structure (LSS) RX J0910 at $z =1.105$. Using multi-wavelength data, including 102 spectral members drawn from the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey and precise photometric redshifts, we calculate star formation rates and map the specific star formation rate density of the LSS galaxies. These analyses along with an investigation of the color-magnitude properties of LSS galaxies indicate lower levels of star formation activity in the region between the merging clusters relative to the outskirts of the system. We suggest that gravitational tidal forces due to the potential of the merging halos may be the physical mechanism responsible for the observed suppression of star formation in galaxies caught between the merging clusters. △ Less

Submitted 12 May, 2017; v1 submitted 9 May, 2017; originally announced May 2017.

Comments: 6 pages, 2 figures, 1 table, accepted in MNRAS Letters, v2 corrected a few typographical errors

Efficient Computation of Sparse Matrix Functions for Large-Scale Electronic Structure Calculations: The CheSS Library

Authors: Stephan Mohr, William Dawson, Michael Wagner, Damien Caliste, Takahito Nakajima, Luigi Genovese

Abstract: We present CheSS, the "Chebyshev Sparse Solvers" library, which has been designed to solve typical problems arising in large-scale electronic structure calculations using localized basis sets. The library is based on a flexible and efficient expansion in terms of Chebyshev polynomials and presently features the calculation of the density matrix, the calculation of matrix powers for arbitrary power… ▽ More We present CheSS, the "Chebyshev Sparse Solvers" library, which has been designed to solve typical problems arising in large-scale electronic structure calculations using localized basis sets. The library is based on a flexible and efficient expansion in terms of Chebyshev polynomials and presently features the calculation of the density matrix, the calculation of matrix powers for arbitrary powers, and the extraction of eigenvalues in a selected interval. CheSS is able to exploit the sparsity of the matrices and scales linearly with respect to the number of nonzero entries, making it well-suited for large-scale calculations. The approach is particularly adapted for setups leading to small spectral widths of the involved matrices and outperforms alternative methods in this regime. By coupling CheSS to the DFT code BigDFT, we show that such a favorable setup is indeed possible in practice. In addition, the approach based on Chebyshev polynomials can be massively parallelized, and CheSS exhibits excellent scaling up to thousands of cores even for relatively small matrix sizes. △ Less

Submitted 10 October, 2017; v1 submitted 3 April, 2017; originally announced April 2017.