Search

The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow-up kilonova observations have shown that NSMs produce copious amounts of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of double-neutron stars (DNSs) with the levels of r-process enrichment seen in ultra-faint dwarf (UFD) galaxies such as Reticulum II and Tucana III. Assuming that such dwarf systems have lost a significant fraction of their stellar mass through tidal stripping, we conclude that contrary to most current models, it is the DNSs with rather large natal kicks but very short merging timescales that can enrich UFD-type galaxies. These binaries are either on highly eccentric orbits or form with very short separations due to an additional mass transfer between the first-born neutron star and a naked helium star, the progenitor of the second neutron star. These DNSs are born with a frequency that agrees with the statistics of the r-process UFDs, and merge well within the virial radius of their host halos, therefore contributing significantly to their r-process enrichment. [ABSTRACT FROM AUTHOR]

We analyze data from several studies of metal-poor stars in the Milky Way, focusing individually on the main r-process elements (Eu) as well as the lighter neutron-capture element Sr, at the neutron-magic peak N = 50. Because these elements were injected in an explosion, we calculate the mass swept up when the blast wave first becomes radiative, yielding a lower limit for the dilution of such elements and hence a lower limit on the ejecta mass that is incorporated into the next generation of stars. Our study demonstrates that in order to explain the largest enhancements in [Eu/Fe] observed in stars at low [Fe/H] metallicities, individual r-process production events must synthesize a minimum of roughly 10−3M⊙ of r-process material. This provides a critical constraint on galactic chemical evolution models. We also show independently that if the site of Mg production is the same as that of Eu, individual injection events must synthesize up to ∼10−3M⊙ of r-process material. On the other hand, demanding that Sr traces Mg production results in r-process masses per event of ∼10−5M⊙ . This suggests that the astrophysical sites responsible for the genesis of the main r-process elements need to operate at a drastically reduced rate when compared to standard core-collapse supernovae. [ABSTRACT FROM AUTHOR]

One of the main challenges of current tidal disruption events (TDEs) studies is that emission arising from active galactic nucleus (AGN) activity may potentially mimic the expected X-ray emission of a TDE. Here we compare the X-ray properties of TDEs and AGNs to determine a set of characteristics that would allow us to discriminate between flares arising from these two objects. We find that at peak, TDEs are brighter than AGNs found at similar redshifts. However, compared to pre-flare upper limits, highly variable AGNs can produce flares of a similar order of magnitude as those seen from X-ray TDEs. Nevertheless, TDEs decay significantly more monotonically, and their emission exhibits little variation in spectral hardness as a function of time. We also find that X-ray TDEs are less absorbed, and their emission is much softer than the emission detected from AGNs found at similar redshifts. We derive the X-ray luminosity function (LF) for X-ray TDEs using the events from Auchettl et al. Interestingly, our X-ray LF closely matches the theoretically derived LF by Milosavljević et al., which assumes a higher TDE rate currently estimated from observations. Using our results and the results of Stone & Metzger, we estimate a TDE rate of (0.7–4.7) × 10−4 year−1 per galaxy, higher than current observational estimates. We find that TDEs can contribute significantly to the LF of AGNs for z ≲ 0.4, while there is no evidence that TDEs influence the growth of 106–7M⊙ BHs. However, BHs < 106M⊙ can grow from TDEs arising from super-Eddington accretion without contributing significantly to the observed AGN LF at z = 0. [ABSTRACT FROM AUTHOR]

