Showing 1–50 of 89 results for author: Fragos, T

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!

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

Computing Challenges for the Einstein Telescope project

Authors: Stefano Bagnasco, Antonella Bozzi, Tassos Fragos, Alba Gonzalvez, Steffen Hahn, Gary Hemming, Lia Lavezzi, Paul Laycock, Gonzalo Merino, Silvio Pardi, Steven Schramm, Achim Stahl, Andres Tanasijczuk, Nadia Tonello, Sara Vallero, John Veitch, Patrice Verdier

Abstract: The discovery of gravitational waves, first observed in September 2015 following the merger of a binary black hole system, has already revolutionised our understanding of the Universe. This was further enhanced in August 2017, when the coalescence of a binary neutron star system was observed both with gravitational waves and a variety of electromagnetic counterparts; this joint observation marked… ▽ More The discovery of gravitational waves, first observed in September 2015 following the merger of a binary black hole system, has already revolutionised our understanding of the Universe. This was further enhanced in August 2017, when the coalescence of a binary neutron star system was observed both with gravitational waves and a variety of electromagnetic counterparts; this joint observation marked the beginning of gravitational multimessenger astronomy. The Einstein Telescope, a proposed next-generation ground-based gravitational-wave observatory, will dramatically increase the sensitivity to sources: the number of observations of gravitational waves is expected to increase from roughly 100 per year to roughly 100'000 per year, and signals may be visible for hours at a time, given the low frequency cutoff of the planned instrument. This increase in the number of observed events, and the duration with which they are observed, is hugely beneficial to the scientific goals of the community but poses a number of significant computing challenges. Moreover, the currently used computing algorithms do not scale to this new environment, both in terms of the amount of resources required and the speed with which each signal must be characterised. This contribution will discuss the Einstein Telescope's computing challenges, and the activities that are underway to prepare for them. Available computing resources and technologies will greatly evolve in the years ahead, and those working to develop the Einstein Telescope data analysis algorithms will need to take this into account. It will also be important to factor into the initial development of the experiment's computing model the availability of huge parallel HPC systems and ubiquitous Cloud computing; the design of the model will also, for the first time, include the environmental impact as one of the optimisation metrics. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: Submitted for the proceedings of CHEP2023

Dynamical tides in binaries: Inconsistencies in the implementation of Zahn's prescription

Authors: Luca Sciarini, Sylvia Ekström, Patrick Eggenberger, Georges Meynet, Tassos Fragos, Han Feng Song

Abstract: Binary evolution codes are essential tools to help in understanding the evolution of binary systems. They contain a great deal of physics, for example stellar evolution, stellar interactions, mass transfer, tides, orbital evolution. Since many of these processes are difficult to account for in detail, we often rely on prescriptions obtained in earlier studies. We highlight that the impact of the d… ▽ More Binary evolution codes are essential tools to help in understanding the evolution of binary systems. They contain a great deal of physics, for example stellar evolution, stellar interactions, mass transfer, tides, orbital evolution. Since many of these processes are difficult to account for in detail, we often rely on prescriptions obtained in earlier studies. We highlight that the impact of the dynamical tides with radiative damping has been implemented inconsistently with respect to its original theoretical formulation in many studies. We derive a new analytical solution for the evolution toward synchronization in the case of circular orbits and propose turnkey equations for the case of eccentric orbits that can be used in population synthesis studies. We compare the strength of the tidal torque obtained with this new formula with respect to that obtained with the formula generally used in literature by studying how the evolution toward synchronization of main sequence stellar models is affected. We conclude that by using an incorrect formula for the tidal torque, as has been done in many binary codes, the strength of the dynamical tides with radiative damping is over- or underestimated depending on whether the star is close to or far from synchronization. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Report number: aa48424-23

Multi-messenger prospects for black hole - neutron star mergers in the O4 and O5 runs

Authors: Alberto Colombo, Raphaël Duqué, Om Sharan Salafia, Floor S. Broekgaarden, Francesco Iacovelli, Michele Mancarella, Igor Andreoni, Francesco Gabrielli, Fabio Ragosta, Giancarlo Ghirlanda, Tassos Fragos, Andrew J. Levan, Silvia Piranomonte, Andrea Melandri, Bruno Giacomazzo, Monica Colpi

Abstract: The existence of merging black hole-neutron star (BHNS) binaries has been ascertained through the observation of their gravitational wave (GW) signals. However, to date, no definitive electromagnetic (EM) emission has been confidently associated with these mergers. Such an association could help unravel crucial information on these systems, for example, their BH spin distribution, the equation of… ▽ More The existence of merging black hole-neutron star (BHNS) binaries has been ascertained through the observation of their gravitational wave (GW) signals. However, to date, no definitive electromagnetic (EM) emission has been confidently associated with these mergers. Such an association could help unravel crucial information on these systems, for example, their BH spin distribution, the equation of state (EoS) of NS and the rate of heavy element production. We model the multi-messenger (MM) emission from BHNS mergers detectable during the fourth (O4) and fifth (O5) observing runs of the LIGO-Virgo-KAGRA GW detector network, in order to provide detailed predictions that can help enhance the effectiveness of observational efforts and extract the highest possible scientific information from such remarkable events. Our methodology is based on a population synthesis-approach, which includes the modelling of the signal-to-noise ratio of the GW signal in the detectors, the GW-inferred sky localization of the source, the kilonova (KN) optical and near-infrared light curves, the relativistic jet gamma-ray burst (GRB) prompt emission peak photon flux, and the GRB afterglow light curves in the radio, optical and X-ray bands. The resulting prospects for BHNS MM detections during O4 are not promising, with a GW detection rate of $15.0^{+15.4}_{-8.8}$ yr$^{-1}$, but joint MM rates of $\sim 10^{-1}$ yr$^{-1}$ for the KN and $\sim 10^{-2}$ yr$^{-1}$ for the jet-related emission. In O5 we find an overall increase in expected detection rates by around an order of magnitude, owing to both the enhanced sensitivity of the GW detector network, and the coming online of future EM facilities. Finally, we discuss direct searches for the GRB radio afterglow with large-field-of-view instruments as a new possible follow-up strategy in the context of ever-dimming prospects for KN detection. △ Less

Submitted 12 April, 2024; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: Submitted to A&A. 17 pages, 11 figures, 2 tables. Comments are welcome!

Journal ref: A&A 686, A265 (2024)

Detecting the Undetected: Overcoming Biases in Gravitational-Wave Population Studies

Authors: Ryan Raikman, Simone Bavera, Tassos Fragos

Abstract: In the flourishing field of gravitational-wave astronomy, accurately inferring binary black hole merger formation channels is paramount. The Bayesian hierarchical model selection analysis offers a promising methodology (see, e.g., One Channel to Rule Them All, Zevin et al. 2021). However, recently, Cheng et al. (2023) highlighted a critical caveat: observed channels absent in known models can bias… ▽ More In the flourishing field of gravitational-wave astronomy, accurately inferring binary black hole merger formation channels is paramount. The Bayesian hierarchical model selection analysis offers a promising methodology (see, e.g., One Channel to Rule Them All, Zevin et al. 2021). However, recently, Cheng et al. (2023) highlighted a critical caveat: observed channels absent in known models can bias branching fraction estimates. In this research note, we introduce a test to detect missing channels in such analyses. Our findings show a commendable success rate in identifying these elusive channels. Yet, in scenarios where missing channels closely overlap with recognized ones, discerning the difference remains challenging. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 4 pages, 1 figure, submitted to research notes AAS

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

N-enhancement in GN-z11: First evidence for supermassive stars nucleosynthesis in proto-globular clusters-like conditions at high redshift ?

