Common envelope evolution, the key orbital tightening phase of the traditional formation channel for close binaries, is a multistage process that presents many challenges to the establishment of a fully descriptive, predictive theoretical framework. In an approach complementary to global 3D hydrodynamical modeling, we explore the range of applicability for a simplified drag formalism that incorporates the results of local hydrodynamic "wind tunnel" simulations into a semi-analytical framework in the treatment of the common envelope dynamical inspiral phase using a library of realistic giant branch stellar models across the low, intermediate, and high mass regimes. In terms of a small number of key dimensionless parameters, we characterize a wide range of common envelope events, revealing the broad range of applicability of the drag formalism as well its self-similar nature across mass regimes and ages. Limitations arising from global binary properties and local structural quantities are discussed together with the opportunity for a general prescriptive application for this formalism.

Tidal disruption events (TDEs) are among the brightest transients in the optical, ultraviolet, and X-ray sky. These flares are set into motion when a star is torn apart by the tidal field of a massive black hole, triggering a chain of events which is -- so far -- incompletely understood. However, the disruption process has been studied extensively for almost half a century, and unlike the later stages of a TDE, our understanding of the disruption itself is reasonably well converged. In this Chapter, we review both analytical and numerical models for stellar tidal disruption. Starting with relatively simple, order-of-magnitude physics, we review models of increasing sophistication, the semi-analytic ``affine formalism,'' hydrodynamic simulations of the disruption of polytropic stars, and the most recent hydrodynamic results concerning the disruption of realistic stellar models. Our review surveys the immediate aftermath of disruption in both typical and more unusual TDEs, exploring how the fate of the tidal debris changes if one considers non-main sequence stars, deeply penetrating tidal encounters, binary star systems, and sub-parabolic orbits. The stellar tidal disruption process provides the initial conditions needed to model the formation of accretion flows around quiescent massive black holes, and in some cases may also lead to directly observable emission, for example via shock breakout, gravitational waves or runaway nuclear fusion in deeply plunging TDEs.

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

White dwarf stars that enter the tidal radius of black holes with masses $\lesssim$ $10^5$ M$_{\odot}$ are doomed to be ripped apart by tidal forces. Black holes in this mass range between stellar black holes and supermassive black holes have not been conclusively identified so the detection of a tidal disruption of a white dwarf would provide clear evidence for the existence of intermediate-mass black holes. In this review, we present a theoretical and observational overview of the transient events that result from the tidal disruptions of white dwarfs by intermediate-mass black holes. This includes discussion of the latest simulations and predicted properties, the results of observational searches, as well as a summary of the potential for gravitational wave emission to be detected with upcoming missions.

Binary systems undergoing unstable Roche Lobe overflow spill gas into their circumbinary environment as their orbits decay toward coalescence. In this paper, we use a suite of hydrodynamic models of coalescing binaries involving an extended donor and a more compact accretor. We focus on the period of unstable Roche Lobe overflow that ends as the accretor plunges within the envelope of the donor at the onset of a common envelope phase. During this stage, mass is removed from the donor and flung into the circumbinary environment. Across a wide range of binary mass ratios, we find that the mass expelled as the separation decreases from the Roche limit to the donor's original radius is of the order of 25% of the accretor's mass. We study the kinematics of this ejecta and its dependencies on binary properties and find that it assembles into a toroidal circumbinary distribution. These circumbinary tori have approximately constant specific angular momentum due to momentum transport by spiral shocks launched from the orbiting binary. We show that an analytic model with these torus properties captures many of the main features of the azimuthally-averaged profiles of our hydrodynamic simulations. Our results, in particular the simple relationship between accretor mass and expelled mass and its spatial distribution, may be useful in interpreting stellar coalescence transients like luminous red novae, and in initializing hydrodynamic simulations of the subsequent common envelope phase.

Planetary engulfment events have long been proposed as a lithium (Li) enrichment mechanism contributing to the population of Li-rich giants $(A(\mathrm{Li}) \geq 1.5$ dex). Using $GALAH$ survey data and MESA stellar models, we calculate the strength and duration of the Li enrichment signature produced in the convective envelope of a host star that has engulfed a hot Jupiter (HJ) companion. We consider solar-metallicity stars in the mass range of $1.0-2.0~\mathrm{M_{\odot}}$ and the Li supplied by a HJ of $1.0~\mathrm{M_{J}}$. We explore engulfment events that occur near the main sequence turn-off (MSTO) and out to orbital separations of $R_{\star}{\sim}~0.1~\mathrm{AU}= 22~\mathrm{R_{\odot}}$. We map our results onto the Hertzsprung-Russell (H-R) Diagram, revealing the parameter space where planetary engulfment events produce significant Li enrichment signatures. We also map the associated survival times of these signatures, which range across 9 orders of magnitude. Our calculations indicate that if the HJ engulfment event occurs near the MSTO or on the subgiant branch, Li enrichment can be measured at a $5\sigma$ confidence level and with strengths that exceed meteoritic abundance measurements. Moreover, for stars of $1.4~\mathrm{M_{\odot}}$, these signatures are predicted to survive for up to 1 Gyr. We determine that Li enrichment beyond the subgiant branch must be produced by other mechanisms, such as the Cameron-Fowler process or accretion of material from an AGB companion.

