Search for: All records

Abstract This review is focused on tests of Einstein’s theory of general relativity with gravitational waves that are detectable by ground-based interferometers and pulsar-timing experiments. Einstein’s theory has been greatly constrained in the quasi-linear, quasi-stationary regime, where gravity is weak and velocities are small. Gravitational waves are allowing us to probe a complimentary, yet previously unexplored regime: the non-linear and dynamicalextreme gravity regime. Such a regime is, for example, applicable to compact binaries coalescing, where characteristic velocities can reach fifty percent the speed of light and gravitational fields are large and dynamical. This review begins with the theoretical basis and the predicted gravitational-wave observables of modified gravity theories. The review continues with a brief description of the detectors, including both gravitational-wave interferometers and pulsar-timing arrays, leading to a discussion of the data analysis formalism that is applicable for such tests. The review then discusses gravitational-wave tests using compact binary systems, and ends with a description of the first gravitational wave observations by advanced LIGO, the stochastic gravitational wave background observations by pulsar timing arrays, and the tests that can be performed with them. 

Abstract While supermassive black hole (SMBH) binaries are not the only viable source for the low-frequency gravitational wave background (GWB) signal evidenced by the most recent pulsar timing array (PTA) data sets, they are expected to be the most likely. Thus, connecting the measured PTA GWB spectrum and the underlying physics governing the demographics and dynamics of SMBH binaries is extremely important. Previously, Gaussian processes (GPs) and dense neural networks have been used to make such a connection by being built as conditional emulators; their input is some selected evolution or environmental SMBH binary parameters and their output is the emulated mean and standard deviation of the GWB strain ensemble distribution over many Universes. In this paper, we use a normalizing flow (NF) emulator that is trained on the entirety of the GWB strain ensemble distribution, rather than only mean and standard deviation. As a result, we can predict strain distributions that mirror underlying simulations very closely while also capturing frequency covariances in the strain distributions as well as statistical complexities such as tails, non-Gaussianities, and multimodalities that are otherwise not learnable by existing techniques. In particular, we feature various comparisons between the NF-based emulator and the GP approach used extensively in past efforts. Our analyses conclude that the NF-based emulator not only outperforms GPs in the ease and computational cost of training but also outperforms in the fidelity of the emulated GWB strain ensemble distributions. 

Abstract Massive black hole binaries (MBHBs) produce gravitational waves (GWs) that are detectable with pulsar timing arrays. We determine the properties of the host galaxies of simulated MBHBs at the time they are producing detectable GW signals. The population of MBHB systems we evaluate is from theIllustriscosmological simulations taken in tandem with post processing semi-analytic models of environmental factors in the evolution of binaries. Upon evolving to the GW frequency regime accessible by pulsar timing arrays, we calculate the detection probability of each system using a variety of different values for pulsar noise characteristics in a plausible near-future International Pulsar Timing Array dataset. We find that detectable systems have host galaxies that are clearly distinct from the overall binary population and from most galaxies in general. With conservative noise factors, we find that host stellar metallicity, for example, peaks at $\sim 2Z_{\odot }$as opposed to the total population of galaxies which peaks at $\sim 0.6Z_{\odot }$. Additionally, the most detectable systems are much brighter in magnitude and more red in color than the overall population, indicating their likely identity as large ellipticals with diminished star formation. These results can be used to develop effective search strategies for identifying host galaxies and electromagnetic counterparts following GW detection by pulsar timing arrays. 

Abstract Pulsar timing arrays (PTAs) are sensitive to low-frequency gravitational waves (GWs), which induce correlated changes in millisecond pulsars’ timing residuals. PTA collaborations around the world have recently announced evidence of a nanohertz gravitational wave background (GWB), which may be produced by a population of supermassive black hole binaries (SMBHBs). The GWB is often modeled as following a power-law power spectral density (PSD); however, a GWB produced by a cosmological population of SMBHBs is expected to have a more complex power spectrum due to the discrete nature of the sources. In this paper, we investigate using at-process PSD to model the GWB, which allows us to fit for both the underlying power-law amplitude and spectral index as well as deviations from that power law, which may be produced by individual nearby binaries. We create simulated data sets based on the properties of the NANOGrav 15 yr data set, and we demonstrate that thet-process PSD can accurately recover the PSD when deviations from a power law are present. With longer timed data sets and more pulsars, we expect the sensitivity of our PTAs to improve, which will allow us to precisely measure the PSD of the GWB and study the sources producing it. 