Authors: C. Charbonnel, D. Schaerer, N. Prantzos, L. Ramírez-Galeano, T. Fragos, A. Kuruvandothi, R. Marques-Chaves, M. Gieles

Abstract: Unusually high N/O abundance ratios were recently reported for a very compact, intensively star-forming object GN-z11 at $z=10.6$ from JWST/NIRSpec observations. We present an empirical comparison with the C, N, and O abundance ratios in Galactic globular clusters (GCs) over a large metallicity range. We show that hot hydrogen-burning nucleosynthesis within supermassive stars (SMS) formed through… ▽ More Unusually high N/O abundance ratios were recently reported for a very compact, intensively star-forming object GN-z11 at $z=10.6$ from JWST/NIRSpec observations. We present an empirical comparison with the C, N, and O abundance ratios in Galactic globular clusters (GCs) over a large metallicity range. We show that hot hydrogen-burning nucleosynthesis within supermassive stars (SMS) formed through runaway collisions can consistently explain the observed abundances ratio in GN-z11 and in GCs. This suggests that a proto-globular cluster hosting a SMS could be at the origin of the strong N-enrichment in GN-z11. Our model predicts the behavior of N/O, C/O, and Ne/O ratios as a function of metallicity, which can be tested if high-$z$ objects similar to GN-z11 are detected with JWST in the future. Further studies and statistics will help differentiate the proto-GC scenario from the Wolf-Rayet scenario that we quantify with a population synthesis model, and shed more light on this peculiar object. △ Less

Submitted 17 April, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: Accepted by A&A Letters (revised version; in press)

Journal ref: A&A 673, L7 (2023)

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

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)

The ionizing and heating power of ultraluminous X-ray sources under the geometrical beaming model

Authors: K. Kovlakas, T. Fragos, D. Schaerer, A. Mesinger

Abstract: While there is now a consensus that X-ray binaries (XRBs) are the dominant X-ray sources in the early Universe and play a significant role during the epoch of heating of the intergalactic medium (IGM), recent studies report contradicting results regarding their contribution in the nebular emission of local Universe galaxies. Ultraluminous X-ray sources (ULXs), which dominate the X-ray budget of no… ▽ More While there is now a consensus that X-ray binaries (XRBs) are the dominant X-ray sources in the early Universe and play a significant role during the epoch of heating of the intergalactic medium (IGM), recent studies report contradicting results regarding their contribution in the nebular emission of local Universe galaxies. Ultraluminous X-ray sources (ULXs), which dominate the X-ray budget of normal galaxies, may be important interstellar-medium (ISM) ionizing sources. However, their output in the extreme UV (EUV) and soft--X-ray part of the spectrum remains observationally unconstrained. In this paper, we predict the ionizing and heating power from ULX populations under the geometrical beaming scenario, and three models describing the emission from super-critical accretion disks. We find that our theoretical spectra for ULX populations cannot (can) explain the HeII (NeV) emission observed in some galaxies, with their contribution being less (more) important than the underlying stellar population. Stochastic fluctuations in the number of ULXs may allow for equal contributions in the HeII emission, in a fraction of galaxies. We provide average spectra of ULX populations as an input to local, and early-Universe studies. We find that the soft--X-ray emission arising from super-critical accretion is significant for the heating of the IGM, and consistent with recent constraints from the 21-cm cosmic signal. Based on the dependence on the adopted compact-object (CO) mass and accretion model, we encourage efforts in modeling ULX spectra via simulations, and their combination with detailed binary population synthesis models. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: 8 pages. 6 figures, 1 table. This article has been accepted for publication in the Astronomy & Astrophysics journal. The full Table 1 will be made available via CDS

Journal ref: A&A 665, A28 (2022)

The $χ_\mathrm{eff}-z$ correlation of field binary black hole mergers and how 3G gravitational-wave detectors can constrain it

Authors: Simone S. Bavera, Maya Fishbach, Michael Zevin, Emmanouil Zapartas, Tassos Fragos

Abstract: Understanding the origin of merging binary black holes is currently one of the most pressing quests in astrophysics. We show that if isolated binary evolution dominates the formation mechanism of merging binary black holes, one should expect a correlation between the effective spin parameter, $χ_\mathrm{eff}$, and the redshift of the merger, $z$, of binary black holes. This correlation comes from… ▽ More Understanding the origin of merging binary black holes is currently one of the most pressing quests in astrophysics. We show that if isolated binary evolution dominates the formation mechanism of merging binary black holes, one should expect a correlation between the effective spin parameter, $χ_\mathrm{eff}$, and the redshift of the merger, $z$, of binary black holes. This correlation comes from tidal spin-up systems preferentially forming and merging at higher redshifts due to the combination of weaker orbital expansion from low metallicity stars given their reduced wind mass loss rate, delayed expansion and have smaller maximal radii during the supergiant phase compared to stars at higher metallicity. As a result, these tightly bound systems merge with short inspiral times. Given our fiducial model of isolated binary evolution, we show that the origin of a $χ_\mathrm{eff}-z$ correlation in the detectable LIGO--Virgo binary black hole population is different from the intrinsic population, which will become accessible only in the future by third-generation gravitational-wave detectors such as Einstein Telescope and Cosmic Explorer. Finally, we compare our model predictions with population predictions based on the current catalog of binary black hole mergers and find that current data favor a positive correlation of $χ_\mathrm{eff}-z$ as predicted by our model of isolated binary evolution. △ Less

Submitted 27 June, 2022; v1 submitted 6 April, 2022; originally announced April 2022.

Comments: 14 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 665, A59 (2022)

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

The 21-cm signal from the Cosmic Dawn: metallicity dependence of high mass X-ray binaries

Authors: Harman Deep Kaur, Yuxiang Qin, Andrei Mesinger, Andrea Pallottini, Tassos Fragos, Antara Basu-Zych

Abstract: X-rays from High-Mass X-ray Binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; $z\sim 10-15$) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift… ▽ More X-rays from High-Mass X-ray Binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; $z\sim 10-15$) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift observations, as well as theoretical models, imply that the integrated X-ray luminosity to star formation rate of HMXBs ($L_{\rm X}/{\rm SFR}$) should increase in metal-poor environments, typical of early galaxies. Here we study the impact of the metallicity ($Z$) dependence of $L_{\rm X}/{\rm SFR}$ during the EoH. For our fiducial models, galaxies with star formation rates of order $10^{-3} - 10^{-1}$ $M_\odot$ yr$^{-1}$ and metallicities of order $10^{-3} - 10^{-2}$ $Z_\odot$ are the dominant contributors to the X-ray background (XRB) during this period. Different $L_{\rm X}/{\rm SFR}$-$Z$ relations result in factors of $\sim$ 3 differences in these ranges, as well as in the mean IGM temperature and the large-scale 21-cm power, at a given redshift. We compute mock 21-cm observations adopting as a baseline a 1000h integration with the upcoming Square Kilometer Array (SKA), for two different $L_{\rm X}/{\rm SFR}$-$Z$ relations. We perform inference on these mock observations using the common simplification of a constant $L_{\rm X}/{\rm SFR}$, finding that constant $L_{\rm X}/{\rm SFR}$ models can recover the IGM evolution of the more complicated $L_{\rm X}/{\rm SFR}$-$Z$ simulations only during the EoH. At $z<10$, where the typical galaxies are more polluted, constant $L_{\rm X}/{\rm SFR}$ models over-predict the XRB and its relative contribution to the early stages of the reionization. △ Less

Submitted 3 May, 2022; v1 submitted 21 March, 2022; originally announced March 2022.

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

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