Close Double Neutron Stars (DNSs) have been observed as Galactic radio pulsars, while their mergers have been detected as gamma-ray bursts and gravitational-wave sources. They are believed to have experienced at least one common-envelope episode (CEE) during their evolution prior to DNS formation. In the last decades there have been numerous efforts to understand the details of the common-envelope phase, but its computational modelling remains challenging. We present and discuss the properties of the donor and the binary at the onset of the Roche-lobe overflow (RLOF) leading to these CEEs as predicted by rapid binary population synthesis models. These properties can be used as initial conditions for detailed simulations of the common-envelope phase. There are three distinctive populations, classified by the evolutionary stage of the donor at the moment of the onset of the RLOF: giant donors with fully-convective envelopes, cool donors with partially-convective envelopes, and hot donors with radiative envelopes. We also estimate that, for standard assumptions, tides would not circularise a large fraction of these systems by the onset of RLOF. This makes the study and understanding of eccentric mass transferring systems relevant for DNS populations.

Galactic center black holes appear to be nearly universally surrounded by dense stellar clusters. When these black holes go through an active accretion phase, the multiple components of the accretion disk, stellar cluster, and black hole system all coexist. We analyze the effect of drag forces on highly eccentric stellar orbits incurred as stars puncture through the disk plane. Disk crossings dissipate orbital energy, drawing eccentric stars into more circular orbits. For high surface density disks, such as those found around black holes accreting near the Eddington mass accretion limit, the magnitude of this energy dissipation can be larger than the mean scatterings that stars receive by two body relaxation. One implication of this is the presence of a disk "loss cone" for highly eccentric stellar orbits where the dissipation from disk interaction outweighs scatter via two body relaxation. The disk loss cone is larger than the tidal disruption loss cone for near-Eddington black hole accretion rates. Stellar orbits within the disk loss cone are lost from the overall cluster as stellar orbits are circularized and stars are potentially ablated by their high-velocity impacts with the disk. We find, however, that the presence of the disk loss cone has a minimal effect on the overall rate of stellar tidal disruptions. Stars are still efficiently fed to the black hole from more-distant stellar orbits that receive large-enough per orbit scatter to jump over the disk loss cone and end up tidally disrupted.

In Kawahara et al. (2018) and Masuda et al. (2019), we reported the discovery of four self-lensing binaries consisting of F/G-type stars and (most likely) white dwarfs whose masses range from 0.2 to 0.6 solar masses. Here we present their updated system parameters based on new radial velocity data from the Tillinghast Reflector Echelle Spectrograph at the Fred Lawrence Whipple Observatory, and the Gaia parallaxes and spectroscopic parameters of the primary stars. We also briefly discuss the astrophysical implications of these findings.

We study the process of runaway, unstable Roche lobe overflow in coalescing binary systems and its dependence on the properties of the binary involved. We create three-dimensional hydrodynamic models of binary coalescences, and follow them through a phase of increasing Roche lobe overflow until the accretor is engulfed by the donor at the onset of a common envelope phase. In these models, we vary binary properties of mass ratio, donor structure and spin, and equation of state through the gas adiabatic index. We compare the numerical results to semi-analytic models of binary orbit evolution based on mass and angular momentum exchange between two point masses. Using our hydrodynamic simulations, we measure the key parameters: the donor mass loss rate and the angular momentum exchanged per unit mass loss from the donor. Using these calibrations, the semi-analytic model closely reproduces the escalating mass loss and runaway orbital decay observed in the hydrodynamic models. The semi-analytic model accurately reproduces the major differences in orbit evolution that arise with varying mass ratio and donor structure. We encapsulate the semi-analytic model in a publicly-released python package RLOF. We apply this model to the observed period decay and subsequent merger of the binary V1309 Sco, and find that it can simultaneously reproduce the observed orbital decay and time of outburst. We further demonstrate that there is a relationship between period derivative and second derivative that can be a useful metric for evaluating candidate merging binaries.

We present three-dimensional local hydrodynamic simulations of flows around objects embedded within stellar envelopes using a "wind tunnel" formalism. Our simulations model the common envelope dynamical inspiral phase in binary star systems in terms of dimensionless flow characteristics. We present suites of simulations that study the effects of varying the binary mass ratio, stellar structure, equation of state, relative Mach number of the object's motion through the gas, and density gradients across the gravitational focusing scale. For each model, we measure coefficients of accretion and drag experienced by the embedded object. These coefficients regulate the coupled evolution of the object's masses and orbital tightening during the dynamical inspiral phase of the common envelope. We extrapolate our simulation results to accreting black holes with masses comparable to that of the population of LIGO black holes. We demonstrate that the mass and spin accrued by these black holes per unit orbital tightening are directly related to the ratio of accretion to drag coefficients. We thus infer that the mass and dimensionless spin of initially non-rotating black holes change by of order $1\%$ and 0.05, respectively, in a typical example scenario. Our prediction that the masses and spins of black holes remain largely unmodified by a common envelope phase aids in the interpretation of the properties of the growing observed population of merging binary black holes. Even if these black holes passed through a common envelope phase during their assembly, features of mass and spin imparted by previous evolutionary epochs should be preserved.

Binary stars are common. While only those with small separations may exchange gas with one another, even the widest binaries interact with their gaseous surroundings. Drag forces and accretion rates dictate how these systems are transformed by these interactions. We perform three-dimensional hydrodynamic simulations of Bondi-Hoyle-Lyttleton flows, in which a binary moves supersonically relative to a homogeneous medium, using the adaptive mesh refinement code FLASH. We simulate a range of values of the initial semi-major axis of the orbit relative to the gravitational focusing impact parameter of the pair. When the binary separation is less than the gravitational focusing impact parameter, the pair orbits within a shared bow shock. When the pair is wider, each object has an individual bow-shock structure. The long-term evolution of the binary is determined by the timescales for accretion, slowing of the center of mass, and orbital inspiral. We find a clear hierarchy of these timescales; a binary's center-of-mass motion is slowed over a shorter timescale than the pair inspirals or accretes. In contrast to previous analytic predictions, which assume an unperturbed background medium, we find that the timescale for orbital inspiral is proportional to the semi-major axis to the $0.19 \pm 0.01$ power. This positive scaling indicates that gaseous drag forces can drive binaries either to coalescence or to the critical separation at which gravitational radiation dominates their further evolution. We discuss the implications of our results for binaries embedded in the interstellar medium, active galactic nuclei disks, and common envelope phases.

