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.

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.

Ultra-luminous X-ray sources (ULXs) are X-ray sources located off-galactic centre and have luminosities exceeding the Eddington limit of a stellar-mass black hole ($L_X>10^{39}\;{\rm erg\,s}^{-1}$). With the discovery of X-ray pulsations in some of these objects (e.g. M82 X-2), we know a part of the ULX population has a neutron star (NS) accretor. We present systematic modelling of intermediate-mass X-ray binaries (IMXBs; donor-star mass range $2.0$ to $8.0\;$M$_{\odot}$) to explain the formation of this sub-population of ULXs. Using MESA, we explore the allowed initial parameter space of binary systems consisting of a neutron star and an intermediate-mass donor star, that could explain the observed properties of ULXs. Our simulations take into account beaming effects and also include stellar rotation, tides, general angular momentum losses, and a detailed and self-consistent calculation of the mass-transfer rate. Exploring the initial parameters that lead to the formation of neutron-star ULXs, 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 NS mass ($1.3\;$M$_{\odot}$ and $2.0\;$M$_{\odot}$) we find that IMXBs can produce NS-ULXs with typical time-averaged isotropic-equivalent X-ray luminosities between $10^{39}$--$10^{41}\;{\rm erg\,s}^{-1}$ on a timescale of up to $\sim\!1.0\;$Myr for the lower luminosities, assuming geometrical beaming. Finally, we estimate their likelihood of detection, the types of white-dwarf remnants left behind by the donors, as well as the total amount of mass accreted by the neutron stars. We also compare our results to the observed pulsating ULXs and infer their initial parameters. Our results suggest that a large subset of the observed pulsating ULX population can be explained by IMXBs in a super-Eddington mass-transfer phase.

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 $\chi_{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 $\chi_{eff} < 0.1$, while the remaining 20% of events with $\chi_{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.

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.

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 $\alpha_{\rm CE}$-equivalent efficiency of $\approx$5.

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.

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.

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.

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\beta)$ 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.

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.

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.

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 $\eta$, such that $\frac{\dot{M}_2}{\dot{M}_1}=q^\eta$ at any time, to determine the binary mass ratio evolution. Depending on the choice of $\eta$, 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 $\eta$ 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. $\eta=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.

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.

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, $\chi_{\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 $\chi_{\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 $\chi_{\rm eff}$. originating from lower mass black hole mergers born at low redshift.

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.

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.

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).

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.

We have endeavoured to understand the formation and evolution of the black hole (BH) X-ray binary LMC X-3. We estimate the properties of the system at 4 evolutionary stages: 1) at the Zero Age Main Sequence (ZAMS), 2) just prior to the supernova (SN) explosion of the primary, 3) just after the SN, and 4) at the moment of RLO onset.We use a hybrid approach, combining detailed stellar structure and binary evolution calculations with approximate population synthesis models. This allows us to estimate potential natal kicks and the evolution of the BH spin. In the whole analysis we incorporate as model constraints the most up-to-date observational information, encompassing the binary's orbital properties, the companion star mass, effective temperature, surface gravity and radius, as well as the black hole's mass and spin. We find that LMC X-3 began as a ZAMS system with the mass of the primary star in the range $M_{\rm{1,ZAMS}}$ = 22-31 $\rm{M_{\odot}}$ and a secondary star of $M_{\rm{2,ZAMS}} = 5.0-8.3M_{\odot}$, in a wide ($P_{ZAMS} \gtrsim 2.000\, \rm days$) and eccentric ($e_{\rm{ZAMS}} \gtrsim 0.23$) orbit. Just prior to the SN, the primary has a mass of $M_{\mathrm{1,preSN}} = 11.1-18.0\,\rm M_{\odot}$, with the secondary star largely unaffected. The orbital period decreases to $0.6-1.7\, \rm days$, and is still eccentric $0 \leq e_{\rm{preSN}} \leq 0.44$. We find that a symmetric SN explosion with no or small natal kicks (a few tens of $\rm km\, s^{-1}$) imparted on the BH cannot be formally excluded, however, large natal kicks in excess of $\gtrsim 120 \,\rm km\, s^{-1}$ increase the estimated formation rate by an order of magnitude. Following the SN, the system has a BH $M_{\mathrm{BH,postSN}} = 6.4-8.2\,\rm M_{\odot}$ and is put into an eccentric orbit. At the RLO onset the orbit is circularised and it has an orbital period of $P_{\rm{RLO}} = 0.8-1.4\,\rm days$.

We present measurements of the evolution of normal-galaxy X-ray emission from $z \approx$ 0-7 using local galaxies and galaxy samples in the 6 Ms Chandra Deep Field-South (CDF-S) survey. The majority of the CDF-S galaxies are observed at rest-frame energies above 2 keV, where the emission is expected to be dominated by X-ray binary (XRB) populations; however, hot gas is expected to provide small contributions to the observed- frame < 1 keV emission at $z < 1$. We show that a single scaling relation between X-ray luminosity ($L_{\rm X}$) and star-formation rate (SFR) is insufficient for characterizing the average X-ray emission at all redshifts. We establish that scaling relations involving not only SFR, but also stellar mass ($M_\star$) and redshift, provide significantly improved characterizations of the average X-ray emission from normal galaxy populations at $z \approx$ 0-7. We further provide the first empirical constraints on the redshift evolution of X-ray emission from both low-mass XRB (LMXB) and high-mass XRB (HMXB) populations and their scalings with $M_\star$ and SFR, respectively. We find $L_{\rm 2-10 keV}({\rm LMXB})/M_\star \propto (1 + z)^{2-3}$ and $L_{\rm 2-10 keV}({\rm HMXB})$/SFR $\propto (1 + z)$, and show that these relations are consistent with XRB population-synthesis model predictions, which attribute the increase in LMXB and HMXB scaling relations with redshift as being due to declining host galaxy stellar ages and metallicities, respectively. We discuss how emission from XRBs could provide an important source of heating to the intergalactic medium in the early Universe, exceeding that of active galactic nuclei.