The Stellar Age Dependence of X-ray Emission from Normal Star-Forming Galaxies in the GOODS Fields

Authors: Woodrow Gilbertson, Bret Lehmer, Keith Doore, Rafael Eufrasio, Antara Basu-Zych, William Brandt, Tassos Fragos, Kristen Garofali, Konstantinos Kovlakas, Bin Luo, Paolo Tozzi, Fabio Vito, Benjamin F. Williams, Yongquan Xue

Abstract: The Chandra Deep Field-South and North surveys (CDFs) provide unique windows into the cosmic history of X-ray emission from normal (non-active) galaxies. Scaling relations of normal galaxy X-ray luminosity (L_X) with star formation rate (SFR) and stellar mass (M_star) have been used to show that the formation rates of low-mass and high-mass X-ray binaries (LMXBs and HMXBs, respectively) evolve wit… ▽ More The Chandra Deep Field-South and North surveys (CDFs) provide unique windows into the cosmic history of X-ray emission from normal (non-active) galaxies. Scaling relations of normal galaxy X-ray luminosity (L_X) with star formation rate (SFR) and stellar mass (M_star) have been used to show that the formation rates of low-mass and high-mass X-ray binaries (LMXBs and HMXBs, respectively) evolve with redshift across z = 0-2 following L_HMXB/SFR ~ 1 + z and L_LMXB/M_star ~ (1 + z)^{2-3}. However, these measurements alone do not directly reveal the physical mechanisms behind the redshift evolution of X-ray binaries (XRBs). We derive star-formation histories for a sample of 344 normal galaxies in the CDFs, using spectral energy distribution (SED) fitting of FUV-to-FIR photometric data, and construct a self-consistent, age-dependent model of the X-ray emission from the galaxies. Our model quantifies how X-ray emission from hot gas and XRB populations vary as functions of host stellar-population age. We find that (1) the ratio L_X/M_star declines by a factor of ~1000 from 0-10 Gyr and (2) the X-ray SED becomes harder with increasing age, consistent with a scenario in which the hot gas contribution to the X-ray SED declines quickly for ages above 10 Myr. When dividing our sample into subsets based on metallicity, we find some indication that L_X/M_star is elevated for low-metallicity galaxies, consistent with recent studies of X-ray scaling relations. However, additional statistical constraints are required to quantify both the age and metallicity dependence of X-ray emission from star-forming galaxies. △ Less

Submitted 6 December, 2021; originally announced December 2021.

Stochastic gravitational-wave background as a tool to investigate multi-channel astrophysical and primordial black-hole mergers

Authors: Simone S. Bavera, Gabriele Franciolini, Giulia Cusin, Antonio Riotto, Michael Zevin, Tassos Fragos

Abstract: The formation of merging binary black holes can occur through multiple astrophysical channels such as, e.g., isolated binary evolution and dynamical formation or, alternatively, have a primordial origin. Increasingly large gravitational-wave catalogs of binary black-hole mergers have allowed for the first model selection studies between different theoretical predictions to constrain some of their… ▽ More The formation of merging binary black holes can occur through multiple astrophysical channels such as, e.g., isolated binary evolution and dynamical formation or, alternatively, have a primordial origin. Increasingly large gravitational-wave catalogs of binary black-hole mergers have allowed for the first model selection studies between different theoretical predictions to constrain some of their model uncertainties and branching ratios. In this work, we show how one could add an additional and independent constraint to model selection by using the stochastic gravitational-wave background. In contrast to model selection analyses that have discriminating power only up to the gravitational-wave detector horizons (currently at redshifts $z\lesssim 1$ for LIGO-Virgo), the stochastic gravitational-wave background accounts for the redshift integration of all gravitational-wave signals in the Universe. As a working example, we consider the branching ratio results from a model selection study that includes potential contribution from astrophysical and primordial channels. We renormalize the relative contribution of each channel to the detected event rate to compute the total stochastic gravitational-wave background energy density. The predicted amplitude lies below the current observational upper limits of GWTC-2 by LIGO-Virgo, indicating that the results of the model selection analysis are not ruled out by current background limits. Furthermore, given the set of population models and inferred branching ratios, we find that, even though the predicted background will not be detectable by current generation gravitational-wave detectors, it will be accessible by third-generation detectors such as the Einstein Telescope and space-based detectors such as LISA. △ Less

Submitted 21 December, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 12 pages, 4 figures, accepted in A&A

The Evolution of Binaries under the Influence of Radiation-Driven Winds from a Stellar Companion

Authors: Sophie Lund Schrøder, Morgan MacLeod, Enrico Ramirez-Ruiz, Ilya Mandel, Tassos Fragos, Abraham Loeb, Rosa Wallace Everson

Abstract: Interacting binaries are of general interest as laboratories for investigating the physics of accretion, which gives rise to the bulk of high-energy radiation in the Galaxy. They allow us to probe stellar evolution processes that cannot be studied in single stars. Understanding the orbital evolution of binaries is essential in order to model the formation of compact binaries. Here we focus our att… ▽ More Interacting binaries are of general interest as laboratories for investigating the physics of accretion, which gives rise to the bulk of high-energy radiation in the Galaxy. They allow us to probe stellar evolution processes that cannot be studied in single stars. Understanding the orbital evolution of binaries is essential in order to model the formation of compact binaries. Here we focus our attention on studying orbital evolution driven by angular momentum loss through stellar winds in massive binaries. We run a suite of hydrodynamical simulations of binary stars hosting one mass losing star with varying wind velocity, mass ratio, wind velocity profile and adiabatic index, and compare our results to analytic estimates for drag and angular momentum loss. We find that, at leading order, orbital evolution is determined by the wind velocity and the binary mass ratio. Small ratios of wind to orbital velocities and large accreting companion masses result in high angular momentum loss and a shrinking of the orbit. For wider binaries and binaries hosting lighter mass-capturing companions, the wind mass-loss becomes more symmetric, which results in a widening of the orbit. We present a simple analytic formula that can accurately account for angular momentum losses and changes in the orbit, which depends on the wind velocity and mass ratio. As an example of our formalism, we compare the effects of tides and winds in driving the orbital evolution of high mass X-ray binaries, focusing on Vela X-1 and Cygnus X-1 as examples. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: 22 pages, 16 figures, submitted to ApJ, 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

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)

Approximations to the spin of close Black-hole-Wolf-Rayet binaries

Authors: Simone S. Bavera, Michael Zevin, Tassos Fragos

Abstract: Population synthesis studies of binary black-hole mergers often lack robust black-hole spin estimates as they cannot accurately follow tidal spin-up during the late black-hole-Wolf-Rayet evolutionary phase. We provide an analytical approximation of the dimensionless second-born black-hole spin given the binary orbital period and Wolf-Rayet stellar mass at helium depletion or carbon depletion. Thes… ▽ More Population synthesis studies of binary black-hole mergers often lack robust black-hole spin estimates as they cannot accurately follow tidal spin-up during the late black-hole-Wolf-Rayet evolutionary phase. We provide an analytical approximation of the dimensionless second-born black-hole spin given the binary orbital period and Wolf-Rayet stellar mass at helium depletion or carbon depletion. These approximations are obtained from fitting a sample of around $10^5$ detailed MESA simulations that follow the evolution and spin up of close black-hole--Wolf-Rayet systems with metallicities in the range $[10^{-4},1.5Z_\odot]$. Following the potential spin up of the Wolf-Rayet progenitor, the second-born black-hole spin is calculated using up-to-date core collapse prescriptions that account for any potential disk formation in the collapsing Wolf-Rayet star. The fits for second-born black hole spin provided in this work can be readily applied to any astrophysical modeling that relies on rapid population synthesis, and will be useful for the interpretation of gravitational-wave sources using such models. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 3 pages, 1 figure, submitted to AAS journal

