Showing 1–50 of 165 results for author: Perets, H B

New Rotation Period Measurements for Kepler Stars Using Deep Learning: The 100K Sample

Authors: Ilay Kamai, Hagai B. Perets

Abstract: We propose a new framework to predict stellar properties from light curves. We analyze the light-curve data from the Kepler space mission and develop a novel tool for deriving the stellar rotation periods for main-sequence stars. Using this tool, we provide the largest (108785 stars) and most accurate (an average error of $1.6$ Days) sample of stellar rotations to date. Our model, LightPred, is a… ▽ More We propose a new framework to predict stellar properties from light curves. We analyze the light-curve data from the Kepler space mission and develop a novel tool for deriving the stellar rotation periods for main-sequence stars. Using this tool, we provide the largest (108785 stars) and most accurate (an average error of $1.6$ Days) sample of stellar rotations to date. Our model, LightPred, is a novel deep-learning model designed to extract stellar rotation periods from light curves. The model utilizes a dual-branch architecture combining Long Short-Term Memory (LSTM) and Transformer components to capture both temporal and global features within the data. We train LightPred on a dataset of simulated light curves generated using a realistic spot model and enhance its performance through self-supervised contrastive pre-training on Kepler light curves. Our evaluation demonstrates that LightPred outperforms classical methods like the Autocorrelation Function (ACF) in terms of accuracy and robustness. We apply LightPred to the Kepler dataset, generating the largest catalog to date of stellar rotation periods for main-sequence stars. Our analysis reveals a systematic shift towards shorter periods compared to previous studies, suggesting a potential revision of stellar age estimates. We also investigate the impact of stellar activity on period determination and find evidence for a distinct period-activity relation. Additionally, we confirm tidal synchronization in eclipsing binaries with orbital periods shorter than 10 days. Our findings highlight the potential of deep learning in extracting fundamental stellar properties from light curves, opening new avenues for understanding stellar evolution and population demographics. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Adaptive Habitability of Exoplanets: Thriving Under Extreme Environmental Change

Authors: Itay Weintraub, Hagai B. Perets

Abstract: The dynamic nature of life's ability to thrive in diverse and changing planetary environments suggests that habitability and survival depend on the evolutionary path and life adaptation to environmental conditions. Here we explore such "adaptive habitability" through astro-ecological models. We study the interplay between temperature adaptation and environmental fluctuations, particularly those in… ▽ More The dynamic nature of life's ability to thrive in diverse and changing planetary environments suggests that habitability and survival depend on the evolutionary path and life adaptation to environmental conditions. Here we explore such "adaptive habitability" through astro-ecological models. We study the interplay between temperature adaptation and environmental fluctuations, particularly those induced by solar activity and orbital dynamics. We present a simplified ecological-evolutionary model to investigate the limits of life's adaptability on a planetary scale. By incorporating complexities such as multiple niches, migration, species interactions, and realistic temperature variations, we demonstrate the potential for adaptive habitability in the face of both gradual and abrupt environmental changes. Through simulations encompassing monotonic, periodic, and secular dynamical evolution-induced temperature profiles, we identify critical thresholds for survival and extinction, highlighting the importance of phenotypic variance and dispersal rates in adapting to varying environmental conditions. These findings underscore the significance of considering temporal variations in assessing exoplanet habitability and expanding the search space for potentially habitable worlds. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ. 10 figures,16 pages

Three-Body Binary Formation in Clusters: Analytical Theory

Authors: Yonadav Barry Ginat, Hagai B. Perets

Abstract: Binary formation in clusters through triple encounters between three unbound stars, 'three-body' binary formation, is one of the main dynamical formation processes of binary systems in dense environments. In this paper, we use an analytical probabilistic approach to study the process for the equal mass case and calculate a probability distribution for the orbital parameters of three-body-formed bi… ▽ More Binary formation in clusters through triple encounters between three unbound stars, 'three-body' binary formation, is one of the main dynamical formation processes of binary systems in dense environments. In this paper, we use an analytical probabilistic approach to study the process for the equal mass case and calculate a probability distribution for the orbital parameters of three-body-formed binaries, as well as their formation rate. For the first time, we give closed-form analytical expressions to the full orbital parameter distribution, accounting for both energy and angular momentum conservation. This calculation relies on the sensitive dependence of the outcomes of three-body scatterings on the initial conditions: here we compute the rate of three-body binaries from ergodic interactions, which allow for an analytical derivation of the distribution of orbital parameters of the binaries thus created. We find that soft binaries are highly favoured in this process and that these binaries have a super-thermal eccentricity distribution, while the few hard three-body binaries have an eccentricity distribution much closer to thermal. The analytical results predict and reproduce simulation results of three-body scattering experiments in the literature well. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Comments welcome

Soft no more: gas shielding protects soft binaries from disruption in gas-rich environments

Authors: Mor Rozner, Hagai B. Perets

Abstract: Binaries in dense environments are traditionally classified as soft or hard based on their binding energy relative to the kinetic energy of surrounding stars. Heggie's law suggests that stellar encounters tend to soften soft binaries and harden hard binaries, altering their separations. However, interactions with gas in such environments can significantly modify this behavior. This study investiga… ▽ More Binaries in dense environments are traditionally classified as soft or hard based on their binding energy relative to the kinetic energy of surrounding stars. Heggie's law suggests that stellar encounters tend to soften soft binaries and harden hard binaries, altering their separations. However, interactions with gas in such environments can significantly modify this behavior. This study investigates the impact of gas on binary softening and its consequences. We find that gas interactions can actually harden binaries, extending the soft-hard boundary to larger separations. This introduces a "shielding radius" within which binaries are likely to harden due to gas interactions, surpassing the traditional soft-hard limit. Consequently, a notable portion of binaries initially classified as "soft" may become "hard" when both gas and stars are considered. We propose a two-stage formation process for hard binaries: initial soft binary formation, either dynamically or through gas-assisted capture, followed by gas-induced hardening before eventual disruption. In environments with low gas density but high gas content, the shielding radius could exceed the typical hard-soft limit by an order of magnitude, leading to a significant fraction of originally soft binaries effectively becoming hard. Conversely, in high gas-density environments, gas-induced hardening may dominate, potentially rendering the entire binary population hard. Gas hardening emerges as a crucial factor in shaping binary populations in gas-rich settings, such as clusters, star-forming regions, and possibly AGN disks. This highlights the complex interplay between gas dynamics and stellar interactions in binary evolution within dense environments. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: Comments are welcome!

Upper Limits of $^{44}$Ti Decay Emission in Four Nearby Thermonuclear Supernova Remnants

Authors: Jianbin Weng, Ping Zhou, Hagai B. Perets, Daniel R. Wik, Yang Chen

Abstract: To identify progenitors and investigate evidence of He burning, we searched for decay radiation of freshly synthesized $^{44}$Ti in four young nearby thermonuclear supernova remnants: Kepler, SN 1885, G1.9+0.3 and SN 1006, by analysing the up-to-date NuSTAR archival data. No apparent flux excess from the 68 and 78 keV line emissions accompanying decay was detected above the power law continuum app… ▽ More To identify progenitors and investigate evidence of He burning, we searched for decay radiation of freshly synthesized $^{44}$Ti in four young nearby thermonuclear supernova remnants: Kepler, SN 1885, G1.9+0.3 and SN 1006, by analysing the up-to-date NuSTAR archival data. No apparent flux excess from the 68 and 78 keV line emissions accompanying decay was detected above the power law continuum applied for the remnants and the absorbed stray light. By comparing the inferred upper limits of the line flux and the initial $^{44}$Ti masses with a wide variety of supernova nucleosynthesis models, we placed constraints on the supernova progenitors. We derived the first NuSTAR line flux upper limit for Kepler and ruled out most of the double-detonation scenarios with a thick He layer under low density. We estimated, for the first time, the upper limit for SN 1885, which is high because of the large distance yet still remains consistent with the He shell detonation. The new flux and mass limit of G1.9+0.3 derived from a longer total exposure is lower than the results from previous studies and evidently excludes explosive burning of He-rich matter. The relatively advanced age and the large spatial extent of SN 1006 have prevented meaningful constraints. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: 14 pages, 12 figures. Accepted to MNRAS

Morphology and Mach Number Dependence of Subsonic Bondi-Hoyle Accretion

Authors: Logan J. Prust, Hila Glanz, Lars Bildsten, Hagai B. Perets, Friedrich K. Roepke

Abstract: We carry out three-dimensional computations of the accretion rate onto an object (of size $R_{\rm sink}$ and mass $m$) as it moves through a uniform medium at a subsonic speed $v_{\infty}$. The object is treated as a fully-absorbing boundary (e.g. a black hole). In contrast to early conjectures, we show that when $R_{\rm sink}\ll R_{A}=2Gm/v^2$ the accretion rate is independent of $v_{\infty}$ and… ▽ More We carry out three-dimensional computations of the accretion rate onto an object (of size $R_{\rm sink}$ and mass $m$) as it moves through a uniform medium at a subsonic speed $v_{\infty}$. The object is treated as a fully-absorbing boundary (e.g. a black hole). In contrast to early conjectures, we show that when $R_{\rm sink}\ll R_{A}=2Gm/v^2$ the accretion rate is independent of $v_{\infty}$ and only depends on the entropy of the ambient medium, its adiabatic index, and $m$. Our numerical simulations are conducted using two different numerical schemes via the Athena++ and Arepo hydrodynamics solvers, which reach nearly identical steady-state solutions. We find that pressure gradients generated by the isentropic compression of the flow near the accretor are sufficient to suspend much of the surrounding gas in a near-hydrostatic equilibrium, just as predicted from the spherical Bondi-Hoyle calculation. Indeed, the accretion rates for steady flow match the Bondi-Hoyle rate, and are indicative of isentropic flow for subsonic motion where no shocks occur. We also find that the accretion drag may be predicted using the Safronov number, $Θ=R_{A}/R_{\rm sink}$, and is much less than the dynamical friction for sufficiently small accretors ($R_{\rm sink}\ll R_{A}$). △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: 14 pages, 7 figures. Submitted to ApJ. Comments welcome!

A Triple Scenario for the Formation of Wide Black Hole Binaries Such As Gaia BH1

Authors: A. Generozov, H. B. Perets

Abstract: Recently, several non-interacting black hole-stellar binaries have been identified in Gaia data. For example, Gaia BH1, where a Sun-like star is in moderate eccentricity (e=0.45), 185-day orbit around a black hole. This orbit is difficult to explain through binary evolution. The present-day separation suggests the progenitor system would have undergone an episode of common envelope evolution, but… ▽ More Recently, several non-interacting black hole-stellar binaries have been identified in Gaia data. For example, Gaia BH1, where a Sun-like star is in moderate eccentricity (e=0.45), 185-day orbit around a black hole. This orbit is difficult to explain through binary evolution. The present-day separation suggests the progenitor system would have undergone an episode of common envelope evolution, but a common envelope should shrink the period below the observed one. Since the majority of massive stars form in higher multiplicity systems, a triple evolution scenario is more likely for the progenitors of BH binaries. Here we show that such systems can indeed be more easily explained via evolution in hierarchical triple systems. von Zeipel--Lidov--Kozai oscillations or instabilities can delay the onset of the common envelope phase in the inner binary of the triple, so that the black hole progenitor and low-mass star are more widely separated when it begins, leading to the formation of wider binaries. There are also systems with similar periods but larger eccentricities, where the BH progenitor is a merger product of the inner binary in the triple. Such mergers lead to a more top-heavy black hole mass function. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 14 pages, 10 figures. To be submitted to ApJ

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.

Comments: Accepted for publication by A&A (with some revisions)

Journal ref: A&A 685, A45 (2024)