We study the properties of tidal disruption event (TDE) host galaxies in the context of a catalog of ∼500,000 galaxies from the Sloan Digital Sky Survey. We explore whether selection effects can account for the overrepresentation of TDEs in E+A/post-starburst galaxies by creating matched galaxy samples. Accounting for possible selection effects due to black hole (BH) mass, redshift completeness, strong active galactic nucleus presence, bulge colors, and surface brightness can reduce the apparent overrepresentation of TDEs in E+A host galaxies by a factor of ∼4 (from ∼×100–190 to ∼×25–48), but cannot fully explain the preference. We find that TDE host galaxies have atypical photometric properties compared to similar, “typical” galaxies. In particular, TDE host galaxies tend to live in or near the “green valley” between star-forming and passive galaxies, and have bluer bulge colors ( mag), lower half-light surface brightnesses (by ∼1 mag/arcsec2), higher Sérsic indices (), and higher bulge-to-total-light ratios () than galaxies with matched BH masses. We find that TDE host galaxies appear more centrally concentrated and that all have high galaxy Sérsic indices and B/T fractions—on average in the top 10% of galaxies of the same BH mass—suggesting a higher nuclear stellar density. We identify a region in the Sérsic index and BH mass parameter space that contains ∼2% of our reference catalog galaxies but of TDE host galaxies. The unique photometric properties of TDE host galaxies may be useful for selecting candidate TDEs for spectroscopic follow-up observations in large transient surveys. [ABSTRACT FROM AUTHOR]

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 of fast ejecta driven by shocks at the onset of common envelope, and later, ∼0.3 M⊙ 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⊙ subgiant-branch primary star with a radius of 30–40 R⊙. 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 10 orbits of the original binary, which had a pre-merger period of ∼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.6 M⊙. This analysis represents a promising step toward a more detailed understanding of flows in common envelope episodes through direct observational constraints. [ABSTRACT FROM AUTHOR]

The r-process nuclei are robustly synthesized in the material ejected during neutron star binary mergers (NSBMs). If NSBMs are indeed solely responsible for the solar system r-process abundances, a galaxy like our own would be required to host a few NSBMs per million years, with each event ejecting, on average, about 5 × 10−2M⊙ of r-process material. Because the ejecta velocities in the tidal tail are significantly larger than those in ordinary supernovae, NSBMs deposit a comparable amount of energy into the ISM. In contrast to extensive efforts studying spherical models for supernova remnant evolution, calculations quantifying the impact of NSBM ejecta in the ISM have been lacking. To better understand their evolution, we perform a suite of three-dimensional hydrodynamic simulations of isolated NSBM ejecta expanding in environments with conditions adopted from Milky-Way-like galaxy simulations. Although the remnant morphology is highly complex at early times, the subsequent radiative evolution is remarkably similar to that of a standard supernova. This implies that sub-resolution supernova feedback models can be used in galaxy-scale simulations that are unable to resolve the key evolutionary phases of NSBMs. Among other quantities, we examine the radius, mass, and kinetic energy content of the remnant at shell formation. We find that the shell formation epoch is attained when the swept-up mass is about 103(nH/1 cm−3)−2/7M⊙; at this point, the mass fraction of r-process material is enhanced up to two orders of magnitude in relation to a solar metallicity ISM. [ABSTRACT FROM AUTHOR]

Main-sequence disruptions of stars by supermassive black holes result in the production of an extended, geometrically thin debris stream winding repeatedly around the black hole. In the absence of black hole spin, in-plane relativistic precession causes this stream to intersect with itself after a single winding. In this paper we show that relativistic precessions arising from black hole spin can induce deflections out of the original orbital plane that prevent the stream from self-intersecting even after many windings. This naturally leads to a “dark period” in which the flare is not observable for some time, persisting for up to a dozen orbital periods of the most bound material, which translates to years for disruptions around black holes with masses . When the stream eventually self-intersects, the distance from the black hole and the angle at which this collision occurs determine the rate of energy dissipation. We find that more-massive black holes () tend to have more violent stream self-intersections, resulting in prompt accretion. For these tidal disruption events (TDEs), the accretion rate onto the black hole should still closely follow the original fallback rate after a fixed delay time , . For lower black hole masses (), we find that flares are typically slowed down by about an order of magnitude, resulting in the majority of TDEs being sub-Eddington at peak. This also implies that current searches for TDEs are biased toward prompt flares, with slowed flares likely having been unidentified. [ABSTRACT FROM AUTHOR]

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 potentially 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 three-dimensional 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 one to two 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. [ABSTRACT FROM AUTHOR]

