Showing 1–41 of 41 results for author: Golovich, N

Investigating the mixing between two black hole populations in LIGO-Virgo-KAGRA GWTC-3

Authors: Ming-Feng Ho, Scott Ellis Perkins, Simeon Bird, William Dawson, Nathan Golovich, Jessica R. Lu, Peter McGill

Abstract: We introduce a population model to analyze the mixing between hypothesised power-law and $\sim 35 M_\odot$ Gaussian bump black hole populations in the latest gravitational wave catalog, GWTC-3, estimating their co-location and separation. We find a relatively low level of mixing, $3.1^{+5.0}_{-3.1}\%$, between the power-law and Gaussian populations, compared to the percentage of mergers containing… ▽ More We introduce a population model to analyze the mixing between hypothesised power-law and $\sim 35 M_\odot$ Gaussian bump black hole populations in the latest gravitational wave catalog, GWTC-3, estimating their co-location and separation. We find a relatively low level of mixing, $3.1^{+5.0}_{-3.1}\%$, between the power-law and Gaussian populations, compared to the percentage of mergers containing two Gaussian bump black holes, $5.0^{+3.2}_{-1.7}\%$. Our analysis indicates that black holes within the Gaussian bump are generally separate from the power-law population, with only a minor fraction engaging in mixing and contributing to the $\mathcal{M} \sim 14 M_\odot$ peak in the chirp mass. This leads us to identify a distinct population of Binary Gaussian Black Holes (BGBHs) that arise from mergers within the Gaussian bump. We suggest that current theories for the formation of the massive $35 M_\odot$ Gaussian bump population may need to reevaluate the underlying mechanisms that drive the preference for BGBHs. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Comments: 19 pages, 11 figures. Submitted to PRD

High-quality Extragalactic Legacy-field Monitoring (HELM) with DECam

Authors: Ming-Yang Zhuang, Qian Yang, Yue Shen, Monika Adamow, Douglas N. Friedel, R. A. Gruendl, Xin Liu, Paul Martini, Timothy M. C. Abbott, Scott F. Anderson, Roberto J. Assef, Franz E. Bauer, Rich Bielby, W. N. Brandt, Colin J. Burke, Jorge Casares, Yu-Ching Chen, Gisella De Rosa, Alex Drlica-Wagner, Tom Dwelly, Alice Eltvedt, Gloria Fonseca Alvarez, Jianyang Fu, Cesar Fuentes, Melissa L. Graham , et al. (23 additional authors not shown)

Abstract: High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the… ▽ More High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the $(u)gri(z)$ bands, over a collective sky area of $\sim 38$ deg${\rm ^2}$. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependences on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first $\sim 3.5$ years of light curves during 2019A--2022A. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 11 pages, 8 figures, 4 tables, submitted to ApJS. Median source catalogs and light curves of individual objects are publicly available at https://ariel.astro.illinois.edu/helm/

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

Distribution of H-alpha Emitters in Merging Galaxy Clusters

Authors: David Wittman, Dariush Imani, Rutger Hartmann Olden, Nathan Golovich

Abstract: Studies of star formation in various galaxy cluster mergers have reached apparently contradictory conclusions regarding whether mergers stimulate star formation, quench it, or have no effect. Because the mergers studied span a range of time since pericenter (TSP), it is possible that the apparent effect on star formation is a function of TSP. We use a sample of 12 bimodal mergers to assess the sta… ▽ More Studies of star formation in various galaxy cluster mergers have reached apparently contradictory conclusions regarding whether mergers stimulate star formation, quench it, or have no effect. Because the mergers studied span a range of time since pericenter (TSP), it is possible that the apparent effect on star formation is a function of TSP. We use a sample of 12 bimodal mergers to assess the star formation as a function of TSP. We measure the equivalent width of the H-alpha emission line in ${\sim}100$ member galaxies in each merger, classify galaxies as emitters or non-emitters, and then classify emitters as star-forming galaxies (SFG) or active galactic nucleus (AGN) based on the [NII] $\lambda6583$ line. We quantify the distribution of SFG and AGN relative to non-emitters along the spatial axis defined by the subcluster separation. The SFG and AGN fractions vary from merger to merger, but show no trend with TSP. The spatial distribution of SFG is consistent with that of non-emitters in eight mergers, but show significant avoidance of the system center in the remaining four mergers, including the three with the lowest TSP. If there is a connection between star formation activity and TSP, probing it further will require more precise TSP estimates and more mergers with TSP in the range of 0-400 Myr. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: accepted to AJ

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