Dynamical evolution of second-generation circumstellar/protoplanetary disks in evolved wide binary systems

Authors: Raffaele S. Cattolico, Hagai B. Perets

Abstract: In mass-transferring wide binary stellar systems, the companion star can capture some of the mass released in wind by the primary evolved star, and form an accretion disk. Such accretion disks could evolve to form disks of comparable properties to protoplanetary disks and may enable the formation of new planets and/or the interactions, re-growth, and re-migration of pre-existing planets in the new… ▽ More In mass-transferring wide binary stellar systems, the companion star can capture some of the mass released in wind by the primary evolved star, and form an accretion disk. Such accretion disks could evolve to form disks of comparable properties to protoplanetary disks and may enable the formation of new planets and/or the interactions, re-growth, and re-migration of pre-existing planets in the newly formed disks. We study the formation and the dynamical evolution of such "second generation" (SG) protoplanetary disks in evolved wide binary systems, with primaries in the mass range of 1-8 Solar mass. We follow their evolution from the asymptotic giant branch (AGB) phase of the red giant stellar donor until its evolution to become a white dwarf. We perform 1D semi-analytical numerical simulations for several binary systems varying the mass of the evolved stellar donor and the initial orbital distance, taking into account the changing mass-loss rates and the binary orbital expansion due to mass-loss. We calculate the radial density profile of the formed SG accretion disk and its temperature profile, considering a non-stationary viscosity profile which depends on the radial temperature profile. We find that SG circumstellar disks evolve to form a long-lived stable structure over the lifetime of the donor star, and we show that we can consistently produce the observed conditions and accretion rate inferred for the Mira evolved wide-binary system. The quasi-steady state radial surface density profiles are comparable with the typical range of masses and densities of observed (regular "first-generation") protoplanetary disks. This suggests that realistic SG disks can give rise to a second phase of planet formation and dynamics in old wide binary systems. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Thermonuclear explosion criteria for direct and indirect collisions of CO white dwarfs: a study of the impact-parameter threshold for detonation

Authors: Hila Glanz, Hagai B. Perets, Ruediger Pakmor

Abstract: The physical collisions of two white dwarfs (WDs) (i.e. not slow mergers) have been shown to produce type-Ia-like supernovae (SNe) explosions. Most studies of WD-collisions have focused on zero impact-parameter (direct) collisions, which can also be studied in 2D. However, the vast majority of WD collisions arising from any evolutionary channels suggested to date, are expected to be indirect, i.e.… ▽ More The physical collisions of two white dwarfs (WDs) (i.e. not slow mergers) have been shown to produce type-Ia-like supernovae (SNe) explosions. Most studies of WD-collisions have focused on zero impact-parameter (direct) collisions, which can also be studied in 2D. However, the vast majority of WD collisions arising from any evolutionary channels suggested to date, are expected to be indirect, i.e. have a non-negligible impact parameter upon collision. Here we use the highest resolution 3D simulations to date (making use of the AREPO code), in order to explore both direct and indirect collisions and the conditions in which they give rise to a detonation and the production of a luminous SNe. Using our simulations, we find a detonation criterion that can provide the critical impact parameter for an explosion to occur, depending on the density profile of the colliding WDs, their composition, and their collision velocities. We find that the initial velocity has a significant impact on the amount of 56Ni production from the explosion. Furthermore, the 56Ni production is also strongly dependent on the numerical resolution. While our results from the head-on collision with large initial velocities produce more 56Ni than previous simulations, those with small and comparable velocities produce significantly less 56Ni than previous works. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: Comments are welcome

ULTRASAT: A wide-field time-domain UV space telescope

Authors: Y. Shvartzvald, E. Waxman, A. Gal-Yam, E. O. Ofek, S. Ben-Ami, D. Berge, M. Kowalski, R. Bühler, S. Worm, J. E. Rhoads, I. Arcavi, D. Maoz, D. Polishook, N. Stone, B. Trakhtenbrot, M. Ackermann, O. Aharonson, O. Birnholtz, D. Chelouche, D. Guetta, N. Hallakoun, A. Horesh, D. Kushnir, T. Mazeh, J. Nordin , et al. (19 additional authors not shown)

Abstract: The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is scheduled to be launched to geostationary orbit in 2026. It will carry a telescope with an unprecedentedly large field of view (204 deg$^2$) and NUV (230-290nm) sensitivity (22.5 mag, 5$σ$, at 900s). ULTRASAT will conduct the first wide-field survey of transient and variable NUV sources and will revolutionize our ability to study the hot… ▽ More The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is scheduled to be launched to geostationary orbit in 2026. It will carry a telescope with an unprecedentedly large field of view (204 deg$^2$) and NUV (230-290nm) sensitivity (22.5 mag, 5$σ$, at 900s). ULTRASAT will conduct the first wide-field survey of transient and variable NUV sources and will revolutionize our ability to study the hot transient universe: It will explore a new parameter space in energy and time-scale (months long light-curves with minutes cadence), with an extra-Galactic volume accessible for the discovery of transient sources that is $>$300 times larger than that of GALEX and comparable to that of LSST. ULTRASAT data will be transmitted to the ground in real-time, and transient alerts will be distributed to the community in $<$15 min, enabling a vigorous ground-based follow-up of ULTRASAT sources. ULTRASAT will also provide an all-sky NUV image to $>$23.5 AB mag, over 10 times deeper than the GALEX map. Two key science goals of ULTRASAT are the study of mergers of binaries involving neutron stars, and supernovae: With a large fraction ($>$50%) of the sky instantaneously accessible, fast (minutes) slewing capability and a field-of-view that covers the error ellipses expected from GW detectors beyond 2025, ULTRASAT will rapidly detect the electromagnetic emission following BNS/NS-BH mergers identified by GW detectors, and will provide continuous NUV light-curves of the events; ULTRASAT will provide early (hour) detection and continuous high (minutes) cadence NUV light curves for hundreds of core-collapse supernovae, including for rarer supernova progenitor types. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: 40 pages, 16 figures, 3 tables. Submitted to the AAS journals

Born to be wide: the distribution of wide binaries in the field and soft binaries in clusters

Authors: Mor Rozner, Hagai B. Perets

Abstract: Most stars, binaries, and higher multiplicity systems are thought to form in stellar clusters and associations, which later dissociate. Very wide binaries can be easily disrupted in clusters due to dynamical evaporation (soft binaries) and/or due to tidal disruption by the gravitational potential of the cluster. Nevertheless, wide binaries are quite frequent in the field, where they can sometimes… ▽ More Most stars, binaries, and higher multiplicity systems are thought to form in stellar clusters and associations, which later dissociate. Very wide binaries can be easily disrupted in clusters due to dynamical evaporation (soft binaries) and/or due to tidal disruption by the gravitational potential of the cluster. Nevertheless, wide binaries are quite frequent in the field, where they can sometimes play a key role in the formation of compact binaries, and serve as tools to study key physical processes. Here we use analytic tools to study the dynamical formation of soft binaries in clusters, and their survival as field binaries following cluster dispersion. We derive the expected properties of very wide binaries both in clusters and in the field. We analytically derive their detailed distributions, including wide-binary fraction as a function of mass in different cluster environments, binaries mass functions and mass ratios, and the distribution of their orbital properties. We show that our calculations agree well on most aspects with the results of N-body simulations, but show some different binary-fraction dependence on the cluster mass. We find that the overall fraction of wide binaries scales as $\propto N_\star^{-1}$ where $N_\star$ is the size of the cluster, even for non-equal mass stars. More massive stars are more likely to capture wide companions, with most stars above five solar mass likely to capture at least one stellar companion, and triples formation is found to be frequent. △ Less

Submitted 28 June, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: Accepted to publication in APJ

Capture of stars into gaseous discs around massive black holes: Alignment, circularization and growth

Authors: Aleksey Generozov, Hagai B. Perets

Abstract: The majority of massive black holes (MBHs) likely hosted gas discs during their lifetimes. These could either be long-lived active galactic nuclei (AGN) discs, or shorter-lived discs formed following singular gas infall events, as was likely the case in our own Galactic Center. Stars and compact objects in such environments are therefore expected to interact with the gaseous disc as they go throug… ▽ More The majority of massive black holes (MBHs) likely hosted gas discs during their lifetimes. These could either be long-lived active galactic nuclei (AGN) discs, or shorter-lived discs formed following singular gas infall events, as was likely the case in our own Galactic Center. Stars and compact objects in such environments are therefore expected to interact with the gaseous disc as they go through it, and potentially become aligned and fully embedded within it. The interactions of embedded stars with the gas could give rise to a plethora of physical processes affecting the stars, including growth through accretion of gas, migration in the disc, stellar captures, and mergers with other stars. The impact of such processes strongly depends on the population of stars that eventually align with the disc and become embedded in it. Here we make use of analytic tools to analyze the alignment process, accounting for both geometric drag and gas dynamical friction. We find that up to $\sim 50\%$ of main sequence stars and stellar mass black holes in the central 0.1 pc can align with AGN discs in the Galactic Center and similar galactic nuclei. The orbits of aligned stars are typically circularized and are prograde with respect to the AGN disc. Furthermore, alignment and accretion are intimately linked, and the capture of stars by an AGN disc can potentially explain the origin of the young stellar disc in the Galactic Center with a top-heavy mass function, even without the need for a star-formation event. △ Less

Submitted 11 April, 2024; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: ~10 pages, 8 figures + Appendices. This is an author-produced PDF of an article accepted for publication in MNRAS following peer-review. The version of record is available at DOI:10.1093/mnras/stad1016

Binary formation through gas-assisted capture and the implications for stellar, planetary and compact-object evolution

Authors: Mor Rozner, Aleksey Generozov, Hagai B. Perets

Abstract: Binary systems are ubiquitous and their formation requires two-body interaction and dissipation. In gaseous media, interactions between two initially unbound objects could result in gas-assisted binary formation, induced by a loss of kinetic energy to the ambient gas medium. Here we derive analytically the criteria for gas-assisted binary capture through gas dynamical friction dissipation. We vali… ▽ More Binary systems are ubiquitous and their formation requires two-body interaction and dissipation. In gaseous media, interactions between two initially unbound objects could result in gas-assisted binary formation, induced by a loss of kinetic energy to the ambient gas medium. Here we derive analytically the criteria for gas-assisted binary capture through gas dynamical friction dissipation. We validate them with few-body simulations and explore this process in different gas-rich environments: gas-embedded star-forming regions (SFR), gas-enriched globular clusters, AGN disks and protoplanetary-disks. We find that gas-assisted binary capture is highly efficient in SFRs, potentially providing a main channel for the formation of binaries. It could also operate under certain conditions in gas-enriched globular clusters. Thin AGN disks could also provide a fertile ground for gas-assisted binary capture and in particular the formation of black-hole/other compact object binaries, the production of gravitational-wave (GW) and other high-energy transients. Large-scale gaseous disks might be too thick to enable gas-assisted binary capture and previous estimates of the production of GW-sources could be overestimated, and sensitive to specific conditions and the structure of the disks. In protoplanetary-disks, while gas-assisted binary capture can produce binary KBOs, dynamical friction by small planetsimals is likely to be more efficient. Overall, we show that gas-assisted binary formation is robust and can contribute significantly to the binary formation rate in many environments. In fact, the gas-assisted binary capture rates are sufficiently high such that they will lead to multicaptures, and the formation of higher multiplicity systems. △ Less

Submitted 24 April, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: Accepted to MNRAS

Resonant Dynamical Friction Around a Super-Massive Black Hole: Analytical Description

Authors: Yonadav Barry Ginat, Taras Panamarev, Bence Kocsis, Hagai B. Perets