A neo-neutron star is a hot neutron star that has just become transparent to neutrinos. In a core-collapse supernova or accretion-induced collapse of a white dwarf, the neo-neutron star phase directly follows the proto-neutron star phase, about 30–60 s after the initial collapse. It will also be present in a binary neutron star merger in the case where the “born-again” hot massive compact star does not immediately collapse into a black hole. Eddington or even super-Eddington luminosities are present for some time. A neo-neutron star produced in a core-collapse supernova is not directly observable, but the one produced by a binary merger, likely associated with an off-axis short gamma-ray burst, may be observable for some time as well as when produced in the accretion-induced collapse of a white dwarf. We present a first step in the study of this neo-neutron star phase in a spherically symmetric configuration, thus ignoring fast rotation and also ignoring the effect of strong magnetic fields. We put particular emphasis on determining how long the star can sustain a near-Eddington luminosity and also show the importance of positrons and contraction energy during the neo-neutron star phase. We finally discuss the observational prospects for neutron star mergers triggered by LIGO and for accretion-induced collapse transients. [ABSTRACT FROM AUTHOR]

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 semimajor 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 inspiraling. 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 inspiraling is proportional to the semimajor axis to the 0.19 ± 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. [ABSTRACT FROM AUTHOR]

An isolated star moving supersonically through a uniform gas accretes material from its gravitationally induced wake. The rate of accretion is set by the accretion radius of the star and is well described by classical Bondi–Hoyle–Lyttleton theory. Stars, however, are not born in isolation. They form in clusters where they accrete material that is influenced by all the stars in the cluster. We perform three-dimensional hydrodynamic simulations of clusters of individual accretors embedded in a uniform-density wind in order to study how the accretion rates experienced by individual cluster members are altered by the properties of the ambient gas and the cluster itself. We study accretion as a function of number of cluster members, mean separation between them, and size of their individual accretion radii. We determine the effect of these key parameters on the aggregate and individual accretion rates, which we compare to analytic predictions. We show that when the accretion radii of the individual objects in the cluster substantially overlap, the surrounding gas is effectively accreted into the collective potential of the cluster prior to being accreted onto the individual stars. We find that individual cluster members can accrete drastically more than they would in isolation, in particular when the flow is able to cool efficiently. This effect could potentially modify the luminosity of accreting compact objects in star clusters and could lead to the rejuvenation of young star clusters as well as globular clusters with low inclination and low eccentricity. [ABSTRACT FROM AUTHOR]

While once rare, observations of stars being tidally disrupted by supermassive black holes are quickly becoming commonplace. To continue to learn from these events, it is necessary to robustly and systematically compare our growing number of observations with theory. We present a tidal disruption module for the Modular Open Source Fitter for Transients (MOSFiT) and the results from fitting 14 tidal disruption events (TDEs). Our model uses FLASH simulations of TDEs to generate bolometric luminosities and passes these luminosities through viscosity and reprocessing transformation functions to create multiwavelength light curves. It then uses an MCMC fitting routine to compare these theoretical light curves with observations. We find that none of the events show evidence for viscous delays exceeding a few days, supporting the theory that our current observing strategies in the optical/UV are missing a significant number of viscously delayed flares. We find that the events have black hole masses of 106–108 M⊙ and that the masses we predict are as reliable as those based on bulk galaxy properties. We also find that there is a preference for stars with mass <1 M⊙, as expected when low-mass stars greatly outnumber high-mass stars. [ABSTRACT FROM AUTHOR]

We use a simple framework to calculate the time evolution of the composition of the fallback material onto a supermassive black hole arising from the tidal disruption of main-sequence stars. We study stars with masses between 0.8 and 3.0 M⊙, at evolutionary stages from zero-age main sequence to terminal-age main sequence, built using the Modules for Experiments in Stellar Astrophysics code. We show that most stars develop enhancements in nitrogen (14N) and depletions in carbon (12C) and oxygen (16O) over their lifetimes, and that these features are more pronounced for higher mass stars. We find that, in an accretion-powered tidal disruption flare, these features become prominent only after the time of peak of the fallback rate and appear at earlier times for stars of increasing mass. We postulate that no severe compositional changes resulting from the fallback material should be expected near peak for a wide range of stellar masses and, as such, are unable to explain the extreme helium-to-hydrogen line ratios observed in some TDEs. On the other hand, the resulting compositional changes could help explain the presence of nitrogen-rich features, which are currently only detected after peak. When combined with the shape of the light curve, the time evolution of the composition of the fallback material provides a clear method to help constrain the nature of the disrupted star. This will enable a better characterization of the event by helping break the degeneracy between the mass of the star and the mass of the black hole when fitting tidal disruption light curves. [ABSTRACT FROM AUTHOR]