We report the discovery of the fifth self-lensing binary in which a low-mass white dwarf (WD) gravitationally magnifies its 15th magnitude G-star companion, KIC 8145411, during eclipses. The system was identified from a pair of such self-lensing events in the Kepler photometry, and was followed up with the Tillinghast Reflector Echelle Spectrograph (TRES) on the 1.5m telescope at the Fred Lawrence Whipple Observatory and the High-Dispersion Spectrograph (HDS) on the Subaru 8.2m telescope. A joint analysis of the TRES radial velocities, the HDS spectrum, and the Kepler photometry of the primary star determines the WD mass $0.20\pm0.01\,M_\odot$, orbital semi-major axis $1.28\pm0.03\,\mathrm{au}$, and orbital eccentricity $0.14\pm0.02$. Because such extremely low-mass WDs cannot be formed in isolation within the age of the Galaxy, their formation is believed to involve binary interactions that truncated evolution of the WD progenitor. However, the observed orbit of the KIC 8145411 system is at least ten times wider than required for this scenario to work. The presence of this system in the Kepler sample, along with its similarities to field blue straggler binaries presumably containing WDs, may suggest that some 10% of post-AGB binaries with Sun-like primaries contain such anomalous WDs.

We model explosions driven by the coalescence of a black hole or neutron star with the core of its massive-star companion. Upon entering a common envelope phase, a compact object may spiral all the way to the core. The concurrent release of energy is likely to be deposited into the surrounding common envelope, powering a merger-driven explosion. We use hydrodynamic models of binary coalescence to model the common envelope density distribution at the time of coalescence. We find toroidal profiles of material, concentrated in the binary's equatorial plane and extending to many times the massive star's original radius. We use the spherically-averaged properties of this circumstellar material (CSM) to estimate the emergent light curves that result from the interaction between the blast wave and the CSM. We find that typical merger-driven explosions are brightened by up to three magnitudes by CSM interaction. From population synthesis models we discover that the brightest merger-driven explosions, $M_V \sim -18$ to $-19$, are those involving black holes because they have the most massive and extended CSM. Black hole coalescence events are also common; they represent about 50% of all merger-driven explosions and approximately 0.3% of the core-collapse rate. Merger-driven explosions offer a window into the highly-uncertain physics of common envelope interactions in binary systems by probing the properties of systems that merge rather than eject their envelopes.

The level of mass segregation in the core of globular clusters has been previously proposed as a potential indicator of the dynamical constituents of the system, such as presence of a significant population of stellar-mass black holes (BHs), or even a central intermediate-mass black hole (IMBH). However, its measurement is limited to clusters with high-quality Hubble Space Telescope data. Thanks to a set of state-of-the-art direct N-body simulations with up to 200k particles inclusive of stellar evolution, primordial binaries, and varying BH/neutron stars, we highlight for the first time the existence of a clear and tight linear relation between the degree of mass segregation and the cluster structural concentration index. The latter is defined as the ratio of the radii containing 5% and 50% of the integrated light ($R_5/R_\mathrm{50}$), making it robustly measurable without the need to individually resolve low-mass stars. Our simulations indicate that given $R_5/R_\mathrm{50}$, the mass segregation $\Delta m$ (defined as the difference in main sequence median mass between center and half-light radius) is expressed as $\Delta m/M_\odot = -1.166 R_5/R_\mathrm{h} + 0.3246$, with a root-mean-square error of $0.0148$. In addition, we can explain its physical origin and the values of the fitted parameters through basic analytical modeling. Such correlation is remarkably robust against a variety of initial conditions (including presence of primordial binaries and IMBHs) and cluster ages, with a slight dependence in best-fit parameters on the prescriptions used to measure the quantities involved. Therefore, this study highlights the potential to develop a new observational tool to gain insight on the dynamical status of globular clusters and on its dark remnants.

The tidal disruption of stars by (super-)massive black holes in galactic nuclei has been discussed in theoretical terms for about 30 years but only in the past decade have we been able to detect such events in substantial numbers. Thus, we are now starting to carry out observational tests of models for the disruption. We are also formulating expectations for the inspiral and disruption of white dwarfs by intermediate-mass black holes with masses $< 10^5\;{\rm M}_\odot$. Such events are very rich with information and open a new window to intermediate-mass black holes, thought to live in dwarf galaxies and star clusters. They can inform us of the demographics of intermediate-mass black holes, stellar populations and dynamics in their immediate vicinity, and the physics of accretion of hydrogen-deficient material. The combination of upcoming transient surveys using ground-based, electromagnetic observatories and low-frequency gravitational wave observations is ideal for exploiting tidal disruptions of white dwarfs. The detection rate of gravitational wave signals, optimistically, may reach a few dozen per year in a volume up to $z\approx 0.1$. Gravitational wave observations are particularly useful because they yield the masses of the objects involved and allow determination of the spin of the black hole, affording tests of physical models for black hole formation and growth. They also give us advance warning of the electromagnetic flares by weeks or more. The right computing infrastructure for modern models for the disruption process and event rates will allow us to make the most of the upcoming observing facilities.

