Showing 1–50 of 54 results for author: Andrews, J J

Formation of Wind-Fed Black Hole High-mass X-ray Binaries: The Role of Roche-lobe-Overflow Post Black-Hole Formation

Authors: Zepei Xing, Tassos Fragos, Emmanouil Zapartas, Tom M. Kwan, Lixin Dai, Ilya Mandel, Matthias U. Kruckow, Max Briel, Jeff J. Andrews, Simone S. Bavera, Seth Gossage, Konstantinos Kovlakas, Kyle A. Rocha, Meng Sun, Philipp M. Srivastava

Abstract: The three dynamically confirmed wind-fed black hole high-mass X-ray binaries (BH-HMXBs) are suggested to all contain a highly spinning black hole (BH). However, based on the theories of efficient angular momentum transport inside the stars, we expect that the first-born BHs in binary systems should have low spins, which is consistent with gravitational-wave observations. As a result, the origin of… ▽ More The three dynamically confirmed wind-fed black hole high-mass X-ray binaries (BH-HMXBs) are suggested to all contain a highly spinning black hole (BH). However, based on the theories of efficient angular momentum transport inside the stars, we expect that the first-born BHs in binary systems should have low spins, which is consistent with gravitational-wave observations. As a result, the origin of the high BH spins measured in wind-fed BH-HMXBs remains a mystery. In this paper, we conduct a binary population synthesis study on wind-fed BH-HMXBs at solar metallicity with the use of the newly developed code POSYDON, considering three scenarios for BH accretion: Eddington-limited, moderately super-Eddington, and fully conservative accretion. Taking into account the conditions for accretion-disk formation, we find that regardless of the accretion model, these systems are more likely to have already experienced a phase of Roche-lobe overflow after the BH formation. To account for the extreme BH spins, highly conservative accretion onto BHs is required, when assuming the accreted material carries the specific angular momentum at the innermost stable orbit. Besides, in our simulations we found that the systems with donor stars within the mass range of $10-20\,M_{\odot}$ are prevalent, posing a challenge in explaining simultaneously all observed properties of the BH-HMXB in our Galaxy, Cygnus X-1, and potentially hinting that the accretion efficiency onto non-degenerate stars, before the formation of the BH, is also more conservative than assumed in our simulations. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: Submitted to A&A. Comments are welcome!

A Lack of Mass-Gap Compact Object Binaries in APOGEE

Authors: Meir Schochet, Jamie Tayar, Jeff J. Andrews

Abstract: Depending principally on mass, the compact object remnant left behind after a star's life may be a white dwarf (WD), neutron star (NS), or black hole (BH). While we have large samples of each of these remnants, we lack knowledge of the exact conditions separating these outcomes. The boundary between low-mass BHs and massive NSs is particularly poorly understood, as few objects between 2-5… ▽ More Depending principally on mass, the compact object remnant left behind after a star's life may be a white dwarf (WD), neutron star (NS), or black hole (BH). While we have large samples of each of these remnants, we lack knowledge of the exact conditions separating these outcomes. The boundary between low-mass BHs and massive NSs is particularly poorly understood, as few objects between 2-5 $M_\odot$ are known. To probe this regime, we search the APOGEE DR17 dataset of 657,000 unique stars for binary systems with one stellar component that exhibit large radial velocity shifts over multiple observations. We identify 4751 likely binary systems, and estimate a minimum mass for each system's "invisible companion" under the assumption of tidal synchronization. Two systems have companion masses $\gtrsim$ 2 $M_\odot$, although we conclude that neither are good candidates for possessing a mass-gap NS or BH companions. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 4 pages, 6 figures, accepted for publication in RNAAS

Theia 456: Tidally Shredding an Open Cluster

Authors: Kyle R. Tregoning, Jeff J. Andrews, Marcel A. Agüeros, Phillip A. Cargile, Julio Chanamé, Jason L. Curtis, Simon C. Schuler

Abstract: The application of clustering algorithms to the Gaia astrometric catalog has revolutionized our census of stellar populations in the Milky Way, including the discovery of many new, dispersed structures. We focus on one such structure, Theia 456 (COIN-Gaia-13), a loosely bound collection of ~320 stars spanning ~120 pc that has previously been shown to exhibit kinematic, chemical, and gyrochronal co… ▽ More The application of clustering algorithms to the Gaia astrometric catalog has revolutionized our census of stellar populations in the Milky Way, including the discovery of many new, dispersed structures. We focus on one such structure, Theia 456 (COIN-Gaia-13), a loosely bound collection of ~320 stars spanning ~120 pc that has previously been shown to exhibit kinematic, chemical, and gyrochronal coherency, indicating a common origin. We obtain follow-up radial velocities and supplement these with Gaia astrometry to perform an in-depth dynamical analysis of Theia 456. By integrating stellar orbits through a Milky Way potential, we find the currently dispersed structure coalesced into a small cluster in the past. Via Bayesian modeling, we derive a kinematic age of 245 +/- 3 Myr (statistical), a half-mass radius of 9 +/- 2 pc, and an initial velocity dispersion of 0.14 +/- 0.02 km/s. Our results are entirely independent of model isochrones, details of stellar evolution, and internal cluster dynamics, and the statistical precision in our age derivation rivals that of the most precise age-dating techniques known today, though our imperfect knowledge of the Milky Way potential and simple spherical model for Theia 456 at birth add additional uncertainties. Using posterior predictive checking, we confirm these results are robust under reasonable variations to the Milky Way potential. Such low density structures that are disrupted by the Galactic tides before virializing may be ubiquitous, signifying that Theia 456 is a valuable benchmark for studying the dynamical history of stellar populations in the Milky Way. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 22 pages, 15 figures, 2 tables. Submitted to AAS journals. Comments welcome

To Be or not to Be: the role of rotation in modeling Galactic Be X-ray Binaries

Authors: Kyle Akira Rocha, Vicky Kalogera, Zoheyr Doctor, Jeff J. Andrews, Meng Sun, Seth Gossage, Simone S. Bavera, Tassos Fragos, Konstantinos Kovlakas, Matthias U. Kruckow, Devina Misra, Philipp M. Srivastava, Zepei Xing, Emmanouil Zapartas

Abstract: Be X-ray binaries (Be-XRBs) are crucial in understanding high-mass X-ray binaries, featuring a rapidly rotating Be star and a neutron star companion in an eccentric orbit, intermittently accreting material from the Be star's decretion disk. Originating from binary stellar evolution, Be-XRBs are of significant interest to binary population synthesis (BPS) studies, encapsulating the physics of super… ▽ More Be X-ray binaries (Be-XRBs) are crucial in understanding high-mass X-ray binaries, featuring a rapidly rotating Be star and a neutron star companion in an eccentric orbit, intermittently accreting material from the Be star's decretion disk. Originating from binary stellar evolution, Be-XRBs are of significant interest to binary population synthesis (BPS) studies, encapsulating the physics of supernovae, common envelope, and mass transfer (MT). Using the POSYDON BPS code, employing pre-computed grids of detailed binary stellar evolution models, we investigate the Galactic Be-XRB population. POSYDON incorporates stellar rotation self-consistently during MT phases, enabling a detailed examination of the rotational distribution of Be stars. Our fiducial BPS and Be-XRB model align well with the orbital properties of Galactic Be-XRBs, emphasizing the role of rotational constraints. Our modeling reveals a bimodal rotational distribution of Be-XRB-like systems, in excellent agreement with literature values. All Be-XRBs undergo an MT phase before the first compact object forms, with over half experiencing a second MT phase from a stripped helium companion (Case BB). Computing rotationally-limited MT efficiencies and applying them to our population, we find that the majority of Be-XRBs have undergone highly non-conservative MT (beta ~ 0.15). Our study underscores the importance of detailed angular momentum modeling during MT in interpreting Be-XRB populations, emphasizing this population as a key probe for the stability and efficiency of MT in interacting binaries. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 20 pages, 10 figures, Submitted to ApJ

Total r-process Yields of Milky Way Neutron Star Mergers

Authors: Erika M. Holmbeck, Jeff J. Andrews

Abstract: While it is now known that double neutron star binary systems (DNSs) are copious producers of heavy elements, there remains much speculation about whether they are the sole or even principal site of rapid neutron-capture (r-process) nucleosynthesis, one of the primary ways in which heavy elements are produced. The occurrence rates, delay times, and galactic environments of DNSs hold sway over esti… ▽ More While it is now known that double neutron star binary systems (DNSs) are copious producers of heavy elements, there remains much speculation about whether they are the sole or even principal site of rapid neutron-capture (r-process) nucleosynthesis, one of the primary ways in which heavy elements are produced. The occurrence rates, delay times, and galactic environments of DNSs hold sway over estimating their total contribution to the elemental abundances in the Solar system and the Galaxy. Furthermore, the expected elemental yield for DNSs may depend on the merger parameters themselves -- such as their stellar masses and radii -- which is not currently considered in many galactic chemical evolution models. Using the characteristics of the observed sample of DNSs in the Milky Way as a guide, we predict the expected nucleosynthetic yields that a population of DNSs would produce upon merger, and we compare that nucleosynthetic signature to the heavy-element abundance pattern of the Solar system elements. We find that with our current models, the present DNS population favors production of the lighter r-process elements, while underproducing the heaviest elements relative to the Solar system. This inconsistency could imply an additional site for the heaviest elements or a population of DNSs much different from that observed today. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 12 pages, 6 figures, 2 tables

Exploring the nature of ultra-luminous X-ray sources across stellar population ages using detailed binary evolution calculations

Authors: Devina Misra, Konstantinos Kovlakas, Tassos Fragos, Jeff J. Andrews, Simone S. Bavera, Emmanouil Zapartas, Zepei Xing, Aaron Dotter, Kyle Akira Rocha, Philipp M. Srivastava, Meng Sun

Abstract: Ultra-luminous X-ray sources (ULXs) are sources observed to exceed the Eddington limit of a stellar-mass black hole (BH). A fraction of ULX sources show X-ray pulses which are evidence for accreting neutron stars (NSs). Theoretical studies have suggested that NSs dominate the compact objects of intrinsic ULXs, even though the majority of observed sample is non-pulsating, implying that X-ray pulses… ▽ More Ultra-luminous X-ray sources (ULXs) are sources observed to exceed the Eddington limit of a stellar-mass black hole (BH). A fraction of ULX sources show X-ray pulses which are evidence for accreting neutron stars (NSs). Theoretical studies have suggested that NSs dominate the compact objects of intrinsic ULXs, even though the majority of observed sample is non-pulsating, implying that X-ray pulses from many NS ULXs are unobservable. We use POSYDON to generate and study X-ray binary populations spanning starburst ages 5 to 1000Myr. Following theoretical predictions for the alignment of the NS spin axis with the accretion disc, we estimate the required accreted mass in ULXs so that the alignment suppresses observable X-ray pulses. While the properties of ULXs are sensitive to model assumptions, there are certain trends that the populations follow. Young and old stellar populations are dominated by BH and NS accretors, respectively. The donors go from massive H-rich main-sequence (MS) stars in young populations (<100Myr) to low-mass post-MS H-rich stars in older populations (>100Myr), with stripped He-rich giant stars dominating the populations at around 100Myr. In addition, we find that NS ULXs exhibit stronger geometrical beaming than BH ULXs, leading to an under-representation of NS accretors in observed populations. Coupled with our finding that X-ray pulses are suppressed in at least 60% of the NS ULXs, we suggest that the observed fraction of ULXs with detectable X-ray pulses is very small, in agreement with observations. This study investigates the effects of age on ULXs as well as the effects of different model assumptions on ULX demographics. We show that geometrical beaming and the mass-accretion phase are critical aspects of understanding ULX observations. Our results suggest that even though most ULXs have accreting NSs, those with observable X-ray pulses would be very few. △ Less