We have studied high mass X-ray binary (HMXB) populations within two low-metallicity, starburst galaxies, Haro 11 and VV 114. These galaxies serve as analogs to high-redshift (z>2) Lyman break galaxies, and within the larger sample of Lyman break analogs (LBAs) are sufficiently nearby (<87 Mpc) to be spatially-resolved by Chandra. Previous studies of the X-ray emission in LBAs have found that the 2-10 keV luminosity per star formation rate (SFR) in these galaxies is elevated, potentially because of their low metallicities (12+log[O/H]= 8.3-8.4). Theoretically, the progenitors of XRBs forming in lower metallicity environments lose less mass from stellar winds over their lifetimes, producing more massive compact objects (i.e., neutron stars and black holes), and thus resulting in more numerous and luminous HMXBs per SFR. We have performed an in-depth study of the only two LBAs that have spatially-resolved 2-10 keV emission with Chandra to present the bright end of the X-ray luminosity distribution of HMXBs (L$_X>10^{39}$ erg/s; ultraluminous X-ray sources, ULXs) in these low-metallicity galaxies, based on 8 detected ULXs. Comparing with the star-forming galaxy X-ray luminosity function (XLF) presented by Mineo et al. (2012), Haro 11 and VV 114 host ~4 times more L$_X>10^{40}$ erg/s sources than expected given their SFRs. We simulate the effects of source blending from crowded lower luminosity HMXBs using the star-forming galaxy XLF and then vary the XLF shapes until we reproduce the observed point source luminosity distributions. We find that these LBAs have a shallower bright end slope than the standard XLF. If we conservatively assume that the brightest X-ray source from each galaxy is powered by an AGN rather than a HMXB and eliminate these sources from consideration, the luminosity distribution becomes poorly constrained but does appear to be consistent with a standard XLF. [Abridged]

The formation of Black Hole (BH) Low-Mass X-ray Binaries (LMXB) poses a theoretical challenge, as low-mass companions are not expected to survive the common-envelope scenario with the BH progenitor. Here we propose a formation mechanism that skips the common-envelope scenario and relies on triple-body dynamics. We study the evolution of hierarchical triples, following the secular dynamical evolution up to the octupole-level of approximation, including general relativity, tidal effects and post-main-sequence evolution, such as mass loss, changes to stellar radii and supernovae. During the dynamical evolution of the triple system, the "eccentric Kozai-Lidov" mechanism can cause large eccentricity excitations in the LMXB progenitor, resulting in three main BH-LMXB formation channels. Here we define BH-LMXB candidates as systems where the inner BH companion star crosses its Roche limit. In the "eccentric" channel (~ 81% of the LMXBs in our simulations), the donor star crosses its Roche limit during an extreme eccentricity excitation, while still on a wide orbit. Second, we find a "giant" LMXB channel (~ 11%), where a system undergoes only moderate eccentricity excitations, but the donor star fills its Roche lobe after evolving toward the giant branch. Third, we identify a "classical" channel (~8%), where tidal forces and magnetic braking shrink and circularize the orbit to short periods, triggering mass transfer. Finally, for the giant channel, we predict an eccentric ($\sim 0.3-0.6$), preferably inclined (~40, ~140 degreed) tertiary, typically on a wide enough orbit (~10^4AU), to potentially become unbound later in the triple evolution.

The recent discovery of a neutron star accretor in the ultra-luminous X-ray source M82 X-2 challenges our understanding of high-mass X-ray binary formation and evolution. By combining binary population synthesis and detailed mass-transfer models, however, we show that the binary parameters of M82 X-2 are not surprising provided non-conservative mass transfer is allowed. Specifically, the donor-mass lower limit and orbital period measured for M82 X-2 lie near the most probable values predicted by population synthesis models, and systems such as M82 X-2 should exist in approximately 13% of the galaxies with a star-formation history similar to M82. We conclude that the binary system that formed M82 X-2 is most likely less than 50 Myr old and contains a donor star which had an initial mass of approximately 8-10 M$_\odot$, while the NS's progenitor star had an initial mass in the $8-25\,\rm M_{\odot}$ range. The donor star still currently resides on the main sequence, and is capable of continued MT on the thermal timescale, while in the ultra-luminous X-ray regime, for as long as 400,000 years.

Galactic field black hole (BH) low-mass X-ray binaries (LMXBs) are believed to form in situ via the evolution of isolated binaries. In the standard formation channel, these systems survived a common envelope phase, after which the remaining helium core of the primary star and the subsequently formed BH are not expected to be highly spinning. However, the measured spins of BHs in LMXBs cover the whole range of spin parameters. We propose here that the BH spin in LMXBs is acquired through accretion onto the BH after its formation. In order to test this hypothesis, we calculated extensive grids of detailed binary mass-transfer sequences. For each sequence, we examined whether, at any point in time, the calculated binary properties are in agreement with their observationally inferred counterparts of 16 Galactic LMXBs. The "successful" sequences give estimates of the mass that the BH has accreted since the onset of Roche-Lobe overflow. We find that in all Galactic LMXBs with measured BH spin, the origin of the spin can be accounted for by the accreted matter, and we make predictions about the maximum BH spin in LMXBs where no measurement is yet available. Furthermore, we derive limits on the maximum spin that any BH can have depending on current properties of the binary it resides in. Finally we discuss the implication that our findings have on the BH birth-mass distribution, which is shifted by ~1.5 Msolar towards lower masses, compared to the currently observed one.

We present direct constraints on how the formation of low-mass X-ray binary (LMXB) populations in galactic fields depends on stellar age. In this pilot study, we utilize Chandra and Hubble Space Telescope (HST) data to detect and characterize the X-ray point source populations of three nearby early-type galaxies: NGC 3115, 3379, and 3384. The luminosity-weighted stellar ages of our sample span 3-10 Gyr. X-ray binary population synthesis models predict that the field LMXBs associated with younger stellar populations should be more numerous and luminous per unit stellar mass than older populations due to the evolution of LMXB donor star masses. Crucially, the combination of deep Chandra and HST observations allows us to test directly this prediction by identifying and removing counterparts to X-ray point sources that are unrelated to the field LMXB populations, including LMXBs that are formed dynamically in globular clusters, Galactic stars, and background AGN/galaxies. We find that the "young" early-type galaxy NGC 3384 (~2-5 Gyr) has an excess of luminous field LMXBs (L_X > (5-10) x 10^37 erg/s) per unit K-band luminosity (L_K; a proxy for stellar mass) than the "old" early-type galaxies NGC 3115 and 3379 (~8-10 Gyr), which results in a factor of ~2-3 excess of LX/LK for NGC 3384. This result is consistent with the X-ray binary population synthesis model predictions; however, our small galaxy sample size does not allow us to draw definitive conclusions on the evolution field LMXBs in general. We discuss how future surveys of larger galaxy samples that combine deep Chandra and HST data could provide a powerful new benchmark for calibrating X-ray binary population synthesis models.

