Search

Observations show that globular clusters might be among the best places to find millisecond pulsars. However, the globular cluster Terzan 6 seems to be an exception without any pulsar discovered, although its high stellar encounter rate suggests that it harbors dozens of them. We report the discovery of the first radio pulsar, PSR J1750-3116A, likely associated with Terzan 6 in a search of C-band (4-8 GHz) data from the Robert C. Byrd Green Bank Telescope with a spin period of 5.33 ms and dispersion measure, DM$\simeq$ 383 ${\rm pc\,cm^{-3}}$. The mean flux density of this pulsar is approximately 3 $\mu$Jy. The DM agrees well with predictions from the Galactic free electron density model, assuming a distance of 6.7 kpc for Terzan 6. PSR J1750-3116A is likely an isolated millisecond pulsar, potentially formed through dynamical interactions, considering the core-collapsed classification and the exceptionally high stellar encounter rate of Terzan 6. This is the highest radio frequency observation that has led to the discovery of a pulsar in a globular cluster to date. While L-band (1-2GHz) observations of this cluster are unlikely to yield significant returns due to propagation effects, we predict that further pulsar discoveries in Terzan 6 will be made by existing radio telescopes at higher frequencies., Comment: Updated to match the published version, correcting the PSR name to PSR J1750-3116A without rounding

The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index \gamma. In this paper, we consider the next logical step beyond this minimal spectral model, allowing for a running (i.e., logarithmic frequency dependence) of the spectral index, \gamma_run(f) = \gamma + \beta \ln(f/f_ref). We fit this running-power-law (RPL) model to the NANOGrav 15-year data and perform a Bayesian model comparison with the minimal constant-power-law (CPL) model, which results in a 95% credible interval for the parameter \beta consistent with no running, \beta \in [-0.80,2.96], and an inconclusive Bayes factor, B(RPL vs. CPL) = 0.69 +- 0.01. We thus conclude that, at present, the minimal CPL model still suffices to adequately describe the NANOGrav signal; however, future data sets may well lead to a measurement of nonzero \beta. Finally, we interpret the RPL model as a description of primordial GWs generated during cosmic inflation, which allows us to combine our results with upper limits from big-bang nucleosynthesis, the cosmic microwave background, and LIGO-Virgo-KAGRA., Comment: 17 pages, 4 figures, 2 tables

Radio-frequency interference (RFI) is becoming an increasingly significant problem for most radio telescopes. Working with Green Bank Telescope data from PSR J1730+0747 in the form of complex-valued channelized voltages and their respective high-resolution power spectral densities, we evaluate a variety of statistical measures to characterize RFI. As a baseline for performance comparison, we use median absolute deviation (MAD) in complex channelized voltage data and spectral kurtosis (SK) in power spectral density data to characterize and filter out RFI. From a new perspective, we implement the Shapiro-Wilks (SW) test for normality and two information theoretical measures, spectral entropy (SE) and spectral relative entropy (SRE), and apply them to mitigate RFI. The baseline RFI mitigation algorithms are compared against our novel RFI detection algorithms to determine how effective and robust the performance is. Except for MAD, we find significant improvements in signal-to-noise ratio through the application of SE, symmetrical SRE, asymmetrical SRE, SK, and SW. These algorithms also do a good job of characterizing broadband RFI. Time- and frequency-variable RFI signals are best detected by SK and SW tests., Comment: 13 pages, 26 figures, accepted for publication in RAS Techniques and Instruments

Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions are idealizations that may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks: after fitting our fiducial model to real data, we generate a population of simulated data-set replications, and use them to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients assume extreme values for the real data when compared to the replications. We confirm that the NANOGrav 15 yr data set is consistent with power-law and Hellings--Downs assumptions. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2\ $\sigma$ and 3\ $\sigma$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background as reported by NANOGrav., Comment: 20 pages, 14 Figures

Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4\sigma$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in \textit{Gaia} DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6\%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 80 kyr, a companion mass of 0.3 $M_\odot$ (indicative of a white dwarf or low-mass star), and a semi-major axis of 2000 AU. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios.