Submitted 20 December, 2023; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 22 pages, 12 figures, Accepted by A&A. Fixed typos and updated references. Referee's comments were addressed

Journal ref: A&A 682, A69 (2024)

From ZAMS to Merger: Detailed Binary Evolution Models of Coalescing Neutron Star-Black Hole Systems at Solar Metallicity

Authors: Zepei Xing, Simone S. Bavera, Tassos Fragos, Matthias U. Kruckow, Jaim Román-Garza, Jeff J. Andrews, Aaron Dotter, Konstantinos Kovlakas, Devina Misra, Philipp M. Srivastava, Kyle A. Rocha, Meng Sun, Emmanouil Zapartas

Abstract: Neutron star $-$ black hole (NSBH) merger events bring us new opportunities to constrain theories of stellar and binary evolution, and understand the nature of compact objects. In this work, we investigate the formation of merging NSBH binaries at solar metallicity by performing a binary population synthesis study of merging NSBH binaries with the newly developed code POSYDON. The latter incorpora… ▽ More Neutron star $-$ black hole (NSBH) merger events bring us new opportunities to constrain theories of stellar and binary evolution, and understand the nature of compact objects. In this work, we investigate the formation of merging NSBH binaries at solar metallicity by performing a binary population synthesis study of merging NSBH binaries with the newly developed code POSYDON. The latter incorporates extensive grids of detailed single and binary evolution models, covering the entire evolution of a double compact object progenitor. We explore the evolution of NSBHs originating from different formation channels, which in some cases differ from earlier studies performed with rapid binary population synthesis codes. Then, we present the population properties of merging NSBH systems and their progenitors such as component masses, orbital features, and BH spins, and investigate the model uncertainties in our treatment of common envelope (CE) evolution and core-collapse process. We find that at solar metallicity, under the default model assumptions, most of the merging NSBHs have BH masses in a range of $3-11\,M{_\odot}$ and chirp masses within $1.5-4\,M{_\odot}$. Independently of our model variations, the BH always forms first with dimensionless spin parameter $\lesssim 0.2$, which is correlated to the initial binary orbital period. Some BHs can subsequently spin up moderately ($χ_{\rm BH} \lesssim 0.4$) due to mass transfer, which we assume to be Eddington limited. Binaries that experienced CE evolution rarely demonstrate large tilt angles. Conversely, approximately $40\%$ of the binaries that undergo only stable mass transfer without CE evolution contain an anti-aligned BH. Finally, accounting for uncertainties in both the population modeling and the NS equation of state, we find that $0-18.6\%$ of NSBH mergers may be accompanied by an electromagnetic counterpart. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 18 pages, 10 figures, submitted to A&A, comments are welcome

The formation of $30\,M_\odot$ merging black holes at solar metallicity

Authors: Simone S. Bavera, Tassos Fragos, Emmanouil Zapartas, Jeff J. Andrews, Vicky Kalogera, Christopher P. L. Berry, Matthias Kruckow, Aaron Dotter, Konstantinos Kovlakas, Devina Misra, Kyle A. Rocha, Philipp M. Srivastava, Meng Sun, Zepei Xing

Abstract: The maximum mass of black holes formed in isolated binaries is determined by stellar winds and the interactions between the binary components. We consider for the first time fully self-consistent detailed stellar structure and binary evolution calculations in population-synthesis models and a new, qualitatively different picture emerges for the formation of black-hole binaries, compared to studies… ▽ More The maximum mass of black holes formed in isolated binaries is determined by stellar winds and the interactions between the binary components. We consider for the first time fully self-consistent detailed stellar structure and binary evolution calculations in population-synthesis models and a new, qualitatively different picture emerges for the formation of black-hole binaries, compared to studies employing rapid population synthesis models. We find merging binary black holes can form with a non-negligible rate ($\sim 4\times10^{-7}\,M_\odot^{-1}$) at solar metallicity. Their progenitor stars with initial masses $\gtrsim 50\,M_\odot$ do not expand to supergiant radii, mostly avoiding significant dust-driven or luminous blue variable winds. Overall, the progenitor stars lose less mass in stellar winds, resulting in black holes as massive as $\sim 30\,M_\odot$, and, approximately half of them avoid a mass-transfer episode before forming the first-born black hole. Finally, binaries with initial periods of a few days, some of which may undergo episodes of Roche-lobe overflow mass transfer, result in mildly spinning first-born black holes, $χ_\mathrm{BH1} \lesssim 0.2$, assuming efficient angular-momentum transport. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 14 pages, 6 figures, manuscript submitted for publication

Direct Statistical Constraints on the Natal Kick velocity of a Black Hole in an X-ray Quiet Binary

Authors: Sharan Banagiri, Zoheyr Doctor, Vicky Kalogera, Chase Kimball, Jeff J. Andrews

Abstract: In recent years, a handful of ``dark" binaries have been discovered with a non-luminous compact object. Astrometry and radial velocity measurements of the bright companion allow us to measure the post-supernova orbital elements of such a binary. In this paper, we develop a statistical formalism to use such measurements to infer the pre-supernova orbital elements, and the natal kick imparted by the… ▽ More In recent years, a handful of ``dark" binaries have been discovered with a non-luminous compact object. Astrometry and radial velocity measurements of the bright companion allow us to measure the post-supernova orbital elements of such a binary. In this paper, we develop a statistical formalism to use such measurements to infer the pre-supernova orbital elements, and the natal kick imparted by the supernova (SN). We apply this formalism to the recent discovery of an X-ray quiet binary with a black hole, VFTS 243, in the Large Magellanic Cloud. Assuming an isotropic, Maxwellian distribution on natal kicks and using broad agnostic mass priors, we find that kick velocity can be constrained to $V_k < 72$ km/s and the dispersion of the kick distribution to $σ_k < 68 $ km/s at 90 % confidence. We find that a Blaauw kick cannot be ruled out and }that at least about $0.6 M_{\odot}$ was lost during the supernova with 90 % confidence. The pre-SN orbital separation is found to be robustly constrained to be around $0.3$ AU. △ Less

Submitted 20 January, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 13 pages, 5 figured

A Black Hole Kicked At Birth: MAXI J1305-704

Authors: Chase Kimball, Sam Imperato, Vicky Kalogera, Kyle A. Rocha, Zoheyr Doctor, Jeff J. Andrews, Aaron Dotter, Emmanouil Zapartas, Simone S. Bavera, Konstantinos Kovlakas, Tassos Fragos, Phillip M. Srivastava, Devina Misra, Meng Sun, Zepei Xing

Abstract: When a compact object is formed in a binary, any mass lost during core collapse will impart a kick on the binary's center of mass. Asymmetries in this mass loss or neutrino emission would impart an additional natal kick on the remnant black hole or neutron star, whether it was formed in a binary or in isolation. While it is well established that neutron stars receive natal kicks upon formation, it… ▽ More When a compact object is formed in a binary, any mass lost during core collapse will impart a kick on the binary's center of mass. Asymmetries in this mass loss or neutrino emission would impart an additional natal kick on the remnant black hole or neutron star, whether it was formed in a binary or in isolation. While it is well established that neutron stars receive natal kicks upon formation, it is unclear whether black holes do as well. Here, we consider the low-mass X-ray binary MAXI J1305-704, which has been reported to have a space velocity $\gtrsim$ 200 km/s. In addition to integrating its trajectory to infer its velocity upon formation of its black hole, we account for recent estimates of its period, black hole mass, mass ratio, and donor effective temperature from photometric and spectroscopic observations. We find that if MAXI J1305-704 formed via isolated binary evolution in the thick Galactic disk, then the supernova that formed its black hole imparted a natal kick of at least 70 km/s while ejecting less than $\simeq 1$ M$_\odot$ with 95% confidence assuming uninformative priors on mass loss and natal kick velocity. △ Less

Submitted 19 July, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: Accepted Astrophysical Journal Letters; 9 pages, 5 figures

Investigating the Lower Mass Gap with Low Mass X-ray Binary Population Synthesis

Authors: Jared C. Siegel, Ilia Kiato, Vicky Kalogera, Christopher P. L. Berry, Thomas J. Maccarone, Katelyn Breivik, Jeff J. Andrews, Simone S. Bavera, Aaron Dotter, Tassos Fragos, Konstantinos Kovlakas, Devina Misra, Kyle A. Rocha, Philipp M. Srivastava, Meng Sun, Zepei Xing, Emmanouil Zapartas

Abstract: Mass measurements from low-mass black hole X-ray binaries (LMXBs) and radio pulsars have been used to identify a gap between the most massive neutron stars (NSs) and the least massive black holes (BHs). BH mass measurements in LMXBs are typically only possible for transient systems: outburst periods enable detection via all-sky X-ray monitors, while quiescent periods enable radial-velocity measure… ▽ More Mass measurements from low-mass black hole X-ray binaries (LMXBs) and radio pulsars have been used to identify a gap between the most massive neutron stars (NSs) and the least massive black holes (BHs). BH mass measurements in LMXBs are typically only possible for transient systems: outburst periods enable detection via all-sky X-ray monitors, while quiescent periods enable radial-velocity measurements of the low-mass donor. We quantitatively study selection biases due to the requirement of transient behavior for BH mass measurements. Using rapid population synthesis simulations (COSMIC), detailed binary stellar-evolution models (MESA), and the disk instability model of transient behavior, we demonstrate that transient-LMXB selection effects introduce observational biases, and can suppress mass-gap BHs in the observed sample. However, we find a population of transient LMXBs with mass-gap BHs form through accretion-induced collapse of a NS during the LMXB phase, which is inconsistent with observations. These results are robust against variations of binary evolution prescriptions. The significance of this accretion-induced collapse population depends upon the maximum NS birth mass $M_\mathrm{ NS, birth-max}$. To reflect the observed dearth of low-mass BHs, COSMIC and MESA models favor $M_\mathrm{ NS, birth-max} \lesssim2M_{\odot}$. In the absence of further observational biases against LMXBs with mass-gap BHs, our results indicate the need for additional physics connected to the modeling of LMXB formation and evolution. △ Less

Submitted 25 July, 2023; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: 21 pages, accepted to ApJ

X-ray luminosity function of high-mass X-ray binaries: Studying the signatures of different physical processes using detailed binary evolution calculations

Authors: Devina Misra, Konstantinos Kovlakas, Tassos Fragos, Margaret Lazzarini, Simone S. Bavera, Bret D. Lehmer, Andreas Zezas, Emmanouil Zapartas, Zepei Xing, Jeff J. Andrews, Aaron Dotter, Kyle A. Rocha, Philipp M. Srivastava, Meng Sun

Abstract: The ever-expanding observational sample of X-ray binaries (XRBs) makes them excellent laboratories for constraining binary evolution theory. Such constraints can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing to observed XLFs. In this work, we focus on high-mass XRBs (HMXBs) and study the effects on the XLF of variou… ▽ More The ever-expanding observational sample of X-ray binaries (XRBs) makes them excellent laboratories for constraining binary evolution theory. Such constraints can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing to observed XLFs. In this work, we focus on high-mass XRBs (HMXBs) and study the effects on the XLF of various, poorly-constrained assumptions regarding physical processes such as the common-envelope phase, the core-collapse, and wind-fed accretion. We use the new binary population synthesis code POSYDON, which employs extensive pre-computed grids of detailed stellar structure and binary evolution models, to simulate the evolution of binaries. We generate 96 synthetic XRB populations corresponding to different combinations of model assumptions. The generated HMXB XLFs are feature-rich, deviating from the commonly assumed single-power law. We find a break in our synthetic XLF at luminosity $\sim 10^{38}$ erg s$^{-1}$, similar to observed XLFs. However, we find also a general overabundance of XRBs (up to a factor of $\sim$10 for certain model parameter combinations) driven primarily by XRBs with black hole accretors. Assumptions about the transient behavior of Be-XRBs, asymmetric supernova kicks, and common-envelope physics can significantly affect the shape and normalization of our synthetic XLFs. We find that less well-studied assumptions regarding the circularization of the orbit at the onset of Roche-lobe overflow and criteria for the formation of an X-ray emitting accretion disk around wind-accreting black holes can also impact our synthetic XLFs. Our study reveals the importance of large-scale parameter studies, highlighting the power of XRBs in constraining binary evolution theory. △ Less

