Search

Type IIb supernovae (SNe IIb) present a unique opportunity for investigating the evolutionary channels and mechanisms governing the evolution of stripped-envelope SN progenitors due to a variety of observational constraints available. Comparison of these constraints with the full distribution of theoretical properties not only help ascertain the prevalence of observed properties in nature, but can also reveal currently unobserved populations. In this follow-up paper, we use the large grid of models presented in Sravan et al. 2019 to derive distributions of single and binary SNe IIb progenitor properties and compare them to constraints from three independent observational probes: multi-band SN light-curves, direct progenitor detections, and X-ray/radio observations. Consistent with previous work, we find that while current observations exclude single stars as SN IIb progenitors, SN IIb progenitors in binaries can account for them. We also find that the distributions indicate the existence of an unobserved dominant population of binary SNe IIb at low metallicity that arise due to mass transfer initiated on the Hertzsprung Gap. In particular, our models indicate the existence of a group of highly stripped (envelope mass ~0.1-0.2 M_sun) progenitors that are compact (<50 R_sun) and blue (T_eff <~ 10^5K) with ~10^4.5-10^5.5 L_sun and low density circumstellar mediums. As discussed in Sravan et al. 2019, this group is necessary to account for SN IIb fractions and likely exist regardless of metallicity. The detection of the unobserved populations indicated by our models would support weak stellar winds and inefficient mass transfer in SN IIb progenitors., Comment: Resubmitted to the Astrophysical Journal after incorporating suggestions from the referee

GW170817 began gravitational-wave multimessenger astronomy. However, GW170817 will not be representative of detections in the coming years -- typical gravitational-wave sources will be closer the detection horizon, have larger localization regions, and (when present) will have correspondingly weaker electromagnetic emission. In its design state, the gravitational-wave detector network in the mid-2020s will consist of up to five similar-sensitivity second-generation interferometers. The instantaneous sky-coverage by the full network is nearly isotropic, in contrast to the configuration during the first \change{three} observing runs. Along with the coverage of the sky, there are also commensurate increases in the average horizon for a given binary mass. We present a realistic set of localizations for binary neutron stars and neutron star--black hole binaries, incorporating intra-network duty cycles and selection effects on the astrophysical distributions. Based on the assumption of an $80\%$ duty cycle, and that two instruments observe a signal above the detection threshold, we anticipate a median of $28$ sq.\ deg.\ for binary neutron stars, and $50$--$120$ sq.\ deg.\ for neutron star--black hole (depending on the population assumed). These distributions have a wide spread, and the best localizations, even for networks with fewer instruments, will have localizations of $1$--$10$ sq.\ deg.\ range. The full five instrument network reduces localization regions to a few tens of degrees at worst., Comment: 20 pages, 8 figures, 3 tables, accepted in ApJ

Observations with next-generation ground-based detectors further enhanced with multi-messenger (electromagnetic and neutrino) detections will allow us to probe new extreme astrophysics. Target sources included: core-collapse supernovae, continuous emission from isolated or accreting neutron stars, and bursts from magnetars and other pulsars., Comment: 13 pages, 1 figure, White Paper Submitted to Astro2020 (2020 Astronomy and Astrophysics Decadal Survey) by GWIC 3G Science Case Team (GWIC: Gravitational Wave International Committee)

Next-generation observations will revolutionize our understanding of binary black holes and will detect new sources, such as intermediate-mass black holes. Primary science goals include: Discover binary black holes throughout the observable Universe; Reveal the fundamental properties of black holes; Uncover the seeds of supermassive black holes., Comment: 14 pages, 3 figures, White Paper Submitted to Astro2020 (2020 Astronomy and Astrophysics Decadal Survey) by GWIC 3G Science Case Team (GWIC: Gravitational Wave International Committee)

Pair-instability and pulsational pair-instability supernovae (PPISN) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. If these BBHs are the product of binary evolution, then PPISNe could occur in very close binaries. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion, and what impact these events have on the formation of merging BBHs through binary evolution. For this, we have computed a set of models of metal-poor ($Z_\odot/10$) single helium stars using the \texttt{MESA} software instrument. For PPISN progenitors with pre-pulse masses $>50M_\odot$ we find that, after a pulse, heat deposited throughout the layers of the star that remain bound cause it to expand to more than $100R_\odot$ for periods of $10^2-10^4\;$~yrs depending on the mass of the progenitor. This results in long-lived phases of Roche-lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated to PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final black hole spin by $\sim 30\%$, induces eccentricities below the threshold of detectability of the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors., Comment: 22 pages, 17 figures, including 4 pages of appendices

Type IIb supernovae (SNe) are important candidates to understand mechanisms that drive the stripping of stripped-envelope (SE) supernova (SN) progenitors. While binary interactions and their high incidence are generally cited to favor them as Type IIb SN progenitors, this idea has not been tested using models covering a broad parameter space. In this paper, we use non-rotating single- and binary-star models at solar and low metallicities spanning a wide parameter space in primary mass, mass ratio, orbital period, and mass transfer efficiencies. We find that our single- and binary-star models contribute to roughly equal, however small, numbers of Type IIb SNe at solar metallicity. Binaries only dominate as progenitors at low metallicity. We also find that our models can account for less than half the observationally inferred rate for Type IIb SNe at solar metallicity, with computed rates ~<4% of core-collapse (CC) SNe. On the other hand, our models can account for the rates currently indicated by observations at low metallicity, with computed rates as high as 15% of CC SNe. However, this requires low mass transfer efficiencies (~<0.1) to prevent most progenitors from entering contact. We suggest that the stellar wind mass-loss rates at solar metallicity used in our models are too high. Lower mass-loss rates would widen the parameter space for binary Type IIb SNe at solar metallicity by allowing stars that initiate mass transfer earlier in their evolution to reach CC without getting fully stripped., Comment: Accepted to The Astrophysical Journal