The extragalactic X-ray binary IC 10 X-1 has attracted attention as it is possibly the host of the most massive stellar-mass black-hole (BH) known to date. Here we consider all available observational constraints and construct its evolutionary history up to the instant just before the formation of the BH. Our analysis accounts for the simplest possible history that includes three evolutionary phases: binary orbital dynamics at core collapse, common envelope (CE) evolution, and evolution of the BH--helium star binary progenitor of the observed system. We derive the complete set of constraints on the progenitor system at various evolutionary stages. Specifically: right before the core collapse event, we find the mass of the BH immediate progenitor to be > 31 Msun (at 95% of confidence, same hereafter). The magnitude of the natal kick imparted to the BH is constrained to be < 130 km/s. Furthermore, we find that the "enthalpy" formalism recently suggested by Ivanova & Chaichenets is able to explain the existence of IC 10 X-1 without the need of invoking unreasonably high CE efficiencies. With this physically motivated formalism, we find that the CE efficiency required to explain the system is in the range of 0.6--1.

We report significant anisotropies in the projected two-dimensional (2D) spatial distributions of Globular Clusters (GCs) of the giant Virgo elliptical galaxy NGC4649 (M60). Similar features are found in the 2D distribution of low-mass X-ray binaries (LMXBs), both associated with GCs and in the stellar field. Deviations from azimuthal symmetry suggest an arc-like excess of GCs extending north at 4-15 kpc galactocentric radii in the eastern side of major axis of NGC4649. This feature is more prominent for red GCs, but still persists in the 2D distribution of blue GCs. High and low luminosity GCs also show some segregation along this arc, with high-luminosity GCs preferentially located in the southern end and low-luminosity GCs in the northern section of the arc. GC-LMXBs follow the anisotropy of red-GCs, where most of them reside; however, a significant overdensity of (high-luminosity) field LMXBs is present to the south of the GC arc. These results suggest that NGC4649 has experienced mergers and/or multiple accretions of less massive satellite galaxies during its evolution, of which the GCs in the arc may be the fossil remnant. We speculate that the observed anisotropy in the field LMXB spatial distribution indicates that these X-ray binaries may be the remnants of a star formation event connected with the merger, or maybe be ejected from the parent red GCs, if the bulk motion of these clusters is significantly affected by dynamical friction. We also detect a luminosity enhancement in the X-ray source population of the companion spiral galaxy NGC4647. We suggest that these may be younger high mass X-ray binaries formed as a result of the tidal interaction of this galaxy with NGC4649.

We investigate the radial distribution of the low-mass X-ray binary (LMXB) population in the elliptical galaxy NGC 4649, using Chandra and Hubble data to separate the field and globular cluster (GC) populations. GCs with LMXBs have the same radial distribution as the parent red and blue GCs. The radial profile of field LMXBs follows the V-band profile within the D25 of NGC 4649. Using the spatial information provided by our data, we find that the global galaxy-wide relations between cumulative number and luminosity of LMXBs and the integrated stellar mass hold on local scales within D25. An excess of field LMXBs with respect to the V-light is observed in the galaxy's outskirt, which may be partially due to unidentified GC sources or to a rejuvenated field LMXB population caused by past merging interactions.

X-ray photons, because of their long mean-free paths, can easily escape the galactic environments where they are produced, and interact at long distances with the inter-galactic medium, potentially having a significant contribution to the heating and reionization of the early Universe. The two most important sources of X-ray photons in the Universe are active galactic nuclei (AGN) and X-ray binaries (XRBs). In this Letter we use results from detailed, large scale population synthesis simulations to study the energy feedback of XRBs, from the first galaxies (z~ 20) until today. We estimate that X-ray emission from XRBs dominates over AGN at z>6-8. The shape of the spectral energy distribution of the emission from XRBs shows little change with redshift, in contrast to its normalization which evolves by ~4 orders of magnitude, primarily due to the evolution of the cosmic star-formation rate. However, the metallicity and the mean stellar age of a given XRB population affect significantly its X-ray output. Specifically, the X-ray luminosity from high-mass XRBs per unit of star-formation rate varies an order of magnitude going from solar metallicity to less than 10% solar, and the X-ray luminosity from low-mass XRBs per unit of stellar mass peaks at an age of ~300 Myr and then decreases gradually at later times, showing little variation for mean stellar ages > 3 Gyr. Finally, we provide analytical and tabulated prescriptions for the energy output of XRBs, that can be directly incorporated in cosmological simulations.

We present the largest-scale comparison to date between observed extragalactic X-ray binary (XRB) populations and theoretical models of their production. We construct observational X-ray luminosity functions (oXLFs) using Chandra observations of 12 late-type galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). For each galaxy, we obtain theoretical XLFs (tXLFs) by combining XRB synthetic models, constructed with the population synthesis code StarTrack, with observational star formation histories (SFHs). We identify highest-likelihood models both for individual galaxies and globally, averaged over the full galaxy sample. Individual tXLFs successfully reproduce about half of oXLFs, but for some galaxies we are unable to find underlying source populations, indicating that galaxy SFHs and metallicities are not well matched and/or XRB modeling requires calibration on larger observational samples. Given these limitations, we find that best models are consistent with a product of common envelope ejection efficiency and central donor concentration ~=0.1, and a 50% uniform -- 50% "twins" initial mass-ratio distribution. We present and discuss constituent subpopulations of tXLFs according to donor, accretor and stellar population characteristics. The galaxy-wide X-ray luminosity due to low-mass and high-mass XRBs, estimated via our best global model tXLF, follows the general trend expected from the L_x - star formation rate and L_x - stellar mass relations of Lehmer et al 2010. Our best models are also in agreement with modeling of the evolution both of XRBs over cosmic time and of the galaxy X-ray luminosity with redshift.