Submitted 14 March, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: 31 pages, 32 figures, Accepted by A&A. Fixed typos and updated references. Referee's comments were addressed

Journal ref: A&A 672, A99 (2023)

A Sample of Neutron Star and Black Hole Binaries Detected through Gaia DR3 Astrometry

Authors: Jeff J. Andrews, Kirsty Taggart, Ryan Foley

Abstract: With its exquisite astrometric precision, the latest Gaia data release includes $\sim$$10^5$ astrometric binaries, each of which have measured orbital periods, eccentricities, and the Thiele-Innes orbital parameters. Using these and an estimate of the luminous stars' masses, we derive the companion stars' masses, from which we identify a sample of 24 binaries in long period orbits (… ▽ More With its exquisite astrometric precision, the latest Gaia data release includes $\sim$$10^5$ astrometric binaries, each of which have measured orbital periods, eccentricities, and the Thiele-Innes orbital parameters. Using these and an estimate of the luminous stars' masses, we derive the companion stars' masses, from which we identify a sample of 24 binaries in long period orbits ($P_{\rm orb}\sim{\rm yrs}$) with a high probability of hosting a massive ($>$1.4 $M_{\odot}$), dark companion: a neutron star (NS) or black hole (BH). The luminous stars in these binaries tend to be F-, G-, and K-dwarfs with the notable exception of one hot subdwarf. Follow-up spectroscopy of eight of these stars shows no evidence for contamination by white dwarfs or other luminous stars. The dark companions in these binaries span a mass range of 1.35-2.7 $M_{\odot}$ and therefore likely includes both NSs and BHs without a significant mass gap in between. Furthermore, the masses of several of these objects are $\simeq$1.7 $M_{\odot}$, similar to the mass of at least one of the merging compact objects in GW190425. Given that these orbits are too wide for significant mass accretion to have occurred, this sample implies that some NSs are born heavy ($\gtrsim$1.5 $M_{\odot}$). Additionally, the low orbital velocities ($\lesssim$20 km s$^{-1}$) of these binaries requires that at least some heavy NSs receive low natal kicks, otherwise they would have been disrupted during core collapse. Although none will become gravitational wave sources within a Hubble time, these systems will be exceptionally useful for testing binary evolution theory. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Comments: 12 pages, 6 figures, 2 tables, submitted AAS Journals, comments welcome

Weighing the Darkness III: How Gaia Could, but Probably Won't, Astrometrically Detect Free-Floating Black Holes

Authors: Jeff J. Andrews

Abstract: The gravitational pull of an unseen companion to a luminous star is well-known to cause deviations to the parallax and proper motion of a star. In a previous paper in this series, we argue that the astrometric mission Gaia can identify long-period binaries by precisely measuring these arcs. An arc in a star's path can also be caused by a fly-by -- the hyperbolic encounter with another massive obje… ▽ More The gravitational pull of an unseen companion to a luminous star is well-known to cause deviations to the parallax and proper motion of a star. In a previous paper in this series, we argue that the astrometric mission Gaia can identify long-period binaries by precisely measuring these arcs. An arc in a star's path can also be caused by a fly-by -- the hyperbolic encounter with another massive object. We quantify the apparent acceleration over time induced by a companion star as a function of the impact parameter, velocity of interaction, and companion mass. In principle, Gaia could be used to astrometrically identify the contribution of massive compact halo objects to the local dark matter potential of the Milky Way. However, after quantifying their rate and Gaia's sensitivity, we find that fly-bys are so rare that Gaia will probably never observe one. Therefore every star in the Gaia database exhibiting astrometric acceleration is likely in a long-period binary with another object. Nevertheless, we show how intermediate mass black holes, if they exist in the Solar Neighborhood, could be detected by the anomalously large accelerations they induce on nearby stars. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 8 pages, 5 figures, submitted to ApJ. Comments welcome

Active Learning for Computationally Efficient Distribution of Binary Evolution Simulations

Authors: Kyle Akira Rocha, Jeff J. Andrews, Christopher P. L. Berry, Zoheyr Doctor, Aggelos K. Katsaggelos, Juan Gabriel Serra Pérez, Pablo Marchant, Vicky Kalogera, Scott Coughlin, Simone S. Bavera, Aaron Dotter, Tassos Fragos, Konstantinos Kovlakas, Devina Misra, Zepei Xing, Emmanouil Zapartas

Abstract: Binary stars undergo a variety of interactions and evolutionary phases, critical for predicting and explaining observed properties. Binary population synthesis with full stellar-structure and evolution simulations are computationally expensive requiring a large number of mass-transfer sequences. The recently developed binary population synthesis code POSYDON incorporates grids of MESA binary star… ▽ More Binary stars undergo a variety of interactions and evolutionary phases, critical for predicting and explaining observed properties. Binary population synthesis with full stellar-structure and evolution simulations are computationally expensive requiring a large number of mass-transfer sequences. The recently developed binary population synthesis code POSYDON incorporates grids of MESA binary star simulations which are then interpolated to model large-scale populations of massive binaries. The traditional method of computing a high-density rectilinear grid of simulations is not scalable for higher-dimension grids, accounting for a range of metallicities, rotation, and eccentricity. We present a new active learning algorithm, psy-cris, which uses machine learning in the data-gathering process to adaptively and iteratively select targeted simulations to run, resulting in a custom, high-performance training set. We test psy-cris on a toy problem and find the resulting training sets require fewer simulations for accurate classification and regression than either regular or randomly sampled grids. We further apply psy-cris to the target problem of building a dynamic grid of MESA simulations, and we demonstrate that, even without fine tuning, a simulation set of only $\sim 1/4$ the size of a rectilinear grid is sufficient to achieve the same classification accuracy. We anticipate further gains when algorithmic parameters are optimized for the targeted application. We find that optimizing for classification only may lead to performance losses in regression, and vice versa. Lowering the computational cost of producing grids will enable future versions of POSYDON to cover more input parameters while preserving interpolation accuracies. △ Less

Submitted 16 September, 2022; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: 21 pages, 10 figures, ApJ in press

Journal ref: Astrophysical Journal; 938(1):64(15); 2022

Constraining Black Hole Natal Kicks with Astrometric Microlensing

Authors: Jeff J. Andrews, Vicky Kalogera

Abstract: Multiple pieces of evidence suggest that neutron stars receive large kicks when formed from the remnant of a collapsing star. However, the evidence for whether black holes (BH) receive natal kicks is less clear, reliant on weak constraints from the analysis of BH X-ray binaries and massive runaway and walkaway stars. Here we show for the first time that recent microlensing detections offer a new m… ▽ More Multiple pieces of evidence suggest that neutron stars receive large kicks when formed from the remnant of a collapsing star. However, the evidence for whether black holes (BH) receive natal kicks is less clear, reliant on weak constraints from the analysis of BH X-ray binaries and massive runaway and walkaway stars. Here we show for the first time that recent microlensing detections offer a new method for measuring the kicks BHs receive at birth. When a BH is identified through both photometric and astrometric microlensing and when the lensed star has a known distance and proper motion, the mass, distance and proper motion of the BH can be determined. We study the runaway velocities for components of eccentric binaries disrupted during a supernova, finding the peculiar velocity correlates strongly with the kick a BH received a birth, typically within 20\%, even when the natal kick is smaller than the orbital velocity. Therefore, by measuring the peculiar velocity of a BH or other compact object that formed from a binary that disrupted during core collapse, we are in effect measuring the natal kick that object received. We focus on MOA-2011-BLG-191/OGLE-2011-BLG-0462, an isolated, single BH detected by microlensing and consider a range of possible formation scenarios, including its formation from the disruption of a binary during a supernova event. We determine that MOA-2011-BLG-191/OGLE-2011-BLG-0462 has a Milky Way orbit consistent with a thick disk population, but if it were formed within the kinematic thin disk it received a natal kick $\lesssim$100 km s$^{-1}$. △ Less

Submitted 15 April, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 14 pages, 8 figures, accepted for publication in ApJ

POSYDON: A General-Purpose Population Synthesis Code with Detailed Binary-Evolution Simulations

Authors: Tassos Fragos, Jeff J. Andrews, Simone S. Bavera, Christopher P. L. Berry, Scott Coughlin, Aaron Dotter, Prabin Giri, Vicky Kalogera, Aggelos Katsaggelos, Konstantinos Kovlakas, Shamal Lalvani, Devina Misra, Philipp M. Srivastava, Ying Qin, Kyle A. Rocha, Jaime Roman-Garza, Juan Gabriel Serra, Petter Stahle, Meng Sun, Xu Teng, Goce Trajcevski, Nam Hai Tran, Zepei Xing, Emmanouil Zapartas, Michael Zevin

Abstract: Most massive stars are members of a binary or a higher-order stellar systems, where the presence of a binary companion can decisively alter their evolution via binary interactions. Interacting binaries are also important astrophysical laboratories for the study of compact objects. Binary population synthesis studies have been used extensively over the last two decades to interpret observations of… ▽ More Most massive stars are members of a binary or a higher-order stellar systems, where the presence of a binary companion can decisively alter their evolution via binary interactions. Interacting binaries are also important astrophysical laboratories for the study of compact objects. Binary population synthesis studies have been used extensively over the last two decades to interpret observations of compact-object binaries and to decipher the physical processes that lead to their formation. Here, we present POSYDON, a novel, binary population synthesis code that incorporates full stellar-structure and binary-evolution modeling, using the MESA code, throughout the whole evolution of the binaries. The use of POSYDON enables the self-consistent treatment of physical processes in stellar and binary evolution, including: realistic mass-transfer calculations and assessment of stability, internal angular-momentum transport and tides, stellar core sizes, mass-transfer rates and orbital periods. This paper describes the detailed methodology and implementation of POSYDON, including the assumed physics of stellar- and binary-evolution, the extensive grids of detailed single- and binary-star models, the post-processing, classification and interpolation methods we developed for use with the grids, and the treatment of evolutionary phases that are not based on pre-calculated grids. The first version of POSYDON targets binaries with massive primary stars (potential progenitors of neutron stars or black holes) at solar metallicity. △ Less

Submitted 7 August, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: 60 pages, 33 figures, 8 tables, referee's comments addressed. The code and the accompanying documentations and data products are available at https:\\posydon.org

A Young, Low-Density Stellar Stream in the Milky Way Disk: Theia 456

Authors: Jeff J. Andrews, Jason L. Curtis, Julio Chanamé, Marcel A. Agüeros, Simon C. Schuler, Marina Kounkel, Kevin R. Covey

Abstract: Our view of the variety of stellar structures pervading the local Milky Way has been transformed by the application of clustering algorithms to the Gaia catalog. In particular, several stellar streams have been recently discovered that are comprised of hundreds to thousands of stars and span several hundred parsecs. We analyze one such structure, Theia 456, a low-density stellar stream extending n… ▽ More Our view of the variety of stellar structures pervading the local Milky Way has been transformed by the application of clustering algorithms to the Gaia catalog. In particular, several stellar streams have been recently discovered that are comprised of hundreds to thousands of stars and span several hundred parsecs. We analyze one such structure, Theia 456, a low-density stellar stream extending nearly 200 pc and 20$^{\circ}$ across the sky. By supplementing Gaia astrometric data with spectroscopic metallicities from LAMOST and photometric rotation periods from the Zwicky Transient Facility (ZTF) and the Transiting Exoplanet Survey Satellite (TESS), we establish Theia 456's radial velocity coherence, and we find strong evidence that members of Theia 456 have a common age ($\simeq$175 Myr), common dynamical origin, and formed from chemically homogeneous pre-stellar material ([Fe/H] = $-$0.07 dex). Unlike well-known stellar streams in the Milky Way, which are in its halo, Theia 456 is firmly part of the thin disk. If our conclusions about Theia 456 can be applied to even a small fraction of the remaining $\simeq$8300 independent structures in the Theia catalog, such low-density stellar streams may be ubiquitous. We comment on the implications this has for the nature of star-formation throughout the Galaxy. △ Less

Submitted 15 April, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 12 pages, 10 figures, 1 table; accepted for publication in AJ

Gaia May Detect Hundreds of Well-characterised Stellar Black Holes

Authors: Chirag Chawla, Sourav Chatterjee, Katelyn Breivik, Chaithanya Krishna Moorthy, Jeff J. Andrews, Robyn E. Sanderson

Abstract: Detection of black holes (BHs) with detached luminous companions (LCs) can be instrumental in connecting the BH properties with their progenitors' since the latter can be inferred from the observable properties of the LC. Past studies showed the promise of Gaia astrometry in detecting BH-LC binaries. We build upon these studies by: 1) initialising the zero-age binary properties based on realistic,… ▽ More Detection of black holes (BHs) with detached luminous companions (LCs) can be instrumental in connecting the BH properties with their progenitors' since the latter can be inferred from the observable properties of the LC. Past studies showed the promise of Gaia astrometry in detecting BH-LC binaries. We build upon these studies by: 1) initialising the zero-age binary properties based on realistic, metallicity-dependent star-formation history in the Milky Way (MW), 2) evolving these binaries to current epoch to generate realistic MW populations of BH-LC binaries, 3) distributing these binaries in the MW preserving the complex age-metallicity-Galactic position correlations, 4) accounting for extinction and reddening using three-dimensional dust maps, 5) examining the extended Gaia mission's ability to resolve BH-LC binaries. We restrict ourselves to detached BH-LC binaries with orbital period <10 yr such that Gaia can observe at least one full orbit. We find: 1) the extended Gaia mission can astrometrically resolve 30-300 detached BH-LC binaries depending on our assumptions of supernova physics and astrometric detection threshold; 2) Gaia's astrometry alone can indicate BH candidates for 10-100 BH-LC binaries by constraining the dark primary mass >3 Msun; 3) distributions of observables including orbital periods, eccentricities, and component masses are sensitive to the adopted binary evolution model, hence can directly inform binary evolution models. Finally, we comment on the potential to further characterise these BH binaries through radial velocity measurements and observation of X-ray counterparts. △ Less