Abstract: We derive an analytical model for the so-called phenomenon of `resonant dynamical friction', where a disc of stars around a super-massive black hole interacts with a massive perturber, so as to align its inclination with the disc's orientation. We show that it stems from a singular behaviour of the orbit-averaged equations of motion, which leads to a rapid alignment of the argument of the ascendin… ▽ More We derive an analytical model for the so-called phenomenon of `resonant dynamical friction', where a disc of stars around a super-massive black hole interacts with a massive perturber, so as to align its inclination with the disc's orientation. We show that it stems from a singular behaviour of the orbit-averaged equations of motion, which leads to a rapid alignment of the argument of the ascending node $Ω$ of each of the disc stars, with that of the perturber, $Ω_{\rm p}$, with a phase-difference of $90^\circ$. This phenomenon occurs for all stars whose maximum possible $\dotΩ$ (maximised over all values of $Ω$ for all the disc stars), is greater than $\dotΩ_{\rm p}$; this corresponds approximately to all stars whose semi-major axes are less than twice that of the perturber. The rate at which the perturber's inclination decreases with time is proportional to its mass and is shown to be much faster than Chandrasekhar's dynamical friction. We find that the total alignment time is inversely proportional to the root of the perturber's mass. This persists until the perturber enters the disc. The predictions of this model agree with a suite of numerical $N$-body simulations which we perform to explore this phenomenon, for a wide range of initial conditions, masses, \emph{etc.}, and are an instance of a general phenomenon. Similar effects could occur in the context of planetary systems, too. △ Less

Submitted 7 August, 2023; v1 submitted 27 November, 2022; originally announced November 2022.

Comments: (v1) Submitted for publication. Comments welcome. (v2) Updated to match version accepted to MNRAS

Using $^{44}$Ti Emission to Differentiate Between Thermonuclear Supernova Progenitors

Authors: Daniel Kosakowski, Mark Ivan Ugalino, Robert Fisher, Or Graur, Alexey Bobrick, Hagai B. Perets

Abstract: The radiosotope $^{44}$Ti is produced through $α$-rich freezeout and explosive helium burning in type Ia supernovae (SNe Ia). In this paper, we discuss how the detection of $^{44}$Ti, either through late-time light curves of SNe Ia, or directly via gamma rays, can uniquely constrain the origin of SNe Ia. In particular, building upon recent advances in the hydrodynamical simulation of helium-ignite… ▽ More The radiosotope $^{44}$Ti is produced through $α$-rich freezeout and explosive helium burning in type Ia supernovae (SNe Ia). In this paper, we discuss how the detection of $^{44}$Ti, either through late-time light curves of SNe Ia, or directly via gamma rays, can uniquely constrain the origin of SNe Ia. In particular, building upon recent advances in the hydrodynamical simulation of helium-ignited double white dwarf binaries, we demonstrate that the detection of $^{44}$Ti in a nearby SN Ia or in a young galactic supernova remnant (SNR) can discriminate between the double-detonation and double-degenerate channels of sub-Chandrasekhar (sub-$M_{\rm Ch}$) and near-Chandrasekhar (near-$M_{\rm Ch}$) SNe Ia. In addition, we predict that the late-time light curves of calcium-rich transients are entirely dominated by $^{44}$Ti. △ Less

Submitted 19 October, 2022; originally announced October 2022.

The origins of Calcium-rich supernovae from disruptions of CO white-dwarfs by hybrid He-CO white-dwarfs

Authors: Yossef Zenati, Hagai B. Perets, Luc Dessart, Wynn V. Jacobson-Gal'an, Silvia Toonen, Armin Rest