The prevalence of binary stars at close separations implies that many of these systems will interact or merge during the binary's lifetime. This paper presents hydrodynamic simulations of the scenario of binary coalescence through unstable mass transfer, which drives the pair to closer separations. When the donor star does not rotate synchronously with respect to the orbit, dynamical tidal waves are excited in its envelope. We show that resonance crossings with high azimuthal-order $(m\sim3$ to $6$) fundamental modes induce a visible "polygram" distortion to the star. As the binary orbit tightens, the system sweeps through resonance with modes of decreasing azimuthal order, which are selectively excited. We compare our hydrodynamic simulations to predictions from linear theory of resonant mode excitation. The linear theory provides an estimate of mode amplitudes to within a factor of two, even as the oscillations become quite non-linear as the stars coalesce. We estimate that a wave with 10% radial amplitude generates approximately 1% photometric variability; this may be detectible if such a binary coalescence is caught in action by future photometric all-sky surveys.

Recent observations have revealed that the remnants of stellar-coalescence transients are bipolar. This raises the questions of how these bipolar morphologies arise and what they teach us about the mechanisms of mass ejection during stellar mergers and common-envelope phases. In this paper, we analyze hydrodynamic simulations of the lead-in to binary coalescence, a phase of unstable Roche lobe overflow that takes the binary from the Roche limit separation to the engulfment of the more compact accretor within the envelope of the extended donor. As mass transfer runs away at increasing rates, gas trails away from the binary. Contrary to previous expectations, early mass loss from the system remains bound to the binary and forms a circumbinary torus. Later ejecta, generated as the accretor grazes the surface of the donor, have very different morphologies and are unbound. These two components of mass loss from the binary interact as later, higher-velocity ejecta collide with the circumbinary torus formed by earlier mass loss. Unbound ejecta are redirected toward the poles, and escaping material creates a bipolar outflow. Our findings show that the transition from bound to unbound ejecta from coalescing binaries can explain the bipolar nature of their remnants, with implications for our understanding of the origin of bipolar remnants of stellar-coalescence transients and, perhaps, some preplanetary nebulae.

Luminous red nova transients, presumably from stellar coalescence, exhibit long-term precursor emission over hundreds of binary orbits, leading to impulsive outbursts with durations similar to a single orbital period. In an effort to understand these signatures, we present and analyze a hydrodynamic model of unstable mass transfer from a giant-star donor onto a more compact accretor in a binary system. Our simulation begins with mass transfer at the Roche limit separation and traces a phase of runaway decay leading up to the plunge of the accretor within the envelope of the donor. We characterize the fluxes of mass and angular momentum through the system and show that the orbital evolution can be reconstructed from measurements of these quantities. The morphology of outflow from the binary changes significantly as the binary orbit tightens. At wide separations, a thin stream of relatively high-entropy gas trails from the outer Lagrange points. As the orbit tightens, the orbital motion desynchronizes from the donor's rotation, and low-entropy ejecta trace a broad fan of largely ballistic trajectories. An order-of-magnitude increase in mass ejection rate accompanies the plunge of the accretor with the envelope of the donor. We argue that this transition marks the precursor-to-outburst transition observed in stellar coalescence transients.

Despite recent observational efforts, unequivocal signs for the presence of intermediate-mass black holes (IMBHs) in globular clusters (GCs) have not been found yet. Especially when the presence of IMBHs is constrained through dynamical modeling of stellar kinematics, it is fundamental to account for the displacement that the IMBH might have with respect to the GC centre. In this paper we analyse the IMBH wandering around the stellar density centre using a set of realistic direct N-body simulations of star cluster evolution. Guided by the simulation results, we develop a basic yet accurate model that can be used to estimate the average IMBH radial displacement ($\left<r_\mathrm{bh}\right>$) in terms of structural quantities as the core radius ($r_\mathrm{c}$), mass ($M_\mathrm{c}$), and velocity dispersion ($\sigma_\mathrm{c}$), in addition to the average stellar mass ($m_\mathrm{c}$) and the IMBH mass ($M_\mathrm{bh}$). The model can be expressed by the equation $\left<r_\mathrm{bh}\right>/r_\mathrm{c}=A(m_\mathrm{c}/M_\mathrm{bh})^\alpha[\sigma_\mathrm{c}^2r_\mathrm{c}/(GM_\mathrm{c})]^\beta$, in which the free parameters $A,\alpha,\beta$ are calculated through comparison with the numerical results on the IMBH displacement. The model is then applied to Galactic GCs, finding that for an IMBH mass equal to 0.1% of the GC mass, the typical expected displacement of a putative IMBH is around $1''$ for most Galactic GCs, but IMBHs can wander to larger angular distances in some objects, including a prediction of a $2.5''$ displacement for NGC 5139 ($\omega$ Cen), and $>10''$ for NGC 5053, NGC 6366 and ARP2.

We report the discovery of three edge-on binaries with white dwarf companions that gravitationally magnify (instead of eclipsing) the light of their stellar primaries, as revealed by a systematic search for pulses with long periods in the Kepler photometry. We jointly model the self-lensing light curves and radial-velocity orbits to derive the white dwarf masses, all of which are close to 0.6 Solar masses. The orbital periods are long, ranging from 419 to 728 days, and the eccentricities are low, all less than 0.2. These characteristics are reminiscent of the orbits found for many blue stragglers in open clusters and the field, for which stable mass transfer due to Roche-lobe overflow from an evolving primary (now a white dwarf) has been proposed as the formation mechanism. Because the actual masses for our three white dwarf companions have been accurately determined, these self-lensing systems would provide excellent tests for models of interacting binaries.