A milestone of multi-messenger astronomy has been achieved with the detection of gravitational waves from a binary neutron star merger accompanied by observations of several associated electromagnetic counterparts. Joint observations can reveal details of the engines that drive the electromagnetic and gravitational-wave emission. However, locating and identify an electromagnetic counterparts to a gravitational-wave event is heavily reliant on localization of the source through gravitational-wave information. We explore the sky localization of a simulated set of neutron star mergers as the worldwide network of gravitational-wave detectors evolves through the next decade, performing the first such study for neutron star -- black hole binary sources. Currently, three detectors are observing with additional detectors in Japan and India expected to become operational in the coming years. With three detectors, we recover a median neutron star -- black hole binary sky localization of 60 deg$^2$ at the 90\% credible level. As all five detectors become operational, sources can be localized to a median of 11 deg$^2$ on the sky., Comment: 8 pages, 5 figures, accepted in ApJL

We demonstrate a method to fully characterize mass-transferring double white dwarf (DWD) systems with a helium-rich (He) WD donor based on the mass--radius relationship for He WDs. Using a simulated Galactic population of DWDs, we show that donor and accretor masses can be inferred for up to $\sim\, 60$ systems observed by both Laser Interferometer Space Antenna (LISA) and Gaia. Half of these systems will have mass constraints $\Delta\,M_{\rm{D}}\lesssim0.2M_{\odot}$ and $\Delta\,M_{\rm{A}}\lesssim2.3\,M_{\odot}$. We also show how the orbital frequency evolution due to astrophysical processes and gravitational radiation can be decoupled from the total orbital frequency evolution for up to $\sim 50$ of these systems., Comment: Updated to match publication in ApJ Letters; 8 pages, 4 figures

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible., Comment: 312 pages, 90 figures. Browse the current version at https://github.com/LSSTScienceCollaborations/ObservingStrategy, new contributions welcome!

Type IIb supernovae (SNe) present a unique opportunity for understanding the progenitors of stripped-envelope (SE) SNe as the stellar progenitor of several Type IIb SNe have been identified in pre-explosion images. In this paper, we use Bayesian inference and a large grid of non-rotating solar-metallicity single and binary stellar models to derive the associated probability distributions of single and binary progenitors of the Type IIb SN 2016gkg using existing observational constraints. We find that potential binary star progenitors have smaller pre-SN hydrogen-envelope and helium-core masses than potential single-star progenitors typically by 0.1 Msun and 2 Msun, respectively. We find that, a binary companion, if present, is a main-sequence or red-giant star. Apart from this, we do not find strong constraints on the nature of the companion star. We demonstrate that the range of progenitor helium-core mass inferred from observations could help improve constraints on the progenitor. We find that the probability that the progenitor of SN 2016gkg was a binary is 22% when we use constraints only on the progenitor luminosity and effective temperature. Imposing the range of pre-SN progenitor hydrogen-envelope mass and radius inferred from SN light-curves the probability the progenitor is a binary increases to 44%. However, there is no clear preference for a binary progenitor. This is in contrast to binaries being the currently favored formation channel for Type IIb SNe. Our analysis demonstrates the importance of statistical inference methods to constrain progenitor channels., Comment: Accepted to The Astrophysical Journal Letters

We explore the long-term evolution of mass-transferring white dwarf binaries undergoing both direct-impact and disk accretion and explore implications of such systems to gravitational wave astronomy. We cover a broad range of initial component masses and show that these systems, the majority of which lie within the LISA sensitivity range, exhibit prominent negative orbital frequency evolution (chirp) for a significant fraction of their lifetimes. Using a galactic population synthesis, we predict ~$2700$ double white dwarfs will be observable by LISA with negative chirps less than $-0.1\ yr^{-2}$. We also show that detections of mass-transferring double white dwarf systems by LISA may provide astronomers with unique ways of probing the physics governing close compact object binaries., Comment: Accepted in The Astrophysical Journal

Many giant exoplanets are found near their Roche limit and in mildly eccentric orbits. In this study we examine the fate of such planets through Roche-lobe overflow as a function of the physical properties of the binary components, including the eccentricity and the asynchronicity of the rotating planet. We use a direct three-body integrator to compute the trajectories of the lost mass in the ballistic limit and investigate the possible outcomes. We find three different outcomes for the mass transferred through the Lagrangian point $L_{1}$: (i) self-accretion by the planet, (ii) direct impact on the stellar surface, (iii) disk formation around the star. We explore the parameter space of the three different regimes and find that at low eccentricities, $e\lesssim 0.2$, mass overflow leads to disk formation for most systems, while for higher eccentricities or retrograde orbits self-accretion is the only possible outcome. We conclude that the assumption often made in previous work that when a planet overflows its Roche lobe it is quickly disrupted and accreted by the star is not always valid., Comment: 6 pages, 4 Figures, Accepted for publication to ApJ