Pulsar timing array experiments have recently uncovered evidence for a nanohertz gravitational wave background by precisely timing an ensemble of millisecond pulsars. The next significant milestones for these experiments include characterizing the detected background with greater precision, identifying its source(s), and detecting continuous gravitational waves from individual supermassive black hole binaries. To achieve these objectives, generating accurate and precise times of arrival of pulses from pulsar observations is crucial. Incorrect polarization calibration of the observed pulsar profiles may introduce errors in the measured times of arrival. Further, previous studies (e.g., van Straten 2013; Manchester et al. 2013) have demonstrated that robust polarization calibration of pulsar profiles can reduce noise in the pulsar timing data and improve timing solutions. In this paper, we investigate and compare the impact of different polarization calibration methods on pulsar timing precision using three distinct calibration techniques: the Ideal Feed Assumption (IFA), Measurement Equation Modeling (MEM), and Measurement Equation Template Matching (METM). Three NANOGrav pulsars-PSRs J1643$-$1224, J1744$-$1134, and J1909$-$3744-observed with the 800 MHz and 1.5 GHz receivers at the Green Bank Telescope (GBT) are utilized for our analysis. Our findings reveal that all three calibration methods enhance timing precision compared to scenarios where no polarization calibration is performed. Additionally, among the three calibration methods, the IFA approach generally provides the best results for timing analysis of pulsars observed with the GBT receiver system. We attribute the comparatively poorer performance of the MEM and METM methods to potential instabilities in the reference noise diode coupled to the receiver and temporal variations in the profile of the reference pulsar, respectively., Comment: 21 pages, 6 figures, 1 table, Accepted for Publication in Astrophysical Journal

In 2007 we were part of a team that discovered the so-called ``Lorimer Burst'', the first example of a new class of objects now known as fast radio bursts (FRBs). These enigmatic events are only a few ms in duration and occur at random locations on the sky at a rate of a few thousand per day. Several thousand FRBs are currently known. While it is now well established that they have a cosmological origin, and about 10\% of all currently known sources have been seen to exhibit multiple bursts, the origins of these enigmatic sources are currently poorly understood. In this article, we review the discovery of FRBs and present some of the highlights from the vast body of work by an international community. Following a brief overview of the scale of the visible Universe in \S 1, we describe the key moments in radio astronomy (\S 2) that led up to the discovery of the Lorimer burst (\S 3). Early efforts to find more FRBs are described in \S 4 which led to the discovery of the first repeating source (\S 5). In \S 6, as we close out on the second decade of FRBs, we outline some of the many open questions in the field and look ahead to the coming years where many surprises are surely in store., Comment: 23 pages, 9 figures, accepted for publication in Astrophysics and Space Science; manuscript based on the 2023 Shaw Prize Lecture in Astronomy delivered by the authors

We present the discovery of 8 isolated millisecond pulsars in Globular Cluster (GC) NGC 6517 using the Five-Hundred-meter Aperture Spherical radio Telescope (FAST). The spin periods of those pulsars (namely PSR J1801-0857K to R, or, NGC 6517K to R) are all shorter than 10 ms. With these discoveries, NGC 6517 is currently the GC with the most known pulsars in the FAST sky. The largest difference in dispersion measure of the pulsars in NGC 6517 is 11.2 cm$^{-3}$ pc, the second among all GCs. The fraction of isolated pulsars in this GC (16 of 17, 94$\%$) is consistent with previous studies indicating an overabundance of isolated pulsars in the densest GCs, especially in those undergoing cluster core collapse. Considering the FAST GC pulsar discoveries, we modeled the GC pulsar population using the empirical Bayesian method described by Turk and Lorimer with the recent counts. Using this approach, we find that the expected number of potential pulsars in GCs seems to be correlated with the central escape velocity, hence, the GCs Liller 1, NGC 6441, M54 (NGC 6715), and $\omega$-Cen (NGC 5139) are expected to host the largest numbers of pulsars., Comment: 21 pages, 2 figures, accepted for publication in The Astrophysical Journal Letters

Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect achromatic noise properties in each pulsar. The models we compare include existing DM models used by NANOGrav and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar dataset. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from $\log_{10}A_{\text{RN}} = -14.1^{+0.1}_{-0.1}$ to $-14.7^{+0.3}_{-0.5}$, and the spectral index broadens from $\gamma_{\text{RN}} = 2.6^{+0.5}_{-0.4}$ to $\gamma_{\text{RN}} = 3.5^{+1.2}_{-0.9}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs.