The role of core-collapse physics in the observability of black-hole neutron-star mergers as multi-messenger sources

Authors: Jaime Román-Garza, Simone S. Bavera, Tassos Fragos, Emmanouil Zapartas, Devina Misra, Jeff Andrews, Scotty Coughlin, Aaron Dotter, Konstantinos Kovlakas, Juan Gabriel Serra, Ying Qin, Kyle A. Rocha, Nam Hai Tran

Abstract: Recent detailed 1D core-collapse simulations have brought new insights on the final fate of massive stars, which are in contrast to commonly used parametric prescriptions. In this work, we explore the implications of these results to the formation of coalescing black-hole (BH) - neutron-star (NS) binaries, such as the candidate event GW190426_152155 reported in GWTC-2. Furthermore, we investigate… ▽ More Recent detailed 1D core-collapse simulations have brought new insights on the final fate of massive stars, which are in contrast to commonly used parametric prescriptions. In this work, we explore the implications of these results to the formation of coalescing black-hole (BH) - neutron-star (NS) binaries, such as the candidate event GW190426_152155 reported in GWTC-2. Furthermore, we investigate the effects of natal kicks and the NS's radius on the synthesis of such systems and potential electromagnetic counterparts linked to them. Synthetic models based on detailed core-collapse simulations result in an increased merger detection rate of BH-NS systems ($\sim 2.3$ yr$^{-1}$), 5 to 10 times larger than the predictions of "standard" parametric prescriptions. This is primarily due to the formation of low-mass BH via direct collapse, and hence no natal kicks, favored by the detailed simulations. The fraction of observed systems that will produce an electromagnetic counterpart, with the detailed supernova engine, ranges from $2$-$25$%, depending on uncertainties in the NS equation of state. Notably, in most merging systems with electromagnetic counterparts, the NS is the first-born compact object, as long as the NS's radius is $\lesssim 12\,\mathrm{km}$. Furthermore, core-collapse models that predict the formation of low-mass BHs with negligible natal kicks increase the detection rate of GW190426_152155-like events to $\sim 0.6 \, $yr$^{-1}$; with an associated probability of electromagnetic counterpart $\leq 10$% for all supernova engines. However, increasing the production of direct-collapse low-mass BHs also increases the synthesis of binary BHs, over-predicting their measured local merger density rate. In all cases, models based on detailed core-collapse simulation predict a ratio of BH-NSs to binary BHs merger rate density that is at least twice as high as other prescriptions. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 13 pages, 5 figures, 3 tables

Journal ref: The Astrophysical Journal Letters, 2021, vol. 912, no 2, p. L23

One Channel to Rule Them All? Constraining the Origins of Binary Black Holes using Multiple Formation Pathways

Authors: Michael Zevin, Simone S. Bavera, Christopher P. L. Berry, Vicky Kalogera, Tassos Fragos, Pablo Marchant, Carl L. Rodriguez, Fabio Antonini, Daniel E. Holz, Chris Pankow

Abstract: The second LIGO-Virgo catalog of gravitational wave transients has more than quadrupled the observational sample of binary black holes. We analyze this catalog using a suite of five state-of-the-art binary black hole population models covering a range of isolated and dynamical formation channels and infer branching fractions between channels as well as constraints on uncertain physical processes t… ▽ More The second LIGO-Virgo catalog of gravitational wave transients has more than quadrupled the observational sample of binary black holes. We analyze this catalog using a suite of five state-of-the-art binary black hole population models covering a range of isolated and dynamical formation channels and infer branching fractions between channels as well as constraints on uncertain physical processes that impact the observational properties of mergers. Given our set of formation models, we find significant differences between the branching fractions of the underlying and detectable populations, and that the diversity of detections suggests that multiple formation channels are at play. A mixture of channels is strongly preferred over any single channel dominating the detected population: an individual channel does not contribute to more than $\simeq 70\%$ of the observational sample of binary black holes. We calculate the preference between the natal spin assumptions and common envelope efficiencies in our models, favoring natal spins of isolated black holes of $\lesssim 0.1$, and marginally preferring common envelope efficiencies of $\gtrsim 2.0$ while strongly disfavoring highly inefficient common envelopes. We show that it is essential to consider multiple channels when interpreting gravitational wave catalogs, as inference on branching fractions and physical prescriptions becomes biased when contributing formation scenarios are not considered or incorrect physical prescriptions are assumed. Although our quantitative results can be affected by uncertain assumptions in model predictions, our methodology is capable of including models with updated theoretical considerations and additional formation channels. △ Less

Submitted 6 April, 2021; v1 submitted 19 November, 2020; originally announced November 2020.

Comments: 27 pages (14 pages main text + 13 pages appendices/references), 8 figures, 1 table, published in ApJ

Journal ref: The Astrophysical Journal 910, 152 (2021)

The Metallicity Dependence of the High-Mass X-ray Binary Luminosity Function

Authors: Bret D. Lehmer, Rafael T. Eufrasio, Antara Basu-Zych, Keith Doore, Tassos Fragos, Kristen Garofali, Konstantinos Kovlakas, Benjamin J. Williams, Andreas Zezas, Luidhy Santana-Silva

Abstract: We present detailed constraints on the metallicity dependence of the high mass X-ray binary (HMXB) X-ray luminosity function (XLF). We analyze ~5 Ms of Chandra data for 55 actively star-forming galaxies at D < 30 Mpc with gas-phase metallicities spanning 12 + log(O/H) = 7-9.2. Within the galactic footprints, our sample contains a total of 1311 X-ray point sources, of which ~49% are expected to be… ▽ More We present detailed constraints on the metallicity dependence of the high mass X-ray binary (HMXB) X-ray luminosity function (XLF). We analyze ~5 Ms of Chandra data for 55 actively star-forming galaxies at D < 30 Mpc with gas-phase metallicities spanning 12 + log(O/H) = 7-9.2. Within the galactic footprints, our sample contains a total of 1311 X-ray point sources, of which ~49% are expected to be HMXBs, with the remaining sources likely to be low-mass X-ray binaries (LMXBs; ~22%) and unrelated background sources (~29%). We construct a model that successfully characterizes the average HMXB XLF over the full metallicity range. We demonstrate that the SFR-normalized HMXB XLF shows clear trends with metallicity, with steadily increasing numbers of luminous and ultraluminous X-ray sources (logL(erg/s) = 38-40.5) with declining metallicity. However, we find that the low-luminosity (logL(erg/s) = 36-38) HMXB XLF appears to show a nearly constant SFR scaling and slope with metallicity. Our model provides a revised scaling relation of integrated LX/SFR versus 12 + log(O/H) and a new characterization of its the SFR-dependent stochastic scatter. The general trend of this relation is broadly consistent with past studies based on integrated galaxy emission; however, our model suggests that this relation is driven primarily by the high-luminosity end of the HMXB XLF. Our results have implications for binary population synthesis models, the nature of super-Eddington accreting objects (e.g., ultraluminous X-ray sources), recent efforts to identify active galactic nucleus candidates in dwarf galaxies, and the X-ray radiation fields in the early Universe during the epoch of cosmic heating at z > 10. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted for publication in ApJ

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

X-ray Binary Luminosity Function Scaling Relations in Elliptical Galaxies: Evidence for Globular Cluster Seeding of Low-Mass X-ray Binaries in Galactic Fields