Gravitational waves (GWs) from binary black hole (BBH) mergers provide a new probe of massive-star evolution and the formation channels of binary compact objects. By coupling the growing sample of BBH systems with population synthesis models, we can begin to constrain the parameters of such models and glean unprecedented knowledge about the inherent physical processes that underpin binary stellar evolution. In this study, we apply a hierarchical Bayesian model to mass measurements from a synthetic GW sample to constrain the physical prescriptions in population models and the relative fraction of systems generated from various channels. We employ population models of two canonical formation scenarios in our analysis --- isolated binary evolution involving a common-envelope phase and dynamical formation within globular clusters --- with model variations for different black hole natal kick prescriptions. We show that solely with chirp mass measurements, it is possible to constrain natal kick prescriptions and the relative fraction of systems originating from each formation channel with $\mathcal{O}(100)$ of confident detections. This framework can be extended to include additional formation scenarios, model parameters, and measured properties of the compact binary., Comment: 12 pages, 8 figures

The detection of electromagnetic counterparts to gravitational waves has great promise for the investigation of many scientific questions. It has long been hoped that in addition to providing extra, non-gravitational information about the sources of these signals, the detection of an electromagnetic signal in conjunction with a gravitational wave could aid in the analysis of the gravitational signal itself. That is, knowledge of the sky location, inclination, and redshift of a binary could break degeneracies between these extrinsic, coordinate-dependent parameters and the physical parameters, such as mass and spin, that are intrinsic to the binary. In this paper, we investigate this issue by assuming a perfect knowledge of extrinsic parameters, and assessing the maximal impact of this knowledge on our ability to extract intrinsic parameters. However, we find only modest improvements in a few parameters --- namely the primary component's spin --- and conclude that, even in the best case, the use of additional information from electromagnetic observations does not improve the measurement of the intrinsic parameters significantly., Comment: 7 pages, 4 figures, to appear in ApJ

The recent detections of GW150914 and GW151226 imply an abundance of stellar-mass binary-black-hole mergers in the local universe. While ground-based gravitational-wave detectors are limited to observing the final moments before a binary merges, space-based detectors, such as the Laser Interferometer Space Antenna (LISA), can observe binaries at lower orbital frequencies where such systems may still encode information about their formation histories. In particular, the orbital eccentricity and mass of binary black holes in the LISA frequency band can be used together to discriminate between binaries formed in isolation in galactic fields and those formed in dense stellar environments such as globular clusters. In this letter, we explore the orbital eccentricity and mass of binary-black-hole populations as they evolve through the LISA frequency band. Overall we find that there are two distinct populations discernible by LISA. We show that up to ~90% of binaries formed either dynamically or in isolation have eccentricities measurable by LISA. Finally, we note how measured eccentricities of low-mass binary black holes evolved in isolation could provide detailed constraints on the physics of black-hole natal kicks and common-envelope evolution., Comment: 7 pages, 4 figures, accepted for publication in ApJL

Massive stars are usually found in binaries, and binaries with periods less than 10 days may have a preference for near equal component masses. In this paper we investigate the evolution of these binaries all the way to contact and the possibility that these systems can be progenitors of double neutron star binaries. The small orbital separations of observed double neutron star binaries suggest that the progenitor systems underwent a common envelope phase at least once during their evolution. Bethe & Brown (1998) proposed that massive binary twins will undergo a common envelope evolution while both components are ascending the red giant branch or asymptotic giant branch simultaneously, also known as double-core evolution. Using models generated from the stellar evolution code Evolve Zero Age Main Sequence, we determine the range of mass ratios resulting in both components simultaneously ascending the RGB or AGB as a function of the difference in birth times, t. We find that, even for a generous t=5 Myr, the minimum mass ratio qmin=0.933 for an 8 Solar Mass primary and increases for larger primaries. We use a hydrodynamics code, StarSmasher, to study specifically the evolution of q=1 common envelope systems as a function of initial component mass, age, and orbital separation. We find the dynamical stability limit, the largest orbital separation where the binary becomes dynamically unstable, as a function of the component mass and age. Finally, we calculate the efficiency of ejecting matter during the inspiral phase to extrapolate the properties of the remnant binary from our numerical results, assuming the common envelope is completely ejected. We find that for the nominal core masses, there is a minimum orbital separation for a given component mass such that the helium cores survive common envelope evolution in a tightly bound binary and are viable progenitors for double neutron stars., Comment: 51 pages, 26 figures, accepted for publication in ApJ

We calculate the long-term evolution of angular momentum in double white dwarf binaries undergoing direct impact accretion over a broad range of parameter space. We allow the rotation rate of both components to vary, and account for the exchange of angular momentum between the spins of the white dwarfs and the orbit, while conserving the total angular momentum. We include gravitational, tidal, and mass transfer effects in the orbital evolution, and allow the Roche radius of the donor star to vary with both the stellar mass and the rotation rate. We examine the long-term stability of these systems, focusing in particular on those systems that may be progenitors of AM CVn or Type Ia Supernovae. We find that our analysis yields an increase in the predicted number of stable systems compared to that in previous studies. Additionally, we find that by properly accounting for the effects of asynchronism between the donor and the orbit on the Roche-lobe size, we eliminate oscillations in the orbital parameters which are found in previous studies. Removing these oscillations can reduce the peak mass transfer rate in some systems, keeping them from entering an unstable mass transfer phase., Comment: Manuscript accepted for publication by The Astrophysical Journal. Comments welcome