We present a new method for the determination of the two-dimensional (2D) projected spatial distribution of globular clusters (GCs) in external galaxies. This method is based on the K-Nearest Neighbor density estimator of Dressler (1980), complemented by MonteCarlo simulations to establish the statistical significance of the results. We apply this method to NGC4261, a "test galaxy" where significant 2D anisotropy in the GC distribution has been reported. We confirm that the 2D distribution of GC is not azimuthally isotropic. Moreover, we demonstrate that the 2D distribution departures from the average GC radial distribution results in highly significant spiral-like or broken shell features. Overall, the same perturbations are found in "red" and "blue" GCs, but with some differences. In particular, we observe a central feature, roughly aligned with the minor axis of NGC4261, composed of red and most luminous GCs. Blue and fainter GCs are more frequent at large radial distances and follow the spiral-like features of the overall density structure. These results suggest a complex merging history for NGC4261.

In this paper, we study the relationship between the 2-10 keV X-ray luminosity (L_X), assumed to originate from X-ray binaries (XRBs), and star formation rate (SFR) in UV-selected z<0.1 Lyman break analogs (LBAs). We present Chandra observations for four new GALEX-selected LBAs. Including previously studied LBAs, Haro 11 and VV 114, we find that LBAs demonstrate L_X/SFR ratios that are elevated by ~1.5sigma compared to local galaxies, similar to the ratios found for stacked LBGs in the early Universe (z>2). We show that these LBAs are unlikely to harbor AGN, based on their optical and X-ray spectra and the spatial distribution of the X-rays in three spatially extended cases. We expect that high-mass X-ray binaries (HMXBs) dominate the X-ray emission in these galaxies, based on their high specific SFRs (sSFRs=SFR/M* > 10^{-9}/yr), which suggest the prevalence of young stellar populations. Since both LBAs and LBGs have lower dust attenuations and metallicities compared to similar samples of more typical local galaxies, we investigate the effects of dust extinction and metallicity on the L_X/SFR for the broader population of galaxies with high sSFRs (>10^{-10}/yr). The estimated dust extinctions (corresponding to column densities of N_H<10^{22}/cm^2) are expected to have insignificant effects on observed L_X/SFR ratio for the majority of galaxy samples. We find that the observed relationship between L_X/SFR and metallicity appears consistent with theoretical expectations from X-ray binary population synthesis models. Therefore, we conclude that lower metallicities, related to more luminous HMXBs such as ultraluminous X-ray sources (ULXs), drive the elevated L_X/SFR observed in our sample of z<0.1 LBAs. The relatively metal-poor, active mode of star formation in LBAs and distant z>2 LBGs may yield higher total HMXB luminosity than found in typical galaxies in the local Universe.

Emission from X-ray binaries (XRBs) is a major component of the total X-ray luminosity of normal galaxies, so X-ray studies of high redshift galaxies allow us to probe the formation and evolution of X-ray binaries on very long timescales. In this paper, we present results from large-scale population synthesis models of binary populations in galaxies from z = 0 to 20. We use as input into our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history (SFH) and metallicity evolution of each galaxy. We run a grid of 192 models, varying all the parameters known from previous studies to affect the evolution of XRBs. We use our models and observationally derived prescriptions for hot gas emission to create theoretical galaxy X-ray luminosity functions (XLFs) for several redshift bins. Models with low CE efficiencies, a 50% twins mass ratio distribution, a steeper IMF exponent, and high stellar wind mass loss rates best match observational results from Tzanavaris & Georgantopoulos (2008), though they significantly underproduce bright early-type and very bright (Lx > 10d41) late-type galaxies. These discrepancies are likely caused by uncertainties in hot gas emission and SFHs, AGN contamination, and a lack of dynamically formed Low-mass XRBs. In our highest likelihood models, we find that hot gas emission dominates the emission for most bright galaxies. We also find that the evolution of the normal galaxy X-ray luminosity density out to z = 4 is driven largely by XRBs in galaxies with X-ray luminosities between 10d40 and 10d41 erg/s.

We present comprehensive observational results of the metallicity effect on the fraction of globular clusters (GC) that contain low-mass X-ray binaries (LMXB), by utilizing all available data obtained with Chandra for LMXBs and HST ACS for GCs. Our primary sample consists of old elliptical galaxies selected from the ACS Virgo and Fornax surveys. To improve statistics at both the lowest and highest X-ray luminosity, we also use previously reported results from other galaxies. It is well known that the LMXB fraction is considerably higher in red, metal-rich, than in blue, metal-poor GCs. In this paper, we test whether this metallicity effect is X-ray luminosity-dependent, and find that the effect holds uniformly in a wide luminosity range. This result is statistically significant (at >= 3 sigma) in LMXBs with luminosities in the range LX = 2 x 10^37 - 5 x 10^38 erg s-1, where the ratio of LMXB fractions in metal-rich to metal-poor GCs is R = 3.4 +- 0.5. A similar ratio is also found at lower (down to 10^36 erg s-1) and higher luminosities (up to the ULX regime), but with less significance (~2 sigma confidence). Because different types of LMXBs dominate in different luminosities, our finding requires a new explanation for the metallicity effect in dynamically formed LMXBs. We confirm that the metallicity effect is not affected by other factors such as stellar age, GC mass, stellar encounter rate, and galacto-centric distance.

High redshift galaxies permit the study of the formation and evolution of X-ray binary populations on cosmological timescales, probing a wide range of metallicities and star-formation rates. In this paper, we present results from a large scale population synthesis study that models the X-ray binary populations from the first galaxies of the universe until today. We use as input to our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history and metallicity evolution of the universe. Our modeling, which is constrained by the observed X-ray properties of local galaxies, gives predictions about the global scaling of emission from X-ray binary populations with properties such as star-formation rate and stellar mass, and the evolution of these relations with redshift. Our simulations show that the X-ray luminosity density (X-ray luminosity per unit volume) from X-ray binaries in our Universe today is dominated by low-mass X-ray binaries, and it is only at z>2.5 that high-mass X-ray binaries become dominant. We also find that there is a delay of ~1.1 Gyr between the peak of X-ray emissivity from low-mass Xray binaries (at z~2.1) and the peak of star-formation rate density (at z~3.1). The peak of the X-ray luminosity from high-mass X-ray binaries (at z~3.9), happens ~0.8 Gyr before the peak of the star-formation rate density, which is due to the metallicity evolution of the Universe.

