Search

In the last few years, a mysterious new class of astrophysical objects has been uncovered. These are spatially coincident with the nuclei of external galaxies and show X-ray variations that repeat on timescales of minutes to a month. They manifest in three different ways in the data: stable quasi-periodic oscillations (QPOs), quasi-periodic eruptions (QPEs) and quasi-periodic outflows (QPOuts). QPOs are systems that show smooth recurrent X-ray brightness variations while QPEs are sudden changes that appear like eruptions. QPOuts represent systems that exhibit repeating outflows moving at mildly-relativistic velocities of about 0.1-0.3c, where c is the speed of light. Their underlying physical mechanism is a topic of heated debate, with most models proposing that they originate either from instabilities within the inner accretion flow or from orbiting objects. There is a huge excitement especially from the latter class of models as it has been argued that some repeating systems could host extreme mass-ratio inspirals, potentially detectable with upcoming space-based gravitational wave interferometers. Consequently, paving the path for an era of "persistent" multi-messenger astronomy. Here we summarize the recent findings on the topics, including the newest observational data, various physical models and their numerical implementation., Comment: 22 pages, 3 figures, summary of the Session MA3 submitted to Proceedings of the Seventeenth Marcel Grossmann Meeting, Pescara, Italy, 2024

Periodic collisions between a star on an inclined orbit around a supermassive black hole and its accretion disk offers a promising explanation for X-ray "quasi-periodic eruptions" (QPEs). Each passage through the disk shocks and compresses gas ahead of the star, which subsequently re-expands above the disk as a quasi-spherical cloud. We present spherically symmetric Monte Carlo radiation transport simulations which follow the production of photons behind the radiation-mediated shock, Comptonization by hot electrons, and the eventual escape of the radiation through the expanding debris. Such one-dimension calculations are approximately justified for thin disks, through which the star of radius $R_{\star}$ passes faster than the shocked gas can flow around the star. For collision speeds $v_{\rm coll} \gtrsim 0.15 c$ and disk surface densities $\Sigma \sim 10^{3}$ g cm$^{-2}$ characteristic of those encountered by stellar orbits consistent with QPE recurrence times, the predicted transient light curves exhibit peak luminosities $\gtrsim 10^{42}$ erg s$^{-1}$ and Comptonized quasi-thermal (Wien-like) spectra which peak at energies $h\nu \sim 100$ eV, broadly consistent with QPE properties. For these conditions, gas and radiation are out of equilibrium and the emission temperature is harder than the blackbody value due to inefficient photon production behind the shock. Alternatively, for higher disk densities and/or lower shock velocities, QPE emission could instead represent the comparatively brief phase shortly after shock break-out, though in this case the bulk of the radiation is thermalized and occurs in the ultraviolet instead of the X-ray band. In either scenario, reproducing the observed eruption properties (duration, luminosity, temperature) requires a large radius $R_{\star} \gtrsim 10R_{\odot}$, which may point to inflation of the star's atmosphere from repeated collisions., Comment: 34 pages, 20 figures. Submitted to ApJ

We consider Roche lobe overflow (RLO) from a low-mass star on a nearly circular orbit, onto a supermassive black hole (SMBH). If mass transfer is unstable, its rate accelerates in a runaway process, resulting in highly super-Eddington mass accretion rates, accompanied by an optically-thick outflow emanating from the SMBH vicinity. This produces a week-month long, bright optical/Ultraviolet flare, accompanied by a year-decade long X-ray precursor and post-cursor emitted from the accretion flow onto the SMBH. Such ``Circular Tidal Disruption Events (TDEs)" represent a new class of nuclear transients, occurring at up to $1-10\%$ of the canonical parabolic tidal disruption event rate. Near breakup rotation and strong tidal deformation of the star prior to disruption could lead to strong magnetic fields, making circular-TDEs possible progenitors of jetted TDEs. Outflows prior to the final stellar disruption produce a circum-nuclear environment (CNM) with $\sim \rm 10^{-2} \, M_\odot$ at distances of $\sim 0.01-0.1 \, \rm pc$, likely leading to bright radio emission, and also similar to the CNM inferred for jetted TDEs. We discuss broader connections between circular TDEs and other recently identified classes of transients associated with galactic nuclei, such as repeating-TDEs and Quasi-Periodic X-ray Eruptions, as well as possible connections to luminous fast blue optical transients such as AT2018cow. We also discuss observational signatures of the analogous RLO of a white dwarf around an intermediate mass BH, which may be a multi-messenger source in the LISA era., Comment: 22 pages, 7 figures. Submitted to ApJ. Comments are welcome