We present a study of binary–single interactions with energy-loss terms such as tidal dissipation and gravitational-wave (GW) 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 relatively 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 depend on the initial orbital parameters, the mass hierarchy, the type of interacting object, and the energy dissipation mechanism. This formalism is applied here to study the formation of two-body GW captures, for which we estimate absolute and relative rates relevant to Advanced LIGO detections. It is shown that two-body GW captures should have compelling observational implications if a sizable fraction of detected compact binaries are formed via dynamical interactions. [ABSTRACT FROM AUTHOR]

We perform the first systematic study of 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 in 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 semimajor 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 (NSs). In this case, we find the rate of highly eccentric white dwarf—NS 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. [ABSTRACT FROM AUTHOR]

During a common envelope (CE) 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 CE, 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 CEs using the “wind tunnel” formalism developed by MacLeod et al. 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 CE is rare, and if it occurs, transitory. The implication for LIGO black hole binary assembly is that by avoiding strong accretion feedback, CE interactions should still result in the substantial orbital tightening needed to produce merging binaries. [ABSTRACT FROM AUTHOR]

We perform a comprehensive study of the X-ray emission from 70 transient sources that have been classified as tidal disruption events (TDEs) in the literature. We explore the properties of these candidates, using nearly three decades of X-ray observations to quantify their properties and characteristics. We find that the emission from X-ray TDEs increase by two to three orders of magnitude, compared to pre-flare constraints. These emissions evolve significantly with time, and decay with power-law indices that are typically shallower than the canonical t−5/3 decay law, implying that X-ray TDEs are viscously delayed. These events exhibit enhanced (relative to galactic) column densities and are quite soft in nature, with no strong correlation between the amount of detected soft and hard emission. At their peak, jetted events have an X-ray to optical ratio ≫1, whereas non-jetted events have a ratio ∼1, which suggests that these events undergo reprocessing at different rates. X-ray TDEs have long T90 values, consistent with what would be expected from a viscously driven accretion disk formed by the disruption of a main-sequence star by a black hole with a mass <107M⊙. The isotropic luminosities of X-ray TDEs are bimodal, such that jetted and non-jetted events are separated by a “reprocessing valley” that we suggest is naturally populated by optical/UV TDEs that most likely produce X-rays, but this emission is “veiled” from observations due to reprocessing. Our results suggest that non-jetted X-ray TDEs likely originate from partial disruptions and/or disruptions of low-mass stars. [ABSTRACT FROM AUTHOR]

In this paper, we model the observable signatures of tidal disruptions of white dwarf (WD) stars using massive black holes (MBHs) of moderate mass, ≈103–105M⊙. 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⊙ 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 light curves and spectra reminiscent of Type I supernovae. The properties are strongly viewing angle dependent, and key spectral signatures are ≈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 multi-wavelength 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 BHs, 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 nondetection 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. [ABSTRACT FROM AUTHOR]

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, but also 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 simulations. We study clusters initially composed of 105 or 2 × 105 stars with IMBHs of 75 and 150 solar masses, and we 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 an orbital semimajor axis at least three times tighter than the second-most-bound object over 90% of the time. These companionships have typical periods on the order of years and are subject to cycles of exchange and destruction. The most frequently observed, long-lived pairings persist for ∼107 years. 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. The 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. The 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 occur with promising frequency. [ABSTRACT FROM AUTHOR]

