Showing 1–46 of 46 results for author: Faber, J

A Heavily Scattered Fast Radio Burst Is Viewed Through Multiple Galaxy Halos

Authors: Jakob T. Faber, Vikram Ravi, Stella Koch Ocker, Myles B. Sherman, Kritti Sharma, Liam Connor, Casey Law, Nikita Kosogorov, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, James W. Lamb, Paul Rasmussen, Jean J. Somalwar, Sander Weinreb, David P. Woody

Abstract: We present a multi-wavelength study of the apparently non-repeating, heavily scattered fast radio burst, FRB 20221219A, detected by the Deep Synoptic Array 110 (DSA-110). The burst exhibits a moderate dispersion measure (DM) of $706.7^{+0.6}_{-0.6}$ $\mathrm{pc}~\mathrm{cm}^{-3}$ and an unusually high scattering timescale of $τ_{\mathrm{obs}} = 19.2_{-2.7}^{+2.7}$ ms at 1.4 GHz. We associate the F… ▽ More We present a multi-wavelength study of the apparently non-repeating, heavily scattered fast radio burst, FRB 20221219A, detected by the Deep Synoptic Array 110 (DSA-110). The burst exhibits a moderate dispersion measure (DM) of $706.7^{+0.6}_{-0.6}$ $\mathrm{pc}~\mathrm{cm}^{-3}$ and an unusually high scattering timescale of $τ_{\mathrm{obs}} = 19.2_{-2.7}^{+2.7}$ ms at 1.4 GHz. We associate the FRB with a Milky Way-like host galaxy at $z_{\mathrm{host}} = 0.554$ of stellar mass $\mathrm{log}_{10}(M_{\star, \mathrm{host}}) = 10.20^{+0.04}_{-0.03} ~M_\odot$. We identify two intervening galaxy halos at redshifts $z_{\mathrm{igh1}} = 0.492$ and $z_{\mathrm{igh2}} = 0.438$, with low impact parameters, $b_{\mathrm{igh1}} = 43.0_{-11.3}^{+11.3}$ kpc and $b_{\mathrm{igh2}} = 36.1_{-11.3}^{+11.3}$ kpc, and intermediate stellar masses, $\mathrm{log}_{10}(M_{\star, \mathrm{igh1}}) = 10.01^{+0.02}_{-0.02} ~M_\odot$ and $\mathrm{log}_{10}(M_{\star, \mathrm{igh2}}) = 10.60^{+0.02}_{-0.02} ~M_\odot$. The presence of two such galaxies suggests that the sightline is significantly overcrowded compared to the median sightline to this redshift, as inferred from the halo mass function. We perform a detailed analysis of the sightline toward FRB 20221219A, constructing both DM and scattering budgets. Our results suggest that, unlike most well-localized sources, the host galaxy does not dominate the observed scattering. Instead, we posit that an intersection with a single partially ionized cloudlet in the circumgalactic medium of an intervening galaxy could account for the substantial scattering in FRB 20221219A and remain in agreement with typical electron densities inferred for extra-planar dense cloud-like structures in the Galactic and extragalactic halos (e.g., high-velocity clouds). △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 18 pages, 6 figures, submitted to ApJ, comments appreciated

Searching for Magnetar Binaries Disrupted by Core-Collapse Supernovae

Authors: Myles B. Sherman, Vikram Ravi, Kareem El-Badry, Kritti Sharma, Stella Koch Ocker, Nikita Kosogorov, Liam Connor, Jakob T. Faber

Abstract: Core-collapse Supernovae (CCSNe) are considered the primary magnetar formation channel, with 15 magnetars associated with supernova remnants (SNRs). A large fraction of these should occur in massive stellar binaries that are disrupted by the explosion, meaning that $\sim45\%$ of magnetars should be nearby high-velocity stars. Here we conduct a multi-wavelength search for unbound stars, magnetar bi… ▽ More Core-collapse Supernovae (CCSNe) are considered the primary magnetar formation channel, with 15 magnetars associated with supernova remnants (SNRs). A large fraction of these should occur in massive stellar binaries that are disrupted by the explosion, meaning that $\sim45\%$ of magnetars should be nearby high-velocity stars. Here we conduct a multi-wavelength search for unbound stars, magnetar binaries, and SNR shells using public optical ($uvgrizy-$bands), infrared ($J-$, $H-$, $K-$, and $K_s-$bands), and radio ($888$ MHz, $1.4$ GHz, and $3$ GHz) catalogs. We use Monte Carlo analyses of candidates to estimate the probability of association with a given magnetar based on their proximity, distance, proper motion, and magnitude. In addition to recovering a proposed magnetar binary, a proposed unbound binary, and 13 of 15 magnetar SNRs, we identify two new candidate unbound systems: an OB star from the Gaia catalog we associate with SGR J1822.3-1606, and an X-ray pulsar we associate with 3XMM J185246.6+003317. Using a Markov-Chain Monte Carlo simulation that assumes all magnetars descend from CCSNe, we constrain the fraction of magnetars with unbound companions to $5\lesssim f_u \lesssim 24\%$, which disagrees with neutron star population synthesis results. Alternate formation channels are unlikely to wholly account for the lack of unbound binaries as this would require $31\lesssim f_{nc} \lesssim 66\%$ of magnetars to descend from such channels. Our results support a high fraction ($48\lesssim f_m \lesssim 86\%$) of pre-CCSN mergers, which can amplify fossil magnetic fields to preferentially form magnetars. △ Less

Submitted 15 May, 2024; v1 submitted 7 April, 2024; originally announced April 2024.

Comments: 35 pages, 8 figures, 5 tables, accepted to MNRAS after revision; updated XTE J1810 data, pulsar test analysis, MW neutron star census

Morphologies of Bright Complex Fast Radio Bursts with CHIME/FRB Voltage Data

Authors: Jakob T. Faber, Daniele Michilli, Ryan Mckinven, Jianing Su, Aaron B. Pearlman, Kenzie Nimmo, Robert A. Main, Victoria Kaspi, Mohit Bhardwaj, Shami Chatterjee, Alice P. Curtin, Matt Dobbs, Gwendolyn Eadie, B. M. Gaensler, Zarif Kader, Calvin Leung, Kiyoshi W. Masui, Ayush Pandhi, Emily Petroff, Ziggy Pleunis, Masoud Rafiei-Ravandi, Ketan R. Sand, Paul Scholz, Kaitlyn Shin, Kendrick Smith , et al. (1 additional authors not shown)

Abstract: We present the discovery of twelve thus far non-repeating fast radio burst (FRB) sources, detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources were selected from a database comprising of order $10^3$ CHIME/FRB full-array raw voltage data recordings, based on their exceptionally high brightness and complex morphology. Our study examines the time-frequency… ▽ More We present the discovery of twelve thus far non-repeating fast radio burst (FRB) sources, detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources were selected from a database comprising of order $10^3$ CHIME/FRB full-array raw voltage data recordings, based on their exceptionally high brightness and complex morphology. Our study examines the time-frequency characteristics of these bursts, including drifting, microstructure, and periodicities. The events in this sample display a variety of unique drifting phenomenologies that deviate from the linear negative drifting phenomenon seen in many repeating FRBs, and motivate a possible new framework for classifying drifting archetypes. Additionally, we detect microstructure features of duration $\lesssim$ 50 $μs$ in seven events, with some as narrow as $\approx$ 7 $μs$. We find no evidence of significant periodicities. Furthermore, we report the polarization characteristics of seven events, including their polarization fractions and Faraday rotation measures (RMs). The observed $|\mathrm{RM}|$ values span a wide range of $17.24(2)$ - $328.06(2) \mathrm{~rad~m}^{-2}$, with linear polarization fractions between $0.340(1)$ - $0.946(3)$. The morphological properties of the bursts in our sample appear broadly consistent with predictions from both relativistic shock and magnetospheric models of FRB emission, as well as propagation through discrete ionized plasma structures. We address these models and discuss how they can be tested using our improved understanding of morphological archetypes. △ Less

Submitted 26 December, 2023; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 29 pages, 14 figures, submitted to ApJ, comments appreciated

Deep Synoptic Array Science: Implications of Faraday Rotation Measures of Localized Fast Radio Bursts

Authors: Myles B. Sherman, Liam Connor, Vikram Ravi, Casey Law, Ge Chen, Kritti Sharma, Morgan Catha, Jakob T. Faber, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, James W. Lamb, Paul Rasmussen, Jun Shi, Dana Simard, Jean Somalwar, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: Faraday rotation measures (RMs) of fast radio bursts (FRBs) offer the prospect of directly measuring extragalactic magnetic fields. We present an analysis of the RMs of ten as yet non-repeating FRBs detected and localized to host galaxies by the 110-antenna Deep Synoptic Array (DSA-110). We combine this sample with published RMs of 15 localized FRBs, nine of which are repeating sources. For each F… ▽ More Faraday rotation measures (RMs) of fast radio bursts (FRBs) offer the prospect of directly measuring extragalactic magnetic fields. We present an analysis of the RMs of ten as yet non-repeating FRBs detected and localized to host galaxies by the 110-antenna Deep Synoptic Array (DSA-110). We combine this sample with published RMs of 15 localized FRBs, nine of which are repeating sources. For each FRB in the combined sample, we estimate the host-galaxy dispersion measure (DM) contributions and extragalactic RM. We find compelling evidence that the extragalactic components of FRB RMs are often dominated by contributions from the host-galaxy interstellar medium (ISM). Specifically, we find that both repeating and as yet non-repeating FRBs show a correlation between the host-DM and host-RM in the rest frame, and we find an anti-correlation between extragalactic RM (in the observer frame) and redshift for non-repeaters, as expected if the magnetized plasma is in the host galaxy. Important exceptions to the ISM origin include a dense, magnetized circum-burst medium in some repeating FRBs, and the intra-cluster medium (ICM) of host or intervening galaxy clusters. We find that the estimated ISM magnetic-field strengths, $\bar{B}_{||}$, are characteristically larger than those inferred from Galactic radio pulsars. This suggests either increased ISM magnetization in FRB hosts in comparison with the Milky Way, or that FRBs preferentially reside in regions of increased magnetic-field strength within their hosts. △ Less