Materials used to study nuclear fusion can retain atmospheric helium unless pretreated before an experiment. Understanding helium outgassing is important for accurate diagnostics in experiments surrounding nuclear fusion. The presence of helium is often cited as the primary evidence that a nuclear reaction has occurred, so it is imperative that known sources of helium are mitigated prior to proceeding with novel nuclear experiments. It is also necessary to ensure hermiticity when transferring gas aliquots from an experiment to a mass spectrometer. In this article, we present studies of detecting helium leak rates in systems used in novel nuclear experiments. We also present studies of helium retention in materials subjected to various heating profiles and atmospheric concentrations. Without pretreatment, stainless-steel 316 retains between 15 $\unicode{x2013}$ 240 pmol of $^{ 4}$He or an areal outgassing amount of 0.07 $\unicode{x2013}$ 1.20 pmol/$cm^{ 2}$. It also may reabsorb $^{ 4}$He from the atmosphere in time. These studies also demonstrate that it is necessary to pretreat most materials prior to performing experiments where the presence of $^{ 4}$He is being used as an indicator for novel nuclear reactions., Comment: 10 pages, 4 figures, 5 tables

Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis., Comment: 11 pages, 5 figures

The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $\eta>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $\eta>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations., Comment: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

Laboratories and Universities regularly apply for approval from the United States Nuclear Regulatory Commission (NRC) to use neutron generators for experimental research. To comply with the regulations set by the NRC, adequate shielding is necessary to ensure that the radiation dose rates experienced by an operator, and outside the walls of any containment buildings, are below the prescribed levels. Typically, the neutron source needs to be shielded such that the radiation dose rate experienced by any user is less than 0.25 mRem/hr (500mrem/2000-hour work-year). To address this requirement, we investigate the effectiveness of boronated concrete, boronated water, and light water shielding materials and their applicability to three neutron sources. We present our findings on the radiation shielding design and calculations for three neutron sources situated inside shielding layers. Our modeling utilized the Monte Carlo n-Particle transport codes (MCNP6.2) to simulate neutron attenuation of the shielding. The simulation results reveal that a light water shielding can sufficiently reduce the dose rate for an individual located as close as two meters from the source. Therefore, this shield design can effectively decrease the radiation dose below the maximum recommended limit., Comment: 9 pages, 6 figures, 2 tables

In recent years, the CHIME (Canadian Hydrogen Intensity Mapping Experiment) interferometer has revealed a large number of Fast Radio Bursts (FRBs), including a sizable population that demonstrates repeating behavior. This transit facility, employing a real-time FRB search pipeline, continually scans the sky with declinations between $-10^{\circ}$ and $90^{\circ}$ for events with fluences $\gtrapprox 0.4$ Jy ms. We simulate a population of repeating FRBs by performing Monte Carlo simulations of underlying source populations processed through a mock CHIME/FRB observing pipeline. Assuming intrinsic repeater rates follow a Poisson distribution, we test assumptions about the burst populations of the repeater sample, and construct models of the FRB sample assuming various cosmological distributions. We infer the completeness of CHIME/FRB observations as a function of observing cadence and redshifts out to 0.5. We find that, if all simulated bursts have a fixed Poisson probability of repetition over their integrated time of observation, repeating burst detections across comoving volume should continue to grow near linearly on the order of decades. We predict that around 170 of the current CHIME/FRB one-off sources will ultimately repeat. We also make projections for FRB repeaters by future facilities and demonstrate that the number of repeaters they find could saturate on a $\sim$3 yr timescale., Comment: 10 pages, 7 figures

Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them., Comment: 14 pages, 8 figures, version matching published paper

X-ray diffraction (XRD) is an essential technique to determine a material's crystal structure in high-throughput experimentation, and has recently been incorporated in artificially intelligent agents in autonomous scientific discovery processes. However, rapid, automated and reliable analysis method of XRD data matching the incoming data rate remains a major challenge. To address these issues, we present CrystalShift, an efficient algorithm for probabilistic XRD phase labeling that employs symmetry-constrained pseudo-refinement optimization, best-first tree search, and Bayesian model comparison to estimate probabilities for phase combinations without requiring phase space information or training. We demonstrate that CrystalShift provides robust probability estimates, outperforming existing methods on synthetic and experimental datasets, and can be readily integrated into high-throughput experimental workflows. In addition to efficient phase-mapping, CrystalShift offers quantitative insights into materials' structural parameters, which facilitate both expert evaluation and AI-based modeling of the phase space, ultimately accelerating materials identification and discovery., Comment: 13 pages, 6 figures

