Recent numerical simulations of globular clusters (GCs) have shown that stellar-mass black holes (BHs) play a fundamental role in driving cluster evolution and shaping their present-day structure. Rapidly mass-segregating to the center of GCs, BHs act as a dynamical energy source via repeated super-elastic scattering, delaying onset of core collapse and limiting mass segregation for visible stars. While recent discoveries of BH candidates in Galactic and extragalactic GCs have further piqued interest in BH-mediated cluster dynamics, numerical models show that even if significant BH populations remain in today's GCs, they are typically in configurations that are not directly detectable. We demonstrated in Weatherford et al. (2018) that an anti-correlation between a suitable measure of mass segregation ($\Delta$) in observable stellar populations and the number of retained BHs in GC models can be applied to indirectly probe BH populations in real GCs. Here, we estimate the number and total mass of BHs presently retained in 50 Milky Way GCs from the ACS Globular Cluster Survey by measuring $\Delta$ between populations of main sequence stars, using correlations found between $\Delta$ and BH retention in the CMC Cluster Catalog models. We demonstrate that the range in $\Delta$'s distribution from our models matches that for observed GCs to a remarkable degree. Our results further provide the narrowest constraints to-date on the retained BH populations in the GCs analyzed. Of these 50 GCs, we identify NGCs 2808, 5927, 5986, 6101, and 6205 to presently contain especially large BH populations, each with total BH mass exceeding $10^3\,\rm{M_{\odot}}$.

Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae. However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called "upper-mass gap" would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through pre-collapse stellar collisions in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as $\sim20\%$ of all BH progenitors undergo one or more collisions prior to stellar collapse and up to $\sim1\%$ of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the "multiple-generation" merger channel described in other analyses. Finally, we describe how stellar collisions in clusters may provide a unique pathway to pair-instability supernovae and briefly discuss the expected rate of these events and other electromagnetic transients.

We study the formation rate of binary black hole mergers formed through gravitational-wave emission between unbound, single black holes in globular clusters. While the formation of these binaries in very dense systems such as galactic nuclei has been well studied, we show here that this process can operate in lower-density stellar systems as well, forming binaries at a rate similar to other proposed pathways for creating eccentric mergers. Recent advances in post-Newtonian cluster dynamics indicate that a large fraction of dynamically-assembled binary black holes merge inside their host clusters during weak and strong binary-single and binary-binary interactions, and that these systems may retain measurable eccentricities as they travel through the LIGO and LISA sensitivity bands. Using an analytic approach to modeling binary black holes from globular clusters, we show that the formation of merging binaries from previously unbound black holes can operate at a similar rate to mergers forming during strong binary encounters, and that these binaries inhabit a unique region of the gravitational-wave frequency space which can be identified by proposed deci-Hertz space-based detectors.

We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics.

Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt, despite their potential role as massive seeds for quasars and sources of tidal disruption events, ultra-luminous X-ray sources, dwarf galaxy feedback, and hypervelocity stars. Gravitational wave (GW) observations can help to find and confirm the existence of IMBHs. Current and upcoming detectors, such as LIGO, Virgo, KAGRA, LISA, ET, and DECIGO promise to identify the full range from stellar-mass to supermassive black holes (SMBHs). In this paper, we address the question of whether IMBHs can produce GWs in galactic nuclei. We consider the possibility that stellar black holes (SBHs) form bound systems and later coalesce with an IMBH through gravitational captures in the dense nucleus. We show that this mechanism is efficient for IMBH masses in the range $\sim 3\times 10^3$ M$_\odot$--$2\times 10^4$ M$_\odot$. We find that the typical distributions of peak frequencies and merger timescales depend mainly on the IMBH mass. In particular, the typical peak frequency is about $0.2\,$Hz, $0.1\,$Hz, $0.09\,$Hz, and $0.05\,$Hz for $M_{\rm IMBH}=5\times 10^3$ M$_\odot$, $8\times 10^3$ M$_\odot$, $1\times 10^4$ M$_\odot$, and $2\times 10^4$ M$_\odot$, respectively. Our results show that, at design sensitivity, both DECIGO and ET should be able to detect these IMBH--SBH mergers. Furthermore, most of the mergers will appear eccentric ($e \gtrsim 0.1$), providing an indication of their dynamical origin.

We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters.