Abstract Low-frequency solar radio emission is sourced in the solar corona, with sub-100 MHz radio emission largely originating from the ∼10⁵K plasma around 2 optical radii. However, the region of emission has yet to be constrained at 35–45 MHz due to both instrumentation limitations and the rarity of astronomical events, such as total solar eclipses, which allow for direct observational approaches. In this work, we present the results from a student-led project to commission a low-frequency radio telescope array situated in the path of totality of the 2024 total solar eclipse in an effort to probe the middle corona. The Deployable Low-Band Ionosphere and Transient Experiment (DLITE) is a low-frequency radio array comprised of four dipole antennas, optimized to observe at 35–45 MHz, and capable of resolving the brightest radio sources in the sky. We constructed a DLITE station in Observatory Park, a dark-sky park in Montville, Ohio. Results of observations during the total solar eclipse demonstrate that DLITE stations can be quickly deployed for observations and provide constraints on the radius of solar emission at our center observing frequency of 42 MHz. In this work, we outline the construction of DLITE Ohio and the solar observation results from the total solar eclipse that transversed North America in 2024 April. 

Abstract Finding and characterizing gravitational waves from individual supermassive black hole binaries is a central goal of pulsar timing array experiments, which will require analysis methods that can be efficient on our rapidly growing datasets. Here we present a novel approach built on three key elements: (i) precalculating and interpolating expensive matrix operations; (ii) semi-analytically marginalizing over the gravitational-wave phase at the pulsars; (iii) numerically marginalizing over the pulsar distance uncertainties. With these improvements the recent NANOGrav 15 yr dataset can be analyzed in minutes after an $O\left(1\mathrm{h}\right)$setup phase, instead of an analysis taking days–weeks with previous methods. The same setup can be used to efficiently analyze the dataset under any sinusoidal deterministic model. In particular, this will aid testing the binary hypothesis by allowing for efficient analysis of competing models (e.g. incoherent, monopolar, or dipolar sine wave model) and scrambled datasets for false alarm studies. The same setup can be updated in minutes for new realizations of the data, which enables large simulation studies. 

Abstract A pulsar’s scintillation bandwidth is inversely proportional to the scattering delay, making accurate measurements of scintillation bandwidth critical to characterize unmitigated delays in efforts to measure low-frequency gravitational waves with pulsar timing arrays. In this pilot work, we searched for a subset of known pulsars within ∼97% of the data taken with the Puerto Rico Ultimate Pulsar Processing Instrument for the AO327 survey with the Arecibo telescope, attempting to measure the scintillation bandwidths in the data set by fitting to the 2D autocorrelation function of their dynamic spectra. We successfully measured 38 bandwidths from 23 pulsars (six without prior literature values), finding that: almost all of the measurements are larger than the predictions from NE2001 and YMW16 (two popular galactic models); NE2001 is more consistent with our measurements than YMW16; Gaussian fits to the bandwidth are more consistent with both electron density models than Lorentzian ones; and for the 17 pulsars with prior literature values, the measurements between various sources often vary by a few factors. The success of Gaussian fits may be due to the use of Gaussian fits to train models in previous work. The variance of literature values over time could relate to the scaling factor used to compare measurements, but also seems consistent with time-varying interstellar medium parameters. This work can be extended to the rest of AO327 to further investigate these trends, highlighting the continuing importance of large archival data sets for projects beyond their initial conception. 