Supermassive black holes and their surrounding dense stellar environments nourish a variety of astrophysical phenomena. We focus on the distribution of stellar-mass black holes around the supermassive black hole and the consequent formation of extreme mass-ratio inspirals (EMRIs). We derive a steady-state distribution, considering the effects of two-body scatterings and gravitational wave emission, and calculate the EMRI formation rate, eccentricity distribution and EMRI-to-plunge ratio. Our model predicts: (I) A stronger segregation than previously estimated at the outskirts of the sphere of influence (at $\sim0.01\rm pc$ to $2\rm pc$ for a Milky Way-like galaxy). (II) An increased EMRI-to-plunge ratio, favoring EMRIs at galaxies where stellar-mass black holes are scarce. (III) A detection of about $2\times10^3$ resolvable EMRIs, with a signal-to-noise ratio above $20$, along a $4$ year LISA mission time. (IV) A confusion noise, induced by a cosmological population of unresolved EMRIs, reducing LISA sensitivity in the $1-5\ \rm mHz$ frequency range by up to a factor of $\approx2$, relative to the instrumental noise., Comment: 15 pages, 6 figures. Submitted to ApJ

A modest fraction of the stars in galactic nuclei fed towards the central supermassive black hole (SMBH) approach on low-eccentricity orbits driven by gravitational-wave radiation (extreme mass ratio inspiral, EMRI). In the likely event that a gaseous accretion disk is created in the nucleus during this slow inspiral (e.g., via an independent tidal-disruption event; TDE), star-disk collisions generate regular short-lived flares consistent with the observed quasi-periodic eruption (QPE) sources. We present a model for the coupled star-disk evolution which self-consistently accounts for mass and thermal energy injected into the disk from stellar collisions and associated mass ablation. For weak collision/ablation heating, the disk is thermally-unstable and undergoes limit-cycle oscillations which modulate its properties and lead to accretion-powered outbursts on timescales of years to decades, with a time-averaged accretion rate $\sim 0.1 \dot{M}_{\rm Edd}$. Stronger collision/ablation heating acts to stabilize the disk, enabling roughly steady accretion at the EMRI-stripping rate. In either case, the stellar destruction time through ablation, and hence the maximum QPE lifetime, is $\sim 10^{2}-10^{3}$ yr, far longer than fall-back accretion after a TDE. The quiescent accretion disks in QPE sources may at the present epoch be self-sustaining and fed primarily by EMRI ablation. Indeed, the observed range of secular variability broadly match those predicted for collision-fed disks. Changes in the QPE recurrence pattern following such outbursts, similar to that observed in GSN 069, could arise from temporary misalignment between the EMRI-fed disk and the SMBH equatorial plane as the former regrows its mass after a state transition., Comment: 24 pages, 6 figures. Submitted to ApJ. Comments welcome!

Quasi-periodic Eruptions (QPEs) represent a novel class of extragalactic X-ray transients that are known to repeat at roughly regular intervals of a few hours to days. Their underlying physical mechanism is a topic of heated debate, with most models proposing that they originate either from instabilities within the inner accretion flow or from orbiting objects. At present, our knowledge of how QPEs evolve over an extended timescale of multiple years is limited, except for the unique QPE source GSN 069. In this study, we present results from strategically designed Swift observing programs spanning the past three years, aimed at tracking eruptions from eRO-QPE1. Our main results are: 1) the recurrence time of eruptions can vary between 0.6 and 1.2 days, 2) there is no detectable secular trend in evolution of the recurrence times, 3) consistent with prior studies, their eruption profiles can have complex shapes, and 4) the peak flux of the eruptions has been declining over the past 3 years with the eruptions barely detected in the most recent Swift dataset taken in June of 2023. This trend of weakening eruptions has been reported recently in GSN 069. However, because the background luminosity of eRO-QPE1 is below our detection limit, we cannot verify if the weakening is correlated with the background luminosity (as is claimed to be the case for GSN 069). We discuss these findings within the context of various proposed QPE models., Comment: Resubmitted to ApJ Letters after implementing referee comments