Submitted 7 September, 2023; v1 submitted 13 August, 2023; originally announced August 2023.

Comments: 13 pages, 4 figures, 1 table ; added references to section 1, 2, and 3; added a paragraph in Section 3 to discuss simulations addressing interplay of local and ISM magnetic field; added bar to B|| to indicate average over line of sight

Deep Synoptic Array Science: Polarimetry of 25 New Fast Radio Bursts Provides Insights into their Origins

Authors: Myles B. Sherman, Liam Connor, Vikram Ravi, Casey Law, Ge Chen, Morgan Catha, Jakob T. Faber, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, James W. Lamb, Paul Rasmussen, Kritti Sharma, Jun Shi, Dana Simard, Jean Somalwar, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: We report on a full-polarization analysis of the first 25 as yet non-repeating FRBs detected at 1.4 GHz by the 110-antenna Deep Synoptic Array (DSA-110) during commissioning observations. We present details of the data-reduction, calibration, and analysis procedures developed for this novel instrument. Faraday rotation measures (RMs) are searched between $\pm10^6$ rad m$^{-2}$ and detected for 20… ▽ More We report on a full-polarization analysis of the first 25 as yet non-repeating FRBs detected at 1.4 GHz by the 110-antenna Deep Synoptic Array (DSA-110) during commissioning observations. We present details of the data-reduction, calibration, and analysis procedures developed for this novel instrument. Faraday rotation measures (RMs) are searched between $\pm10^6$ rad m$^{-2}$ and detected for 20 FRBs with magnitudes ranging from $4-4670$ rad m$^{-2}$. $15/25$ FRBs are consistent with 100% polarization, 10 of which have high ($\ge70\%$) linear-polarization fractions and 2 of which have high ($\ge30\%$) circular-polarization fractions. Our results disfavor multipath RM scattering as a dominant depolarization mechanism. Polarization-state and possible RM variations are observed in the four FRBs with multiple sub-components. We combine the DSA-110 sample with polarimetry of previously published FRBs, and compare the polarization properties of FRB sub-populations and FRBs with Galactic pulsars. Although FRB polarization fractions are typically higher than those of Galactic pulsars, and cover a wider range than those of pulsar single pulses, they resemble those of the youngest (characteristic ages $<10^{5}$ yr) pulsars. Our results support a scenario wherein FRB emission is intrinsically highly linearly polarized, and propagation effects can result in conversion to circular polarization and depolarization. Young pulsar emission and magnetospheric-propagation geometries may form a useful analogy for the origin of FRB polarization. △ Less

Submitted 15 February, 2024; v1 submitted 13 August, 2023; originally announced August 2023.

Comments: 57 pages, 24 figures, 5 tables ; version accepted by ApJ ; revised classification to include polarized signal-to-noise; added appendices for PA error/variability, J1935 polarization ; corrected literature FRB sample and supplemented with omitted data; added discussion of biases, limitations, and selection criteria; revised sigmaRM analysis ; additional plots added ; revised FRB20220424E analysis

Deep Synoptic Array Science: First FRB and Host Galaxy Catalog

Authors: C. J. Law, K. Sharma, V. Ravi, G. Chen, M. Catha, L. Connor, J. T. Faber, G. Hallinan, C. Harnach, G. Hellbourg, R. Hobbs, D. Hodge, M. Hodges, J. W. Lamb, P. Rasmussen, M. B. Sherman, J. Shi, D. Simard, R. Squillace, S. Weinreb, D. P. Woody, N. Yadlapalli

Abstract: Fast Radio Bursts (FRBs) are a powerful and mysterious new class of transient that are luminous enough to be detected at cosmological distances. By associating FRBs to host galaxies, we can measure intrinsic and environmental properties that test FRB origin models, in addition to using them as precise probes of distant cosmic gas. The Deep Synoptic Array (DSA-110) is a radio interferometer built t… ▽ More Fast Radio Bursts (FRBs) are a powerful and mysterious new class of transient that are luminous enough to be detected at cosmological distances. By associating FRBs to host galaxies, we can measure intrinsic and environmental properties that test FRB origin models, in addition to using them as precise probes of distant cosmic gas. The Deep Synoptic Array (DSA-110) is a radio interferometer built to maximize the rate at which it can simultaneously detect and localize FRBs. Here, we present the first sample of FRBs and host galaxies discovered by the DSA-110. This sample of 11 FRBs is the largest, most uniform sample of localized FRBs to date, as it is selected based on association to host galaxies identified in optical imaging by Pan-STARRS1. These FRBs have not been observed to repeat and their radio properties (dispersion, temporal scattering, energy) are similar to that of the known non-repeating FRB population. Most host galaxies have ongoing star formation, as has been identified before for FRB hosts. Two hosts of the new sample are massive, quiescent galaxies. The distribution of star-formation history across this host-galaxy sample shows that the delay-time distribution is wide, with a powerlaw model that spans from $\sim100$\,Myr to $\gtrsim2$\,Gyr. This requires the existence of one or more progenitor formation channels associated with old stellar populations, such as the binary evolution of compact objects. △ Less

Submitted 24 January, 2024; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 21 pages. Submitted to AAS Journals. Includes changes based on referee comments and improved host galaxy analysis

Detection of ultra-fast radio bursts from FRB 20121102A

Authors: M. P. Snelders, K. Nimmo, J. W. T. Hessels, Z. Bensellam, L. P. Zwaan, P. Chawla, O. S. Ould-Boukattine, F. Kirsten, J. T. Faber, V. Gajjar

Abstract: Fast radio bursts (FRBs) are extragalactic transients with typical durations of milliseconds. FRBs have been shown, however, to fluctuate on a wide range of timescales: some show sub-microsecond sub-bursts while others last up to a few seconds in total. Probing FRBs on a range of timescales is crucial for understanding their emission physics, how to detect them effectively, and how to maximize the… ▽ More Fast radio bursts (FRBs) are extragalactic transients with typical durations of milliseconds. FRBs have been shown, however, to fluctuate on a wide range of timescales: some show sub-microsecond sub-bursts while others last up to a few seconds in total. Probing FRBs on a range of timescales is crucial for understanding their emission physics, how to detect them effectively, and how to maximize their utility as astrophysical probes. FRB 20121102A is the first-known repeating FRB source. Here we show that FRB 20121102A is able to produce isolated microsecond-duration bursts whose total durations are more than ten times shorter than all other known FRBs to date. The polarimetric properties of these micro-bursts resemble those of the longer-lasting bursts, suggesting a common emission mechanism producing FRBs spanning a factor of 1,000 in duration. Furthermore, this work shows that there exists a population of ultra-fast radio bursts that current wide-field FRB searches are missing due to insufficient time-resolution. These results indicate that FRBs occur more frequently and with greater diversity than initially thought. This could also influence our understanding of energy, wait time, and burst rate distributions. △ Less

Submitted 7 November, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: Published in Nature Astronomy, doi.org/10.1038/s41550-023-02101-x

Deep Synoptic Array science: Two fast radio burst sources in massive galaxy clusters

Authors: Liam Connor, Vikram Ravi, Morgan Catha, Ge Chen, Jakob T. Faber, James W. Lamb, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, Casey Law, Paul Rasmussen, Jack Sayers, Kritti Sharma, Myles B. Sherman, Jun Shi, Dana Simard, Jean Somalwar, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: The hot gas that constitutes the intracluster medium (ICM) has been studied at X-ray and millimeter/sub-millimeter wavelengths (Sunyaev-Zeldovich effect) for decades. Fast radio bursts (FRBs) offer an additional method of directly measuring the ICM and gas surrounding clusters, via observables such as dispersion measure (DM) and Faraday rotation measure (RM). We report the discovery of two FRB sou… ▽ More The hot gas that constitutes the intracluster medium (ICM) has been studied at X-ray and millimeter/sub-millimeter wavelengths (Sunyaev-Zeldovich effect) for decades. Fast radio bursts (FRBs) offer an additional method of directly measuring the ICM and gas surrounding clusters, via observables such as dispersion measure (DM) and Faraday rotation measure (RM). We report the discovery of two FRB sources detected with the Deep Synoptic Array (DSA-110) whose host galaxies belong to massive galaxy clusters. In both cases, the FRBs exhibit excess extragalactic DM, some of which likely originates in the ICM of their respective clusters. FRB 20220914A resides in the galaxy cluster Abell 2310 at z=0.1125 with a projected offset from the cluster center of 520 kpc. The host of a second source, FRB 20220509G, is an elliptical galaxy at z=0.0894 that belongs to the galaxy cluster Abell 2311 at projected offset of 870 kpc. These sources represent the first time an FRB has been localized to a galaxy cluster. We combine our FRB data with archival X-ray, SZ, and optical observations of these clusters in order to infer properties of the ICM, including a measurement of gas temperature from DM and ySZ of 0.8-3.9 keV. We then compare our results to massive cluster halos from the IllustrisTNG simulation. Finally, we describe how large samples of localized FRBs from future surveys will constrain the ICM, particularly beyond the virial radius of clusters. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Deep Synoptic Array science: A massive elliptical host among two galaxy-cluster fast radio bursts

Authors: Kritti Sharma, Jean Somalwar, Casey Law, Vikram Ravi, Morgan Catha, Ge Chen, Liam Connor, Jakob T. Faber, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, James W. Lamb, Paul Rasmussen, Myles B. Sherman, Jun Shi, Dana Simard, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: The stellar population environments associated with fast radio burst (FRB) sources provide important insights for developing their progenitor theories. We expand the diversity of known FRB host environments by reporting two FRBs in massive galaxy clusters discovered by the Deep Synoptic Array (DSA-110) during its commissioning observations. FRB 20220914A has been localized to a star-forming, late-… ▽ More The stellar population environments associated with fast radio burst (FRB) sources provide important insights for developing their progenitor theories. We expand the diversity of known FRB host environments by reporting two FRBs in massive galaxy clusters discovered by the Deep Synoptic Array (DSA-110) during its commissioning observations. FRB 20220914A has been localized to a star-forming, late-type galaxy at a redshift of 0.1139 with multiple starbursts at lookback times less than $\sim$3.5 Gyr in the Abell 2310 galaxy cluster. Although the host galaxy of FRB 20220914A is similar to typical FRB hosts, the FRB 20220509G host stands out as a quiescent, early-type galaxy at a redshift of 0.0894 in the Abell 2311 galaxy cluster. The discovery of FRBs in both late and early-type galaxies adds to the body of evidence that the FRB sources have multiple formation channels. Therefore, even though FRB hosts are typically star-forming, there must exist formation channels consistent with old stellar population in galaxies. The varied star formation histories of the two FRB hosts we report indicate a wide delay-time distribution of FRB progenitors. Future work in constraining the FRB delay-time distribution, using methods we develop herein, will prove crucial in determining the evolutionary histories of FRB sources. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 18 pages, 9 figures

