Search

The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentially explained by a progenitor system that is asymmetric in terms of the stellar component masses, i.e. with a mass ratio $q$ of 0.7-0.8, the known Galactic population of merging double neutron star (DNS) systems (i.e. those that will coalesce within billions of years or less) has, until now, only consisted of nearly equal-mass ($q > 0.9$) binaries. PSR J1913+1102 is a DNS system in a 5-hour, low-eccentricity ($e = 0.09$) orbit, implying an orbital separation of 1.8 solar radii, with the two neutron stars predicted to coalesce in 470 million years due to gravitational-wave emission. Here we report that the masses of the two neutron stars, as measured by a dedicated pulsar timing campaign, are $1.62 \pm 0.03$ and $1.27 \pm 0.03$ solar masses for the pulsar and companion neutron star, respectively; with a measured mass ratio $q = 0.78 \pm 0.03$, it is the most asymmetric DNS among known merging systems. Based on this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30% (90% confidence) of the total population of merging DNS binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event., 17 pages, 3 figures, published in Nature on 9 July 2020 (https://www.nature.com/articles/s41586-020-2439-x); this is the authors' final submitted version before review. Includes methods and supplementary information

Statistical anisotropy in the nanohertz-frequency gravitational-wave background (GWB) is expected to be detected by pulsar timing arrays (PTAs) in the near future. By developing a frequentist statistical framework that intrinsically restricts the GWB power to be positive, we establish scaling relations for multipole-dependent anisotropy decision thresholds that are a function of the noise properties, timing baselines, and cadences of the pulsars in a PTA. We verify that $(i)$ a larger number of pulsars, and $(ii)$ factors that lead to lower uncertainty on the cross-correlation measurements between pulsars, lead to a higher overall GWB signal-to-noise ratio, and lower anisotropy decision thresholds with which to reject the null hypothesis of isotropy. Using conservative simulations of realistic NANOGrav datasets, we predict that an anisotropic GWB with angular power $C_{l=1} > 0.3\,C_{l=0}$ may be sufficient to produce tension with isotropy at the $p = 3\times10^{-3}$ ($\sim3\sigma$) level in near-future NANOGrav data with a $20$~yr baseline. We present ready-to-use scaling relationships that can map these thresholds to any number of pulsars, configuration of pulsar noise properties, and sky coverage. We discuss how PTAs can improve the detection prospects for anisotropy, as well as how our methods can be adapted for more versatile searches., Comment: Submitted to ApJ. Comments welcome

The characterization of nanohertz-frequency gravitational waves (GWs) with pulsar-timing arrays requires a continual expansion of datasets and monitored pulsars. Whereas detection of the stochastic GW background is predicated on measuring a distinctive pattern of inter-pulsar correlations, characterizing the background's spectrum is driven by information encoded in the power spectra of the individual pulsars' time series. We propose a new technique for rapid Bayesian characterization of the stochastic GW background that is fully parallelized over pulsar datasets. This Factorized Likelihood (FL) technique empowers a modular approach to parameter estimation of the GW background, multi-stage model selection of a spectrally-common stochastic process and quadrupolar inter-pulsar correlations, and statistical cross-validation of measured signals between independent pulsar sub-arrays. We demonstrate the equivalence of this technique's efficacy with the full pulsar-timing array likelihood, yet at a fraction of the required time. Our technique is fast, easily implemented, and trivially allows for new data and pulsars to be combined with legacy datasets without re-analysis of the latter., Comment: 14 pages, 6 figures. Matches version accepted by PRD

We report results from continued timing observations of PSR J0740+6620, a high-mass, 2.8-ms radio pulsar in orbit with a likely ultra-cool white dwarf companion. Our data set consists of combined pulse arrival-time measurements made with the 100-m Green Bank Telescope and the Canadian Hydrogen Intensity Mapping Experiment telescope. We explore the significance of timing-based phenomena arising from general-relativistic dynamics and variations in pulse dispersion. When using various statistical methods, we find that combining $\sim 1.5$ years of additional, high-cadence timing data with previous measurements confirms and improves upon previous estimates of relativistic effects within the PSR J0740+6620 system, with the pulsar mass $m_{\rm p} = 2.08^{+0.07}_{-0.07}$ M$_\odot$ (68.3\% credibility) determined by the relativistic Shapiro time delay. For the first time, we measure secular variation in the orbital period and argue that this effect arises from apparent acceleration due to significant transverse motion. After incorporating contributions from Galactic differential rotation and off-plane acceleration in the Galactic potential, we obtain a model-dependent distance of $d = 1.14^{+0.17}_{-0.15}$ kpc (68.3\% credibility). This improved distance confirms the ultra-cool nature of the white dwarf companion determined from recent optical observations. We discuss the prospects for future observations with next-generation facilities, which will likely improve the precision on $m_{\rm p}$ for J0740+6620 by an order of magnitude within the next few years., Final version after minor corrections during referee process. Published in the Astrophysical Journal Letters on 1 July 2021

The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its 12.5-yr pulsar timing array dataset, with median characteristic strain amplitude at periods of a year of $A_{\rm yr} = 1.92^{+0.75}_{-0.55} \times 10^{-15}$. However, evidence for the quadrupolar Hellings \& Downs interpulsar correlations, which are characteristic of gravitational wave signals, was not yet significant. We emulate and extend the NANOGrav dataset, injecting a wide range of stochastic gravitational wave background (GWB) signals that encompass a variety of amplitudes and spectral shapes, and quantify three key milestones: (I) Given the amplitude measured in the 12.5 yr analysis and assuming this signal is a GWB, we expect to accumulate robust evidence of an interpulsar-correlated GWB signal with 15--17 yrs of data, i.e., an additional 2--5 yrs from the 12.5 yr dataset; (II) At the initial detection, we expect a fractional uncertainty of $40\%$ on the power-law strain spectrum slope, which is sufficient to distinguish a GWB of supermassive black-hole binary origin from some models predicting more exotic origins;(III) Similarly, the measured GWB amplitude will have an uncertainty of $44\%$ upon initial detection, allowing us to arbitrate between some population models of supermassive black-hole binaries. In addition, power-law models are distinguishable from those having low-frequency spectral turnovers once 20~yrs of data are reached. Even though our study is based on the NANOGrav data, we also derive relations that allow for a generalization to other pulsar-timing array datasets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier., 15 pages, 7 figures

Using neural networks, we integrate the ability to account for Doppler smearing due to a pulsar's orbital motion with the pulsar population synthesis package \psrpoppy\ to develop accurate modeling of the observed binary pulsar population. As a first application, we show that binary neutron star systems where the two components have highly unequal mass are, on average, easier to detect than systems which are symmetric in mass. We then investigate the population of ultra-compact ($1.5 \, {\rm min} \leq P_{\rm b} \leq 15\,\rm min$) neutron star--white dwarf (NS--WD) and double neutron star (DNS) systems which are promising sources for the Laser Interferometer Space Antenna gravitational-wave detector. Given the non-detection of these systems in radio surveys thus far, we estimate a 95\% confidence upper limit of $\sim$1450 and $\sim$1100 ultra-compact NS--WD and DNS systems in the Milky Way that are beaming towards the Earth respectively. We also show that using survey integration times in the range 20~s to 200~s with time-domain resampling will maximize the signal-to-noise ratio as well as the probability of detection of these ultra-compact binary systems. Among all the large scale radio pulsar surveys, those that are currently being carried out at the Arecibo radio telescope have $\sim$50--80\% chance of detecting at least one of these systems using current integration integration times and $\sim$80--95\% using optimal integration times in the next several years., 10 pages, 7 figures, submitted to ApJ. An older version of this analysis appears as a dissertation chapter in arXiv:2008.03842

When galaxies merge, the supermassive black holes in their centers may form binaries and emit low-frequency gravitational radiation in the process. In this paper, we consider the galaxy 3C 66B, which was used as the target of the first multimessenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational-wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C 66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C 66B to less than (1.65 ± 0.02) × 109 M ⊙ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data over “blind” pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

In this dissertation, we estimate the population of different classes of BNS systems that are visible to gravitational-wave observatories. Given that no ultra-compact BNS systems have been discovered in pulsar radio surveys, we place a 95\% confidence upper limit of $\sim$850 and $\sim$1100 ultra-compact neutron star--white dwarf and double neutron star (DNS) systems that are beaming towards the Earth, respectively. We show that among all of the current radio pulsar surveys, the ones at the Arecibo radio telescope have the best chance of detecting an ultra-compact BNS system. We also show that adopting a survey integration time of $t_{\rm int} \sim 1$~min will maximize the signal-to-noise ratio, and thus, the probability of detecting an ultra-compact BNS system. Similarly, we use the sample of nine observed DNS systems to derive a Galactic DNS merger rate of $\mathcal{R}_{\rm MW} = 37^{+24}_{-11}$~Myr$^{-1}$, where the errors represent 90\% confidence intervals. Extrapolating this rate to the observable volume for LIGO, we derive a merger detection rate of $\mathcal{R} = 1.9^{+1.2}_{-0.6} \times \left(D_{\rm r}/100 \ \rm Mpc \right)^3 \rm yr^{-1}$, where $D_{\rm r}$ is the range distance for LIGO. This rate is consistent with that derived using the DNS mergers observed by LIGO. Finally, we measure the sense of rotation of the older millisecond pulsar, pulsar A, in the DNS J0737--3039 system and find that it rotates prograde with respect to its orbit. This is the first direct measurement of the sense of rotation of a pulsar and a direct confirmation of the rotating lighthouse model for pulsars. This result confirms that the spin angular momentum vector is closely aligned with the orbital angular momentum, suggesting that kick of the supernova producing the second born pulsar J0737--3039B was small.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects., 14 pages, 13 figures, fixed typo in abstract of earlier version

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $, Comment: 10 pages, 6 figures, submitted to ApJ

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al. 2021a; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision, multi-band measurements of several pulsars indicate frequency-dependent DMs. Compared to the narrowband analysis in NG12.5, the TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses., Comment: 62 pages, 55 figures, 5 tables, 3 appendices. Data available at http://nanograv.org/data/ and via DOI 10.5281/zenodo.4312887

We present time-of-arrival (TOA) measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of TOA measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for 14 pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications., Comment: 54 pages, 52 figures, 7 tables, 1 appendix. Data are available at http://nanograv.org/data/ and via the DOI 10.5281/zenodo.4312297

We calculate the dispersion measure (DM) contributed by the intergalactic medium (IGM) to the total measured DM for fast radio bursts (FRBs). We use the MareNostrum Instituto de Ciencias del Espacio (MICE) Onion Universe simulation (Fosalba et al. 2008) to track the evolution of the dark matter particle density over a large range of redshifts. We convert this dark matter particle number density to the corresponding free electron density and then integrate it to find the DM as a function of redshift. This approach yields an intergalactic DM of $\rm DM_{\rm IGM}(z = 1) = 800^{+7000}_{-170}$ pc cm$^{-3}$, with the large errors representative of the structure in the IGM. We place limits on the redshifts of the current population of observed FRBs. We also use our results to estimate the host galaxy contribution to the DM for the first repeater, FRB 121102, and show that the most probable host galaxy DM contribution, $\rm DM_{\rm host} \approx 310$ pc cm$^{-3}$, is consistent with the estimate made using the Balmer emission lines in the spectrum of the host galaxy, $\rm DM_{\rm Balmer} = 324$ pc cm$^{-3}$ (Tendulkar et al. 2017). We also compare our predictions for the host galaxy contribution to the DM for the observations of FRB 180924 (Bannister et al. 2019) and FRB 190523 (Ravi et al. 2019), both of which have been localized to a host galaxy., Resubmitted to ApJ following second round of referee comments. Expanded sub-section comparing our results to those from other works (Sec. 3.3)

We present the discovery and timing of the young (age $\sim 28.6$ kyr) pulsar PSR J0837$-$2454. Based on its high latitude ($b = 9.8^{\circ}$) and dispersion measure (DM $ = 143$~pc~cm$^{-3}$), the pulsar appears to be at a $z$-height of $>$1 kpc above the Galactic plane, but near the edge of our Galaxy. This is many times the observed scale height of the canonical pulsar population, which suggests this pulsar may have been born far out of the plane. If accurate, the young age and high $z$-height imply that this is the first pulsar known to be born from a runaway O/B star. In follow-up imaging with the Australia Telescope Compact Array (ATCA), we detect the pulsar with a flux density $S_{1400} = 0.18 \pm 0.05$ mJy. We do not detect an obvious supernova remnant around the pulsar in our ATCA data, but we detect a co-located, low-surface-brightness region of $\sim$1.5$^\circ$ extent in archival Galactic and Extragalactic All-sky MWA Survey data. We also detect co-located H$\alpha$ emission from the Southern H$\alpha$ Sky Survey Atlas. Distance estimates based on these two detections come out to $\sim$0.9 kpc and $\sim$0.2 kpc respectively, both of which are much smaller than the distance predicted by the NE2001 model ($6.3$ kpc) and YMW model ($>25$ kpc) and place the pulsar much closer to the plane of the Galaxy. If the pulsar/remnant association holds, this result also highlights the inherent difficulty in the classification of transients as "Galactic" (pulsar) or "extragalactic" (fast radio burst) toward the Galactic anti-center based solely on the modeled Galactic electron contribution to a detection., Comment: Published in ApJ. 12 pages, 9 figures, 2 tables

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature., Comment: 25 pages, 14 figures, 5 tables, 3 appendices. Published in The Astrophysical Journal Letters. Please send any comments/questions to Joseph Simon (joe.simon@nanograv.org). Jupyter notebook tutorials and some MCMC chain files are available at https://github.com/nanograv/12p5yr_stochastic_analysis

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}, 21 pages, 12 figures, 9 tables. Published in The Astrophysical Journal. Please send any comments/questions to S. R. Taylor (srtaylor@caltech.edu)