Wide-field survery have recently detected recurring optical and X-ray sources near galactic nuclei, with period spanning hours to years. These phenomena could result from repeated partial tidal disruptions of stars by supermassive black holes (SMBHs) or by interaction between star and SMBH-accretion discs. We study the physical processes that produce period changes in such sources, highlighting the key role of the interaction between the orbiting star and the accretion disc. We focus on ASASSN-14ko - a repeatedly flaring optical source with a mean period $P_0 = 115 \, \rm d$ and a detected period decay $\dot{P} = -2.6\times 10^{-3}$ (Payne et al. 2022). We argue that the system's $\dot{P}$ is most compatible with true orbital decay produced by hydrodynamical drag as a star passes through the accretion disc on an inclined orbit, twice per orbit. The star is likely a sun-like star whose envelope is somewhat inflated, possibly due to tidal heating. Star-disc interaction inevitably leads to drag-induced stripping of mass from the star, which may be the dominant component in powering the observed flares. We discuss ASASSN-14ko's possible formation history and observational tests of our interpretation of the measured $\dot P$. Our results imply that partial tidal disruption events manifesting as repeating nuclear transients cannot be modeled without accounting for the cumulative impact of tidal heating over many orbits. We discuss the implications of our results for other repeating transients, and predict that the recurrence time of Quasi-Periodic Eruptions is expected to decay at a rate of order $|\dot{P}| \approx 10^{-6}-10^{-5}$., Comment: 13 pages, 4 figures. Accepted for publication in MNRAS

We examine the effects of two-body interactions in a nuclear star cluster surrounding a supermassive black hole. We evaluate the energy flux, analogously to the particle flux calculation of Bahcall and Wolf (1976). We show that there are two types of power-law steady-state solutions: one with zero energy flux and constant particle flux and the other with constant energy flux and zero particle flux. We therefore prove that a zero particle flux solution, which corresponds to the case of an accreting supermassive black hole, can be obtained by requiring a constant energy flux. Consequently, this solution can be derived by simple dimensional analysis, bypassing the need for detailed calculation. Finally, we show that this characteristic, of zero particle flux for constant energy flux and vice versa, is not unique to the Keplerian potential of a supermassive black hole but holds for any central potential of the form $\phi\propto r^{-\beta}$., Comment: 6 pages, Accepted for publication in ApJ

Like most galaxies, the Milky Way harbors a supermassive black hole (SMBH) at its center, surrounded by a nuclear star cluster. In this dense star cluster, direct collisions can occur between stars before they evolve off the main-sequence. Using a statistical approach, we characterize the outcomes of these stellar collisions within the inner parsec of the Galactic Center (GC). Close to the SMBH, where the velocity dispersion is larger than the escape speed from a Sun-like star, collisions lead to mass loss. We find that the stellar population within $0.01$ pc is halved within about a Gyr because of destructive collisions. Additionally, we predict a diffuse population of peculiar low-mass stars in the GC. These stars have been divested of their outer layers in the inner $0.01$ pc before migrating to larger distances from the SMBH. Between $0.01$ and $0.1$ pc from the SMBH, collisions can result in mergers. Our results suggest that repeated collisions between lower mass stars can produce massive ($\gtrsim 10$ M$_\odot$) stars, and there may be $\sim 100$ of them residing in this region. We provide predictions on the number of G objects, dust and gas enshrouded stellar objects, that may result from main-sequence stellar collisions. Lastly, we comment on uncertainties in our model and possible connections between stellar collisions and the missing red giants in the GC., Comment: Submitted to ApJ. Comments welcome. 24 pages, 13 figures