Although the details of the core-collapse supernova mechanism are not fully understood, it is generally accepted that the energy released in the collapse produces a shock that disrupts the star and produces the explosion. Some of the stellar material does not receive enough energy to escape the potential well of the newly formed neutron star and it falls back on to the core. This fallback plays an important role in setting compact remnant masses and has been invoked to explain a number of observed phenomena: late-time neutrino emission, r-process element production, peculiar supernovae, and long-duration gamma-ray bursts. To examine the link between fallback and these observed phenomena, we must understand the nature of fallback and currently, multiple mechanisms have been proposed that lead to different characteristics of the fallback. In this paper, we model a series of three-dimensional explosions by following the fallback for two different massive progenitors to study both the nature of the fallback and the role of asymmetric explosions on this fallback. One of the most-cited mechanisms behind fallback argues that it is caused by the reverse shock produced when the outward moving supernova shock decelerates in the relatively flat density profile of the stellar envelope. We show that this mechanism plays a minor, if any, role in driving fallback and most supernova fallback is produced by prompt mechanisms., Comment: 37 pages, 9 figures, submitted to ApJ

Although not nearly as numerous as binaries with two white dwarfs, eccentric neutron star-white dwarf (NS-WD) binaries are important gravitational-wave (GW) sources for the next generation of space-based detectors sensitive to low frequency waves. Here we investigate periastron precession in these sources as a result of general relativistic, tidal, and rotational effects; such precession is expected to be detectable for at least some of the detected binaries of this type. Currently, two eccentric NS-WD binaries are known in the galactic field, PSR J1141-6545 and PSR B2303+46, both of which have orbits too wide to be relevant in their current state to GW observations. However, population synthesis studies predict the existence of a significant Galactic population of such systems. Though small in most of these systems, we find that tidally induced periastron precession becomes important when tides contribute to more than 3% of the total precession rate. For these systems, accounting for tides when analyzing periastron precession rate measurements can improve estimates of the WD component mass inferred and, in some cases, will prevent us from misclassifying the object. However, such systems are rare due to rapid orbital decay. To aid the inclusion of tidal effects when using periastron precession as a mass measurement tool, we derive a function that relates the WD radius and periastron precession constant to the WD mass., Comment: Published in The Astrophysical Journal

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localisation have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high-latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyse methods for assessing the self-consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localisation. We apply these methods to two existing sky localisation techniques representing the two above-mentioned categories, using a set of simulated inspiral-only signals from compact binary systems with total mass $\le 20\,M_\odot$ and non-spinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor $\approx 20$ smaller for fully coherent techniques than for the specific variant of the "triangulation-based" technique used during the last science runs, at the expense of a factor $\approx 1000$ longer processing time.

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., Comment: 15 pages, 9 figures, submitted to ApJ

Here we present CAFein, a new computational tool for investigating radiative dissipation of dynamic tides in close binaries and of non-adiabatic, non-radial stellar oscillations in isolated stars in the linear regime. For the latter, CAFein computes the non-adiabatic eigenfrequencies and eigenfunctions of detailed stellar models. The code is based on the so-called Riccati method, a numerical algorithm that has been successfully applied to a variety of stellar pulsators, and which doesn't suffer of the major drawbacks of commonly used shooting and relaxation schemes. Here we present an extension of the Riccati method to investigate dynamic tides in close binaries. We demonstrate CAFein's capabilities as a stellar pulsation code both in the adiabatic and non-adiabatic regime, by reproducing previously published eigenfrequencies of a polytrope, and by successfully identifying the unstable modes of a stellar model in the $\beta$ Cephei/SPB region of the Hertzsprung-Russell diagram. Finally, we verify CAFein's behavior in the dynamic tides regime by investigating the effects of dynamic tides on the eigenfunctions and orbital and spin evolution of massive Main Sequence stars in eccentric binaries, and of hot Jupiter host stars. The plethora of asteroseismic data provided by the NASA's Kepler satellite, some of which include the direct detection of tidally excited stellar oscillations, make CAFein quite timely. Furthermore, the increasing number of observed short-period detached double white dwarfs (WD) and the observed orbital decay in the tightest of such binaries open up a new possibility of investigating WD interiors through the effects of tides on their orbital evolution, Comment: 17 pages, 9 figures, 6 tables. Accepted for publication in ApJ. Revised to match the accepted version. A new section on dynamic tides in a binary hosting a hot Jupiter has been added

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., Comment: ApJ in press. For redacted long table 1, see: http://www.pa.msu.edu/~strader/4649/table.tex

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

We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parallel random number generation scheme, as well as a parallel sorting algorithm, required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. The implementation uses the Message Passing Interface (MPI) library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude, from 10^5 to 10^7. We find that our results are in good agreement with self-similar core-collapse solutions, and the core collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within less than 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N=10^5, 128 for N=10^6 and 256 for N=10^7. The runtime reaches a saturation with the addition of more processors beyond these limits which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60x, 100x, and 220x, respectively., Comment: 53 pages, 13 figures, accepted for publication in ApJ Supplements

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., Comment: 14 pages, 8 figures, Accepted for publication in ApJ. Main changes from last revision: updated Figure 2, added section 4.3

We explore possible systematic errors in the mass measurements of stellar mass black holes. We find that significant errors can arise from the assumption of zero or constant emission from the accretion flow, which is commonly used when determining orbital inclination by modelling ellipsoidal variations. For A0620-00, the system with the best available data, we show that typical data sets and analysis procedures can lead to systematic underestimates of the inclination by ten degrees or more. A careful examination of the available data for the 15 other X-ray transients with low-mass donors suggests that this effect may significantly reduce the black hole mass estimates in several other cases, most notably that of GRO J0422+32. With these revisions, our analysis of the black hole mass distribution in soft X-ray transients does not suggest any "mass gap" between the low end of the distribution and the maximum theoretical neutron star mass, as has been identified in previous studies. Nevertheless, we find that the mass distribution retains other previously identified characteristics, namely a peak around 8M\odot, a paucity of sources with masses below 5M\odot, and a sharp drop-off above 10M\odot., Comment: 18 pages, 10 figures; submitted to ApJ

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., Comment: 11 pages, 11 figures, accepted by ApJ

