Recent studies suggest that globular clusters (GCs) may retain a substantial population of stellar-mass black holes (BHs), in contrast to the long-held belief of a few to zero BHs. We model the population of BH low-mass X-ray binaries (BH-LMXBs), an ideal observable proxy for elusive single BHs, produced from a representative group of Milky Way GCs with variable BH populations. We simulate the formation of BH-binaries in GCs through exchange interactions between binary and single stars in the company of tens to hundreds of BHs. Additionally, we consider the impact of the BH population on the rate of compact binaries undergoing gravitational wave driven mergers. The characteristics of the BH-LMXB population and binary properties are sensitive to the GCs structural parameters as well as its unobservable BH population. We find that GCs retaining $\sim 1000$ BHs produce a galactic population of $\sim 150$ ejected BH-LMXBs whereas GCs retaining only $\sim20$ BHs produce zero ejected BH-LMXBs. Moreover, we explore the possibility that some of the presently known BH-LMXBs might have originated in GCs and identify five candidate systems.

The recent LIGO detection of gravitational waves (GW150914), likely originating from the merger of two $\sim 30 M_\odot$ black holes suggests progenitor stars of low metallicity ($[Z/Z_\odot] \lesssim 0.3$), constraining when and where the progenitor of GW150914 may have formed. We combine estimates of galaxy properties (metallicity, star formation rate and merger rate) across cosmic time to predict the low redshift black hole - black hole merger rate as a function of present day host galaxy mass, $M_\mathrm{gal}$, and the formation redshift of the progenitor system $z_\mathrm{form}$ for different progenitor metallicities $Z_\mathrm{c}$. At $Z_\mathrm{c}=0.1 Z_\odot$, the signal is dominated by binaries in massive galaxies with $z_\mathrm{form}\simeq 2$, with a small contribution from binaries formed around $z_\mathrm{form}\simeq 0.5$ in dwarf galaxies. For $Z_\mathrm{c}=0.01Z_\odot$, fast mergers are possible and very recent star formation in dwarfs likely dominates. Additional gravitational wave detections from merging massive black holes will provide constraints on the mass-metallicity relation and massive star formation at high redshifts.

We present the SuperNova Explosion Code (SNEC), an open-source Lagrangian code for the hydrodynamics and equilibrium-diffusion radiation transport in the expanding envelopes of supernovae. Given a model of a progenitor star, an explosion energy, and an amount and distribution of radioactive nickel, SNEC generates the bolometric light curve, as well as the light curves in different broad bands assuming black body emission. As a first application of SNEC, we consider the explosions of a grid of 15 Msun (at zero-age main sequence) stars whose hydrogen envelopes are stripped to different extents and at different points in their evolution. The resulting light curves exhibit plateaus with durations of ~20-100 days if >~1.5-2 Msun of hydrogen-rich material is left and no plateau if less hydrogen-rich material is left. If these shorter plateau lengths are not seen for Type IIP supernovae in nature, it suggests that, at least for zero-age main sequence masses <~ 20 Msun, hydrogen mass loss occurs as an all or nothing process. This perhaps points to the important role binary interactions play in generating the observed mass-stripped supernovae (i.e., Type Ib/c events). These light curves are also unlike what is typically seen for Type IIL supernovae, arguing that simply varying the amount of mass loss cannot explain these events. The most stripped models begin to show double-peaked light curves similar to what is often seen for Type IIb supernovae, confirming previous work that these supernovae can come from progenitors that have a small amount of hydrogen and a radius of ~500 Rsun.

A longstanding question in stellar evolution is which massive stars produce black holes (BHs) rather than neutron stars (NSs) upon death. It has been common practice to assume that a given zero-age main sequence (ZAMS) mass star (and perhaps a given metallicity) simply produces either an NS or a BH, but this fails to account for a myriad of other variables that may effect this outcome, such as spin, binarity, or even stochastic differences in the stellar structure near core collapse. We argue that instead a probabilistic description of NS versus BH formation may be better suited to account for the current uncertainties in understanding how massive stars die. We present an initial exploration of the probability that a star will make a BH as a function of its ZAMS mass, $P_{\rm BH}(M_{\rm ZAMS})$. Although we find that it is difficult to derive a unique $P_{\rm BH}(M_{\rm ZAMS})$ using current measurements of both the BH mass distribution and the degree of chemical enrichment by massive stars, we demonstrate how $P_{\rm BH}(M_{\rm ZAMS})$ changes with these various observational and theoretical uncertainties. We anticipate that future studies of Galactic BHs and theoretical studies of core collapse will refine $P_{\rm BH}(M_{\rm ZAMS})$ and argue that this framework is an important new step toward better understanding BH formation. A probabilistic description of BH formation will be useful as input for future population synthesis studies that are interested in the formation of X-ray binaries, the nature and event rate of gravitational wave sources, and answering questions about chemical enrichment.