Roughly half of the quasi-periodic eruption (QPE) sources in galactic nuclei exhibit a remarkably regular alternating "long-short'' pattern of recurrence times between consecutive flares. We show that a main-sequence star (brought into the nucleus as an extreme mass-ratio inspiral; EMRI) which passes twice per orbit through the accretion disk of the supermassive black-hole (SMBH) on a mildly eccentric inclined orbit, each time shocking and ejecting optically-thick gas clouds above and below the midplane, naturally reproduces observed properties of QPE flares. Inefficient photon production in the ejecta renders the QPE emission much harder than the blackbody temperature, enabling the flares to stick out from the softer quiescent disk spectrum. Destruction of the star via mass ablation limits the QPE lifetime to decades, precluding a long-lived AGN as the gaseous disk. By contrast, a tidal disruption event (TDE) naturally provides a transient gaseous disk on the requisite radial scale, with a rate exceeding the EMRI inward migration rate, suggesting that many TDEs should host a QPE. This picture is consistent with the X-ray TDE observed several years prior to the QPE appearance from GSN 069. Remarkably, a second TDE-like flare was observed from this event, starting immediately after detectable QPE activity ceased; this event could plausibly result from the (partial or complete) destruction of the QPE-generating star triggered by runaway mass-loss, though other explanations cannot be excluded. Our model can also be applied to black hole-disk collisions, such as those invoked in the context of the candidate SMBH binary OJ 287., Comment: 21 pages, 4 figures, 2 tables, Accepted to ApJ

We discuss the formation and evolution of systems comprised of a low-mass ($M_\star \lesssim 4 \, \rm M_\odot$) main sequence star, orbiting a $10^5-10^7 \, \rm M_\odot$ supermassive black hole with an orbital period of order $\sim$hours, and a mild eccentricity ($e\approx0.1-0.2$), episodically shedding mass at each pericenter passage. We argue that the resulting mass transfer is likely unstable, with Roche lobe overflow initially driven by gravitational wave emission, but then being accelerated by the star's expansion in response to its mass loss, undergoing a runaway process. We show that such systems are naturally produced by two-body gravitational encounters within the inner parsec of a galaxy, followed by gravitational wave circularization and inspiral from initially highly eccentric orbits. We argue that such systems can produce recurring flares similar to the recently identified class of X-ray transients known as Quasi-Periodic Eruptions, observed at the centers of a few distant galaxies., Comment: 13 pages, 2 figures. Submitted to ApJ. Comments are welcome

Quasi-Periodic Erupters (QPEs) are a remarkable class of objects exhibiting very large amplitude quasi-periodic X-ray flares. Although numerous dynamical models have been proposed to explain them, relatively little attention has been given to using the properties of their radiation to constrain their dynamics. Here we show that the observed luminosity, spectrum, repetition period, duty cycle, and fluctuations in the latter two quantities point toward a model in which: a main sequence star on a moderately eccentric orbit around a supermassive black hole periodically transfers mass to the Roche lobe of the black hole; orbital dynamics lead to mildly-relativistic shocks near the black hole; and thermal X-rays at the observed temperature are emitted by the gas as it flows away from the shock. Strong X-ray irradiation of the star by the flare itself augments the mass transfer, creates fluctuations in flare timing, and stirs turbulence in the stellar atmosphere that amplifies magnetic field to a level at which magnetic stresses can accelerate infall of the transferred mass toward the black hole., Comment: 15 pgs., submitted to Ap. J