Submitted 17 March, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 26 pages, 18 figures, and 1 table; accepted for publication in The Astrophysical Journal

Weighing the Darkness II: Astrometric Measurement of Partial Orbits with Gaia

Authors: Jeff J. Andrews, Katelyn Breivik, Chirag Chawla, Carl Rodriguez, Sourav Chatterjee

Abstract: Over the course of several years, stars trace helical trajectories as they traverse across the sky due to the combined effects of proper motion and parallax. It is well known that the gravitational pull of an unseen companion can cause deviations to these tracks. Several studies have pointed out that the astrometric mission Gaia will be able to identify a slew of new exoplanets, stellar binaries,… ▽ More Over the course of several years, stars trace helical trajectories as they traverse across the sky due to the combined effects of proper motion and parallax. It is well known that the gravitational pull of an unseen companion can cause deviations to these tracks. Several studies have pointed out that the astrometric mission Gaia will be able to identify a slew of new exoplanets, stellar binaries, and compact object companions with orbital periods as short as tens of days to as long as Gaia's lifetime. Here, we use mock astrometric observations to demonstrate that Gaia can identify and characterize black hole companions to luminous stars with orbital periods longer than Gaia's lifetime. Such astrometric binaries have orbital periods too long to exhibit complete orbits, and instead are identified through curvature in their characteristic helical paths. By simultaneously measuring the radius of this curvature and the orbital velocity, constraints can be placed on the underlying orbit. We quantify the precision with which Gaia can measure orbital accelerations and apply that to model predictions for the population of black holes orbiting stars in the stellar neighborhood. Although orbital degeneracies imply that many of the accelerations induced by hidden black holes could also be explained by faint low-mass stars, we discuss how the nature of certain putative black hole companions can be confirmed with high confidence using Gaia data alone. △ Less

Submitted 11 October, 2021; originally announced October 2021.

Comments: 12 pages, 11 figures; submitted to AAS Journals; comments welcome

Probing the progenitors of spinning binary black-hole mergers with long gamma-ray bursts

Authors: Simone S. Bavera, Tassos Fragos, Emmanouil Zapartas, Enrico Ramirez-Ruiz, Pablo Marchant, Luke Z. Kelley, Michael Zevin, Jeff J. Andrews, Scott Coughlin, Aaron Dotter, Konstantinos Kovlakas, Devina Misra, Juan G. Serra-Perez, Ying Qin, Kyle A. Rocha, Jaime Román-Garza, Nam H. Tran, Zepei Xing

Abstract: Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and popul… ▽ More Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed sub-population of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ~10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ~85% and ~60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction $f_\mathrm{B}\sim 0.05$ our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift $z \in [0,2.5]$, a tentative comparison with the metallicity distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances. △ Less

Submitted 3 December, 2021; v1 submitted 30 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A Letters, 12 pages, 6 figures, 1 table

Combining Astrometry and Elemental Abundances: The Case of the Candidate Pre-Gaia Halo Moving Groups G03-37, G18-39, and G21-22

Authors: S. C. Schuler, J. J. Andrews, V. R. Clanzy II, M. Mourabit, J. Chanamé, M. A. Agüeros

Abstract: While most moving groups are young and nearby, a small number have been identified in the Galactic halo. Understanding the origin and evolution of these groups is an important piece of reconstructing the formation history of the halo. Here we report on our analysis of three putative halo moving groups: G03-37, G18-39, and G21-22. Based on Gaia EDR3 data, the stars associated with each group show s… ▽ More While most moving groups are young and nearby, a small number have been identified in the Galactic halo. Understanding the origin and evolution of these groups is an important piece of reconstructing the formation history of the halo. Here we report on our analysis of three putative halo moving groups: G03-37, G18-39, and G21-22. Based on Gaia EDR3 data, the stars associated with each group show some scatter in velocity (e.g., Toomre diagram) and integrals of motion (energy, angular momentum) spaces, counter to expectations of moving-group stars. We choose the best candidate of the three groups, G21-22, for follow-up chemical analysis based on high-resolution spectroscopy of six presumptive members. Using a new Python code that uses a Bayesian method to self-consistently propagate uncertainties from stellar atmosphere solutions in calculating individual abundances and spectral synthesis, we derive the abundances of $α$- (Mg, Si, Ca, Ti), Fe-peak (Cr, Sc, Mn, Fe, Ni), odd-$Z$ (Na, Al, V), and neutron-capture (Ba, Eu) elements for each star. We find that the G21-22 stars are not chemically homogeneous. Based on the kinematic analysis for all three groups and the chemical analysis for G21-22, we conclude the three are not genuine moving groups. The case for G21-22 demonstrates the benefit of combining kinematic and chemical information in identifying conatal populations when either alone may be insufficient. Comparing the integrals of motion and velocities of the six G21-22 stars with those of known structures in the halo, we tentatively associate them with the Gaia-Enceladus accretion event. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: 21 pages, 12 figures, seven tables; accepted to AJ

The Heraklion Extragalactic Catalogue (HECATE): a value-added galaxy catalogue for multi-messenger astrophysics

Authors: Konstantinos Kovlakas, Andreas Zezas, Jeff J. Andrews, Antara Basu-Zych, Tassos Fragos, Ann Hornschemeier, Konstantinos Kouroumpatzakis, Bret Lehmer, Andrew Ptak

Abstract: We present the Heraklion Extragalactic Catalogue, or HECATE, an all-sky value-added galaxy catalogue, aiming to facilitate present and future multi-wavelength and multi-messenger studies in the local Universe. It contains 204,733 galaxies up to a redshift of 0.047 (D<200 Mpc), and it is >50% complete in terms of the B-band luminosity density at distances in the 0-170 Mpc range. By incorporating an… ▽ More We present the Heraklion Extragalactic Catalogue, or HECATE, an all-sky value-added galaxy catalogue, aiming to facilitate present and future multi-wavelength and multi-messenger studies in the local Universe. It contains 204,733 galaxies up to a redshift of 0.047 (D<200 Mpc), and it is >50% complete in terms of the B-band luminosity density at distances in the 0-170 Mpc range. By incorporating and homogenising data from astronomical databases and multi-wavelength surveys, the catalogue offers positions, sizes, distances, morphological classifications, star-formation rates, stellar masses, metallicities, and nuclear activity classifications. This wealth of information can enable a wide-range of applications, such as: (i) demographic studies of extragalactic sources, (ii) initial characterisation of transient events, and (iii) searches for electromagnetic counterparts of gravitational-wave events. The catalogue is publicly available to the community at a dedicated portal, which will also host future extensions in terms of the covered volume, and data products. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 21 pages. 11 figures, 5 tables. This article has been accepted for publication in the Monthly Notices of the Royal Astronomical Society, published by Oxford University Press on behalf of the Royal Astronomical Society

Revisiting the explodability of single massive star progenitors of stripped-envelope supernovae

Authors: E. Zapartas, M. Renzo, T. Fragos, A. Dotter, J. J. Andrews, S. S. Bavera, S. Coughlin, D. Misra, K. Kovlakas, J. Román-Garza, J. G. Serra, Y. Qin, K. A. Rocha, N. H. Tran, Z. P. Xing

Abstract: Stripped-envelope supernovae (Types IIb, Ib, and Ic) that show little or no hydrogen comprise roughly one-third of the observed explosions of massive stars. Their origin and the evolution of their progenitors are not yet fully understood. Very massive single stars stripped by their own winds ($\gtrsim 25-30 M_{\odot}$ at solar metallicity) are considered viable progenitors of these events. However… ▽ More Stripped-envelope supernovae (Types IIb, Ib, and Ic) that show little or no hydrogen comprise roughly one-third of the observed explosions of massive stars. Their origin and the evolution of their progenitors are not yet fully understood. Very massive single stars stripped by their own winds ($\gtrsim 25-30 M_{\odot}$ at solar metallicity) are considered viable progenitors of these events. However, recent 1D core-collapse simulations show that some massive stars may collapse directly into black holes after a failed explosion, with a weak or no visible transient. In this letter, we estimate the effect of direct collapse into a black hole on the rates of stripped-envelope supernovae that arise from single stars. For this, we compute single-star MESA models at solar metallicity and map their final state to their core-collapse outcome following prescriptions commonly used in population synthesis. According to our models, no single stars that have lost their entire hydrogen-rich envelope are able to explode, and only a fraction of progenitors left with a thin hydrogen envelope do (IIb progenitor candidates), unless we use a prescription that takes the effect of turbulence into account or invoke increased wind mass-loss rates. This result increases the existing tension between the single-star paradigm to explain most stripped-envelope supernovae and their observed rates and properties. At face value, our results point toward an even higher contribution of binary progenitors to stripped-envelope supernovae. Alternatively, they may suggest inconsistencies in the common practice of mapping different stellar models to core-collapse outcomes and/or higher overall mass loss in massive stars. △ Less

