Modern large-scale physics experiments create datasets with sizes and streaming rates that can exceed those from industry leaders such as Google Cloud and Netflix. Fully processing these datasets ...requires both sufficient compute power and efficient workflows. Recent advances in Machine Learning (ML) and Artificial Intelligence (AI) can either improve or replace existing domain-specific algorithms to increase workflow efficiency. Not only can these algorithms improve the physics performance of current algorithms, but they can often be executed more quickly, especially when run on coprocessors such as GPUs or FPGAs. In the winter of 2023, MIT hosted the Accelerating Physics with ML at MIT workshop, which brought together researchers from gravitational-wave physics, multi-messenger astrophysics, and particle physics to discuss and share current efforts to integrate ML tools into their workflows. The following white paper highlights examples of algorithms and computing frameworks discussed during this workshop and summarizes the expected computing needs for the immediate future of the involved fields.
We present new proper motion measurements of optically emitting oxygen-rich knots of supernova remnant 1E 0102.2-7219 (E0102), which are used to estimate the remnant's center of expansion and age. ...Four epochs of high resolution Hubble Space Telescope images spanning 19 yr were retrieved and analyzed. We found a robust center of expansion of alpha=1:04:02.48 and delta=-72:01:53.92 (J2000) with 1-sigma uncertainty of 1.77 arcseconds using 45 knots from images obtained with the Advanced Camera for Surveys using the F475W filter in 2003 and 2013 having the highest signal-to-noise ratio. We also estimate an upper limit explosion age of 1738 +/- 175 yr by selecting knots with the highest proper motions, that are assumed to be the least decelerated. We find evidence of an asymmetry in the proper motions of the knots as a function of position angle. We conclude that these asymmetries were most likely caused by interaction between E0102's original supernova blast wave and an inhomogeneous surrounding environment, as opposed to intrinsic explosion asymmetry. The observed non-homologous expansion suggests that the use of a free expansion model inaccurately offsets the center of expansion and leads to an overestimated explosion age. We discuss our findings as they compare to previous age and center of expansion estimates of E0102 and their relevance to a recently identified candidate central compact object.
Narrow transient emission lines (flash-ionization features) in early supernova (SN) spectra trace the presence of circumstellar material (CSM) around the massive progenitor stars of core-collapse ...SNe. The lines disappear within days after the SN explosion, suggesting that this material is spatially confined, and originates from enhanced mass loss shortly (months to a few years) prior to explosion. We performed a systematic survey of H-rich (Type II) SNe discovered within less than two days from explosion during the first phase of the Zwicky Transient Facility (ZTF) survey (2018-2020), finding thirty events for which a first spectrum was obtained within \(< 2\) days from explosion. The measured fraction of events showing flash ionisation features (\(>36\%\) at \(95\%\) confidence level) confirms that elevated mass loss in massive stars prior to SN explosion is common. We find that SNe II showing flash ionisation features are not significantly brighter, nor bluer, nor more slowly rising than those without. This implies that CSM interaction does not contribute significantly to their early continuum emission, and that the CSM is likely optically thin. We measured the persistence duration of flash ionisation emission and find that most SNe show flash features for \(\approx 5 \) days. Rarer events, with persistence timescales \(>10\) days, are brighter and rise longer, suggesting these may be intermediate between regular SNe II and strongly-interacting SNe IIn.
We present a three-dimensional kinematic reconstruction of the optically-emitting, oxygen-rich ejecta of supernova remnant N132D in the Large Magellanic Cloud. Data were obtained with the 6.5 m ...Magellan telescope in combination with the IMACS+GISMO instrument and survey O III \(\lambda\lambda\)4959,5007 line emission in a \({\sim}\)3\(^{\prime}~\times\) 3\(^{\prime}\) region centered on N132D. The spatial and spectral resolution of our data enable detailed examination of the optical ejecta structure. The majority of N132D's optically bright oxygen ejecta are arranged in a torus-like geometry tilted approximately 28\(^{\circ}\) with respect to the plane of the sky. The torus has a radius of 4.4 pc (\(D_{\rm LMC}\)/50 kpc), exhibits a blue-shifted radial velocity asymmetry of \(-3000\) to \(+2300\) km s\(^{-1}\), and has a conspicuous break in its circumference. Assuming homologous expansion from the geometric center of O-rich filaments, the average expansion velocity of 1745 km s\(^{-1}\) translates to an age since explosion of 2450 \(\pm\) 195 yr. A faint, spatially-separated "runaway knot" (RK) with total space velocity of 3650 km s\(^{-1}\) is nearly perpendicular to the torus plane and coincident with X-ray emission that is substantially enhanced in Si relative to the LMC and N132D's bulk ejecta. These kinematic and chemical signatures suggest that the RK may have had its origin deep within the progenitor star. Overall, the main shell morphology and high-velocity, Si-enriched components of N132D have remarkable similarity with that of Cassiopeia A, which was the result of a Type IIb supernova explosion. Our results underscore the need for further observations and simulations that can robustly reconcile whether the observed morphology is dominated by explosion dynamics or shaped by interaction with the environment.
Although not nearly as numerous as binaries with two white dwarfs, eccentric neutron star-white dwarf (NS-WD) binaries are important gravitational-wave (GW) sources for the next generation of ...space-based detectors sensitive to low frequency waves. Here we investigate periastron precession in these sources as a result of general relativistic, tidal, and rotational effects; such precession is expected to be detectable for at least some of the detected binaries of this type. Currently, two eccentric NS-WD binaries are known in the galactic field, PSR J1141-6545 and PSR B2303+46, both of which have orbits too wide to be relevant in their current state to GW observations. However, population synthesis studies predict the existence of a significant Galactic population of such systems. Though small in most of these systems, we find that tidally induced periastron precession becomes important when tides contribute to more than 3% of the total precession rate. For these systems, accounting for tides when analyzing periastron precession rate measurements can improve estimates of the WD component mass inferred and, in some cases, will prevent us from misclassifying the object. However, such systems are rare due to rapid orbital decay. To aid the inclusion of tidal effects when using periastron precession as a mass measurement tool, we derive a function that relates the WD radius and periastron precession constant to the WD mass.
The Large Synoptic Survey Telescope (LSST) will open a discovery frontier for faint and fast transients with its ability to detect variable flux components down to \(\sim\)24.5 mag in a \(\sim\)30 ...second exposure. Spectroscopic follow-up of such phenomena - which are necessary for understanding the physics of stellar explosions - can require a rapid response and several hours with a 8-10m telescope, making it both expensive and difficult to acquire. The future Extremely Large Telescopes (ELTs) would be able to provide not only spectroscopy but capabilities such as spectropolarimetry and high-resolution diffraction-limited imaging that would contribute to future advances in our physical understanding of stellar explosions. In this white paper we focus on several specific scientific impacts in the field of explosive transient astrophysics that will be generated by the combination of LSST's discovery abilities and ELTs' follow-up capacities. First, we map the uncharted frontier of discovery phase-space in terms of intrinsic luminosity and timescales for explosive transients, where we expect the unexpected. We then focus on six areas with open science questions for known transients: the progenitors of thermonuclear supernovae (SNe), mass loss prior to core collapse, asymmetry in stellar explosions, light echoes, high-\(z\) transients, and strongly lensed SNe. We conclude with a brief discussion of the practical aspects of ELT & LSST synergy.
Whether supernovae are a significant source of dust has been a long-standing debate. The large quantities of dust observed in high-redshift galaxies raise a fundamental question as to the origin of ...dust in the Universe since stars cannot have evolved to the AGB dust-producing phase in high-redshift galaxies. In contrast, supernovae occur within several millions of years after the onset of star formation. This white paper focuses on dust formation in supernova ejecta with US-Extremely Large Telescope (ELT) perspective during the era of JWST and LSST.
The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered ...a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the NMMA (Nuclear-physics and Multi-Messenger Astrophysics) framework. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. It also enables us to classify electromagnetic observations, e.g., to distinguish between supernovae and kilonovae. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions. The extension of the NMMA code presented here is the first attempt of analysing the gravitational-wave signal, the kilonovae, and the GRB afterglow simultaneously, which reduces the uncertainty of our constraints. Incorporating all available information, we estimate the radius of a 1.4 solar mass neutron star to be \(R=11.98^{+0.35}_{-0.40}\) km.
The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and ...neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.
We perform a Bayesian analysis of the mass distribution of stellar-mass black holes using the observed masses of 15 low-mass X-ray binary systems undergoing Roche lobe overflow and five high-mass, ...wind-fed X-ray binary systems. Using Markov Chain Monte Carlo calculations, we model the mass distribution both parametrically---as a power law, exponential, gaussian, combination of two gaussians, or log-normal distribution---and non-parametrically---as histograms with varying numbers of bins. We provide confidence bounds on the shape of the mass distribution in the context of each model and compare the models with each other by calculating their relative Bayesian evidence as supported by the measurements, taking into account the number of degrees of freedom of each model. The mass distribution of the low-mass systems is best fit by a power-law, while the distribution of the combined sample is best fit by the exponential model. We examine the existence of a "gap" between the most massive neutron stars and the least massive black holes by considering the value, M_1%, of the 1% quantile from each black hole mass distribution as the lower bound of black hole masses. The best model (the power law) fitted to the low-mass systems has a distribution of lower-bounds with M_1% > 4.3 Msun with 90% confidence, while the best model (the exponential) fitted to all 20 systems has M_1% > 4.5 Msun with 90% confidence. We conclude that our sample of black hole masses provides strong evidence of a gap between the maximum neutron star mass and the lower bound on black hole masses. Our results on the low-mass sample are in qualitative agreement with those of Ozel, et al (2010).