MEGASIM: Distribution and Detection of Earth Trojan Asteroids

Authors: Travis Yeager, Nathan Golovich, Kerianne Pruett

Abstract: Using N-body simulation results from the MEGASIM dataset, we present spatial distributions of Earth Trojan Asteroids and assess the detectability of the population in current and next-generation ground-based astronomical surveys (Yeager & Golovich 2022). Our high-fidelity Earth Trojan Asteroid (ETA) distribution maps show never-before-seen high-resolution spatial features that evolve over timescal… ▽ More Using N-body simulation results from the MEGASIM dataset, we present spatial distributions of Earth Trojan Asteroids and assess the detectability of the population in current and next-generation ground-based astronomical surveys (Yeager & Golovich 2022). Our high-fidelity Earth Trojan Asteroid (ETA) distribution maps show never-before-seen high-resolution spatial features that evolve over timescales up to 1 Gyr. The simulation was synchronized to start times and timelines of two observational astronomy surveys, 1) the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) and 2) the Zwicky Transient Facility (ZTF). We calculate upper limits for the number of ETAs potentially observable with both the ZTF and LSST surveys. Due to the Yarkovsky Effect, we find no stable ETAs on billion year timescales likely to be detected by any ETA survey, as no C-type or S-type ETAs (with $H<22$ and $H<24$, respectively) are likely to be stable on billion year timescales, and ETAs large enough to remain stable on billion year timescales are very rare relative to the rest of the ETA population. We find that a twilight ETA survey will not drastically increase the likelihood of individual ETA detection, but would provide orders of magnitude more observations of select ETA populations. The null detection to date from ZTF restricts the potential ETA population to hundreds of objects larger than 100 meters (at $H\approx22$, Lifset et al. 2021), while a null detection by LSST will further restrict the ETA population to tens of objects larger than 100 meters. △ Less

Submitted 7 November, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

On the Need for a Near-Earth Object Characterization Constellation in Low-Earth Orbit

Authors: Nathan Golovich

Abstract: In 2005, the United States Congress passed a bill mandating the detection, tracking, cataloguing and characterization of 90\% of the 140 meter and larger near-Earth objects (NEOs) by 2020. At the deadline $\sim$35\% were detected, tracked and catalogued, but only a small fraction were characterized. At the present rate, it will take 40 years to meet the detection mandate, and there are insufficien… ▽ More In 2005, the United States Congress passed a bill mandating the detection, tracking, cataloguing and characterization of 90\% of the 140 meter and larger near-Earth objects (NEOs) by 2020. At the deadline $\sim$35\% were detected, tracked and catalogued, but only a small fraction were characterized. At the present rate, it will take 40 years to meet the detection mandate, and there are insufficient global facilities dedicated to NEO characterization to come close to the characterization threshold. The major surveys focus mainly on detection and initial orbit determination, which must be refined in order to fully be tracked and catalogued. Characterization requires observations spanning multiple wavelengths, cadences, and instruments, so it is challenging for observers to acquire the requisite data in a timely manner for planetary defense. Two upcoming surveys will easily meet the 90\% threshold for detection, but each will require separate facilities to tip and queue to refine orbits and characterize new discoveries, and they will provide too many discoveries for ground and space-based assets to keep up with. Here, I argue for a constellation of proliferating small satellites carrying visible and infrared sensors that would offer the needed coverage and flexibility to follow up detections from current and upcoming surveys in a timely manner. Such a constellation would enable NASA to move beyond the detection focused investments and fully meet the 2005 Congressional mandate. △ Less

Submitted 26 January, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: Submitted to Icarus, 27 pages, 13 figures, silly title typo fixed

Report number: LLNL-JRNL-843984

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.

MEGASIM: Lifetimes and Resonances of Earth Trojan Asteroids -- The Death of Primordial ETAs?

Authors: Travis Yeager, Nathan Golovich