Neutron star (binary neutron star and neutron star–black hole) mergers are believed to produce short-duration gamma-ray bursts (GRBs). They are also believed to be the dominant source of gravitational waves to be detected by the advanced LIGO and advanced VIRGO and the dominant source of the heavy r-process elements in the universe. Whether or not these mergers produce short-duration GRBs depends sensitively on the fate of the core of the remnant (whether, and how quickly, it forms a black hole). In this paper, we combine the results of Newtonian merger calculations and equation of state studies to determine the fate of the cores of neutron star mergers. Using population studies, we can determine the distribution of these fates to compare to observations. We find that black hole cores form quickly only for equations of state that predict maximum non-rotating neutron star masses below 2.3–2.4 solar masses. If quick black hole formation is essential in producing GRBs, LIGO/Virgo observed rates compared to GRB rates could be used to constrain the equation of state for dense nuclear matter. [ABSTRACT FROM AUTHOR]

We investigate the production sites and the enrichment history of r-process elements in the Galaxy, as traced by the [Eu/Fe] ratio, using the high resolution, cosmological zoom-in simulation “Eris.” At z = 0, Eris represents a close analog to the Milky Way, making it the ideal laboratory to understand the chemical evolution of our Galaxy. Eris formally traces the production of oxygen and iron due to supernovae (SNe) Ia and SNe II. We include in post-processing the production of r-process elements from compact binary mergers. Unlike previous studies, we find that the nucleosynthetic products from compact binary mergers can be incorporated into stars of very low metallicity and at early times, even with a minimum delay time of 100 Myr. This conclusion is relatively insensitive to modest variations in the merger rate, minimum delay time, and the delay time distribution. By implementing a first-order prescription for metal mixing, we can further improve the agreement between our model and the data for the chemical evolution of both [α/Fe] and [Eu/Fe]. We argue that compact binary mergers could be the dominant source of r-process nucleosynthesis in the Galaxy. [ABSTRACT FROM AUTHOR]

We present light curves and spectra of the tidal disruption event (TDE) ASASSN-18pg/AT 2018dyb spanning a period of one year. The event shows a plethora of strong emission lines, including the Balmer series, He ii, He i, and metal lines of O iii λ3760 and N iii λλ4100, 4640 (blended with He ii). The latter lines are consistent with originating from the Bowen fluorescence mechanism. By analyzing literature spectra of past events, we conclude that these lines are common in TDEs. The spectral diversity of optical TDEs is thus larger than previously thought and includes N-rich events besides H- and He-rich events. We study how the spectral lines evolve with time, by means of their width, relative strength, and velocity offsets. The velocity width of the lines starts at ∼13,000 km s−1 and decreases with time. The ratio of He ii to N iii increases with time. The same is true for ASASSN-14li, which has a very similar spectrum to AT 2018dyb but its lines are narrower by a factor of >2. We estimate a black hole mass of MBH = M⊙ by using the M–σ relation. This is consistent with the black hole mass derived using the MOSFiT transient fitting code. The detection of strong Bowen lines in the optical spectrum is an indirect proof for extreme ultraviolet and (reprocessed) X-ray radiation and favors an accretion origin for the TDE optical luminosity. A model where photons escape after multiple scatterings through a super-Eddington thick disk and its optically thick wind, viewed at an angle close to the disk plane, is consistent with the observations. [ABSTRACT FROM AUTHOR]

The physics of core-collapse (CC) supernovae (SNe) and how the explosions depend on progenitor properties are central questions in astronomy. For only a handful of SNe, the progenitor star has been identified in pre-explosion images. Supernova remnants (SNRs), which are observed long after the original SN event, provide a unique opportunity to increase the number of progenitor measurements. Here we systematically examine the stellar populations in the vicinities of 23 known SNRs in the Small Magellanic Cloud (SMC) using the star formation history (SFH) maps of Harris & Zaritsky. We combine the results with constraints on the SNR metal abundances and environment from X-ray and optical observations. We find that 22 SNRs in the SMC have local SFHs and properties consistent with a CC explosion, several of which are likely to have been high-mass progenitors. This result supports recent theoretical findings that high-mass progenitors can produce successful explosions. We estimate the mass distribution of the CC progenitors and find that this distribution is similar to a Salpeter IMF (within the uncertainties), while this result is shallower than the mass distribution found in M31 and M33 by Jennings et al. and Díaz-Rodríguez et al. using a similar approach. Additionally, we find that a number of the SMC SNRs exhibit a burst of star formation between 50 and 200 Myr ago. As these sources are likely CC, this signature may be indicative of massive stars undergoing delayed CC as a consequence of binary interaction, rapid rotation, or low metallicity. In addition, the lack of Type Ia SNRs in the SMC is possibly a result of the short visibility times of these sources, as they may fall below the sensitivity limits of current radio observations. [ABSTRACT FROM AUTHOR]