Abstract We report on findings from scintillation analyses using high-cadence observations of eight canonical pulsars with observing baselines ranging from 1–3 yr. We obtain scintillation bandwidth and timescale measurements for all pulsars in our survey and scintillation arc curvature measurements for four, and we detect multiple arcs for two. We find evidence of a previously undocumented scattering screen along the line of sight (LOS) to PSR J1645−0317, as well as evidence that a scattering screen along the LOS to PSR J2313+4253 may reside somewhere within the Milky Way’s Orion–Cygnus arm. We report evidence of a significant change in the scintillation pattern in PSR J2022+5154 from the previous two decades of literature, wherein both the scintillation bandwidth and timescale decreased by an order of magnitude relative to earlier observations at the same frequencies, potentially as a result of a different screen dominating the observed scattering. By augmenting the results of previous studies, we find general agreement with estimations of scattering delays from pulsar observations and predictions by the NE2001 electron density model but not for the newest data we have collected, providing some evidence of changes in the ISM along various LOSs over the timespans considered. In a similar manner, we find additional evidence of a correlation between a pulsar’s dispersion measure and the overall variability of its scattering delays over time. The plethora of interesting science obtained through these observations demonstrates the capabilities of the Green Bank Observatory’s 20 m telescope to contribute to pulsar-based studies of the interstellar medium. 

Abstract Nonradiative shocks accelerate particles and heat astrophysical plasmas. While supernova remnants are the most well-studied example, neutron star (NS) bow shocks are also nonradiative and Balmer dominated. NS bow shocks are likely ubiquitous in the interstellar medium due to their large speeds imparted at birth, and they are thought to be a discrete source population contributing to the Galactic cosmic-ray spectrum. To date, nine NS bow shocks have been directly observed in Hαimages. Most of these shocks have been characterized using narrowband Hαimaging and slit spectroscopy, which do not resolve the multicomponent velocity structure of the shocks and their spatial geometry. Here we present integral field spectroscopy of three NS bow shocks: J0742−2822, J1741−2054, and J2225+6535 (the Guitar Nebula). We measure the shock properties simultaneously in four dimensions: the 2D projected shock morphology, the radial velocity structure, and the Hαflux. The broad-to-narrow line ratio (I_b/I_n) is inferred from radial velocity profiles, and for J1741−2054, the narrow line is detected in multiple regions of the shock. The inferred line ratios and widths suggest that NS bow shocks represent a low-shock velocity regime (V≲ 200 km s⁻¹) in whichI_b/I_nis high, distinct from the shock regime probed by supernova remnants. Our results illustrate a need for nonradiative shock models at velocities lower than previously considered, which will reveal the electron–ion temperature ratios and particle acceleration efficiencies of these bow shocks. 

Abstract Supermassive black hole binaries (SMBHBs) are thought to form in galaxy mergers, possessing the potential to produce electromagnetic (EM) radiation as well as gravitational waves (GWs) detectable with pulsar timing arrays (PTAs). Once GWs from individually resolved SMBHBs are detected, the identification of the host galaxy will be a major challenge due to the ambiguity in possible EM signatures and the poor localization capability of PTAs. To aid EM observations in choosing follow-up sources, we use NANOGrav’s galaxy catalog to quantify the number of plausible hosts in both realistic and idealistic scenarios. We outline a host identification pipeline that injects a single-source GW signal into a simulated PTA data set, recovers the signal using production-level techniques, quantifies the localization region and number of galaxies contained therein, and finally imposes cuts on the galaxies using parameter estimates from the GW search. In an ideal case, the 90% credible areas span 29–241 deg², containing about 14–341 galaxies. After cuts, the number of galaxies remaining ranges from 22 at worst to one true host at best. In a realistic case, these areas range from 287 to 530 deg²and enclose about 285–1238 galaxies. After cuts, the number of galaxies is 397 at worst and 27 at best. While the signal-to-noise ratio is the primary determinant of the localization area of a given source, we find that the area is also influenced by the proximity to nearby pulsars on the sky and the binary chirp mass.