Authors: Bret D. Lehmer, Andrew P. Ferrell, Keith Doore, Rafael T. Eufrasio, Erik B. Monson, David M. Alexander, Antara Basu-Zych, William N. Brandt, Greg Sivakoff, Panayiotis Tzanavaris, Mihoko Yukita, Tassos Fragos, Andrew Ptak

Abstract: We investigate X-ray binary (XRB) luminosity function (XLF) scaling relations for Chandra detected populations of low-mass XRBs (LMXBs) within the footprints of 24 early-type galaxies. Our sample includes Chandra and HST observed galaxies at D < 25 Mpc that have estimates of the globular cluster (GC) specific frequency (SN) reported in the literature. As such, we are able to directly classify X-ra… ▽ More We investigate X-ray binary (XRB) luminosity function (XLF) scaling relations for Chandra detected populations of low-mass XRBs (LMXBs) within the footprints of 24 early-type galaxies. Our sample includes Chandra and HST observed galaxies at D < 25 Mpc that have estimates of the globular cluster (GC) specific frequency (SN) reported in the literature. As such, we are able to directly classify X-ray-detected sources as being either coincident with unrelated background/foreground objects, GCs, or sources that are within the fields of the galaxy targets. We model the GC and field LMXB population XLFs for all galaxies separately, and then construct global models characterizing how the LMXB XLFs vary with galaxy stellar mass and SN. We find that our field LMXB XLF models require a component that scales with SN, and has a shape consistent with that found for the GC LMXB XLF. We take this to indicate that GCs are "seeding" the galactic field LMXB population, through the ejection of GC-LMXBs and/or the diffusion of the GCs in the galactic fields themselves. However, we also find that an important LMXB XLF component is required for all galaxies that scales with stellar mass, implying that a substantial population of LMXBs are formed "in situ," which dominates the LMXB population emission for galaxies with SN < 2. For the first time, we provide a framework quantifying how directly-associated GC LMXBs, GC-seeded LMXBs, and in-situ LMXBs contribute to LMXB XLFs in the broader early-type galaxy population. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Comments: Accepted for publication in ApJS. Data products and catalogs available at https://lehmer.uark.edu/downloads/

An alternative interpretation of GW190412 as a binary black hole merger with a rapidly spinning secondary

Authors: Ilya Mandel, Tassos Fragos

Abstract: The LIGO-Virgo collaboration recently reported the properties of GW190412, a binary black hole merger with unequal component masses (mass ratio $0.25^{+0.06}_{-0.04}$ when using the EOBNR PHM approximant) and a non-vanishing effective spin aligned with the orbital angular momentum. They used uninformative priors to infer that the more massive black hole had a dimensionless spin magnitude between 0… ▽ More The LIGO-Virgo collaboration recently reported the properties of GW190412, a binary black hole merger with unequal component masses (mass ratio $0.25^{+0.06}_{-0.04}$ when using the EOBNR PHM approximant) and a non-vanishing effective spin aligned with the orbital angular momentum. They used uninformative priors to infer that the more massive black hole had a dimensionless spin magnitude between 0.17 and 0.59 at 90\% confidence. We argue that, within the context of isolated binary evolution, it is more natural to assume a priori that the first-born, more massive black hole has a negligible spin, while the spin of the less massive black hole is preferentially aligned with the orbital angular momentum if it is spun up by tides. Under this astrophysically motivated prior, we conclude that the lower mass black hole had a dimensionless spin component between 0.64 and 0.99 along the orbital angular momentum. △ Less

Submitted 29 April, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Comments: Updating to version accepted by ApJ Letters

The origin of pulsating ultra-luminous X-ray sources: Low- and intermediate-mass X-ray binaries containing neutron star accretors

Authors: Devina Misra, Tassos Fragos, Thomas Tauris, Emmanouil Zapartas, David R. Aguilera-Dena

Abstract: Ultra-luminous X-ray sources (ULXs) are those X-ray sources located away from the centre of their host galaxy with luminosities exceeding the Eddington limit of a stellar-mass black hole ($L_X>10^{39}\;{\rm erg\,s}^{-1}$). The discovery of X-ray pulsations in some of these objects (e.g. M82~X-2) suggests that a certain fraction of the ULX population may have a neutron star accretor. We present sys… ▽ More Ultra-luminous X-ray sources (ULXs) are those X-ray sources located away from the centre of their host galaxy with luminosities exceeding the Eddington limit of a stellar-mass black hole ($L_X>10^{39}\;{\rm erg\,s}^{-1}$). The discovery of X-ray pulsations in some of these objects (e.g. M82~X-2) suggests that a certain fraction of the ULX population may have a neutron star accretor. We present systematic modelling of low- and intermediate-mass X-ray binaries (LMXBs and IMXBs; donor-star mass range $0.92$--$8.0$~M$_{\odot}$ and neutron-star accretors) to explain the formation of this sub-population of ULXs. Using MESA, we explored the allowed initial parameter space of binary systems consisting of a neutron star and a low- or intermediate-mass donor star that could explain the observed properties of ULXs. Our simulations take into account beaming effects, stellar rotation, general angular momentum losses, and a detailed and self-consistent calculation of the mass-transfer rate. We study the conditions that lead to dynamical stability of these systems, which depends strongly on the response of the donor star to mass loss. Using two values for the initial neutron star mass ($1.3$~M$_{\odot}$ and $2.0$~M$_{\odot}$), we present two sets of mass-transfer calculation grids. We find that LMXBs/IMXBs can produce NS-ULXs with typical time-averaged isotropic-equivalent X-ray luminosities of $10^{39}$--$10^{41}\;{\rm erg\,s}^{-1}$ on a timescale up to $\sim\!1.0\;{\rm Myr}$ for the lower luminosities. We also estimate their likelihood of detection, the types of white-dwarf remnants left behind by the donors, and the total amount of mass accreted by the neutron stars. We also compare our results to the observed pulsating ULXs. Our results suggest that a large subset of the observed pulsating ULX population can be explained by LMXBs/IMXBs undergoing a super-Eddington mass-transfer phase. △ Less

Submitted 25 August, 2020; v1 submitted 2 April, 2020; originally announced April 2020.

Comments: 19 pages, 13 figures, Accepted by A&A. Parameter space was increased to include low-mass XRBs and corresponding changes made to the text (including the title) and figures 4, 6-11. Changed axes for figures 1 and 2. Fixed typos and updated references. Added arguments about why spin period is not an accurate reflection of mass accretion rate in the introduction

Journal ref: A&A 642, A174 (2020)

New insight on the far-UV SED and HeII emission from low metallicity galaxies

Authors: D. Schaerer, Y. Izotov, T. Fragos

Abstract: Understanding the ionizing spectrum of low-metallicity galaxies is of great importance for modeling and interpreting emission line observations of early/distant galaxies. Although a wide suite of stellar evolution, atmosphere, population synthesis, and photoionization models, taking many physical processes into account now exist, all models face a common problem: the inability to explain the pre… ▽ More Understanding the ionizing spectrum of low-metallicity galaxies is of great importance for modeling and interpreting emission line observations of early/distant galaxies. Although a wide suite of stellar evolution, atmosphere, population synthesis, and photoionization models, taking many physical processes into account now exist, all models face a common problem: the inability to explain the presence of nebular HeII emission, which is observed in many low metallicity galaxies, both in UV and optical spectra. Several possible explanations have been proposed in the literature, including Wolf-Rayet (WR) stars, binaries, very massive stars, X-ray sources, or shocks. However, none has so far been able to explain the major observations. We briefly discuss the HeII problem, available empirical data, and observed trends combining X-ray, optical and other studies. We present a simple and consistent physical model showing that X-ray binaries could explain the long-standing nebular HeII problem. Our model, described in Schaerer et al. (2019), successfully explains the observed trends and strength of nebular HeII emission in large samples of low metallicity galaxies and in individual galaxies, which have been studied in detail and with multi-wavelength observations. Our results have in particular important implications for the interpretation of galaxy spectra in the early Universe, which will be obtained with upcoming and future facilities. △ Less