Abstract: We present an analysis of lifetimes and resonances of Earth Trojan Asteroids (ETAs) in the MEGASIM data set (Yeager & Golovich 2022). Trojan asteroids co-orbit the Sun with a planet but remain bound to the Lagrange points, L4 (60° leading the planet) or L5 (60° trailing). In the circular three-body approximation, the stability of a Trojan asteroid depends on the ratio of the host planet mass and t… ▽ More We present an analysis of lifetimes and resonances of Earth Trojan Asteroids (ETAs) in the MEGASIM data set (Yeager & Golovich 2022). Trojan asteroids co-orbit the Sun with a planet but remain bound to the Lagrange points, L4 (60° leading the planet) or L5 (60° trailing). In the circular three-body approximation, the stability of a Trojan asteroid depends on the ratio of the host planet mass and the central mass. For the inner planets, the range of stability becomes increasingly small, so perturbations from the planets have made primordial Trojans rare. To date there have been just two ETAs (2010 TK7 and 2020 XL5), several Mars Trojans, and a Venus Trojan discovered. The estimated lifetimes of the known inner system Trojans are less than a million years, suggesting they are interlopers rather than members of a stable and long-lasting population. With the largest ETA n-body simulation to date, we are able to track their survival across a wide initialized parameter space. We find the remaining fraction of ETAs over time is well fit with a stretched exponential function that when extrapolated beyond our simulation run time predicts zero ETAs by 2.33 Gyr. We also show correlations between ETA ejections and the periods of the Milankovitch cycles. Though Earth's orbital dynamics dominate the instabilities of ETAs, we provide evidence that ETA ejections are linked to resonances found in the variation of the orbital elements of many, if not all of the planets. △ Less

Submitted 3 June, 2022; originally announced June 2022.

A New Blind Asteroid Detection Scheme

Authors: Nathan Golovich, Noah Lifset, Robert Armstrong, Eric Green, Michael D. Schneider, Roger Pearce

Abstract: As astronomical photometric surveys continue to tile the sky repeatedly, the potential to pushdetection thresholds to fainter limits increases; however, traditional digital-tracking methods cannotachieve this efficiently beyond time scales where motion is approximately linear. In this paper weprototype an optimal detection scheme that samples under a user defined prior on a parameterizationof the… ▽ More As astronomical photometric surveys continue to tile the sky repeatedly, the potential to pushdetection thresholds to fainter limits increases; however, traditional digital-tracking methods cannotachieve this efficiently beyond time scales where motion is approximately linear. In this paper weprototype an optimal detection scheme that samples under a user defined prior on a parameterizationof the motion space, maps these sampled trajectories to the data space, and computes an optimalsignal-matched filter for computing the signal to noise ratio of trial trajectories. We demonstrate thecapability of this method on a small test data set from the Dark Energy Camera. We recover themajority of asteroids expected to appear and also discover hundreds of new asteroids with only a fewhours of observations. We conclude by exploring the potential for extending this scheme to larger datasets that cover larger areas of the sky over longer time baselines. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 28 pages, 12 figures

Time-series and Phasecurve Photometry of Episodically-Active Asteroid (6478) Gault in a Quiescent State Using APO, GROWTH, P200 and ZTF

Authors: Josiah N. Purdum, Zhong-Yi Lin, Bryce T. Bolin, Kritti Sharma, Philip I. Choi, Varun Bhalerao, Harsh Kumar, Robert Quimby, Joannes C. Van Roestel, Chengxing Zhai, Yanga R. Fernandez, Josef Hanuš, Carey M. Lisse, Dennis Bodewits, Christoffer Fremling, Nathan Ryan Golovich, Chen-Yen Hsu, Wing-Huen Ip, Chow-Choong Ngeow, Navtej S. Saini, Michael Shao, Yuhan Yao, Tomás Ahumada, Shreya Anand, Igor Andreoni , et al. (27 additional authors not shown)

Abstract: We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations.… ▽ More We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations. In addition, we were able to constrain Gault's rotation period using time-series photometric lightcurves taken over 17 hours on multiple days in 2020 August, September and October. The photometric lightcurves have a repeating $\lesssim$0.05 magnitude feature suggesting that (6478) Gault has a rotation period of ~2.5 hours and may have a semi-spherical or top-like shape, much like Near-Earth Asteroids Ryugu and Bennu. The rotation period of ~2.5 hours is near to the expected critical rotation period for an asteroid with the physical properties of (6478) Gault suggesting that its activity observed over multiple epochs is due to surface mass shedding from its fast rotation spun up by the Yarkovsky-O'Keefe-Radzievskii-Paddack effect. △ Less