One of the goals of gravitational-wave astronomy is simultaneous detection of gravitational-wave signals from merging compact-object binaries and the electromagnetic transients from these mergers. With the next generation of advanced ground-based gravitational wave detectors under construction, we examine the benefits of the proposed extension of the detector network to include a fourth site in Australia in addition to the network of Hanford, Livingston and Cascina sites. Using Bayesian parameter-estimation analyses of simulated gravitational-wave signals from a range of coalescing-binary locations and orientations, we study the improvement in parameter estimation. We find that an Australian detector can break degeneracies in several parameters; in particular, the localization of the source on the sky is improved by a factor of ~4, with more modest improvements in distance and binary inclination estimates. This enhanced ability to localize sources on the sky will be crucial in any search for electromagnetic counterparts to detected gravitational-wave signals., Comment: This paper has been withdrawn by the authors due to very extensive changes from the original version. Please see the paper "Estimating parameters of coalescing compact binaries with proposed advanced detector networks"

Galactic interacting double white dwarfs (DWD) are guaranteed gravitational wave (GW) sources for the GW detector LISA, with more than 10^4 binaries expected to be detected over the mission's lifetime. Part of this population is expected to be eccentric, and here we investigate the potential for constraining the white dwarf (WD) properties through apsidal precession in these binaries. We analyze the tidal, rotational, and general relativistic contributions to apsidal precession by using detailed He WD models, where the evolution of the star's interior is followed throughout the cooling phase. In agreement with previous studies of zero-temperature WDs, we find that apsidal precession in eccentric DWDs can lead to a detectable shift in the emitted GW signal when binaries with cool (old) components are considered. This shift increases significantly for hot (young) WDs. We find that apsidal motion in hot (cool) DWDs is dominated by tides at orbital frequencies above ~10^{-4}Hz (10^{- 3}$Hz). The analysis of apsidal precession in these sources while ignoring the tidal component would lead to an extreme bias in the mass determination, and could lead us to misidentify WDs as neutron stars or black holes. We use the detailed WD models to show that for older, cold WDs, there is a unique relationship that ties the radius and apsidal precession constant to the WD masses, therefore allowing tides to be used as a tool to constrain the source masses., Comment: 23 pages, 7 figures, revised to match accepted ApJ version

The system PSR J0737-3039 is the only binary pulsar known to consist of two radio pulsars (PSR J0737-3039 A and PSR J0737-3039 B). This unique configuration allows measurements of spin orientation for both pulsars: pulsar A's spin is tilted from the orbital angular momentum by no more than 14 degrees at 95% confidence; pulsar B's by 130 +/- 1 degrees at 99.7% confidence. This spin-spin misalignment requires that the origin of most of B's present-day spin is connected to the supernova that formed pulsar B. Under the simplified assumption of a single, instantaneous kick during the supernova, the spin could be thought of as originating from the off-center nature of the kick, causing pulsar B to tumble to its misaligned state. With this assumption, and using current constraints on the kick magnitude, we find that pulsar B's instantaneous kick must have been displaced from the center of mass of the exploding star by at least 1 km and probably 5--10 km. Regardless of the details of the kick mechanism and the process that produced pulsar B's current spin, the measured spin-spin misalignment in the double pulsar system provides an empirical, direct constraint on the angular momentum production in this supernova. This constraint can be used to guide core-collapse simulations and the quest for understanding the spins and kicks of compact objects., Comment: 12 pages, 3 figures, as accepted by ApJ

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., Comment: Manuscript: 18 pages, 2 tables, 2 figure. Supplementary Information: 34 pages, 6 figures. Advance Online Publication (AOP) on http://www.nature.com/nature on October 20, 2010. To Appear in Nature on November 4, 2010

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., Comment: 6 pages, 3 figures, to be published in the proceedings of the conference 'Binary star evolution: mass loss, accretion, and mergers', June 22-25, Mykonos, Greece (2010)

The common-envelope (CE) phase is an important stage in the evolution of binary stellar populations. The most common way to compute the change in orbital period during a CE is to relate the binding energy of the envelope of the Roche-lobe filling giant to the change in orbital energy. Especially in population-synthesis codes, where the evolution of millions of stars must be computed and detailed evolutionary models are too expensive computationally, simple approximations are made for the envelope binding energy. In this study, we present accurate analytic prescriptions based on detailed stellar-evolution models that provide the envelope binding energy for giants with metallicities between Z = 10-4 and Z = 0.03 and masses between 0.8 Msun and 100 Msun, as a function of the metallicity, mass, radius and evolutionary phase of the star. Our results are also presented in the form of electronic data tables and Fortran routines that use them. We find that the accuracy of our fits is better than 15% for 90% of our model data points in all cases, and better than 10% for 90% of our data points in all cases except the asymptotic giant branches for three of the six metallicities we consider. For very massive stars (M > 50 Msun), when stars lose more than ~20% of their initial mass due to stellar winds, our fits do not describe the models as accurately. Our results are more widely applicable - covering wider ranges of metallicity and mass - and are of higher accuracy than those of previous studies.