Planets accompany most sun-like stars. The orbits of many are sufficiently close that they will be engulfed when their host stars ascend the giant branch. This Letter compares the power generated by orbital decay of an engulfed planet to the intrinsic stellar luminosity. Orbital decay power is generated by drag on the engulfed companion by the surrounding envelope. As stars ascend the giant branch their envelope density drops and so does the power injected through orbital decay, scaling approximately as $L_{\rm decay} \propto R_*^{-9/2}$. Their luminosity, however, increases along the giant branch. These opposed scalings indicate a crossing, where $L_{\rm decay}= L_*$. We consider the engulfment of planets along isochrones in the Hertzsprung-Russell (H-R) diagram. We find that the conditions for such a crossing occur around $L_*\approx 10^2$~$L_\odot$ (or $a\approx 0.1$~au) for Jovian planetary companions. The consumption of closer-in giant planets, such as hot Jupiters, leads to $L_{\rm decay}\gg L_*$, while more distant planets such as warm Jupiters, $a\approx 0.5$~au, lead to minor perturbations of their host stars with $L_{\rm decay} \ll L_*$. Our results map out the parameter space along the giant branch in the H-R Diagram where interaction with planetary companions leads to significant energetic disturbance of host stars.

In this paper, we model the observable signatures of tidal disruptions of white dwarf (WD) stars by massive black holes (MBHs) of moderate mass, $\approx 10^3 - 10^5 M_\odot$. When the WD passes deep enough within the MBH's tidal field, these signatures include thermonuclear transients from burning during maximum compression. We combine a hydrodynamic simulation that includes nuclear burning of the disruption of a $0.6 M_\odot$ C/O WD with a Monte Carlo radiative transfer calculation to synthesize the properties of a representative transient. The transient's emission emerges in the optical, with lightcurves and spectra reminiscent of type I SNe. The properties are strongly viewing-angle dependent, and key spectral signatures are $\approx 10,000$ km s$^{-1}$ Doppler shifts due to the orbital motion of the unbound ejecta. Disruptions of He WDs likely produce large quantities of intermediate-mass elements, offering a possible production mechanism for Ca-rich transients. Accompanying multiwavelength transients are fueled by accretion and arise from the nascent accretion disk and relativistic jet. If MBHs of moderate mass exist with number densities similar to those of supermassive MBHs, both high energy wide-field monitors and upcoming optical surveys should detect tens to hundreds of WD tidal disruptions per year. The current best strategy for their detection may therefore be deep optical follow up of high-energy transients of unusually-long duration. The detection rate or the non-detection of these transients by current and upcoming surveys can thus be used to place meaningful constraints on the extrapolation of the MBH mass function to moderate masses.

We present a study on binary-single interactions with energy loss terms such as tidal dissipation and gravitational wave emission added to the equation-of-motion. The inclusion of such terms leads to the formation of compact binaries that form during the three-body interaction through two-body captures. These binaries predominantly merge relative promptly at high eccentricity, with several observable and dynamical consequences to follow. Despite their possibility for being observed in both present and upcoming transient surveys, their outcomes are not firmly constrained. In this paper we present an analytical framework that allows to estimate the cross section of such two-body captures, which permits us to study how the corresponding rates depends on the initial orbital parameters, the mass hierarchy, the type of interacting objects, and the energy dissipation mechanism. This formalism is applied here to study the formation of two-body gravitational wave captures, for which we estimate absolute and relative rates relevant to Advanced LIGO detections. It is shown that two-body gravitational wave captures should have compelling observational implications if a sizable fraction of detected compact binaries are formed via dynamical interactions.

We construct a menu of objects that can give rise to bright flares when disrupted by massive black holes (BHs), ranging from planets to evolved stars. Through their tidal disruption, main sequence and evolved stars can effectively probe the existence of otherwise quiescent supermassive BHs and white dwarfs can probe intermediate mass BHs. Many low-mass white dwarfs possess extended hydrogen envelopes, which allow for the production of prompt flares in disruptive encounters with moderately massive BHs of $10^5$ to $10^7~M_\odot$--masses that may constitute the majority of massive BHs by number. These objects are a missing link in two ways: (1) for probing moderately massive BHs and (2) for understanding the hydrodynamics of the disruption of objects with tenuous envelopes. A flare arising from the tidal disruption of a $0.17~M_\odot$ white dwarf by a $10^5~M_\odot$ BH reaches a maximum between 0.6 and 11 days, with a peak fallback rate that is usually super-Eddington and results in a flare that is likely brighter than a typical tidal disruption event. Encounters stripping only the envelope can provide hydrogen-only fallback, while encounters disrupting the core evolve from H- to He-rich fallback. While most tidal disruption candidates observed thus far are consistent with the disruptions of main sequence stars, the rapid timescales of nuclear transients such as Dougie and PTF10iya are naturally explained by the disruption of low-mass white dwarfs. As the number of observed flares continues to increase, the menu presented here will be essential for characterizing nuclear BHs and their environments through tidal disruptions.