We present the first systematic study of strong binary–single and binary–binary black hole (BH) interactions with the inclusion of general relativity. By including general relativistic effects in the equations of motion during strong encounters, the dissipation of orbital energy from the emission of gravitational waves (GWs) can lead to captures and subsequent inspirals with appreciable eccentricities when entering the sensitive frequency ranges of the LIGO and Virgo GW detectors. It has been shown that binary–single interactions significantly contribute to the rate of eccentric mergers, but no studies have looked exclusively into the contribution from binary–binary interactions. To this end, we perform binary–binary and binary–single scattering experiments with general relativistic dynamics up through the 2.5 post-Newtonian order included, both in a controlled setting to gauge the importance of non-dissipative post-Newtonian terms and derive scaling relations for the cross section of GW captures, as well as experiments tuned to the strong interactions from state-of-the art globular cluster (GC) models to assess the relative importance of the binary–binary channel in facilitating GW captures and the resultant eccentricity distributions of inspiral from channel. Although binary–binary interactions are 10–100 times less frequent in GCs than binary–single interactions, their longer lifetime and more complex dynamics leads to a higher probability for GW captures to occur during the encounter. We find that binary–binary interactions contribute 25%–45% of the eccentric mergers that occur during strong BH encounters in GCs, regardless of the properties of the cluster environment. The inclusion of higher multiplicity encounters in dense star clusters therefore have major implications on the predicted rates of highly eccentric binaries potentially detectable by the LIGO/Virgo network. Because gravitational waveforms of eccentric inspirals are distinct from those generated by merging binaries that have circularized, measurements of eccentricity in such systems would highly constrain their formation scenario. [ABSTRACT FROM AUTHOR]

Galaxy formation simulations demonstrate that cosmic-ray (CR) feedback may be important in the launching of galactic-scale winds. CR protons dominate the bulk of the CR population, yet most observational constraints of CR feedback come from synchrotron emission of CR electrons. In this paper, we analyze 105 months of Fermi Gamma-ray Space Telescope observations of the Small Magellanic Cloud (SMC), with the aim of exploring CR feedback and transport in an external galaxy. We produce maps of the 2–300 GeV emission and detect statistically significant extended emission along the “Bar” and the “Wing,” where active star formation is occurring. Gamma-ray emission is not detected above 13 GeV, and we set stringent upper limits on the flux above this energy. We find the best fit to the gamma-ray spectrum is a single-component model with a power law of index and an exponential cutoff energy of GeV. We assess the relative contribution of pulsars and CRs to the emission, and we find that pulsars may produce up to 10% of the flux above 100 MeV. Thus, we attribute most of the gamma-ray emission (based on its spectrum and morphology) to CR interactions with the interstellar medium. We show that the gamma-ray emissivity of the SMC is at least five times smaller than that of the Milky Way and that the SMC is far below the “calorimetric limit,” where all CR protons experience pion losses. We interpret these findings as evidence that CRs are escaping the SMC via advection and diffusion. [ABSTRACT FROM AUTHOR]