There is increasing evidence that globular clusters retain sizeable black hole populations at present day. This is supported by dynamical simulations of cluster evolution, which have unveiled the spatial distribution and mass spectrum of black holes in clusters across cosmic age. However, the black hole populations of young, high metallicity clusters remain unconstrained. The black holes hosted by these clusters are expected to rapidly mass segregate early in their evolutionary history, forming central cusps of hundreds to thousands of black holes. Here, we argue that the host young cluster can accumulate gas from its dense surroundings, from which the black hole cusp members can accrete at highly enhanced rates. Their collective accretion luminosity can be substantial and provides a novel observational constraint for nearby young massive clusters. We test this hypothesis by performing 3D hydrodynamic simulations where we embed discretized potentials, representing our black holes, within the core potential of a massive cluster. This system moves supersonically with respect to a gaseous medium from which it accretes. We study the accretion of this black hole cusp for different cusp populations and determine the integrated accretion luminosity of the black hole cusp. We apply our results to the young massive clusters of the Antennae Galaxies and find that a typical cusp accretion luminosity should be in excess of $\approx 10^{40}\,{\rm ergs\,\,s^{-1}}$. We argue that no strong candidates of this luminous signal have been observed, and constrain the cusp population of a typical cluster in the Antennae Galaxies to $\lesssim10-2\times10^2$ $10\,M_\odot$ black holes. This is in conflict with predictions of dynamical simulations, which suggest that these systems should harbor several hundred to thousands of black holes, and we explore the reasons for this discrepancy in our conclusions.

The first detection of gravitational waves from a neutron star - neutron star (NS-NS) merger, GW170817, and the increasing number of observations of short gamma-ray bursts (SGRBs) have greatly motivated studies of the origins of NS-NS and neutron star - black hole (NS-BH) binaries. We calculate the merger rates of NS-NS and NS-BH binaries from globular clusters (GCs) using realistic GC simulations with the \texttt{CMC} cluster catalog. We use a large sample of models with a range of initial numbers of stars, metallicities, virial radii and galactocentric distances, representative of the present-day Milky Way GCs, to quantify the inspiral times and volumetric merger rates as a function of redshift, both inside and ejected from clusters. We find that over the complete lifetime of most GCs, stellar BHs dominate the cluster cores and prevent the mass segregation of NSs, thereby reducing the dynamical interaction rates of NSs so that at most a few NS binary mergers are ever produced. We estimate the merger rate in the local universe to be $\sim\rm{0.02\,Gpc^{-3}\,yr^{-1}}$ for both NS-NS and NS-BH binaries, or a total of $\sim 0.04$~Gpc$^{-3}$~yr$^{-1}$ for both populations. These rates are about 5 orders of magnitude below the current empirical merger rate from LIGO/Virgo. We conclude that dynamical interactions in GCs do not play a significant role in enhancing the NS-NS and NS-BH merger rates.

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

Star-to-star dispersion of r-process elements has been observed in a significant number of old, metal-poor globular clusters. We investigate early-time neutron-star mergers as the mechanism for this enrichment. Through both numerical modeling and analytical arguments, we show that neutron-star mergers cannot be induced through dynamical interactions early in the history of the cluster, even when the most liberal assumptions about neutron-star segregation are assumed. Therefore, if neutron-star mergers are the primary mechanism for r-process dispersion in globular clusters, they likely result from the evolution of isolated, primordial binaries in the clusters. Through population modeling, we find that moderate fractions of GCs with enrichment are only possible when a significant number of double neutron-star progenitors proceed through Case BB mass transfer --- under various assumptions for the initial properties of globular clusters, a neutron-star merger with the potential for enrichment will occur in ~15-60% (~30-90%) of globular clusters if this mass transfer proceeds stably (unstably). The strong anti-correlation between the pre-supernova orbital separation and post-supernova systemic velocity due to mass loss in the supernova leads to efficient ejection of most enrichment candidates from their host clusters. Thus, most enrichment events occur shortly after the double neutron stars are born. This requires star-forming gas that can absorb the r-process ejecta to be present in the globular cluster 30-50 Myr after the initial burst of star formation. If scenarios for redistributing gas in globular clusters cannot act on these timescales, the number of neutron-star merger enrichment candidates drops severely, and it is likely that another mechanism, such as r-process enrichment from collapsars, is at play.