Submitted 27 March, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: 23 pages, 18 figures; Accepted by ApJ Letters

A Search for L4 Earth Trojan Asteroids Using a Novel Track-Before-Detect Multi-Epoch Pipeline

Authors: Noah Lifset, Nathan Golovich, Eric Green, Robert Armstrong, Travis Yeager

Abstract: Earth Trojan Asteroids are an important but elusive population that co-orbit with Earth at the L4 and L5 Lagrange points. There is only one known, but a large population is theoretically stable and could provide insight into our solar system's past and present as well as planetary defense. In this paper, we present the results of an Earth Trojan survey that uses a novel shift-and-stack detection m… ▽ More Earth Trojan Asteroids are an important but elusive population that co-orbit with Earth at the L4 and L5 Lagrange points. There is only one known, but a large population is theoretically stable and could provide insight into our solar system's past and present as well as planetary defense. In this paper, we present the results of an Earth Trojan survey that uses a novel shift-and-stack detection method on two nights of data from the Dark Energy Camera. We find no new Earth Trojan Asteroids. We calculate an upper limit on the population that is consistent with previous searches despite much less sky coverage. Additionally, we elaborate on previous upper limit calculations using current asteroid population statistics and an extensive asteroid simulation to provide the most up to date population constraints. We find an L4 Earth Trojan population of NET < 1 for H = 13.93, NET < 7 for H = 16, and NET < 938 for H = 22. △ Less

Submitted 24 June, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Journal ref: The Astronomical Journal, 161:282 (11pp), 2021 June

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

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

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

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

Particle re-acceleration and Faraday-complex structures in the RXC J1314.4-2515 galaxy cluster

Authors: C. Stuardi, A. Bonafede, D. Wittor, F. Vazza, A. Botteon, N. Locatelli, D. Dallacasa, N. Golovich, M. Hoeft, R. J. van Weeren, M. Brüggen, F. de Gasperin

Abstract: Radio relics are sites of electron (re)acceleration in merging galaxy clusters but the mechanism of acceleration and the topology of the magnetic field in and near relics are yet to be understood. We are carrying out an observational campaign on double relic galaxy clusters starting with RXC J1314.4-2515. With $Jansky Very Large Array$ multi-configuration observations in the frequency range 1-4 GH… ▽ More Radio relics are sites of electron (re)acceleration in merging galaxy clusters but the mechanism of acceleration and the topology of the magnetic field in and near relics are yet to be understood. We are carrying out an observational campaign on double relic galaxy clusters starting with RXC J1314.4-2515. With $Jansky Very Large Array$ multi-configuration observations in the frequency range 1-4 GHz, we perform both spectral and polarization analyses, using the Rotation Measure synthesis technique. We use archival $XMM-Newton$ observations to constrain the properties of the shocked region. We discover a possible connection between the activity of a radio galaxy and the emission of the eastern radio relic. In the northern elongated arc of the western radio relic, we detect polarized emission with an average polarization fraction of $31 \ \%$ at 3 GHz and we derive the Mach number of the underlying X-ray shock. Our observations reveal low levels of fractional polarization and Faraday-complex structures in the southern region of the relic, which point to the presence of thermal gas and filamentary magnetic field morphology inside the radio emitting volume. We measured largely different Rotation Measure dispersion from the two relics. Finally, we use cosmological magneto-hydrodynamical simulations to constrain the magnetic field, viewing angle, and to derive the acceleration efficiency of the shock. We find that the polarization properties of RXC J1314.4-2515 are consistent with a radio relic observed at $70^{\circ}$ with respect to the line of sight and that efficient re-acceleration of fossil electrons has taken place. △ Less

Submitted 30 August, 2019; originally announced August 2019.

Comments: 24 pages, 19 figures, 5 tables, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 489, Issue 3, November 2019, Pages 3905-3926

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

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

Brightest Cluster Galaxy Alignments in Merging Clusters

Authors: David Wittman, Drake Foote, Nathan Golovich