CHIME/FRB Discovery of 25 Repeating Fast Radio Burst Sources

Authors: The CHIME/FRB Collaboration, :, Bridget C. Andersen, Kevin Bandura, Mohit Bhardwaj, P. J. Boyle, Charanjot Brar, Tomas Cassanelli, S. Chatterjee, Pragya Chawla, Amanda M. Cook, Alice P. Curtin, Matt Dobbs, Fengqiu Adam Dong, Jakob T. Faber, Mateus Fandino, Emmanuel Fonseca, B. M. Gaensler, Utkarsh Giri, Antonio Herrera-Martin, Alex S. Hill, Adaeze Ibik, Alexander Josephy, Jane F. Kaczmarek, Zarif Kader , et al. (35 additional authors not shown)

Abstract: We present the discovery of 25 new repeating fast radio burst (FRB) sources found among CHIME/FRB events detected between 2019 September 30 and 2021 May 1. The sources were found using a new clustering algorithm that looks for multiple events co-located on the sky having similar dispersion measures (DMs). The new repeaters have DMs ranging from $\sim$220 pc cm$^{-3}$ to $\sim$1700 pc cm$^{-3}$, an… ▽ More We present the discovery of 25 new repeating fast radio burst (FRB) sources found among CHIME/FRB events detected between 2019 September 30 and 2021 May 1. The sources were found using a new clustering algorithm that looks for multiple events co-located on the sky having similar dispersion measures (DMs). The new repeaters have DMs ranging from $\sim$220 pc cm$^{-3}$ to $\sim$1700 pc cm$^{-3}$, and include sources having exhibited as few as two bursts to as many as twelve. We report a statistically significant difference in both the DM and extragalactic DM (eDM) distributions between repeating and apparently nonrepeating sources, with repeaters having lower mean DM and eDM, and we discuss the implications. We find no clear bimodality between the repetition rates of repeaters and upper limits on repetition from apparently nonrepeating sources after correcting for sensitivity and exposure effects, although some active repeating sources stand out as anomalous. We measure the repeater fraction over time and find that it tends to an equilibrium of $2.6_{-2.6}^{+2.9}$% over our total time-on-sky thus far. We also report on 14 more sources which are promising repeating FRB candidates and which merit follow-up observations for confirmation. △ Less

Submitted 15 March, 2023; v1 submitted 20 January, 2023; originally announced January 2023.

Comments: ApJ in press. Comments are still welcome and follow-up observations are encouraged!

Deep Synoptic Array science: a 50 Mpc fast radio burst constrains the mass of the Milky Way circumgalactic medium

Authors: Vikram Ravi, Morgan Catha, Ge Chen, Liam Connor, James M. Cordes, Jakob T. Faber, James W. Lamb, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, Casey Law, Paul Rasmussen, Kritti Sharma, Myles B. Sherman, Jun Shi, Dana Simard, Jean J. Somalwar, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: We present the Deep Synoptic Array (DSA-110) discovery and interferometric localization of the so far non-repeating FRB 20220319D. The FRB originates in a young, rapidly star-forming barred spiral galaxy, IRAS 02044$+$7048, at just 50 Mpc. Although the NE2001 and YMW16 models for the Galactic interstellar-medium (ISM) contribution to the DM of FRB 20220319D exceed its total observed DM, we show th… ▽ More We present the Deep Synoptic Array (DSA-110) discovery and interferometric localization of the so far non-repeating FRB 20220319D. The FRB originates in a young, rapidly star-forming barred spiral galaxy, IRAS 02044$+$7048, at just 50 Mpc. Although the NE2001 and YMW16 models for the Galactic interstellar-medium (ISM) contribution to the DM of FRB 20220319D exceed its total observed DM, we show that uncertainties in these models accommodate an extragalactic origin for the burst. We derive a conservative upper limit on the DM contributed by the circumgalactic medium (CGM) of the Milky Way: the limit is either 28.7 pc cm$^{-3}$ and 47.3 pc cm$^{-3}$, depending on which of two pulsars nearby on the sky to FRB 20220319D is used to estimate the ISM DM. These limits both imply that the total Galactic CGM mass is $<10^{11}M_{\odot}$, and that the baryonic mass of the Milky Way is $\lesssim60\%$ of the cosmological average given the total halo mass. More stringent albeit less conservative constraints are possible when the DMs of pulsars in the distant globular cluster M53 are additionally considered. Although our constraints are sensitive to possible anisotropy in the CGM and to the assumed form of the radial-density profile, they are not subject to uncertainties in the chemical and thermal properties of the CGM. Our results strongly support scenarios commonly predicted by galaxy-formation simulations wherein feedback processes expel baryonic matter from the halos of galaxies like the Milky Way. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Comments: 17 pages, 8 figures, 3 tables, submitted to AAS Journals

Deep Synoptic Array science I: discovery of the host galaxy of FRB 20220912A

Authors: Vikram Ravi, Morgan Catha, Ge Chen, Liam Connor, Jakob T. Faber, James W. Lamb, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, Casey Law, Paul Rasmussen, Kritti Sharma, Myles B. Sherman, Jun Shi, Dana Simard, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli, Tomas Ahumada, Dillon Dong, Christoffer Fremling , et al. (3 additional authors not shown)

Abstract: We report the detection and interferometric localization of the repeating fast radio burst (FRB) source FRB 20220912A during commissioning observations with the Deep Synoptic Array (DSA-110). Two bursts were detected from FRB 20220912A, one each on 2022 October 18 and 2022 October 25. The best-fit position is (R.A. J2000, decl. J2000) = (23:09:04.9, +48:42:25.4), with a 90% confidence error ellips… ▽ More We report the detection and interferometric localization of the repeating fast radio burst (FRB) source FRB 20220912A during commissioning observations with the Deep Synoptic Array (DSA-110). Two bursts were detected from FRB 20220912A, one each on 2022 October 18 and 2022 October 25. The best-fit position is (R.A. J2000, decl. J2000) = (23:09:04.9, +48:42:25.4), with a 90% confidence error ellipse of $\pm2$ arcsec and $\pm1$ arcsec in right ascension and declination respectively. The two bursts have disparate polarization properties and temporal profiles. We find a Faraday rotation measure that is consistent with the low value of $+0.6$ rad m$^{-2}$ reported by CHIME/FRB. The DSA-110 localization overlaps with the galaxy PSO J347.2702+48.7066 at a redshift $z=0.0771$, which we identify as the likely host. PSO J347.2702$+$48.7066 has a stellar mass of approximately $10^{10}M_{\odot}$, modest internal dust extinction, and a star-formation rate likely in excess of $0.1\,M_{\odot}$ yr$^{-1}$. The host-galaxy contribution to the dispersion measure is likely $\lesssim50$ pc cm$^{-3}$. The FRB 20220912A source is therefore likely viewed along a tenuous plasma column through the host galaxy. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 10 pages, 7 figures, 2 tables, submitted to AAS Journals

Addition of tabulated equation of state and neutrino leakage support to IllinoisGRMHD

Authors: Leonardo R. Werneck, Zachariah B. Etienne, Ariadna Murguia-Berthier, Roland Haas, Federico Cipolletta, Scott C. Noble, Lorenzo Ennoggi, Federico G. Lopez Armengol, Bruno Giacomazzo, Thiago Assumpção, Joshua Faber, Tanmayee Gupte, Bernard J. Kelly, Julian H. Krolik

Abstract: We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which imp… ▽ More We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which implements a new, AMR-capable neutrino leakage scheme in the Einstein Toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by Carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of IllinoisGRMHD, as well as NRPyEOS and NRPyLeakage, is pedagogically documented in Jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the Einstein Toolkit. △ Less

Submitted 14 December, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 20 pages, 9 figures. v3 erratum: fixed typo in Eq. (21)

Journal ref: Phys. Rev. D 107, 044037 (2023)

Handing off the outcome of binary neutron star mergers for accurate and long-term post-merger simulations

Authors: Federico G. Lopez Armengol, Zachariah B. Etienne, Scott C. Noble, Bernard J. Kelly, Leonardo R. Werneck, Brendan Drachler, Manuela Campanelli, Federico Cipolletta, Yosef Zlochower, Ariadna Murguia-Berthier, Lorenzo Ennoggi, Mark Avara, Riccardo Ciolfi, Joshua Faber, Grace Fiacco, Bruno Giacomazzo, Tanmayee Gupte, Trung Ha, Julian H. Krolik, Vassilios Mewes, Richard O'Shaughnessy, Jesús M. Rueda-Becerril, Jeremy Schnittman

Abstract: We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at larg… ▽ More We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the Handoff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to HARM3D, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the Handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the \handoff approach. △ Less

Submitted 31 October, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Comments: 27 pages, 22 figures

Journal ref: Phys. Rev. D 106, 083015 (2022)

Multiband Detection of Repeating FRB 20180916B

Authors: Ketan R. Sand, Jakob Faber, Vishal Gajjar, Daniele Michilli, Bridget C. Andersen, Bhal Chandra Joshi, Sanjay Kudale, Maura Pilia, Bryan Brzycki, Tomas Cassanelli, Steve Croft, Biprateep Dey, Hoang John, Calvin Leung, Ryan Mckinven, Cherry Ng, Aaron B. Pearlman, Emily Petroff, Danny C. Price, Andrew Siemion, Kendrick Smith, Shriharsh P. Tendulkar