Among the growing list of confirmed exoplanets, the number of planets identified in dense star clusters remains sparse. Previous analyses have suggested this may be due in part to dynamical interactions that unbind large fractions of planets from their host stars, limiting the survival of planetary systems in clusters. Thus, alternative detection strategies may be necessary to study planets in clusters that may no longer be bound to a host star. Here, we use the cluster Monte Carlo code CMC to explore the evolution of planetary systems in dense star clusters. Depending on a number of initial conditions, we show that $10-50\%$ of primordial planetary systems are broken through dynamical encounters over a cluster's full lifetime, populating clusters with "free-floating" planets. Furthermore, a large number ($30-80\%$) of planets are ejected from their host cluster through strong dynamical encounters and/or tidal loss. Additionally, we show that planets naturally mix with stellar-mass black holes (BHs) in the central regions of their host cluster. As a consequence, up to a few hundreds of planets will be tidally disrupted through close passages of BHs. We show these BH-planet tidal disruption events (TDEs) occur in clusters at a rate of up to $10^{-5}\,\rm{yr}^{-1}$ in a Milky Way-type galaxy. We predict BH-planet TDEs should be detected by upcoming transient surveys such as LSST at a rate of roughly a few events per year. The observed rate of BH-planet TDEs would place new constraints upon the formation and survival of both planetary systems and BHs in dense star clusters. Additionally, depending on various assumptions including the initial number of planets and their orbital properties, we predict that typical globular clusters may contain a few dynamically-formed NS-planet systems at present as well as up to roughly 100 dynamically-formed WD-planet systems.

When two black holes merge in a dense star cluster, they form a new black hole with a well-defined mass and spin. If that "second-generation" black hole remains in the cluster, it will continue to participate in dynamical encounters, form binaries, and potentially merge again. Using a grid of 96 dynamical models of dense star clusters and a cosmological model of cluster formation, we explore the production of binary black hole mergers where at least one component of the binary was forged in a previous merger. We create four hypothetical universes where every black hole born in the collapse of a massive star has a dimensionless Kerr spin parameter, $\chi_{\rm birth}$, of 0.0, 0.1, 0.2, or 0.5. We show that if all stellar-born black holes are non-spinning ($\chi_{\rm birth} = 0.0$), then more than 10% of merging binary black holes from clusters have components formed from previous mergers, accounting for more than 20% of the mergers from globular clusters detectable by LIGO/Virgo. Furthermore, nearly 7% of detectable mergers would have a component with a mass $\gtrsim 55M_{\odot}$, placing them clearly in the mass "gap" region where black holes cannot form from isolated collapsing stars due to the pulsational-pair instability mechanism. On the other hand, if black holes are born spinning, then the contribution from these second-generation mergers decreases, making up as little as 1% of all detections from globular clusters when $\chi_{\rm birth} = 0.5$. We make quantitative predictions for the detected masses, mass ratios, and spin properties of first- and second-generation mergers from dense star clusters, and show how these distributions are highly sensitive to the birth spins of black holes.

As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, which leads to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward so-called "core collapse," the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the core-collapsed and non-core-collapsed clusters. This suggests an internal energy source is at work, delaying the onset of core collapse in many clusters. Primordial binary stars have been thought for a long time to provide this energy source, but recent analyses have cast doubt upon the corresponding "binary-burning" mechanism as a viable explanation. Over the past decade, a large amount of both observational and theoretical work has suggested that many stellar-mass black holes (BHs) are retained in typical clusters today and that they play a dynamically-significant role in these clusters throughout their entire lifetimes. Here we review our latest understanding of the formation and evolution of BH populations in GCs and demonstrate that, through their dynamical interaction with their host cluster, BHs can naturally explain the distinction between core-collapsed and non-core-collapsed clusters through a process we call "black hole burning."

In a dense stellar environment, such as the core of a globular cluster (GC), dynamical interactions with black holes (BHs) are expected to lead to a variety of astrophysical transients. Here we explore tidal disruption events (TDEs) of stars by stellar-mass BHs through collisions and close encounters. Using state-of-the-art $N$-body simulations, we show that these TDEs occur at significant rates throughout the evolution of typical GCs and we study how their relative rates relate to cluster parameters such mass and size. By incorporating a realistic cosmological model of GC formation, we predict a BH - main-sequence-star TDE rate of approximately $3\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ in the local universe ($z<0.1$) and a cosmological rate that peaks at roughly $25\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ for redshift 3. Furthermore, we show that the ejected mass associated with these TDEs could produce optical transients of luminosity $\sim 10^{41} - 10^{44}\rm\,erg\,s^{-1}$ with timescales of about a day to a month. These should be readily detectable by optical transient surveys such as the Zwicky Transient Facility. Finally, we comment briefly on BH - giant encounters and discuss how these events may contribute to the formation of BH - white-dwarf binaries.