We report the discovery and initial follow-up of a double neutron star (DNS) system, PSR J1946$+$2052, with the Arecibo L-Band Feed Array pulsar (PALFA) survey. PSR J1946$+$2052 is a 17-ms pulsar in a 1.88-hour, eccentric ($e \, =\, 0.06$) orbit with a $\gtrsim 1.2 \, M_\odot$ companion. We have used the Jansky Very Large Array to localize PSR J1946$+$2052 to a precision of 0.09 arcseconds using a new phase binning mode. We have searched multiwavelength catalogs for coincident sources but did not find any counterparts. The improved position enabled a measurement of the spin period derivative of the pulsar ($\dot{P} \, = \, 9\,\pm \, 2 \,\times 10^{-19}$); the small inferred magnetic field strength at the surface ($B_S \, = \, 4 \, \times \, 10^9 \, \rm G$) indicates that this pulsar has been recycled. This and the orbital eccentricity lead to the conclusion that PSR J1946$+$2052 is in a DNS system. Among all known radio pulsars in DNS systems, PSR J1946$+$2052 has the shortest orbital period and the shortest estimated merger timescale, 46 Myr; at that time it will display the largest spin effects on gravitational wave waveforms of any such system discovered to date. We have measured the advance of periastron passage for this system, $\dot{\omega} \, = \, 25.6 \, \pm \, 0.3\, \deg \rm yr^{-1}$, implying a total system mass of only 2.50 $\pm$ 0.04 $M_\odot$, so it is among the lowest mass DNS systems. This total mass measurement combined with the minimum companion mass constrains the pulsar mass to $\lesssim 1.3 \, M_\odot$., Comment: Accepted for publication by ApJL, 8 pages, 3 figures

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs., 10 pages, 11 figures. Accepted by Astrophysical Journal. Please send any comments/questions to S. J. Vigeland (vigeland@uwm.edu)