Abstract: The orientations of brightest cluster galaxies (BCGs) and their host clusters tend to be aligned, but the mechanism driving this is not clear. To probe the role of cluster mergers in this process, we quantify alignments of 38 BCGs in 22 clusters undergoing major mergers (up to $\sim 1$ Gyr after first pericenter). We find alignments entirely consistent with those of clusters in general. This sugge… ▽ More The orientations of brightest cluster galaxies (BCGs) and their host clusters tend to be aligned, but the mechanism driving this is not clear. To probe the role of cluster mergers in this process, we quantify alignments of 38 BCGs in 22 clusters undergoing major mergers (up to $\sim 1$ Gyr after first pericenter). We find alignments entirely consistent with those of clusters in general. This suggests that alignments are robust against major cluster mergers. If, conversely, major cluster mergers actually help orient the BCG, such a process is acting quickly because the orientation is in place within $\sim 1$ Gyr after first pericenter. △ Less

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

Comments: accepted to 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

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.

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

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

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

On the absence of radio halos in clusters with double relics

Authors: A. Bonafede, R. Cassano, M. Brüggen, G. A. Ogrean, C. J. Riseley, V. Cuciti, F. de Gasperin, N. Golovich, R. Kale, T. Venturi, R. J van Weeren, D. R. Wik, D. Wittman

Abstract: Pairs of radio relics are believed to form during cluster mergers, and are best observed when the merger occurs in the plane of the sky. Mergers can also produce radio halos, through complex processes likely linked to turbulent re-acceleration of cosmic-ray electrons. However, only some clusters with double relics also show a radio halo. Here, we present a novel method to derive upper limits on th… ▽ More Pairs of radio relics are believed to form during cluster mergers, and are best observed when the merger occurs in the plane of the sky. Mergers can also produce radio halos, through complex processes likely linked to turbulent re-acceleration of cosmic-ray electrons. However, only some clusters with double relics also show a radio halo. Here, we present a novel method to derive upper limits on the radio halo emission, and analyse archival X-ray Chandra data, as well as galaxy velocity dispersions and lensing data, in order to understand the key parameter that switches on radio halo emission. We place upper limits on the halo power below the $P_{\rm 1.4 \, GHz}\, M_{500}$ correlation for some clusters, confirming that clusters with double relics have different radio properties. Computing X-ray morphological indicators, we find that clusters with double relics are associated with the most disturbed clusters. We also investigate the role of different mass-ratios and time-since-merger. Data do not indicate that the merger mass ratio has an impact on the presence or absence of radio halos (the null hypothesis that the clusters belong to the same group cannot be rejected). However, the data suggests that the absence of radio halos could be associated with early and late mergers, but the sample is too small to perform a statistical test. Our study is limited by the small number of clusters with double relics. Future surveys with LOFAR, ASKAP, MeerKat and SKA will provide larger samples to better address this issue. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: 12 pages, 7 figures, MNRAS accepted

MC$^2$: Multi-wavelength and dynamical analysis of the merging galaxy cluster ZwCl 0008.8+5215: An older and less massive Bullet Cluster

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

Abstract: We analyze a rich dataset including Subaru/SuprimeCam, HST/ACS and WFC3, Keck/DEIMOS, Chandra/ACIS-I, and JVLA/C and D array for the merging galaxy cluster ZwCl 0008.8+5215. With a joint Subaru/HST weak gravitational lensing analysis, we identify two dominant subclusters and estimate the masses to be M$_{200}=\text{5.7}^{+\text{2.8}}_{-\text{1.8}}\times\text{10}^{\text{14}}\,\text{M}_{\odot}$ and… ▽ More We analyze a rich dataset including Subaru/SuprimeCam, HST/ACS and WFC3, Keck/DEIMOS, Chandra/ACIS-I, and JVLA/C and D array for the merging galaxy cluster ZwCl 0008.8+5215. With a joint Subaru/HST weak gravitational lensing analysis, we identify two dominant subclusters and estimate the masses to be M$_{200}=\text{5.7}^{+\text{2.8}}_{-\text{1.8}}\times\text{10}^{\text{14}}\,\text{M}_{\odot}$ and 1.2$^{+\text{1.4}}_{-\text{0.6}}\times10^{14}$ M$_{\odot}$. We estimate the projected separation between the two subclusters to be 924$^{+\text{243}}_{-\text{206}}$ kpc. We perform a clustering analysis on confirmed cluster member galaxies and estimate the line of sight velocity difference between the two subclusters to be 92$\pm$164 km s$^{-\text{1}}$. We further motivate, discuss, and analyze the merger scenario through an analysis of the 42 ks of Chandra/ACIS-I and JVLA/C and D polarization data. The X-ray surface brightness profile reveals a remnant core reminiscent of the Bullet Cluster. The X-ray luminosity in the 0.5-7.0 keV band is 1.7$\pm$0.1$\times$10$^{\text{44}}$ erg s$^{-\text{1}}$ and the X-ray temperature is 4.90$\pm$0.13 keV. The radio relics are polarized up to 40$\%$. We implement a Monte Carlo dynamical analysis and estimate the merger velocity at pericenter to be 1800$^{+\text{400}}_{-\text{300}}$ km s$^{-\text{1}}$. ZwCl 0008.8+5215 is a low-mass version of the Bullet Cluster and therefore may prove useful in testing alternative models of dark matter. We do not find significant offsets between dark matter and galaxies, as the uncertainties are large with the current lensing data. Furthermore, in the east, the BCG is offset from other luminous cluster galaxies, which poses a puzzle for defining dark matter -- galaxy offsets. △ Less