Submitted 17 December, 2021; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: Published in Astronomy & Astrophysics Letters; One main enhancement: added Couch et al. (2020) in the list of supernova engines

Journal ref: A&A 656, L19 (2021)

The hydrodynamic evolution of binary black holes embedded within the vertically stratified disks of active galactic nuclei

Authors: Nicholas Kaaz, Sophie Lund Schrøder, Jeff J. Andrews, Andrea Antoni, Enrico Ramirez-Ruiz

Abstract: Stellar-mass black holes can become embedded within the gaseous disks of active galactic nuclei (AGNs). Afterwards, their interactions are mediated by their gaseous surroundings. In this work, we study the evolution of stellar-mass binary black holes (BBHs) embedded within AGN disks using a combination of three-dimensional hydrodynamic simulations and analytic methods, focusing on environments in… ▽ More Stellar-mass black holes can become embedded within the gaseous disks of active galactic nuclei (AGNs). Afterwards, their interactions are mediated by their gaseous surroundings. In this work, we study the evolution of stellar-mass binary black holes (BBHs) embedded within AGN disks using a combination of three-dimensional hydrodynamic simulations and analytic methods, focusing on environments in which the AGN disk scale height $H$ is $\gtrsim$ the BBH sphere of influence. We model the local surroundings of the embedded BBHs using a wind tunnel formalism and characterize different accretion regimes based on the local properties of the disk, which range from wind-dominated to quasi-spherical. We use our simulations to develop prescriptions for mass accretion and drag for embedded BBHs. We use these prescriptions, along with AGN disk models that can represent the Toomre-unstable outer regions of AGN disks, to study the long-term evolution of the BBHs as they migrate through the disk. We find that BBHs typically merge within $\lesssim 5-30\,{\rm Myr}$, increasing their mass significantly in the process, allowing BBHs to enter (or cross) the pair-instability supernova mass gap. The rate at which gas is supplied to these BBHs often exceeds the Eddington limit, sometimes by several orders of magnitude. We conclude that most embedded BBHs will merge before migrating significantly in the disk. Depending on the conditions of the ambient gas and the distance to the system, LISA can detect the transition between the gas-dominated and gravitational wave dominated regime for inspiraling BBHs that are formed sufficiently close to the AGN ($\lesssim$ 0.1 pc). We also discuss possible electromagnetic signatures during and following the inspiral, finding that it is generally unlikely but not inconceivable for the bolometric luminosity of the BBH to exceed that of the host AGN. △ Less

Submitted 25 January, 2023; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: 19 pages, 8 figures, accepted to ApJ

Targeted modeling of GW150914's binary black hole source with dart_board

Authors: Jeff J. Andrews, Julianne Cronin, Vicky Kalogera, Christopher Berry, Andreas Zezas

Abstract: We present a new method to extract statistical constraints on the progenitor properties and formation channels of individual gravitational-wave sources. Although many different models have been proposed to explain the binary black holes detected by the LIGO Scientific and Virgo Collaboration (LVC), formation through isolated binary evolution remains the best explored channel. Under the assumption… ▽ More We present a new method to extract statistical constraints on the progenitor properties and formation channels of individual gravitational-wave sources. Although many different models have been proposed to explain the binary black holes detected by the LIGO Scientific and Virgo Collaboration (LVC), formation through isolated binary evolution remains the best explored channel. Under the assumption of formation through binary evolution, we use the statistical wrapper dart_board coupled with the rapid binary evolution code COSMIC to model the progenitor of GW150914, the first gravitational-wave signal detected by the LVC. Our Bayesian method combines the likelihood generated from the gravitational-wave signal with a prior describing the population of stellar binaries, and the Universe's star-formation and metallicity evolution. We find that the dominant evolutionary channel for GW150914 did not involve a common-envelope phase, but instead the system most probably (70%-90%) formed through stable mass transfer. This result is robust against variations of various model parameters, and it is reversed only when dynamical instability in binaries becomes more likely when a strict condition favoring common envelopes is adopted. Our analysis additionally provides a quantitative description of the progenitors relevant to each channel. △ Less

Submitted 16 May, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 14 pages, 5 figures, 1 table, accepted for publication in ApJL

The impact of mass-transfer physics on the observable properties of field binary black hole populations

Authors: Simone S. Bavera, Tassos Fragos, Michael Zevin, Christopher P. L. Berry, Pablo Marchant, Jeff J. Andrews, Scott Coughlin, Aaron Dotter, Konstantinos Kovlakas, Devina Misra, Juan G. Serra-Perez, Ying Qin, Kyle A. Rocha, Jaime Román-Garza, Nam H. Tran, Emmanouil Zapartas

Abstract: We study the impact of mass-transfer physics on the observable properties of binary black hole populations formed through isolated binary evolution. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of $χ_{eff}$, $M_{chirp}$ and $q$. We find that low common envelope efficiency translates to tighter orbits post c… ▽ More We study the impact of mass-transfer physics on the observable properties of binary black hole populations formed through isolated binary evolution. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of $χ_{eff}$, $M_{chirp}$ and $q$. We find that low common envelope efficiency translates to tighter orbits post common envelope and therefore more tidally spun up second-born black holes. However, these systems have short merger timescales and are only marginally detectable by current gravitational-waves detectors as they form and merge at high redshifts ($z\sim 2$), outside current detector horizons. Assuming Eddington-limited accretion efficiency and that the first-born black hole is formed with a negligible spin, we find that all non-zero $χ_{eff}$ systems in the detectable population can come only from the common envelope channel as the stable mass-transfer channel cannot shrink the orbits enough for efficient tidal spin-up to take place. We find the local rate density ($z\simeq 0.01$) for the common envelope channel is in the range $\sim 17-113~Gpc^{-3}yr^{-1}$ considering a range of $α_{CE} \in [0.2,5.0]$ while for the stable mass transfer channel the rate density is $\sim 25~Gpc^{-3}yr^{-1}$. The latter drops by two orders of magnitude if the mass accretion onto the black hole is not Eddington limited because conservative mass transfer does not shrink the orbit as efficiently as non-conservative mass transfer does. Finally, using GWTC-2 events, we constrain the lower bound of branching fraction from other formation channels in the detected population to be $\sim 0.2$. Assuming all remaining events to be formed through either stable mass transfer or common envelope channels, we find moderate to strong evidence in favour of models with inefficient common envelopes. △ Less

Submitted 15 February, 2021; v1 submitted 30 October, 2020; originally announced October 2020.

Comments: 26 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 647, A153 (2021)

A census of ultraluminous X-ray sources in the local Universe

Authors: Konstantinos Kovlakas, Andreas Zezas, Jeff J. Andrews, Antara Basu-Zych, Tassos Fragos, Ann Hornschemeier, Bret Lehmer, Andrew Ptak

Abstract: Using the Chandra Source Catalog 2.0 and a newly compiled catalogue of galaxies in the local Universe, we deliver a census of ultraluminous X-ray source (ULX) populations in nearby galaxies. We find 629 ULX candidates in 309 galaxies with distance smaller than 40\,Mpc. The foreground/background contamination is ${\sim}20\%$. The ULX populations in bona-fide star-forming galaxies scale on average w… ▽ More Using the Chandra Source Catalog 2.0 and a newly compiled catalogue of galaxies in the local Universe, we deliver a census of ultraluminous X-ray source (ULX) populations in nearby galaxies. We find 629 ULX candidates in 309 galaxies with distance smaller than 40\,Mpc. The foreground/background contamination is ${\sim}20\%$. The ULX populations in bona-fide star-forming galaxies scale on average with star-formation rate (SFR) and stellar mass ($M_\star$) such that the number of ULXs per galaxy is $0.45^{+0.06}_{-0.09}\times\frac{\rm SFR}{\rm M_\odot\,yr^{-1}}{+}3.3^{+3.8}_{-3.2}\times\frac{M_\star}{\rm M_\odot}$. The scaling depends strongly on the morphological type. This analysis shows that early spiral galaxies contain an additional population of ULXs that scales with $M_\star$. We also confirm the strong anti-correlation of the ULX rate with the host galaxy's metallicity. In the case of early-type galaxies we find that there is a non-linear dependence of the number of ULXs with $M_\star$, which is interpreted as the result of star-formation history differences. Taking into account age and metallicity effects, we find that the predictions from X-ray binary population synthesis models are consistent with the observed ULX rates in early-type galaxies, as well as, spiral/irregular galaxies. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: 21 pages, 19 figures, 6 tables. Accepted for publication in MNRAS

Mass Ratios of Merging Double Neutron Stars as Implied by the Milky Way Population

Authors: Jeff J. Andrews

Abstract: Of the seven known double neutron stars (DNS) with precisely measure masses in the Milky Way that will merge within a Hubble time, all but one has a mass ratio, $q$, close to unity. Recently, precise measurements of three post-Keplerian parameters in the DNS J1913$+$1102 constrain this system to have a significantly non-unity mass ratio of 0.78$\pm$0.03. One may be tempted to conclude that approxi… ▽ More Of the seven known double neutron stars (DNS) with precisely measure masses in the Milky Way that will merge within a Hubble time, all but one has a mass ratio, $q$, close to unity. Recently, precise measurements of three post-Keplerian parameters in the DNS J1913$+$1102 constrain this system to have a significantly non-unity mass ratio of 0.78$\pm$0.03. One may be tempted to conclude that approximately one out of seven (14\%) DNS mergers detected by gravitational wave observatories will have mass ratios significantly different from unity. However J1913$+$1102 has a relatively long merger time of 470 Myr. We show that when merger times and observational biases are taken into account, the population of Galactic DNSs imply that $\simeq98\%$ of all merging DNSs will have $q>$0.9. We then apply two separate fitting formulas informed by 3D hydrodynamic simulations of DNS mergers to our results on Galactic DNS masses, finding that either $\simeq$0.004 ${M_{\odot}}$ or $\simeq$0.010 ${M_{\odot}}$ of material will be ejected at merger, depending on which formula is used. These ejecta masses have implications for both the peak bolometric luminosities of electromagnetic counterparts (which we find to be $\sim$10$^{41}$ erg s$^{-1}$) as well as the $r$-process enrichment of the Milky Way. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Comments: 7 pages, 3 figures, 1 table, submitted to ApJL. Comments welcome

The ELM Survey. VIII. 98 Double White Dwarf Binaries

Authors: Warren R. Brown, Mukremin Kilic, Alekzander Kosakowski, Jeff J. Andrews, Craig O. Heinke, Marcel A. Agueros, Fernando Camilo, A. Gianninas, J. J. Hermes, Scott J. Kenyon

Abstract: We present the final sample of 98 detached double white dwarf (WD) binaries found in the Extremely Low Mass (ELM) Survey, a spectroscopic survey targeting <0.3 Msun He-core WDs completed in the Sloan Digital Sky Survey footprint. Over the course of the survey we observed ancillary low mass WD candidates like GD278, which we show is a P=0.19 d double WD binary, as well as candidates that turn out t… ▽ More We present the final sample of 98 detached double white dwarf (WD) binaries found in the Extremely Low Mass (ELM) Survey, a spectroscopic survey targeting <0.3 Msun He-core WDs completed in the Sloan Digital Sky Survey footprint. Over the course of the survey we observed ancillary low mass WD candidates like GD278, which we show is a P=0.19 d double WD binary, as well as candidates that turn out to be field blue straggler/subdwarf A-type stars with luminosities too large to be WDs given their Gaia parallaxes. Here, we define a clean sample of ELM WDs that is complete within our target selection and magnitude range 15<g_0<20 mag. The measurements are consistent with 100% of ELM WDs being 0.0089 < P < 1.5 d double WD binaries, 35% of which belong to the Galactic halo. We infer these are mostly He+CO WD binaries given the measurement constraints. The merger rate of the observed He+CO WD binaries exceeds the formation rate of stable mass transfer AM CVn binaries by a factor of 25, and so the majority of He+CO WD binaries must experience unstable mass transfer and merge. The shortest-period systems like J0651+2844 are signature {\it LISA} verification binaries that can be studied with gravitational waves and light. △ Less