We present the initial–final mass relation (IFMR) based on the self-consistent analysis of Sirius B and 79 white dwarfs from 13 star clusters. We have also acquired additional signal on eight white dwarfs previously analyzed in the NGC 2099 cluster field, four of which are consistent with membership. These re-observed white dwarfs have masses ranging from 0.72 to 0.97 M⊙, with initial masses from 3.0 to 3.65 M⊙, where the IFMR has an important change in slope that these new data help to observationally confirm. In total, this directly measured IFMR has small scatter (σ = 0.06 M⊙) and spans from progenitors of 0.85 to 7.5 M⊙. Applying two different stellar evolutionary models to infer two different sets of white dwarf progenitor masses shows that, when the same model is also used to derive the cluster ages, the resulting IFMR has weak sensitivity to the adopted model at all but the highest initial masses (>5.5 M⊙). The nonlinearity of the IFMR is also clearly observed with moderate slopes at lower masses (0.08 Mfinal/Minitial) and higher masses (0.11 Mfinal/Minitial) that are broken up by a steep slope (0.19 Mfinal/Minitial) between progenitors from 2.85 to 3.6 M⊙. This IFMR shows total stellar mass loss ranges from 33% of Minitial at 0.83 M⊙ to 83% of Minitial at 7.5 M⊙. Testing this total mass loss for dependence on progenitor metallicity, however, finds no detectable sensitivity across the moderate range of −0.15 < [Fe/H] < +0.15. [ABSTRACT FROM AUTHOR]

Supernova explosions are inherently asymmetric and can accelerate new-born neutron stars (NSs) to hundreds of km s−1. Two prevailing theories to explain NS kicks are ejecta asymmetries (e.g., conservation of momentum between NS and ejecta) and anisotropic neutrino emission. Observations of supernova remnants (SNRs) can give us insights into the mechanism that generates these NS kicks. In this paper, we investigate the relationship between NS kick velocities and the X-ray morphologies of 18 SNRs observed with the Chandra X-ray Observatory and the Röntgen Satellite (ROSAT). We measure SNR asymmetries using the power-ratio method (a multipole expansion technique), focusing on the dipole, quadrupole, and octupole power ratios. Our results show no correlation between the magnitude of the power ratios and NS kick velocities, but we find that for Cas A and G292.0+1.8, whose emission traces the ejecta distribution, their NSs are preferentially moving opposite to the bulk of the X-ray emission. In addition, we find a similar result for PKS 1209–51, CTB 109, and Puppis A; however, their emission is dominated by circumstellar/interstellar material, so their asymmetries may not reflect their ejecta distributions. Our results are consistent with the theory that NS kicks are a consequence of ejecta asymmetries as opposed to anisotropic neutrino emission. In the future, additional observations to measure NS proper motions within ejecta-dominated SNRs are necessary to robustly constrain the NS kick mechanism. [ABSTRACT FROM AUTHOR]

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 105–—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 white dwarf by a 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. [ABSTRACT FROM AUTHOR]

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 , in the other, . 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, flow. We illustrate the impact of these newly derived coefficients by integrating the inspiral of a secondary object through the envelopes of () and () 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 focused 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. [ABSTRACT FROM AUTHOR]

We study the emission from tidal disruption events (TDEs) produced as radiation from black hole accretion propagates through an extended, optically thick envelope formed from stellar debris. We analytically describe key physics controlling spectrum formation, and present detailed radiative transfer calculations that model the spectral energy distribution and optical line strengths of TDEs near peak brightness. The steady-state transfer is coupled to a solver for the excitation and ionization states of hydrogen, helium, and oxygen (as a representative metal), without assuming local thermodynamic equilibrium. Our calculations show how an extended envelope can reprocess a fraction of soft X-rays and produce the observed optical fluxes of the order of 1043 erg s−1, with an optical/UV continuum that is not described by a single blackbody. Variations in the mass or size of the envelope may help explain how the optical flux changes over time with roughly constant color. For high enough accretion luminosities, X-rays can escape to be observed simultaneously with the optical flux. Due to optical depth effects, hydrogen Balmer line emission is often strongly suppressed relative to helium line emission (with He ii-to-H line ratios of at least 5:1 in some cases) even in the disruption of a solar-composition star. We discuss the implications of our results to understanding the type of stars destroyed in TDEs and the physical processes responsible for producing the observed flares. [ABSTRACT FROM AUTHOR]