We present the X-ray source catalog for the Chandra monitoring observations of the elliptical galaxy, NGC 4649. The galaxy has been observed with Chandra ACIS-S3 in six separate pointings, reaching a total exposure of 299 ks. There are 501 X-ray sources detected in the 0.3-8.0 keV band in the merged observation or in one of the six individual observations; 399 sources are located within the D_25 ellipse. The observed 0.3-8.0 keV luminosities of these 501 sources range from 9.3x10^{36} erg/s to 5.4x10^{39} erg/s. The 90% detection completeness limit within the D_25 ellipse is 5.5x10^{37} erg/s. Based on the surface density of background active galactic nuclei (AGNs) and detection completeness, we expect ~45 background AGNs among the catalog sources (~15 within the D_25 ellipse). There are nine sources with luminosities greater than 10^{39} erg/s, which are candidates for ultraluminous X-ray sources. The nuclear source of NGC 4649 is a low-luminosity AGN, with an intrinsic 2.0-8.0 keV X-ray luminosity of 1.5x10^{38} erg/s. The X-ray colors suggest that the majority of the catalog sources are low-mass X-ray binaries (LMXBs). We find that 164 of the 501 X-ray sources show long-term variability, indicating that they are accreting compact objects. We discover four transient candidates and another four potential transients. We also identify 173 X-ray sources (141 within the D_25 ellipse) that are associated with globular clusters (GCs) based on Hubble Space Telescope and ground-based data; these LMXBs tend to be hosted by red GCs. Although NGC 4649 has a much larger population of X-ray sources than the structurally similar early-type galaxies, NGC 3379 and NGC 4278, yet the X-ray source properties are comparable in all three systems.

We report the discovery of a new ultraluminous X-ray source (ULX) associated with a globular cluster in the elliptical galaxy NGC 4649. The X-ray source was initially detected with a luminosity below 5 x 10^38 erg/s, but in subsequent observations 7 and 11 years later it had brightened substantially to 2 - 3 x 10^39 erg/s. Over the course of six separate observations it displayed significant spectral variability, in both continuum slope and absorption column. Short-term variability in the X-ray flux was also present in at least one observation. The properties of this object appear consistent with a stellar-mass black hole accreting at super-Eddington rates (i.e. in the ultraluminous accretion state), although a highly super-Eddington neutron star cannot be excluded. The coincidence of an increase in absorption column with a possible enhancement in short-term variability in at least one observation is suggestive of a clumpy radiatively-driven wind crossing our line-of-sight to the object

We present g and z photometry and size estimates for globular clusters (GCs) in the massive Virgo elliptical NGC 4649 (M60) using a five-pointing Hubble Space Telescope/Advanced Camera for Surveys mosaic. The metal-poor GCs show a monotonic negative metallicity gradient of (-0.43 +/- 0.10) dex per dex in radius over the full radial range of the data, out to ~ 24 kpc. There is evidence for substantial color substructure among the metal-rich GCs. The metal-poor GCs have typical sizes ~ 0.4 pc larger than the metal-rich GCs out to large galactocentric distances (~> 20 kpc), favoring an intrinsic explanation for the size difference rather than projection effects. There is no clear relation between half-light radius and galactocentric distance beyond ~ 15 kpc, suggesting that the sizes of GCs are not generically set by tidal limitation. Finally, we identify ~ 20 candidate ultra-compact dwarfs that extend down to surprisingly faint absolute magnitudes (M_z ~ -8.5), and may bridge the gap between this class and "extended clusters" in the Local Group. Three of the brighter candidates have published radial velocities and can be confirmed as bona fide ultra-compact dwarfs; follow-up spectroscopy will determine the nature of the remainder of the candidates.

We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7, we expect that the observed properties of these LBGs are dominated by young stellar populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z<4 LBGs (covering ~90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with redshift (z) and SFR as log LX = 0.93 log (1+z) + 0.65 log SFR + 39.80. By comparing our observations with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution of LX/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 2-10 keV luminosity per Mpc^3) and find excellent agreement. While there are no significant stacked detections at z>5, we use our upper limits from 5<z<8 LBGs to constrain the SMBH accretion history of the Universe around the epoch of reionization.

Globular clusters (GCs) effectively produce dynamically-formed low-mass X-ray binaries (LMXBs). Observers detect ~100 times more LMXBs per stellar mass in GCs compared to stars in the fields of galaxies. It has also been observationally established that metal-rich GCs are about 3 times more likely to contain an X-ray source than their metal-poor counterparts. Recent observations have shown that this ratio holds in extragalactic GCs for all bright X-ray sources with Lx between 2X10^{37} and 5X10^{38} erg/s. In this Letter, we propose that the observed metallicity dependence of LMXBs in extragalactic GCs can be explained by the differences in the number densities and average masses of red giants in populations of different metallicities. Red giants serve as seeds for the dynamical production of bright LMXBs via two channels - binary exchanges and physical collisions - and the increase of the number densities and masses of red giants boost LMXB production, leading to the observed difference. We also discuss a possible effect of the age difference in stellar populations of different metallicities.

We report the spectral and temporal variability properties of 18 candidate transient and potential transient (TC and PTC) sources detected in deep multi-epoch Chandra observation of the nearby elliptical galaxies, NGC 3379, NGC 4278 and NGC 4697. Only one source can be identified with a background counterpart, leaving 17 TCs + PTCs in the galaxies. Of these, 14 are in the galaxy field, supporting the theoretical picture that the majority of field X-ray binaries (XRBs) will exhibit transient accretion for >75% of their lifetime. Three sources are coincident with globular clusters (GCs), including two high-luminosity candidate black hole (BH) XRBs, with Lx=5.4E38 erg/s, and Lx=2.8E39 erg/s, respectively. The spectra, luminosities and temporal behavior of these 17 sources suggest that the transient population is heterogeneous, including neutron star (NS) and BH XRBs in both normal and high-rate accretion modes, and super soft sources containing white dwarf binaries. Our TC and PTC detections are noticeably fewer that the number expected from the populations synthesis (PS) models of Fragos et al. (2009), tailored to our new Chandra pointings of NGC 4278. We attribute this discrepancy to the PS assumption that the transient population is composed of NS XRBs, as well as differences between the statistical analysis and error estimates used in the model and our observations.