Abstract: We present a multiband study of FRB 20180916B, a repeating source with a 16.3 day periodicity. We report the detection of 4, 1 and 7 bursts from observations spanning 3 days using upgraded Giant Metrewave Radio Telescope (300-500 MHz), Canadian Hydrogen Intensity Mapping Experiment (400-800 MHz) and Green Bank Telescope (600-1000 MHz), respectively. We report the first-ever detection of the source… ▽ More We present a multiband study of FRB 20180916B, a repeating source with a 16.3 day periodicity. We report the detection of 4, 1 and 7 bursts from observations spanning 3 days using upgraded Giant Metrewave Radio Telescope (300-500 MHz), Canadian Hydrogen Intensity Mapping Experiment (400-800 MHz) and Green Bank Telescope (600-1000 MHz), respectively. We report the first-ever detection of the source in the 800-1000 MHz range along with one of the widest instantaneous bandwidth detection (200 MHz) at lower frequencies. We identify 30 $μ$s wide structures in one of the bursts at 800 MHz, making it the lowest frequency detection of such structures for this FRB thus far. There is also a clear indication of high activity of the source at a higher frequency during earlier phases of the activity cycle. We identify a gradual decrease in the rotation measure over two years and no significant variations in the dispersion measure. We derive useful conclusions about progenitor scenarios, energy distribution, emission mechanisms, and variation of downward drift rate of emission with frequency. Our results reinforce that multiband observations are an effective approach to study repeaters and even one-off events to better understand their varying activity and spectral anomalies. △ Less

Submitted 20 June, 2022; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: 20 pages, 10 figures

Journal ref: ApJ 932 98 (2022)

Absence of bursts between 4-8 GHz from FRB20200120E located in an M81 Globular Cluster

Authors: Vishal Gajjar, Daniele Michilli, Jakob T. Faber, Sabrina Berger, Steve Croft, Aaron B. Pearlman, Ketan R. Sand, Paul Scholz, Andrew P. V. Siemion

Abstract: We report the non-detection of dispersed bursts between 4 - 8 GHz from 2.5 hours of observations of FRB20200120E at 6 GHz using the Robert C. Byrd Green Bank Telescope. Our fluence limits are several times lower than the average burst fluences reported at 600 and 1400 MHz. We conclude that these non-detections are either due to high-frequency bursts being weaker and/or scintillation-induced modula… ▽ More We report the non-detection of dispersed bursts between 4 - 8 GHz from 2.5 hours of observations of FRB20200120E at 6 GHz using the Robert C. Byrd Green Bank Telescope. Our fluence limits are several times lower than the average burst fluences reported at 600 and 1400 MHz. We conclude that these non-detections are either due to high-frequency bursts being weaker and/or scintillation-induced modulated. It is also likely that our observations were non-concurrent with any activity window of FRB20200120E. △ Less

Submitted 16 July, 2021; originally announced July 2021.

Comments: To be published in RNAAS

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

Authors: Ariadna Murguia-Berthier, Scott C. Noble, Luke F. Roberts, Enrico Ramirez-Ruiz, Leonardo R. Werneck, Michael Kolacki, Zachariah B. Etienne, Mark Avara, Manuela Campanelli, Riccardo Ciolfi, Federico Cipolletta, Brendan Drachler, Lorenzo Ennoggi, Joshua Faber, Grace Fiacco, Bruno Giacomazzo, Tanmayee Gupte, Trung Ha, Bernard J. Kelly, Julian H. Krolik, Federico G. Lopez Armengol, Ben Margalit, Tim Moon, Richard O'Shaughnessy, Jesús M. Rueda-Becerril , et al. (3 additional authors not shown)

Abstract: The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merge… ▽ More The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind. △ Less

Submitted 31 January, 2022; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: 19 pages, 14 figures, published in ApJ

Re-Analysis of Breakthrough Listen Observations of FRB 121102: Polarization Properties of Eight New Spectrally Narrow Bursts

Authors: Jakob T. Faber, Vishal Gajjar, Andrew P. V. Siemion, Steve Croft, Daniel Czech, David DeBoer, Julia DeMarines, Jamie Drew, Howard Isaacson, Brian C. Lacki, Matt Lebofsky, David H. E. MacMahon, Cherry Ng, Imke de Pater, Danny C. Price, Sofia Z. Sheikh, Claire Webb, S. Pete Worden

Abstract: We report polarization properties for eight narrowband bursts from FRB 121102 that have been re-detected in a high-frequency (4-8 GHz) Breakthrough Listen observation with the Green Bank Telescope, originally taken on 2017 August 26. The bursts were found to exhibit nearly 100% linear polarization, Faraday rotation measures (RM) bordering 9.3$\times$10$^4$ rad-m$^{-2}$, and stable polarization pos… ▽ More We report polarization properties for eight narrowband bursts from FRB 121102 that have been re-detected in a high-frequency (4-8 GHz) Breakthrough Listen observation with the Green Bank Telescope, originally taken on 2017 August 26. The bursts were found to exhibit nearly 100% linear polarization, Faraday rotation measures (RM) bordering 9.3$\times$10$^4$ rad-m$^{-2}$, and stable polarization position angles (PA), all of which agree with burst properties previously reported for FRB 121102 at the same epoch. We confirm that these detections are indeed physical bursts with limited spectral occupancies and further support the use of sub-banded search techniques in FRB detection. △ Less

Submitted 19 January, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 4 pages, 1 table. Accepted to RNAAS

Numerical generation of vector potentials from specified magnetic fields

Authors: Zachary J. Silberman, Thomas R. Adams, Joshua A. Faber, Zachariah B. Etienne, Ian Ruchlin

Abstract: Many codes have been developed to study highly relativistic, magnetized flows around and inside compact objects. Depending on the adopted formalism, some of these codes evolve the vector potential $\mathbf{A}$, and others evolve the magnetic field $\mathbf{B}=\nabla\times\mathbf{A}$ directly. Given that these codes possess unique strengths, sometimes it is desirable to start a simulation using a c… ▽ More Many codes have been developed to study highly relativistic, magnetized flows around and inside compact objects. Depending on the adopted formalism, some of these codes evolve the vector potential $\mathbf{A}$, and others evolve the magnetic field $\mathbf{B}=\nabla\times\mathbf{A}$ directly. Given that these codes possess unique strengths, sometimes it is desirable to start a simulation using a code that evolves $\mathbf{B}$ and complete it using a code that evolves $\mathbf{A}$. Thus transferring the data from one code to another would require an inverse curl algorithm. This paper describes two new inverse curl techniques in the context of Cartesian numerical grids: a cell-by-cell method, which scales approximately linearly with the numerical grid, and a global linear algebra approach, which has worse scaling properties but is generally more robust (e.g., in the context of a magnetic field possessing some nonzero divergence). We demonstrate these algorithms successfully generate smooth vector potential configurations in challenging special and general relativistic contexts. △ Less

Submitted 27 March, 2018; originally announced March 2018.

Comments: 20 pages, 9 figures

The Influence of Neutrinos on r-Process Nucleosynthesis in the Ejecta of Black Hole-Neutron Star Mergers

Authors: Luke F. Roberts, Jonas Lippuner, Matthew D. Duez, Joshua A. Faber, Francois Foucart, James C. Lombardi Jr., Sandra Ning, Christian D. Ott, Marcelo Ponce

Abstract: During the merger of a black hole and a neutron star, baryonic mass can become unbound from the system. Because the ejected material is extremely neutron-rich, the r-process rapidly synthesizes heavy nuclides as the material expands and cools. In this work, we map general relativistic models of black hole-neutron star (BHNS) mergers into a Newtonian smoothed particle hydrodynamics (SPH) code and f… ▽ More During the merger of a black hole and a neutron star, baryonic mass can become unbound from the system. Because the ejected material is extremely neutron-rich, the r-process rapidly synthesizes heavy nuclides as the material expands and cools. In this work, we map general relativistic models of black hole-neutron star (BHNS) mergers into a Newtonian smoothed particle hydrodynamics (SPH) code and follow the evolution of the thermodynamics and morphology of the ejecta until the outflows become homologous. We investigate how the subsequent evolution depends on our mapping procedure and find that the results are robust. Using thermodynamic histories from the SPH particles, we then calculate the expected nucleosynthesis in these outflows while varying the level of neutrino irradiation coming from the postmerger accretion disk. We find that the ejected material robustly produces r-process nucleosynthesis even for unrealistically high neutrino luminosities, due to the rapid velocities of the outflow. Nonetheless, we find that neutrinos can have an impact on the detailed pat- tern of the r-process nucleosynthesis. Electron neutrinos are captured by neutrons to produce protons while neutron capture is occurring. The produced protons rapidly form low mass seed nuclei for the r-process. These low mass seeds are eventually incorporated into the first r-process peak at A~78, producing mainly Ge and Se. We consider the mechanism of this process in detail and discuss if it can impact galactic chemical evolution of the first peak r-process nuclei. △ Less

Submitted 6 April, 2017; v1 submitted 28 January, 2016; originally announced January 2016.

Comments: 14 pages, 12 figures, accepted for publication in MNRAS

An Efficient Radiative Cooling Approximation for Use in Hydrodynamic Simulations

Authors: James C. Lombardi Jr., William G. McInally, Joshua A. Faber

Abstract: To make relevant predictions about observable emission, hydrodynamical simulation codes must employ schemes that account for radiative losses, but the large dimensionality of accurate radiative transfer schemes is often prohibitive. Stamatellos and collaborators introduced a scheme for smoothed particle hydrodynamics (SPH) simulations based on the notion of polytropic pseudo-clouds that uses only… ▽ More To make relevant predictions about observable emission, hydrodynamical simulation codes must employ schemes that account for radiative losses, but the large dimensionality of accurate radiative transfer schemes is often prohibitive. Stamatellos and collaborators introduced a scheme for smoothed particle hydrodynamics (SPH) simulations based on the notion of polytropic pseudo-clouds that uses only local quantities to estimate cooling rates. The computational approach is extremely efficient and works well in cases close to spherical symmetry, such as in star formation problems. Unfortunately, the method, which takes the local gravitational potential as an input, can be inaccurate when applied to non-spherical configurations, limiting its usefulness when studying disks or stellar collisions, among other situations of interest. Here, we introduce the "pressure scale height method," which incorporates the fluid pressure scale height into the determination of column densities and cooling rates, and show that it produces more accurate results across a wide range of physical scenarios while retaining the computational efficiency of the original method. The tested models include spherical polytropes as well as disks with specified density and temperature profiles. We focus on applying our techniques within an SPH code, although our method can be implemented within any particle-based Lagrangian or grid-based Eulerian hydrodynamic scheme. Our new method may be applied in a broad range of situations, including within the realm of stellar interactions, collisions, and mergers. △ Less