Following the procedure described in Pol et al., we update the Galactic double neutron star (DNS) merger rate by including the new, highly eccentric DNS system J0509+3801 (Lynch et al.). This leads to a new Galactic DNS merger rate of $\mathcal{R}_{\rm MW} = 37^{+24}_{-11}$ Myr$^{-1}$, where the errors represent 90% confidence intervals. The corresponding DNS merger detection rate for Advanced LIGO is $\mathcal{R} = 1.9^{+1.2}_{-0.6} \times \left( D_{\rm r} / 100 \ \rm Mpc \right)^3 \rm yr^{-1},$ where $D_{\rm r}$ is the range distance. Using the LIGO O3 range distance of 130 Mpc (Abbott et al.), we predict that LIGO will detect anywhere between three and seven DNS mergers per year of observing at O3 sensitivity., Published in Research Notes of the AAS. The data for figure 1 is available with the RNAAS publication

We apply the algorithm published by Liang et al. (2014) to describe the Double Pulsar system J0737$-$3039 and extract the sense of rotation of first born recycled pulsar PSR J0737$-$3039A. We find that this pulsar is rotating prograde in its orbit. This is the first direct measurement of the sense of rotation of a pulsar with respect to its orbit and a direct confirmation of the rotating lighthouse model for pulsars. This result confirms that the spin angular momentum vector is closely aligned with the orbital angular momentum, suggesting that kick of the supernova producing the second born pulsar J0737$-$3039B was small., Comment: Accepted for publication in ApJ