Abstract: Calcium-rich (SN 2005E-like) explosions are very faint (typical -15.5, type I supernovae (SNe) showing strong Ca-lines, mostly observed in old stellar environments. Several models for such SNe had been explored and debated, but non were able to consistently reproduce the observed properties of Ca-rich SNe, nor their rates and host-galaxy distributions. Here we show that the disruptions of low-mass… ▽ More Calcium-rich (SN 2005E-like) explosions are very faint (typical -15.5, type I supernovae (SNe) showing strong Ca-lines, mostly observed in old stellar environments. Several models for such SNe had been explored and debated, but non were able to consistently reproduce the observed properties of Ca-rich SNe, nor their rates and host-galaxy distributions. Here we show that the disruptions of low-mass carbon-oxygen (CO) white dwarfs (WDs) by hybrid Helium-CO (HeCO) WDs during their merger could explain the origin and properties of such SNe. We make use of detailed multi-dimensional hydrodynamical-thermonuclear (FLASH) simulations models of HeCO-CO WD-WD mergers to characterize such explosions. We find that the accretion of CO material onto a HeCO-WD heats its He-shell and eventually leads to its "weak" detonation and ejection and the production of a sub-energetic $\sim10^{49}$erg Ca-rich SN, while leaving the CO core of the HeCO-WD intact as a hot remnant white dwarf, possibly giving rise to X-ray emission as it cools down. We model the detailed light-curve and spectra of such explosions to find an excellent agreement with observations of Ia/c Ca-rich SNe. In these models most of the He is burned, producing Ia/c Ca-rich SNe but higher He-enrichment levels might potentially also explain the Ib Ca-rich SNe. We thereby provide a viable, consistent model for the origins of Ca-rich. These findings can shed new light on the role of Ca-rich in the chemical evolution of galaxies and the intra-cluster medium, and their contribution to the observed 511 kev signal in the Galaxy originating from positrons produced from ${}^{44}\textrm{Ti}$ decay. Finally, the origins of such SNe points to the key-role of HeCO WDs as SN progenitors, and their potential role as progenitors of other thermonuclear SNe including normal type Ia SNe. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Submitted to ApJ, comments are welcome

Report number: AAS40760

The Late-Time Light Curves of Type Ia Supernovae: Confronting Models with Observations

Authors: Vishal Tiwari, Or Graur, Robert Fisher, Ivo Seitenzahl, Shing-Chi Leung, Ken'ichi Nomoto, Hagai Binyamin Perets, Ken Shen

Abstract: Type Ia supernovae (SNe Ia) play a crucial role as standardizable candles in measurements of the Hubble constant and dark energy. Increasing evidence points towards multiple possible explosion channels as the origin of normal SNe Ia, with possible systematic effects on the determination of cosmological parameters. We present, for the first time, a comprehensive comparison of publicly-available SN… ▽ More Type Ia supernovae (SNe Ia) play a crucial role as standardizable candles in measurements of the Hubble constant and dark energy. Increasing evidence points towards multiple possible explosion channels as the origin of normal SNe Ia, with possible systematic effects on the determination of cosmological parameters. We present, for the first time, a comprehensive comparison of publicly-available SN Ia model nucleosynthetic data with observations of late-time light curve observations of SN Ia events. These models span a wide range of white dwarf (WD) progenitor masses, metallicities, explosion channels, and numerical methodologies. We focus on the influence of $^{57}$Ni and its isobaric decay product $^{57}$Co in powering the late-time ($t > 1000$ d) light curves of SNe Ia. $^{57}$Ni and $^{57}$Co are neutron-rich relative to the more abundant radioisotope $^{56}$Ni, and are consequently a sensitive probe of neutronization at the higher densities of near-Chandrashekhar (near-$M_{\rm Ch}$) progenitor WDs. We demonstrate that observations of one SN Ia event, SN 2015F is only consistent with a sub-$M_{\rm Ch}$ WD progenitor. Observations of four other events (SN 2011fe, SN 2012cg, SN 2014J, SN2013aa) are consistent with both near-$M_{\rm Ch}$ and sub-$M_{\rm Ch}$ progenitors. Continued observations of late-time light curves of nearby SNe Ia will provide crucial information on the nature of the SN Ia progenitors. △ Less

Submitted 29 June, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: 13 pages, 5 figures, 3 tables. Accepted for publication in MNRAS. Comments welcome!

Are there any pristine comets? Constraints from pebble structure

Authors: Uri Malamud, Wolf A. Landeck, Dorothea Bischoff, Christopher Kreuzig, Hagai B. Perets, Bastian Gundlach, Jurgen Blum

Abstract: We show that if comets (or any small icy planetesimals such as Kuiper belt objects) are composed of pebble piles, their internal radiogenic as well as geochemical heating results in considerably different evolutionary outcomes compared to similar past studies. We utilize a 1D thermo-physical evolution code, modified to include state-of-the-art empirical measurements of pebble thermal conductivity… ▽ More We show that if comets (or any small icy planetesimals such as Kuiper belt objects) are composed of pebble piles, their internal radiogenic as well as geochemical heating results in considerably different evolutionary outcomes compared to similar past studies. We utilize a 1D thermo-physical evolution code, modified to include state-of-the-art empirical measurements of pebble thermal conductivity and compression, the latter obtained through a new laboratory experiment presented here for the first time. Results indicate that due to the low pebble thermal conductivity, the peak temperatures attained during evolution are much higher than in any previous study given the same formation time. Assuming meteoritic radiogenic abundances, we find that only extremely small, sub-kilometre comets have the potential to retain the primordial, uniform and thermally unprocessed composition from which they formed. Comets with radii in excess of about 20 km are typically swept by rapid and energetically powerful aqueous hydration reactions. Across the full range of comet sizes and formation times, evolutions result in the processing and differentiation of various volatile species, and a radially heterogeneous nucleus stucture. Our computations however also indicate that the assumed fraction of radionuclides is a pivotal free parameter, because isotopic analyses of the only available cometary samples suggest that no 26Al was ever present in comet 81P/Wild 2. We show that if comets formed early in the protoplanetary disc (within 1-3 Myr), the radionuclide abundances indeed must be much smaller than those typically assumed based on meteoritic samples. We discuss the importance of our findings for the formation, present-day attributes and future research of comets. △ Less

Submitted 31 May, 2022; originally announced June 2022.

Detailed properties of gravitational-wave mergers from flyby perturbations of wide binary black holes in the field

Authors: Yael Raveh, Erez Michaely, Hagai Binyamin Perets

Abstract: Wide black hole binaries (wide-BBHs; $\geqslant 10^3$ AU) in the field can be perturbed by random stellar flybys that excite their eccentricities. Once a wide binary is driven to a sufficiently small pericenter approach, gravitational wave (GW) emission becomes significant, and the binary inspirals and merges. In our previous study, using simplified models for wide-BBHs, we found that successive f… ▽ More Wide black hole binaries (wide-BBHs; $\geqslant 10^3$ AU) in the field can be perturbed by random stellar flybys that excite their eccentricities. Once a wide binary is driven to a sufficiently small pericenter approach, gravitational wave (GW) emission becomes significant, and the binary inspirals and merges. In our previous study, using simplified models for wide-BBHs, we found that successive flybys lead to significant merger fractions of wide-BBHs in less than Hubble time, making the flyby perturbation mechanism a relevant contributor to the production rate of GW-sources. However, the exact rates and detailed properties of the resulting GW sources depend on the wide binary progenitors. In this paper we use detailed population synthesis models for the initial wide-BBH population, considering several populations corresponding to different natal-kick models and metallicities, and then follow the wide-BBHs evolution due to flyby perturbations and GW-emission. We show that the cumulative effect of flybys is conductive for the production of GW sources in non-negligible rates of $1-20$ Gpc$^{-3}$ yr$^{-1}$, which are sensitive to the natal kicks model. Such rates are relevant to the observationally inferred rate. Our models, now derived from detailed population of binaries, provide the detailed properties of the produced GW-sources, including mass-functions and delay times. The GW mergers are circularized when enter the aLIGO band; have a preference for high velocity dispersion host galaxies (in particular ellipticals); have a relatively uniform delay-time distribution; and likely have mildly correlated (less than isolated evolution channels and more than dynamical channels) prograde spin-spin and spin-orbits. △ Less

Submitted 26 April, 2022; originally announced April 2022.

Comments: 13 pages, 8 figures

3D Hydrodynamical Simulations of Helium-Ignited Double-degenerate White Dwarf Mergers

Authors: Niranjan Roy, Vishal Tiwari, Alexey Bobrick, Daniel Kosakowski, Robert Fisher, Hagai B. Perets, Rahul Kashyap, Pablo Lorén-Aguilar, Enrique García-Berro

Abstract: The origins of type Ia supernovae (SNe Ia) are still debated. Some of the leading scenarios involve a double detonation in double white dwarf (WD) systems. In these scenarios, helium shell detonation occurs on top of a carbon-oxygen (CO) WD, which then drives the detonation of the CO-core, producing a SN Ia. Extensive studies have been done on the possibility of a double helium detonation, followi… ▽ More The origins of type Ia supernovae (SNe Ia) are still debated. Some of the leading scenarios involve a double detonation in double white dwarf (WD) systems. In these scenarios, helium shell detonation occurs on top of a carbon-oxygen (CO) WD, which then drives the detonation of the CO-core, producing a SN Ia. Extensive studies have been done on the possibility of a double helium detonation, following a dynamical helium mass-transfer phase onto a CO-WD. However, 3D self-consistent modeling of the double-WD system, the mass transfer, and the helium shell detonation have been little studied. Here we use 3D hydrodynamical simulations to explore this case in which a helium detonation occurs near the point of Roche lobe overflow of the donor WD and may lead to an SN Ia through the dynamically driven double-degenerate double-detonation (D6) mechanism. We find that the helium layer of the accreting primary WD does undergo a detonation, while the underlying carbon-oxygen core does not, leading to an extremely rapid and faint nova-like transient instead of a luminous SN Ia event. This failed core detonation suggests that D6 SNe Ia may be restricted to the most massive carbon-oxygen primary WDs. We highlight the nucleosynthesis of the long-lived radioisotope $^{44}$Ti during explosive helium burning, which may serve as a hallmark both of successful as well as failed D6 events which subsequently detonate as classical double-degenerate mergers. △ Less

Submitted 25 July, 2022; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: 11 pages, 4 figures, 1 table. Published in the Astrophysical Journal Letters

Journal ref: 2022 ApJL 932 L24

The Circumstellar Environments of Double-Peaked, Calcium-strong Supernovae 2021gno and 2021inl

Authors: Wynn Jacobson-Galán, Padma Venkatraman, Raffaella Margutti, David Khatami, Giacomo Terreran, Ryan J. Foley, Rodrigo Angulo, Charlotte R. Angus, Katie Auchettl, Peter K. Blanchard, Alexey Bobrick, Joe S. Bright, Cirilla D. Couch, David A. Coulter, Karoli Clever, Kyle W. Davis, Thomas de Boer, Lindsay DeMarchi, Sierra A. Dodd, David O. Jones, Jessica Johnson, Charles D. Kilpatrick, Nandita Khetan, Zhisen Lai, Danial Langeroodi , et al. (20 additional authors not shown)

Abstract: We present panchromatic observations and modeling of calcium-strong supernovae (SNe) 2021gno in the star-forming host galaxy NGC 4165 (D = 30.5 Mpc) and 2021inl in the outskirts of elliptical galaxy NGC 4923 (D = 80 Mpc), both monitored through the Young Supernova Experiment (YSE) transient survey. The multi-color light curves of both SNe show two peaks, the former peak being derived from shock co… ▽ More We present panchromatic observations and modeling of calcium-strong supernovae (SNe) 2021gno in the star-forming host galaxy NGC 4165 (D = 30.5 Mpc) and 2021inl in the outskirts of elliptical galaxy NGC 4923 (D = 80 Mpc), both monitored through the Young Supernova Experiment (YSE) transient survey. The multi-color light curves of both SNe show two peaks, the former peak being derived from shock cooling emission (SCE) and/or shock interaction with circumstellar material (CSM). The primary peak in SN 2021gno is coincident with luminous, rapidly decaying X-ray emission ($L_x = 5 \times 10^{41}$ erg s$^{-1}$) detected by Swift-XRT at $δt = 1$ day after explosion, this observation being the second ever detection of X-rays from a calcium-strong transient. We interpret the X-ray emission from SN 2021gno in the context of shock interaction with dense CSM that extends to $r < 3 \times 10^{14}$ cm. Based on modeling of the SN 2021gno X-ray spectrum, we calculate a CSM mass range of $M_{\rm CSM} = (0.3 - 1.6) \times 10^{-3}$ M$_{\odot}$ and particle densities of $n = (1-4) \times 10^{10}$ cm$^{-3}$. Radio non-detections of SN 2021gno indicate a low-density environment at larger radii ($r > 10^{16}$ cm) and a progenitor mass loss rate of $\dot{M} < 10^{-4}$ M$_{\odot}$ yr$^{-1}$, for $v_w = 500$ km s$^{-1}$. For radiation derived from SCE, modeling of the primary light curve peak in both SNe indicates an extended progenitor envelope mass and radius of $M_e = 0.02 - 0.05$ M$_{\odot}$ and $R_e = 30 - 230$ R$_{\odot}$. The explosion properties of SNe 2021gno and 2021inl suggest progenitor systems containing either a low-mass massive star or a white dwarf (WD), the former being unlikely for either object given the lack of star formation at both explosion sites. Furthermore, the progenitor environments of both SNe are consistent with explosion models for low-mass hybrid He/C/O WD + C/O WD binaries. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 33 pages, 19 figures. Submitted to ApJ. Comments welcome!

Binary evolution, gravitational-wave mergers and explosive transients in multiple-populations gas-enriched globular-clusters

Authors: Mor Rozner, Hagai B. Perets

Abstract: Most globular clusters (GCs) show evidence for multiple stellar populations, suggesting the occurrence of several distinct star-formation episodes. The large fraction of second population (2P) stars observed requires a very large 2P gaseous mass to have accumulated in the cluster core to form these stars. Hence the first population of stars (1P) in the cluster core has had to become embedded in 2P… ▽ More Most globular clusters (GCs) show evidence for multiple stellar populations, suggesting the occurrence of several distinct star-formation episodes. The large fraction of second population (2P) stars observed requires a very large 2P gaseous mass to have accumulated in the cluster core to form these stars. Hence the first population of stars (1P) in the cluster core has had to become embedded in 2P gas, just prior to the formation of later populations. Here we explore the evolution of binaries in ambient 2P gaseous media of multiple-population GCs. We mostly focus on black hole binaries and follow their evolution as they evolve from wide binaries towards short periods through interaction with ambient gas, followed by gravitational-wave (GW) dominated inspiral and merger. We show this novel GW-merger channel could provide a major contribution to the production of GW-sources. We consider various assumptions and initial conditions and calculate the resulting gas-mediated change in the population of binaries and the expected merger rates due to gas-catalyzed GW-inspirals. For plausible conditions and assumptions, we find an expected GW merger rate observable by aLIGO of the order of up to a few tens of $\rm{Gpc^{-3} yr^{-1}}$, and an overall range for our various models of $0.08-25.51 \ \rm{Gpc^{-3} yr^{-1}}$. Finally, our results suggest that the conditions and binary properties in the early stage of GCs could be critically affected by gas-interactions and may require a major revision in the current modeling of the evolution of GCs. △ Less

Submitted 31 May, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: accepted to APJ

Compact objects formation, retention, and growth through accretion onto gas-embedded white-dwarfs/neutron-stars in gas-enriched globular-clusters

Authors: Hagai B. Perets

Abstract: Observations of pulsars in globular clusters (GCs) give evidence that more >10-20% of neutron stars (NSs) ever formed in GCs were retained there. However, the velocity distribution of field pulsars peaks at 5-10 times the escape velocities of GCs. Consequently, only a small fraction of GC-NSs should have been retained, even accounting for low-velocity NSs formed through electron-capture supernovae… ▽ More Observations of pulsars in globular clusters (GCs) give evidence that more >10-20% of neutron stars (NSs) ever formed in GCs were retained there. However, the velocity distribution of field pulsars peaks at 5-10 times the escape velocities of GCs. Consequently, only a small fraction of GC-NSs should have been retained, even accounting for low-velocity NSs formed through electron-capture supernovae. Thus, too few low-velocity NSs should have been retained in GCs, giving rise to the NS retention problem in GCs. Here we suggest a novel solution, in which the progenitors of most GC-NSs were ONe white-dwarfs (WDs) that accreted ambient intra-cluster gas and formed low-velocity NSs through accretion induced collapse (AIC). The existence of an early gas-enriched environment in GCs is supported by observations of GC multiple stellar populations. It is thought that 10s-100s of Myrs after the formation of the first generation of stars, and after ONe-WDs were already formed, GCs were replenished with gas which formed a second generation of stars. Accretion of such replenished gas onto the ONe-WDs catalyzed the AIC processes. The number of AIC-formed NSs is then sufficient to explain the large number of NSs retained in GCs. Similar processes might also drive CO-WDs to produce type-Ia SNe or to merge and form NSs, and similarly drive NSs to AIC and mergers producing BHs. Moreover, the wide variety of gas-catalyzed binary mergers and explosive transients suggest to occur in the gas-rich environments of AGN disk could similarly, and even more efficiently, occur in second-generation gas in GCs. △ Less

Submitted 10 February, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: Revised following comments from the community. Submitted

Constraints on the origins of hypervelocity stars: velocity distribution, mergers and star-formation history

Authors: Aleksey Generozov, Hagai B. Perets

Abstract: In recent years surveys have identified several dozen B stars in the Milky Way halo moving faster than the local escape speed. The origin of most of these hypervelocity stars (HVSs) is still poorly constrained. Here we show that the velocity distribution, and in particular the deficiency in >700 km/s HVSs is inconsistent with binary disruptions by the massive black hole (MBH) in the Galactic Centr… ▽ More In recent years surveys have identified several dozen B stars in the Milky Way halo moving faster than the local escape speed. The origin of most of these hypervelocity stars (HVSs) is still poorly constrained. Here we show that the velocity distribution, and in particular the deficiency in >700 km/s HVSs is inconsistent with binary disruptions by the massive black hole (MBH) in the Galactic Centre. This conclusion holds in the full and empty loss cone regime, and for secular instabilities in eccentric disks. Accounting for multiple close encounters between binaries and the MBH, does not qualitatively change the results. Moreover, there is no observed counterpart population in the Galactic Centre that is consistent with the HVSs. The star-formation history could be tuned explain the HVS velocity distribution, but this tuning would produce a mismatch with the observed HVS flight times. Frequent stellar collisions of the binary components due to interactions with the MBH do not significantly impact the velocity distribution in the Galactic halo. Such collisions, however, can leave observable remnants in the Galactic Centre, and potentially explain the origins of G2-like dust clouds. △ Less

Submitted 30 December, 2022; v1 submitted 27 December, 2021; originally announced December 2021.

Comments: 11 pages, 10 figures, accepted to MNRAS. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. The version of record [MNRAS, 513, 4257] is available online at DOI:10.1093/mnras/stac1108

Chaotic dynamics of wide triples induced by galactic tides: a novel channel for producing compact binaries, mergers, and collisions

Authors: Evgeni Grishin, Hagai B. Perets

Abstract: Recent surveys show that wide ($>10^4$ AU) binaries and triples are abundant in the field. We study the long-term evolution of wide hierarchical triple systems and the role played by the Galactic tidal (GT) field. We find that when the timescales of the secular von-Ziepel-Lidov-Kozai and the GT oscillations are comparable, triple evolution becomes chaotic which leads to extreme eccentricities. Con… ▽ More Recent surveys show that wide ($>10^4$ AU) binaries and triples are abundant in the field. We study the long-term evolution of wide hierarchical triple systems and the role played by the Galactic tidal (GT) field. We find that when the timescales of the secular von-Ziepel-Lidov-Kozai and the GT oscillations are comparable, triple evolution becomes chaotic which leads to extreme eccentricities. Consequently, the close pericentre approaches of the inner-binary components lead to strong interactions, mergers and collisions. We use a novel secular evolution code to quantify the key parameters and carry out a population-synthesis study of low and intermediate-mass wide-orbit triples. We find that in $\sim9\%$ of low-mass wide-triples the inner main-sequence binaries collide or tidally-inspiral within $10\ \rm Gyr$, with direct collisions are $6$ times more likely to occur. For the intermediate-mass sample, $\sim7.6\%$ of the systems merge or inspiral with roughly equal probabilities. We discuss the relative fractions of different stellar merger/inspiral outcomes as a function of their evolutionary stage (Main-Sequence, MS; Red-Giant, RG; or White-Dwarf, WD), their transient electromagnetic signatures and the final products of the merger/inspiral. In particular, the rate of WD-WD direct-collisions that lead to type-Ia Supernovae is comparable to other dynamical channels and accounts for at most $0.1\%$ of the observed rate. RG inspirals provide a novel channel for the formation of eccentric common-envelope-evolution binaries. The catalysis of mergers/collisions in triples due to GT could explain a significant fraction, or even the vast majority, of blue-stragglers in the field, produce progenitors for cataclysmic-variables, and give-rise to mergers/collisions of double-RG binaries. △ Less

Submitted 10 March, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: Accepted to MNRAS

Inflated Eccentric Migration of Evolving Gas-Giants I: Accelerated Formation and Destruction of Hot and Warm Jupiters

Authors: Mor Rozner, Hila Glanz, Hagai B. Perets, Evgeni Grishin

Abstract: Hot and warm Jupiters (HJs and WJs correspondingly) are gas-giants orbiting their host stars at very short orbital periods ($P_{HJ}<10$ days; $10<P_{WJ}<200$ days). HJs and a significant fraction of WJs are thought to have migrated from an initially farther-out birth locations. While such migration processes have been extensively studied, the thermal evolution of gas-giants and its coupling to the… ▽ More Hot and warm Jupiters (HJs and WJs correspondingly) are gas-giants orbiting their host stars at very short orbital periods ($P_{HJ}<10$ days; $10<P_{WJ}<200$ days). HJs and a significant fraction of WJs are thought to have migrated from an initially farther-out birth locations. While such migration processes have been extensively studied, the thermal evolution of gas-giants and its coupling to the the migration processes are usually overlooked. In particular, gas-giants end their core-accretion phase with large radii and then contract slowly to their final radii. Moreover, intensive heating can slow the contraction at various evolutionary stages. The initial inflated large radii lead to faster tidal migration due to the strong dependence of tides on the radius. Here we explore this accelerated migration channel, which we term inflated eccentric migration, using a semi-analytical self-consistent modeling of the thermal-dynamical evolution of the migrating gas-giants, later validated by our numerical model (see a companion paper, paper II). We demonstrate our model for specific examples and carry a population synthesis study. Our results provide a general picture of the properties of the formed HJs\&WJs via inflated migration, and the dependence on the initial parameters/distributions. We show that tidal migration of gas-giants could occur far more rapidly then previously thought and lead to accelerated destruction and formation of HJs and enhanced formation rate of WJs. Accounting for the coupled thermal-dynamical evolution is therefore critical to the understanding of HJs/WJs formation, evolution and final properties of the population and play a key role in their migration process. △ Less

Submitted 20 April, 2022; v1 submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted to APJ

Inflated Eccentric Migration of evolving gas giants II: Numerical methodology and basic concepts

Authors: Hila Glanz, Mor Rozner, Hagai B. Perets, Evgeni Grishin

Abstract: Hot and Warm Jupiters (HJs&WJs) are gas-giant planets orbiting their host stars at short orbital periods, posing a challenge to their efficient in-situ formation. Therefore, most of the HJs&WJs are thought to have migrated from an initially farther-out birth locations. Current migration models, i.e disc-migration (gas-dissipation driven) and eccentric-migration (tidal evolution driven), fail to pr… ▽ More Hot and Warm Jupiters (HJs&WJs) are gas-giant planets orbiting their host stars at short orbital periods, posing a challenge to their efficient in-situ formation. Therefore, most of the HJs&WJs are thought to have migrated from an initially farther-out birth locations. Current migration models, i.e disc-migration (gas-dissipation driven) and eccentric-migration (tidal evolution driven), fail to produce the occurrence rate and orbital properties of HJs&WJs. Here we study the role of the thermal evolution and its coupling to tidal evolution. We use the AMUSE, numerical environment, and MESA, planetary evolution modeling, to model in detail the coupled internal and orbital evolution of gas-giants during their eccentric-migration. In a companion paper, we use a simple semi-analytic model, validated by our numerical model, and run a population-synthesis study. We consider the initially inflated radii of gas-giants (expected following their formation), as well study the effects of the potential slowed contraction and even re-inflation of gas-giants (due to tidal and radiative heating) on the eccentric-migration. Tidal forces that drive eccentric-migration are highly sensitive to the planetary structure and radius. Consequently, we find that this form of inflated eccentric-migration operates on significantly (up to an order of magnitude) shorter timescales than previously studied eccentric-migration models. Thereby, inflated eccentric-migration gives rise to more rapid formation of HJs&WJs, higher occurrence rates of WJs, and higher rates of tidal disruptions, compared with previous eccentric migration models which consider constant ~Jupiter radii for HJ&WJ progenitors. Coupled thermal-dynamical evolution of eccentric gas-giants can therefore play a key-role in their evolution. △ Less

Submitted 18 April, 2022; v1 submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJ

Spatially Resolved X-ray Study of Supernova Remnant G306.3$-$0.9 with Unusually High Calcium Abundance

Authors: Jianbin Weng, Ping Zhou, Yang Chen, Shing-Chi Leung, Silvia Toonen, Hagai B. Perets, Ken'ichi Nomoto, Yossef Zenati, Jacco Vink

Abstract: G306.3$-$0.9 is an asymmetric Galactic supernova remnant (SNR), whose progenitor has been thought to be a Type Ia supernova (SN), but its high Ca abundance appears inconsistent with the Type Ia origin. Hoping to uncover the reason for its asymmetry and the origin of this SNR, we performed a spatially resolved X-ray spectroscopic analysis of XMM-Newton and Chandra observation data. We divided the S… ▽ More G306.3$-$0.9 is an asymmetric Galactic supernova remnant (SNR), whose progenitor has been thought to be a Type Ia supernova (SN), but its high Ca abundance appears inconsistent with the Type Ia origin. Hoping to uncover the reason for its asymmetry and the origin of this SNR, we performed a spatially resolved X-ray spectroscopic analysis of XMM-Newton and Chandra observation data. We divided the SNR into 13 regions and analyzed the spectra using two-temperature models (0.2 keV + 1 keV). Compared to the southwestern regions, the northeastern regions have higher metal abundances and a lower gas density. This suggests that the asymmetric morphology results from the non-uniform ambient environment. We found that neither Type Ia nor core-collapse SN models can account for the abnormally high abundance ratios of Ar/Si, Ca/Si, or the shape of the abundance curve. A comparison with the Ca-rich transient models shows that G306.3$-$0.9 is likely to be the first identified Galactic "Ca-rich transient" remnant, although the theoretical production of element S is lower. We also note that the conclusion for the SNR's origin relies on the measured abundance ratios and existing nucleosynthesis models. Between two groups of Ca-rich transient explosion models, we prefer the He shell detonation for an accreting WD, rather than the merger of a white dwarf and a neutron star. △ Less

Submitted 22 October, 2021; v1 submitted 19 October, 2021; originally announced October 2021.

Comments: 21 pages, 10 figures. Accepted to ApJ

Binaries are softer than they seem: Effects of an external potential on the scattering dynamics of binaries

Authors: Yonadav Barry Ginat, Hagai B. Perets

Abstract: Binary evolution is influenced by dynamical scattering with other stars in dense environments. Heggie's law states that, due to their environments, hard binaries (whose orbital energy surpasses the typical energy of single stars) tend to harden (increase their orbital energy), while soft binaries typically soften. Here, we show that Heggie's law sometimes needs to be revised, by accounting for an… ▽ More Binary evolution is influenced by dynamical scattering with other stars in dense environments. Heggie's law states that, due to their environments, hard binaries (whose orbital energy surpasses the typical energy of single stars) tend to harden (increase their orbital energy), while soft binaries typically soften. Here, we show that Heggie's law sometimes needs to be revised, by accounting for an external potential, for example, for binaries in nuclear stellar discs or AGN discs, that are affected by the central massive black hole, or binary planetesimals in proto-planetary discs, affected by the host star. We find that in such environments, where the Hill radius is finite, binary-single scattering can have different outcomes. In particular, a three-body encounter could be cut short due to stars being ejected beyond the Hill radius, thereby ceasing to participate in further close interactions. This leads to a systematic difference in the energy changes brought about by the encounter, and in particular slows binary hardening, and even causes some hard binaries to soften, on average, rather than harden. We use our previously derived analytical, statistical solution to the bound chaotic three-body problem to quantitatively characterise the revision of the hardening-softening phase transition and evolution of binaries. We also provide an analytical calculation of the mean hardening rate of binaries in any environment (also reproducing the results of detailed N-body simulations). We show that the latter exhibits a non-trivial dependence on the Hill radius induced by the environment. △ Less

Submitted 5 October, 2021; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: Updated to match version accepted to MNRAS. Comments Welcome

Return of the TEDI: revisiting the Triple Evolution Dynamical Instability channel in triple stars

Authors: Adrian S. Hamers, Hagai B. Perets, Todd A. Thompson, Patrick Neunteufel

Abstract: Triple-star systems exhibit a phenomenon known as the Triple Evolution Dynamical Instability (TEDI), in which mass loss in evolving triples triggers short-term dynamical instabilities, potentially leading to collisions of stars, exchanges, and ejections. Previous work has shown that the TEDI is an important pathway to head-on stellar collisions in the Galaxy, significantly exceeding the rate of co… ▽ More Triple-star systems exhibit a phenomenon known as the Triple Evolution Dynamical Instability (TEDI), in which mass loss in evolving triples triggers short-term dynamical instabilities, potentially leading to collisions of stars, exchanges, and ejections. Previous work has shown that the TEDI is an important pathway to head-on stellar collisions in the Galaxy, significantly exceeding the rate of collisions due to random encounters in globular clusters. Here, we revisit the TEDI evolutionary pathway using state-of-the-art population synthesis methods that self-consistently take into account stellar evolution and binary interactions, as well as gravitational dynamics and perturbations from passing stars in the field. We find Galactic TEDI-induced collision rates on the order of 1e-4/yr, consistent with previous studies which were based on more simplified methods. The majority of TEDI-induced collisions involve main sequence stars, potentially producing blue straggler stars. Collisions are also possible involving more evolved stars, potentially producing eccentric post-common-envelope systems, and white dwarfs collisions leading to Type Ia supernovae (although the latter with a negligible contribution to the Galactic rate). In our simulations, the TEDI is not only triggered by adiabatic wind mass loss, but also by Roche lobe overflow in the inner binary: when the donor star becomes less massive than the accretor, the inner binary orbit widens, triggering triple dynamical instability. We find that collision rates are increased by ~17% when fly-bys in the field are taken into account. In addition to collisions, we find that the TEDI produces ~1e-4/yr of unbound stars, although none with escape speeds in excess of 1e3 km/s. △ Less

Submitted 2 December, 2021; v1 submitted 28 July, 2021; originally announced July 2021.

Comments: Accepted for publication in ApJ. 22 pages, 16 figures

Supernova explosions in active galactic nuclear discs

Authors: Evgeni Grishin, Alexey Bobrick, Ryosuke Hirai, Ilya Mandel, Hagai B. Perets

Abstract: Active galactic nuclei (AGN) are prominent environments for stellar capture, growth and formation. These environments may catalyze stellar mergers and explosive transients, such as thermonuclear and core-collapse supernovae (SNe). SN explosions in AGN discs generate strong shocks, leading to unique observable signatures. We develop an analytical model which follows the evolution of the shock propa… ▽ More Active galactic nuclei (AGN) are prominent environments for stellar capture, growth and formation. These environments may catalyze stellar mergers and explosive transients, such as thermonuclear and core-collapse supernovae (SNe). SN explosions in AGN discs generate strong shocks, leading to unique observable signatures. We develop an analytical model which follows the evolution of the shock propagating in the disc until it eventually breaks out. We derive the peak luminosity, bolometric lightcurve, and breakout time. The peak luminosities may exceed $10^{45}$ erg s$^{-1}$ and last from hours to days. The brightest explosions occur in regions of reduced density; either off-plane, or in discs around low-mass central black holes ($\sim 10^6\ M_\odot$), or in starved subluminous AGNs. Explosions in the latter two sites are easier to observe due to a reduced AGN background luminosity. We perform suites of 1D Lagrangian radiative hydrodynamics SNEC code simulations to validate our results and obtain the luminosity in different bands, and 2D axisymmetric Eulerian hydrodynamics code HORMONE simulations to study the morphology of the ejecta and its deviation from spherical symmetry. The observed signature is expected to be a bright blue, UV, or X-ray flare on top of the AGN luminosity from the initial shock breakout, while the subsequent red part of the lightcurve will largely be unobservable. We estimate the upper limit for the total event rate to be $\mathcal{R}\lesssim 100\ \rm yr^{-1}\ Gpc^{-3}$ for optimal conditions and discuss the large uncertainties in this estimate. Future high-cadence transient searches may reveal these events. Some existing tidal disruption event candidates may originate from AGN supernovae. △ Less

Submitted 12 July, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: Accepted to MNRAS

Common envelope evolution of eccentric binaries

Authors: Hila Glanz, Hagai B. Perets

Abstract: Common envelope evolution (CEE) is believed to be an important stage in the evolution of binary/multiple stellar systems. Following this stage, the CE is thought to be ejected, leaving behind a compact binary (or a merger product). Although extensively studied, the CEE process is still little understood, and although most binaries have non-negligible eccentricity, the effect of initial eccentricit… ▽ More Common envelope evolution (CEE) is believed to be an important stage in the evolution of binary/multiple stellar systems. Following this stage, the CE is thought to be ejected, leaving behind a compact binary (or a merger product). Although extensively studied, the CEE process is still little understood, and although most binaries have non-negligible eccentricity, the effect of initial eccentricity on the CEE has been little explored. Moreover, most studies assume a complete circularization of the orbit by the CE onset, while observationally such eccentricities are detected in many post-CE binaries. Here we use smoothed particle hydro-dynamical simulations (SPH) to study the evolution of initially eccentric ($0\le e\le0.95$) CE-systems. We find that initially eccentric binaries only partially circularize. In addition, higher initial eccentricity leads to a higher eccentricity following the end of the inspiral phase, with eccentricities as high as 0.18 in the most eccentric cases, and even higher if the initial peri-center of the orbit is located inside the star (e.g. following a kick into an eccentric orbit, rather than a smooth transition). CEE of more eccentric binaries leads to enhanced dynamical mass-loss of the CE compared with more circular binaries, and depends on the initial closest approach of the binary. We show that our results and the observed eccentricities of post-CE binaries suggest that the typical assumptions of circular orbits following CEE might potentially be revised. We expect post-CE eccentricities to affect the delay time distributions of various transients such as supernovae, gamma-ray bursts and gravitational-wave sources by up to tens of percents. △ Less

Submitted 8 August, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: Accepted to MNRAS

Transients from ONe White-Dwarf - Neutron-Star/Black-Hole Mergers

Authors: Alexey Bobrick, Yossef Zenati, Hagai B. Perets, Melvyn B. Davies, Ross Church

Abstract: We conduct the first 3D hydrodynamic simulations of oxygen-neon white dwarf-neutron star/black hole mergers (ONe WD-NS/BH mergers). Such mergers constitute a significant fraction, and may even dominate, the inspiral rates of all WD-NS binaries. We post-process our simulations to obtain the nuclear evolution of these systems and couple the results to a supernova spectral synthesis code to obtain th… ▽ More We conduct the first 3D hydrodynamic simulations of oxygen-neon white dwarf-neutron star/black hole mergers (ONe WD-NS/BH mergers). Such mergers constitute a significant fraction, and may even dominate, the inspiral rates of all WD-NS binaries. We post-process our simulations to obtain the nuclear evolution of these systems and couple the results to a supernova spectral synthesis code to obtain the first light curves and spectra for these transients. We find that the amount of $^{56}$Ni synthesised in these mergers grows as a strong function of the WD mass, reaching typically $0.05$ and up to $0.1\,{\rm M}_\odot$ per merger. Photodisintegration leads to similar amounts of $^4$He and about a ten times smaller amount of $^1$H. The nuclear yields from these mergers, in particular those of $^{55}$Mn, may contribute significantly to Galactic chemical evolution. The transients expected from ONe WD-NS mergers are dominantly red/infrared, evolve on month-long timescales and reach bolometric magnitudes of up to -16.5. The current surveys must have already detected these transients or are, alternatively, putting strong constraints on merger scenarios. The properties of the expected transients from WD-NS mergers best agree with faint type Iax supernovae. The Vera Rubin Observatory (LSST) will be detecting up to thousands of merging ONe WD-NS systems per year. We simulate a subset of our models with 2D axisymmetric FLASH code to investigate why they have been challenging for previous studies. We find that the likely main challenge has been effectively modelling the nuclear statistical equilibrium regime in such mergers. △ Less

Submitted 13 January, 2024; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: Accepted version, comments welcome

Probing supermassive stars and massive black hole seeds through gravitational wave inspirals

Authors: Yael Raveh, Yonadav Barry Ginat, Hagai B. Perets, Tyrone E. Woods

Abstract: We propose a novel source of gravitational wave emission: the inspirals of compact fragments inside primordial supermassive stars (SMSs). Such systems are thought to be an essential channel in the as-yet little understood formation of supermassive black holes (SMBHs). One model suggests that high accretion rates of $0.1$-1 M$_\odot$/yr attainable in atomically-cooled primordial halos can lead to t… ▽ More We propose a novel source of gravitational wave emission: the inspirals of compact fragments inside primordial supermassive stars (SMSs). Such systems are thought to be an essential channel in the as-yet little understood formation of supermassive black holes (SMBHs). One model suggests that high accretion rates of $0.1$-1 M$_\odot$/yr attainable in atomically-cooled primordial halos can lead to the formation of a nuclear-burning SMS. This will ultimately undergo collapse through a relativistic instability, leaving a massive BH remnant. Recent simulations suggest that supermassive stars rarely form in isolation, and that companion stars and even black holes formed may be captured/accreted and inspiral to the SMS core due to gas dynamical friction. Here, we explore the GW emission produced from such inspirals, which could probe the formation and evolution of SMS and seeds of the first supermassive black holes. We use a semi-analytic gas-dynamical friction model of the inspirals in the SMS to characterize their properties. We find such sources could potentially be observable by upcoming space-born GW-detectors at their formation redshifts with the benefit of gravitational lensing. Mergers within closely-related quasi-stars may produce a much stronger signal, though disambiguating such events from other high-z events may prove challenging. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 6 pages, 5 figures

Distinguishing Tidal Disruption Events from Impostors

Authors: Ann Zabludoff, Iair Arcavi, Stephanie La Massa, Hagai B. Perets, Benny Trakhtenbrot, B. Ashley Zauderer, Katie Auchettl, Jane L. Dai, K. Decker French, Tiara Hung, Erin Kara, Giuseppe Lodato, W. Peter Maksym, Yujing Qin, Enrico Ramirez-Ruiz, Nathaniel Roth, Jessie C. Runnoe, Thomas Wevers

Abstract: Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transi… ▽ More Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transient types, such as active galactic nuclei (AGN) and supernovae (SNe), which can make TDE classification ambiguous. A further complication is that "TDE signatures" have not been uniformly observed to similar sensitivities or even targeted across all candidates. This chapter reviews those events that are unusual relative to other TDEs, including the possibility of TDEs in pre-existing AGN, and summarises those characteristics thought to best distinguish TDEs from continuously accreting AGN, strongly flaring AGN, SNe, and Gamma-Ray Bursts (GRBs), as well as other potential impostors like stellar collisions, "micro-TDEs," and circumbinary accretion flows. We conclude that multiple observables should be used to classify any one event as a TDE. We also consider the TDE candidate population as a whole, which, for certain host galaxy or SMBH characteristics, is distinguishable statistically from non-TDEs, suggesting that at least some TDE candidates do in fact arise from SMBH-disrupted stars. △ Less

Submitted 8 April, 2021; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: 57 pages, 17 figures. Accepted for publication in Springer Space Science Reviews. Chapter in ISSI review "The Tidal Disruption of Stars by Massive Black Holes" vol. 79. Corrected several typos from previous submission

Thermonuclear explosion of a massive hybrid HeCO white-dwarf triggered by a He-detonation on a companion

Authors: Ruediger Pakmor, Yossef Zenati, Hagai B. Perets, Silvia Toonen

Abstract: Normal type Ia supernovae (SNe) are thought to arise from the thermonuclear explosion of massive ($>0.8$ M$_\odot$) carbon-oxygen white dwarfs (WDs), although the exact mechanism is debated. In some models helium accretion onto a carbon-oxygen (CO) WD from a companion was suggested to dynamically trigger a detonation of the accreted helium shell. The helium detonation then produces a shock that af… ▽ More Normal type Ia supernovae (SNe) are thought to arise from the thermonuclear explosion of massive ($>0.8$ M$_\odot$) carbon-oxygen white dwarfs (WDs), although the exact mechanism is debated. In some models helium accretion onto a carbon-oxygen (CO) WD from a companion was suggested to dynamically trigger a detonation of the accreted helium shell. The helium detonation then produces a shock that after converging on itself close to the core of the CO-WD, triggers a secondary carbon detonation and gives rise to an energetic explosion. However, most studies of such scenarios have been done in one or two dimensions, and/or did not consider self-consistent models for the accretion and the He-donor. Here we make use of detailed 3D simulation to study the interaction of a He-rich hybrid $0.69\,\mathrm{M_\odot}$ HeCO WD with a more massive $0.8\,\mathrm{M_\odot}$ CO~WD. We find that accretion from the hybrid WD onto the CO~WD gives rise to a helium detonation. However, the helium detonation does not trigger a carbon detonation in the CO~WD. Instead, the helium detonation burns through the accretion stream to also burn the helium shell of the donor hybrid HeCO-WD. The detonation of its massive helium shell then compresses its CO core, and triggers its detonation and full destruction. The explosion gives rise to a faint, likely highly reddened transient, potentially observable by the Vera Rubin survey, and the high-velocity ($\sim 1000\,\mathrm{km s^{-1}}$) ejection of the heated surviving CO~WD companion. Pending on uncertainties in stellar evolution we estimate the rate of such transient to be up to $\sim10\%$ of the rate of type Ia SNe. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: 14 pages, 10 figures, accepted by MNRAS, comments welcome

A Measurement of Stellar Surface Gravity Hidden in Radial Velocity Differences of Co-moving Stars

Authors: Matthew Moschella, Oren Slone, Jeff A. Dror, Matteo Cantiello, Hagai B. Perets

Abstract: The gravitational redshift induced by stellar surface gravity is notoriously difficult to measure for non-degenerate stars, since its amplitude is small in comparison with the typical Doppler shift induced by stellar radial velocity. In this study, we make use of the large observational data set of the Gaia mission to achieve a significant reduction of noise caused by these random stellar motions.… ▽ More The gravitational redshift induced by stellar surface gravity is notoriously difficult to measure for non-degenerate stars, since its amplitude is small in comparison with the typical Doppler shift induced by stellar radial velocity. In this study, we make use of the large observational data set of the Gaia mission to achieve a significant reduction of noise caused by these random stellar motions. By measuring the differences in velocities between the components of pairs of co-moving stars and wide binaries, we are able to statistically measure the combined effects of gravitational redshift and convective blueshifting of spectral lines, and nullify the effect of the peculiar motions of the stars. For the subset of stars considered in this study, we find a positive correlation between the observed differences in Gaia radial velocities and the differences in surface gravity and convective blueshift inferred from effective temperature and luminosity measurements. The results rule out a null signal at the $5σ$ level for our full data-set. Additionally, we study the sub-dominant effects of binary motion, and possible systematic errors in radial velocity measurements within Gaia. Results from the technique presented in this study are expected to improve significantly with data from the next Gaia data release. Such improvements could be used to constrain the mass-luminosity relation and stellar models which predict the magnitude of convective blueshift. △ Less

Submitted 20 June, 2022; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 7 pages, 4 figures MNRAS accepted version

Journal ref: MNRAS 514 (2022) 1071-1076

No velocity-kicks are required to explain large-distance offsets of Ca-rich supernovae and short-GRBs

Authors: Hagai B. Perets, Paz Beniamini

Abstract: The environments of explosive transients link their progenitors to the underlying stellar population, providing critical clues for their origins. However, some Ca-rich supernovae (SNe) and short gamma ray burst (sGRBs) appear to be located at remote locations, far from the stellar population of their host galaxy, challenging our understanding of their origin and/or physical evolution. These findin… ▽ More The environments of explosive transients link their progenitors to the underlying stellar population, providing critical clues for their origins. However, some Ca-rich supernovae (SNe) and short gamma ray burst (sGRBs) appear to be located at remote locations, far from the stellar population of their host galaxy, challenging our understanding of their origin and/or physical evolution. These findings instigated models suggesting that either large velocity kicks were imparted to their progenitors, allowing them to propagate to large distances and attain their remote locations; or that they formed in dense globular clusters residing in the halos. Here we show that instead, the large spatial-offsets of these transients are naturally explained by the observations of highly extended underlying stellar populations in (mostly early type) galaxy halos, typically missed since they can only be identified through ultra-deep/stacked images. Consequently, no large velocity kicks, nor halo globular cluster environments are required in order to explain the origin of these transients. These findings support thermonuclear explosions on white-dwarfs, for the origins of Ca-rich SNe progenitors, and no/small-natal-kicks given to sGRB progenitors. Since early-type galaxies contain older stellar populations, transient arising from older stellar populations would have larger fractions of early-type galaxy hosts, and consequently larger fractions of transients at large offsets. This is verified by our results for sGRBs and Ca-rich SNe showing different offset distributions in early vs. late-type galaxies. Furthermore, once divided between early and late type galaxies, the offsets' distributions of the different transients are consistent with each other. Finally, we point out that studies of other transients' offset distribution (e.g. Ia-SNe or FRBs) should similarly consider the host galaxy-type. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: Comments are welcome

Extreme mass-ratio gravitational-wave sources: Mass segregation and post binary tidal-disruption captures

Authors: Yael Raveh, Hagai B. Perets

Abstract: The gravitational-wave (GW) inspirals of stellar-mass compact objects onto a supermassive black hole (MBH), are some of the most promising GW sources detectable by next-generation space-born GW-detectors. The rates and characteristics of such extreme mass ratio inspirals (EMRIs) sources are highly uncertain. They are determined by the dynamics of stars near MBHs, and the rate at which compacts obj… ▽ More The gravitational-wave (GW) inspirals of stellar-mass compact objects onto a supermassive black hole (MBH), are some of the most promising GW sources detectable by next-generation space-born GW-detectors. The rates and characteristics of such extreme mass ratio inspirals (EMRIs) sources are highly uncertain. They are determined by the dynamics of stars near MBHs, and the rate at which compacts objects are driven to the close proximity of the MBH. Here we consider weakly and strongly mass-segregated nuclear clusters, and the evolution of stars captured into highly eccentric orbits following binary disruptions by the MBH. We make use of a Monte-Carlo approach to model the diffusion of both captured objects, and compact-objects brought through two-body relaxation processes. We calculate the rates of GW-inspirals resulting from relaxation-driven objects, and characterize EMRIs properties. We correct previous studies and show that relaxation-driven sources produce GW-sources with lower-eccentricity than previously found, and provide the detailed EMRI eccentricity distribution in the weak and strong mass-segregation regimes. We also show that binary-disruption captured-stars could introduce low-eccentricity GW-sources of stellar black-hole EMRIs in mass-segregated clusters. The eccentricities of the GW-sources from the capture channel, however, are strongly affected by relaxation processes, and are significantly higher than previously suggested. We find that both the rate and eccentricity distribution of EMRIs could probe the dynamics near MBHs, and the contribution of captured stars, characterize the mass-function of stellar compact objects, and verify whether weak or strong mass-segregation processes take place near MBHs. △ Less

Submitted 4 January, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 9 pages, 6 figures. Accepted for publication in MNRAS

Rapid destruction of planetary debris around white dwarfs through aeolian erosion

Authors: Mor Rozner, Dimitri Veras, Hagai B. Perets

Abstract: The discovery of numerous debris disks around white dwarfs (WDs), gave rise to extensive study of such disks and their role in polluting WDs, but the formation and evolution of these disks is not yet well understood. Here we study the role of aeolian (wind) erosion in the evolution of solids in WD debris disks. Aeolian erosion is a destructive process that plays a key role in shaping the propertie… ▽ More The discovery of numerous debris disks around white dwarfs (WDs), gave rise to extensive study of such disks and their role in polluting WDs, but the formation and evolution of these disks is not yet well understood. Here we study the role of aeolian (wind) erosion in the evolution of solids in WD debris disks. Aeolian erosion is a destructive process that plays a key role in shaping the properties and size-distribution of planetesimals, boulders and pebbles in gaseous protoplanetary disks. Our analysis of aeolian erosion in WD debris disks shows it can also play an important role in these environments. We study the effects of aeolian erosion under different conditions of the disk, and its erosive effect on planetesimals and boulders of different sizes. We find that solid bodies smaller than $\sim 5 \rm{km}$ will be eroded within the short disk lifetime. We compare the role of aeolian erosion in respect to other destructive processes such as collisional fragmentation and thermal ablation. We find that aeolian erosion is the dominant destructive process for objects with radius $\lesssim 10^3 \rm{cm}$ and at distances $\lesssim 0.6 R_\odot$ from the WD. Thereby, aeolian erosion constitutes the main destructive pathway linking fragmentational collisions operating on large objects with sublimation of the smallest objects and Poynting-Robertson drag, which leads to the accretion of the smallest particles onto the photosphere of WDs, and the production of polluted WDs. △ Less

Submitted 3 February, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: Accepted to MNRAS

Linking globular cluster structural parameters and their evolution: multiple stellar populations

Authors: Alessandra Mastrobuono-Battisti, Hagai B. Perets

Abstract: Globular clusters (GCs) are known to host multiple stellar populations showing chemical anomalies in the content of light elements. The origin of such anomalies observed in Galactic GCs is still debated. Here we analyse data compiled from the Hubble Space Telescope, ground-based surveys and Gaia DR2 and explore relationships between the structural properties of GCs and the fraction of second popul… ▽ More Globular clusters (GCs) are known to host multiple stellar populations showing chemical anomalies in the content of light elements. The origin of such anomalies observed in Galactic GCs is still debated. Here we analyse data compiled from the Hubble Space Telescope, ground-based surveys and Gaia DR2 and explore relationships between the structural properties of GCs and the fraction of second population (2P) stars. Given the correlations we find, we conclude that the main factor driving the formation/evolution of 2P stars is the cluster mass. The existing strong correlations between the 2P fraction and the rotational velocity and concentration parameter could derive from their correlation with the cluster mass. Furthermore, we observe that increasing cluster escape velocity corresponds to higher 2P fractions. Each of the correlations found is bimodal, with a different behaviour detected for low and high mass (or escape velocity) clusters. These correlations could be consistent with an initial formation of more centrally concentrated 2P stars in deeper cluster potentials, followed by a long-term tidal stripping of stars from clusters outskirts. The latter are dominated by the more extended distributed first population (1P) stars, and therefore stronger tidal stripping would preferentially deplete the 1P population, raising the cluster 2P fraction. This also suggests a tighter distribution of initial 2P fractions than observed today. In addition, higher escape velocities allow better retention of low-velocity material ejected from 1P stars, providing an alternative/additional origin for the observed differences and the distributions of 2P fractions amongst GCs. △ Less

Submitted 18 May, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: 18 pages, 8 figures (5 in the main text, 3 in the appendix), 7 tables (3 in the main text, 4 in the appendix). Accepted for publication in MNRAS

An Analytical, Statistical Approximate Solution for Dissipative and non-Dissipative Binary-Single Stellar Encounters

Authors: Yonadav Barry Ginat, Hagai B. Perets

Abstract: We present a statistical approximate solution of the bound, non-hierarchical three-body problem, and extend it to a general analysis of encounters between hard binary systems and single stars. Any such encounter terminates when one of the three stars is ejected to infinity, leaving behind a remnant binary; the problem of binary-single star-scattering consists of finding the probability distributio… ▽ More We present a statistical approximate solution of the bound, non-hierarchical three-body problem, and extend it to a general analysis of encounters between hard binary systems and single stars. Any such encounter terminates when one of the three stars is ejected to infinity, leaving behind a remnant binary; the problem of binary-single star-scattering consists of finding the probability distribution of the orbital parameters of the remnant binary, as a function of the total energy and the total angular momentum. Here, we model the encounter as a series of close, non-hierarchical, triple approaches, interspersed with hierarchical phases, in which the system consists of an inner binary and a star that orbits it -- this turns the evolution of the entire encounter to a random walk between consecutive hierarchical phases. We use the solution of the bound, non-hierarchical three-body problem to find the walker's transition probabilities, which we generalise to situations in which tidal interactions are important. Besides tides, any dissipative process may be incorporated into the random walk model, as it is completely general. Our approximate solution can reproduce the results of the extensive body of past numerical simulations, and can account for different environments and different dissipative effects. Therefore, this model can effectively replace the need for direct few-body integrations for the study of binary-single encounters in any environment. Furthermore, it allows for a simply inclusion of dissipative forces typically not accounted for in full N-body integration schemes. △ Less

Submitted 18 June, 2021; v1 submitted 30 October, 2020; originally announced November 2020.

Comments: (v1) Comments welcome. (v3) Updated to match version accepted to PRX

Journal ref: Phys. Rev. X 11, 031020 (2021)

Nuclear Burning in Collapsar Accretion Disks

Authors: Yossef Zenati, Daniel M. Siegel, Brian D. Metzger, Hagai B. Perets

Abstract: The core collapse of massive, rapidly-rotating stars are thought to be the progenitors of long-duration gamma-ray bursts (GRB) and their associated hyper-energetic supernovae (SNe). At early times after the collapse, relatively low angular momentum material from the infalling stellar envelope will circularize into an accretion disk located just outside the black hole horizon, resulting in high acc… ▽ More The core collapse of massive, rapidly-rotating stars are thought to be the progenitors of long-duration gamma-ray bursts (GRB) and their associated hyper-energetic supernovae (SNe). At early times after the collapse, relatively low angular momentum material from the infalling stellar envelope will circularize into an accretion disk located just outside the black hole horizon, resulting in high accretion rates necessary to power a GRB jet. Temperatures in the disk midplane at these small radii are sufficiently high to dissociate nuclei, while outflows from the disk can be neutron-rich and may synthesize r-process nuclei. However, at later times, and for high progenitor angular momentum, the outer layers of the stellar envelope can circularize at larger radii $\gtrsim 10^{7}$ cm, where nuclear reaction can take place in the disk midplane ((e.g.~$^{4}$He + $^{16}$O $\rightarrow$ $^{20}$Ne + $γ$).. Here we explore the effects of nuclear burning on collapsar accretion disks and their outflows by means of hydrodynamical $α$-viscosity torus simulations coupled to a 19-isotope nuclear reaction network, which are designed to mimic the late infall epochs in collapsar evolution when the viscous time of the torus has become comparable to the envelope fall-back time. Our results address several key questions, such as the conditions for quiescent burning and accretion versus detonation and the generation of $^{56}$Ni in disk outflows, which we show could contribute significantly to powering GRB supernovae. Being located in the slowest, innermost layers of the ejecta, the latter could provide the radioactive heating source necessary to make the spectral signatures of r-process elements visible in late-time GRB-SNe spectra. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: Submitted to MNRAS, comments are welcome

The Wide-Binary Origin of The Pluto-Charon System

Authors: Mor Rozner, Evgeni Grishin, Hagai B. Perets

Abstract: The Pluto-Charon binary system is the best-studied representative of the binary Kuiper-belt population. Its origins are vital to understanding the formation of other Kupier-belt objects (KBO) and binaries, and the evolution of the outer solar-system. The Pluto-Charon system is believed to form following a giant impact between two massive KBOs at relatively low velocities. However, the likelihood o… ▽ More The Pluto-Charon binary system is the best-studied representative of the binary Kuiper-belt population. Its origins are vital to understanding the formation of other Kupier-belt objects (KBO) and binaries, and the evolution of the outer solar-system. The Pluto-Charon system is believed to form following a giant impact between two massive KBOs at relatively low velocities. However, the likelihood of a random direct collision between two of the most massive KBOs is low, and is further constrained by the requirement of a low-velocity collision, making this a potentially fine-tuned scenario. Here we expand our previous studies and suggest that the proto-Pluto-Charon system was formed as a highly inclined wide-binary, which was then driven through secular/quasi-secular evolution into a direct impact. Since wide-binaries are ubiquitous in the Kuiper-belt with many expected to be highly inclined, our scenario is expected to be robust. We use analytic tools and few-body simulations of the triple Sun-(proto-)Pluto-Charon system to show that a large parameter-space of initial conditions leads to such collisions. The velocity of such an impact is the escape velocity of a bound system, which naturally explains the low-velocity impact. The dynamical evolution and the origins of the Pluto-Charon system could therefore be traced to similar secular origins as those of other binaries and contact-binaries (e.g. Arrokoth), and suggest they play a key role in the evolution of KBOs. △ Less

Submitted 12 August, 2020; v1 submitted 20 July, 2020; originally announced July 2020.

The former companion of the hyper-velocity star S5-HVS1

Authors: Wenbin Lu, Jim Fuller, Yael Raveh, Hagai B. Perets, Ting S. Li, Matthew W. Hosek Jr., Tuan Do

Abstract: The hyper-velocity star S5-HVS1, ejected 5 Myr ago from the Galactic Center at 1800 km/s, was most likely produced by tidal break-up of a tight binary by the supermassive black hole SgrA*. Taking a Monte Carlo approach, we show that the former companion of S5-HVS1 was likely a main-sequence star between 1.2 and 6 solar masses and was captured into a highly eccentric orbit with pericenter distance… ▽ More The hyper-velocity star S5-HVS1, ejected 5 Myr ago from the Galactic Center at 1800 km/s, was most likely produced by tidal break-up of a tight binary by the supermassive black hole SgrA*. Taking a Monte Carlo approach, we show that the former companion of S5-HVS1 was likely a main-sequence star between 1.2 and 6 solar masses and was captured into a highly eccentric orbit with pericenter distance in the range 1-10 AU and semimajor axis about $10^3$ AU. We then explore the fate of the captured star. We find that the heat deposited by tidally excited stellar oscillation modes leads to runaway disruption if the pericenter distance is smaller than about 3 AU. Over the past 5 Myr, its angular momentum has been significantly modified by orbital relaxation, which may stochastically drive the pericenter inwards below 3 AU and cause tidal disruption. We find an overall survival probability in the range 5% to 50%, depending on the local relaxation time in the close environment of the captured star, and the initial pericenter at capture. The pericenter distance of the surviving star has migrated to 10-100 AU, making it potentially the most extreme member of the S-star cluster. From the ejection rate of S5-HVS1-like stars, we estimate that there may currently be a few stars in such highly eccentric orbits. They should be detectable (typically Ks < 18.5 mag) by the GRAVITY instrument and by future Extremely Large Telescopes and hence provide an extraordinary probe of the spin of SgrA*. △ Less

Submitted 13 February, 2021; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: MNRAS accepted after minor revision

Erosion driven size-redistribution of protoplanetary disk solids and the onset of streaming Instability and Pebble Accretion

Authors: Evgeni Grishin, Mor Rozner, Hagai B. Perets

Abstract: The formation of the first planetesimals and the final growth of planetary cores relies on the abundance of small pebbles. The efficiencies of both the streaming instability (SI) process, suggested to catalyze the early growth of planetesimals, and the pebble-accretion process, suggested to accelerate the growth of planetary cores, depend on the sizes of solids residing in the disk. In particular,… ▽ More The formation of the first planetesimals and the final growth of planetary cores relies on the abundance of small pebbles. The efficiencies of both the streaming instability (SI) process, suggested to catalyze the early growth of planetesimals, and the pebble-accretion process, suggested to accelerate the growth of planetary cores, depend on the sizes of solids residing in the disk. In particular, these processes were found to be sensitive to size distribution of solids, and efficient planetesimal formation and growth through these channels require a limited pebble size distribution. Here we show that aeolian erosion, a process that efficiently grinds down boulders into a mono-sized distribution of pebbles, provides a natural upper limit for the maximal pebble sizes (in terms of their Stokes number). We find the dependence of this upper limit on the radial separation, disk age, turbulence strength, and the grain-size composition of the boulders in the disk. SI is favorable in areas with a Stokes number less than 0.1, which is found in the inner sub-astronomical-unit regions of the disk. This upper limit shapes the size distribution of small pebbles and thereby catalyzes the early onset of planetesimal formation due to SI, and the later core accretion growth through pebble accretion. △ Less

Submitted 21 July, 2020; v1 submitted 7 April, 2020; originally announced April 2020.

Comments: Published in ApJL

Journal ref: The Astrophysical Journal Letters, 898:L13 (6pp), 2020 July 20

Simulations of common envelope evolution in triple systems: Circumstellar case

Authors: Hila Glanz, Hagai B. Perets

Abstract: The dynamical evolution of triple stellar systems could induce the formation of compact binaries and binary mergers. Common envelope (CE) evolution, which plays a major role in the evolution of compact binary systems, can similarly play a key role in the evolution of triples. Here we use hydrodynamical simulations coupled with few-body dynamics to provide the first detailed models of triple common… ▽ More The dynamical evolution of triple stellar systems could induce the formation of compact binaries and binary mergers. Common envelope (CE) evolution, which plays a major role in the evolution of compact binary systems, can similarly play a key role in the evolution of triples. Here we use hydrodynamical simulations coupled with few-body dynamics to provide the first detailed models of triple common envelope (TCE) evolution. We focus on the circumstellar case, where the envelope of an evolved giant engulfs a compact binary orbiting the giant, which then in-spirals into the core of the evolved star. Through our exploratory modeling we find several possible outcomes of such TCE: (1) The merger of the binary inside the third star's envelope; (2) The disruption of the in-spiraling binary following its plunge, leading to a chaotic triple dynamics of the stellar-core and the two components of the former disrupted binary. The chaotic evolution typically leads to the in-spiral and merger of at least one of the former binary components with the core, and sometimes to the ejection of the second, or alternatively its further now-binary common-envelope evolution. The in-spiral in TCE leads to overall slower in-spiral, larger mass ejection and the production of more aspherical remnant, compared with a corresponding binary case of similar masses, due to the energy/momentum extraction from the inner-binary. We expect TCE to play a key role in producing various types of stellar-mergers and unique compact binary systems, and potentially induce transient electromagnetic and gravitational-wave sources. △ Less

Submitted 15 October, 2020; v1 submitted 31 March, 2020; originally announced April 2020.

Comments: Accepted to MNRAS. Comments are welcomed

Rates of Stellar Tidal Disruption

Authors: Nicholas C. Stone, Eugene Vasiliev, Michael Kesden, Elena M. Rossi, Hagai B. Perets, Pau Amaro-Seoane

Abstract: Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astr… ▽ More Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: Accepted for publication in Springer Space Science Reviews. Chapter in ISSI review, "The Tidal Disruption of Stars by Massive Black Holes," vol. 79

Origin of (2014) MU69-like Kuiper-belt contact binaries from wide binaries

Authors: Evgeni Grishin, Uri Malamud, Hagai B. Perets, Oliver Wandel, Christoph M. Schaefer

Abstract: Following its flyby and first imaging the Pluto-Charon binary, the New Horizons spacecraft visited the Kuiper-Belt-Object (KBO) (486958) 2014 MU69 (Arrokoth). Imaging showed MU69 to be a contact-binary, made of two individual lobes connected by a narrow neck, rotating at low spin period (15.92 h), and having high obliquity (~98 deg), similar to other KBO contact-binaries inferred through photometr… ▽ More Following its flyby and first imaging the Pluto-Charon binary, the New Horizons spacecraft visited the Kuiper-Belt-Object (KBO) (486958) 2014 MU69 (Arrokoth). Imaging showed MU69 to be a contact-binary, made of two individual lobes connected by a narrow neck, rotating at low spin period (15.92 h), and having high obliquity (~98 deg), similar to other KBO contact-binaries inferred through photometric observations. The origin of such peculiar configurations is puzzling, and all scenarios suggested for the origins of contact-binaries fail to reproduce such properties and their likely high frequency. Here we show that semi-secular perturbations operating only on ultra-wide (~0.1-0.4 Hill-radius) KBO-binaries can robustly lead to gentle, slow-speed binary mergers at arbitrarily high obliquities, but low rotational velocities, that can reproduce MU69's (and similar oblique contact binaries) characteristics. Using N-body simulations, we find that ~15% of all ultra-wide binaries with cosine-uniform inclination distribution are likely to merge through this process. Moreover, we find that such mergers are sufficiently gentle as to only slightly deform the KBO shape, and can produce the measured rotation speed of MU69. The semi-secular contact-binary formation channel not only explains the observed properties of MU69, but could also apply for other Kuiper/asteroid belt binaries, and for Solar/extra-solar moon systems. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: Submitted on 27 June 2019

Journal ref: Nature, 580, 463 (2020)