In an influential recent paper, Harvey et al. derive an upper limit to the self-interaction cross section of dark matter (DM) ( DM m < 0.47 cm2 g−1 at 95% confidence) by averaging the DM-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 DM m to 2 cm2 g−1, 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 DM 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.
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 supernova 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 (tE > 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 tE > 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 criterion in order to statistically constrain the abundance of BHs in the Milky Way. The future Wide Field Infrared Survey Telescope (WFIRST) microlensing survey provides both precise photometry and astrometry and will yield individual masses of BHs, 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 BH present-day mass function, BH multiplicity, and BH kick velocity distributions.
Abstract
Using
N
-body simulation results from the MEGASIM data set, we present spatial distributions of Earth Trojan Asteroids and assess the detectability of the population in current and ...next-generation ground-based astronomical surveys. 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 it 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 m (at
H
≈ 22), while a null detection by LSST will further restrict the ETA population to tens of objects larger than 100 m.
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 the 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
α
emission line in ∼100 member galaxies in each merger, classify galaxies as emitters or nonemitters, and then classify emitters as star-forming galaxies (SFGs) or active galactic nucleus (AGN) based on the N
ii
λ
6583 line. We quantify the distribution of SFG and AGN relative to nonemitters 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 nonemitters 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.
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 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 Optical Gravitational Lensing Experiment-IV (OGLE-IV) and Nancy Grace Roman Space Telescope (Roman) microlensing surveys. We provide a proof-of-concept analysis for adding PBHs into PopSyCLE , and thus include many simplifying assumptions, such as f DM , the fraction of dark matter composed of PBHs, and m ¯ PBH , mean PBH mass. Assuming m ¯ PBH = 30 M ⊙ , we find ∼3.6 f DM times as many PBH microlensing events than stellar evolved black hole events, a PBH average peak Einstein crossing time of ∼91.5 days, estimate on order of 10 2 f DM PBH events within the 8 yr OGLE-IV results, and estimate Roman to detect ∼1000 f DM PBH microlensing events throughout its planned microlensing survey.
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 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 MACS J1752.0+4440 (
z
= 0.365) and ZWCL 1856.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 MACS J1752.0+4440 (ZWCL 1856.8+6616) is
M
200
=
14.7
−
3.3
+
3.8
×
10
14
M
⊙
(
M
200
=
2.4
−
0.7
+
0.9
×
10
14
M
⊙
). 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 ZWCL 1856.8+6618 is likely at a later merger phase than MACS J1752.0+4440. We postulate that the SW radio relic in MACS J1752.0+4440 is a result of particle reacceleration.
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, 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
⊙
. 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 ≈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.
Abstract
Modern surveys of gravitational microlensing events have progressed to detecting thousands per year, and 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 the microlensing Einstein radius crossing time (
t
E
) distribution. However, systematics in individual light curves as well as oversimplistic 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
t
E
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 the previous analyses overestimated the number of
t
E
> 100 day events due to their oversimplistic model ignoring parallax effects. We use our fits to identify the hundreds filling a regime in the microlensing parameter space that are 50% pure of black holes. Finally, we have released the largest-ever catalog of Markov Chain Monte Carlo parameter estimates for microlensing events.
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, including the Galactic plane, every few nights. The ZTF observes ∼109 stars in the g and r bands and can significantly contribute to the observed microlensing population. We predict that the ZTF will observe ∼1100 microlensing events in 3 yr of observing within 10° latitude of the Galactic plane, with ∼500 events in the outer Galaxy ( ≥ 10°). This yield increases to ∼1400 (∼800) events by combining every three ZTF exposures, ∼1800 (∼900) events if the ZTF observes for a total of 5 yr, and ∼2400 (∼1300) events for a 5 yr survey with postprocessing 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. The ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing.
Abstract
We present an analysis of lifetimes and resonances of Earth Trojan Asteroids (ETAs) in the MEGASIM data set. 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 TK
7
and 2020 XL
5
), several Mars Trojans, and a Venus Trojan discovered. The estimated lifetimes of the known inner system Trojans are shorter 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 that 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.