Motivated by suggestions that binaries with almost equal-mass components ("twins") play an important role in the formation of double neutron stars and may be rather abundant among binaries, we study the stability of synchronized close and contact binaries with identical components in circular orbits. In particular, we investigate the dependency of the innermost stable circular orbit on the core mass, and we study the coalescence of the binary that occurs at smaller separations. For twin binaries composed of convective main-sequence stars, subgiants, or giants with low mass cores (M_c <~0.15M, where M is the mass of a component), a secular instability is reached during the contact phase, accompanied by a dynamical mass transfer instability at the same or at a slightly smaller orbital separation. Binaries that come inside this instability limit transfer mass gradually from one component to the other and then coalesce quickly as mass is lost through the outer Lagrangian points. For twin giant binaries with moderate to massive cores (M_c >~0.15M), we find that stable contact configurations exist at all separations down to the Roche limit, when mass shedding through the outer Lagrangian points triggers a coalescence of the envelopes and leaves the cores orbiting in a central tight binary. In addition to the formation of binary neutron stars, we also discuss the implications of our results for the production of planetary nebulae with double degenerate central binaries., Comment: 17 pages, accepted to ApJ, final version includes discussion of planetary nebulae with central binaries and a new figure about shock heating, visualizations at http://webpub.allegheny.edu/employee/j/jalombar/movies/

We model a circular mass-transferring binary system to calculate the exchange of angular momentum between stellar spins and the orbit due to direct impact of the mass transfer stream onto the surface of the accretor. We simulate mass transfer by calculating the ballistic motion of a point mass ejected from the $L_1$ point of the donor star, conserving the total linear and angular momentum of the system, and treating the stars as uniform density spheres with main sequence radii determined by their masses. We show that, contrary to previous assumptions in the literature, direct impact does not always act as a sink of orbital angular momentum and may in fact increase it by facilitating the transfer of angular momentum from the spin of the donor to the orbit. Here, we show an example of the exchange of angular momentum, as well as a measure of the orbital angular momentum changes for a variety of binary star systems with main sequence components., Comment: 2 pages, 2 figures, Conference Proceedings for the International Conference on Binaries, Mykonos, Greece. Updated Version of Fig. 1b, correcting a scaling error. Results remain unchanged, but the numerical scaling factors have been decreased

Since the discovery of the first double neutron star (DNS) system in 1975 by Hulse and Taylor, there are currently 8 confirmed DNS in our galaxy. For every system, the masses of both neutron stars, the orbital semi- major axis and eccentricity are measured, and proper motion is known for half of the systems. Using the orbital parameters and kinematic information, if available, as constraints for all system, we investigate the immediate progenitor mass of the second-born neutron star and the magnitude of the supernova kick it received at birth, with the primary goal to understand the core collapse mechanism leading to neutron star formation. Compared to earlier studies, we use a novel method to address the uncertainty related to the unknown radial velocity of the observed systems. For PSR B1534+12 and PSR B1913+16, the kick magnitudes are 150 - 270 km/s and 190 - 450 km/s (with 95% confidence) respectively, and the progenitor masses of the 2nd born neutron stars are 1.3 - 3.4 Msun and 1.4 - 5.0 Msun (95%), respectively. These suggest that the 2nd born neutron star was formed by an iron core collapse supernova in both systems. For PSR J0737-3039, on the other hand, the kick magnitude is only 5 - 120 km/s (95%), and the progenitor mass of the 2nd born neutron star is 1.3 - 1.9 Msun (95%). Because of the relatively low progenitor mass and kick magnitude, the formation of the 2nd born neutron star in PSR J0737-3039 is potentially connected to an electron capture supernova of a massive O - Ne - Mg white dwarf. For the remaining 5 Galactic DNS, the kick magnitude ranges from several tens to several hundreds of km/s, and the progenitor mass of the 2nd formed neutron star can be as low as ~1.5 Msun, or as high as ~8 Msun. Therefore in these systems, it is not clear which type of supernova is more likely to form the 2nd neutron star., Comment: 13 pages, 16 figures, accepted for publication in ApJ

One of the primary goals when studying stellar systems with neutron stars has been to reveal the physical properties of progenitors and understand how neutron star spins and birth kicks are determined. Over the years a consensus understanding had been developed, but recently some of the basic elements of this understanding are being challenged by current observations of some binary systems and their theoretical interpretation. In what follows we review such recent developments and highlight how they are interconnected; we particularly emphasize some of the assumptions and caveats of theoretical interpretations and examine their validity (e.g., in connection to the unknown radial velocities of pulsars or the nuances of multi-dimensional statistical analysis). The emerging picture does not erase our earlier understanding; instead it broadens it as it reveals additional pathways for neutron star formation and evolution: neutron stars probably form at the end of both core collapse of Fe cores of massive stars and electron-capture supernovae of ONeMg cores of lower-mass stars; birth kicks are required to be high (well in excess of 100 km/s) for some neutron stars and low (< 100 km/s) for others depending on the formation process; and spin up may occur not just through Roche-lobe overflow but also through wind accretion or phases of hypercritical accretion during common envelope evolution., Comment: 9 pages,4 figures, proceedings paper for 40 Years of Pulsars Conference

Double black hole binaries are among the most important sources of gravitational radiation for ground-based detectors such as LIGO or VIRGO. Even if formed with lower efficiency than double neutron star binaries, they could dominate the predicted detection rates, since black holes are more massive than neutron stars and therefore could be detected at greater distances. Here we discuss an evolutionary process that can very significantly limit the formation of close double black hole binaries: the vast majority of their potential progenitors undergo a common envelope (CE) phase while the donor, one of the massive binary components, is evolving through the Hertzsprung gap. Our latest theoretical understanding of the CE process suggests that this will probably lead to a merger, inhibiting double black hole formation. Barring uncertainties in the physics of CE evolution, we use population synthesis calculations, and find that the corresponding reduction in the merger rate of double black holes formed in galactic fields is so great (by ~500) that their contribution to inspiral detection rates for ground-based detectors could become relatively small (1 in 10) compared to double neutron star binaries. A similar process also reduces the merger rates for double neutron stars, by factor of ~5, eliminating most of the previously predicted ultracompact NS-NS systems. Our predicted detection rates for Advanced LIGO are now much lower for double black holes (~2/yr), but are still quite high for double neutron stars (~20/yr). If double black holes were found to be dominant in the detected inspiral signals, this could indicate that they mainly originate from dense star clusters (not included here) or that our theoretical understanding of the CE phase requires significant revision., Comment: 13 pages, significant additions/changes (ApJ)