Submitted 14 March, 2017; originally announced March 2017.

Comments: 22 pages, 19 figures, accepted for publication in the Astrophysical Journal on March 13, 2017

The Mismeasure of Mergers: Revised Limits on Self-interacting Dark Matter in Merging Galaxy Clusters

Authors: David Wittman, Nathan Golovich, William A. Dawson

Abstract: In an influential recent paper, Harvey et al (2015) derive an upper limit to the self-interaction cross section of dark matter ($σ_{\rm DM} < 0.47$ cm$^2$/g at 95\% confidence) by averaging the dark matter-galaxy offsets in a sample of merging galaxy clusters. Using much more comprehensive data on the same clusters, we identify several substantial errors in their offset measurements. Correcting th… ▽ More In an influential recent paper, Harvey et al (2015) derive an upper limit to the self-interaction cross section of dark matter ($σ_{\rm DM} < 0.47$ cm$^2$/g at 95\% confidence) by averaging the dark matter-galaxy offsets in a sample of merging galaxy clusters. Using much more comprehensive data on the same clusters, we identify several substantial errors in their offset measurements. Correcting these errors relaxes the upper limit on $σ_{\rm DM}$ to $\lesssim 2$ cm$^2$/g, following the Harvey et al prescription for relating offsets to cross sections in a simple solid body scattering model. Furthermore, many clusters in the sample violate the assumptions behind this prescription, so even this revised upper limit should be used with caution. Although this particular sample does not tightly constrain self-interacting dark matter models when analyzed this way, we discuss how merger ensembles may be used more effectively in the future. We conclude that errors inherent in using single-band imaging to identify mass and light peaks do not necessarily average out in a sample of this size, particularly when a handful of substructures constitute a majority of the weight in the ensemble. △ Less

Submitted 11 December, 2018; v1 submitted 20 January, 2017; originally announced January 2017.

Comments: this version accepted to ApJ; various improvements in presentation relative to v1

Journal ref: ApJ 869, 104 (2018)

The Case for Electron Re-Acceleration at Galaxy Cluster Shocks

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

Abstract: On the largest scales, the Universe consists of voids and filaments making up the cosmic web. Galaxy clusters are located at the knots in this web, at the intersection of filaments. Clusters grow through accretion from these large-scale filaments and by mergers with other clusters and groups. In a growing number of galaxy clusters, elongated Mpc-size radio sources have been found, so-called radio… ▽ More On the largest scales, the Universe consists of voids and filaments making up the cosmic web. Galaxy clusters are located at the knots in this web, at the intersection of filaments. Clusters grow through accretion from these large-scale filaments and by mergers with other clusters and groups. In a growing number of galaxy clusters, elongated Mpc-size radio sources have been found, so-called radio relics. These relics are thought to trace relativistic electrons in the intracluster plasma accelerated by low-Mach number collisionless shocks generated by cluster-cluster merger events. A long-standing problem is how low-Mach number shocks can accelerate electrons so efficiently to explain the observed radio relics. Here we report on the discovery of a direct connection between a radio relic and a radio galaxy in the merging galaxy cluster Abell 3411-3412. This discovery indicates that fossil relativistic electrons from active galactic nuclei are re-accelerated at cluster shocks. It also implies that radio galaxies play an important role in governing the non-thermal component of the intracluster medium in merging clusters. △ Less