We present a new analytical solution to the steady-state distribution of stars close to a central supermassive black hole of mass $M_{\bullet}$ in the center of a galaxy. Assuming a continuous mass function of the form $dN/dm \propto m^{\gamma}$, stars with a specific orbital energy $x = GM_{\bullet}/r - v^2/2$ are scattered primarily by stars of mass $m_{\rm d}(x) \propto x^{-5/(4\gamma+10)}$ that dominate the scattering of both lighter and heavier species at that energy. Stars of mass $m_{\rm d}(x)$ are exponentially rare at energies lower than $x$, and follow a density profile $n(x') \propto x'^{3/2}$ at energies $x' \gg x$. Our solution predicts a negligible flow of stars through energy space for all mass species, similarly to the conclusions of Bahcall & Wolf (1977), but in contrast to the assumptions of Alexander & Hopman (2009). This is the first analytic solution which smoothly transitions between regimes where different stellar masses dominate the scattering., Comment: 8 pages, 1 figure. Submitted to ApJ

Most stellar evolution models predict that black holes (BHs) should not exist above approximately $50-70$ M$_\odot$, the lower limit of the pair-instability mass gap. However, recent LIGO/Virgo detections indicate the existence of BHs with masses at and above this threshold. We suggest that massive BHs, including intermediate mass black holes (IMBHs), can form in galactic nuclei through collisions between stellar-mass black holes and the surrounding main-sequence stars. Considering dynamical processes such as collisions, mass segregation, and relaxation, we find that this channel can be quite efficient, forming IMBHs as massive as $10^4$ M$_\odot$. This upper limit assumes that (1) the BHs accrete a substantial fraction of the stellar mass captured during each collision and (2) that the rate at which new stars are introduced into the region near the SMBH is high enough to offset depletion by stellar disruptions and star-star collisions. We discuss deviations from these key assumptions in the text. Our results suggest that BHs in the pair-instability mass gap and IMBHs may be ubiquitous in galactic centers. This formation channel has implications for observations. Collisions between stars and BHs can produce electromagnetic signatures, for example, from x-ray binaries and tidal disruption events. Additionally, formed through this channel, both black holes in the mass gap and IMBHs can merge with the supermassive black hole at the center of a galactic nucleus through gravitational waves. These gravitational wave events are extreme and intermediate mass ratio inspirals (EMRIs and IMRIs, respectively).

The shock breakout emission is the first light that emerges from a supernova. In the spherical case it is characterized by a brief UV flash. In an axisymmetric, non-spherical prolate explosion, the shock first breaches the surface along the symmetry axis, then peels around to larger angles, producing a breakout light curve which may differ substantially from the spherically symmetric case. We study the emergence of a non-relativistic, bipolar shock from a spherical star, and estimate the basic properties of the associated bolometric shock breakout signal. We identify four possible classes of breakout light curves, depending on the degree of asphericity. Compared to spherical breakouts, we find that the main distinguishing features of significantly aspherical breakouts are 1) a longer and fainter initial breakout flash and 2) an extended phase of slowly-declining, or even rising, emission which is produced as ejecta flung out by the oblique breakout expand and cool. We find that the breakout flash has a maximum duration of roughly $\sim R_*/v_{\rm bo}$ where $R_*$ is the stellar radius and $v_{\rm bo}$ is the velocity of the fastest-moving ejecta. For a standard Wolf--Rayet progenitor, the duration of the X-ray flash seen in SN 2008D exceeds this limit, and the same holds true for the prompt X-ray emission of low-luminosity GRBs such as GRB 060218. This suggests that these events cannot be explained by an aspherical explosion within a typical Wolf--Rayet star, implying that they originate from non-standard progenitors with larger breakout radii., Comment: 21 pages, 9 figures, accepted for publication in MNRAS

Many massive stars appear to undergo enhanced mass loss during late stages of their evolution. In some cases, the ejected mass likely originates from non-terminal explosive outbursts, rather than continuous winds. Here we study the dependence of the ejecta mass, $m_{\rm ej}$, on the energy budget $E$ of an explosion deep within the star, using both analytical arguments and numerical hydrodynamics simulations. Focusing on polytropic stellar models, we find that for explosion energies smaller than the stellar binding energy, the ejected mass scales as $m_{\rm ej} \propto E^{\varepsilon_{m}}$, where $\varepsilon_m = 2.4-3.0$ depending on the polytropic index. The loss of energy due to shock breakout emission near the stellar edge leads to the existence of a minimal mass-shedding explosion energy, corresponding to a minimal ejecta mass. For a wide range of progenitors, from Wolf-Rayet stars to red supergiants, we find a similar limiting energy of $E_{\rm min} \approx 10^{46}-10^{47} \rm \, erg$, almost independent of the stellar radius. The corresponding minimal ejecta mass varies considerably across different progenitors, ranging from $\sim \! 10^{-8} \, \rm M_\odot$ in compact stars, up to $\sim \! 10^{-2} \, \rm M_\odot$ in red supergiants. We discuss implications of our results for pre-supernova outbursts driven by wave heating, and complications caused by the non-constant opacity and adiabatic index of realistic stars., Comment: 9 pages, 5 figures. Submitted to MNRAS