We study the process of mass segregation through 2-body relaxation in galactic nuclei with a central massive black hole (MBH). This study has bearing on a variety of astrophysical questions, from the distribution of X-ray binaries at the Galactic centre, to tidal disruptions of main-sequence and giant stars, to inspirals of compact objects into the MBH, an important category of events for the future space borne gravitational wave interferometer LISA. In relatively small galactic nuclei, typical hosts of MBHs with masses in the range 1e4-1e7 Msun, the relaxation induces the formation of a steep density cusp around the MBH and strong mass segregation. Using a spherical stellar dynamical Monte-Carlo code, we simulate the long-term relaxational evolution of galactic nucleus models with a spectrum of stellar masses. Our focus is the concentration of stellar black holes to the immediate vicinity of the MBH. Special attention is given to models developed to match the conditions in the Milky Way nucleus., Comment: 7 pages, 2 figures. Proceedings for the Galactic Centre Workshop 2006. Essentially a digest of the work reported in astro-ph/0603280

The study of how stars distribute themselves around a massive black hole (MBH) in the center of a galaxy is an important prerequisite for the understanding of many galactic-center processes. These include the observed overabundance of point X-ray sources at the Galactic center, the prediction of rates and characteristics of tidal disruptions of extended stars by the MBH and of inspirals of compact stars into the MBH, the latter being events of high importance for the future space borne gravitational wave interferometer LISA. In relatively small galactic nuclei, hosting MBHs with masses in the range 10^5-10^7 Msun, the single most important dynamical process is 2-body relaxation. It induces the formation of a steep density cusp around the MBH and strong mass segregation, as more massive stars lose energy to lighter ones and drift to the central regions. Using a spherical stellar dynamical Monte-Carlo code, we simulate the long-term relaxational evolution of galactic nucleus models with a spectrum of stellar masses. Our focus is the concentration of stellar black holes to the immediate vicinity of the MBH. We quantify this mass segregation for a variety of galactic nucleus models and discuss its astrophysical implications. Special attention is given to models developed to match the conditions in the Milky Way nucleus; we examine the presence of compact objects in connection to recent high-resolution X-ray observations., Comment: 28 pages, 24 figures, ApJ accepted. Small changes to follow referee's suggestions

We present a theoretical study of double compact objects as potential short/hard gamma-ray burst (GRB) progenitors. An updated population synthesis code StarTrack is used to calculate properties of double neutron stars and black-hole neutron star binaries. We obtain their formation rates, estimate merger times and finally predict their most likely merger locations and afterglow properties for different types of host galaxies. Our results serve for a direct comparison with the recent HETE-II and SWIFT observations of several short bursts, for which afterglows and host galaxies were detected. We also discuss the possible constraints these observations put on the evolutionary models of double compact object formation. We emphasize that our double compact object models can successfully reproduce at the same time short GRBs within both young, star-forming galaxies (e.g., GRB 050709 and GRB 051221A), as well as within old, elliptical hosts (e.g., GRB 050724 and probably GRB 050509B), Comment: ApJ, in press

Recent pulsar surveys have increased the number of observed double neutron stars (DNS) in our galaxy enough so that observable trends in their properties are starting to emerge. In particular, it has been noted that the majority of DNS have eccentricities less than 0.3, which are surprisingly low for binaries that survive a supernova explosion that we believe imparts a significant kick to the neutron star. To investigate this trend, we generate many different theoretical distributions of DNS eccentricities using Monte Carlo population synthesis methods. We determine which eccentricity distributions are most consistent with the observed sample of DNS binaries. In agreement with Chaurasia & Bailes (2005), assuming all double neutron stars are equally as probable to be discovered as binary pulsars, we find that highly eccentric, coalescing DNS are less likely to be observed because of their accelerated orbital evolution due to gravitational wave emission and possible early mergers. Based on our results for coalescing DNS, we also find that models with vanishingly or moderately small kicks (sigma < about 50 km/s) are inconsistent with the current observed sample of such DNS. We discuss the implications of our conclusions for DNS merger rate estimates of interest to ground-based gravitational-wave interferometers. We find that, although orbital evolution due to gravitational radiation affects the eccentricity distribution of the observed sample, the associated upwards correction factor to merger rate estimates is rather small (typically 10-40%)., Comment: 9 pages, 8 figures, accepted by ApJ. Figures reduced and some content changed, references added

In this paper we consider the population of eccentric binaries with a neutron star and a white dwarf that has been revealed in our galaxy in recent years through binary pulsar observations. We apply our statistical analysis method (Kim, Kalogera, & Lorimer 2003)and calculate the Galactic formation rate of these binaries empirically. We then compare our results with rate predictions based on binary population synthesis from various research groups and for various ranges of model input parameters. For our reference moel, we find the Galactic formation rate of these eccentric systems to be ~7 per Myr, about an order of magnitude smaller than results from binary evolution estimations. However, the empirical estimates are calculated with no correction for pulsar beaming, and therefore they should be taken as lower limits. Despite uncertainties that exceed an order of magnitude, there is significant overlap of the various rate calculations. This consistency lends confidence that our current understanding of the formation of these eccentric NS-WD binaries is reasonable., Comment: 7 pages, 2 figures, 2 tables, to appear in proceedings of Aspen Center for Physics Conference on ``Binary Radio Pulsars'' Eds. F. Rasio and I. Stairs