Submitted 9 January, 2017; v1 submitted 5 January, 2017; originally announced January 2017.

Comments: Published in Nature Astronomy on Jan 4, 2017

MC$^2$: Dynamical Analysis of the Merging Galaxy Cluster MACS J1149.5+2223

Authors: Nathan Golovich, William A. Dawson, David Wittman, Georgiana A. Ogrean, Reinout J. van Weeren, Annalisa Bonafede

Abstract: We present an analysis of the merging cluster MACS J1149.5+2223 using archival imaging from Subaru/Suprime-Cam and multi-object spectroscopy from Keck/DEIMOS and Gemini/GMOS. We employ two and three dimensional substructure tests and determine that MACS J1149.5+2223 is composed of two separate mergers between three subclusters occurring $\sim$1 Gyr apart. The primary merger gives rise to elongated… ▽ More We present an analysis of the merging cluster MACS J1149.5+2223 using archival imaging from Subaru/Suprime-Cam and multi-object spectroscopy from Keck/DEIMOS and Gemini/GMOS. We employ two and three dimensional substructure tests and determine that MACS J1149.5+2223 is composed of two separate mergers between three subclusters occurring $\sim$1 Gyr apart. The primary merger gives rise to elongated X-ray morphology and a radio relic in the southeast. The brightest cluster galaxy is a member of the northern subcluster of the primary merger. This subcluster is very massive (16.7$^{+\text{1.25}}_{-\text{1.60}}\times\text{10}^{\text{14}}$ M$_{\odot}$). The southern subcluster is also very massive (10.8$^{+\text{3.37}}_{-\text{3.54}}\times\text{10}^{\text{14}}$ M$_{\odot}$), yet it lacks an associated X-ray surface brightness peak, and it has been unidentified previously despite the detailed study of this \emph{Frontier Field} cluster. A secondary merger is occurring in the north along the line of sight with a third, less massive, subcluster (1.20$^{+\text{0.19}}_{-\text{0.34}}\times\text{10}^{\text{14}}$ M$_{\odot}$). We perform a Monte Carlo dynamical analysis on the main merger and estimate a collision speed at pericenter of 2770$^{+\text{610}}_{-\text{310}}$ km s$^{-\text{1}}$. We show the merger to be returning from apocenter with core passage occurring 1.16$^{+\text{0.50}}_{-\text{0.25}}$ Gyr before the observed state. We identify the line of sight merging subcluster in a strong lensing analysis in the literature and show that it is likely bound to MACS J1149 despite having reached an extreme collision velocity of $\sim$4000 km s$^{-\text{1}}$. △ Less

Submitted 3 August, 2016; originally announced August 2016.

Comments: 17 pages, 12 figures

Frontier Fields Clusters: Deep Chandra Observations of the Complex Merger MACS J1149.6+2223

Authors: G. A. Ogrean, R. J. van Weeren, C. Jones, W. Forman, W. A. Dawson, N. Golovich, F. Andrade-Santos, S. S. Murray, P. Nulsen, E. Roediger, A. Zitrin, E. Bulbul, R. Kraft, A. Goulding, K. Umetsu, T. Mroczkowski, A. Bonafede, S. Randall, J. Sayers, E. Churazov, L. David, J. Merten, M. Donahue, B. Mason, P. Rosati , et al. (2 additional authors not shown)

Abstract: The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock precedin… ▽ More The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock preceding the cold front. The substructure analysis of the cluster identified several components with large relative radial velocities, thus indicating that at least some collisions occur almost along the line of sight. The inclination of the mergers with respect to the plane of the sky poses significant observational challenges at X-ray wavelengths. MACS J1149.6+2223 possibly hosts a steep-spectrum radio halo. If the steepness of the radio halo is confirmed, then the radio spectrum, combined with the relatively regular ICM morphology, could indicate that MACS J1149.6+2223 is an old merging cluster. △ Less

Submitted 18 March, 2016; originally announced March 2016.

Comments: Accepted to ApJ

A powerful double radio relic system discovered in PSZ1 G108.18-11.53: evidence for a shock with non-uniform Mach number?