Submitted 31 January, 2020; originally announced February 2020.

Comments: 22 pages, ApJ accepted

COSMIC Variance in Binary Population Synthesis

Authors: Katelyn Breivik, Scott Coughlin, Michael Zevin, Carl L. Rodriguez, Kyle Kremer, Claire S. Ye, Jeff J. Andrews, Michael Kurkowski, Matthew C. Digman, Shane L. Larson, Frederic A. Rasio

Abstract: The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observat… ▽ More The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA. △ Less

Submitted 10 September, 2020; v1 submitted 3 November, 2019; originally announced November 2019.

Comments: reflects the version submitted to ApJ; 17 pages, 2 Tables, 4 Figures; corrects typo in SNR calculation

LISA and the Existence of a Fast-Merging Double Neutron Star Formation Channel

Authors: Jeff J. Andrews, Katelyn Breivik, Chris Pankow, Daniel J. D'Orazio, Mohammadtaher Safarzadeh

Abstract: Using a Milky Way double neutron star (DNS) merger rate of 210 Myr$^{-1}$, as derived by the Laser Interferometer Gravitational-Wave Observatory (LIGO), we demonstrate that the Laser Interferometer Space Antenna (LISA) will detect on average 240 (330) DNSs within the Milky Way for a 4-year (8-year) mission with a signal-to-noise ratio greater than 7. Even adopting a more pessimistic rate of 42 Myr… ▽ More Using a Milky Way double neutron star (DNS) merger rate of 210 Myr$^{-1}$, as derived by the Laser Interferometer Gravitational-Wave Observatory (LIGO), we demonstrate that the Laser Interferometer Space Antenna (LISA) will detect on average 240 (330) DNSs within the Milky Way for a 4-year (8-year) mission with a signal-to-noise ratio greater than 7. Even adopting a more pessimistic rate of 42 Myr$^{-1}$, as derived by the population of Galactic DNSs, we find a significant detection of 46 (65) Milky Way DNSs. These DNSs can be leveraged to constrain formation scenarios. In particular, traditional NS-discovery methods using radio telescopes are unable to detect DNSs with $P_{\rm orb}$ $\lesssim$1 hour (merger times $\lesssim$10 Myr). If a fast-merging channel exists that forms DNSs at these short orbital periods, LISA affords, perhaps, the only opportunity to observationally characterize these systems; we show that toy models for possible formation scenarios leave unique imprints on DNS orbital eccentricities, which may be measured by LISA for values as small as $\sim$10$^{-2}$. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: 10 pages, 4 figures, submitted to ApJL

Weighing the Darkness: Astrometric Mass Measurement of Hidden Stellar Companions using Gaia

Authors: Jeff J. Andrews, Katelyn Breivik, Sourav Chatterjee

Abstract: In astrometric binaries, the presence of a dark, unseen star can be inferred from the gravitational pull it induces on its luminous binary companion. While the orbit of such binaries can be characterized with precise astrometric measurements, constraints made from astrometry alone are not enough to measure the component masses. In this work, we determine the precision with which Gaia can astrometr… ▽ More In astrometric binaries, the presence of a dark, unseen star can be inferred from the gravitational pull it induces on its luminous binary companion. While the orbit of such binaries can be characterized with precise astrometric measurements, constraints made from astrometry alone are not enough to measure the component masses. In this work, we determine the precision with which Gaia can astrometrically measure the orbits and -- with additional observations -- the component masses, for luminous stars hosting hidden companions. Using realistic mock Gaia observations, we find that Gaia can precisely measure the orbits of binaries hosting hidden brown-dwarfs out to tens of pc and hidden white dwarf and neutron star companions at distances as far as several hundred pc. Heavier black hole companions may be measured out to 1 kpc or farther. We further determine how orbital period affects this precision, finding that Gaia can characterize orbits with periods as short as 10 days and as long as a few 10$^3$ days, with the best measured orbits having periods just short of Gaia's mission lifetime. Extending Gaia's nominal five-year mission lifetime by an additional five years not only allows for the measurement of longer period orbits, but those longer period binaries can be seen at even greater distances. △ Less

Submitted 12 September, 2019; originally announced September 2019.

Comments: 14 pages, 7 figures, 1 table, accepted for publication in ApJ

The Complete Evolution of a Neutron-Star Binary through a Common Envelope Phase Using 1D Hydrodynamic Simulations

Authors: Tassos Fragos, Jeff J. Andrews, Enrico Ramirez-Ruiz, Georges Meynet, Vicky Kalogera, Ronald E. Taam, Andreas Zezas

Abstract: Over forty years of research suggests that the common envelope phase, in which an evolved star engulfs its companion upon expansion, is the critical evolutionary stage forming short-period, compact-object binary systems, such as coalescing double compact objects, X-ray binaries, and cataclysmic variables. In this work, we adapt the one-dimensional hydrodynamic stellar evolution code, MESA, to mode… ▽ More Over forty years of research suggests that the common envelope phase, in which an evolved star engulfs its companion upon expansion, is the critical evolutionary stage forming short-period, compact-object binary systems, such as coalescing double compact objects, X-ray binaries, and cataclysmic variables. In this work, we adapt the one-dimensional hydrodynamic stellar evolution code, MESA, to model the inspiral of a 1.4M$_{\odot}$ neutron star (NS) inside the envelope of a 12$M_{\odot}$ red supergiant star. We self-consistently calculate the drag force experienced by the NS as well as the back-reaction onto the expanding envelope as the NS spirals in. Nearly all of the hydrogen envelope escapes, expanding to large radii ($\sim$10$^2$ AU) where it forms an optically thick envelope with temperatures low enough that dust formation occurs. We simulate the NS orbit until only 0.8M$_{\odot}$ of the hydrogen envelope remains around the giant star's core. Our results suggest that the inspiral will continue until another $\approx$0.3M$_{\odot}$ are removed, at which point the remaining envelope will retract. Upon separation, a phase of dynamically stable mass transfer onto the NS accretor is likely to ensue, which may be observable as an ultraluminous X-ray source. The resulting binary, comprised of a detached 2.6M$_{\odot}$ helium-star and a NS with a separation of 3.3-5.7R$_{\odot}$, is expected to evolve into a merging double neutron-star, analogous to those recently detected by LIGO/Virgo. For our chosen combination of binary parameters, our estimated final separation (including the phase of stable mass transfer) suggests a very high $α_{\rm CE}$-equivalent efficiency of $\approx$5. △ Less

Submitted 29 July, 2019; originally announced July 2019.

Comments: 10 pages, 5 figures, submitted to The Astrophysical Journal Letters. Comments welcome

Voyage 2050 White Paper: All-Sky Visible and Near Infrared Space Astrometry

Authors: David Hobbs, Anthony Brown, Erik Høg, Carme Jordi, Daisuke Kawata, Paolo Tanga, Sergei Klioner, Alessandro Sozzetti, Łukasz Wyrzykowski, Nic Walton, Antonella Vallenari, Valeri Makarov, Jan Rybizki, Fran Jiménez-Esteban, José A. Caballero, Paul J. McMillan, Nathan Secrest, Roger Mor, Jeff J. Andrews, Tomaž Zwitter, Cristina Chiappini, Johan P. U. Fynbo, Yuan-Sen Ting, Daniel Hestroffer, Lennart Lindegren , et al. (5 additional authors not shown)

Abstract: A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost… ▽ More A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR. △ Less

Submitted 26 July, 2019; originally announced July 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1609.07325

Did GW170817 harbor a pulsar?

Authors: Enrico Ramirez-Ruiz, Jeff J. Andrews, Sophie L. Schrøder

Abstract: If the progenitor of GW170817 harbored a pulsar, then a Poynting flux dominated bow-shock cavity would have been expected to form around the traveling binary. The characteristic size of this evacuated region depends strongly on the spin-down evolution of the pulsar companion, which in turn depends on the merging timescale of the system. If this evacuated region is able to grow to a sufficiently la… ▽ More If the progenitor of GW170817 harbored a pulsar, then a Poynting flux dominated bow-shock cavity would have been expected to form around the traveling binary. The characteristic size of this evacuated region depends strongly on the spin-down evolution of the pulsar companion, which in turn depends on the merging timescale of the system. If this evacuated region is able to grow to a sufficiently large scale, then the deceleration of the jet, and thus the onset of the afterglow, would be noticeably delayed. The first detection of afterglow emission, which was uncovered 9.2 days after the $γ$-ray burst trigger, can thus be used to constrain the size of a pre-existing pulsar-wind cavity. We use this information, together with a model of the jet to place limits on the presence of a pulsar in GW170817 and discuss the derived constraints in the context of the observed double neutron star binary population. We find that the majority of Galactic systems that are close enough to merge within a Hubble time would have carved a discernibly large pulsar-wind cavity, inconsistent with the onset timescale of the X-ray afterglow of GW170817. Conversely, the recently detected system J1913+1102, which host a low-luminosity pulsar, provides a congruous Milky Way analog of GW170817's progenitor model. This study highlights the potential of the proposed observational test for gaining insight into the origin of double neutron star binaries, in particular if the properties of Galactic systems are representative of the overall merging population. △ Less

Submitted 28 August, 2019; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: Accepted for publication in ApJL, 6 pages, 5 figures

Double Neutron Star Populations and Formation Channels

Authors: Jeff J. Andrews, Ilya Mandel

Abstract: In the past five years, the number of known double neutron stars (DNS) in the Milky Way has roughly doubled. We argue that the observed sample can be split into three distinct sub-populations based on their orbital characteristics: (i) short-period, low-eccentricity binaries; (ii) wide binaries; and (iii) short-period, high-eccentricity binaries. These sub-populations also exhibit distinct spin pe… ▽ More In the past five years, the number of known double neutron stars (DNS) in the Milky Way has roughly doubled. We argue that the observed sample can be split into three distinct sub-populations based on their orbital characteristics: (i) short-period, low-eccentricity binaries; (ii) wide binaries; and (iii) short-period, high-eccentricity binaries. These sub-populations also exhibit distinct spin period and spindown rate properties. We focus on sub-population (iii), which contains the Hulse-Taylor binary. Contrary to previous analysis, we demonstrate that, if they are the product of primordial binary evolution, the $P_{\rm orb}$ and $e$ distribution of these systems requires that the second-born NSs must have been formed with small natal kicks ($\lesssim$25 km s$^{-1}$) and have pre-SN masses narrowly distributed around 3.2 M$_{\odot}$. These constraints challenge binary evolution theory and further predict closely aligned spin and orbital axes, inconsistent with the Hulse-Taylor binary's measured spin-orbit misalignment angle of $\approx$20$^{\circ}$. Motivated by the similarity of these DNSs to B2127+11C, a DNS residing in the globular cluster M15, we argue that this sub-population is consistent with being formed in, and then ejected from, globular clusters. This scenario provides a pathway for the formation and merger of DNSs in stellar environments without recent star formation, as observed in the host galaxy population of short gamma ray bursts and the recent detection by LIGO of a merging DNS in an old stellar population. △ Less

Submitted 8 August, 2019; v1 submitted 29 April, 2019; originally announced April 2019.

Comments: 8 pages, 4 figures, 1 table, accepted for publication in ApJL

Journal ref: 2019ApJ...880L...8A

Double Neutron Star Formation: Merger Times, Systemic Velocities, and Travel Distances