Submitted 13 November, 2014; originally announced November 2014.

Comments: 12 pages, 7 figures, submitted to MNRAS

Binary Neutron Star Mergers

Authors: Joshua A. Faber, Frederic A. Rasio

Abstract: We review the current status of studies of the coalescence of binary neutron star systems. We begin with a discussion of the formation channels of merging binaries and we discuss the most recent theoretical predictions for merger rates. Next, we turn to the quasi-equilibrium formalisms that are used to study binaries prior to the merger phase and to generate initial data for fully dynamical simula… ▽ More We review the current status of studies of the coalescence of binary neutron star systems. We begin with a discussion of the formation channels of merging binaries and we discuss the most recent theoretical predictions for merger rates. Next, we turn to the quasi-equilibrium formalisms that are used to study binaries prior to the merger phase and to generate initial data for fully dynamical simulations. The quasi-equilibrium approximation has played a key role in developing our understanding of the physics of binary coalescence and, in particular, of the orbital instability processes that can drive binaries to merger at the end of their lifetimes. We then turn to the numerical techniques used in dynamical simulations, including relativistic formalisms, (magneto-)hydrodynamics, gravitational-wave extraction techniques, and nuclear microphysics treatments. This is followed by a summary of the simulations performed across the field to date, including the most recent results from both fully relativistic and microphysically detailed simulations. Finally, we discuss the likely directions for the field as we transition from the first to the second generation of gravitational-wave interferometers and while supercomputers reach the petascale frontier. △ Less

Submitted 17 April, 2012; originally announced April 2012.

Comments: 82 pages, 18 figures (reproduced with permission of the original authors); to appear in Living Reviews in Relativity

The Einstein Toolkit: A Community Computational Infrastructure for Relativistic Astrophysics

Authors: Frank Löffler, Joshua Faber, Eloisa Bentivegna, Tanja Bode, Peter Diener, Roland Haas, Ian Hinder, Bruno C. Mundim, Christian D. Ott, Erik Schnetter, Gabrielle Allen, Manuela Campanelli, Pablo Laguna

Abstract: We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, c… ▽ More We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 62 pages, 20 figures

Accretion disks around kicked black holes: Post-kick Dynamics

Authors: Marcelo Ponce, Joshua A. Faber, James C. Lombardi Jr

Abstract: Numerical calculations of merging black hole binaries indicate that asymmetric emission of gravitational radiation can kick the merged black hole at up to thousands of km/s, and a number of systems have been observed recently whose properties are consistent with an active galactic nucleus containing a supermassive black hole moving with substantial velocity with respect to its broader accretion di… ▽ More Numerical calculations of merging black hole binaries indicate that asymmetric emission of gravitational radiation can kick the merged black hole at up to thousands of km/s, and a number of systems have been observed recently whose properties are consistent with an active galactic nucleus containing a supermassive black hole moving with substantial velocity with respect to its broader accretion disk. We study here the effect of an impulsive kick delivered to a black hole on the dynamical evolution of its accretion disk using a smoothed particle hydrodynamics code, focusing attention on the role played by the kick angle with respect to the orbital angular momentum vector of the pre-kicked disk. We find that for more vertical kicks, for which the angle between the kick and the normal vector to the disk $θ\lesssim 30^\circ$, a gap remains present in the inner disk, in accordance with the prediction from an analytic collisionless Keplerian disk model, while for more oblique kicks with $θ\gtrsim 45^\circ$, matter rapidly accretes toward the black hole. There is a systematic trend for higher potential luminosities for more oblique kick angles for a given black hole mass, disk mass and kick velocity, and we find large amplitude oscillations in time in the case of a kick oriented $60^\circ$ from the vertical. △ Less

Submitted 19 December, 2011; v1 submitted 8 July, 2011; originally announced July 2011.

Comments: Movies available as 'ancillary files'. Accepted for publication in ApJ

Journal ref: ApJ 745 (2012) 71

Tidal break-up of binary stars at the Galactic center and its consequences

Authors: Fabio Antonini, Joshua Faber, Alessia Gualandris, David Merritt

Abstract: The tidal breakup of binary star systems by the supermassive black hole (SMBH) in the center of the galaxy has been suggested as the source of both the observed sample of hypervelocity stars (HVSs) in the halo of the Galaxy and the S-stars that remain in tight orbits around Sgr A*. Here, we use a post-Newtonian N-body code to study the dynamics of main-sequence binaries on highly elliptical boun… ▽ More The tidal breakup of binary star systems by the supermassive black hole (SMBH) in the center of the galaxy has been suggested as the source of both the observed sample of hypervelocity stars (HVSs) in the halo of the Galaxy and the S-stars that remain in tight orbits around Sgr A*. Here, we use a post-Newtonian N-body code to study the dynamics of main-sequence binaries on highly elliptical bound orbits whose periapses lie close to the SMBH, determining the properties of ejected and bound stars as well as collision products. Unlike previous studies, we follow binaries that remain bound for several revolutions around the SMBH, finding that in the case of relatively large periapses and highly inclined binaries the Kozai resonance can lead to large periodic oscillations in the internal binary eccentricity and inclination. Collisions and mergers of the binary elements are found to increase significantly for multiple orbits around the SMBH, while HVSs are primarily produced during a binary's first passage. This process can lead to stellar coalescence and eventually serve as an important source of young stars at the galactic center. △ Less

Submitted 28 December, 2009; v1 submitted 10 September, 2009; originally announced September 2009.

Comments: accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.713:90-104,2010

Fully General Relativistic Simulations of Black Hole-Neutron Star Mergers

Authors: Zachariah B. Etienne, Joshua A. Faber, Yuk Tung Liu, Stuart L. Shapiro, Keisuke Taniguchi, Thomas W. Baumgarte

Abstract: Black hole-neutron star (BHNS) binaries are expected to be among the leading sources of gravitational waves observable by ground-based detectors, and may be the progenitors of short-hard gamma ray bursts (SGRBs) as well. Here, we discuss our new fully general relativistic calculations of merging BHNS binaries, which use high-accuracy, low-eccentricity, conformal thin-sandwich configurations as i… ▽ More Black hole-neutron star (BHNS) binaries are expected to be among the leading sources of gravitational waves observable by ground-based detectors, and may be the progenitors of short-hard gamma ray bursts (SGRBs) as well. Here, we discuss our new fully general relativistic calculations of merging BHNS binaries, which use high-accuracy, low-eccentricity, conformal thin-sandwich configurations as initial data. Our evolutions are performed using the moving puncture method and include a fully relativistic, high-resolution shock-capturing hydrodynamics treatment. Focusing on systems in which the neutron star is irrotational and the black hole is nonspinning with a 3:1 mass ratio, we investigate the inspiral, merger, and disk formation in the system. We find that the vast majority of material is promptly accreted and no more than 3% of the neutron star's rest mass is ejected into a tenuous, gravitationally bound disk. We find similar results for mass ratios of 2:1 and 1:1, even when we reduce the NS compaction in the 2:1 mass ratio case. These ambient disks reach temperatures suitable for triggering SGRBs, but their masses may be too small to produce the required total energy output. We measure gravitational waveforms and compute the effective strain in frequency space, finding measurable differences between our waveforms and those produced by binary black hole mergers within the advanced LIGO band. These differences appear at frequencies corresponding to the emission that occurs when the NS is tidally disrupted and accreted by the black hole. The resulting information about the radius of the neutron star may be used to constrain the neutron star equation of state. △ Less

Submitted 16 July, 2008; v1 submitted 14 December, 2007; originally announced December 2007.

Comments: 22 pages, 14 figures, fixed a few typos

Journal ref: Phys.Rev.D77:084002,2008

Relativistic black hole-neutron star binaries in quasiequilibrium: effects of the black hole excision boundary condition

Authors: Keisuke Taniguchi, Thomas W. Baumgarte, Joshua A. Faber, Stuart L. Shapiro

Abstract: We construct new models of black hole-neutron star binaries in quasiequilibrium circular orbits by solving Einstein's constraint equations in the conformal thin-sandwich decomposition together with the relativistic equations of hydrostationary equilibrium. We adopt maximal slicing, assume spatial conformal flatness, and impose equilibrium boundary conditions on an excision surface (i.e., the app… ▽ More We construct new models of black hole-neutron star binaries in quasiequilibrium circular orbits by solving Einstein's constraint equations in the conformal thin-sandwich decomposition together with the relativistic equations of hydrostationary equilibrium. We adopt maximal slicing, assume spatial conformal flatness, and impose equilibrium boundary conditions on an excision surface (i.e., the apparent horizon) to model the black hole. In our previous treatment we adopted a "leading-order" approximation for a parameter related to the black-hole spin in these boundary conditions to construct approximately nonspinning black holes. Here we improve on the models by computing the black hole's quasilocal spin angular momentum and setting it to zero. As before, we adopt a polytropic equation of state with adiabatic index Gamma=2 and assume the neutron star to be irrotational. In addition to recomputing several sequences for comparison with our earlier results, we study a wider range of neutron star masses and binary mass ratios. To locate the innermost stable circular orbit we search for turning points along both the binding energy and total angular momentum curves for these sequences. Unlike for our previous approximate boundary condition, these two minima now coincide. We also identify the formation of cusps on the neutron star surface, indicating the onset of tidal disruption. Comparing these two critical binary separations for different mass ratios and neutron star compactions we distinguish those regions that will lead to a tidal disruption of the neutron star from those that will result in the plunge into the black hole of a neutron star more or less intact, albeit distorted by tidal forces. △ Less

Submitted 1 February, 2008; v1 submitted 26 October, 2007; originally announced October 2007.

Comments: Minor corrections, Fig.8 revised, 15 pages, 15 figures, published in Phys. Rev. D

Journal ref: Phys.Rev.D77:044003,2008