An electron beam may be used to generate bremsstrahlung photons that go on to create photoneutrons within metals. This serves as a low-energy neutron source for irradiation experiments [1-3]. In this article, we present simulation results for optimizing photoneutron yield for a 10-MeV electron beam on tantalum foil and erbium deuteride (ErD$_3$). The thickness of the metal layers was varied. A tantalum foil thickness of 1.5 mm resulted in the most photons reaching the second metal layer. When a second metal layer of ErD$_3$ was included, the photoneutron yield increased with the thickness of the secondary layer. When the electron beam was directly incident upon a layer of ErD$_3$, the photoneutron yield did not differ significantly from the yield when a layer of tantalum was included. The directional photoneutron yield reached a maximum level when the thickness of the ErD$_3$ layer was around 12 cm. About 1 neutron was generated per $10^4$ source electrons. When using a 2-mA beam current, it is possible to generate up to $10^{ 12}$ neutrons per second, making this combination a relatively-inexpensive neutron generator., Comment: 15 pages, 9 figures, 4 tables

In this paper we discuss the advantages and challenges of utilizing Fission Fragment Rocket Engines (FFREs) to dramatically reduce transit time in space travel, for example, traveling to Mars. We discuss methods to decrease the size and weight of FFREs. These include utilizing metallic deuterides as moderators, driving the engines with electron beam bremsstrahlung, and operating the FFREs as subcritical assemblies, not as nuclear reactors. We discuss these and other new innovations based upon improved materials and technology that may be integrated into a revolutionary nuclear rocket technology., Comment: 10 pages, 2 figures, 2 tables

Future advanced nuclear rocket propulsion, and the availability of new nuclear power cycle designs, will benefit substantially from the large current investment in alternative nuclear energy that is underway today. We propose a new nuclear cycle which includes the primary fission of lithium-6, followed by secondary fusion of deuterium and tritium, and a secondary fission of lithium-7 by tritium. This cycle does not produce nuclear waste from its nuclear fuel, since all byproducts of these cascade reactions are stable, provided that the triton production during the primary reaction is fully consumed in the secondary reactions. This cycle may, however, activate surrounding technical materials from its neutron flux. This light-element nuclear fuel is readily obtained through the ongoing expansion of the lithium mining industry and electric vehicle (EV) battery recycling industries., Comment: 10 pages, 2 figures

We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately $3.3 \times 10^{-14}$. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613$-$0200. This suggests that this pulsar has some excess noise which can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array., Comment: 29 pages, 5 figures

Noise characterization for pulsar-timing applications accounts for interstellar dispersion by assuming a known frequency-dependence of the delay it introduces in the times of arrival (TOAs). However, calculations of this delay suffer from mis-estimations due to other chromatic effects in the observations. The precision in modeling dispersion is dependent on the observed bandwidth. In this work, we calculate the offsets in infinite-frequency TOAs due to mis-estimations in the modeling of dispersion when using varying bandwidths at the Green Bank Telescope. We use a set of broadband observations of PSR J1643-1224, a pulsar with an excess of chromatic noise in its timing residuals. We artificially restricted these observations to a narrowband frequency range, then used both data sets to calculate residuals with a timing model that does not include short-scale dispersion variations. By fitting the resulting residuals to a dispersion model, and comparing the ensuing fitted parameters, we quantify the dispersion mis-estimations. Moreover, by calculating the autocovariance function of the parameters we obtained a characteristic timescale over which the dispersion mis-estimations are correlated. For PSR J1643-1224, which has one of the highest dispersion measures (DM) in the NANOGrav pulsar timing array, we find that the infinite-frequency TOAs suffer from a systematic offset of ~22 microseconds due to DM mis-estimations, with correlations over ~1 month. For lower-DM pulsars, the offset is ~7 microseconds. This error quantification can be used to provide more robust noise modeling in NANOGrav's data, thereby increasing sensitivity and improving parameter estimation in gravitational wave searches.

Kirigami-inspired designs can enable self-folding three-dimensional materials from flat, two-dimensional sheets. Hierarchical designs of connected levels increase the diversity of possible target structures, yet they can lead to longer folding times in the presence of fluctuations. Here, we study the effect of rotational coupling between levels on the self-folding of two-level kirigami designs driven by thermal noise in a fluid. Naturally present due to hydrodynamic resistance, we find that optimization of this coupling as control parameter can significantly improve a structure's self-folding rate and yield.

This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data., Comment: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email comments@nanograv.org

Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems., Comment: 23 pages, 13 figures, 2 tables. Accepted for publication in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data., Comment: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org