Over a hundred millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. Here we study MSPs in GCs using realistic $N$-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae (including both accretion-induced and merger-induced collapse of white dwarfs) can be spun up through mass transfer to form MSPs. Both NS formation and spin-up through accretion are greatly enhanced through dynamical interaction processes. We find that our models for average GCs at the present day with masses $\approx 2 \times 10^5\,M_\odot$ can produce up to $10-20$ MSPs, while a very massive GC model with mass $\approx 10^6\,M_\odot$ can produce close to $100$. We show that the number of MSPs is anti-correlated with the total number of stellar-mass black holes (BHs) retained in the host cluster. The radial distributions are also affected: MSPs are more concentrated towards the center in a host cluster with a smaller number of retained BHs. As a result, the number of MSPs in a GC could be used to place constraints on its BH population. Some intrinsic properties of MSP systems in our models (such as the magnetic fields and spin periods) are in good overall agreement with observations, while others (such as the distribution of binary companion types) less so, and we discuss the possible reasons for such discrepancies. Interestingly, our models also demonstrate the possibility of dynamically forming NS--NS and NS--BH binaries in GCs, although the predicted numbers are very small.

Ultra-compact binaries (UCBs) are systems containing compact or degenerate stars with orbital periods less than one hour. Tens of millions of UCBs are predicted to exist within theGalaxy emitting gravitational waves (GWs) at mHz frequencies. Combining GW searches with electromagnetic (EM) surveys like Gaia and LSST will yield a comprehensive, multimessenger catalog of UCBs in the galaxy. Joint EM and GW observations enable measurements of masses, radii, and orbital dynamics far beyond what can be achieved by independent EM or GW studies. GW+EM surveys of UCBs in the galaxy will yield a trove of unique insights into the nature of white dwarfs, the formation of compact objects, dynamical interactions in binaries, and energetic, accretion-driven phenomena like Type Ia superonovae.

The dynamical processing of black holes in the dense cores of globular clusters (GCs), makes them efficient factories for producing binary black holes (BBHs). Here we explore the population of BBHs that form dynamically in GCs and may be observable at mHz frequencies or higher with LISA. We use our Monte Carlo stellar dynamics code, which includes gravitational radiation reaction effects for all BH encounters. By creating a representative local universe of GCs, we show that up to dozens of these systems may be resolvable by LISA with signal-to-noise ratios of at least 5. Approximately one third of these binaries will have measurable eccentricities ($e > 10^{-3}$) in the LISA band and a small number ($\lesssim 5$) may evolve from the LISA band to the LIGO band during the LISA mission.

About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary stars. However, observations suggest that the eccentricity distributions of most observed binaries, and particularly those with masses < 5 Msun, are flatter than thermal and follow more closely to a uniform distribution. Nonetheless, it is often argued that dynamical interactions in a star cluster would quickly thermalize the binaries, which could justify imposing a thermal eccentricity distribution at birth for all binaries. In this paper we investigate the validity of this assumption. We develop our own rapid semi-analytic model for binary evolution in star clusters, and also compare with detailed N-body and Monte Carlo star cluster models. We show that, for nearly all binaries, dynamical encounters fail to convert an initially uniform eccentricity distribution to thermal within a star cluster's lifetime. Thus, if a thermal eccentricity distribution is observed, it is likely imprinted upon formation rather than through subsequent long-term dynamical processing. Theoretical investigations that initialize all binaries with a thermal distribution will make incorrect predictions for the evolution of the binary population. Such models may overpredict the merger rate for binaries with modest orbital separations by a factor of about two.

The dynamical assembly of binary black holes (BBHs) in dense star clusters (SCs) is one of the most promising pathways for producing observable gravitational wave (GW) sources, however several other formation scenarios likely operate as well. One of the current outstanding questions is how these different pathways may be distinguished apart. In this paper we suggest a new multi-messenger observable that can be used to constrain the formation of BBH mergers originating from SCs: the electromagnetic signal from tidal disruptions (TDs) of stars by BBHs. Such TDs will show variability in their light curve from the orbital motion of the disruptive BBHs, and can therefore be used to map the BBH orbital period distribution, and thereby also the dynamical mechanisms that eventually drive the BBHs to merger. Using an analytical approach including General Relativistic effects, we find that the orbital period distribution of BBHs within globular clusters peaks on timescales of days, which we argue is unique to this assembly pathway. We propose that the search for variable TDs in current and future EM transient surveys might be used to constrain the merger history of BBHs in SCs.

Globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the so-called "core-collapsed" and "non-core-collapsed" clusters. Here, we use our H\'enon-type Monte Carlo code, CMC, to explore initial cluster parameters that map into this bimodality. Remarkably, we find that by varying the initial size of clusters (specified in our initial conditions in terms of the initial virial radius, $r_v$) within a relatively narrow range consistent with the measured radii of young star clusters in the local universe ($r_v \approx 0.5-5$ pc), our models reproduce the variety of present-day cluster properties. Furthermore, we show that stellar-mass black holes (BHs) play an intimate role in this mapping from initial conditions to the present-day structural features of GCs. We identify "best-fit" models for three GCs with known observed BH candidates, NGC 3201, M22, and M10, and show that these clusters harbor populations of $\sim 50-100$ stellar-mass BHs at present. As an alternative case, we also compare our models to the core-collapsed cluster NGC 6752 and show that this cluster likely contains few BHs at present. Additionally, we explore the formation of BH binaries in GCs and demonstrate that these systems form naturally in our models in both detached and mass-transferring configurations with a variety of companion stellar types, including low-mass main sequence stars, white dwarfs, and sub-subgiants.

Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic population of dynamically-formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission between unbound black holes during chaotic resonant encounters, with roughly half of those events having eccentricities detectable by current ground-based gravitational-wave detectors. The mergers that occur inside clusters typically have lower masses than binaries that were ejected from the cluster many Gyrs ago. Gravitational-wave captures from globular clusters contribute 1-2 Gpc^-3 yr^-1 to the binary merger rate in the local universe, increasing to ~10 Gpc^-3 yr^-1 at z~3. Finally, we discuss some of the technical difficulties associated with post-Newtonian scattering encounters, and how care must be taken when measuring the binary parameters during a dynamical capture.

We explore the formation of double-compact-object binaries in Milky Way (MW) globular clusters (GCs) that may be detectable by the Laser Interferometer Space Antenna (LISA). We use a set of 137 fully evolved GC models that, overall, effectively match the properties of the observed GCs in the MW. We estimate that, in total, the MW GCs contain $\sim21$ sources that will be detectable by LISA. These detectable sources contain all combinations of black hole (BH), neutron star, and white dwarf components. We predict $\sim7$ of these sources will be BH-BH binaries. Furthermore, we show that some of these BH-BH binaries can have signal-to-noise ratios large enough to be detectable at the distance of the Andromeda galaxy or even the Virgo cluster.

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.

Numerical simulations have shown that black holes (BHs) can strongly influence the evolution and present-day observational properties of globular clusters (GCs). Using a Monte Carlo code, we construct GC models that match the Milky Way (MW) cluster NGC 3201, the first cluster in which a stellar-mass BH was identified through radial-velocity measurements. We predict that NGC 3201 contains $\gtrsim 200$ stellar-mass BHs. Furthermore, we explore the dynamical formation of main sequence-BH binaries and demonstrate that systems similar to the observed BH binary in NGC 3201 are produced naturally. Additionally, our models predict the existence of bright blue-straggler-BH binaries unique to core-collapsed clusters, which otherwise retain few BHs.

We explore the population of mass-transferring binaries ejected from globular clusters (GCs) with both black hole (BH) and neutron star (NS) accretors. We use a set of 137 fully evolved globular cluster models which span a large range in cluster properties and, overall, match very well the properties of old GCs observed in the Milky Way. We identify all binaries ejected from our set of models that eventually undergo mass-transfer. These binaries are ejected from their host clusters over a wide range of ejection times and include white dwarf, giant, and main sequence donors. We calculate the orbits of these ejected systems in the Galactic potential to determine their present-day positions in the Galaxy and compare to the distribution of observed low-mass X-ray binaries (XRBs) in the Milky Way. We estimate $\sim 300$ mass-transferring NS binaries and $\sim 180$ mass-transferring BH binaries may currently be present in the Milky Way that originated from within GCs. Of these, we estimate, based on mass-transfer rates and duty cycles at the present time, at most a few would be observable as BH--XRBs and NS--XRBs at the present day. Based on our results, XRBs that originated from GCs are unlikely to contribute significantly to the total population of low-mass XRBs in the Galactic field.

We explore the formation of mass-transferring binary systems containing black holes within globular clusters. We show that it is possible to form mass-transferring black hole binaries with main sequence, giant, and white dwarf companions with a variety of orbital parameters in globular clusters spanning a large range in present-day properties. All mass-transferring black hole binaries found in our models at late times are dynamically created. The black holes in these systems experienced a median of $\sim 30$ dynamical encounters within the cluster before and after acquiring the donor. Furthermore, we show that the presence of mass-transferring black hole systems has little correlation with the total number of black holes within the cluster at any time. This is because the net rate of formation of black hole-non-black hole binaries in a cluster is largely independent of the total number of retained black holes. Our results suggest that the detection of a mass-transferring black hole binary in a globular cluster does not necessarily indicate that the host cluster contains a large black hole population.

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.

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.

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.