Relativistic hydrodynamics in the presence of puncture black holes

Authors: Joshua A. Faber, Thomas W. Baumgarte, Zachariah B. Etienne, Stuart L. Shapiro, Keisuke Taniguchi

Abstract: Many of the recent numerical simulations of binary black holes in vacuum adopt the moving puncture approach. This successful approach avoids the need to impose numerical excision of the black hole interior and is easy to implement. Here we wish to explore how well the same approach can be applied to moving black hole punctures in the presence of relativistic hydrodynamic matter. First, we evolve… ▽ More Many of the recent numerical simulations of binary black holes in vacuum adopt the moving puncture approach. This successful approach avoids the need to impose numerical excision of the black hole interior and is easy to implement. Here we wish to explore how well the same approach can be applied to moving black hole punctures in the presence of relativistic hydrodynamic matter. First, we evolve single black hole punctures in vacuum to calibrate our BSSN (Baumgarte-Shapiro-Shibata-Nakamura) implementation and to confirm that the numerical solution for the exterior spacetime is invariant to any ``junk'' (i.e., constraint-violating) initial data employed in the black hole interior. Then we focus on relativistic Bondi accretion onto a moving puncture Schwarzschild black hole as a numerical testbed for our high-resolution shock-capturing relativistic hydrodynamics scheme. We find that the hydrodynamical equations can be evolved successfully in the interior without imposing numerical excision. These results help motivate the adoption of the moving puncture approach to treat the binary black hole-neutron star problem using conformal thin-sandwich initial data. △ Less

Submitted 15 October, 2007; v1 submitted 17 August, 2007; originally announced August 2007.

Comments: 21 pages, 21 figures, RevTex, accepted by PRD; minor updates

Journal ref: Phys.Rev.D76:104021,2007

Filling the holes: Evolving excised binary black hole initial data with puncture techniques

Authors: Zachariah B. Etienne, Joshua A. Faber, Yuk Tung Liu, Stuart L. Shapiro, Thomas W. Baumgarte

Abstract: We follow the inspiral and merger of equal-mass black holes (BHs) by the moving puncture technique and demonstrate that both the exterior solution and the asymptotic gravitational waveforms are unchanged when the initial interior solution is replaced by constraint-violating ``junk'' initial data. We apply this result to evolve conformal thin-sandwich (CTS) binary BH initial data by filling their… ▽ More We follow the inspiral and merger of equal-mass black holes (BHs) by the moving puncture technique and demonstrate that both the exterior solution and the asymptotic gravitational waveforms are unchanged when the initial interior solution is replaced by constraint-violating ``junk'' initial data. We apply this result to evolve conformal thin-sandwich (CTS) binary BH initial data by filling their excised interiors with arbitrary, but smooth, initial data and evolving with standard puncture gauge choices. The waveforms generated for both puncture and filled-CTS initial data are remarkably similar, and there are only minor differences between irrotational and corotational CTS BH binaries. Even the interior solutions appear to evolve to the same constraint-satisfying solution at late times, independent of the initial data. △ Less

Submitted 27 September, 2007; v1 submitted 16 July, 2007; originally announced July 2007.

Comments: 5 pages, 5 figures, accepted by PRD Rapid Communications, RevTeX

Journal ref: Phys.Rev.D76:101503,2007

Quasiequilibrium black hole-neutron star binaries in general relativity

Authors: Keisuke Taniguchi, Thomas W. Baumgarte, Joshua A. Faber, Stuart L. Shapiro

Abstract: We construct quasiequilibrium sequences of black hole-neutron star binaries in general relativity. We solve Einstein's constraint equations in the conformal thin-sandwich formalism, subject to black hole boundary conditions imposed on the surface of an excised sphere, together with the relativistic equations of hydrostatic equilibrium. In contrast to our previous calculations we adopt a flat spa… ▽ More We construct quasiequilibrium sequences of black hole-neutron star binaries in general relativity. We solve Einstein's constraint equations in the conformal thin-sandwich formalism, subject to black hole boundary conditions imposed on the surface of an excised sphere, together with the relativistic equations of hydrostatic equilibrium. In contrast to our previous calculations we adopt a flat spatial background geometry and do not assume extreme mass ratios. We adopt a Gamma=2 polytropic equation of state and focus on irrotational neutron star configurations as well as approximately nonspinning black holes. We present numerical results for ratios of the black hole's irreducible mass to the neutron star's ADM mass in isolation of M_{irr}^{BH}/M_{ADM,0}^{NS} = 1, 2, 3, 5, and 10. We consider neutron stars of baryon rest mass M_B^{NS}/M_B^{max} = 83% and 56%, where M_B^{max} is the maximum allowed rest mass of a spherical star in isolation for our equation of state. For these sequences, we locate the onset of tidal disruption and, in cases with sufficiently large mass ratios and neutron star compactions, the innermost stable circular orbit. We compare with previous results for black hole-neutron star binaries and find excellent agreement with third-order post-Newtonian results, especially for large binary separations. We also use our results to estimate the energy spectrum of the outgoing gravitational radiation emitted during the inspiral phase for these binaries. △ Less

Submitted 24 October, 2007; v1 submitted 19 January, 2007; originally announced January 2007.

Comments: 17 pages, 15 figures, published in Phys. Rev. D

Journal ref: Phys.Rev.D75:084005,2007

Quasiequilibrium sequences of black-hole--neutron-star binaries in general relativity

Authors: Keisuke Taniguchi, Thomas W. Baumgarte, Joshua A. Faber, Stuart L. Shapiro

Abstract: We construct quasiequilibrium sequences of black hole-neutron star binaries for arbitrary mass ratios by solving the constraint equations of general relativity in the conformal thin-sandwich decomposition. We model the neutron star as a stationary polytrope satisfying the relativistic equations of hydrodynamics, and account for the black hole by imposing equilibrium boundary conditions on the su… ▽ More We construct quasiequilibrium sequences of black hole-neutron star binaries for arbitrary mass ratios by solving the constraint equations of general relativity in the conformal thin-sandwich decomposition. We model the neutron star as a stationary polytrope satisfying the relativistic equations of hydrodynamics, and account for the black hole by imposing equilibrium boundary conditions on the surface of an excised sphere (the apparent horizon). In this paper we focus on irrotational configurations, meaning that both the neutron star and the black hole are approximately nonspinning in an inertial frame. We present results for a binary with polytropic index n=1, mass ratio M_{irr}^{BH}/M_{B}^{NS}=5 and neutron star compaction M_{ADM,0}^{NS}/R_0=0.0879, where M_{irr}^{BH} is the irreducible mass of the black hole, M_{B}^{NS} the neutron star baryon rest-mass, and M_{ADM,0}^{NS} and R_0 the neutron star Arnowitt-Deser-Misner mass and areal radius in isolation, respectively. Our models represent valid solutions to Einstein's constraint equations and may therefore be employed as initial data for dynamical simulations of black hole-neutron star binaries. △ Less

Submitted 15 September, 2006; originally announced September 2006.

Comments: 5 pages, 1 figure, revtex4, published in Phys.Rev.D

Journal ref: Phys.Rev. D74 (2006) 041502

Black Hole-Neutron Star Binary Merger Calculations: GRB Progenitors and the Stability of Mass Transfer

Authors: Joshua A. Faber, Thomas W. Baumgarte, Stuart L. Shapiro, Keisuke Taniguchi, Frederic A. Rasio

Abstract: We have calculated the first dynamical evolutions of merging black hole-neutron star binaries that treat the combined spacetime in a nonperturbative general relativistic framework. Using the conformal flatness approximation, we have studied how the location of the tidal disruption radius with respect to the the black hole horizon and innermost stable circular orbit (ISCO) affects the qualitative… ▽ More We have calculated the first dynamical evolutions of merging black hole-neutron star binaries that treat the combined spacetime in a nonperturbative general relativistic framework. Using the conformal flatness approximation, we have studied how the location of the tidal disruption radius with respect to the the black hole horizon and innermost stable circular orbit (ISCO) affects the qualitative evolution of the system. Based on simple arguments, we show that for a binary mass ratio q>~0.24, tidal disruption occurs outside the ISCO, while the opposite is true for q<~0.24. When tidal disruption occurs sufficiently far outside the ISCO, mass is transferred unstably from the neutron star to the black hole, resulting in the complete disruption of the neutron star. When tidal disruption occurs slightly within the ISCO, we find that some of the mass forms an extremely hot disk around the black hole. The resulting configurations in this case are excellent candidates for the progenitors of short-hard gamma ray bursts. △ Less

Submitted 19 May, 2006; originally announced May 2006.

Comments: 8 pages, 5 figures, in Proceedings of the Albert Einstein Century International Conference, Paris, France, 2005, edited by Jean-Michel Alimi and Andre Fuzfa

General Relativistic Binary Merger Simulations and Short Gamma Ray Bursts

Authors: Joshua A. Faber, Thomas W. Baumgarte, Stuart L. Shapiro, Keisuke Taniguchi

Abstract: The recent localization of some short-hard gamma ray bursts (GRBs) in galaxies with low star formation rates has lent support to the suggestion that these events result from compact object binary mergers. We discuss how new simulations in general relativity are helping to identify the central engine of short-hard GRBs. Motivated by our latest relativistic black hole-neutron star merger calculati… ▽ More The recent localization of some short-hard gamma ray bursts (GRBs) in galaxies with low star formation rates has lent support to the suggestion that these events result from compact object binary mergers. We discuss how new simulations in general relativity are helping to identify the central engine of short-hard GRBs. Motivated by our latest relativistic black hole-neutron star merger calculations, we discuss a scenario in which these events may trigger short-hard GRBs, and compare this model to competing relativistic models involving binary neutron star mergers and the delayed collapse of hypermassive neutron stars. Distinguishing features of these models may help guide future GRB and gravitational wave observations to identify the nature of the sources. △ Less

Submitted 10 March, 2006; originally announced March 2006.

Comments: 5 Pages, 3 figures, accepted by ApJ Lett

Journal ref: Astrophys.J. 641 (2006) L93-L96

The Dynamical Evolution of Black Hole-Neutron Star Binaries in General Relativity: Simulations of Tidal Disruption