We present Chandra studies of the X-ray binary (XRB) populations in the bulge and ring regions of the ring galaxy NGC 1291. We detect 169 X-ray point sources in the galaxy, 75 in the bulge and 71 in the ring, utilizing the four available Chandra observations totaling an effective exposure of 179 ks. We report photometric properties of these sources in a point-source catalog. There are ~40% of the bulge sources and ~25% of the ring sources showing >3\sigma long-term variability in their X-ray count rate. The X-ray colors suggest that a significant fraction of the bulge (~75%) and ring (~65%) sources are likely low-mass X-ray binaries (LMXBs). The spectra of the nuclear source indicate that it is a low-luminosity AGN with moderate obscuration; spectral variability is observed between individual observations. We construct 0.3-8.0 keV X-ray luminosity functions (XLFs) for the bulge and ring XRB populations, taking into account the detection incompleteness and background AGN contamination. We reach 90% completeness limits of ~1.5\times10^{37} and ~2.2\times10^{37} erg/s for the bulge and ring populations, respectively. Both XLFs can be fit with a broken power-law model, and the shapes are consistent with those expected for populations dominated by LMXBs. We perform detailed population synthesis modeling of the XRB populations in NGC 1291, which suggests that the observed combined XLF is dominated by an old LMXB population. We compare the bulge and ring XRB populations, and argue that the ring XRBs are associated with a younger stellar population than the bulge sources, based on the relative overdensity of X-ray sources in the ring, the generally harder X-ray color of the ring sources, the overabundance of luminous sources in the combined XLF, and the flatter shape of the ring XLF.

In recent years, accurate observational constraints become available for an increasing number of Galactic X-ray binaries. Together with proper motion measurements, we could reconstruct the full evolutionary history of X-ray binaries back to the time of compact object formation. In this paper, we present the first study of the persistent X-ray source Cygnus\;X-1 that takes into account of all available observational constraints. Our analysis accounts for three evolutionary phases: orbital evolution and motion through the Galactic potential after the formation of black hole (BH), and binary orbital dynamics at the time of core collapse. We find that the mass of the BH immediate progenitor is $15.0 - 20.0$ M$_\sun$, and at the time of core collapse, the BH has potentially received a small kick velocity of $\le 77$ km s$^{-1}$ at 95% confidence. If the BH progenitor mass is less than $\sim 17$ M$_\sun$, a non zero natal kick velocity is required to explain the currently observed properties of Cygnus\;X-1. Since the BH has only accreted mass from its companion's stellar wind, the negligible amount of accreted mass is impossible to explain the observationally inferred BH spin of $a_* > 0.95$, and the origin of this extreme BH spin must be connected to the BH formation itself. Right after the BH formation, we find that the BH companion is a $19.8 - 22.6$ M$_\sun$ main sequence star, orbiting the BH at a period of $4.7 - 5.2$ days. Furthermore, recent observations show that the BH companion is currently super-synchronized. This super-synchronism indicates that the strength of tides exerted on the BH companion should be weaker by a factor of at least two compared to the usually adopted strength.

Black hole (BH) X-ray binaries (XRBs) are X-ray luminous binary systems comprising a BH accreting matter from a companion star. Understanding their origins sheds light on the still not well understood physics of BH formation. M33 X-7 hosts one of the most massive stellar-mass BH among all XRBs known to date, a 15.65 Msun BH orbiting a 70 Msun companion star in a 3.45 day orbit. The high masses of the two components and the tight orbit relative to the large H-rich stellar component challenge our understanding of the typically invoked BH-XRBs formation channels. The measured underluminosity of the optical component further complicates the picture. A solution to the evolutionary history of this system that can account for all its observed properties has yet to be presented, and here we propose the first scenario that is consistent with the complete set of current observational constraints. In our model, M33 X-7 started its life hosting a 85-99 Msun primary and a 28-32 Msun companion in a Keplerian orbit of 2.8-3.1 days. In order to form a BH of 15.65 Msun, the initially most massive component transferred part of its envelope to the companion star and lost the rest in a strong stellar wind. During this dynamically stable mass transfer phase the companion accreted matter, to become the presently underluminous 70 Msun star.

M33 X-7 is among the most massive X-Ray binary stellar systems known, hosting a rapidly spinning 15.65 Msun black hole orbiting an underluminous 70 Msun Main Sequence companion in a slightly eccentric 3.45 day orbit. Although post-main-sequence mass transfer explains the masses and tight orbit, it leaves unexplained the observed X-Ray luminosity, star's underluminosity, black hole's spin, and eccentricity. A common envelope phase, or rotational mixing, could explain the orbit, but the former would lead to a merger and the latter to an overluminous companion. A merger would also ensue if mass transfer to the black hole were invoked for its spin-up. Here we report that, if M33 X-7 started as a primary of 85-99 Msun and a secondary of 28-32 Msun, in a 2.8-3.1 day orbit, its observed properties can be consistently explained. In this model, the Main Sequence primary transferred part of its envelope to the secondary and lost the rest in a wind; it ended its life as a ~16 Msun He star with a Fe-Ni core which collapsed to a black hole (with or without an accompanying supernova). The release of binding energy and, possibly, collapse asymmetries "kicked" the nascent black hole into an eccentric orbit. Wind accretion explains the X-Ray luminosity, while the black hole spin can be natal.

In black hole X-ray binaries, a misalignment between the spin axis of the black hole and the orbital angular momentum can occur during the supernova explosion that forms the compact object. In this letter we present population synthesis models of Galactic black hole X-ray binaries, and study the probability density function of the misalignment angle, and its dependence on our model parameters. In our modeling, we also take into account the evolution of misalignment angle due to accretion of material onto the black hole during the X-ray binary phase. The major factor that sets the misalignment angle for X-ray binaries is the natal kick that the black hole may receive at its formation. However, large kicks tend to disrupt binaries, while small kicks allow the formation of XRBs and naturally select systems with small misalignment angles. Our calculations predict that the majority (>67%) of Galactic field BH XRBs have rather small (>10 degrees) misalignment angles, while some systems may reach misalignment angles as high as ~90 degrees and even higher. This results is robust among all population synthesis models. The assumption of small small misalignment angles is extensively used to observationally estimate black hole spin magnitudes, and for the first time we are able to confirm this assumption using detailed population synthesis calculations.