We report deep, single-dish 21 cm observations of NGC 1052-DF2, taken with the Green Bank Telescope. NGC 1052-DF2, proposed to be lacking in dark matter, is currently classified as an ultra-diffuse galaxy in the NGC 1052 group. We do not detect the galaxy, and derive an upper limit on the HI, mass. The galaxy is extremely gas-poor, and we find that a $3��\, M_{HI}$ detection at a distance of 19 Mpc and using a line width of 3.2 km $\rm s^{-1}$ would have an upper limit of $M_{HI,lim} < 5.5 \times 10^5$ M$_{\odot}$. At this mass limit, the gas fraction of neutral gas mass to stellar mass is extremely low, at $M_{HI}$/M$_{*}$ $\, < \, 0.0027$. This extremely low gas fraction, comparable to Galactic dwarf spheroidals and gas-poor dwarf ellipticals, implies that either the galaxy is within the virial radius of NGC1052, where its gas has been stripped due to its proximity to the central galaxy, or that NGC 1052-DF2 is at distance large enough to inhibit detection of its gas. We also estimated the upper limit of the HI mass of NGC 1052-DF2 resided at 13 Mpc. This would give an HI mass of $M_{HI,lim} < 2.5 \times 10^5$ M$_{\odot}$, and an HI gas fraction of $M_{HI}$/M$_{*}$r $\, < \, 0.0012$, becoming even more extreme. While the dark matter fraction would be less extreme at this distance, the neutral gas fraction would be unprecedented for an object in a low density environment., Accepted for publication in ApJL