Submitted 22 November, 2019; originally announced November 2019.

Comments: 5 pages, 2 figures. Proceedings of IAU Symposium No. 352, 2019 "Uncovering early galaxy evolution in the ALMA and JWST era". Eds. E. da Cunha, J. Hodge, J. Afonso, L. Pentericci

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

The origin of spin in binary black holes: Predicting the distributions of the main observables of Advanced LIGO

Authors: Simone S. Bavera, Tassos Fragos, Ying Qin, Emmanouil Zapartas, Coenraad J. Neijssel, Ilya Mandel, Aldo Batta, Sebastian M. Gaebel, Chase Kimball, Simon Stevenson

Abstract: We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of the main observables of Advanced LIGO. We used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code.… ▽ More We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of the main observables of Advanced LIGO. We used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. By assuming efficient angular momentum transport, we are able to present a model that is capable of simultaneously predicting the three main gravitational-wave observables: the effective inspiral spin parameter $χ_{eff}$, the chirp mass $M_{chirp}$, and the cosmological redshift of merger $z_{merger}$. We find an excellent agreement between our model and the ten events from the first two advanced detector observing runs. We make predictions for the third observing run O3 and for Advanced LIGO design sensitivity. We expect around 80% of events with $χ_{eff} < 0.1$, while the remaining 20% of events with $χ_{eff} \ge 0.1$ are split into ~10% with $M_{chirp} < 15$ M$_\odot$ and ~10% with $M_{chirp} \ge 15$ M$_\odot$. In conclusion, the favorable comparison of the existing LIGO/Virgo observations with our model predictions gives support to the idea that the majority, if not all of the observed mergers, originate from the evolution of isolated binaries. The first-born black hole has negligible spin because it lost its envelope after it expanded to become a giant star, while the spin of the second-born black hole is determined by the tidal spin up of its naked helium star progenitor by the first-born black hole companion after the binary finished the common-envelope phase. △ Less

Submitted 24 January, 2020; v1 submitted 28 June, 2019; originally announced June 2019.

Comments: Accepted for publication in A&A on January 21th 2020, 19 pages, 9 figures

Journal ref: A&A 635, A97 (2020)

X-ray Binary Luminosity Function Scaling Relations for Local Galaxies Based on Subgalactic Modeling

Authors: Bret D. Lehmer, Rafael T. Eufrasio, Panayiotis Tzanavaris, Antara Basu-Zych, Tassos Fragos, Andrea Prestwich, Mihoko Yukita, Andreas Zezas, Ann E. Hornschemeier, Andrew Ptak