(abridged) We present a statistical study of the low-mass X-ray binary (LMXB) populations of three nearby, old elliptical galaxies: NGC 3379, NGC 4278, and NGC 4697. With a cumulative ~1 Ms Chandra ACIS observing time, we detect 90-170 LMXBs within the D25 ellipse of each galaxy. Cross-correlating Chandra X-ray sources and HST optical sources, we identify 75 globular cluster (GC) LMXBs and 112 field LMXBs. In the low luminosity range allowed by our deeper data (LX < 5 x 1037 erg s-1), we find a significant relative lack of GC-LMXBs, when compared with field sources. Using the co-added sample from the three galaxies, we find that the incompleteness-corrected X-ray luminosity functions (XLFs) of GC and field LMXBs differ at ~4# significance at LX < 5 x 1037 erg s-1. As previously reported, these XLFs are consistent at higher luminosities. Our observations may indicate a potential predominance of GC-LMXBs with donors evolved beyond the main sequence, when compared to current models, but their efficient formation requires relatively high initial binary fractions in clusters. The field LMXB XLF can be fitted with either a single power-law model plus a localized excess at a luminosity of 5-6 x 1037 erg s-1, or a broken power-law with a similar low-luminosity break. This XLF may be explained with NS-red-giant LMXBs, contributing to ~15% of total LMXBs population at ~5x1037 erg s-1. The difference in the GC and field XLFs is consistent with different origins and/or evolutionary paths between the two LMXB populations, although a fraction of the field sources are likely to have originated in GCs.

We explore the formation and evolution of the black hole X-ray binary system M33 X-7. In particular, we examine whether accounting for systematic errors in the stellar parameters inherent to single star models, as well as the uncertainty in the distance to M33, can explain the discrepancy between the observed and expected luminosity of the ~70 solar masses companion star. Our analysis assumes no prior interactions between the companion star and the black hole progenitor. We use four different stellar evolution codes, modified to include a variety of current stellar wind prescriptions. For the models satisfying the observational constraints on the donor star's effective temperature and luminosity, we recalculate the black hole mass, the orbital separation, and the mean X-ray luminosity. Our best model, satisfying simultaneously all observational constraints except the observationally inferred companion mass, consists of a ~13 solar masses black hole and a ~54 solar masses companion star. We conclude that a star with the observed mass and luminosity can not be explained via single star evolution models, and that a prior interaction between the companion star and the black hole progenitor should be taken into account.

We present the properties of the discrete X-ray sources detected in our monitoring program of the globular cluster (GC) rich elliptical galaxy, NGC 4278, observed with Chandra ACIS-S in six separate pointings, resulting in a co-added exposure of 458-ks. From this deep observation, 236 sources have been detected within the region overlapped by all observations, 180 of which lie within the D25 ellipse of the galaxy. These 236 sources range in Lx from 3.5E36 erg/s (with 3sigma upper limit <1E37 erg/s) to ~2E40 erg/s, including the central nuclear source which has been classified as a LINER. From optical data, 39 X-ray sources have been determined to be coincident with a globular cluster, these sources tend to have high X-ray luminosity, with ten of these sources exhibiting Lx>1E38 erg/s. From X-ray source photometry, it has been determined that the majority of the 236 point sources that have well constrained colors, have values that are consistent with typical LMXB spectra, with 29 of these sources expected to be background objects from the logN-logS relation. There are 103 sources in this population that exhibit long-term variability, indicating that they are accreting compact objects. 3 of these sources have been identified as transient candidates, with a further 3 possible transients. Spectral variations have also been identified in the majority of the source population, where a diverse range of variability has been identified, indicating that there are many different source classes located within this galaxy.

We propose a physically motivated and self-consistent prescription for the modeling of transient neutron star (NS) low-mass X-ray binary (LMXB) properties, such as duty cycle (DC), outburst duration and recurrence time. We apply this prescription to the population synthesis (PS) models of field LMXBs presented by Fragos et al. (2008), and compare the transient LMXB population to the Chandra X-ray survey of the two elliptical galaxies NGC 3379 and NGC 4278, which revealed several transient sources (Brassington et al., 2008, 2009). We are able to exclude models with a constant DC for all transient systems, while models with a variable DC based on the properties of each system are consistent with the observed transient populations. We predict that the majority of the observed transient sources in these two galaxies are LMXBs with red giant donors. Our comparison suggests that LMXBs formed through evolution of primordial field binaries are dominant in globular cluster (GC) poor elliptical galaxies, while they still have a significant contribution in GC rich ones.

In recent years, an increasing number of proper motions have been measured for Galactic X-ray binaries. When supplemented with accurate determinations of the component masses, orbital period, and donor effective temperature, these kinematical constraints harbor a wealth of information on the system's past evolution. Here, we consider all this available information to reconstruct the full evolutionary history of the black hole X-ray binary XTE J1118+480, assuming that the system originated in the Galactic disk and the donor has solar metallicity. This analysis accounts for four evolutionary phases: mass transfer through the ongoing X-ray phase, tidal evolution before the onset of Roche-lobe overflow, motion through the Galactic potential after the formation of the black hole, and binary orbital dynamics due to explosive mass loss and possibly a black hole natal kick at the time of core collapse. We find that right after black hole formation, the system consists of a ~6.0-10.0 solar masses black hole and a ~1.0-1.6 solar masses main-sequence star. We also find that that an asymmetric natal kick is not only plausible but required for the formation of this system, and derive a lower and upper limit on the black hole natal kick velocity magnitude of 80 km/s and 310 km/s, respectively.