Authors: Joshua A. Faber, Thomas W. Baumgarte, Stuart L. Shapiro, Keisuke Taniguchi, Frederic A. Rasio

Abstract: We calculate the first dynamical evolutions of merging black hole-neutron star binaries that construct the combined black hole-neutron star spacetime in a general relativistic framework. We treat the metric in the conformal flatness approximation, and assume that the black hole mass is sufficiently large compared to that of the neutron star so that the black hole remains fixed in space. Using a… ▽ More We calculate the first dynamical evolutions of merging black hole-neutron star binaries that construct the combined black hole-neutron star spacetime in a general relativistic framework. We treat the metric in the conformal flatness approximation, and assume that the black hole mass is sufficiently large compared to that of the neutron star so that the black hole remains fixed in space. Using a spheroidal spectral methods solver, we solve the resulting field equations for a neutron star orbiting a Schwarzschild black hole. The matter is evolved using a relativistic, Lagrangian, smoothed particle hydrodynamics (SPH) treatment. We take as our initial data recent quasiequilibrium models for synchronized neutron star polytropes generated as solutions of the conformal thin-sandwich (CTS) decomposition of the Einstein field equations. We are able to construct from these models relaxed SPH configurations whose profiles show good agreement with CTS solutions. Our adiabatic evolution calculations for neutron stars with low compactness show that mass transfer, when it begins while the neutron star orbit is still outside the innermost stable circular orbit, is more unstable than is typically predicted by analytical formalisms. This dynamical mass loss is found to be the driving force in determining the subsequent evolution of the binary orbit and the neutron star, which typically disrupts completely within a few orbital periods. The majority of the mass transferred onto the black hole is accreted promptly; a significant fraction (~30%) of the mass is shed outward as well, some of which will become gravitationally unbound and ejected completely from the system. The remaining portion forms an accretion disk around the black hole, and could provide the energy source for short-duration gamma ray bursts. △ Less

Submitted 11 November, 2005; originally announced November 2005.

Comments: 32 pages, 16 figures, 2 tables, RevTeX, accepted by PRD

Journal ref: Phys.Rev.D73:024012,2006

Black hole-neutron star binaries in general relativity: effects of neutron star spin

Authors: Keisuke Taniguchi, Thomas W. Baumgarte, Joshua A. Faber, Stuart L. Shapiro

Abstract: We present new sequences of general relativistic, quasiequilibrium black hole-neutron star binaries. We solve for the gravitational field in the conformal thin-sandwich decomposition of Einstein's field equations, coupled to the equations of relativistic hydrostatic equilibrium for a perfect fluid. We account for the black hole by solving these equations in the background metric of a Schwarzschi… ▽ More We present new sequences of general relativistic, quasiequilibrium black hole-neutron star binaries. We solve for the gravitational field in the conformal thin-sandwich decomposition of Einstein's field equations, coupled to the equations of relativistic hydrostatic equilibrium for a perfect fluid. We account for the black hole by solving these equations in the background metric of a Schwarzschild black hole whose mass is much larger than that of the neutron star. The background metric is treated in Kerr-Schild as well as isotropic coordinates. For the neutron star, we assume a polytropic equation of state with adiabatic index Gamma=2, and solve for both irrotational and corotational configurations. By comparing the results of irrotational and synchronized configurations with the same background metric, we conclude that the effect of the rotation on the location of tidal break-up is only on the order of a few percent. The different choices in the background also lead to differences of order a few percent, which may be an indication of the level to which these configurations approximate quasiequilibrium. △ Less

Submitted 24 October, 2007; v1 submitted 20 May, 2005; originally announced May 2005.

Comments: 27 pages, 18 figures, published in Phys. Rev. D

Journal ref: Phys.Rev. D72 (2005) 044008

Relativistic SPH Calculations of Compact Binary Mergers

Authors: Frederic Rasio, Joshua Faber, Shiho Kobayashi, Pablo Laguna

Abstract: We review recent progress in understanding the hydrodynamics of compact binary mergers using relativistic smoothed particle hydrodynamics (SPH) codes. Recent results are discussed for both double neutron stars and black hole - neutron star binaries. The gravitational wave signals from these mergers contain a lot of information about the masses and radii of the stars and the equation of state of… ▽ More We review recent progress in understanding the hydrodynamics of compact binary mergers using relativistic smoothed particle hydrodynamics (SPH) codes. Recent results are discussed for both double neutron stars and black hole - neutron star binaries. The gravitational wave signals from these mergers contain a lot of information about the masses and radii of the stars and the equation of state of the neutron star fluid. For black hole - neutron star mergers, the final fate of the fluid and the gravitational wave emission depend strongly on the mass ratio and the black hole spin. △ Less

Submitted 28 February, 2005; originally announced March 2005.

Comments: 12 pages, to appear in proceedings of JGRG14 (Yukawa Institute for Theoretical Physics, Kyoto, Japan, Nov 29 - Dec 3, 2004)

Tidal Interactions and Disruptions of Giant Planets on Highly Eccentric Orbits

Authors: Joshua A. Faber, Frederic A. Rasio, Bart Willems

Abstract: We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be strippe… ▽ More We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be stripped from it, but in no case do we find enough energy transferred to the planet to lead to complete disruption. In light of the three new extrasolar planets discovered with periods shorter than two days, we argue that the shortest-period cases observed in the period-mass relation may be explained by a model whereby planets undergo strong tidal encounters with stars, after either being scattered by dynamical interactions into highly eccentric orbits, or tidally captured from nearly parabolic orbits. Although this scenario does provide a natural explanation for the edge found for planets at twice the Roche limit, it does not explain how such planets will survive the inevitable expansion that results from energy injection during tidal circularization. △ Less

Submitted 15 November, 2004; v1 submitted 15 July, 2004; originally announced July 2004.

Comments: 32 pages, 10 figures, accepted by Icarus; updated for referees suggestions and recent observations

Journal ref: Icarus 175 (2005) 248

Relativistic Calculations of Coalescing Binary Neutron Stars

Authors: Joshua Faber, Philippe Grandclément, Frederic Rasio

Abstract: We have designed and tested a new relativistic Lagrangian hydrodynamics code, which treats gravity in the conformally flat approximation to general relativity. We have tested the resulting code extensively, finding that it performs well for calculations of equilibrium single-star models, collapsing relativistic dust clouds, and quasi-circular orbits of equilibrium solutions. By adding in a radia… ▽ More We have designed and tested a new relativistic Lagrangian hydrodynamics code, which treats gravity in the conformally flat approximation to general relativity. We have tested the resulting code extensively, finding that it performs well for calculations of equilibrium single-star models, collapsing relativistic dust clouds, and quasi-circular orbits of equilibrium solutions. By adding in a radiation reaction treatment, we compute the full evolution of a coalescing binary neutron star system. We find that the amount of mass ejected from the system, much less than a percent, is greatly reduced by the inclusion of relativistic gravitation. The gravity wave energy spectrum shows a clear divergence away from the Newtonian point-mass form, consistent with the form derived from relativistic quasi-equilibrium fluid sequences. △ Less

Submitted 2 March, 2004; originally announced March 2004.

Comments: 7 pages, proceedings of the ICGC 2004 meeting, to appear in Pramana

Journal ref: Pramana63:837-844,2004

Mergers of Irrotational Neutron Star Binaries in Conformally Flat Gravity

Authors: Joshua A. Faber, Philippe Grandclément, Frederic A. Rasio

Abstract: We present the first results from our new general relativistic, Lagrangian hydrodynamics code, which treats gravity in the conformally flat (CF) limit. The evolution of fluid configurations is described using smoothed particle hydrodynamics (SPH), and the elliptic field equations of the CF formalism are solved using spectral methodsin spherical coordinates. The code was tested on models for whic… ▽ More We present the first results from our new general relativistic, Lagrangian hydrodynamics code, which treats gravity in the conformally flat (CF) limit. The evolution of fluid configurations is described using smoothed particle hydrodynamics (SPH), and the elliptic field equations of the CF formalism are solved using spectral methodsin spherical coordinates. The code was tested on models for which the CF limit is exact, finding good agreement with the classical Oppenheimer-Volkov solution for a relativistic static spherical star as well as the exact semi-analytic solution for a collapsing spherical dust cloud. By computing the evolution of quasi-equilibrium neutron star binary configurations in the absence of gravitational radiation backreaction, we have confirmed that these configurations can remain dynamically stable all the way to the development of a cusp. With an approximate treatment of radiation reaction, we have calculated the complete merger of an irrotational binary configuration from the innermost point on an equilibrium sequence through merger and remnant formation and ringdown, finding good agreement withprevious relativistic calculations. In particular, we find that mass loss is highly suppressed by relativistic effects, but that, for a reasonably stiff neutron star equation of state, the remnant is initially stable against gravitational collapse because of its strong differential rotation. The gravity wave signal derived from our numerical calculation has an energy spectrum which matches extremely well with estimates based solely on quasi-equilibrium results, deviating from the Newtonian power-law form at frequencies below 1 kHz, i.e., within the reach of advanced interferometric detectors. △ Less

Submitted 2 April, 2004; v1 submitted 21 December, 2003; originally announced December 2003.

Comments: 44 pages, 21 figures, 2 tables, RevTeX, submitted to PRD; revised and updated

Journal ref: Phys.Rev. D69 (2004) 124036

The speed of gravity has not been measured from time delays

Authors: Joshua A. Faber

Abstract: The recent passage of Jupiter by the quasar QSO J0842+1835 at a separation of 3.7 arcminutes on September 8, 2002, combined with recent advances in interferometric radio timing, has allowed for the first measurement of higher-order post-Newtonian terms in the Shapiro time delay which depend linearly on the velocity of the gravitating body. Claims have been made that these measurements also allow… ▽ More The recent passage of Jupiter by the quasar QSO J0842+1835 at a separation of 3.7 arcminutes on September 8, 2002, combined with recent advances in interferometric radio timing, has allowed for the first measurement of higher-order post-Newtonian terms in the Shapiro time delay which depend linearly on the velocity of the gravitating body. Claims have been made that these measurements also allow for the measurement of the propagation speed of the gravitational force. This conclusion disagrees with recent calculations done in the parameterized post-Newtonian (PPN) model, which find no dependence of the velocity-dependent terms in the time delay on the speed of gravity to the stated order. Here, to test out these claims and counterclaims, we calculate the time delay in the limit of an instantaneous gravitational force, and find that the velocity-dependent terms are in complete agreement with previous PPN calculations, with no dependence on the speed of gravity. We conclude that the speed of gravity cannot be determined by measuring these terms in the Shapiro time delay, and suggest a reason why other groups mistakenly came to the opposite conclusion. △ Less

