Engineering topics

McGill, Peter (ORCID:0000000210526749)

Publications and source records attributed to McGill, Peter (ORCID:0000000210526749).

5 records

DOE OSTI2024

Primordial Black Hole Dark Matter Simulations Using PopSyCLE

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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Astrometric Microlensing by Primordial Black Holes with the Roman Space Telescope

Abstract Primordial black holes (PBHs) could explain some fraction of dark matter and shed light on many areas of early-Universe physics. Despite over half a century of research interest, a PBH population has so far eluded detection. The most competitive constraints on the fraction of dark matter comprised of PBHs ( f DM ) in the (10 −9 –10) M ⊙ mass ranges come from photometric microlensing and bound f DM ≲ 10 −2 –10 −1 . With the advent of the Roman Space Telescope with its submilliarcsecond astrometric capabilities and its planned Galactic Bulge Time Domain Survey (GBTDS), detecting astrometric microlensing signatures will become routine. Compared with photometric microlensing, astrometric microlensing signals are sensitive to different lens masses–distance configurations and contain different information, making it a complimentary lensing probe. At submilliarcsecond astrometric precision, astrometric microlensing signals are typically detectable at larger lens–source separations than photometric signals, suggesting a microlensing detection channel of pure astrometric events. We use a Galactic simulation to predict the number of detectable microlensing events during the GBTDS via this pure astrometric microlensing channel. Assuming an absolute astrometric precision floor for bright stars of 0.1 mas for the GBTDS, we find that the number of detectable events peaks at ≈10 3 f DM for a population of 1 M ⊙ PBHs and tapers to ≈10 f DM and ≈100 f DM at 10 −4 M ⊙ and 10 3 M ⊙ , respectively. Accounting for the distinguishability of PBHs from stellar lenses, we conclude the GBTDS will be sensitive to a PBH population at f DM down to ≈10 −1 –10 −3 for (10 −1 –10 2 ) M ⊙ likely yielding novel PBH constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Spatially resolved microlensing time-scale distributions across the Galactic bulge with the VVV survey

ABSTRACT We analyse 1602 microlensing events found in the VISTA Variables in the Via Lactea (VVV) near-infrared (NIR) survey data. We obtain spatially resolved, efficiency-corrected time-scale distributions across the Galactic bulge (|ℓ| < 10°, |b| < 5°), using a Bayesian hierarchical model. Spatially resolved peaks and means of the time-scale 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 time-scales in the central bulge fields (|ℓ| < 5°) are on average shorter than the non-central (|ℓ| > 5°) fields, with the average peak of the lognormal time-scale distribution at 23.6 ± 1.9 d for the central fields and 29.0 ± 3.0 d 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 time-scale distributions. We propose two pairs of Roman fields, centred on (ℓ = ±9, 5°, b = −0.125°) and (ℓ = −5°, b = ±1.375°) as good targets to measure the asymmetry in longitude and latitude, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Disentangling the Black Hole Mass Spectrum with Photometric Microlensing Surveys

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.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Bringing 2D Eclipse Mapping out of the Shadows with Leave-one-out Cross Validation

Abstract Eclipse mapping is a technique for inferring 2D brightness maps of transiting exoplanets from the shape of an eclipse light curve. With JWST’s unmatched precision, eclipse mapping is now possible for a large number of exoplanets. However, eclipse mapping has only been applied to two planets, and the nuances of fitting eclipse maps are not yet fully understood. Here, we use Leave-one-out Cross Validation (LOO-CV) to investigate eclipse mapping, with application to a JWST NIRISS/SOSS observation of the ultrahot Jupiter WASP-18b. LOO-CV is a technique that provides insight into the out-of-sample predictive power of models on a data-point-by-data-point basis. We show that constraints on planetary brightness patterns behave as expected, with large-scale variations driven by the phase-curve variation in the light curve and smaller-scale structures constrained by the eclipse ingress and egress. For WASP-18b we show that the need for higher model complexity (smaller-scale features) is driven exclusively by the shape of the eclipse ingress and egress. We use LOO-CV to investigate the relationship between planetary brightness map components when mapping under a positive-flux constraint to better understand the need for complex models. Finally, we use LOO-CV to understand the degeneracy between the competing “hot spot” and “plateau” brightness map models of WASP-18b, showing that the plateau model is driven by the ingress shape and the hot spot model is driven by the egress shape, but preference for neither model is due to outliers or unmodeled signals. Based on this analysis, we make recommendations for the use of LOO-CV in future eclipse-mapping studies.

79 ASTRONOMY AND ASTROPHYSICS↗