We present the properties of the discrete X-ray sources detected in our monitoring program of the `typical' elliptical galaxy, NGC 3379, observed with Chandra ACIS-S in five separate pointings, resulting in a co-added exposure of 324-ks. From this deep observation, 132 sources have been detected within the region overlapped by all observations, 98 of which lie within the D25 ellipse of the galaxy. These 132 sources range in Lx from 6E35 erg/s (with 3sigma upper limit <4E36 erg/s) to ~2E39 erg/s, including one source with Lx>1E39 erg/s, which has been classified as a ULX. From optical data, 10 X-ray sources have been determined to be coincident with a globular cluster, these sources tend to have high X-ray luminosity, with three of these sources exhibiting Lx> 1E38 erg/s. From X-ray source photometry, it has been determined that the majority of the 132 sources that have well constrained colors, have values that are consistent with typical LMXB spectra. Additionally to this, a sub-population of 10 sources has been found to exhibit very hard spectra and it is expected that most of these sources are absorbed background AGN. There are 64 sources in this population that exhibit long-term variability, indicating that they are accreting compact objects. Five of these sources have been identified as transient candidates, with a further 3 possible transients. Spectral variations have also been identified in the majority of the source population, where a diverse range of variability has been identified, indicating that there are many different source classes located within this galaxy.

We present theoretical models for the formation and evolution of populations of low-mass X-ray binaries (LMXB) in the two elliptical galaxies NGC 3379 and NGC 4278. The models are calculated with the recently updated StarTrack code (Belczynski et al., 2006), assuming only a primordial galactic field LMXB population. StarTrack is an advanced population synthesis code that has been tested and calibrated using detailed binary star calculations and incorporates all the important physical processes of binary evolution. The simulations are targeted to modeling and understanding the origin of the X-ray luminosity functions (XLF) of point sources in these galaxies. For the first time we explore the population XLF down to luminosities of 3X10^36 erg/s, as probed by the most recent observational results (Kim et al., 2006). We consider models for the formation and evolution of LMXBs in galactic fields with different CE efficiencies, stellar wind prescriptions, magnetic braking laws and initial mass functions. We identify models that produce an XLF in excellent agreement with the observations both in shape and absolute normalization. We also find that the treatment of the outburst luminosity of transient systems remains a crucial factor for the determination of the XLF since the modeled populations are dominated by transient X-ray systems.

In recent years, an increasing number of proper motions have been measured for Galactic X-ray binaries. When supplemented with accurate determinations of the component masses, orbital period, and donor luminosity and effective temperature, these kinematical constraints harbor a wealth of information on the systems' past evolution. The constraints on compact object progenitors and kicks derived from this are of immense value for understanding compact object formation and exposing common threads and fundamental differences between black hole and neutron star formation. Here, we present the results of such an analysis for the black hole X-ray binary XTE J1118+480. We present results from modeling the mass transfer phase, following the motion in the Galaxy back to the birth site of the black hole, and examining the dynamics of symmetric and asymmetric core-collapses of the black hole progenitor.

We revisit the formation of PSR J0737-3039, taking into account the most recent observational constraints. We show that the most likely kick velocity and progenitor parameters depend strongly on the consideration of the full five-dimensional PDF for the magnitude and direction of the kick velocity imparted to pulsar B at birth, the mass of pulsar B's pre-supernova helium star progenitor, and the pre-supernova orbital separation, and on the adopted prior assumptions. The priors consist of the transverse systemic velocity, the age of the system, and the treatment of the unknown radial velocity. Since the latter cannot be determined from observation, we adopt a statistical approach and use theoretical radial-velocity distributions obtained from population synthesis calculations for coalescing double neutron stars. We find that the prior assumptions about the pre-supernova helium star mass affect the derived most likely parameters significantly: when the minimum helium star mass required for neutron star formation is assumed to be 2.1Msun, the most likely kick velocity ranges from 70-180km/s; when masses lower than 2.1Msun are assumed to allow neutron star formation, the most likely kick velocity can be as low as a few km/s, although the majority of the considered models still yield most likely kick velocities of 50-170km/s. We also show that the proximity of the double pulsar to the Galactic plane and the small proper motion do not pose stringent constraints on the kick velocity and progenitor mass of pulsar B. Instead, the constraints imposed by the orbital dynamics of asymmetric supernova explosions turn out to be much more restrictive. We conclude that the currently available observational constraints cannot be used to favor a specific core-collapse and neutron star formation mechanism. (abridged)

We report recent Chandra observations of the ULX in the elliptical galaxy NGC3379 that clearly detect two flux variability cycles. Comparing these data with the Chandra observation of ~5yr ago, we measure a flux modulation with a period of ~12.6hr. Moreover, we find that the emission undergoes a correlated spectral modulation, becoming softer at low flux. We argue that our results establish this source as a ULX binary in NGC3379. Given the old stellar population of this galaxy, the ULX is likely to be a soft transient, however historical X-ray sampling suggests that the current on phase has lasted ~10yr. We discuss our results in terms of ADC and wind-feedback models. We constrain the donor mass and orbital period at the onset of mass transfer within 1.15-1.4M and 12.5-16hr, respectively. The duration of the mass-transfer phase so far is probably ~1Gyr and the binary has been a soft X-ray transient throughout this time. These constraints are insensitive to the mass of the accretor.

The formation and evolution of active regions is an inherently complex phenomenon. Magnetic fields generated at the base of the convection zone follow a chaotic evolution before reaching the solar surface. In this article, we use a 2-D probabilistic Cellular Automaton (CA) to model the statistical properties of the magnetic patterns formed on the solar surface and to estimate the magnetic energy released in the interaction of opposite polarities. We assume that newly emerged magnetic flux tubes stimulate the emergence of new magnetic flux in their neighborhood. The flux-tubes move randomly on the surface of the sun, and they cancel and release their magnetic energy when they collide with magnetic flux of opposite polarity, or diffuse into the ``empty'' photosphere. We assume that cancellation of magnetic flux in collisions causes "flares" and determine the released energy as the difference in the square of the magnetic field flux. The statistics of the simulated "flares" follow a power-law distribution in energy, with power-law index a=2.2 +/- 0.1. The size distribution function of the simulated active regions exhibits a power law behavior with index k = 1.93 +/- 0.08, and the fractal dimension of the magnetized areas on the simulated solar surface is close to D = 1.42 +/- 0.12. Both quantities, D and k, are inside the range of the observed values.