We report the discovery of an off-nuclear ultrasoft hyper-luminous X-ray source candidate 3XMM J141711.1+522541 in the inactive S0 galaxy SDSS J141711.07+522540.8 (z = 0.41827, dL = 2.3 Gpc) in the Extended Groth Strip. It is located at a projected offset of ∼1.″0 (5.2 kpc) from the nucleus of the galaxy and was serendipitously detected in five XMM-Newton observations in 2000 July. Two observations have enough counts and can be fitted with a standard thermal disk with an apparent inner disk temperature and a 0.28–14.2 keV unabsorbed luminosity LX ∼ 4 × 1043 erg s−1 in the source rest frame. The source was still detected in three Chandra observations in 2002 August, with similarly ultrasoft but fainter spectra (kTMCD ∼ 0.17 keV, LX ∼ 0.5 × 1043 erg s−1). It was not detected in later observations, including two by Chandra in 2005 October, one by XMM-Newton in 2014 January, and two by Chandra in 2014 September–October, implying a long-term flux variation factor of >14. Therefore the source could be a transient with an outburst in 2000–2002. It has a faint optical counterpart candidate, with apparent magnitudes of mF606W = 26.3 AB mag and mF814W = 25.5 AB mag in 2004 December (implying an absolute V-band magnitude of ∼−15.9 AB mag). We discuss various explanations for the source and find that it is best explained as a massive black hole (BH) embedded in the nucleus of a possibly stripped satellite galaxy, with the X-ray outburst due to tidal disruption of a surrounding star by the BH. The BH mass is ∼105M⊙, assuming the peak X-ray luminosity at around the Eddington limit. [ABSTRACT FROM AUTHOR]

We observed a sample of 10 white dwarf candidates in the rich open cluster NGC 2323 (M50) with the Keck Low-Resolution Imaging Spectrometer. The spectroscopy shows eight to be DA white dwarfs, with six of these having high signal-to-noise ratio appropriate for our analysis. Two of these white dwarfs are consistent with singly evolved cluster membership, and both are high mass ∼1.07 M⊙, and give equivalent progenitor masses of 4.69 M⊙. To supplement these new high-mass white dwarfs and analyze the initial–final mass relation (IFMR), we also looked at 30 white dwarfs from publicly available data that are mostly all high-mass ( M⊙). These original published data exhibited significant scatter, and to test if this scatter is true or simply the result of systematics, we have uniformly analyzed the white dwarf spectra and have adopted thorough photometric techniques to derive uniform cluster parameters for their parent clusters. The resulting IFMR scatter is significantly reduced, arguing that mass-loss rates are not stochastic in nature and that within the ranges of metallicity and mass analyzed in this work mass loss is not highly sensitive to variations in metallicity. Lastly, when adopting cluster ages based on Y2 isochrones, the slope of the high-mass IFMR remains steep and consistent with that found from intermediate-mass white dwarfs, giving a linear IFMR from progenitor masses between 3 and 6.5 M⊙. In contrast, when adopting the slightly younger cluster ages based on PARSEC isochrones, the high-mass IFMR has a moderate turnover near an initial mass of 4 M⊙. [ABSTRACT FROM AUTHOR]

We have expanded the sample of observed white dwarfs in the rich open cluster NGC 2099 (M37) with the Keck Low Resolution Imaging Spectrometer. Of 20 white dwarf candidates, the spectroscopy shows 19 to be true white dwarfs, with 14 of these having high signal-to-noise ratio. We find 11 of these 14 to be consistent with singly evolved cluster members. They span a mass range of ∼0.7–0.95 , excluding a low-mass outlier, corresponding to progenitor masses of ∼3–4 . This region of the initial–final mass relation (IFMR) has large scatter and a slope that remains to be precisely determined. With this large sample of white dwarfs that belong to a single age and metallicity population, we find an IFMR of (0.171 ± 0.057)M, significantly steeper than the linear relation adopted over the full observed white dwarf mass range in many previous studies. Comparison of this new relation from the solar-metallicity NGC 2099 to 18 white dwarfs in the metal-rich Hyades and Praesepe shows that their IFMR also has a consistently steep slope. This strong consistency also suggests that there is no significant metallicity dependence of the IFMR at this mass and metallicity range. As a result, the IFMR can be more reliably determined with this broad sample of 29 total white dwarfs, giving M = (0.163 ± 0.022)M + 0.238 ± 0.071 from of 3–4 . A steep IFMR in this mass range indicates that the full IFMR is nonlinear. [ABSTRACT FROM AUTHOR]