Submitted 14 March, 2003; originally announced March 2003.

Comments: 5 pages, no figures, RevTeX, submitted to ApJ Lett

Measuring Neutron Star Radii with Gravitational Wave Detectors

Authors: Joshua A. Faber, Philippe Grandclement, Frederic A. Rasio, Keisuke Taniguchi

Abstract: Coalescing binary neutron stars (NS) are expected to be an important source of gravitational waves (GW) detectable by laser interferometers. We present here a simple method for determining the compactness ratio M/R of NS based on the observed deviation of the GW energy spectrum from point-mass behavior at the end of an inspiral event. Our method is based on the properties of quasi-equilibrium bi… ▽ More Coalescing binary neutron stars (NS) are expected to be an important source of gravitational waves (GW) detectable by laser interferometers. We present here a simple method for determining the compactness ratio M/R of NS based on the observed deviation of the GW energy spectrum from point-mass behavior at the end of an inspiral event. Our method is based on the properties of quasi-equilibrium binary NS sequences and does not require the computation of the full GW signal h(t). Combined with the measurement of the NS masses from the GW signal during inspiral, the determination of M/R will allow very strong constraints to be placed on the equation of state of nuclear matter at high densities. △ Less

Submitted 10 October, 2002; v1 submitted 23 April, 2002; originally announced April 2002.

Comments: Minor revisions and ref. updates, RevTeX4, 6 pages, 2 figures, accepted by Phys. Rev. Lett

Journal ref: Phys.Rev.Lett. 89 (2002) 231102

Post-Newtonian SPH calculations of binary neutron star coalescence. III. Irrotational systems and gravitational wave spectra

Authors: Joshua A. Faber, Frederic A. Rasio

Abstract: Using our new post-Newtonian (PN) smoothed particle hydrodynamics (SPH) code, we have studied numerically the mergers of neutron star binaries with irrotational initial configurations. Here we describe a new method for constructing numerically accurate initial conditions for irrotational binary systems with circular orbits in PN gravity. We then compute the 3D hydrodynamic evolution of these sys… ▽ More Using our new post-Newtonian (PN) smoothed particle hydrodynamics (SPH) code, we have studied numerically the mergers of neutron star binaries with irrotational initial configurations. Here we describe a new method for constructing numerically accurate initial conditions for irrotational binary systems with circular orbits in PN gravity. We then compute the 3D hydrodynamic evolution of these systems until the two stars have completely merged, and we determine the corresponding GW signals. We present results for systems with different binary mass ratios, and for neutron stars represented by polytropes with $Γ=2$ or $Γ=3$. Compared to mergers of corotating binaries, we find that irrotational binary mergers produce similar peak GW luminosities, but they shed almost no mass at all to large distances. The dependence of the GW signal on numerical resolution for calculations performed with N>10^5 SPH particles is extremely weak, and we find excellent agreement between runs utilizing N=10^5 and N=10^6 SPH particles (the largest SPH calculation ever performed to study such irrotational binary mergers). We also compute GW energy spectra based on all calculations reported here and in our previous works. We find that PN effects lead to clearly identifiable features in the GW energy spectrum of binary neutron star mergers, which may yield important information about the nuclear equation of state at extreme densities. △ Less

Submitted 12 January, 2002; originally announced January 2002.

Comments: 34 pages, 19 figures, submitted to PRD

Journal ref: Phys.Rev. D65 (2002) 084042

Hydrodynamics of Neutron Star Mergers

Authors: Joshua A. Faber, Frederic A. Rasio

Abstract: The final burst of gravitational radiation emitted by coalescing binary neutron stars carries direct information about the neutron star fluid, and, in particular, about the equation of state of nuclear matter at extreme densities. The final merger may also be accompanied by a detectable electromagnetic signal, such as a gamma-ray burst. In this paper, we summarize the results of theoretical work… ▽ More The final burst of gravitational radiation emitted by coalescing binary neutron stars carries direct information about the neutron star fluid, and, in particular, about the equation of state of nuclear matter at extreme densities. The final merger may also be accompanied by a detectable electromagnetic signal, such as a gamma-ray burst. In this paper, we summarize the results of theoretical work done over the past decade that has led to a detailed understanding of this hydrodynamic merger process for two neutron stars, and we discuss the prospects for the detection and physical interpretation of the gravity wave signals by ground-based interferometers such as LIGO. We also present results from our latest post-Newtonian SPH calculations of binary neutron star coalescence, using up to 10^6 SPH particles to compute with higher spatial resolution than ever before the merger of an initially irrotational system. We discuss the detectability of our calculated gravity wave signals based on power spectra. △ Less

Submitted 17 January, 2001; originally announced January 2001.

Comments: 13 pages, 5 figures, to appear in "Astrophysical Sources of Gravitational Radiation", edited by J. M. Centrella (AIP Press)

Evolution of Stellar Collision Products in Globular Clusters - II. Off-axis Collision

Authors: Alison Sills, Joshua A. Faber, James C. Lombardi, Jr., Frederic A. Rasio, Aaron Warren

Abstract: We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitte… ▽ More We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitted much higher numerical resolution of the hydrodynamics. We also present evolutionary tracks for the resulting collision products, which emerge as rapidly rotating blue stragglers. The rotating collision products are brighter, bluer and remain on the main sequence longer than their non-rotating counterparts. In addition, they retain their rapid rotation rates throughout their main sequence lifetime. Rotationally-induced mixing strongly affects the evolution of the collision products, although it is not sufficient to mix the entire star. We discuss the implications of these results for studies of blue straggler populations in clusters. This work shows that off-axis collision products cannot become blue stragglers unless they lose a large fraction of their initial angular momentum. The mechanism for this loss is not apparent, although some possibilities are discussed. △ Less

Submitted 16 August, 2000; originally announced August 2000.

Comments: 25 pages incl. 9 figures (one in colour). Submitted to ApJ

Post-Newtonian SPH calculations of binary neutron star coalescence. II. Binary mass ratio, equation of state, and spin dependence

Authors: Joshua A. Faber, Frederic A. Rasio, Justin B. Manor

Abstract: Using our new Post-Newtonian SPH (smoothed particle hydrodynamics) code, we study the final coalescence and merging of neutron star (NS) binaries. We vary the stiffness of the equation of state (EOS) as well as the initial binary mass ratio and stellar spins. Results are compared to those of Newtonian calculations, with and without the inclusion of the gravitational radiation reaction. We find a… ▽ More Using our new Post-Newtonian SPH (smoothed particle hydrodynamics) code, we study the final coalescence and merging of neutron star (NS) binaries. We vary the stiffness of the equation of state (EOS) as well as the initial binary mass ratio and stellar spins. Results are compared to those of Newtonian calculations, with and without the inclusion of the gravitational radiation reaction. We find a much steeper decrease in the gravity wave peak strain and luminosity with decreasing mass ratio than would be predicted by simple point-mass formulae. For NS with softer EOS (which we model as simple $Γ=2$ polytropes) we find a stronger gravity wave emission, with a different morphology than for stiffer EOS (modeled as $Γ=3$ polytropes as in our previous work). We also calculate the coalescence of NS binaries with an irrotational initial condition, and find that the gravity wave signal is relatively suppressed compared to the synchronized case, but shows a very significant second peak of emission. Mass shedding is also greatly reduced, and occurs via a different mechanism than in the synchronized case. We discuss the implications of our results for gravity wave astronomy with laser interferometers such as LIGO, and for theoretical models of gamma-ray bursts (GRBs) based on NS mergers. △ Less

Submitted 15 October, 2000; v1 submitted 21 June, 2000; originally announced June 2000.

Comments: RevTeX, 38 pages, 24 figures, Minor Corrections, to appear in Phys. Rev. D

Journal ref: Phys.Rev. D63 (2001) 044012

Post-Newtonian SPH calculations of binary neutron star coalescence. I. Method and first results

Authors: Joshua A. Faber, Frederic A. Rasio

Abstract: We present the first results from our Post-Newtonian (PN) Smoothed Particle Hydrodynamics (SPH) code, which has been used to study the coalescence of binary neutron star (NS) systems. The Lagrangian particle-based code incorporates consistently all lowest-order (1PN) relativistic effects, as well as gravitational radiation reaction, the lowest-order dissipative term in general relativity. We tes… ▽ More We present the first results from our Post-Newtonian (PN) Smoothed Particle Hydrodynamics (SPH) code, which has been used to study the coalescence of binary neutron star (NS) systems. The Lagrangian particle-based code incorporates consistently all lowest-order (1PN) relativistic effects, as well as gravitational radiation reaction, the lowest-order dissipative term in general relativity. We test our code on sequences of single NS models of varying compactness, and we discuss ways to make PN simulations more relevant to realistic NS models. We also present a PN SPH relaxation procedure for constructing equilibrium models of synchronized binaries, and we use these equilibrium models as initial conditions for our dynamical calculations of binary coalescence. Though unphysical, since tidal synchronization is not expected in NS binaries, these initial conditions allow us to compare our PN work with previous Newtonian results. We compare calculations with and without 1PN effects, for NS with stiff equations of state, modeled as polytropes with $Γ=3$. We find that 1PN effects can play a major role in the coalescence, accelerating the final inspiral and causing a significant misalignment in the binary just prior to final merging. In addition, the character of the gravitational wave signal is altered dramatically, showing strong modulation of the exponentially decaying waveform near the end of the merger. We also discuss briefly the implications of our results for models of gamma-ray bursts at cosmological distances. △ Less

Submitted 13 April, 2000; v1 submitted 22 December, 1999; originally announced December 1999.

Comments: RevTeX, 37 pages, 17 figures, to appear in Phys. Rev. D, minor corrections only

Journal ref: Phys.Rev. D62 (2000) 064012