We present the Galactic merger rate for double neutron star (DNS) binaries using the observed sample of eight DNS systems merging within a Hubble time. This sample includes the recently discovered, highly relativistic DNS systems J1757$-$1854 and J1946+2052, and is approximately three times the sample size used in previous estimates of the Galactic merger rate by Kim et al. Using this sample, we calculate the vertical scale height for DNS systems in the Galaxy to be $z_0 = 0.4 \pm 0.1$ kpc. We calculate a Galactic DNS merger rate of $\mathcal{R}_{\rm MW} = 42^{+30}_{-14}$ Myr$^{-1}$ at the 90% confidence level. The corresponding DNS merger detection rate for Advanced LIGO is $\mathcal{R}_{\rm LIGO} = 0.18^{+0.13}_{-0.06} \times \left( D_{\rm r}/100 \ \rm Mpc \right)^3 \rm yr^{-1}$, where $D_{\rm r}$ is the range distance. Using this merger detection rate and the predicted range distance of 120$-$170 Mpc for the third observing run of LIGO (Laser Interferometer Gravitational-wave Observatory, Abbott et al.), we predict, accounting for 90% confidence intervals, that LIGO$-$Virgo will detect anywhere between zero and two DNS mergers. We explore the effects of the underlying pulsar population properties on the merger rate and compare our merger detection rate with those estimated using different formation and evolutionary scenario of DNS systems. As we demonstrate, reconciling the rates is sensitive to assumptions about the DNS population, including its radio pulsar luminosity function. Future constraints from further gravitational wave DNS detections and pulsar surveys anticipated in the near future should permit tighter constraints on these assumptions., 12 pages, 6 figures, published in ApJ; Update acknowledgements for authors in this version

We present a rest-frame composite spectrum for Seyfert 1 galaxies using spectra obtained from the DR12 release of the Sloan Digital Sky Survey (SDSS). The spectrum is constructed by combining data from a total of 10,112 galaxies, spanning a redshift range of 0 to 0.793. We produce an electronic table of the median and geometric mean composite Seyfert 1 spectrum. We measure the spectral index of the composite spectrum, and compare it with that of the composite quasar spectrum. We also measure the flux and width of the strong emission lines present in the composite spectrum. We compare the entire spectrum with the quasar spectrum in the context of the AGN unification model. The two composite spectra match extremely well in the blue part of the spectrum, while there is an offset in flux in the red portion of the spectrum., Comment: 11 pages, 15 figures, 4 tables, accepted for publication in MNRAS. The median, geometric mean and ratio spectra are available online at: http://www.ncra.tifr.res.in:8081/~yogesh/seyfert1compositespectrum/ . The median spectra binned in redshift are also made available on the same page