During a common envelope episode in a binary system, the engulfed companion spirals to tighter orbital separations under the influence of drag from the surrounding envelope material. As this object sweeps through material with a steep radial gradient of density, net angular momentum is introduced into the flow, potentially leading to the formation of an accretion disk. The presence of a disk would have dramatic consequences for the outcome of the interaction because accretion might be accompanied by strong, polar outflows with enough energy to unbind the entire envelope. Without a detailed understanding of the necessary conditions for disk formation during common envelope, therefore, it is difficult to accurately predict the population of merging compact binaries. This paper examines the conditions for disk formation around objects embedded within common envelopes using the `wind tunnel' formalism developed by MacLeod et al. (2017). We find that the formation of disks is highly dependent on the compressibility of the envelope material. Disks form only in the most compressible of stellar envelope gas, found in envelopes' outer layers in zones of partial ionization. These zones are largest in low-mass stellar envelopes, but comprise small portions of the envelope mass and radius in all cases. We conclude that disk formation and associated accretion feedback in common envelope is rare, and if it occurs, transitory. The implication for LIGO black hole binary assembly is that by avoiding strong accretion feedback, common envelope interactions should still result in the substantial orbital tightening needed to produce merging binaries.

This paper examines the properties of flows around objects embedded within common envelopes in the simplified context of a "wind tunnel." We establish characteristic relationships between key common envelope flow parameters like the Mach number and density scale height. Our wind tunnel is a three-dimensional, cartesian geometry hydrodynamic simulation setup that includes the gravity of the primary and secondary stars and allows us to study the coefficients of drag and accretion experienced by the embedded object. Accretion and drag lead to a transformation of an embedded object and its orbit during a common envelope phase. We present two suites of simulations spanning a range of density gradients and Mach numbers -- relevant for flow near the limb of a stellar envelope to the deep interior. In one suite, we adopt an ideal gas adiabatic exponent of $\gamma=5/3$, in the other, $\gamma=4/3$. We find that coefficients of drag rise in flows with steeper density gradients and that coefficients of drag and accretion are consistently higher in the more compressible, $\gamma=4/3$ flow. We illustrate the impact of these newly derived coefficients by integrating the inspiral of a secondary object through the envelopes of $3M_\odot$ ($\gamma\approx5/3$) and $80M_\odot$ ($\gamma\approx4/3$) giants. In these examples, we find a relatively rapid initial inspiral because, near the stellar limb, dynamical friction drag is generated mainly from dense gas focussed from deeper within the primary-star's envelope. This rapid initial inspiral timescale carries potential implications for the timescale of transients from early common envelope interaction.

This paper investigates the recent stellar merger transient M31LRN 2015 in the Andromeda galaxy. We analyze published optical photometry and spectroscopy along with a Hubble Space Telescope detection of the color and magnitude of the pre-outburst source. The transient outburst is consistent with dynamically driven ejecta at the onset of a common envelope episode, which eventually leads to the complete merger of a binary system. The light curve appears to contain two components: first $\sim10^{-2} M_\odot$ of fast ejecta driven by shocks at the onset of common envelope, and later, $\sim0.3 M_\odot$ of further ejecta as the secondary becomes more deeply engulfed within the primary. Just prior to merger, we find that the primary star is a $3-5.5 M_\odot$ sub-giant branch primary star with radius of $30-40R_\odot$. Its position in the color-magnitude diagram shows that it is growing in radius, consistent with a picture where it engulfs its companion. By matching the properties of the primary star to the transient outburst, we show that the optical transient lasts less than ten orbits of the original binary, which had a pre-merger period of $\sim 10$ days. We consider the possible orbital dynamics leading up to the merger, and argue that if the system merged due to the Darwin tidal instability it implies a lower mass main sequence companion of $0.1-0.6M_\odot$. This analysis represents a promising step toward a more detailed understanding of flows in common envelope episodes through direct observational constraints.

We perform the first systematic study on how dynamical stellar tides and general relativistic (GR) effects affect the dynamics and outcomes of binary-single interactions. For this, we have constructed an N-body code that includes tides in the affine approximation, where stars are modeled as self-similar ellipsoidal polytropes, and GR corrections using the commonly-used post-Newtonian formalism. Using this numerical formalism, we are able resolve the leading effect from tides and GR across several orders of magnitude in both stellar radius and initial target binary separation. We find that the main effect from tides is the formation of two-body tidal captures that form during the chaotic and resonant evolution of the triple system. The two stars undergoing the capture spiral in and merge. The inclusion of tides can thus lead to an increase on the stellar coalescence rate. We also develop an analytical framework for calculating the cross section of tidal inspirals between any pair of objects with similar mass. From our analytical and numerical estimates we find that the rate of tidal inspirals relative to collisions increases as the initial semi-major axis of the target binary increases and the radius of the interacting tidal objects decreases. The largest effect is therefore found for triple systems hosting white dwarfs and neutron stars. In this case, we find the rate of highly eccentric white dwarf - neutron star mergers to likely be dominated by tidal inspirals. While tidal inspirals occur rarely, we note that they can give rise to a plethora of thermonuclear transients such as Ca-rich transients.

When embedded in dense cluster cores, intermediate mass black holes (IMBHs) acquire close stellar or stellar-remnant companions. These companions are not only gravitationally bound, they tend to hierarchically isolate from other cluster stars through series of multibody encounters. In this paper, we study the demographics of IMBH companions in compact star clusters through direct $N$-body simulation. We study clusters initially composed of $10^5$ or $2\times 10^5$ stars with IMBHs of 75 and 150 solar masses, and follow their evolution for 6-10 Gyr. A tight innermost binary pair of IMBH and stellar object rapidly forms. The IMBH has a companion with orbital semi-major axis at least three times tighter than the second-most bound object over 90% of the time. These companionships have typical periods of order years and are subject to cycles of exchange and destruction. The most frequently observed, long-lived pairings persist for $\sim 10^7$ yr. The demographics of IMBH companions in clusters are diverse; they include both main sequence, giant stars, and stellar remnants. Companion objects may reveal the presence of an IMBH in a cluster in one of several ways. Most-bound companion stars routinely suffer grazing tidal interactions with the IMBH, offering a dynamical mechanism to produce repeated flaring episodes like those seen in the IMBH candidate HLX-1. Stellar winds of companion stars provide a minimum quiescent accretion rate for IMBHs, with implications for radio searches for IMBH accretion in globular clusters. Finally, gravitational wave inspirals of compact objects are found to occur with promising frequency.

Investigations of elemental abundances in the ancient and most metal deficient stars are extremely important because they serve as tests of variable nucleosynthesis pathways and can provide critical inferences of the type of stars that lived and died before them. The presence of r-process elements in a handful of carbon-enhanced metal-poor (CEMP-r) stars, which are assumed to be closely connected to the chemical yield from the first stars, is hard to reconcile with standard neutron star mergers. Here we show that the production rate of dynamically assembled compact binaries in high-z nuclear star clusters can attain a sufficient high value to be a potential viable source of heavy r-process material in CEMP-r stars. The predicted frequency of such events in the early Galaxy, much lower than the frequency of Type II supernovae but with significantly higher mass ejected per event, can naturally lead to a high level of scatter of Eu as observed in CEMP-r stars.

This paper examines flows in the immediate vicinity of stars and compact objects dynamically inspiralling within a common envelope (CE). Flow in the vicinity of the embedded object is gravitationally focused leading to drag and potential to gas accretion. This process has been studied numerically and analytically in the context of Hoyle-Lyttleton accretion (HLA). Yet, within a CE, accretion structures may span a large fraction of the envelope radius, and in so doing sweep across a substantial radial gradient of density. We quantify these gradients using detailed stellar evolution models for a range of CE encounters. We provide estimates of typical scales in CE encounters that involve main sequence stars, white dwarfs, neutron stars, and black holes with giant-branch companions of a wide range of masses. We apply these typical scales to hydrodynamic simulations of 3D HLA with an upstream density gradient. This density gradient breaks the symmetry that defines HLA flow, and imposes an angular momentum barrier to accretion. Material that is focused into the vicinity of the embedded object thus may not be able to accrete. As a result, accretion rates drop dramatically, by 1-2 orders of magnitude, while drag rates are only mildly affected. We provide fitting formulae to the numerically-derived rates of drag and accretion as a function of the density gradient. The reduced ratio of accretion to drag suggests that objects that can efficiently gain mass during CE evolution, such as black holes and neutron stars, may grow less than implied by the HLA formalism.

This paper models the orbital inspiral of a neutron star (NS) through the envelope of its giant-branch companion during a common envelope (CE) episode. These CE episodes are necessary to produce close pairs of NSs that can inspiral and merge due to gravitational wave losses in less than a Hubble time. Because cooling by neutrinos can be very efficient, NSs have been predicted to accumulate significant mass during CE events, perhaps enough to lead them to collapse to black holes. We revisit this conclusion with the additional consideration of CE structure, particularly density gradients across the embedded NS's accretion radius. This work is informed by our recent numerical simulations that find that the presence of a density gradient strongly limits accretion by imposing a net angular momentum to the flow around the NS. Our calculations suggest that NSs should survive CE encounters. They accrete only modest amounts of envelope material, $\lesssim 0.1M_\odot$, which is broadly consistent with mass determinations of double NS binaries. With less mass gain, NSs must spiral deeper to eject their CE, leading to a potential increase in mergers. The survival of NSs in CE events has implications for the formation mechanism of observed double NS binaries, as well as for predicted rates of NS binary gravitational wave inspirals and their electromagnetic counterparts.

White dwarfs (WDs) can be tidally disrupted only by massive black holes (MBHs) with masses less than $\sim10^5 M_\odot$. These tidal interactions feed material to the MBH well above its Eddington limit, with the potential to launch a relativistic jet. The corresponding beamed emission is a promising signpost to an otherwise quiescent MBH of relatively low mass. We show that the mass transfer history, and thus the lightcurve, are quite different when the disruptive orbit is parabolic, eccentric, or circular. The mass lost each orbit exponentiates in the eccentric-orbit case leading to the destruction of the WD after several tens of orbits. We examine the stellar dynamics of clusters surrounding MBHs to show that single-passage WD disruptions are substantially more common than repeating encounters. The $10^{49}$ erg s$^{-1}$ peak luminosity of these events makes them visible to cosmological distances. They may be detectible at rates of as many as tens per year by instruments like Swift. In fact, WD-disruption transients significantly outshine their main-sequence star counterparts, and are the most likely tidal interaction to be detected arising from MBHs with masses less than $10^5 M_\odot$. The detection or non-detection of such WD-disruption transients by Swift is, therefore, a powerful tool to constrain lower end of the MBH mass function. The emerging class of ultra-long gamma ray bursts all have peak luminosities and durations reminiscent of WD disruptions, offering a hint that WD-disruption transients may already be present in existing datasets.

The discovery of the gas cloud G2 on a near-radial orbit about Sgr A* has prompted much speculation on its origin. In this Letter, we propose that G2 formed out of the debris stream produced by the removal of mass from the outer envelope of a nearby giant star. We perform hydrodynamical simulations of the returning tidal debris stream with cooling, and find that the stream condenses into clumps that fall periodically onto Sgr A*. We propose that one of these clumps is the observed G2 cloud, with the rest of the stream being detectable at lower Br-$\gamma$ emissivity along a trajectory that would trace from G2 to the star that was partially disrupted. By simultaneously fitting the orbits of S2, G2, and $\sim$ 2,000 candidate stars, and by fixing the orbital plane of each candidate star to G2 (as is expected for a tidal disruption), we find that several stars have orbits that are compatible with the notion that one of them was tidally disrupted to produce G2. If one of these stars were indeed disrupted, it last encountered Sgr A* hundreds of years ago, and has likely encountered Sgr A* repeatedly. However, while these stars are compatible with the giant disruption scenario given their measured positions and proper motions, their radial velocities are currently unknown. If one of these stars' radial velocity is measured to be compatible with a disruptive orbit, it would strongly suggest its disruption produced G2.

Stars may be tidally disrupted if, in a single orbit, they are scattered too close to a supermassive black hole (SMBH). Tidal disruption events are thought to power luminous but short-lived accretion episodes that can light up otherwise quiescent SMBHs in transient flares. Here we explore a more gradual process of tidal stripping where stars approach the tidal disruption radius by stellar evolution while in an eccentric orbit. After the onset of mass transfer, these stars episodically transfer mass to the SMBH every pericenter passage giving rise to low-level flares that repeat on the orbital timescale. Giant stars, in particular, will exhibit a runaway response to mass loss and "spoon-feed" material to the black hole for tens to hundreds of orbital periods. In contrast to full tidal disruption events, the duty cycle of this feeding mode is of order unity for black holes with mass greater than approximately 10 million solar masses. This mode of quasi-steady SMBH feeding is competitive with indirect SMBH feeding through stellar winds, and spoon-fed giant stars may play a significant role in determining the quiescent luminosity of local SMBHs.

The inspiral and merger of eccentric binaries leads to gravitational waveforms distinct from those generated by circularly merging binaries. Dynamical environments can assemble binaries with high eccentricity and peak frequencies within the {\it LIGO} band. In this paper, we study binary-single stellar scatterings occurring in dense stellar systems as a source of eccentrically-inspiraling binaries. Many interactions between compact binaries and single objects are characterized by chaotic resonances in which the binary-single system undergoes many exchanges before reaching a final state. During these chaotic resonances, a pair of objects has a non-negligible probability of experiencing a very close passage. Significant orbital energy and angular momentum are carried away from the system by gravitational wave (GW) radiation in these close passages and in some cases this implies an inspiral time shorter than the orbital period of the bound third body. We derive the cross section for such dynamical inspiral outcomes through analytical arguments and through numerical scattering experiments including GW losses. We show that the cross section for dynamical inspirals grows with increasing target binary semi-major axis, $a$, and that for equal-mass binaries it scales as $a^{2/7}$. Thus, we expect wide target binaries to predominantly contribute to the production of these relativistic outcomes. We estimate that eccentric inspirals account for approximately one percent of dynamically assembled non-eccentric merging binaries. While these events are rare, we show that binary-single scatterings are a more effective formation channel than single-single captures for the production of eccentrically-inspiraling binaries, even given modest binary fractions.

Sun-like stars are thought to be regularly disrupted by supermassive black holes (SMBHs) within galactic nuclei. Yet, as stars evolve off the main sequence their vulnerability to tidal disruption increases drastically as they develop a bifurcated structure consisting of a dense core and a tenuous envelope. Here we present the first hydrodynamic simulations of the tidal disruption of giant stars and show that the core has a substantial influence on the star's ability to survive the encounter. Stars with more massive cores retain large fractions of their envelope mass, even in deep encounters. Accretion flares resulting from the disruption of giant stars should last for tens to hundreds of years. Their characteristic signature in transient searches would not be the $t^{-5/3}$ decay typically associated with tidal disruption events, but a correlated rise over many orders of magnitude in brightness on months to years timescales. We calculate the relative disruption rates of stars of varying evolutionary stages in typical galactic centers, then use our results to produce Monte Carlo realizations of the expected flaring event populations. We find that the demographics of tidal disruption flares are strongly dependent on both stellar and black hole mass, especially near the limiting SMBH mass scale of $\sim 10^8 M_\odot$. At this black hole mass, we predict a sharp transition in the SMBH flaring diet beyond which all observable disruptions arise from evolved stars, accompanied by a dramatic cutoff in the overall tidal disruption flaring rate. Black holes less massive than this limiting mass scale will show observable flares from both main sequence and evolved stars, with giants contributing up to 10% of the event rate. The relative fractions of stars disrupted at different evolutionary states can constrain the properties and distributions of stars in galactic nuclei other than our own.

White dwarf-neutron star binaries generate detectable gravitational radiation. We construct Newtonian equilibrium models of corotational white dwarf-neutron star (WDNS) binaries in circular orbit and find that these models terminate at the Roche limit. At this point the binary will undergo either stable mass transfer (SMT) and evolve on a secular time scale, or unstable mass transfer (UMT), which results in the tidal disruption of the WD. The path a given binary will follow depends primarily on its mass ratio. We analyze the fate of known WDNS binaries and use population synthesis results to estimate the number of LISA-resolved galactic binaries that will undergo either SMT or UMT. We model the quasistationary SMT epoch by solving a set of simple ordinary differential equations and compute the corresponding gravitational waveforms. Finally, we discuss in general terms the possible fate of binaries that undergo UMT and construct approximate Newtonian equilibrium configurations of merged WDNS remnants. We use these configurations to assess plausible outcomes of our future, fully relativistic simulations of these systems. If sufficient WD debris lands on the NS, the remnant may collapse, whereby the gravitational waves from the inspiral, merger, and collapse phases will sweep from LISA through LIGO frequency bands. If the debris forms a disk about the NS, it may fragment and form planets.