Authors: Jeff J. Andrews, Andreas Zezas

Abstract: The formation and evolution of double neutron stars (DNS) have traditionally been studied using binary population synthesis. In this work, we take an alternative approach by focusing only on the second supernova (SN) forming the DNS and the subsequent orbital decay and merger due to gravitational wave radiation. Using analytic and numerical methods, we explore how different NS natal kick velocity… ▽ More The formation and evolution of double neutron stars (DNS) have traditionally been studied using binary population synthesis. In this work, we take an alternative approach by focusing only on the second supernova (SN) forming the DNS and the subsequent orbital decay and merger due to gravitational wave radiation. Using analytic and numerical methods, we explore how different NS natal kick velocity distributions, pre-SN orbital separations, and progenitor He-star masses affect the post-SN orbital periods, eccentricities, merger times, systemic velocities, and distances traveled by the system before merging. Comparison with the set of 17 known DNSs in the Milky Way shows that DNSs have pre-SN orbital separations ranging between 1 and 44 $R_{\odot}$. Those DNSs with pre-SN separations $\sim$1 $R_{\odot}$ have merger time distributions that peak $\sim$10-100 Myr after formation, regardless of the kick velocity received by the NS. These DNSs are typically formed with systemic velocities $\sim$10$^2$ km s$^{-1}$ and may travel $\sim$1-10 kpc before merging. Depending on progenitor mass of the second-born NS, the short merger time can account for the $r$-process enrichment observed in compact stellar systems such as ultra-faint dwarf galaxies. For Milky Way-mass galaxies only DNSs with the tightest pre-SN orbits and large kick velocities ($\gtrsim$10$^2$ km s$^{-1}$) can escape. However, those DNSs that do escape may travel as far as $\sim$Mpc before merging, which as previous studies have pointed out has implications for identifying the host galaxies to short gamma ray bursts and gravitational wave events. △ Less

Submitted 12 April, 2019; originally announced April 2019.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Using APOGEE Wide Binaries to Test Chemical Tagging with Dwarf Stars

Authors: Jeff J. Andrews, Borja Anguiano, Julio Chanamé, Marcel A. Agüeros, Hannah M. Lewis, Christian R. Hayes, Steven R. Majewski

Abstract: Stars of a common origin are thought to have similar, if not nearly identical, chemistry. Chemical tagging seeks to exploit this fact to identify Milky Way subpopulations through their unique chemical fingerprints. In this work, we compare the chemical abundances of dwarf stars in wide binaries to test the abundance consistency of stars of a common origin. Our sample of 31 wide binaries is identif… ▽ More Stars of a common origin are thought to have similar, if not nearly identical, chemistry. Chemical tagging seeks to exploit this fact to identify Milky Way subpopulations through their unique chemical fingerprints. In this work, we compare the chemical abundances of dwarf stars in wide binaries to test the abundance consistency of stars of a common origin. Our sample of 31 wide binaries is identified from a catalog produced by cross-matching APOGEE stars with UCAC5 astrometry, and we confirm the fidelity of this sample with precision parallaxes from Gaia DR2. For as many as 14 separate elements, we compare the abundances between components of our wide binaries, finding they have very similar chemistry (typically within 0.1 dex). This level of consistency is more similar than can be expected from stars with different origins (which show typical abundance differences of 0.3-0.4 dex within our sample). For the best measured elements, Fe, Si, K, Ca, Mn, and Ni, these differences are reduced to 0.05-0.08 dex when selecting pairs of dwarf stars with similar temperatures. Our results suggest that APOGEE dwarf stars may currently be used for chemical tagging at the level of $\sim$0.1 dex or at the level of $\sim$0.05 dex when restricting for the best-measured elements in stars of similar temperatures. Larger wide binary catalogs may provide calibration sets, in complement to open cluster samples, for on-going spectroscopic surveys. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: 21 pages, 14 figures, accepted for publication in ApJ

Constraining Compact Object Formation with 2M0521

Authors: Katelyn Breivik, Sourav Chatterjee, Jeff J. Andrews

Abstract: We show that the recently discovered binary 2M05215658+4359220 (2M0521), comprised of a giant star (GS) orbiting a suspected black hole (BH) in a ~80 day orbit, may be instrumental in shedding light on uncertain BH-formation physics and can be a test case for studying wind accretion models. Using binary population synthesis with a realistic prescription for the star formation history and metallici… ▽ More We show that the recently discovered binary 2M05215658+4359220 (2M0521), comprised of a giant star (GS) orbiting a suspected black hole (BH) in a ~80 day orbit, may be instrumental in shedding light on uncertain BH-formation physics and can be a test case for studying wind accretion models. Using binary population synthesis with a realistic prescription for the star formation history and metallicity evolution of the Milky Way, we analyze the evolution of binaries containing compact objects (COs) in orbit around GSs with properties similar to 2M0521. We find ~100-1000 CO-GS binaries in the Milky Way observable by Gaia, and 0-12 BH-GS and 0-1 neutron star-GS binaries in the Milky Way with properties similar to 2M0521. We find that all CO-GSs with Porb<5 yr, including 2M0521, go through a common envelope (CE) and hence form a class of higher mass analogs to white dwarf post-CE binaries. We further show how the component masses of 2M0521-like binaries depend strongly on the supernova-engine model we adopt. Thus, an improved measurement of the orbit of 2M0521, imminent with Gaia's third data release, will strongly constrain its component masses and as a result inform supernova-engine models widely used in binary population synthesis studies. These results have widespread implications for the origins and properties of CO binaries, especially those detectable by LIGO and LISA. Finally, we show that the reported X-ray non-detection of 2M0521 is a challenge for wind accretion theory, making 2M0521-like CO-GS binaries a prime target for further study with accretion models. △ Less

Submitted 24 May, 2019; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: 7 pages, 5 figures, Accepted for Publication in ApJL

r-process enrichment of ultra-faint dwarf galaxies by fast merging double neutron stars

Authors: Mohammadtaher Safarzadeh, Enrico Ramirez-Ruiz, Jeff. J. Andrews, Tassos Fragos, Phillip Macias, Evan Scannapieco

Abstract: The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow up kilonova observations have shown that NSMs produce copious amount of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of Double Neutron Stars (DNSs), with the levels of r-process enrichment seen in ultra faint dwarf (UF… ▽ More The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow up kilonova observations have shown that NSMs produce copious amount of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of Double Neutron Stars (DNSs), with the levels of r-process enrichment seen in ultra faint dwarf (UFD) galaxies such as Reticulum II and Tucana III. Assuming that such dwarf systems have lost a significant fraction of their stellar mass through tidal stripping, we conclude that contrary to most current models, it is the DNSs with rather large natal kicks but very short merging timescales that can enrich UFD-type galaxies. These binaries are either on highly eccentric orbits, or form with very short separations due to an additional mass-transfer between the first-born neutron star and a naked helium star, progenitor of the second-born neutron star. These DNSs are born with a frequency that agrees with the statistics of the r-process UFDs, and merge well within the virial radius of their host halos, therefore contributing significantly to their r-process enrichment. △ Less

Submitted 12 January, 2019; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: Accepted for publication in ApJ

A Serendipitous Pulsar Discovery in a Search for a Companion to a Low-Mass White Dwarf

Authors: Jeff J. Andrews, Marcel A. Agüeros, Fernando Camilo, Mukremin Kilic, Alex Gianninas, Warren Brown, Craig Heinke

Abstract: We report the discovery of a previously unidentified pulsar as part of a radio campaign to identify neutron star companions to low-mass white dwarfs (LMWDs) using the Robert C.\ Byrd Green Bank Telescope (GBT). PSR J0802-0955, which is coincident with the position of a WD with a mass of 0.2 solar masses, has a pulse period of 571 ms. Because of its relatively long pulse period, the lack of radial… ▽ More We report the discovery of a previously unidentified pulsar as part of a radio campaign to identify neutron star companions to low-mass white dwarfs (LMWDs) using the Robert C.\ Byrd Green Bank Telescope (GBT). PSR J0802-0955, which is coincident with the position of a WD with a mass of 0.2 solar masses, has a pulse period of 571 ms. Because of its relatively long pulse period, the lack of radial velocity (RV) variations in the radio data, and GBT's large beam size at the observing frequency of 340 MHz, we conclude that PSR J0802-0955 is unassociated with the LMWD at roughly the same position and distance. △ Less

Submitted 15 June, 2018; originally announced June 2018.

Comments: Accepted for publication in Research Notes of the AAS

Validating TGAS wide binaries with Gaia DR2 Radial Velocities and Parallaxes

Authors: Jeff J. Andrews, Julio Chanamé, Marcel A. Agüeros

Abstract: Featuring greatly improved astrometry and precise radial velocities (RVs), the second data release (DR2) from Gaia affords us a unique opportunity to test the validity of samples of stellar wide binaries. In a previous work, we presented a set of wide binaries identified in the joint Tycho-Gaia Astrometric Solution catalog. Here, we use DR2 to confirm the low contamination rate of this sample, the… ▽ More Featuring greatly improved astrometry and precise radial velocities (RVs), the second data release (DR2) from Gaia affords us a unique opportunity to test the validity of samples of stellar wide binaries. In a previous work, we presented a set of wide binaries identified in the joint Tycho-Gaia Astrometric Solution catalog. Here, we use DR2 to confirm the low contamination rate of this sample, thereby verifying the effectiveness of our algorithm for application to DR2 data. △ Less

Submitted 10 May, 2018; originally announced May 2018.

Comments: Accepted for publication in Research Notes of the AAS

Journal ref: RNAAS, 2018, Vol 2, No 2, pg 29

dart_board: Binary Population Synthesis with Markov Chain Monte Carlo

Authors: Jeff J Andrews, Andreas Zezas, Tassos Fragos

Abstract: By employing Monte Carlo random sampling, traditional binary population synthesis (BPS) offers a substantial improvement in efficiency over brute force, grid-based studies. Even so, BPS models typically require a large number of simulation realizations, a computationally expensive endeavor, to generate statistically robust results. Recent advances in statistical methods have led us to revisit the… ▽ More By employing Monte Carlo random sampling, traditional binary population synthesis (BPS) offers a substantial improvement in efficiency over brute force, grid-based studies. Even so, BPS models typically require a large number of simulation realizations, a computationally expensive endeavor, to generate statistically robust results. Recent advances in statistical methods have led us to revisit the traditional approach to BPS. In this work we describe our publicly available code dart_board which combines rapid binary evolution codes, typically used in traditional BPS, with modern Markov chain Monte Carlo methods. dart_board takes a novel approach that treats the initial binary parameters and the supernova kick vector as model parameters. This formulation has several advantages, including the ability to model either populations of systems or individual binaries, the natural inclusion of observational uncertainties, and the flexible addition of new constraints which are problematic to include using traditional BPS. After testing our code with mock systems, we demonstrate the flexibility of dart_board by applying it to three examples: (i) a generic population of high mass X-ray binaries (HMXBs), (ii) the population of HMXBs in the Large Magellanic Cloud (LMC) in which the spatially resolved star formation history is used as a prior, and (iii) one particular HMXB in the LMC, Swift J0513.4-6547, in which we include observations of the system's component masses and orbital period. Although this work focuses on HMXBs, dart_board can be applied to a variety of stellar binaries including the recent detections by gravitational wave observatories of merging compact object binaries. △ Less

Submitted 1 June, 2018; v1 submitted 30 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJS, 39 pages, 21 figures, 5 tables

Wide Binaries in Tycho-Gaia II: Metallicities, Abundances, and Prospects for Chemical Tagging

Authors: Jeff J. Andrews, Julio Chanamé, Marcel A. Agüeros

Abstract: From our recent catalog based on the first Gaia data release (TGAS), we select wide binaries in which both stars have been observed by the Radial Velocity Experiment (RAVE) or the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST). Using RAVE and LAMOST metallicities and RAVE Mg, Al, Si, Ti, and Fe abundances, we find that the differences in the metallicities and elemental abundanc… ▽ More From our recent catalog based on the first Gaia data release (TGAS), we select wide binaries in which both stars have been observed by the Radial Velocity Experiment (RAVE) or the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST). Using RAVE and LAMOST metallicities and RAVE Mg, Al, Si, Ti, and Fe abundances, we find that the differences in the metallicities and elemental abundances of components of wide binaries are consistent with being due to observational uncertainties, in agreement with previous results for smaller and more restricted samples. The metallicity and elemental abundance consistency between wide binary components presented in this work confirms their common origin and bolsters the status of wide binaries as "mini-open clusters." Furthermore, this is evidence that wide binaries are effectively co-eval and co-chemical, supporting their use for e.g., constraining age-activity-rotation relations, the initial-final mass relation for white dwarfs, and M-dwarf metallicity indicators. Additionally, we demonstrate that the common proper motion, common parallax pairs in TGAS with the most extreme separations (s $\gtrsim$ 0.1 pc) typically have inconsistent metallicities, radial velocities or both and are therefore likely to be predominantly comprised of random alignments of unassociated stars with similar astrometry, in agreement with our previous results. Finally, we propose that wide binaries form an ideal data set with which to test chemical tagging as a method to identify stars of common origin, particularly because the stars in wide binaries span a wide range of metallicities, much wider than that spanned by nearby open clusters. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: 15 pages, 11 figures, 3 tables, accepted for publication in MNRAS