The discovery of a binary comprising a black hole (BH) and a millisecond pulsar (MSP) would yield insights into stellar evolution and facilitate exquisitely sensitive tests of general relativity. Globular clusters (GCs) are known to harbor large MSP populations and recent studies suggest that GCs may also retain a substantial population of stellar mass BHs. We modeled the formation of BH+MSP binaries in GCs through exchange interactions between binary and single stars. We found that in dense, massive clusters most of the dynamically formed BH+MSP binaries will have orbital periods of 2 to 10 days, regardless of the mass of the BH, the number of BHs retained by the cluster, and the nature of the GC's binary population. The size of the BH+MSP population is sensitive to several uncertain parameters, including the BH mass function, the BH retention fraction, and the binary fraction in GCs. Based on our models, we estimate that there are $0.6\pm0.2$ dynamically formed BH+MSP binaries in the Milky Way GC system, and place an upper limit on the size of this population of $\sim 10$. Interestingly, we find that BH+MSP binaries will be rare even if GCs retain large BH populations.

We model the formation of black hole-neutron star (BH-NS) binaries via dynamical interactions in globular clusters. We find that in dense, massive clusters, 16-61% of the BH-NS binaries formed by interactions with existing BH binaries will undergo mergers driven by the emission of gravitational radiation. If the BHs are retained by the cluster after merging with a NS, the BHs acquire subsequent NS companions and undergo several mergers. Thus, the merger rate depends critically upon whether or not the BH is retained by the cluster after the merger. Results from numerical relativity suggest that kick imparted to a ~7 M_sun BH after it merges with a NS will greatly exceed the cluster's escape velocity. In this case, the models suggest that the majority of BH-NS mergers in globular clusters occur within 4 Gyrs of the cluster's formation and would be unobservable by Advanced LIGO. For more massive BHs, on the other hand, the post merger kick is suppressed and the BH is retained. Models with 35 M_sun BHs predict Advanced LIGO detection rates in the range 0.04 - 0.7 per year. On the pessimistic end of this range, BH-NS mergers resulting from binary-single star interactions in globular clusters could account for an interesting fraction of all BH-NS mergers. On the optimistic end, this channel may dominate the rate of detectable BH-NS mergers.

We model the emission lines generated in the photoionised debris of a tidally disrupted horizontal branch star. We find that at late times, the brightest optical emission lines are [N II] \lambda\lambda 6548,6583 and [O III] \lambda\lambda 4959,5007. Models of a red clump horizontal branch star undergoing mild disruption by a massive (50 -- 100 M_\sun) black hole yield an emission line spectrum that is in good agreement with that observed in the NGC 1399 globular cluster hosting the ultraluminous X-ray source CXOJ033831.8 - 352604. We make predictions for the UV emission line spectrum that can verify the tidal disruption scenario and constrain the mass of the BH.

Binaries that contain a hot subdwarf (sdB) star and a main sequence companion may have interacted in the past. This binary population has historically helped determine our understanding of binary stellar evolution. We have computed a grid of binary population synthesis models using different assumptions about the minimum core mass for helium ignition, the envelope binding energy, the common envelope ejection efficiency, the amount of mass and angular momentum lost during stable mass transfer, and the criteria for stable mass transfer on the red giant branch and in the Hertzsprung gap. These parameters separately and together can significantly change the entire predicted population of sdBs. Nonetheless, several different parameter sets can reproduce the observed subpopulation of sdB + white dwarf and sdB + M dwarf binaries, which has been used to constrain these parameters in previous studies. The period distribution of sdB + early F dwarf binaries offers a better test of different mass transfer scenarios for stars that fill their Roche lobes on the red giant branch.

Many hot subdwarf B stars (sdBs) are in close binaries, and the favored formation channels for subdwarfs rely on mass transfer in a binary system to strip a core He burning star of its envelope. However, these channels cannot account for sdBs that have been observed in long period binaries nor the narrow mass distribution of isolated (or "singleton") sdBs. We propose a new formation channel involving the merger of a helium white dwarf and a low mass, hydrogen burning star, which addresses these issues. Hierarchical triples whose inner binaries merge and form sdBs by this process could explain the observed long period subdwarf+main sequence binaries. This process would also naturally explain the observed slow rotational speeds of singleton sdBs. We also briefly discuss the implications of this formation channel for extreme horizontal branch morphology in globular clusters and the UV upturn in elliptical galaxies.

We calculate the emission line spectrum produced by the debris released when a white dwarf (WD) is tidally disrupted by an intermediate-mass black hole (IMBH; $M\sim 10^{2}-10^{5}\msun$) and we explore the possibility of using the emission lines to identify such events and constrain the properties of the IMBH. To this end, we adopt and adapt the techniques developed by Strubbe & Quataert to study the optical emission lines produced when a main sequence (MS) star is tidally disrupted by a supermassive black hole. WDs are tidally disrupted outside of the event horizon of a $< 10^{5}\msun$ black hole, which makes these tidal disruption events good signposts of IMBHs. We focus on the optical and UV emission lines produced when the accretion flare photoionizes the stream of debris that remains unbound during the disruption. We find that the spectrum is dominated by lines due to ions of C and O, the strongest of which are \ion{C}{4} $\lambda$1549 at early times and [\ion{O}{3}] $\lambda$5007 at later times. Furthermore, we model the profile of the emission lines in the [\ion{O}{3}] $\lambda\lambda$4959, 5007 doublet and find that it is highly asymmetric with velocity widths of up to $\sim 2500 \rm{\;km\;s^{-1}}$, depending on the properties of the WD-IMBH system and the orientation of the observer. Finally, we compare the models with observations of X-ray flares and optical emission lines in the cores of globular clusters and propose how future observations can test if these features are due to a WD that has been tidally disrupted by an IMBH.