Abstract: We present new Chandra constraints on the X-ray luminosity functions (XLFs) of X-ray binary (XRB) populations, and their scaling relations, for a sample of 38 nearby galaxies (D = 3.4-29 Mpc). Our galaxy sample is drawn primarily from the Spitzer infrared nearby galaxy survey (SINGS), and contains a wealth of Chandra (5.8 Ms total) and multiwavelength data, allowing for star-formation rates (SFRs)… ▽ More We present new Chandra constraints on the X-ray luminosity functions (XLFs) of X-ray binary (XRB) populations, and their scaling relations, for a sample of 38 nearby galaxies (D = 3.4-29 Mpc). Our galaxy sample is drawn primarily from the Spitzer infrared nearby galaxy survey (SINGS), and contains a wealth of Chandra (5.8 Ms total) and multiwavelength data, allowing for star-formation rates (SFRs) and stellar masses (M*) to be measured on subgalactic scales. We divided the 2478 X-ray detected sources into 21 subsamples in bins of specific-SFR (sSFR = SFR/M*) and constructed XLFs. To model the XLF dependence on sSFR, we fit a global XLF model, containing contributions from high-mass XRBs (HMXBs), low-mass XRBs (LMXBs), and background sources from the cosmic X-ray background (CXB) that respectively scale with SFR, M*, and sky area. We find an HMXB XLF that is more complex in shape than previously reported and an LMXB XLF that likely varies with sSFR, potentially due to an age dependence. When applying our global model to XLF data for each individual galaxy, we discover a few galaxy XLFs that significantly deviate from our model beyond statistical scatter. Most notably, relatively low-metallicity galaxies have an excess of HMXBs above ~10^38 erg/s and elliptical galaxies that have relatively rich populations of globular clusters (GCs) show excesses of LMXBs compared to the global model. Additional modeling of how the XRB XLF depends on stellar age, metallicity, and GC specific frequency is required to sufficiently characterize the XLFs of galaxies. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Comments: Accepted for publication in ApJS (37 pages, 16 figures, and 7 tables). Data products available (https://lehmer.uark.edu/downloads/)

Stellar models and isochrones from low-mass to massive stars including pre-main sequence phase with accretion

Authors: L. Haemmerlé, P. Eggenberger, S. Ekström, C. Georgy, G. Meynet, A. Postel, M. Audard, M. Sørensen, T. Fragos

Abstract: Grids of stellar models are useful tools to derive the properties of stellar clusters, in particular young clusters hosting massive stars, and to provide information on the star formation process in various mass ranges. Because of their short evolutionary timescale, massive stars end their life while their low-mass siblings are still on the pre-main sequence (pre-MS) phase. Thus the study of young… ▽ More Grids of stellar models are useful tools to derive the properties of stellar clusters, in particular young clusters hosting massive stars, and to provide information on the star formation process in various mass ranges. Because of their short evolutionary timescale, massive stars end their life while their low-mass siblings are still on the pre-main sequence (pre-MS) phase. Thus the study of young clusters requires consistent consideration of all the phases of stellar evolution. But despite the large number of grids that are available in the literature, a grid accounting for the evolution from the pre-MS accretion phase to the post-MS phase in the whole stellar mass range is still lacking. We build a grid of stellar models at solar metallicity with masses from 0.8 $M_\odot$ to 120 $M_\odot$, including pre-MS phase with accretion. We use the {\sc genec} code to run stellar models on this mass range. The accretion law is chosen to match the observations of pre-MS objects on the Hertzsprung-Russell diagram. We describe the evolutionary tracks and isochrones of our models. The grid is connected to previous MS and post-MS grids computed with the same numerical method and physical assumptions, which provides the widest grid in mass and age to date. Numerical tables of our models and corresponding isochrones are available online. △ Less

Submitted 25 March, 2019; originally announced March 2019.

Journal ref: A&A 624, A137 (2019)

Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s

Authors: R. J. Foley, K. D. Alexander, I. Andreoni, I. Arcavi, K. Auchettl, J. Barnes, G. Baym, E. C. Bellm, A. M. Beloborodov, N. Blagorodnova, J. P. Blakeslee, P. R. Brady, M. Branchesi, J. S. Brown, N. Butler, M. Cantiello, R. Chornock, D. O. Cook, J. Cooke, D. L. Coppejans, A. Corsi, S. M. Couch, M. W. Coughlin, D. A. Coulter, P. S. Cowperthwaite , et al. (88 additional authors not shown)

Abstract: As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse s… ▽ More As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submitted to Astro2020

X-ray binaries as the origin of nebular HeII emission in low-metallicity star-forming galaxies

Authors: D. Schaerer, T. Fragos, Y. Izotov

Abstract: The origin of nebular HeII emission, which is frequently observed in low-metallicity (O/H) star-forming galaxies, remains largely an unsolved question. Using the observed anticorrelation of the integrated X-ray luminosity per unit of star formation rate ($L_X/{\rm SFR}$) of an X-ray binary population with metallicity and other empirical data from the well-studied galaxy I Zw 18, we show that the o… ▽ More The origin of nebular HeII emission, which is frequently observed in low-metallicity (O/H) star-forming galaxies, remains largely an unsolved question. Using the observed anticorrelation of the integrated X-ray luminosity per unit of star formation rate ($L_X/{\rm SFR}$) of an X-ray binary population with metallicity and other empirical data from the well-studied galaxy I Zw 18, we show that the observed HeII 4686 intensity and its trend with metallicity is naturally reproduced if the bulk of He$^+$ ionizing photons are emitted by the X-ray sources. We also show that a combination of X-ray binary population models with normal single and/or binary stellar models reproduces the observed $I(4686)/I(Hβ)$ intensities and its dependency on metallicity and age. We conclude that both empirical data and theoretical models suggest that high-mass X-ray binaries are the main source of nebular HeII emission in low-metallicity star-forming galaxies. △ Less

Submitted 27 February, 2019; originally announced February 2019.

Comments: 5 pages, 3 figures. Accepted for publication in A&A Letters

Journal ref: A&A 622, L10 (2019)

On the Origin of Black Hole Spin in High-mass X-Ray Binaries

Authors: Ying Qin, Pablo Marchant, Tassos Fragos, Georges Meynet, Vicky Kalogera

Abstract: Black hole (BH) spins in low-mass X-ray binaries (LMXBs) cover a range of values that can be explained by accretion after BH birth. In contrast, the three BH spin measurements in high-mass X-ray binaries (HMXBs) show only values near the maximum and likely have a different origin connected to the BH stellar progenitor. We explore here two possible scenarios to explain the high spins of BHs in HMXB… ▽ More Black hole (BH) spins in low-mass X-ray binaries (LMXBs) cover a range of values that can be explained by accretion after BH birth. In contrast, the three BH spin measurements in high-mass X-ray binaries (HMXBs) show only values near the maximum and likely have a different origin connected to the BH stellar progenitor. We explore here two possible scenarios to explain the high spins of BHs in HMXBs: formation in binaries that undergo mass transfer (MT) during the main sequence (MS; Case-A MT), and very close binaries undergoing chemically homogeneous evolution (CHE). We find that both scenarios are able to produce high-spin BHs if internal angular momentum (AM) transport in the progenitor star after its MS evolution is not too strong (i.e., weak coupling between the stellar core and its envelope). If instead efficient AM transport is assumed, we find that the resulting BH spins are always too low with respect to observations. The Case-A MT model provides a good fit for the BH spins, the masses of the two components, and the final orbital periods for two of the three BHs in HMXBs with measured spins. For one of them, the mass predicted for the BH companion is significantly lower than observed, but this depends strongly on the assumed efficiency of mass transfer. The CHE models predict orbital periods that are too large for all three cases considered here. We expect the Case-A MT to be much more frequent at the metallicities involved, so we conclude that the Case-A MT scenario is preferred. Finally, we predict that the stellar companions of HMXBs formed through the Case-A MT have enhanced nitrogen surface abundances, which can be tested by future observations. △ Less

Submitted 20 January, 2019; v1 submitted 30 October, 2018; originally announced October 2018.

Comments: 8 pages, 5 figures and 1 table, 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

Formation of massive close binaries I. Detached evolution

Authors: Mads Sørensen, Tassos Fragos, Georges Meynet, Lionel Haemmerlé

Abstract: We explore the initial conditions that lead to detached binary evolution of massive pre-MS binaries and ask how large a fraction of the observed binary systems may have been initially formed as low-mass protobinaries and later undergone a significant accretion phase while remaining detached. We develop a family of analytic models to describe the orbital separation, $a$, and mass ratio, $q$, evolut… ▽ More We explore the initial conditions that lead to detached binary evolution of massive pre-MS binaries and ask how large a fraction of the observed binary systems may have been initially formed as low-mass protobinaries and later undergone a significant accretion phase while remaining detached. We develop a family of analytic models to describe the orbital separation, $a$, and mass ratio, $q$, evolution. For a given mass accretion rate onto the binary system, we define a recipe for distributing this mass between the two components. For this we introduce a parameter $η$, such that $\frac{\dot{M}_2}{\dot{M}_1}=q^η$ at any time, to determine the binary mass ratio evolution. Depending on the choice of $η$, any type of mass ratio evolution is possible. Furthermore, we use MESA, a detailed stellar structure code, to calculate an extensive grid of binary sequences where a protobinary undergoes accretion, and we identify the initial conditions that separates detached from non-detached pre-MS binary evolution. A value of $η$ around 2 allows accretion growth in detached systems to form close massive binaries on the Zero Age Main-Sequence with minimum orbital periods down to about 1.2 days for $M_{\rm 1, ZAMS}=20-30M_{\odot}$ twin-binaries. $η=2$ can also reproduce the observed population of binary systems with primary stars above $6M_{\odot}$. The whole observed range of massive close binaries can form via accretion growth in detached systems, making the binary formation channel of accretion growth a strong contender to explain the formation of massive close binaries, including progenitors of coalescing binary black holes. △ Less

Submitted 20 August, 2018; originally announced August 2018.

Comments: 27 pages, 19 figures, 3 tables, submitted to A&A. Abstract context removed to to fit arXiv

The spin of the second-born black hole in coalescing binary black holes

Authors: Y. Qin, T. Fragos, G. Meynet, J. Andrews, M. Sørensen, H. F. Song

Abstract: Various binary black hole formation channels have been proposed since the first gravitational event GW150914 was discovered by the Advanced Laser Interferometer Gravitational-Wave Observatory (AdLIGO). For all evolutionary channels based on the evolution of isolated binaries, the immediate progenitor of the binary black hole is a close binary system composed of a black hole and a helium star. We p… ▽ More Various binary black hole formation channels have been proposed since the first gravitational event GW150914 was discovered by the Advanced Laser Interferometer Gravitational-Wave Observatory (AdLIGO). For all evolutionary channels based on the evolution of isolated binaries, the immediate progenitor of the binary black hole is a close binary system composed of a black hole and a helium star. We perform detailed binary evolution and systematically explore the parameter space of initial binary properties, including initial black hole and helium star masses, initial rotation of the helium star as well as metallicity. We argue that the spin of the first-born black hole at its birth is negligible ($\lesssim 0.1$), hence the second-born black hole's spin dominates the measured effective spin, $χ_{\rm eff}$, from gravitational wave events of double black hole mergers. We find that tides can be important only when orbital periods are shorter than 2 days. Upon core collapse, the helium star produces a black hole (the second-born black hole in the system) with a spin that can span the entire range from zero to maximally spinning. We show that the bimodal distribution of the spin of the second-born black hole obtained in recent papers is mainly due to oversimplifying assumptions. We find an anti-correlation between the merging timescale of the two black holes, T$_{\rm merger}$, and the effective spin $χ_{\rm eff}$. Finally, we provide new prescriptions for the tidal coefficient E$_2$ for both H-rich and the helium-rich stars. We predict that, with future improvements to AdLIGO's sensitivity, the sample of merging binary black hole systems will show an overdensity of sources with positive but small $χ_{\rm eff}$. originating from lower mass black hole mergers born at low redshift. △ Less

Submitted 15 February, 2018; originally announced February 2018.

Comments: 20 pages, 15 figures; submitted to A&A

Journal ref: A&A 616, A28 (2018)

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

On the Spatially Resolved Star-Formation History in M51 II: X-ray Binary Population Evolution

Authors: B. D. Lehmer, R. T Eufrasio, L. Markwardt, A. Zezas, A. Basu-Zych, T. Fragos, A. E. Hornschemeier, A. Ptak, P. Tzanavaris, M. Yukita

Abstract: We present a new technique for empirically calibrating how the X-ray luminosity function (XLF) of X-ray binary (XRB) populations evolves following a star-formation event. We first utilize detailed stellar population synthesis modeling of far-UV to far-IR photometry of the nearby face-on spiral galaxy M51 to construct maps of the star-formation histories (SFHs) on subgalactic (~400 pc) scales. Next… ▽ More We present a new technique for empirically calibrating how the X-ray luminosity function (XLF) of X-ray binary (XRB) populations evolves following a star-formation event. We first utilize detailed stellar population synthesis modeling of far-UV to far-IR photometry of the nearby face-on spiral galaxy M51 to construct maps of the star-formation histories (SFHs) on subgalactic (~400 pc) scales. Next, we use the ~850 ks cumulative Chandra exposure of M51 to identify and isolate 2-7 keV detected point sources within the galaxy, and we use our SFH maps to recover the local properties of the stellar populations in which each X-ray source is located. We then divide the galaxy into various subregions based on their SFH properties (e.g., star-formation rate [SFR] per stellar mass [M*] and mass-weighted stellar age) and group the X-ray point sources according to the characteristics of the regions in which they are found. Finally, we construct and fit a parameterized XLF model that quantifies how the XLF shape and normalization evolves as a function of the XRB population age. Our best-fit model indicates the XRB XLF per unit stellar mass declines in normalization, by ~3-3.5 dex, and steepens in slope from ~10 Myr to ~10 Gyr. We find that our technique recovers results from past studies of how XRB XLFs and XRB luminosity scaling relations vary with age and provides a self-consistent picture for how the XRB XLF evolves with age. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ (19 pages, 15 figures, and 2 tables)

On the Spatially Resolved Star Formation History in M51 I: Hybrid UV+IR Star Formation Laws and IR Emission from Dust Heated by Old Stars

Authors: R. T. Eufrasio, B. D. Lehmer, A. Zezas, E. Dwek, R. G. Arendt, A. Basu-Zych, T. Wiklind, M. Yukita, T. Fragos, A. E. Hornschemeier, L. Markwardt, A. Ptak, P. Tzanavaris

Abstract: We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determin… ▽ More We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determine SFH step intensities and combines this with a grid-based approach to determine three extinction parameters. We apply our procedure to the extensive FUV-to-FIR photometric data of M51, convolved to a common spatial resolution and pixel scale, and make the resulting maps publicly available. We recover, for M51a, a peak star formation rate (SFR) between 0.1 and 5 Gyr ago, with much lower star formation activity over the last 100 Myr. For M51b, we find a declining SFR toward the present day. In the outskirt regions of M51a, which includes regions between M51a and M51b, we recover a SFR peak between 0.1 and 1 Gyr ago, which corresponds to the effects of the interaction between M51a and M51b. We utilize our results to (1) illustrate how UV+IR hybrid SFR laws vary across M51, and (2) provide first-order estimates for how the IR luminosity per unit stellar mass varies as a function of the stellar age. From the latter result, we find that IR emission from dust heated by stars is not always associated with young stars, and that the IR emission from M51b is primarily powered by stars older than 5 Gyr. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ

Close binary evolution. III. Impact of tides, wind magnetic braking, and internal angular momentum transport

Authors: H. F. Song, G. Meynet, A. Maeder, S. Ekstrom, P. Eggenberger, C. Georgy, Y. Qin, T. Fragos, M. Soerensen, F. Barblan, G. A. Wade

Abstract: Massive stars with solar metallicity lose important amounts of rotational angular momentum through their winds. When a magnetic field is present at the surface of a star, efficient angular momentum losses can still be achieved even when the mass-loss rate is very modest, at lower metallicities, or for lower-initial-mass stars. In a close binary system, the effect of wind magnetic braking also inte… ▽ More Massive stars with solar metallicity lose important amounts of rotational angular momentum through their winds. When a magnetic field is present at the surface of a star, efficient angular momentum losses can still be achieved even when the mass-loss rate is very modest, at lower metallicities, or for lower-initial-mass stars. In a close binary system, the effect of wind magnetic braking also interacts with the influence of tides, resulting in a complex evolution of rotation. We study the interactions between the process of wind magnetic braking and tides in close binary systems. We discuss the evolution of a 10 M$_\odot$ star in a close binary system with a 7 M$_\odot$ companion using the Geneva stellar evolution code. The initial orbital period is 1.2 days. The 10 M$_\odot$ star has a surface magnetic field of 1 kG. Various initial rotations are considered. We use two different approaches for the internal angular momentum transport. In one of them, angular momentum is transported by shear and meridional currents. In the other, a strong internal magnetic field imposes nearly perfect solid-body rotation. The evolution of the primary is computed until the first mass-transfer episode occurs. The cases of different values for the magnetic fields and for various orbital periods and mass ratios are briefly discussed. We show that, independently of the initial rotation rate of the primary and the efficiency of the internal angular momentum transport, the surface rotation of the primary will converge, in a time that is short with respect to the main-sequence lifetime, towards a slowly evolving velocity that is different from the synchronization velocity. (abridged). △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: 11 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 609, A3 (2018)

Radiation Backgrounds at Cosmic Dawn: X-Rays from Compact Binaries

Authors: Piero Madau, Tassos Fragos

Abstract: We compute the expected X-ray diffuse background and radiative feedback on the intergalactic medium (IGM) from X-ray binaries prior and during the epoch of reionization. The cosmic evolution of compact binaries is followed using a population synthesis technique that treats separately neutron stars and black hole binaries in different spectral states and is calibrated to reproduce the observed X-ra… ▽ More We compute the expected X-ray diffuse background and radiative feedback on the intergalactic medium (IGM) from X-ray binaries prior and during the epoch of reionization. The cosmic evolution of compact binaries is followed using a population synthesis technique that treats separately neutron stars and black hole binaries in different spectral states and is calibrated to reproduce the observed X-ray properties of galaxies at z<4. Together with an updated empirical determination of the cosmic history of star formation, recent modeling of the stellar mass-metallicity relation, and a scheme for absorption by the IGM that accounts for the presence of ionized HII bubbles during the epoch of reionization, our detailed calculations provide refined predictions of the X-ray volume emissivity and filtered radiation background from "normal" galaxies at z>6. Radiative transfer effects modulate the background spectrum, which shows a characteristic peak between 1 and 2 keV. While the filtering of X-ray radiation through the IGM slightly increases the mean excess energy per photoionization, it also weakens the radiation intensity below 1 keV, lowering the mean photoionization and heating rates. Numerical integration of the rate and energy equations shows that the contribution of X-ray binaries to the ionization of the bulk IGM is negligible, with the electron fraction never exceeding 1%. Direct HeI photoionizations are the main source of IGM heating, and the temperature of the largely neutral medium in between HII cavities increases above the temperature of the cosmic microwave background (CMB) only at z<10, when the volume filling factor of HII bubbles is already >0.1. Therefore, in this scenario, it is only at relatively late epochs that the bulk of neutral intergalactic hydrogen may be observable in 21-cm emission against the CMB. △ Less

Submitted 17 April, 2017; v1 submitted 25 June, 2016; originally announced June 2016.

Comments: 14 pages, 8 figures, accepted for publication in The Astrophysical Journal