The emergence of a shock from a medium with a free surface is an important process in various astrophysical phenomena. It generates the first light associated with explosions like supernovae and Gamma Ray Bursts. Most previous works considered planar or spherical geometries, where the shock front is parallel to the surface, and emerges simultaneously from all points. Here we study the hydrodynamics of an oblique planar shock breaking out from the planar surface of a uniform density ideal gas with adiabatic index $\gamma$. We obtain an analytic solution to the flow as a function of the angle between the plane of the shock and the surface $\beta$. We find steady state solutions (in a frame moving with the intersection point of the shock and the surface) up to some critical angle ($\beta_{max}=63.4$ degress for $\gamma=5/3$ and $\beta_{max}=69.3$ degrees for $\gamma=4/3$). We show how this analytic solution can be used in more complicated geometries where the shock is not planar, giving the exact profile of the outermost breakout ejecta. We apply our analytical results to a few realistic problems, such as underwater explosions, detonation under the surface of an asteroid, or off center detonations in a uniform sphere., Comment: 10 pages, 10 figures. Accepted for publication in Physics of Fluids

GW170817 was the first detection of a binary neutron star merger via gravitational waves. The event was observed over a wide range of the electromagnetic spectrum, revealing a thermal kilonova dominating the optical signal during the first $\sim$15 days, and a non-thermal synchrotron emission that has continued to rise $\sim$200 days post-merger, dominating the radio and x-ray emission. At early times, when the kilonova is still dominant, the synchrotron emitting electrons can efficiently cool by up-scattering the kilonova photos through inverse-Compton. Yet, the cooling frequency is not observed up to the X-ray band. This can only be explained if the source is moving at least at a mildly relativistic velocity. We find a lower limit on the source's bulk Lorentz factor of $\Gamma > 2.1$ at $9$ days. This lower limit is model independent, and relies directly on the observed quantities, providing an additional robust evidence to the relativistic motion in this event at early times., Comment: 4 pages, accepted for publishing in MNRAS

Transit timing variations (TTVs) are a powerful tool for characterizing the properties of transiting exoplanets. However, inferring planet properties from the observed timing variations is a challenging task, which is usually addressed by extensive numerical searches. We propose a new, computationally inexpensive method for inverting TTV signals in a planetary system of two transiting planets. To the lowest order in planetary masses and eccentricities, TTVs can be expressed as a linear combination of 3 functions, which we call the \textit{TTV modes}. These functions depend only on the planets' linear ephemerides, and can be either constructed analytically, or by performing 3 orbital integrations of the three-body system. Given a TTV signal, the underlying physical parameters are found by decomposing the data as a sum of the TTV modes. We demonstrate the use of this method by inferring the mass and eccentricity of 6 \textit{Kepler} planets that were previously characterized in other studies. Finally we discuss the implications and future prospects of our new method., Comment: Published on ApJ. Matched to published version

We consider stable mass transfer from the secondary to the primary of an extreme mass ratio binary system. We show that when the mass transfer is sufficiently fast, mass leakage occurs through the outer Lagrange point L2, in addition to the usual transfer through L1. We provide an analytical estimate for the mass leakage rate through L2 and find the conditions in which it is comparable to the mass transfer rate through L1. Focusing on a binary system of a main-sequence star and a super-massive black hole, driven by the emission of gravitational radiation, we show that it may sustain stable mass transfer, along with mass loss through L2. If such a mass-transferring system occurs at our Galactic Centre, it produces a gravitational wave signal detectable by future detectors, such as eLISA. The signal evolves according to the star's adiabatic index and cooling time. For low mass stars, the evolution is faster than the Kelvin-Helmholtz cooling rate driving the star out of the main-sequence. In some cases, the frequency and amplitude of the signal may both decrease with time, contrary to the standard chirp of a coalescing binary. Mass loss through L2, when occurs, decreases the evolution timescale of the emitted gravitational wave signal by up to a few tens of per cent. We conclude that L2 mass ejection is a crucial factor in analyzing gravitational waves signals produced by such systems., Comment: 15 pages, 4 figures, accepted for publication in MNRAS