Wide Binaries in Tycho-{\it Gaia}: Search Method and the Distribution of Orbital Separations

Authors: Jeff J. Andrews, Julio Chanamé, Marcel A. Agüeros

Abstract: We mine the Tycho-{\it Gaia} astrometric solution (TGAS) catalog for wide stellar binaries by matching positions, proper motions, and astrometric parallaxes. We separate genuine binaries from unassociated stellar pairs through a Bayesian formulation that includes correlated uncertainties in the proper motions and parallaxes. Rather than relying on assumptions about the structure of the Galaxy, we… ▽ More We mine the Tycho-{\it Gaia} astrometric solution (TGAS) catalog for wide stellar binaries by matching positions, proper motions, and astrometric parallaxes. We separate genuine binaries from unassociated stellar pairs through a Bayesian formulation that includes correlated uncertainties in the proper motions and parallaxes. Rather than relying on assumptions about the structure of the Galaxy, we calculate Bayesian priors and likelihoods based on the nature of Keplerian orbits and the TGAS catalog itself. We calibrate our method using radial velocity measurements and obtain 6196 high-confidence candidate wide binaries with projected separations $s\lesssim1$ pc. The normalization of this distribution suggests that at least 0.6\% of TGAS stars have an associated, distant TGAS companion in a wide binary. We demonstrate that {\it Gaia}'s astrometry is precise enough that it can detect projected orbital velocities in wide binaries with orbital periods as large as 10$^6$ yr. For pairs with $s\ \lesssim\ 4\times10^4$~AU, characterization of random alignments indicate our contamination to be $\approx$5\%. For $s \lesssim 5\times10^3$~AU, our distribution is consistent with Öpik's Law. At larger separations, the distribution is steeper and consistent with a power-law $P(s)\propto s^{-1.6}$; there is no evidence in our data of any bimodality in this distribution for $s \lesssim$ 1 pc. Using radial velocities, we demonstrate that at large separations, i.e., of order $s \sim$ 1 pc and beyond, any potential sample of genuine wide binaries in TGAS cannot be easily distinguished from ionized former wide binaries, moving groups, or contamination from randomly aligned stars. △ Less

Submitted 25 April, 2017; originally announced April 2017.

Comments: 27 pages, 17 figures, 2 tables, submitted to MNRAS

A closer look at the "characteristic" width of molecular cloud filaments

Authors: G. V. Panopoulou, I. Psaradaki, R. Skalidis, K. Tassis, J. J. Andrews

Abstract: Filaments in Herschel molecular cloud images are found to exhibit a "characteristic width". This finding is in tension with spatial power spectra of the data, which show no indication of this characteristic scale. We demonstrate that this discrepancy is a result of the methodology adopted for measuring filament widths. First, we perform the previously used analysis technique on artificial scale-fr… ▽ More Filaments in Herschel molecular cloud images are found to exhibit a "characteristic width". This finding is in tension with spatial power spectra of the data, which show no indication of this characteristic scale. We demonstrate that this discrepancy is a result of the methodology adopted for measuring filament widths. First, we perform the previously used analysis technique on artificial scale-free data, and obtain a peaked width distribution of filament-like structures. Next, we repeat the analysis on three Herschel maps and reproduce the narrow distribution of widths found in previous studies $-$ when considering the average width of each filament. However, the distribution of widths measured at all points along a filament spine is broader than the distribution of mean filament widths, indicating that the narrow spread (interpreted as a "characteristic" width) results from averaging. Furthermore, the width is found to vary significantly from one end of a filament to the other. Therefore, the previously identified peak at 0.1 pc cannot be understood as representing the typical width of filaments. We find an alternative explanation by modelling the observed width distribution as a truncated power-law distribution, sampled with uncertainties. The position of the peak is connected to the lower truncation scale and is likely set by the choice of parameters used in measuring filament widths. We conclude that a "characteristic" width of filaments is not supported by the available data. △ Less

Submitted 22 November, 2016; originally announced November 2016.

Comments: 15 pages (incl. appendices), 16 figures, accepted to MNRAS

Today a Duo, But Once a Trio? The Double White Dwarf HS 2220$+$2146 May Be A Post-Blue Straggler Binary

Authors: Jeff J. Andrews, Marcel Agüeros, Warren R. Brown, Natalie M. Gosnell, A. Gianninas, Mukremin Kilic, Detlev Koester

Abstract: For sufficiently wide orbital separations {\it a}, the two members of a stellar binary evolve independently. This implies that in a wide double white dwarf (DWD), the more massive WD should always be produced first, when its more massive progenitor ends its main-sequence life, and should therefore be older and cooler than its companion. The bound, wide DWD HS 2220$+$2146 ($a\approx500$ AU) does no… ▽ More For sufficiently wide orbital separations {\it a}, the two members of a stellar binary evolve independently. This implies that in a wide double white dwarf (DWD), the more massive WD should always be produced first, when its more massive progenitor ends its main-sequence life, and should therefore be older and cooler than its companion. The bound, wide DWD HS 2220$+$2146 ($a\approx500$ AU) does not conform to this picture: the more massive WD is the younger, hotter of the pair. We show that this discrepancy is unlikely to be due to past mass-transfer phases or to the presence of an unresolved companion. Instead, we propose that HS 2220$+$2146 formed through a new wide DWD evolutionary channel involving the merger of the inner binary in a hierarchical triple system. The resulting blue straggler and its wide companion then evolved independently, forming the WD pair seen today. Although we cannot rule out other scenarios, the most likely formation channel has the inner binary merging while both stars are still on the main sequence. This provides us with the tantalizing possibility that Kozai-Lidov oscillations may have played a role in the inner binary's merger. {\it Gaia} may uncover hundreds more wide DWDs, leading to the identification of other systems like HS 2220$+$2146. There are already indications that other WD systems may have been formed through different, but related, hierarchical triple evolutionary scenarios. Characterizing these populations may allow for thorough testing of the efficiency with which KL oscillations induce stellar mergers. △ Less

Submitted 1 July, 2016; originally announced July 2016.

Comments: Accepted for publication in ApJ, 12 pages, 6 figures, 2 tables

Constraints on the Initial-Final Mass Relation from Wide Double White Dwarfs

Authors: Jeff J. Andrews, Marcel A. Agüeros, A. Gianninas, Mukremin Kilic, Saurav Dhital, Scott F. Anderson

Abstract: We present observational constraints on the initial-final mass relation (IFMR) using wide double white dwarfs (DWDs). We identify 65 new candidate wide DWDs within the Sloan Digital Sky Survey, bringing the number of candidate wide DWDs to 142. We then engage in a spectroscopic follow-up campaign and collect existing spectra for these objects; using these spectra, we derive masses and cooling ages… ▽ More We present observational constraints on the initial-final mass relation (IFMR) using wide double white dwarfs (DWDs). We identify 65 new candidate wide DWDs within the Sloan Digital Sky Survey, bringing the number of candidate wide DWDs to 142. We then engage in a spectroscopic follow-up campaign and collect existing spectra for these objects; using these spectra, we derive masses and cooling ages for 54 hydrogen (DA) WDs in DWDs. We also identify one new DA/DB pair, four candidate DA/DC pairs, four candidate DA/DAH pairs, and one new candidate triple degenerate system. Because wide DWDs are co-eval and evolve independently, the difference in the pre-WD lifetimes should equal the difference in the WD cooling ages. We use this to develop a Bayesian hierarchical framework and construct a likelihood function to determine the probability that any particular IFMR fits a sample of wide DWDs. We then define a parametric model for the IFMR and find the best parameters indicated by our sample of DWDs. We place robust constraints on the IFMR for initial masses of 2--4 \Msun. The WD masses produced by our model for stars within this mass range differ from those predicted by semi-empirical fits to open cluster WDs. Within this mass range, where there are few constraining open cluster WDs and disagreements in the cluster ages, wide DWDs may provide more reliable constraints on the IFMR. Expanding this method to the many wide DWDs expected to be discovered by Gaia may transform our understanding of the IFMR. △ Less

Submitted 20 October, 2015; originally announced October 2015.

Comments: Accepted for publication in the ApJ. Code to reproduce the IFMR constraints is hosted on GitHub (https://github.com/astroJeff/separated-at-birth). Samples from our posterior distribution of IFMRs are available on figshare at doi:10.6084/m9.figshare.1572148

Ultracool white dwarfs and the age of the Galactic disc

Authors: A. Gianninas, B. Curd, John R. Thorstensen, Mukremin Kilic, P. Bergeron, Jeff J. Andrews, Paul Canton, M. A. Agüeros

Abstract: We present parallax observations and a detailed model atmosphere analysis of 54 cool and ultracool ($T_{\rm eff}$ < 4000 K) white dwarfs (WDs) in the solar neighbourhood. For the first time, a large number of cool and ultracool WDs have distance and tangential velocities measurements available. Our targets have distances ranging from 21 pc to >100 pc, and include five stars within 30 pc. Contrary… ▽ More We present parallax observations and a detailed model atmosphere analysis of 54 cool and ultracool ($T_{\rm eff}$ < 4000 K) white dwarfs (WDs) in the solar neighbourhood. For the first time, a large number of cool and ultracool WDs have distance and tangential velocities measurements available. Our targets have distances ranging from 21 pc to >100 pc, and include five stars within 30 pc. Contrary to expectations, all but two of them have tangential velocities smaller than 150 km s$^{-1}$ thus suggesting Galactic disc membership. The oldest WDs in this sample have WD cooling ages of 10 Gyr, providing a firm lower limit to the age of the thick disc population. Many of our targets have uncharacteristically large radii, indicating that they are low mass WDs. It appears that we have detected the brighter population of cool and ultracool WDs near the Sun. The fainter population of ultracool CO-core WDs remain to be discovered in large numbers. The Large Synoptic Survey Telescope should find these elusive, more massive ultracool WDs in the solar neighbourhood. △ Less

Submitted 27 March, 2015; v1 submitted 10 March, 2015; originally announced March 2015.

Comments: 16 pages, 11 figures, 4 tables, accepted for publication in MNRAS