We summarize our results on the Galactic merger rate of double neutron stars (DNS) in view of the recent discovery of PSR J0737-3039. We also present previously unpublished results for the 'global' probability distribution of merger rate values that incorporate the presently known systematics from the radio pulsar luminosity function. The most likely value obtained from the global distribution is only ~15 per Myr, but a re-analysis of the current pulsar sample and radio luminosities is needed for a reliable assessment of the best fitting distribution. Finally, we use our theoretical understanding of DNS formation to calculate a possible upper limit on the DNS merger rate from current Type Ib/c supernova rate estimates., Comment: 8 pages, 3 figures, to appear in proceedings of Aspen Center for Physics Conference on ``Binary Radio Pulsars'' Eds. F. Rasio and I. Stairs. New results on the systematic uncertainties and possible upper limits on rate estimates are included; Rate calculation is revised

We consider the statistics of pulsar binaries with white dwarf companions (NS-WD). Using the statistical analysis method developed by Kim et al. (2003) we calculate the Galactic coalescence rate of NS-WD binaries due to gravitational-wave emission. We find that the most likely values for the total Galactic coalescence rate (R_tot) of NS-WD binaries lie in the range 0.2--10 per Myr depending on different assumed pulsar population models. For our reference model, we obtain R_tot=4.11_(-2.56)^(+5.25) per Myr at a 68% statistical confidence level. These rate estimates are not corrected for pulsar beaming and as such they are found to be about a factor of 20 smaller than the Galactic coalescence rate estimates for double neutron star systems. Based on our rate estimates, we calculate the gravitational-wave background due to coalescing NS-WD binaries out to extragalactic distances within the frequency band of the Laser Interferometer Space Antenna. We find the contribution from NS-WD binaries to the gravitational-wave background to be negligible., Comment: 20 pages, 2 figures, 2 tables, Accepted for publication in ApJ

Relativistic spin-orbit and spin-spin couplings has been shown to modify the gravitational waveforms expected from inspiraling binaries with a black hole and a neutron star. As a result inspiral signals may be missed due to significant losses in signal-to-noise ratio, if precession effects are ignored in gravitational-wave searches. We examine the sensitivity of the anticipated loss of signal-to-noise ratio on two factors: the accuracy of the precessing waveforms adopted as the true signals and the expected distributions of spin-orbit tilt angles, given the current understanding of their physical origin. We find that the results obtained using signals generated by approximate techniques are in good agreement with the ones obtained by integrating the 2PN equations. This shows that a complete account of all high-order post-Newtonian effects is usually not necessary for the determination of detection efficiencies. Based on our current astrophysical expectations, large tilt angles are not favored and as a result the decrease in detection rate varies rather slowly with respect to the black hole spin magnitude and is within 20--30% of the maximum possible values., Comment: 7 fig., accepted by Phys. Rev. D Minor modifications

We consider the implications of the presence of ~1 stellar-mass black hole (BH) at the center of a dense globular cluster. We show that BH X-ray binaries formed through exchange interactions are likely to have extremely low duty cycles (<0.001), consistent with the absence of BH X-ray transients in Galactic globular clusters. In contrast, we find that BH X-ray binaries formed through tidal capture would be persistent, bright X-ray sources. Given the absence of any such source and the very high interaction rates, we conclude that tidal capture of a main-sequence star by a BH most likely leads to the complete disruption of the star., Comment: 4 pages, to appear in ApJ Letters

We present an analysis method that allows us to estimate the Galactic formation of radio pulsar populations based on their observed properties and our understanding of survey selection effects. More importantly, this method allows us to assign a statistical significance to such rate estimates and calculate the allowed ranges of values at various confidence levels. Here, we apply the method to the question of the double neutron star (NS-NS) coalescence rate using the current observed sample, and we find calculate the most likely value for the total Galactic coalescence rate to lie in the range 3-22 Myr^{-1}, for different pulsar population models. The corresponding range of expected detection rates of NS--NS inspiral are (1-9)x10^{-3} yr^{-1} for the initial LIGO, and 6-50 yr^{-1} for the advanced LIGO. Based on this newly developed statistical method, we also calculate the probability distribution for the expected number of pulsars that could be observed by the Parkes Multibeam survey, when acceleration searches will alleviate the effects of Doppler smearing due to orbital motions. We suggest that the Parkes survey will probably detect 1-2 new binary pulsars like PSRs B1913+16 and/or B1534+12., Comment: 8 pages; to appear in "Radio Pulsars" (ASP Conf. Ser.), eds. M. Bailes, D. Nice, & S. Thorsett

We examine the theoretical implications of a population of low-mass helium-core white dwarfs in globular clusters. In particular, we focus on the observed population in the core of NGC 6397, where several low-mass white dwarf canditates have been identified as ``non-flickerers'' by Cool and collaborators. Age and mass estimates from cooling models, combined with dynamical and evolutionary considerations, lead us to infer that the dark binary companions are C/O white dwarfs rather than neutron stars. Furthermore, we find that the progenitor binaries very likely underwent an exchange interaction within the last 1.e9 years. We examine the prospects for detecting a similar population in other globular clusters, with particular attention to the case of 47 Tuc., Comment: 30 pages, 9 postscript figures, submitted to ApJ