Authors: F. de Gasperin, H. T. Intema, R. J. van Weeren, W. A. Dawson, N. Golovich, D. Wittman, A. Bonafede, M. Bruggen

Abstract: Diffuse radio emission in the form of radio halos and relics has been found in a number of merging galaxy clusters. These structures indicate that shock and turbulence associated with the merger accelerate electrons to relativistic energies. We report the discovery of a radio relic + radio halo system in PSZ1 G108.18-11.53 (z=0.335). This cluster hosts the second most powerful double radio relic s… ▽ More Diffuse radio emission in the form of radio halos and relics has been found in a number of merging galaxy clusters. These structures indicate that shock and turbulence associated with the merger accelerate electrons to relativistic energies. We report the discovery of a radio relic + radio halo system in PSZ1 G108.18-11.53 (z=0.335). This cluster hosts the second most powerful double radio relic system ever discovered. We observed PSZ1 G108.18-11.53 with the Giant Meterwave Radio Telescope (GMRT) and the Westerbork Synthesis Radio Telescope (WSRT). We obtained radio maps at 147, 323, 607 and 1380 MHz. We also observed the cluster with the Keck telescope, obtaining the spectroscopic redshift for 42 cluster members. From the injection index we obtained the Mach number of the shocks generating the two radio relics. For the southern shock we found M = 2.33^{+0.19}_{-0.26}, while the northern shock Mach number goes from M = 2.20^{+0.07}_{-0.14} in the north part down to M = 2.00^{+0.03}_{-0.08} in the southern region. If the relation between the injection index and the Mach number predicted by diffusive shock acceleration (DSA) theory holds, this is the first observational evidence for a gradient in the Mach number along a galaxy cluster merger shock. △ Less

Submitted 12 August, 2015; originally announced August 2015.

Comments: 16 pages, 17 images, accepted MNRAS

MC^2: Galaxy Imaging and Redshift Analysis of the Merging Cluster CIZA J2242.8+5301

Authors: William A. Dawson, M. James Jee, Andra Stroe, Y. Karen Ng, Nathan Golovich, David Wittman, David Sobral, Marcus Bruggen, H. J. A. Rottgering, R. J. van Weeren

Abstract: X-ray and radio observations of CIZA J2242.8+5301 suggest that it is a major cluster merger. Despite being well studied in the X-ray, and radio, little has been presented on the cluster structure and dynamics inferred from its galaxy population. We carried out a deep (i<25) broad band imaging survey of the system with Subaru SuprimeCam (g & i bands) and the Canada France Hawaii Telescope (r band)… ▽ More X-ray and radio observations of CIZA J2242.8+5301 suggest that it is a major cluster merger. Despite being well studied in the X-ray, and radio, little has been presented on the cluster structure and dynamics inferred from its galaxy population. We carried out a deep (i<25) broad band imaging survey of the system with Subaru SuprimeCam (g & i bands) and the Canada France Hawaii Telescope (r band) as well as a comprehensive spectroscopic survey of the cluster area (505 redshifts) using Keck DEIMOS. We use this data to perform a comprehensive galaxy/redshift analysis of the system, which is the first step to a proper understanding the geometry and dynamics of the merger, as well as using the merger to constrain self-interacting dark matter. We find that the system is dominated by two subclusters of comparable richness with a projected separation of 6.9'^{+0.7}_{-0.5} (1.3^{+0.13}_{-0.10} Mpc). We find that the north and south subclusters have similar redshifts of z=0.188 with a relative line-of-sight velocity difference of 69+/-190 km/s. We also find that north and south subclusters have velocity dispersions of 1160^{+100}_{-90} km/s and 1080^{+100}_{-70} km/s, respectively. These correspond to masses of 16.1^{+4.6}_{-3.3}x10^14 M_sun and 13.0^{+4.0}_{-2.5}x10^14 M_sun, respectively. While velocity dispersion measurements of merging clusters can be biased we believe the bias in this system to be minor due to the large projected separation and nearly plane-of-sky merger configuration. CIZA J2242.8+5301 is a relatively clean dissociative cluster merger with near 1:1 mass ratio, which makes it an ideal merger for studying merger associated physical phenomena. △ Less

Submitted 10 October, 2014; originally announced October 2014.

Comments: 17 pages, 14 figures, submitted to ApJ

Report number: LLNL-DOC-661846-DRAFT