-
Conditions for Changing-Look AGNs from Accretion Disk-Induced Tidal Disruption Events
Authors:
Yihan Wang,
Matthew J. Graham,
K. E. Saavik Ford,
Barry McKernan,
Taeho Ryu,
Daniel Stern
Abstract:
The phenomenon of changing-look (CL) behavior in active galactic nuclei (AGN) is characterized by dramatic changes in luminosity and/or emission line profiles over relatively short periods, ranging from months to years. The origin of CL-AGNs remains a mystery, but one proposed explanation involves the response of the inner AGN disk to tidal disruption events (TDEs) around the supermassive black ho…
▽ More
The phenomenon of changing-look (CL) behavior in active galactic nuclei (AGN) is characterized by dramatic changes in luminosity and/or emission line profiles over relatively short periods, ranging from months to years. The origin of CL-AGNs remains a mystery, but one proposed explanation involves the response of the inner AGN disk to tidal disruption events (TDEs) around the supermassive black hole (SMBH). In this Letter, we calculate the predicted frequency of AGN TDEs as a function of SMBH mass and compare the results to the observed CL-AGN distribution. We find that if the fraction of CL-AGNs caused by AGN-TDEs is high, then: (1) most SMBHs in CL-AGN are near maximal spin, with the dimensionless spin parameter $a>0.9$; (2) AGN inner disks have a high surface density ($\geq 10^{7}\, {\rm g\, cm^{-2}}$); (3) typical AGN lifetimes are $\sim 10$-$100$ Myr; and (4) a nuclear star cluster initial mass function (IMF) that scales as $\sim m_*^{-1.6}$ is preferred. Future observations of CL-AGN will help constrain the fraction of CL-AGNs caused by AGN-TDEs, SMBH spins, AGN lifetimes, and the nuclear star cluster IMF.
△ Less
Submitted 17 June, 2024;
originally announced June 2024.
-
Repeating partial disruptions and two-body relaxation
Authors:
Luca Broggi,
Nicholas C. Stone,
Taeho Ryu,
Elisa Bortolas,
Massimo Dotti,
Matteo Bonetti,
Alberto Sesana
Abstract:
Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only…
▽ More
Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only liberate a small amount of stellar mass, but may significantly alter the orbit of the star. We study the phenomenon of repeating partial tidal disruptions (pTDEs) by building a semi-analytical model that accounts for mass loss and tidal excitation. In the empty loss cone regime, where two-body relaxation is weak, we estimate the number of consecutive partial disruptions that a star can undergo, on average, before being significantly affected by two-body encounters. We find that in this empty loss cone regime, a star will be destroyed in a sequence of weak pTDEs, possibly explaining the tension between the low observed TDE rate and its higher theoretical estimates.
△ Less
Submitted 19 June, 2024; v1 submitted 8 April, 2024;
originally announced April 2024.
-
Ephemeral Flaring Transients Following Supernova Explosions in Black-Hole Binary Systems
Authors:
Davide Lazzati,
Rosalba Perna,
Taeho Ryu
Abstract:
The majority of massive stars are members of binary systems, a fraction of which is expected to remain bound after the first star goes off as a supernova. When the second star also explodes, the SN ejecta are bound to interact with the compact object companion. We explore the consequences of this interaction in the case of a black hole (BH) companion. We show that accretion of the SN ejecta by the…
▽ More
The majority of massive stars are members of binary systems, a fraction of which is expected to remain bound after the first star goes off as a supernova. When the second star also explodes, the SN ejecta are bound to interact with the compact object companion. We explore the consequences of this interaction in the case of a black hole (BH) companion. We show that accretion of the SN ejecta by the BH generally occurs with highly super-Eddington rates, which can be strongly modulated if the ejecta are clumpy, as typically observed in supernova remnants. This late accretion produces transient flaring activity with a time delay from the SN explosion of days to months, depending on the orbital separation and the velocity of the ejecta. The flares are expected to have a non-thermal, broad band spectrum, but their high-frequency emission (UV and X-rays) would be absorbed within the remnant. Observed flares should therefore be in the optical and near-infrared. We propose this model as an explanation to the late-time flaring activity observed in the fast blue optical transient AT2022tsd.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
Tidal Disruption Encores
Authors:
Taeho Ryu,
Rosalba Perna,
Matteo Cantiello
Abstract:
Nuclear star clusters (NSCs), made up of a dense concentrations of stars and the compact objects they leave behind, are ubiquitous in the central regions of galaxies, surrounding the central supermassive black hole (SMBH). Close interactions between stars and stellar-mass black holes (sBH) lead to tidal disruption events (TDEs). We uncover an interesting new phenomenon: For a subset of these, the…
▽ More
Nuclear star clusters (NSCs), made up of a dense concentrations of stars and the compact objects they leave behind, are ubiquitous in the central regions of galaxies, surrounding the central supermassive black hole (SMBH). Close interactions between stars and stellar-mass black holes (sBH) lead to tidal disruption events (TDEs). We uncover an interesting new phenomenon: For a subset of these, the unbound debris (to the sBH) remain bound to the SMBH, accreting at a later time, and thus giving rise to a second flare. We compute the rate of such events, and find them ranging within $10^{-6}$ -$10^{-3}$ yr$^{-1}$gal$^{-1}$ for SMBH mass $\simeq 10^{6}-10^{9}M_\odot$. Time delays between the two flares spread over a wide range, from less than a year to hundreds of years. The temporal evolution of the light curves of the second flare can vary between the standard $t^{-5/3}$ power-law to much steeper decays, providing a natural explanation for observed light curves in tension with the classical TDE model. Our predictions have implications for learning about NSC properties and calibrating its sBH population. Some double flares may be electromagnetic counterparts to LISA Extreme-Mass-Ratio-Inspiral (EMRI) sources. Another important implication is the possible existence of TDE-like events in very massive SMBHs, where TDEs are not expected. Such flares can affect spin measurements relying on TDEs in the upper SMBH range.
△ Less
Submitted 27 March, 2024; v1 submitted 23 February, 2024;
originally announced February 2024.
-
Strong lensing of tidal disruption events: Detection rates in imaging surveys
Authors:
K. Szekerczes,
T. Ryu,
S. H. Suyu,
S. Huber,
M. Oguri,
L. Dai
Abstract:
Tidal disruption events (TDEs) are multi-messenger transients in which a star is tidally destroyed by a supermassive black hole at the center of galaxies. The Rubin Observatory Legacy Survey of Space and Time (LSST) is anticipated to annually detect hundreds to thousands of TDEs, such that the first gravitationally lensed TDE may be observed in the coming years. Using Monte-Carlo simulations, we q…
▽ More
Tidal disruption events (TDEs) are multi-messenger transients in which a star is tidally destroyed by a supermassive black hole at the center of galaxies. The Rubin Observatory Legacy Survey of Space and Time (LSST) is anticipated to annually detect hundreds to thousands of TDEs, such that the first gravitationally lensed TDE may be observed in the coming years. Using Monte-Carlo simulations, we quantify the rate of both unlensed and lensed TDEs as a function of limiting magnitudes in four different optical bands ($u$, $g$, $r$, and $i$) for a range of TDE temperatures that match observations. Dependent on the temperature and luminosity model, we find that $g$ and $r$ bands are the most promising bands with unlensed TDE detections that can be as high as ${\sim}10^{4}$ annually. By populating a cosmic volume with realistic distributions of TDEs and galaxies that can act as gravitational lenses, we estimate that a few lensed TDEs (depending on the TDE luminosity model) can be detected annually in $g$ or $r$ bands in the LSST survey, with TDE redshifts in the range of ${\sim}0.5$ to ${\sim}2$. The ratio of lensed to unlensed detections indicates that we may detect ${\sim}1$ lensed event for every $10^{4}$ unlensed events, which is independent of the luminosity model. The number of lensed TDEs decreases as a function of the image separations and time delays, and most of the lensed TDE systems are expected to have image separations below ${\sim}3"$ and time delays within ${\sim}30$ days. At fainter limiting magnitudes, the $i$ band becomes notably more successful. These results suggest that strongly lensed TDEs are likely to be observed within the coming years and such detections will enable us to study the demographics of black holes at higher redshifts through the lensing magnifications.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
Rapid evolution of the recurrence time in the repeating partial tidal disruption event eRASSt J045650.3-203750
Authors:
Zhu Liu,
Taeho Ryu,
A. J. Goodwin,
A. Rau,
D. Homan,
M. Krumpe,
A. Merloni,
I. Grotova,
G. E. Anderson,
A. Malyali,
J. C. A. Miller-Jones
Abstract:
In this letter, we present the results from further X-ray and UV observations of the nuclear transient eRASSt J045650.3-203750 (hereafter J0456-20). We detected five repeating X-ray and UV flares from J0456-20, making it one of the most promising repeating partial tidal disruption event (pTDE) candidates. More importantly, we also found rapid changes in the recurrence time $T_\text{recur}$ of the…
▽ More
In this letter, we present the results from further X-ray and UV observations of the nuclear transient eRASSt J045650.3-203750 (hereafter J0456-20). We detected five repeating X-ray and UV flares from J0456-20, making it one of the most promising repeating partial tidal disruption event (pTDE) candidates. More importantly, we also found rapid changes in the recurrence time $T_\text{recur}$ of the X-ray flares by modelling the long-term X-ray light curve of J0456-20. $T_\text{recur}$ first decreased rapidly from about 300 days to around 230 days. It continued to decrease to around 190 days with an indication of a constant $T_\text{recur}$ evidenced from the latest three cycles. Our hydrodynamic simulations suggest that, in the repeating pTDE scenario, such rapid evolution of $T_\text{recur}$ could be reproduced if the original star is a $1~\mathrm{M}_\odot$ main-sequence star near the terminal age and loses nearly 80-90% of its mass during the initial encounter with a supermassive black hole (SMBH) of mass around $10^5~\mathrm{M}_\odot$. The inferred mass loss of 0.8-0.9 $\mathrm{M}_\odot$ is higher than the estimated value of around 0.13 $\mathrm{M}_\odot$ from observation, which could be explained if the radiation efficiency is low (i.e. $\ll0.1$). Our results indicate that repeating pTDEs could be effective tools to explore the dynamics around supermassive black holes beyond our own Galaxy.
△ Less
Submitted 4 March, 2024; v1 submitted 25 January, 2024;
originally announced January 2024.
-
Solar evolution models with a central black hole
Authors:
Earl P. Bellinger,
Matt E. Caplan,
Taeho Ryu,
Deepika Bollimpalli,
Warrick H. Ball,
Florian Kühnel,
R. Farmer,
S. E. de Mink,
Jørgen Christensen-Dalsgaard
Abstract:
Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here…
▽ More
Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here we compute the evolution of stars having such a BH at their center. We find that such objects can be surprisingly long-lived, with the lightest black holes having no influence over stellar evolution, while more massive ones consume the star over time to produce a range of observable consequences. Models of the Sun born about a BH whose mass has since grown to approximately $10^{-6}~\rm{M}_\odot$ are compatible with current observations. In this scenario, the Sun would first dim to half its current luminosity over a span of 100 Myr as the accretion starts to generate enough energy to quench nuclear reactions. The Sun would then expand into a fully-convective star, where it would shine luminously for potentially several Gyr with an enriched surface helium abundance, first as a sub-subgiant star, and later as a red straggler, before becoming a sub-solar-mass BH. We also present results for a range of stellar masses and metallicities. The unique internal structures of stars harboring BHs may make it possible for asteroseismology to discover them, should they exist. We conclude with a list of open problems and predictions.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Multi-messenger astronomy with black holes: tidal disruption events
Authors:
Thomas Wevers,
Taeho Ryu
Abstract:
This chapter provides an overview of tidal disruption events, aiming to provide an overview of both the theoretical and the observational state of the field, with the overarching goal of introducing them as tools to indirectly observe massive black holes in the Universe. We start by introducing the relevant theoretical concepts, physical scales and timescales with an emphasis on the classical fram…
▽ More
This chapter provides an overview of tidal disruption events, aiming to provide an overview of both the theoretical and the observational state of the field, with the overarching goal of introducing them as tools to indirectly observe massive black holes in the Universe. We start by introducing the relevant theoretical concepts, physical scales and timescales with an emphasis on the classical framework and how this has been (and continues to be) improved since the inception of the field. We then cover the current and future prospects of observing TDEs through a variety of messengers, including photons across the electromagnetic spectrum, as well as gravitational waves and neutrino particles. More recent advancements in the field, including repeating TDEs as well as TDEs by stellar-mass black holes, are also highlighted.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
Simulating the tidal disruption of stars by stellar-mass black holes using moving-mesh hydrodynamics
Authors:
Pavan Vynatheya,
Taeho Ryu,
Ruediger Pakmor,
Selma E. de Mink,
Hagai B. Perets
Abstract:
In the centers of dense star clusters, close encounters between stars and compact objects are likely to occur. We study tidal disruption events of main-sequence (MS) stars by stellar-mass black holes (termed $μ$TDEs), which can shed light on the processes occurring in these clusters, including being an avenue in the mass growth of stellar-mass BHs. Using the moving-mesh hydrodynamics code \texttt{…
▽ More
In the centers of dense star clusters, close encounters between stars and compact objects are likely to occur. We study tidal disruption events of main-sequence (MS) stars by stellar-mass black holes (termed $μ$TDEs), which can shed light on the processes occurring in these clusters, including being an avenue in the mass growth of stellar-mass BHs. Using the moving-mesh hydrodynamics code \texttt{AREPO}, we perform a suite of hydrodynamics simulations of partial $μ$TDEs of realistic, \texttt{MESA}-generated MS stars by varying the initial mass of the star ($0.5\,{\rm M}_{\rm \odot}$ and $1\,{\rm M}_{\rm \odot}$), the age of the star (zero-age, middle-age and terminal-age), the mass of the black hole ($10\,{\rm M}_{\rm \odot}$ and $40\,{\rm M}_{\rm \odot}$) and the impact parameter (yielding almost no mass loss to full disruption). We then examine the dependence of the masses, spins, and orbital parameters of the partially disrupted remnant on the initial encounter parameters. We find that the mass lost from a star decreases exponentially with increasing distance of approach and that a $1\,{\rm M}_{\rm \odot}$ star loses lesser mass than a $0.5\,{\rm M}_{\rm \odot}$. Moreover, a more evolved star is less susceptible to mass loss. Tidal torques at the closest approach spin up the remnant by factors of $10^2$--$10^4$ depending on the impact parameter. The remnant star can be bound (eccentric) or unbound (hyperbolic) to the black hole: hyperbolic orbits occur when the star's central density concentration is relatively low and the black hole-star mass ratio is high, which is the case for the disruption of a $0.5\,{\rm M}_{\rm \odot}$ star. Finally, we provide best-fit analytical formulae for a range of parameters that can be incorporated into cluster codes to model star-black hole interaction more accurately.
△ Less
Submitted 4 March, 2024; v1 submitted 23 October, 2023;
originally announced October 2023.
-
The look of high-velocity red-giant star collisions
Authors:
Luc Dessart,
Taeho Ryu,
Pau Amaro Seoane,
Andrew M. Taylor
Abstract:
High-velocity stellar collisions driven by a supermassive black hole (BH) or BH-driven disruptive collisions, in dense, nuclear clusters can rival the energetics of supergiant star explosions following gravitational collapse of their iron core. Here, starting from a sample of red-giant star collisions simulated with the hydrodynamics code AREPO, we generate photometric and spectroscopic observable…
▽ More
High-velocity stellar collisions driven by a supermassive black hole (BH) or BH-driven disruptive collisions, in dense, nuclear clusters can rival the energetics of supergiant star explosions following gravitational collapse of their iron core. Here, starting from a sample of red-giant star collisions simulated with the hydrodynamics code AREPO, we generate photometric and spectroscopic observables using the nonlocal thermodynamic equilibrium time-dependent radiative transfer code CMFGEN. Collisions from more extended giants or stronger collisions (higher velocity or smaller impact parameter) yield bolometric luminosities on the order of 1e43 erg/s at 1d, evolving on a timescale of a week to a bright plateau at ~1e41 erg/s, before plunging precipitously after 20-40d at the end of the optically-thick phase. This luminosity falls primarily in the UV in the first days, thus when it is at its maximum, and shifts to the optical thereafter. Collisions at lower velocity or from less extended stars produce ejecta that are fainter but may remain optically thick for up to 40d if they have a small expansion rate. These collision debris show a similar spectral evolution as that observed or modeled for blue-supergiant star explosions of massive stars, differing only in the more rapid transition to the nebular phase. Such BH-driven disruptive collisions should be detectable by high-cadence surveys in the UV like ULTRASAT.
△ Less
Submitted 10 October, 2023;
originally announced October 2023.
-
In-plane Tidal Disruption of Stars in Disks of Active Galactic Nuclei
Authors:
Taeho Ryu,
Barry McKernan,
Saavik Ford,
Matteo Cantiello,
Matthew Graham,
Daniel Stern,
Nathan W. C Leigh
Abstract:
Stars embedded in active galactic nucleus (AGN) disks or captured by them may scatter onto the supermassive black hole (SMBH), leading to a tidal disruption event (TDE). Using the moving-mesh hydrodynamics simulations with {\small AREPO}, we investigate the dependence of debris properties in in-plane TDEs in AGN disks on the disk density and the orientation of stellar orbits relative to the disk g…
▽ More
Stars embedded in active galactic nucleus (AGN) disks or captured by them may scatter onto the supermassive black hole (SMBH), leading to a tidal disruption event (TDE). Using the moving-mesh hydrodynamics simulations with {\small AREPO}, we investigate the dependence of debris properties in in-plane TDEs in AGN disks on the disk density and the orientation of stellar orbits relative to the disk gas (pro- and retro-grade). Key findings are: 1) Debris experiences continuous perturbations from the disk gas, which can result in significant and continuous changes in debris energy and angular momentum compared to `naked' TDEs. 2) Above a critical density of a disk around a SMBH with mass $M_{\bullet}$ ($ρ_{\rm crit} \sim 10^{-8}{\rm g~cm^{-3}}(M_{\bullet}/10^{6}{\rm M}_{\odot})^{-2.5}$) for retrograde stars, both bound and unbound debris is fully mixed into the disk. The density threshold for no bound debris return, inhibiting the accretion component of TDEs, is $ρ_{\rm crit,bound} \sim 10^{-9}{\rm g~cm^{-3}}(M_{\bullet}/10^{6}{\rm M}_{\odot})^{-2.5}$. 3) Observationally, AGN-TDEs transition from resembling naked TDEs in the limit of $ρ_{\rm disk}\lesssim 10^{-2}ρ_{\rm crit,bound}$ to fully muffled TDEs with associated inner disk state changes at $ρ_{\rm disk}\gtrsimρ_{\rm crit,bound}$, with a superposition of AGN+TDE in between. Stellar or remnant passages themselves can significantly perturb the inner disk. This can lead to an immediate X-ray signature and optically detectable inner disk state changes, potentially contributing to the changing-look AGN phenomenon. 4) Debris mixing can enriches the average disk metallicity over time if the star's metallicity exceeds that of the disk gas.
△ Less
Submitted 1 October, 2023;
originally announced October 2023.
-
Collisions of red giants in galactic nuclei
Authors:
Taeho Ryu,
Pau Amaro Seoane,
Andrew M. Taylor,
Sebastian T. Ohlmann
Abstract:
In stellar-dense environments, stars can collide with each other. For collisions close to a supermassive black hole (SMBH), the collisional kinetic energy can be so large that the colliding stars can be completely destroyed, potentially releasing an amount of energy comparable to that of a supernova. Such violent collisions, which we call BH-driven disruptive collisions (BDCs), have been examined…
▽ More
In stellar-dense environments, stars can collide with each other. For collisions close to a supermassive black hole (SMBH), the collisional kinetic energy can be so large that the colliding stars can be completely destroyed, potentially releasing an amount of energy comparable to that of a supernova. Such violent collisions, which we call BH-driven disruptive collisions (BDCs), have been examined mostly analytically, with the non-linear hydrodynamical effects being left largely unstudied. Using the moving-mesh hydrodynamics code {\small AREPO}, we investigate high-velocity ($>10^{3}$ km/s) collisions between 1M$_{\odot}$ giants with varying radii, impact parameters, and initial approaching velocities, and estimate their observables. Very strong shocks across the collision surface efficiently convert $\gtrsim10\%$ of the initial kinetic energy into radiation energy. The outcome is a gas cloud expanding supersonically, homologously, and quasi-spherically, generating a flare with a peak luminosity $\simeq 10^{41}-10^{44}$ erg/s in the extreme UV band ($\simeq 10$ eV). The luminosity decreases approximately following a power-law $t^{-0.7}$ initially, then $t^{-0.4}$ after $t\simeq$10 days at which point it would be bright in the optical band ($\lesssim 1$eV). Subsequent, and possibly even brighter, emission would be generated due to the accretion of the gas cloud onto the nearby SMBH, possibly lasting up to multi-year timescales. This inevitable BH-collision product interaction can contribute to the growth of BHs at all mass scales, in particular, seed BHs at high redshifts. Furthermore, the proximity of the events to the central BH makes them a potential tool for probing the existence of dormant BHs, even very massive ones which cannot be probed by tidal disruption events.
△ Less
Submitted 26 January, 2024; v1 submitted 14 July, 2023;
originally announced July 2023.
-
Close Encounters of Star - Black Hole Binaries with Single Stars
Authors:
Taeho Ryu,
Selma de Mink,
Rob Farmer,
Ruediger Pakmor,
Rosalba Perna,
Volker Springel
Abstract:
Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single st…
▽ More
Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single stars, focusing on the identification of final outcomes and their long-term evolution and observational properties, using the moving-mesh hydrodynamics code AREPO. This type of encounters produces five types of outcomes: stellar disruption, stellar collision, weak perturbation of the original binary, binary member exchange, and triple formation. The two decisive parameters are the binary phase angle, which determines which two objects meet at the first closest approach, and the impact parameter, which sets the boundary between violent and non-violent interactions. When the impact parameter is smaller than the semimajor axis of the binary, tidal disruptions and star-BH collisions frequently occur when the BH and the incoming star first meet, while the two stars mostly merge when the two stars meet first instead. In both cases, the BHs accrete from an accretion disk at super-Eddington rates, possibly generating flares luminous enough to be observed. The stellar collision products either form a binary with the BH or remain unbound to the BH. Upon collision, the merged stars are hotter and larger than main sequence stars of the same mass at similar age. Even after recovering their thermal equilibrium state, stellar collision products, if isolated, would remain hotter and brighter than main sequence stars until becoming giants.
△ Less
Submitted 10 July, 2023; v1 submitted 6 July, 2023;
originally announced July 2023.
-
Accretion onto a Supermassive Black Hole Binary Before Merger
Authors:
Mark J. Avara,
Julian H. Krolik,
Manuela Campanelli,
Scott C. Noble,
Dennis Bowen,
Taeho Ryu
Abstract:
While supermassive binary black holes inspiral toward merger they may also experience significant accretion of matter from a surrounding disk. We study the dynamics of this system, simultaneously describing the evolving spacetime and magnetized plasma, and present the first relativistic calculation simulating two equal-mass, non-spinning black holes as they inspiral from a $20M$ ($G=c=1$) initial…
▽ More
While supermassive binary black holes inspiral toward merger they may also experience significant accretion of matter from a surrounding disk. We study the dynamics of this system, simultaneously describing the evolving spacetime and magnetized plasma, and present the first relativistic calculation simulating two equal-mass, non-spinning black holes as they inspiral from a $20M$ ($G=c=1$) initial separation almost to merger, $\simeq 9M$ ($M$=binary mass). Our dynamical results imply important observational consequences: for instance, the accretion rate $\dot M$ onto the black holes first decreases and then reaches a plateau, dropping by only a factor of $\sim 3$ despite the rapid inspiral. An estimated bolometric light curve thus suggests some merging SMBBHs may be quite luminous past the predicted decoupling from the circumbinary disk. The minidisks through which the accretion reaches the black holes are very non-standard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, ``sloshing" streams transfer mass from one minidisk to the other through the L1 point at a rate $\sim 0.1\times$ the accretion rate, carrying kinetic energy at a rate that can be as large as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives minidisks with time-varying tilts. The unsigned poloidal magnetic flux on the black hole event horizon is roughly constant at a dimensionless level $φ\sim 2-3$, but doubles just before merger; if the black holes had significant spin, this flux indicates the potential for powerful jets with variability driven by binary dynamics, another prediction of potentially unique EM signatures. This simulation is the first to employ our multipatch infrastructure \pwmhd, decreasing computational expense to $\sim 3\%$ of conventional single-grid methods' cost.
△ Less
Submitted 16 January, 2024; v1 submitted 29 May, 2023;
originally announced May 2023.
-
Shocks Power Tidal Disruption Events
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran,
Scott Noble,
Mark Avara
Abstract:
Accretion of debris seems to be the natural mechanism to power the radiation emitted during a tidal disruption event (TDE), in which a supermassive black hole tears apart a star. However, this requires the prompt formation of a compact accretion disk. Here, using a fully relativistic global simulation for the long-term evolution of debris in a TDE with realistic initial conditions, we show that at…
▽ More
Accretion of debris seems to be the natural mechanism to power the radiation emitted during a tidal disruption event (TDE), in which a supermassive black hole tears apart a star. However, this requires the prompt formation of a compact accretion disk. Here, using a fully relativistic global simulation for the long-term evolution of debris in a TDE with realistic initial conditions, we show that at most a tiny fraction of the bound mass enters such a disk on the timescale of observed flares. To "circularize" most of the bound mass entails an increase in the binding energy of that mass by a factor $\sim 30$; we find at most an order unity change. Our simulation suggests it would take a time scale comparable to a few tens of the characteristic mass fallback time to dissipate enough energy for "circularization". Instead, the bound debris forms an extended eccentric accretion flow with eccentricity $\simeq 0.4-0.5$ by $\sim 2$ fallback times. Although the energy dissipated in shocks in this large-scale flow is much smaller than the "circularization" energy, it matches the observed radiated energy very well. Nonetheless, the impact of shocks is not strong enough to unbind initially bound debris into an outflow.
△ Less
Submitted 22 August, 2023; v1 submitted 9 May, 2023;
originally announced May 2023.
-
Close encounters of black hole - star binaries with stellar-mass black holes
Authors:
Taeho Ryu,
Ruggero Valli,
Rudiger Pakmor,
Rosalba Perna,
Selma E. de Mink,
Volker Springel
Abstract:
Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics co…
▽ More
Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics code {\small AREPO}. This type of encounters can be divided into two classes depending on whether the final outcome includes BBHs. This outcome is primarily determined by which two objects meet at the first closest approach. BBHs are more likely to form when the star and the incoming BH encounter first with an impact parameter smaller than the binary's semimajor axis. In this case, the star is frequently disrupted. On the other hand, when the two BHs encounter first, frequent consequences are an orbit perturbation of the original binary or a binary member exchange. For the parameters chosen in this study, BBH formation, accompanied by stellar disruption, happens in roughly 1 out of 4 encounters. The close correlation between BBH formation and stellar disruption has possible implications for EM counterparts at the binary's merger. The BH that disrupts the star is promptly surrounded by an optically and geometrically thick disk with accretion rates exceeding the Eddington limit. If the debris disk cools fast enough to become long-lived, EM counterparts can be produced at the time of the BBH merger.
△ Less
Submitted 10 July, 2023; v1 submitted 4 April, 2023;
originally announced April 2023.
-
Quasi-periodic eruptions from impacts between the secondary and a rigidly precessing accretion disc in an extreme mass-ratio inspiral system
Authors:
Alessia Franchini,
Matteo Bonetti,
Alessandro Lupi,
Giovanni Miniutti,
Elisa Bortolas,
Margherita Giustini,
Massimo Dotti,
Alberto Sesana,
Riccardo Arcodia,
Taeho Ryu
Abstract:
X-ray quasi-periodic eruptions (QPEs) represent a recently discovered example of extreme X-ray variability associated with supermassive black holes. These are high-amplitude bursts recurring every few hours that are detected in the soft X-ray band from the nuclei of nearby galaxies whose optical spectra lack the broad emission lines typically observed in unobscured active galaxies. The physical or…
▽ More
X-ray quasi-periodic eruptions (QPEs) represent a recently discovered example of extreme X-ray variability associated with supermassive black holes. These are high-amplitude bursts recurring every few hours that are detected in the soft X-ray band from the nuclei of nearby galaxies whose optical spectra lack the broad emission lines typically observed in unobscured active galaxies. The physical origin of this new X-ray variability phenomenon is still unknown and several theoretical models have been presented. However, no attempt has been made so far to account for the varying QPE recurrence time and luminosity in individual sources, nor for the diversity of the QPE phenomenology in the different known erupters. We present a semi-analytical model based on an extreme mass-ratio inspiral (EMRI) system where the secondary intersects, along its orbit, a rigidly precessing accretion disc surrounding the primary. We assume that QPEs result from emission from an adiabatically expanding, initially optically thick gas cloud expelled from the disc plane at each impact. We produced synthetic X-ray light curves, which we then compared with X-ray data from four QPE sources: GSN 069, eRO-QPE1, eRO-QPE2, and RX J1301.9+2747. Our model aptly reproduces the diversity of QPE properties between the considered objects and it is also able to naturally account for the varying QPE amplitudes and recurrence times in individual sources. Future implementations will enable us to refine the match with the data and to estimate the system parameters precisely, making additional use of multi-epoch QPE data. We briefly discuss the nature of the secondary object, as well as the possible implications of our findings for the EMRI population at large.
△ Less
Submitted 14 June, 2023; v1 submitted 3 April, 2023;
originally announced April 2023.
-
Tidal Peeling Events: low-eccentricity tidal disruption of a star by a stellar-mass black hole
Authors:
Chengcheng Xin,
Zoltan Haiman,
Rosalba Perna,
Yihan Wang,
Taeho Ryu
Abstract:
Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei (AGNs). Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH (micro-tidal disruption event, or micro-TDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. He…
▽ More
Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei (AGNs). Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH (micro-tidal disruption event, or micro-TDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. Here we consider a scenario in which the star might approach the stellar-mass BH in a gradual, nearly circular inspiral, under the influence of dynamical friction on a circum-binary gas disk or three-body interactions in a star cluster. We perform hydro-dynamical simulations of this scenario using the smoothed particle hydrodynamics code PHANTOM. We find that the mass of the star is slowly stripped away by the BH. We call this gradual tidal disruption a "tidal-peeling event", or a TPE. Depending on the initial distance and eccentricity of the encounter, TPEs might exhibit significant accretion rates and orbital evolution distinct from those of a typical (eccentric) micro-TDE.
△ Less
Submitted 22 March, 2023;
originally announced March 2023.
-
Partial stellar tidal disruption events and their rates
Authors:
Elisa Bortolas,
Taeho Ryu,
Luca Broggi,
Alberto Sesana
Abstract:
Tidal disruption events (TDEs) of stars operated by massive black holes (MBHs) will be detected in thousands by upcoming facilities such as the Vera Rubin Observatory. In this work, we assess the rates of standard total TDEs, destroying the entire star, and partial TDEs, in which a stellar remnant survives the interaction, by solving 1-D Fokker-Planck equations. Our rate estimates are based on a n…
▽ More
Tidal disruption events (TDEs) of stars operated by massive black holes (MBHs) will be detected in thousands by upcoming facilities such as the Vera Rubin Observatory. In this work, we assess the rates of standard total TDEs, destroying the entire star, and partial TDEs, in which a stellar remnant survives the interaction, by solving 1-D Fokker-Planck equations. Our rate estimates are based on a novel definition of the loss cone whose size is commensurate to the largest radius at which partial disruptions can occur, as motivated by relativistic hydrodynamical simulations. Our novel approach unveils two important results. First, partial TDEs can be more abundant than total disruptions by a factor of a few to a few tens. Second, the rates of complete stellar disruptions can be overestimated by a factor of a few to a few tens if one neglects partial TDEs, as we find that many of the events classified at total disruptions in the standard framework are in fact partial TDEs. Accounting for partial TDEs is particularly relevant for galaxies harbouring a nuclear stellar cluster featuring many events coming from the empty loss cone. Based on these findings, we stress that partial disruptions should be considered when constraining the luminosity function of TDE flares; accounting for this may reconcile the theoretically estimated TDE rates with the observed ones.
△ Less
Submitted 4 July, 2023; v1 submitted 6 March, 2023;
originally announced March 2023.
-
Close Encounters of Tight Binary Stars with Stellar-mass Black Holes
Authors:
Taeho Ryu,
Rosalba Perna,
Ruediger Pakmor,
Jing-Ze Ma,
Rob Farmer,
Selma E. de Mink
Abstract:
Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of 3D hydrodynamics simulations using the moving-mesh code AREPO, we investigate the formation of transient phenomena and their properties in close encounters between an $2M_{\odot}$ or $20M_{\odot}$ equal-mass…
▽ More
Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of 3D hydrodynamics simulations using the moving-mesh code AREPO, we investigate the formation of transient phenomena and their properties in close encounters between an $2M_{\odot}$ or $20M_{\odot}$ equal-mass circular binary star and single $20M_{\odot}$ black hole (BH). Stars can be disrupted by the BH during dynamical interactions, naturally producing electromagnetic transient phenomena. Encounters with impact parameters smaller than the semimajor axis of the initial binary frequently lead to a variety of transients whose electromagnetic signatures are qualitatively different from those of ordinary disruption events involving just two bodies. These include the simultaneous or successive disruptions of both stars and one full disruption of one star accompanied by successive partial disruptions of the other star. On the other hand, when the impact parameter is larger than the semimajor axis of the initial binary, the binary is either simply tidally perturbed or dissociated into bound and unbound single stars ("micro-Hills" mechanism). The dissociation of $20M_{\odot}$ binaries can produce a runaway star and an active BH moving away from one another. Also, the binary dissociation can either produce an interacting binary with the BH, or a non-interacting, hard binary; both could be candidates of BH high- and low-mass X-ray binaries. Hence our simulations especially confirm that strong encounters can lead to the formation of the (generally difficult to form) BH low-mass X-ray binaries.
△ Less
Submitted 5 January, 2023; v1 submitted 4 November, 2022;
originally announced November 2022.
-
Extremely Relativistic Tidal Disruption Events
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran
Abstract:
Extreme tidal disruption events (eTDEs), which occur when a star passes very close to a supermassive black hole, may provide a way to observe a long-sought general relativistic effect: orbits that wind several times around a black hole and then leave. Through general relativistic hydrodynamics simulations, we show that such eTDEs are easily distinguished from most tidal disruptions, in which stars…
▽ More
Extreme tidal disruption events (eTDEs), which occur when a star passes very close to a supermassive black hole, may provide a way to observe a long-sought general relativistic effect: orbits that wind several times around a black hole and then leave. Through general relativistic hydrodynamics simulations, we show that such eTDEs are easily distinguished from most tidal disruptions, in which stars come close, but not so close, to the black hole. Following the stellar orbit, the debris in eTDEs is initially distributed in a crescent that quickly turns into tight spirals, from which some mass later falls back toward the black hole, while the remainder is ejected. Internal shocks within the infalling debris power the observed emission. The resulting light-curve rises rapidly to roughly the Eddington luminosity, maintains this level for between a few weeks and a year (depending on both the stellar mass and the black hole mass), and then drops. Most of its power is in thermal X-rays at a temperature $\sim (1-2)\times 10^{6}$ K ($\sim 100-200$ eV). The debris evolution and observational features of eTDEs are qualitatively different from ordinary TDEs, making eTDEs a new type of TDE. Although eTDEs are relatively rare for lower-mass black holes, most tidal disruptions around higher-mass black holes are extreme. Their detection offers a view of an exotic relativistic phenomenon previously inaccessible.
△ Less
Submitted 11 March, 2023; v1 submitted 31 October, 2022;
originally announced November 2022.
-
Close Encounters of Stars with Stellar-mass Black Hole Binaries
Authors:
Taeho Ryu,
Rosalba Perna,
Yihan Wang
Abstract:
Many astrophysical environments, from star clusters and globular clusters to the disks of Active Galactic Nuclei, are characterized by frequent interactions between stars and the compact objects that they leave behind. Here, using a suite of $3-D$ hydrodynamics simulations, we explore the outcome of close interactions between $1M_{\odot}$ stars and binary black holes (BBHs) in the gravitational wa…
▽ More
Many astrophysical environments, from star clusters and globular clusters to the disks of Active Galactic Nuclei, are characterized by frequent interactions between stars and the compact objects that they leave behind. Here, using a suite of $3-D$ hydrodynamics simulations, we explore the outcome of close interactions between $1M_{\odot}$ stars and binary black holes (BBHs) in the gravitational wave regime, resulting in a tidal disruption event (TDE) or a pure scattering, focusing on the accretion rates, the back reaction on the BH binary orbital parameters and the increase in the binary BH effective spin. We find that TDEs can make a significant impact on the binary orbit, which is often different from that of pure scattering. Binaries experiencing a prograde (retrograde) TDE tend to be widened (hardened) by up to $\simeq 20\%$. Initially circular binaries become more eccentric by $\lesssim 10\%$ by a prograde or retrograde TDE, whereas the eccentricity of initially eccentric binaries increases (decreases) by a retrograde (prograde) TDE by $\lesssim 5\%$. Overall a single TDE can generally result in changes of the gravitational wave-driven merger time scale by order unity. The accretion rates of both black holes are very highly super-Eddington, showing modulations (preferentially for retrograde TDEs) on a time scale of the orbital period, which can be a characteristic feature of BBH-driven TDEs. Prograde TDEs result in the effective spin parameter $χ$ to vary by $\lesssim 0.02$ while $χ\gtrsim -0.005$ for retrograde TDEs.
△ Less
Submitted 2 August, 2022; v1 submitted 1 June, 2022;
originally announced June 2022.
-
Tidal Disruption Events by Compact Supermassive Black Hole Binaries
Authors:
Taeho Ryu,
Alessandro A. Trani,
Nathan W. C. Leigh
Abstract:
Stars can be tidally destroyed or swallowed by supermassive black hole binaries. Using a large number of accurate few-body simulations, we investigate the enhancement and suppression of full and partial disruption and direct capture events by hard supermassive black hole binaries with a wide ranges of key parameters, i.e., the primary black hole mass ($10^{5}-10^{8}M_{\odot}$), the binary mass rat…
▽ More
Stars can be tidally destroyed or swallowed by supermassive black hole binaries. Using a large number of accurate few-body simulations, we investigate the enhancement and suppression of full and partial disruption and direct capture events by hard supermassive black hole binaries with a wide ranges of key parameters, i.e., the primary black hole mass ($10^{5}-10^{8}M_{\odot}$), the binary mass ratio ($10^{-3}-1$), the ratio of the binary semimajor axis to the hardening radius ($10^{-4}-1$), the binary eccentricity ($0.0-0.9$) and the stellar mass ($0.3-3M_{\odot}$). This is a significant extension of the parameter space compared to previous work. We show that the encounter probabilities of all three events are well-described by the encounter cross section, which is proportional to the pericenter distance. The probability of full disruptions by supermassive black hole binaries can be enhanced by up to a factor of $40-50$ or suppressed by up to a factor of $10$, relative to that by single black holes, depending on the binary parameters. Relativistic effects are not important for the primary black hole mass $\leq 10^{7}M_{\odot}$, but can provide an additional enhancement of the full disruption probability by less than a factor of $2-3$ for higher primary black hole masses. We provide a fitting formula for the full disruption probability by the hard supermassive black hole binaries that works for a wide range of parameters. We also find that partial disruption events can occur multiple times before full disruptions or direct captures, and their probabilities can be greater than that of full disruption events by a factor of three. Because partial disruption events can induce stellar spins and mass loss and change the orbits, it can significantly affect the overall full disruption event rate and the shape of the light curves.
△ Less
Submitted 2 August, 2022; v1 submitted 15 February, 2022;
originally announced February 2022.
-
Starfall: A heavy rain of stars in 'turning on' AGN
Authors:
B. McKernan,
K. E. S. Ford,
M. Cantiello,
M. J. Graham,
A. S. Jermyn,
N. W. C. Leigh,
T. Ryu,
D. Stern
Abstract:
As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and…
▽ More
As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and onto the SMBH. A large AGN loss cone ($θ_{\rm AGN,lc}$) can result from binary scatterings in the inner disk and yield tidal disruption events (TDEs). Signatures of starfall include optical/UV flares that rise in luminosity over time, particularly in the inner disk. If the SMBH mass is $M_{\rm SMBH} \ge 10^{8}M_{\odot}$, flares truncate abruptly and the star is swallowed. If $M_{\rm SMBH}<10^{8}M_{\odot}$, and if the infalling orbit lies within $θ_{\rm AGN,lc}$, the flare is followed by a TDE which can be prograde or retrograde relative to the AGN inner disk. Retrograde AGN TDEs are over-luminous and short-lived as in-plane ejecta collide with the inner disk and a lower AGN state follows. Prograde AGN TDEs add angular momentum to inner disk gas and so start off looking like regular TDEs but are followed by an AGN high state. Searches for such flare signatures test models of AGN `turn on', SMBH mass, as well as disk properties and the embedded population.
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
The Impact of Shocks on the Vertical Structure of Eccentric Disks
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran
Abstract:
Accretion disks whose matter follows eccentric orbits can arise in multiple astrophysical situations. Unlike circular orbit disks, the vertical gravity in eccentric disks varies around the orbit. In this paper, we investigate some of the dynamical effects of this varying gravity on the vertical structure using $1D$ hydrodynamics simulations of individual gas columns assumed to be mutually non-inte…
▽ More
Accretion disks whose matter follows eccentric orbits can arise in multiple astrophysical situations. Unlike circular orbit disks, the vertical gravity in eccentric disks varies around the orbit. In this paper, we investigate some of the dynamical effects of this varying gravity on the vertical structure using $1D$ hydrodynamics simulations of individual gas columns assumed to be mutually non-interacting. We find that time-dependent gravitational pumping generically creates shocks near pericenter; the energy dissipated in the shocks is taken from the orbital energy. Because the kinetic energy per unit mass in vertical motion near pericenter can be large compared to the net orbital energy, the shocked gas can be heated to nearly the virial temperature, and some of it becomes unbound. These shocks affect larger fractions of the disk mass for larger eccentricity and/or disk aspect ratio. If the orbit can be maintained despite orbital energy loss, diverse initial structures evolve in only a few orbits so that they follow a limit-cycle characterized by a low-entropy midplane and a much higher entropy outer layer. In favorable cases (such as the tidal disruption of stars by supermassive black holes), these effects could be a potentially important energy dissipation and mass loss mechanism.
△ Less
Submitted 27 July, 2021; v1 submitted 19 May, 2021;
originally announced May 2021.
-
Measuring stellar and black hole masses of tidal disruption events
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran
Abstract:
The flare produced when a star is tidally disrupted by a supermassive black hole holds potential as a diagnostic of both the black hole mass and the star mass. We propose a new method to realize this potential based upon a physical model of optical/UV light production in which shocks near the apocenters of debris orbits dissipate orbital energy, which is then radiated from that region. Measurement…
▽ More
The flare produced when a star is tidally disrupted by a supermassive black hole holds potential as a diagnostic of both the black hole mass and the star mass. We propose a new method to realize this potential based upon a physical model of optical/UV light production in which shocks near the apocenters of debris orbits dissipate orbital energy, which is then radiated from that region. Measurement of the optical/UV luminosity and color temperature at the peak of the flare leads directly to the two masses. The black hole mass depends mostly on the temperature observed at peak luminosity, while the mass of the disrupted star depends mostly on the peak luminosity. We introduce {\sc TDEmass}, a method to infer the black hole and stellar masses given these two input quantities. Using {\sc TDEmass}, we find, for 21 well-measured events, black hole masses between $5\times 10^5$ and $10^7 M_\odot$ and disrupted stars with initial masses between 0.6 and $13M_\odot$. An open-source {\sc python}-based tool for {\sc TDEmass} is available at https://github.com/taehoryu/TDEmass.git.
△ Less
Submitted 7 October, 2020; v1 submitted 27 July, 2020;
originally announced July 2020.
-
Tidal disruptions of main sequence stars -- IV. Relativistic effects and dependence on black hole mass
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran,
Scott C. Noble
Abstract:
Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass ($10^{5}\leq M_{\rm BH}/\mathrm{M}_{\odot}\leq 5\times 10^{7}$) and stellar mass ($0.3 \leq M_{\star} /\mathrm{M}_{\odot}\leq 3$) as larger $M_{\rm BH}$ leads to stronger relativist…
▽ More
Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass ($10^{5}\leq M_{\rm BH}/\mathrm{M}_{\odot}\leq 5\times 10^{7}$) and stellar mass ($0.3 \leq M_{\star} /\mathrm{M}_{\odot}\leq 3$) as larger $M_{\rm BH}$ leads to stronger relativistic effects. For fixed $M_{\star}$, as $M_{\rm BH}$ increases, the ratio of the maximum pericenter distance yielding full disruptions ($\mathcal{R}_{\rm t}$) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for $M_{\rm BH} = 5\times10^{7}~{\rm M}_\odot$, while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of two correction at $M_{\rm BH} = 5 \times 10^7~{\rm M}_\odot$. We find that for partial disruptions, the fractional remnant mass for a given ratio of the pericenter to $\mathcal{R}_{\rm t}$ is higher for larger $M_{\rm BH}$. These results have several implications. As $M_{\rm BH}$ increases above $\sim 10^7~{\rm M}_\odot$, the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross section ratio for partial to complete disruptions depends only weakly on $M_{\rm BH}$. The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. For $M_{\rm BH} \gtrsim 10^7~{\rm M}_\odot$, the $M_{\rm BH}$-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.
△ Less
Submitted 8 October, 2020; v1 submitted 8 January, 2020;
originally announced January 2020.
-
Tidal disruptions of main sequence stars -- III. Stellar mass dependence of the character of partial disruptions
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran,
Scott C. Noble
Abstract:
In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the stellar mass dependence of the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a resul…
▽ More
In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the stellar mass dependence of the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than main sequence stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris, but of either sign with respect to the black hole potential, while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. The remnants most deeply bound to the black hole go through a second tidal disruption event upon their first return to pericenter; if they have not thermally relaxed, they will be completely disrupted.
△ Less
Submitted 1 September, 2020; v1 submitted 8 January, 2020;
originally announced January 2020.
-
Tidal disruptions of main sequence stars -- II. Simulation methodology and stellar mass dependence of the character of full tidal disruptions
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran,
Scott C. Noble
Abstract:
This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. In the first paper (Paper I), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide detail…
▽ More
This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. In the first paper (Paper I), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions, focusing on the stellar mass dependence of the outcomes for a black hole of mass $10^{6}\rm{M}_{\odot}$. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars ($M_{\star} \lesssim 0.7~ {\rm M}_\odot$) is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor 2--2.5. The traditional estimate of the range of energies found in the debris is $\approx 1.4\times$ too large for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for many stellar vibration times because the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface until the star has reached a distance from the black hole $\sim O(10)r_{\rm t}$.
△ Less
Submitted 8 October, 2020; v1 submitted 8 January, 2020;
originally announced January 2020.
-
Tidal Disruptions of Main Sequence Stars -- I. Observable Quantities and their Dependence on Stellar and Black Hole Mass
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran,
Scott C. Noble
Abstract:
This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses ($0.15 \leq M_*/M_\odot \leq 10$) and six black hole masses (…
▽ More
This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses ($0.15 \leq M_*/M_\odot \leq 10$) and six black hole masses ($10^5 \leq M_{BH}/M_\odot \leq 5 \times 10^7$). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of $M_*$ for $M_*/M_\odot \lesssim 3$. Independent of $M_*$, relativistic effects increase the critical pericenter distance for full disruptions by up to a factor $\sim 3$ relative to the Newtonian prediction. The probability of a direct capture is also independent of $M_*$; at $M_{BH}/M_\odot \simeq 5 \times 10^6$ this probability is equal to that of a complete disruption. The width of the debris energy distribution $ΔE$ can differ from the standard estimate by factors from 0.35 to 2, depending on $M_*$ and $M_{BH}$, implying a corresponding change in the characteristic mass-return timescale. The "frozen-in approximation" is inconsistent with $ΔE$, and mass-loss continues over a long span of time. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly-universal relation between the star's angular momentum and the fraction of $M_*$ remaining. Within the "empty loss-cone" regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.
△ Less
Submitted 1 September, 2020; v1 submitted 8 January, 2020;
originally announced January 2020.
-
Tidal Disruptions of Main Sequence Stars -- V. The Varieties of Disruptions
Authors:
Julian Krolik,
Tsvi Piran,
Taeho Ryu
Abstract:
Tidal disruption events (TDEs), events in which a star passes very close to a supermassive black hole, are generally imagined as leading either to the star's complete disruption or to its passage directly into the black hole. In the former case it is widely believed that in all cases the bound portion of the debris quickly "circularizes" due to relativistic apsidal precession, i.e., forms a compac…
▽ More
Tidal disruption events (TDEs), events in which a star passes very close to a supermassive black hole, are generally imagined as leading either to the star's complete disruption or to its passage directly into the black hole. In the former case it is widely believed that in all cases the bound portion of the debris quickly "circularizes" due to relativistic apsidal precession, i.e., forms a compact accretion disk, and emits a flare of standardized lightcurve and spectrum. We show here that TDEs are more diverse and can be grouped into several distinct categories on the basis of stellar pericenter distance $r_p$; we calculate the relative frequency of these categories. In particular, because rapid circularization requires $r_p \lesssim 10r_g$ ($r_g \equiv GM_{\rm BH}/c^2$), it can happen in only a minority of total disruptions, $\lesssim 1/4$ when the black hole has mass $M_{\rm BH} = 10^6 M_\odot$. For larger pericenter distances, $10 < r_p/r_g < 27$ (for $M_{\rm BH}=10^6M_\odot$), main sequence stars are completely disrupted, but the bound debris orbits are highly eccentric and possess semimajor axes $\sim 100\times$ the scale of the expected compact disk. Partial disruptions with fractional mass-loss $\gtrsim 10\%$ should occur with a rate similar to that of total disruptions; for fractional mass-loss $\gtrsim 50\%$, the rate is $\approx 1/3$ as large. Partial disruptions -- which must precede total disruptions when the stars' angular momenta evolve in the "empty loss-cone" regime -- change the orbital energy by factors $\gtrsim O(1)$. Remnants of partial disruptions are in general far from thermal equilibrium. Depending on the orbital energy of the remnant and conditions within the stellar cluster surrounding the SMBH, it may return after hundreds or thousands of years and be fully disrupted, or it may rejoin the stellar cluster.
△ Less
Submitted 7 October, 2020; v1 submitted 9 January, 2020;
originally announced January 2020.
-
Tidal radii of main sequence stars -- III. Partial disruptions
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran,
Scott C. Noble
Abstract:
In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can d…
▽ More
In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than MS stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris (but of either sign with respect to the black hole potential) while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. Remnants bound to the black hole can possibly go through a second tidal disruption event.
△ Less
Submitted 13 January, 2020; v1 submitted 9 September, 2019;
originally announced September 2019.
-
Tidal radii of main sequence stars -- II. Simulation methodology and the character of full tidal disruptions
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran,
Scott Noble
Abstract:
This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. In the first paper (Ryu et al. 2019a), we gave an overview of this program and discussed the principa…
▽ More
This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. In the first paper (Ryu et al. 2019a), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor $2-2.5$. The traditional estimate of the range of energies found in the debris is approximately accurate for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for a long time because the ${\it instantaneous}$ tidal radius, the distance out to which the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface, stays comparable to the distance to the black hole until the star has reached $O(10)~r_{\rm t}$. These findings indicate significant failings in the popular "frozen-in" approximation.
△ Less
Submitted 13 January, 2020; v1 submitted 27 July, 2019;
originally announced July 2019.
-
Tidal radii of main sequence stars -- I. Physical tidal radius, semi-analytic model and their implications
Authors:
Taeho Ryu,
Julian Krolik,
Tsvi Piran,
Scott Noble
Abstract:
A star is tidally disrupted by a supermassive black hole when their separation is shorter than the "tidal radius". This quantity is often estimated on an order-of-magnitude basis without reference to the star's internal structure. Using MESA models for main sequence stars and fully general relativistic dynamics, we find the physical tidal radius for complete disruption $\cal{R}_t$ for a…
▽ More
A star is tidally disrupted by a supermassive black hole when their separation is shorter than the "tidal radius". This quantity is often estimated on an order-of-magnitude basis without reference to the star's internal structure. Using MESA models for main sequence stars and fully general relativistic dynamics, we find the physical tidal radius for complete disruption $\cal{R}_t$ for a $10^6M_\odot$ black hole (BH). We find that across a factor $\sim20$ in stellar mass $M_*$, i.e., $0.15M_{\odot}\leq M_*\leq3M_\odot$, $\cal{R}_t\sim27\times$(BH's gravitational radius). When comparing $\cal{R}_t$ with the commonly used order-of-magnitude estimate $r_t$, we find that $\cal{R}_t\sim1.05-1.45r_t$ for $0.15M_\odot\leq M_*\leq0.5M_\odot$, but between $0.5 M_\odot$ and $1 M_\odot$, $\cal{R}_t$ drops to $\sim 0.45r_t$, and it remains at this value up to $10 M_\odot$. The near-constancy of $\cal{R}_t$ implies a weaker dependence of the full disruption rate on $M_*$ than when predicted with $r_t$. The characteristic energy width of the debris $ΔE$ ranges from $\sim1.2Δ\cal{E}$ for low-mass stars to $\sim 0.35Δ\cal{E}$ for higher-mass stars, where $Δ\cal{E}=GM_{\rm BH}R_*/\cal{R}_t^{2}$. We present analytic fits for the $M_*$ dependence of $\cal{R}_t$ and $ΔE$; these fits lead to analytic expressions for the time of peak mass fallback rate and the maximal mass fallback rate. Our results also bear on the fraction of events leading to fast or slow circularization, as well as on the character of the tidal event occurring when the remnant of a partial disruption returns to the black hole. Using a semi-analytic model, we show that $\cal{R}_t$ is primarily determined by the star's central density rather than its mean density. For high-mass stars, the full disruption rate is roughly 1/4 the partial disruption rate, while this ratio is close to unity for low-mass stars.
△ Less
Submitted 13 January, 2020; v1 submitted 18 July, 2019;
originally announced July 2019.
-
Turbulence-driven thermal and kinetic energy fluxes in the atmospheres of hot Jupiters
Authors:
Taeho Ryu,
Michael Zingale,
Rosalba Perna
Abstract:
We have performed high resolution $3-$dimensional compressible hydrodynamics simulations to investigate the effects of shocks and turbulence on energy transport into hot Jupiter atmospheres, under a variety of shear gradients. We focus on a local atmospheric region to accurately follow the small-scale structures of turbulence and shocks. We find that the effects of turbulence above and below a she…
▽ More
We have performed high resolution $3-$dimensional compressible hydrodynamics simulations to investigate the effects of shocks and turbulence on energy transport into hot Jupiter atmospheres, under a variety of shear gradients. We focus on a local atmospheric region to accurately follow the small-scale structures of turbulence and shocks. We find that the effects of turbulence above and below a shear layer are different in scale and magnitude: below the shear layer, the effects of turbulence on the vertical energy transfer are local, generally $\lesssim2~\times(\text{scale height})$. However, turbulence can have a spatially and thermally-large influence on almost the entire region above the shear layer. We also find that shock formation is local and transient. Once the atmosphere becomes steady, the time-averaged heat energy flux at $P\sim 1$ bar is insignificant, on the order of 0.001\% of the incoming stellar flux with a shear motion at $P\sim 1$ mbar, and 0.1\% with a deeper shear layer at $P\sim 100$ mbar. Accordingly, the diffusion coefficient is higher for the deeper shear layer. Therefore, our results suggest that turbulence near less dense regions ($P\sim 1$ mbar) does not cause a sufficiently deep and large penetration of thermal energy to account for radius inflation in hot Jupiters, regardless of how violent the turbulence is. However, as the shear layer gets deeper, heat energy transfer becomes more effective throughout the atmosphere (upwards and downwards) due to a larger kinetic energy budget. Therefore, it is more important how deep turbulence occurs in the atmosphere, than how unstable the atmosphere is for effective energy transfer.
△ Less
Submitted 22 September, 2018; v1 submitted 20 June, 2018;
originally announced June 2018.
-
When do star clusters become multiple star systems? II. Toward a half-life formalism with four bodies
Authors:
Timur Ibragimov,
Nathan W. C. Leigh,
Taeho Ryu,
Teresa Panurach,
Rosalba Perna
Abstract:
We present a half-life formalism for describing the disruption of gravitationally-bound few-body systems, with a focus on binary-binary scattering. For negative total encounter energies, the four-body problem has three possible decay products in the point particle limit. For each decay product and a given set of initial conditions, we obtain directly from numerical scattering simulations the half-…
▽ More
We present a half-life formalism for describing the disruption of gravitationally-bound few-body systems, with a focus on binary-binary scattering. For negative total encounter energies, the four-body problem has three possible decay products in the point particle limit. For each decay product and a given set of initial conditions, we obtain directly from numerical scattering simulations the half-life for the distribution of disruption times. As in radioactive decay, the half-lives should provide a direct prediction for the relative fractions of each decay product. We test this prediction with simulated data and find good agreement with our hypothesis. We briefly discuss applications of this feature of the gravitational four-body problem to populations of black holes in globular clusters. This paper, the second in the series, builds on extending the remarkable similarity between gravitational chaos at the macroscopic scale and radioactive decay at the microscopic scale to larger-N systems.
△ Less
Submitted 14 March, 2018;
originally announced March 2018.
-
Population III X-Ray Binaries
Authors:
Taeho Ryu,
Takamitsu L. Tanaka,
Rosalba Perna
Abstract:
Understanding of the role of X-rays for driving the thermal evolution of the intergalactic medium (IGM) at high redshifts is one of important questions in astrophysics. High-mass X-ray binaries (HMXBs) in early stellar populations are prime X-ray source; however, their formation efficiency is not well understood. Using $N$-body simulations, we estimate the HMXB formation rate via mutual gravitatio…
▽ More
Understanding of the role of X-rays for driving the thermal evolution of the intergalactic medium (IGM) at high redshifts is one of important questions in astrophysics. High-mass X-ray binaries (HMXBs) in early stellar populations are prime X-ray source; however, their formation efficiency is not well understood. Using $N$-body simulations, we estimate the HMXB formation rate via mutual gravitational interactions of nascent, small groups of the Population~III stars. We find that HMXBs form at a rate of one per $\gtrsim 10^{4}M_{\odot}$ in newly born stars, and that they emit with a power of $\sim 10^{41} {\rm erg}~{\rm s}^{-1}$ in the $2-10$ keV band per star formation rate (SFR). This value is a factor $\sim 10^{2}$ larger than what is observed in star forming galaxies at lower redshifts; the X-ray production from early HMXBs would have been even more copious, if they also formed \textit{in situ} or via migration in protostellar disks. Combining our results with earlier studies suggests that early HMXBs were highly effective at heating the IGM and leaving a strong 21 cm signature. We discuss broader implications of our results, such as the rate of long gamma-ray bursts from Population~III stars and the direct collapse channel for massive black hole formation.
△ Less
Submitted 24 October, 2017;
originally announced October 2017.
-
Exoplanets around Low-mass Stars Unveiled by K2
Authors:
Teruyuki Hirano,
Fei Dai,
Davide Gandolfi,
Akihiko Fukui,
John H. Livingston,
Kohei Miyakawa,
Michael Endl,
William D. Cochran,
Francisco J. Alonso-Floriano,
Masayuki Kuzuhara,
David Montes,
Tsuguru Ryu,
Simon Albrecht,
Oscar Barragan,
Juan Cabrera,
Szilard Csizmadia,
Hans Deeg,
Philipp Eigmüller,
Anders Erikson,
Malcolm Fridlund,
Sascha Grziwa,
Eike W. Guenther,
Artie P. Hatzes,
Judith Korth,
Tomoyuki Kudo
, et al. (16 additional authors not shown)
Abstract:
We present the detection and follow-up observations of planetary candidates around low-mass stars observed by the K2 mission. Based on light-curve analysis, adaptive-optics imaging, and optical spectroscopy at low and high resolution (including radial velocity measurements), we validate 16 planets around 12 low-mass stars observed during K2 campaigns 5-10. Among the 16 planets, 12 are newly valida…
▽ More
We present the detection and follow-up observations of planetary candidates around low-mass stars observed by the K2 mission. Based on light-curve analysis, adaptive-optics imaging, and optical spectroscopy at low and high resolution (including radial velocity measurements), we validate 16 planets around 12 low-mass stars observed during K2 campaigns 5-10. Among the 16 planets, 12 are newly validated, with orbital periods ranging from 0.96-33 days. For one of the planets (K2-151b) we present ground-based transit photometry, allowing us to refine the ephemerides. Combining our K2 M-dwarf planets together with the validated or confirmed planets found previously, we investigate the dependence of planet radius $R_p$ on stellar insolation and metallicity [Fe/H]. We confirm that for periods $P\lesssim 2$ days, planets with a radius $R_p\gtrsim 2\,R_\oplus$ are less common than planets with a radius between 1-2$\,R_\oplus$. We also see a hint of the "radius valley" between 1.5 and 2$\,R_\oplus$ that has been seen for close-in planets around FGK stars. These features in the radius/period distribution could be attributed to photoevaporation of planetary envelopes by high-energy photons from the host star, as they have for FGK stars. For the M dwarfs, though, the features are not as well defined, and we cannot rule out other explanations such as atmospheric loss from internal planetary heat sources, or truncation of the protoplanetary disk. There also appears to be a relation between planet size and metallicity: those few planets larger than about 3 $R_\oplus$ are found around the most metal-rich M dwarfs.
△ Less
Submitted 18 January, 2018; v1 submitted 9 October, 2017;
originally announced October 2017.
-
Interactions between multiple supermassive black holes in galactic nuclei: a solution to the final parsec problem
Authors:
Taeho Ryu,
Rosalba Perna,
Zoltán Haiman,
Jeremiah P. Ostriker,
Nicholas C. Stone
Abstract:
Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evoluti…
▽ More
Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evolution, and allow us to separately explore the role of dynamical friction and that of multi-body BH interactions on BH mergers. Using the computed merger rates, we infer the stochastic gravitational wave background (GWB). Our dynamical approach is a first attempt to study the dynamical evolution of multiple SMBHs in the host galaxies undergoing mergers with various mass ratios ($10^{-4} < q_{\star} < 1$). Our main result demonstrates that SMBH binaries are able to merge in both scenarios. In the empty loss cone case, we find that BHs merge via multi-body interactions, avoiding the `final parsec' problem, and entering the PTA band with substantial orbital eccentricity. Our full loss cone treatment, albeit more approximate, suggests that the eccentricity becomes even higher when GWs become dominant, leading to rapid coalescences (binary lifetime $\lesssim1 {\rm ~Gyr}$). Despite the lower merger rates in the empty loss cone case, due to their higher mass ratios and lower redshifts, the GWB in the full/empty loss cone models are comparable ($0.70\times10^{-15}$ and $0.53\times10^{-15}$ at a frequency of $1~{\rm yr}^{-1}$, respectively). Finally, we compute the effects of high eccentricities on the GWB spectrum.
△ Less
Submitted 25 September, 2017; v1 submitted 19 September, 2017;
originally announced September 2017.
-
Formation of runaway stars in a star-cluster potential
Authors:
Taeho Ryu,
Nathan W. C. Leigh,
Rosalba Perna
Abstract:
We study the formation of runaway stars due to binary-binary (2+2) interactions in young star-forming clusters and/or associations. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We focus on interactions that produce two single stars and a binary, and study the outcomes as a…
▽ More
We study the formation of runaway stars due to binary-binary (2+2) interactions in young star-forming clusters and/or associations. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We focus on interactions that produce two single stars and a binary, and study the outcomes as a function of the depth of the background potential, within a range typical of cluster cores. As reference parameters for the observational properties, we use those observed for the system of runaway stars AE Aur and $μ$ Col and binary $ι$ Ori. We find that the outcome fractions have no appreciable dependence on the depth of the potential, and neither do the velocities of the ejected single stars. However, as the potential gets deeper and a larger fraction of binaries remain trapped, two binary populations emerge, with the escaped component having higher speeds and shorter semi-major axes than the trapped one. Additionally, we find that the relative angles between the ejected products are generally large. In particular, the angle between the ejected fastest star and the escaped binary is typically $\gtrsim 120-135^{\circ}$, with a peak at around $160^{\circ}$. However, as the potential gets deeper, the angle distribution becomes broader. Finally, we discuss the implications of our results for the interpretation of the properties of the runaway stars AE Aur and $μ$ Col.
△ Less
Submitted 5 June, 2017; v1 submitted 24 March, 2017;
originally announced March 2017.
-
An analytic method for identifying dynamically-formed runaway stars
Authors:
Taeho Ryu,
Nathan W. C. Leigh,
Rosalba Perna
Abstract:
In this paper, we study the 3-body products (two single stars and a binary) of binary-binary (2+2) scattering interactions. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We derive analytically a simple formula relating the angle between the velocity vectors of the two eject…
▽ More
In this paper, we study the 3-body products (two single stars and a binary) of binary-binary (2+2) scattering interactions. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We derive analytically a simple formula relating the angle between the velocity vectors of the two ejected single stars and the orbital separation of the remaining binary. We compare our analytic formulation to numerical scattering simulations, and illustrate that the agreement is excellent, both in isolation and in a homogeneous background potential. Our results are ideally suited for application to the GAIA database, which is expected to identify many hundred runaway stars. The analytic relation presented here has the potential to identify runaway stars formed dynamically with high confidence. Finally, by applying our method to the runaways AE Aur and $μ$ Col, we illustrate that it can be used to constrain the history of the background potential, which was denser than the presently observed density in the case of the Trapezium cluster.
△ Less
Submitted 9 May, 2017; v1 submitted 24 March, 2017;
originally announced March 2017.
-
The remarkable outburst of the highly evolved post-period-minimum dwarf nova SSS J122221.7-311525
Authors:
V. V. Neustroev,
T. R. Marsh,
S. V. Zharikov,
C. Knigge,
E. Kuulkers,
J. P. Osborne,
K. L. Page,
D. Steeghs,
V. F. Suleimanov,
G. Tovmassian,
E. Breedt,
A. Frebel,
Ma. T. Garcia-Diaz,
F. -J. Hambsch,
H. Jacobson,
S. G. Parsons,
T. Ryu,
L. Sabin,
G. Sjoberg,
A. S. Miroshnichenko,
D. E. Reichart,
J. B. Haislip,
K. M. Ivarsen,
A. P. LaCluyze,
J. P. Moore
Abstract:
We report extensive 3-yr multiwavelength observations of the WZ Sge-type dwarf nova SSS J122221.7-311525 during its unusual double superoutburst, the following decline and in quiescence. The second segment of the superoutburst had a long duration of 33 d and a very gentle decline with a rate of 0.02 mag/d, and it displayed an extended post-outburst decline lasting at least 500 d. Simultaneously wi…
▽ More
We report extensive 3-yr multiwavelength observations of the WZ Sge-type dwarf nova SSS J122221.7-311525 during its unusual double superoutburst, the following decline and in quiescence. The second segment of the superoutburst had a long duration of 33 d and a very gentle decline with a rate of 0.02 mag/d, and it displayed an extended post-outburst decline lasting at least 500 d. Simultaneously with the start of the rapid fading from the superoutburst plateau, the system showed the appearance of a strong near-infrared excess resulting in very red colours, which reached extreme values (B-I~1.4) about 20 d later. The colours then became bluer again, but it took at least 250 d to acquire a stable level. Superhumps were clearly visible in the light curve from our very first time-resolved observations until at least 420 d after the rapid fading from the superoutburst. The spectroscopic and photometric data revealed an orbital period of 109.80 min and a fractional superhump period excess <0.8 per cent, indicating a very low mass ratio q<0.045. With such a small mass ratio the donor mass should be below the hydrogen-burning minimum mass limit. The observed infrared flux in quiescence is indeed much lower than is expected from a cataclysmic variable with a near-main-sequence donor star. This strongly suggests a brown-dwarf-like nature for the donor and that SSS J122221.7-311525 has already evolved away from the period minimum towards longer periods, with the donor now extremely dim.
△ Less
Submitted 22 February, 2017; v1 submitted 11 January, 2017;
originally announced January 2017.
-
Multi-Color Simultaneous Photometry of the T-Tauri Star Having A Planetary Candidate CVSO 30
Authors:
Masahiro Onitsuka,
Akihiko Fukui,
Norio Narita,
Teruyuki Hirano,
Nobuhiko Kusakabe,
Tsuguru Ryu,
Motohide Tamura
Abstract:
We present three-band simultaneous observations of a weak-line T-Tauri star CVSO~30 (PTFO~8-8695), which is one of the youngest objects having a candidate transiting planet. The data were obtained with the Multicolor Simultaneous Camera for studying Atmospheres of Transiting exoplanets (MuSCAT) on the 188 cm telescope at Okayama Astrophysical Observatory in Japan. We observed the fading event in t…
▽ More
We present three-band simultaneous observations of a weak-line T-Tauri star CVSO~30 (PTFO~8-8695), which is one of the youngest objects having a candidate transiting planet. The data were obtained with the Multicolor Simultaneous Camera for studying Atmospheres of Transiting exoplanets (MuSCAT) on the 188 cm telescope at Okayama Astrophysical Observatory in Japan. We observed the fading event in the $g^{\prime}_2$-, $r^{\prime}_2$-, and $z_{\rm s,2}$-bands simultaneously. As a result, we find a significant wavelength dependence of fading depths of about 3.1\%, 1.7\%, 1.0\% for the $g^{\prime}_2$-, $r^{\prime}_2$-, and $z_{\rm s,2}$-bands, respectively. A cloudless H/He dominant atmosphere of a hot Jupiter cannot explain this large wavelength dependence. Additionally, we rule out a scenario by the occultation of the gravity-darkened host star. Thus our result is in favor of the fading origin as circumstellar dust clump or occultation of an accretion hotspot.
△ Less
Submitted 6 January, 2017;
originally announced January 2017.
-
The K2-ESPRINT Project VI: K2-105 b, a Hot-Neptune around a Metal-rich G-dwarf
Authors:
Norio Narita,
Teruyuki Hirano,
Akihiko Fukui,
Yasunori Hori,
Fei Dai,
Liang Yu,
John Livingston,
Tsuguru Ryu,
Grzegorz Nowak,
Masayuki Kuzuhara,
Yoichi Takeda,
Simon Albrecht,
Tomoyuki Kudo,
Nobuhiko Kusakabe,
Enric Palle,
Ignasi Ribas,
Motohide Tamura,
Vincent Van Eylen,
Joshua N. Winn
Abstract:
We report on the confirmation that the candidate transits observed for the star EPIC 211525389 are due to a short-period Neptune-sized planet. The host star, located in K2 campaign field 5, is a metal-rich ([Fe/H] = 0.26$\pm$0.05) G-dwarf (T_eff = 5430$\pm$70 K and log g = 4.48$\pm$0.09), based on observations with the High Dispersion Spectrograph (HDS) on the Subaru 8.2m telescope. High-spatial r…
▽ More
We report on the confirmation that the candidate transits observed for the star EPIC 211525389 are due to a short-period Neptune-sized planet. The host star, located in K2 campaign field 5, is a metal-rich ([Fe/H] = 0.26$\pm$0.05) G-dwarf (T_eff = 5430$\pm$70 K and log g = 4.48$\pm$0.09), based on observations with the High Dispersion Spectrograph (HDS) on the Subaru 8.2m telescope. High-spatial resolution AO imaging with HiCIAO on the Subaru telescope excludes faint companions near the host star, and the false positive probability of this target is found to be <$10^{-6}$ using the open source vespa code. A joint analysis of transit light curves from K2 and additional ground-based multi-color transit photometry with MuSCAT on the Okayama 1.88m telescope gives the orbital period of P = 8.266902$\pm$0.000070 days and consistent transit depths of $R_p/R_\star \sim 0.035$ or $(R_p/R_\star)^2 \sim 0.0012$. The transit depth corresponds to a planetary radius of $R_p = 3.59_{-0.39}^{+0.44} R_{\oplus}$, indicating that EPIC 211525389 b is a short-period Neptune-sized planet. Radial velocities of the host star, obtained with the Subaru HDS, lead to a 3σ upper limit of 90 $M_{\oplus} (0.00027 M_{\odot})$ on the mass of EPIC 211525389 b, confirming its planetary nature. We expect this planet, newly named K2-105 b, to be the subject of future studies to characterize its mass, atmosphere, spin-orbit (mis)alignment, as well as investigate the possibility of additional planets in the system.
△ Less
Submitted 5 January, 2017;
originally announced January 2017.
-
On the likelihood of detecting gravitational waves from Population III compact object binaries
Authors:
K. Belczynski,
T. Ryu,
R. Perna,
E. Berti,
T. L. Tanaka,
T. Bulik
Abstract:
We study the contribution of binary black hole (BH-BH) mergers from the first, metal-free stars in the Universe (Pop III) to gravitational wave detection rates. Our study combines initial conditions for the formation of Pop III stars based on N-body simulations of binary formation (including rates, binary fraction, initial mass function, orbital separation and eccentricity distributions) with an u…
▽ More
We study the contribution of binary black hole (BH-BH) mergers from the first, metal-free stars in the Universe (Pop III) to gravitational wave detection rates. Our study combines initial conditions for the formation of Pop III stars based on N-body simulations of binary formation (including rates, binary fraction, initial mass function, orbital separation and eccentricity distributions) with an updated model of stellar evolution specific for Pop III stars. We find that the merger rate of these Pop III BH-BH systems is relatively small (< 0.1 Gpc^-3 yr^-1) at low redshifts (z<2), where it can be compared with the LIGO empirical estimate of 9-240 Gpc^-3 yr^-1 (Abbott et al. 2016). The predicted rates are even smaller for Pop III double neutron star and black hole neutron star mergers. Our rates are compatible with those of Hartwig et al. (2016), but significantly smaller than those found in previous work (Bond & Carr 1984; Belczynski et al. 2004; Kinugawa et al. 2014, 2016). We explain the reasons for this discrepancy by means of detailed model comparisons and point out that (i) identification of Pop III BH-BH mergers may not be possible by advanced LIGO, and (ii) the level of stochastic gravitational wave background from Pop III mergers may be lower than recently estimated (Kowalska et al. 2012; Inayoshi et al. 2016; Dvorkin et al. 2016). We further estimate gravitational wave detection rates for third-generation interferometric detectors. Our calculations are relevant for low to moderately rotating Pop III stars. We can now exclude significant (> 1 per cent) contribution of these stars to low-redshift BH-BH mergers. However, it remains to be tested whether (and at what level) rapidly spinning Pop III stars (homogeneous evolution) can contribute to BH-BH mergers in the local Universe.
△ Less
Submitted 5 December, 2016;
originally announced December 2016.
-
Numerical study of the chaotic N=4 problem in a background potential
Authors:
Taeho Ryu,
Nathan W. C. Leigh,
Rosalba Perna
Abstract:
We perform a large suite of $N=4$ numerical scattering experiments between two identical binaries consisting of identical point particles in a (continuous) background potential. For investigative purposes, albeit without loss of generality, we assume that the potential corresponds to a uniform (natal or star-forming) gas medium. We explore a range of constant gas densities, from…
▽ More
We perform a large suite of $N=4$ numerical scattering experiments between two identical binaries consisting of identical point particles in a (continuous) background potential. For investigative purposes, albeit without loss of generality, we assume that the potential corresponds to a uniform (natal or star-forming) gas medium. We explore a range of constant gas densities, from $n=10~ {\rm cm}^{-3}$ to $10^{5}~ {\rm cm}^{-3}$. These densities are relevant for various astrophysical environments, including molecular clouds (i.e., star-forming regions) and denser, fragmented cores within these clouds. Our primary goal is to characterize the effects of the background potential on the subsequent stellar dynamics. We consider the outcome probabilities as well as the properties of any binaries formed during the binary-binary encounters, such as the distributions of binary binding energies and eccentricities. We also present the final velocity distributions of the ejected single stars. The background potential has two important effects on the stellar dynamics: 1) The potential acts to reset the zero-point of the total system energy, which in turn affects the types and properties of the products of the encounter; 2) For higher $n$ and weakly bound systems (i.e., large semimajor axes), the stellar dynamics are significantly affected when stars become trapped in the potential, oscillating around the system centre of mass (CM). This, in turn, increases the number of scattering events between stars (single, binary or triple) near the CM and makes it harder for single stars to escape to infinity. This ultimately leads to the preferential ionization of triples and wide binaries and the survival of compact binaries, with the single stars escaping at very high ejection velocities.
△ Less
Submitted 13 February, 2017; v1 submitted 28 November, 2016;
originally announced November 2016.
-
Ground-based Transit Observation of the Habitable-zone Super-Earth K2-3d
Authors:
Akihiko Fukui,
John Livingston,
Norio Narita,
Teruyuki Hirano,
Masahiro Onitsuka,
Tsuguru Ryu,
Nobuhiko Kusakabe
Abstract:
We report the first ground-based transit observation of K2-3d, a 1.5 R_Earth planet supposedly within the habitable zone around a bright M-dwarf host star, using the Okayama 188 cm telescope and the multi(grz)-band imager MuSCAT. Although the depth of the transit (0.7 mmag) is smaller than the photometric precisions (1.2, 0.9, and 1.2 mmag per 60 s for the g, r, and z bands, respectively), we marg…
▽ More
We report the first ground-based transit observation of K2-3d, a 1.5 R_Earth planet supposedly within the habitable zone around a bright M-dwarf host star, using the Okayama 188 cm telescope and the multi(grz)-band imager MuSCAT. Although the depth of the transit (0.7 mmag) is smaller than the photometric precisions (1.2, 0.9, and 1.2 mmag per 60 s for the g, r, and z bands, respectively), we marginally but consistently identify the transit signal in all three bands, by taking advantage of the transit parameters from K2, and by introducing a novel technique that leverages multi-band information to reduce the systematics caused by second-order extinction. We also revisit previously analyzed Spitzer transit observations of K2-3d to investigate the possibility of systematic offsets in transit timing, and find that all the timing data can be explained well by a linear ephemeris. We revise the orbital period of K2-3d to be 44.55612 \pm 0.00021 days, which corrects the predicted transit times for 2019, i.e., the era of the James Webb Space Telescope, by \sim80 minutes. Our observation demonstrates that (1) even ground-based, 2 m class telescopes can play an important role in refining the transit ephemeris of small-sized, long-period planets, and that (2) a multi-band imager is useful to reduce the systematics of atmospheric origin, in particular for bluer bands and for observations conducted at low-altitude observatories.
△ Less
Submitted 25 November, 2016; v1 submitted 5 October, 2016;
originally announced October 2016.
-
SEEDS direct imaging of the RV-detected companion to V450 Andromedae, and characterization of the system
Authors:
K. G. Hełminiak,
M. Kuzuhara,
K. Mede,
T. D. Brandt,
R. Kandori,
T. Suenaga,
N. Kusakabe,
N. Narita,
J. C. Carson,
T. Currie,
T. Kudo,
J. Hashimoto,
L. Abe,
E. Akiyama,
W. Brandner,
M. Feldt,
M. Goto,
C. A. Grady,
O. Guyon,
Y. Hayano,
M. Hayashi,
S. S. Hayashi,
T. Henning,
K. W. Hodapp,
M. Ishii
, et al. (27 additional authors not shown)
Abstract:
We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resoluti…
▽ More
We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resolution spectra and radial velocity (RV) measurements, along with RVs from the Lick planet search program. We combined our multi-epoch astrometry with these archival, partially unpublished RVs, and found that the companion is a low-mass star, not a brown dwarf, as previously suggested. We found the best-fitting dynamical masses to be $m_1=1.141_{-0.091}^{+0.037}$ and $m_2=0.279^{+0.023}_{-0.020}$ M$_\odot$. We also performed spectral analysis of the SOPHIE spectra with the iSpec code. The Hipparcos time-series photometry shows a periodicity of $P=5.743$ d, which is also seen in SOPHIE spectra as an RV modulation of the star A. We interpret it as being caused by spots on the stellar surface, and the star to be rotating with the given period. From the rotation and level of activity, we found that the system is $380^{+220}_{-100}$ Myr old, consistent with an isochrone analysis ($220^{+2120}_{-90}$ Myr). This work may serve as a test case for future studies of low-mass stars, brown dwarfs and exoplanets by combination of RV and direct imaging data.
△ Less
Submitted 30 August, 2016;
originally announced August 2016.
-
V5852 Sgr: An Unusual Nova Possibly Associated with the Sagittarius Stream
Authors:
E. Aydi,
P. Mróz,
P. A. Whitelock,
S. Mohamed,
Ł. Wyrzykowski,
A. Udalski,
P. Vaisanen,
T. Nagayama,
M. Dominik,
A. Scholz,
H. Onozato,
R. E. Williams,
S. T. Hodgkin,
S. Nishiyama,
M. Yamagishi,
A. M. S. Smith,
T. Ryu,
A. Iwamatsu,
I. Kawamata
Abstract:
We report spectroscopic and photometric follow-up of the peculiar nova V5852~Sgr (discovered as OGLE-2015-NOVA-01), which exhibits a combination of features from different nova classes. The photometry shows a flat-topped light curve with quasi-periodic oscillations, then a smooth decline followed by two fainter recoveries in brightness. Spectroscopy with the Southern African Large Telescope shows…
▽ More
We report spectroscopic and photometric follow-up of the peculiar nova V5852~Sgr (discovered as OGLE-2015-NOVA-01), which exhibits a combination of features from different nova classes. The photometry shows a flat-topped light curve with quasi-periodic oscillations, then a smooth decline followed by two fainter recoveries in brightness. Spectroscopy with the Southern African Large Telescope shows first a classical nova with an Fe II or Fe IIb spectral type. In the later spectrum, broad emissions from helium, nitrogen and oxygen are prominent and the iron has faded which could be an indication to the start of the nebular phase. The line widths suggest ejection velocities around $1000\,{\rm km\,s^{-1}}$. The nova is in the direction of the Galactic bulge and is heavily reddened by an uncertain amount. The $V$ magnitude 16 days after maximum enables a distance to be estimated and this suggests that the nova may be in the extreme trailing stream of the Sagittarius dwarf spheroidal galaxy. If so it is the first nova to be detected from that, or from any dwarf spheroidal galaxy. Given the uncertainty of the method and the unusual light curve we cannot rule out the possibility that it is in the bulge or even the Galactic disk behind the bulge.
△ Less
Submitted 8 June, 2016;
originally announced June 2016.
-
The K2-ESPRINT Project V: a short-period giant planet orbiting a subgiant star
Authors:
Vincent Van Eylen,
Simon Albrecht,
Davide Gandolfi,
Fei Dai,
Joshua N. Winn,
Teriyuki Hirano,
Norio Narita,
Hans Bruntt,
Jorge Prieto-Arranz,
Victor J. S. Bejar,
Grzegorz Nowak,
Mikkel N. Lund,
Enric Palle,
Ignasi Ribas,
Roberto Sanchis-Ojeda,
Liang Yu,
Pamela Arriagada,
R. Paul Butler,
Jeffrey D. Crane,
Rasmus Handberg,
Hans Deeg,
Jens Jessen-Hansen,
John A. Johnson,
David Nespral,
Leslie Rogers
, et al. (9 additional authors not shown)
Abstract:
We report on the discovery and characterization of the transiting planet K2-39b (EPIC 206247743b). With an orbital period of 4.6 days, it is the shortest-period planet orbiting a subgiant star known to date. Such planets are rare, with only a handful of known cases. The reason for this is poorly understood, but may reflect differences in planet occurrence around the relatively high-mass stars that…
▽ More
We report on the discovery and characterization of the transiting planet K2-39b (EPIC 206247743b). With an orbital period of 4.6 days, it is the shortest-period planet orbiting a subgiant star known to date. Such planets are rare, with only a handful of known cases. The reason for this is poorly understood, but may reflect differences in planet occurrence around the relatively high-mass stars that have been surveyed, or may be the result of tidal destruction of such planets. K2-39 is an evolved star with a spectroscopically derived stellar radius and mass of $3.88^{+0.48}_{-0.42}~\mathrm{R_\odot}$ and $1.53^{+0.13}_{-0.12}~\mathrm{M_\odot}$, respectively, and a very close-in transiting planet, with $a/R_\star = 3.4$. Radial velocity (RV) follow-up using the HARPS, FIES and PFS instruments leads to a planetary mass of $50.3^{+9.7}_{-9.4}~\mathrm{M_\oplus}$. In combination with a radius measurement of $8.3 \pm 1.1~\mathrm{R_\oplus}$, this results in a mean planetary density of $0.50^{+0.29}_{-0.17}$ g~cm$^{-3}$. We furthermore discover a long-term RV trend, which may be caused by a long-period planet or stellar companion. Because K2-39b has a short orbital period, its existence makes it seem unlikely that tidal destruction is wholly responsible for the differences in planet populations around subgiant and main-sequence stars. Future monitoring of the transits of this system may enable the detection of period decay and constrain the tidal dissipation rates of subgiant stars.
△ Less
Submitted 30 May, 2016;
originally announced May 2016.