Showing 1–50 of 259 results for author: Lai, D

Migration of Accreting Planets and Black Holes in Disks

Authors: JT Laune, Rixin Li, Dong Lai

Abstract: Nascent planets are thought to lose angular momentum (AM) to the gaseous protoplanetary disk (PPD) via gravitational interactions, leading to inward migration. A similar migration process also applies to stellar-mass black holes (BHs) embedded in AGN disks. However, AM exchange via accretion onto the planet/BH may strongly influence the migration torque. In this study, we perform 2D global hydrody… ▽ More Nascent planets are thought to lose angular momentum (AM) to the gaseous protoplanetary disk (PPD) via gravitational interactions, leading to inward migration. A similar migration process also applies to stellar-mass black holes (BHs) embedded in AGN disks. However, AM exchange via accretion onto the planet/BH may strongly influence the migration torque. In this study, we perform 2D global hydrodynamic simulations of an accreting planet/BH embedded in a disk, where AM exchange between the planet/BH and disk via gravity, accretion, pressure and viscosity are considered. When accretion is turned off, we recover the linear estimate for Type I migration torque. However, in all of our accreting simulations, we find outward migration due to the positive AM deposited onto the accreting body by the disk gas. Our simulations achieve the global steady state for the transport of mass and AM: The mass and AM fluxes are constant across the disk except for jumps ($Δ\dot M$ and $Δ\dot J$) at the planet's location, and the jumps match the accretion rate and torque on the planet. Our findings suggest that extra caution is needed when applying the standard results of disk migration to accreting planets and BHs. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: Submitted to ApJ

Free-Floating Planets, Survivor Planets, Captured Planets and Binary Planets from Stellar Flybys

Authors: Fangyuan Yu, Dong Lai

Abstract: In star clusters, close stellar encounters can strongly impact the architecture of a planetary system or even destroy it. We present a systematic study on the effects of stellar flybys on two-planet systems. When such a system (with a modest initial planet semi-major axis ratio, $a_1/a_2=0.6-0.8$) experiences stellar flybys, one or both planets can be ejected, forming free-floating planets (FFPs),… ▽ More In star clusters, close stellar encounters can strongly impact the architecture of a planetary system or even destroy it. We present a systematic study on the effects of stellar flybys on two-planet systems. When such a system (with a modest initial planet semi-major axis ratio, $a_1/a_2=0.6-0.8$) experiences stellar flybys, one or both planets can be ejected, forming free-floating planets (FFPs), captured planets (CPs) around the flyby star, and free-floating binary planets (BPs); the remaining single-surviving-planets (SSPs) can have their orbital radii and eccentricities greatly changed. Through numerical experiments, we calculate the formation fractions (or branching ratios) of FFPs, SSPs, CPs and BPs as a function of the pericenter separation of the flyby (in units of the initial planet semi-major axis), and use them to derive the analytical expressions that can be used to compute the formation rates of FFPs, SSPs, CPs and BPs in general cluster environments. We find that the production rates of FFPs and SSPs are similar, while the rate for CPs is a factor of a few smaller. The formation fraction of free-floating BPs depends strongly on $a_1/a_2$ of the initial systems and on the planet masses. For Jupiter-mass planets, the formation fraction of BPs is always less than $1\%$ (for systems with $a_1/a_2=0.8$) and typically much smaller ($\lesssim 0.2\%$ for $a_1/a_2\lesssim 0.7$). The fraction remains less than $1\%$ when considering $4M_{\rm J}$ planets. Overall, when averaging over all possible flybys, the production rate of BPs is less than $0.1\%$ of that for FFPs. We also derive the velocity distribution of FFPs produced by stellar flybys, and the orbital semi-major axis and eccentricity distributions of SSPs, CPs and (rare) free-floating BPs. These results can find applications in future studies of exotic planets and planetary systems. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 19 pages, 30 figures, submitted to ApJ

Extreme resonant eccentricity excitation of stars around merging black-hole binary

Authors: Bin Liu, Dong Lai

Abstract: We study the dynamics of a star orbiting a merging black-hole binary (BHB) in a coplanar triple configuration. During the BHB's orbital decay, the system can be driven across the apsidal precession resonance, where the apsidal precession rate of the stellar orbit matches that of the inner BHB. As a result, the system gets captured into a state of resonance advection until the merger of the BHB, le… ▽ More We study the dynamics of a star orbiting a merging black-hole binary (BHB) in a coplanar triple configuration. During the BHB's orbital decay, the system can be driven across the apsidal precession resonance, where the apsidal precession rate of the stellar orbit matches that of the inner BHB. As a result, the system gets captured into a state of resonance advection until the merger of the BHB, leading to an extreme eccentricity growth of the stellar orbit. This resonance advection occurs when the inner binary has a non-zero eccentricity and unequal masses. The resonant driving of the stellar eccentricity can significantly alter the hardening rate of the inner BHB, and produce observational signatures to uncover the presence of nearby merging or merged BHBs. △ Less

Submitted 7 June, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 5 pages, 5 figures, 1 supplemental material. Accepted for publication in PRL

Mass Ejection Driven by Sudden Energy Deposition in Stellar Envelopes

Authors: Nicholas Corso, Dong Lai

Abstract: A number of stellar astrophysical phenomena, such as tidal novae and planetary engulfment, involve sudden injection of sub-binding energy in a thin layer within the star, leading to mass ejection of the stellar envelope. We use a 1D hydrodynamical model to survey the stellar response and mass loss for various amounts ($E_{\mathrm{dep}}$) and locations of the energy deposition. We find that the tot… ▽ More A number of stellar astrophysical phenomena, such as tidal novae and planetary engulfment, involve sudden injection of sub-binding energy in a thin layer within the star, leading to mass ejection of the stellar envelope. We use a 1D hydrodynamical model to survey the stellar response and mass loss for various amounts ($E_{\mathrm{dep}}$) and locations of the energy deposition. We find that the total mass ejection has a nontrivial dependence on $E_{\mathrm{dep}}$ due to the varying strengths of mass ejection events, which are associated with density/pressure waves breaking out from the stellar surface. The rapid occurrence of multiple breakouts may present a unique observational signature for sudden envelope heating events in stars. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Comments: 8 pages, 6 figures, submitted to ApJ, awaiting review

Potential Melting of Extrasolar Planets by Tidal Dissipation

Authors: Darryl Z. Seligman, Adina D. Feinstein, Dong Lai, Luis Welbanks, Aster G. Taylor, Juliette Becker, Fred C. Adams, Marvin Morgan, Jennifer B. Bergner

Abstract: Tidal heating on Io due to its finite eccentricity was predicted to drive surface volcanic activity, which was subsequently confirmed by the $\textit{Voyager}$ spacecrafts. Although the volcanic activity in Io is more complex, in theory volcanism can be driven by runaway melting in which the tidal heating increases as the mantle thickness decreases. We show that this runaway melting mechanism is g… ▽ More Tidal heating on Io due to its finite eccentricity was predicted to drive surface volcanic activity, which was subsequently confirmed by the $\textit{Voyager}$ spacecrafts. Although the volcanic activity in Io is more complex, in theory volcanism can be driven by runaway melting in which the tidal heating increases as the mantle thickness decreases. We show that this runaway melting mechanism is generic for a composite planetary body with liquid core and solid mantle, provided that (i) the mantle rigidity, $μ$, is comparable to the central pressure, i.e. $μ/ (ρg R_{\rm P})\gtrsim0.1$ for a body with density $ρ$, surface gravitational acceleration $g$, and radius $R_{\rm P}$, (ii) the surface is not molten, (iii) tides deposit sufficient energy, and (iv) the planet has nonzero eccentricity. We calculate the approximate liquid core radius as a function of $μ/ (ρg R_{\rm P})$, and find that more than $90\%$ of the core will melt due to this runaway for $μ/ (ρg R_{\rm P})\gtrsim1$. From all currently confirmed exoplanets, we find that the terrestrial planets in the L98-59 system are the most promising candidates for sustaining active volcanism. However, uncertainties regarding the quality factors and the details of tidal heating and cooling mechanisms prohibit definitive claims of volcanism on any of these planets. We generate synthetic transmission spectra of these planets assuming Venus-like atmospheric compositions with an additional 5, 50, and $98\%$ SO$_2$ component, which is a tracer of volcanic activity. We find a $\gtrsim 3 σ$ preference for a model with SO$_2$ with 5-10 transits with $\textit{JWST}$ for L98-59bcd. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 16 pages, 8 Figures, accepted for publication in ApJ

Dynamical Friction Models for Black-Hole Binary Formation in AGN Disks

Authors: Kecheng Qian, Jiaru Li, Dong Lai

Abstract: Stellar-mass black holes (sBH) embedded in gaseous disks of active galactic nuclei (AGN) can be important sources of detectable gravitational radiation for LIGO/Virgo when they form binaries and coalesce due to orbital decay. In this paper, we study the effect of gas dynamical friction (DF) on the formation of BH binaries in AGN disks using $N$-body simulations. We employ two simplified models of… ▽ More Stellar-mass black holes (sBH) embedded in gaseous disks of active galactic nuclei (AGN) can be important sources of detectable gravitational radiation for LIGO/Virgo when they form binaries and coalesce due to orbital decay. In this paper, we study the effect of gas dynamical friction (DF) on the formation of BH binaries in AGN disks using $N$-body simulations. We employ two simplified models of DF, with the force on the BH depending on $Δ{\bf v}$, the velocity of the sBH relative to the background Keplerian gas. We integrate the motion of two sBH initially on circular orbits around the central supermassive black hole (SMBH), and evaluate the probability of binary formation under various conditions. We find that both models of DF (with different dependence of the frictional coefficient on $|Δ{\bf v}|$) can foster the formation of binaries when the effective friction timescale $τ$ satisfies $Ω_{\rm K}τ\lesssim 20-30$ (where $Ω_{\rm K}$ is the Keplerian frequency around the SMBH): prograde binaries are formed when the DF is stronger (smaller $τ$), while retrograde binaries dominate when the DF is weaker (larger $τ$). We determine the distribution of both prograde and retrograde binaries as a function of initial orbital separation and the DF strength. Using our models of DF, we show that for a given sBH number density in the AGN disk, the formation rate of sBH binaries increases with decreasing $τ$ and can reach a moderate value with a sufficiently strong DF. △ Less

Submitted 12 January, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: 17 pages, 13 figures, submitted to ApJ

Formation of Merging Stellar-Mass Black Hole Binaries by Gravitational Wave Emission in Active Galactic Nucleus Disks

Authors: Barak Rom, Re'em Sari, Dong Lai

Abstract: Many stellar-mass Black Holes (sBHs) are expected to orbit supermassive black holes at galactic centers. For galaxies with Active Galactic Nuclei (AGN), it is likely that the sBHs reside in a disk. We study the formation of sBH binaries via gravitational wave emission in such disks. We examine analytically the dynamics of two sBHs orbiting a supermassive black hole, estimate the capture cross sect… ▽ More Many stellar-mass Black Holes (sBHs) are expected to orbit supermassive black holes at galactic centers. For galaxies with Active Galactic Nuclei (AGN), it is likely that the sBHs reside in a disk. We study the formation of sBH binaries via gravitational wave emission in such disks. We examine analytically the dynamics of two sBHs orbiting a supermassive black hole, estimate the capture cross section, and derive the eccentricity distribution of bound binaries at different frequency bands. We find that the majority of the merging sBH binaries, assembled in this manner, can be measured as highly eccentric, detectable in the LIGO-Virgo-KAGRA (LVK) band from their formation, with $(1-e)\ll1$, through their circularization and up to their merger; the remaining binaries circularize to small eccentricities ($e\lesssim0.3$) before entering the LVK band. More eccentric mergers would be observed for sBHs with higher random velocities, closer to the supermassive black hole, or at lower observing frequency bands, as planned in future gravitational wave detectors such as the Einstein Telescope and LISA. △ Less

Submitted 14 March, 2024; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 11 pages, 6 figures

Planet-Driven Scatterings of Planetesimals Into a Star: Probability, Timescale and Applications

Authors: Laetitia Rodet, Dong Lai

Abstract: A planetary system can undergo multiple episodes of intense dynamical activities throughout its life, resulting in the production of star-grazing planetesimals (or exocomets) and pollution of the host star. Such activity is especially pronounced when giant planets interact with other small bodies during the system's evolution. However, due to the chaotic nature of the dynamics, it is difficult to… ▽ More A planetary system can undergo multiple episodes of intense dynamical activities throughout its life, resulting in the production of star-grazing planetesimals (or exocomets) and pollution of the host star. Such activity is especially pronounced when giant planets interact with other small bodies during the system's evolution. However, due to the chaotic nature of the dynamics, it is difficult to determine the properties of the perturbing planet(s) from the observed planetesimal-disruption activities. In this study, we examine the outcomes of planetesimal-planet scatterings in a general setting. We focus on one-planet systems, and determine the likelihood and timescale of planetesimal disruption by the host star as a function of the planet properties. We obtain a new analytical expression for the minimum distance a scattering body can reach, extending previous results by considering finite planet eccentricity and non-zero planetesimal mass. Through N-body simulations, we derive the distribution of minimum distances and the likelihood and timescales of three possible outcomes of planetesimal-planet scatterings: collision with the planet, ejection, and disruption by the star. For planetesimals with negligible mass, we identify four defining dimensionless parameters (the planet eccentricity, planet-to-star mass ratio, planet radius to semi-major axis ratio, and the stellar disruption radius to planet semi-major axis ratio) that enable us to scale the problem and generalize our findings to a wide range of orbital configurations. Using these results, we explore three applications: falling evaporating bodies in the beta Pictoris system, white dwarf pollution due to planetesimal disruption and planet engulfment by main-sequence stars. △ Less

Submitted 23 December, 2023; v1 submitted 20 August, 2023; originally announced August 2023.

Comments: Accepted to MNRAS (2023)

On the pollution of white dwarfs by exo-Oort cloud comets

Authors: Christopher E. O'Connor, Dong Lai, Darryl Z. Seligman

Abstract: A large fraction of white dwarfs (WDs) have metal-polluted atmospheres, which are produced by accreting material from remnant planetary systems. The composition of the accreted debris broadly resembles that of rocky Solar System objects. Volatile-enriched debris with compositions similar to long-period comets (LPCs) is rarely observed. We attempt to reconcile this dearth of volatiles with the prem… ▽ More A large fraction of white dwarfs (WDs) have metal-polluted atmospheres, which are produced by accreting material from remnant planetary systems. The composition of the accreted debris broadly resembles that of rocky Solar System objects. Volatile-enriched debris with compositions similar to long-period comets (LPCs) is rarely observed. We attempt to reconcile this dearth of volatiles with the premise that exo-Oort clouds (XOCs) occur around a large fraction of planet-hosting stars. We estimate the comet accretion rate from an XOC analytically, adapting the 'loss cone' theory of LPC delivery in the Solar System. We investigate the dynamical evolution of an XOC during late stellar evolution. Using numerical simulations, we show that 1 to 30 per cent of XOC objects remain bound after anisotropic stellar mass loss imparting a WD natal kick of $\sim$1 km/s. We also characterize the surviving comets' distribution function. Surviving planets orbiting a WD can prevent the accretion of XOC comets by the star. A planet's 'dynamical barrier' is effective at preventing comet accretion if the energy kick imparted by the planet exceeds the comet's orbital binding energy. By modifying the loss cone theory, we calculate the amount by which a planet reduces the WD's accretion rate. We suggest that the scarcity of volatile-enriched debris in polluted WDs is caused by an unseen population of 10-100 AU scale giant planets acting as barriers to incoming LPCs. Finally, we constrain the amount of volatiles delivered to a planet in the habitable zone of an old, cool WD. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: 18 pages, 12 figures; submitted to MNRAS. Comments welcome!

Giant planet engulfment by evolved giant stars: light curves, asteroseismology, and survivability

Authors: Christopher E. O'Connor, Lars Bildsten, Matteo Cantiello, Dong Lai

Abstract: About ten percent of Sun-like ($1$-$2 M_\odot$) stars will engulf a $1$-$10 M_{\rm J}$ planet as they expand during the red giant branch (RGB) or asymptotic giant branch (AGB) phase of their evolution. Once engulfed, these planets experience a strong drag force in the star's convective envelope and spiral inward, depositing energy and angular momentum. For these mass ratios, the inspiral takes… ▽ More About ten percent of Sun-like ($1$-$2 M_\odot$) stars will engulf a $1$-$10 M_{\rm J}$ planet as they expand during the red giant branch (RGB) or asymptotic giant branch (AGB) phase of their evolution. Once engulfed, these planets experience a strong drag force in the star's convective envelope and spiral inward, depositing energy and angular momentum. For these mass ratios, the inspiral takes $\sim 10$-$10^{2}$ years ($\sim 10^{2}$-$10^{3}$ orbits); the planet undergoes tidal disruption at a radius of $\sim R_\odot$. We use the Modules for Experiments in Stellar Astrophysics (MESA) software instrument to track the stellar response to the energy deposition while simultaneously evolving the planetary orbit. For RGB stars, as well as AGB stars with $M_{\rm p} \lesssim 5 M_{\rm J}$ planets, the star responds quasistatically but still brightens measurably on a timescale of years. In addition, asteroseismic indicators, such as the frequency spacing or rotational splitting, differ before and after engulfment. For AGB stars, engulfment of a $M_{\rm p} \gtrsim 5 M_{\rm J}$ planet drives supersonic expansion of the envelope, causing a bright, red, dusty eruption similar to a "luminous red nova." Based on the peak luminosity, color, duration, and expected rate of these events, we suggest that engulfment events on the AGB could be a significant fraction of low-luminosity red novae in the Galaxy. We do not find conditions where the envelope is ejected prior to the planet's tidal disruption, complicating the interpretation of short-period giant planets orbiting white dwarfs as survivors of common-envelope evolution. △ Less

Submitted 13 June, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: 24 pages, 11 figures, 1 table. Accepted to ApJ

Eccentricity Growth of Massive Planets inside Cavities of Protoplanetary Discs

Authors: Marina M. Romanova, Alexander V. Koldoba, Galina V. Ustyugova, Dong Lai, Richard V. E. Lovelace

Abstract: We carry out hydrodynamical simulations to study the eccentricity growth of a 1-30 Jupiter mass planet located inside the fixed cavity of a protoplanetary disc. The planet exchanges energy and angular momentum with the disc at resonant locations, and its eccentricity grows due to Lindblad resonances. We observe several phases of eccentricity growth where different eccentric Lindblad resonances dom… ▽ More We carry out hydrodynamical simulations to study the eccentricity growth of a 1-30 Jupiter mass planet located inside the fixed cavity of a protoplanetary disc. The planet exchanges energy and angular momentum with the disc at resonant locations, and its eccentricity grows due to Lindblad resonances. We observe several phases of eccentricity growth where different eccentric Lindblad resonances dominate from 1:3 up to 3:5. The maximum values of eccentricity reached in our simulations are 0.65-0.75. We calculate the eccentricity growth rate for different planet masses and disc parameters and derive analytical dependencies on these parameters. We observe that the growth rate is proportional to both the planet's mass and the characteristic disc mass for a wide range of parameters. In a separate set of simulations, we derived the width of the 1:3 Lindblad resonance. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 19 pages, 18 figures, accepted by MNRAS

Wide binary stars formed in the turbulent interstellar medium

Authors: Siyao Xu, Hsiang-Chih Hwang, Chris Hamilton, Dong Lai

Abstract: The ubiquitous interstellar turbulence regulates star formation and the scaling relations between the initial velocity differences and the initial separations of stars. We propose that the formation of wide binaries with initial separations $r$ in the range $\sim 10^3~\text{AU} \lesssim r \lesssim 10^5$ AU is a natural consequence of star formation in the turbulent interstellar medium. With the de… ▽ More The ubiquitous interstellar turbulence regulates star formation and the scaling relations between the initial velocity differences and the initial separations of stars. We propose that the formation of wide binaries with initial separations $r$ in the range $\sim 10^3~\text{AU} \lesssim r \lesssim 10^5$ AU is a natural consequence of star formation in the turbulent interstellar medium. With the decrease of $r$, the mean turbulent relative velocity $v_\text{tur}$ between a pair of stars decreases, while the largest velocity $v_\text{bon}$ at which they still may be gravitationally bound increases. When $v_\text{tur} < v_\text{bon}$, a wide binary can form. In this formation scenario, we derive the eccentricity distribution $p(e)$ of wide binaries for an arbitrary relative velocity distribution. By adopting a turbulent velocity distribution, we find that wide binaries at a given initial separation generally exhibit a superthermal $p(e)$. This provides a natural explanation for the observed superthermal $p(e)$ of the wide binaries in the Solar neighborhood. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: 7 pages, 4 figures, submitted to ApJL

Hydrodynamical Evolution of Black-Hole Binaries Embedded in AGN Discs: III. The Effects of Viscosity

Authors: Rixin Li, Dong Lai

Abstract: Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs offer a distinct dynamical channel to produce black hole mergers detected in gravitational waves by LIGO/Virgo. To understand their orbital evolution through interactions with the disc gas, we perform a suite of 2D high-resolution, local shearing box, viscous hydrodynamical simulations of equal-mass binaries. We… ▽ More Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs offer a distinct dynamical channel to produce black hole mergers detected in gravitational waves by LIGO/Virgo. To understand their orbital evolution through interactions with the disc gas, we perform a suite of 2D high-resolution, local shearing box, viscous hydrodynamical simulations of equal-mass binaries. We find that viscosity not only smooths the flow structure around prograde circular binaries, but also greatly raises their accretion rates. The torque associated with accretion may be overwhelmingly positive and dominate over the gravitational torque at a high accretion rate. However, the accreted angular momentum per unit mass decreases with increasing viscosity, making it easier to shrink the binary orbit. In addition, retrograde binaries still experience rapid orbital decay, and prograde eccentric binaries still experience eccentricity damping. Our numerical experiments further show that prograde binaries are more likely to be hardened if the physical sizes of the accretors are sufficiently small such that the accretion rate is reduced. The dependency of the binary accretion rate on the accretor size can be weaken through boosted accretion either due to a high viscosity or a more isothermal-like equation of state (EOS). Our results widen the explored parameter space for the hydrodynamics of embedded BBHs and demonstrate that their orbital evolution in AGN discs is a complex, multifaceted problem. △ Less

Submitted 12 January, 2024; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: 13 pages, 8 figures, 1st revision submitted to MNRAS. arXiv admin note: text overlap with arXiv:2207.01125

Hydrodynamical Simulations of Circumbinary Accretion: Balance between Heating and Cooling

Authors: Hai-Yang Wang, Xue-Ning Bai, Dong Lai, Douglas N. C. Lin

Abstract: Hydrodynamical interaction in circumbinary discs (CBDs) plays a crucial role in various astrophysical systems, ranging from young stellar binaries to supermassive black hole binaries in galactic centers. Most previous simulations of binary-disc systems have adopted locally isothermal equation of state. In this study, we use the grid-based code $\texttt{Athena++}$ to conduct a suite of two-dimensio… ▽ More Hydrodynamical interaction in circumbinary discs (CBDs) plays a crucial role in various astrophysical systems, ranging from young stellar binaries to supermassive black hole binaries in galactic centers. Most previous simulations of binary-disc systems have adopted locally isothermal equation of state. In this study, we use the grid-based code $\texttt{Athena++}$ to conduct a suite of two-dimensional viscous hydrodynamical simulations of circumbinary accretion on a cartesian grid, resolving the central cavity of the binary. The gas thermodynamics is treated by thermal relaxation towards an equilibrium temperature (based on the constant$-β$ cooling ansatz, where $β$ is the cooling time in units of the local Keplerian time). Focusing on equal mass, circular binaries in CBDs with (equilibrium) disc aspect ratio $H/R=0.1$, we find that the cooling of the disc gas significantly influences the binary orbital evolution, accretion variability, and CBD morphology, and the effect depends sensitively on the disc viscosity prescriptions. When adopting a constant kinematic viscosity, a finite cooling time ($β\gtrsim 0.1$) leads to binary inspiral as opposed to outspiral and the CBD cavity becomes more symmetric. When adopting a dynamically varying $α-$viscosity, binary inspiral only occurs within a narrow range of cooling time (corresponding to $β$ around 0.5). △ Less

Submitted 19 September, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 20 pages, 20 figures, Accepted for publication in MNRAS

On the Role of Gas Cooling in the Dynamics of Circumbinary Disks

Authors: Hai-Yang Wang, Xue-Ning Bai, Dong Lai

Abstract: Hydrodynamical interactions between binaries and circumbinary disks (CBDs) play an important role in a variety of astrophysical systems, from young stellar binaries to supermassive black hole binaries. Previous simulations of CBDs have mostly employed locally isothermal equation of state. We carry out two-dimensional viscous hydrodynamic simulations of CBDs around equal-mass, circular binaries, tr… ▽ More Hydrodynamical interactions between binaries and circumbinary disks (CBDs) play an important role in a variety of astrophysical systems, from young stellar binaries to supermassive black hole binaries. Previous simulations of CBDs have mostly employed locally isothermal equation of state. We carry out two-dimensional viscous hydrodynamic simulations of CBDs around equal-mass, circular binaries, treating the gas thermodynamics by thermal relaxation towards equilibrium temperature (the constant-$β$ cooling ansatz, where $β$ is the cooling time in units of the local Keplerian time). As an initial study, we use the grid-based code Athena++ on a polar grid, covering an extended disk outside the binary co-orbital region. We find that with a longer cooling time, the accretion variability is gradually suppressed, and the morphology of the CBD becomes more symmetric. The disk also shows evidence of hysterisis behavior depending on the initial conditions. Gas cooling also affects the rate of angular momentum transfer between the binary and the CBD, where given our adopted disk thickness and viscosity ($H/r\sim 0.1$ and $α\sim 0.1$), the binary orbit expands while undergoing accretion for most $β$ values between 0 and 4.0 except over a narrow range of intermediate $β$ values. The validity of using polar grid excising the central domain is also discussed. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: 14 pages, 12 figures, resubmitted to APJ

Hydrodynamical Simulations of Black-Hole Binary Formation in AGN Disks

Authors: Jiaru Li, Adam M. Dempsey, Hui Li, Dong Lai, Shengtai Li

Abstract: We study the close encounters between two single black holes (BHs) embedded in an AGN disk using a series of global 2D hydrodynamics simulations. We find that when the disk density is sufficiently high, bound BH binaries can be formed by the collision of their circum-single disks. Our analysis demonstrates that, after a BH pair passes the pericenter of their relative trajectory, a gas post-collisi… ▽ More We study the close encounters between two single black holes (BHs) embedded in an AGN disk using a series of global 2D hydrodynamics simulations. We find that when the disk density is sufficiently high, bound BH binaries can be formed by the collision of their circum-single disks. Our analysis demonstrates that, after a BH pair passes the pericenter of their relative trajectory, a gas post-collision drag may slow down the BHs, possibly forcing the two BHs to stay tightly bound. A binary formed by a close encounter can have a compact semi-major axis, large eccentricity, and retrograde orbital angular momentum. We provide a fitting formula that can accurately predict whether a close encounter can form a binary based on the gas mass and the incoming energy of the encounter. This fitting formula can be easily implemented in other long-term simulations that study the dynamical evolution of BHs in AGN disks. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 10 pages, 5 figures

High-Eccentricity Migration with Disk-Induced Spin-Orbit Misalignment: a Preference for Perpendicular Hot Jupiters

Authors: Michelle Vick, Yubo Su, Dong Lai

Abstract: High-eccentricity migration is a likely formation mechanism for many observed hot Jupiters, particularly those with a large misalignment between the stellar spin axis and orbital angular momentum axis of the planet. In one version of high-eccentricity migration, an inclined stellar companion excites von Zeipel-Lidov-Kozai (ZLK) eccentricity oscillations of a cold Jupiter, and tidal dissipation cau… ▽ More High-eccentricity migration is a likely formation mechanism for many observed hot Jupiters, particularly those with a large misalignment between the stellar spin axis and orbital angular momentum axis of the planet. In one version of high-eccentricity migration, an inclined stellar companion excites von Zeipel-Lidov-Kozai (ZLK) eccentricity oscillations of a cold Jupiter, and tidal dissipation causes the planet's orbit to shrink and circularize. Throughout this process, the stellar spin can evolve chaotically, resulting in highly misaligned hot Jupiters. Previous population studies of this migration mechanism have assumed that the stellar spin is aligned with the planetary orbital angular momentum when the companion begins to induce ZLK oscillations. However, in the presence of a binary companion, the star's obliquity may be significantly excited during the dissipation of its protoplanetary disk. We calculate the stellar obliquities produced in the protoplanetary disk phase and use these to perform an updated population synthesis of ZLK-driven high-eccentricity migration. We find that the resulting obliquity distribution of HJ systems is predominantly retrograde with a broad peak near 90$^\circ$. The distribution we obtain has intriguing similarities to the recently-observed preponderance of perpendicular planets close to their host stars. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 12 pages, 5 figures

Resonant Excitation of Planetary Eccentricity due to a Dispersing Eccentric Protoplanetary Disk: A New Mechanism of Generating Large Planetary Eccentricities

Authors: Jiaru Li, Dong Lai

Abstract: We present a new mechanism of generating large planetary eccentricities. This mechanism applies to planets within the inner cavities of their companion protoplanetary disks. A massive disk with an inner truncation may become eccentric due to non-adiabatic effects associated with gas cooling, and can retain its eccentricity in long-lived coherently-precessing eccentric modes; as the disk disperses,… ▽ More We present a new mechanism of generating large planetary eccentricities. This mechanism applies to planets within the inner cavities of their companion protoplanetary disks. A massive disk with an inner truncation may become eccentric due to non-adiabatic effects associated with gas cooling, and can retain its eccentricity in long-lived coherently-precessing eccentric modes; as the disk disperses, the inner planet will encounter a secular resonance with the eccentric disk when the planet and the disk have the same apsidal precession rates; the eccentricity of the planet is then excited to a large value as the system goes through the resonance. In this work, we solve the eccentric modes of a model disk for a wide range of masses. We then adopt an approximate secular dynamics model to calculate the long-term evolution of the "planet + dispersing disk" system. The planet attains a large eccentricity (between 0.1 and 0.6) in our calculations, even though the disk eccentricity is quite small ($\lesssim0.05$). This eccentricity excitation can be understood in terms of the mode conversion (``avoided crossing'') phenomenon associated with the evolution of the "planet + disk" eccentricity eigenstates. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 10 pages, 7 figures, submitted to ApJ

Circumbinary Accretion: From Binary Stars to Massive Binary Black Holes

Authors: Dong Lai, Diego J. Muñoz

Abstract: We review recent works on the dynamics of circumbinary accretion, including time variability, angular momentum transfer between the disk and the binary, and the secular evolution of accreting binaries. These dynamics can impact stellar binary formation/evolution, circumbinary planet formation/migration, and the evolution of (super)massive black-hole binaries. We discuss the dynamics and evolution… ▽ More We review recent works on the dynamics of circumbinary accretion, including time variability, angular momentum transfer between the disk and the binary, and the secular evolution of accreting binaries. These dynamics can impact stellar binary formation/evolution, circumbinary planet formation/migration, and the evolution of (super)massive black-hole binaries. We discuss the dynamics and evolution of inclined/warped circumbinary disks and connect with recent observations of protoplanetary disks. A special kind of circumbinary accretion involves binaries embedded in "big" disks, which may contribute to the mergers of stellar-mass black holes in AGN disks. Highlights include: $\bullet$ Circumbinary accretion is highly variable, being modulated at $P_{\rm b}$ (the binary period) or $\sim 5P_{\rm p}$, depending on the binary eccentricity $e_{\rm b}$ and mass ratio $q_{\rm b}$. $\bullet$ The inner region of the circumbinary disk can develop coherent eccentric structure, which may modulate the accretion and affect the physical processes (e.g. planet migration) taking place in the disk. $\bullet$ Over long timescales, circumbinary accretion steers binaries toward equal masses, and it does not always lead to binary orbital decay, as is commonly assumed. The secular orbital evolution depends on the binary parameters ($e_{\rm b}$ and $q_{\rm b}$), and on the thermodynamic properties of the accreting gas. $\bullet$ A misaligned disk around a low-eccentricity binary tends to evolve toward coplanarity due to viscous dissipation. But when $e_{\rm b}$ is significant, the disk can evolve toward "polar alignment", with the disk plane perpendicular to the binary plane. △ Less

Submitted 24 January, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Draft version of 2023 Annual Reviews article, comments are welcome. 39 pages including references, 16 figures

IXPE Detection of Polarized X-rays from Magnetars and Photon Mode Conversion at QED Vacuum Resonance

Authors: Dong Lai

Abstract: The recent observations of the anomalous X-ray pulsars 4U 0142+61 and 1RXS J170849.0-400910 by the Imaging X-ray Polarimetry Explorer (IXPE) opened up a new avenue to study magnetars, neutron stars endowed with superstrong magnetic fields ($B\gtrsim 10^{14}$~G). The detected polarized X-rays from 4U 0142+61 exhibit a 90$^\circ$ linear polarization swing from low photon energies ($E\lesssim 4$~keV)… ▽ More The recent observations of the anomalous X-ray pulsars 4U 0142+61 and 1RXS J170849.0-400910 by the Imaging X-ray Polarimetry Explorer (IXPE) opened up a new avenue to study magnetars, neutron stars endowed with superstrong magnetic fields ($B\gtrsim 10^{14}$~G). The detected polarized X-rays from 4U 0142+61 exhibit a 90$^\circ$ linear polarization swing from low photon energies ($E\lesssim 4$~keV) to high energies ($E\gtrsim 5.5$~keV). We show that this swing can be explained by photon polarization mode conversion at the vacuum resonance in the magnetar atmosphere; the resonance arises from the combined effects of plasma-induced birefringence and QED-induced vacuum birefringence in strong magnetic fields. This explanation suggests that the atmosphere of 4U 0142 be composed of partially ionized heavy elements, and the surface magnetic field be comparable or less than $10^{14}$~G, consistent with the dipole field inferred from the measured spindown. It also implies that the spin axis of 4U 0142+61 is aligned with its velocity direction. The polarized X-rays from 1RXS J170849.0-400910 do not show such $90^\circ$ swing, consistent with magnetar atmospheric emission with $B\gtrsim 5\times 10^{14}$~G. △ Less

Submitted 13 April, 2023; v1 submitted 27 September, 2022; originally announced September 2022.

Comments: 8 pages. Updated to include discussions on both magnetars detected by IXPE. To appear in PNAS in 4/2023

Eccentric debris belts reveal the dynamical history of the companion exoplanet

Authors: Laetitia Rodet, Dong Lai

Abstract: In recent years, a number of eccentric debris belts have been observed in extrasolar systems. The most common explanation for their shape is the presence of a nearby eccentric planetary companion. The gravitational perturbation from such a companion would induce periodic eccentricity variations on the planetesimals in the belt, with a range of precession frequencies. The overall expected shape is… ▽ More In recent years, a number of eccentric debris belts have been observed in extrasolar systems. The most common explanation for their shape is the presence of a nearby eccentric planetary companion. The gravitational perturbation from such a companion would induce periodic eccentricity variations on the planetesimals in the belt, with a range of precession frequencies. The overall expected shape is an eccentric belt with a finite minimum width. However, several observed eccentric debris discs have been found to exhibit a narrower width than the theoretical expectation. In this paper, we study two mechanisms that can produce this small width: (i) the protoplanetary disc can interact with the planet and/or the planetesimals, slowly driving the eccentricity of the former and damping the eccentricities of the latter; (ii) the companion planet could have gained its eccentricity stochastically, through planet-planet scatterings. We show that under appropriate conditions, both of these scenarios offer a plausible way to reduce the minimum width of an eccentric belt exterior to a perturbing planet. However, the effects of protoplanetary discs are diminished at large separations (a > 10 au) due to the scarcity of gas and the limited disc lifetime. These findings suggest that one can use the shape and width of debris discs to shed light on the evolution of extrasolar systems, constraining the protoplanetary disc properties and the prevalence of planet-planet scatterings. Further observations of debris-harbouring systems could confirm whether thin debris belts are a common occurrence, or the results of rare initial conditions or evolutionary processes. △ Less

Submitted 26 September, 2022; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted in MNRAS (2022)

Hydrodynamical Evolution of Black-Hole Binaries Embedded in AGN Discs: II. Dependence on Equation of State, Binary Mass, and Separation Scales

Authors: Rixin Li, Dong Lai

Abstract: Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs offer a promising dynamical channel to produce black-hole mergers that are detectable by LIGO/Virgo. Modeling the interactions between the disc gas and the embedded BBHs is crucial to understand their orbital evolution. Using a suite of 2D high-resolution simulations of prograde equal-mass circular binaries in l… ▽ More Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs offer a promising dynamical channel to produce black-hole mergers that are detectable by LIGO/Virgo. Modeling the interactions between the disc gas and the embedded BBHs is crucial to understand their orbital evolution. Using a suite of 2D high-resolution simulations of prograde equal-mass circular binaries in local disc models, we systematically study how their hydrodynamical evolution depends on the equation of state (EOS; including the $γ$-law and isothermal EOS) and on the binary mass and separation scales (relative to the supermassive BH mass and the Hill radius, respectively). We find that binaries accrete slower and contract in orbit if the EOS is far from isothermal such that the surrounding gas is diffuse, hot, and turbulent. The typical orbital decay rate is of the order of a few times the mass doubling rate. For a fixed EOS, the accretion flows are denser, hotter, and more turbulent around more massive or tighter binaries. The torque associated with accretion is often comparable to the gravitational torque, so both torques are essential in determining the long-term binary orbital evolution. We carry out additional simulations with non-accreting binaries and find that their orbital evolution can be stochastic and is sensitive to the gravitational softening length, and the secular orbital evolution can be very different from those of accreting binaries. Our results indicate that stellar-mass BBHs may be hardened efficiently under ideal conditions, namely less massive and wider binaries embedded in discs with a non-isothermal EOS. △ Less

Submitted 3 July, 2022; originally announced July 2022.

Comments: 14 pages, 9 figures, submitted to MNRAS

Apsidal Alignment and Anti-Alignment of Planets in Mean-Motion Resonance: Disk-Driven Migration and Eccentricity Driving

Authors: JT Laune, Laetitia Rodet, Dong Lai

Abstract: Planets migrating in their natal discs can be captured into mean-motion resonance (MMR), in which the planets' periods are related by integer ratios. Recent observations indicate that planets in MMR can be either apsidally aligned or anti-aligned. How these different configurations arise is unclear. In this paper, we study the MMR capture process of migrating planets, focusing on the property of t… ▽ More Planets migrating in their natal discs can be captured into mean-motion resonance (MMR), in which the planets' periods are related by integer ratios. Recent observations indicate that planets in MMR can be either apsidally aligned or anti-aligned. How these different configurations arise is unclear. In this paper, we study the MMR capture process of migrating planets, focusing on the property of the apsidal angles of the captured planets. We show that the standard picture of MMR capture, in which the planets undergo convergent migration and experience eccentricity damping due to planet-disc interactions, always leads to apsidal anti-alignment of the captured planets. However, when the planets experience eccentricity driving from the disc, apsidally aligned configuration in MMR can be produced. In this configuration, both planets' resonance angles circulate, but a "mixed" resonance angle librates and traps the planets near the nominal resonance location. The MMR capture process in the presence of disc eccentricity driving is generally complex and irregular, and can lead to various outcomes, including apsidal alignment and anti-alignment, as well as the disruption of the resonance. We suggest that the two resonant planets in the K2-19 system, with their moderate eccentricities and aligned apsides, have experienced eccentricity driving from their natal disc in the past. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 19 pages, 16 figures

Dynamical Instability in Multi-Orbiter Systems with Gas Friction

Authors: Jiaru Li, Laetitia Rodet, Dong Lai

Abstract: Closely-packed multi-planet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in AGN discs around supermassive black holes. In this paper, we use $N$-body integrations to evaluate… ▽ More Closely-packed multi-planet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in AGN discs around supermassive black holes. In this paper, we use $N$-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of planetary spacing and planet-to-star mass ratios. We find that the stability of a system depends on the damping timescale $τ$ relative to the zero-friction instability growth timescale $t_{\rm inst}$. In a two-planet system, the frictional damping can stabilise the dynamical evolution if $t_{\rm inst}\gtrsimτ$. With three planets, $t_{\rm inst} \gtrsim 10τ- 100τ$ is needed for stabilisation. When the separations between the planetary orbits are sufficiently small, $t_{\rm inst}$ can be less than the synodic period between the planets, which makes frictional stabilisation unlikely to occur. As the orbital spacing increases, the instability timescale tends to grow exponentially on average, but it can vary by a few orders of magnitude depending on the initial orbital phases of the planets. In general, the stable region (at large orbital spacings) and unstable region (at small orbital spacings) are separated by a transition zone, in which the (in)stability of the system is not guaranteed. We also devise a linear map to analyse the dynamical instability of the "planet + test-mass" system, and we find qualitatively similar results to the $N$-body simulations. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 13 pages, 15 figures, 7 tables, submitted to MNRAS

Long-term Evolution of Tightly-Packed Stellar Black Holes in AGN Disks: Formation of Merging Black-Hole Binaries via Close Encounters

Authors: Jiaru Li, Dong Lai, Laetitia Rodet

Abstract: We study the long-term evolution of two or more stellar black holes (BHs) on initially separated but unstable circular orbits around a supermassive BH (SMBH). Such a close-packed orbital configuration can naturally arise from BH migrations in the AGN disk. Dynamical instability of the orbits leads to recurring close encounters between two BHs, during which the BH separation $r_{\rm p}$ becomes les… ▽ More We study the long-term evolution of two or more stellar black holes (BHs) on initially separated but unstable circular orbits around a supermassive BH (SMBH). Such a close-packed orbital configuration can naturally arise from BH migrations in the AGN disk. Dynamical instability of the orbits leads to recurring close encounters between two BHs, during which the BH separation $r_{\rm p}$ becomes less than the Hill radius $R_{\rm H}$. In the rare very close encounters (with $r_{\rm p}$ several orders of magnitude less than $R_{\rm H}$), a tight merging BH binary can form with the help of gravitational wave emission. We use $N$-body simulations to study the time evolution of close encounters of various degrees of "closeness" and the property of the resulting binary BH mergers. For a typical "SMBH + 2 BHs" system, the averaged cumulative number of close encounters (with $r_{\rm p} \lesssim R_{\rm H}$) scales approximately as $\propto t^{0.5}$. The minimum encounter separation $r_{\rm p}$ follows a linear cumulative distribution $P(<r_{\rm p}) \propto r_{\rm p}$ for $r_{\rm p} \ll R_{\rm H}$. From these, we obtain a semi-analytical expression for the averaged rate of binary captures that lead to BH mergers. Our results suggest that close-packed BHs in AGN disks may take a long time ($\gtrsim 10^7$ orbits around the SMBH) to experience a sufficiently close encounter and form a bound binary, although this time can be shorter if the initial BH orbits are highly aligned. The BH binary mergers produced in this scenario always have high eccentricities when entering the LIGO band, and have a broad distribution of orbital inclinations relative to the original AGN disk. We also explore the effects of the gas disk and find that simple gas drags on the BHs do not necessarily lead to an enhanced BH binary capture rate. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: 15 pages, 17 figures, submitted to ApJ

Hydrodynamical Evolution of Black-Hole Binaries Embedded in AGN Discs

Authors: Rixin Li, Dong Lai

Abstract: Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs are possible progenitors of black-hole mergers detected in gravitational waves by LIGO/VIRGO. To better understand the hydrodynamical evolution of BBHs interacting with the disc gas, we perform a suite of high-resolution 2D simulations of binaries in local disc (shearing-box) models, considering various binary m… ▽ More Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs are possible progenitors of black-hole mergers detected in gravitational waves by LIGO/VIRGO. To better understand the hydrodynamical evolution of BBHs interacting with the disc gas, we perform a suite of high-resolution 2D simulations of binaries in local disc (shearing-box) models, considering various binary mass ratios, eccentricities and background disc properties. We use the $γ$-law equation of state and adopt a robust post-processing treatment to evaluate the mass accretion rate, torque and energy transfer rate on the binary to determine its long-term orbital evolution. We find that circular comparable-mass binaries contract, with an orbital decay rate of a few times the mass doubling rate. Eccentric binaries always experience eccentricity damping. Prograde binaries with higher eccentricities or smaller mass ratios generally have slower orbital decay rates, with some extreme cases exhibiting orbital expansion. The averaged binary mass accretion rate depends on the physical size of the accretor. The accretion flows are highly variable, and the dominant variability frequency is the apparent binary orbital frequency (in the rotating frame around the central massive BH) for circular binaries but gradually shifts to the radial epicyclic frequency as the binary eccentricity increases. Our findings demonstrate that the dynamics of BBHs embedded in AGN discs is quite different from that of isolated binaries in their own circumbinary discs. Furthermore, our results suggest that the hardening timescales of the binaries are much shorter than their migration timescales in the disc, for all reasonable binary and disc parameters. △ Less

Submitted 7 September, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 24 pages, 22 figures, accepted for publication in MNRAS

Evolution of Stellar Orbits Around Merging Massive Black-Hole Binary

Authors: Bin Liu, Dong Lai

Abstract: We study the long-term orbital evolution of stars around a merging massive or supermassive black-hole (BH) binary, taking into account the general relativistic effect induced by the BH spin. When the BH spin is significant compared to and misaligned with the binary orbital angular momentum, the orbital axis ($\hat{\mathbf{l}}$) of the circumbinary star can undergo significant evolution during the… ▽ More We study the long-term orbital evolution of stars around a merging massive or supermassive black-hole (BH) binary, taking into account the general relativistic effect induced by the BH spin. When the BH spin is significant compared to and misaligned with the binary orbital angular momentum, the orbital axis ($\hat{\mathbf{l}}$) of the circumbinary star can undergo significant evolution during the binary orbital decay driven by gravitational radiation. Including the spin effect of the primary (more massive) BH, we find that starting from nearly coplanar orbital orientations, the orbital axes $\hat{\mathbf{l}}$ of circumbinary stars preferentially evolve towards the spin direction after the merger of the BH binary, regardless of the initial BH spin orientation. Such alignment phenomenon, i.e., small final misalignment angle between $\hat{\mathbf{l}}$ and the spin axis of the remanent BH $\hat{\mathbf{S}}$, can be understood analytically using the principle of adiabatic invariance. For the BH binaries with extremely mass ratio ($m_2/m_1\lesssim0.01$), $\hat{\mathbf{l}}$ may experience more complicated evolution as adiabatic invariance breaks down, but the trend of alignment still works reasonably well when the initial binary spin-orbit angle is relatively small. Our result suggests that the correlation between the orientations of stellar orbits and the spin axis of the central BH could provide a potential signature of the merger history of the massive BH. △ Less

Submitted 28 April, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 13 pages, 15 figures, accepted for publication in MNRAS

Non-Trivial Oblique Spin Equilibria of Super-Earths in Multi-planetary Systems

Authors: Yubo Su, Dong Lai

Abstract: Many Sun-like stars are observed to host close-in super-Earths (SEs) as part of a multi-planetary system. In such a system, the spin of the SE evolves due to spin-orbit resonances and tidal dissipation. In the absence of tides, the planet's obliquity can evolve chaotically to large values. However, for close-in SEs, tidal dissipation is significant and suppresses the chaos, instead driving the spi… ▽ More Many Sun-like stars are observed to host close-in super-Earths (SEs) as part of a multi-planetary system. In such a system, the spin of the SE evolves due to spin-orbit resonances and tidal dissipation. In the absence of tides, the planet's obliquity can evolve chaotically to large values. However, for close-in SEs, tidal dissipation is significant and suppresses the chaos, instead driving the spin into various steady states. We find that the attracting steady states of the SE's spin are more numerous than previously thought, due to the discovery of a new class of "mixed-mode" high-obliquity equilibria. These new equilibria arise due to subharmonic responses of the parametrically-driven planetary spin, an unusual phenomenon that arises in nonlinear systems. Many SEs should therefore have significant obliquities, with potentially large impacts on the physical conditions of their surfaces and atmospheres. △ Less

Submitted 17 October, 2022; v1 submitted 5 December, 2021; originally announced December 2021.

Comments: 18 pages, 18 figures, accepted to MNRAS

Secular chaos in white-dwarf planetary systems: origins of metal pollution and short-period planetary companions

Authors: Christopher E. O'Connor, Jean Teyssandier, Dong Lai

Abstract: Secular oscillations in multi-planet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This "secular chaos" readily pushes the inner body to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and test particles using the ring-averaging… ▽ More Secular oscillations in multi-planet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This "secular chaos" readily pushes the inner body to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and test particles using the ring-averaging method, with emphasis on the relationship between the planets' properties and the time-scale and efficiency of chaotic diffusion. We find that secular chaos can excite extreme eccentricities on time-scales spanning several orders of magnitude in a given system. We apply our results to the evolution of planetary systems around white dwarfs (WDs), specifically the tidal disruption and high-eccentricity migration of planetesimals and planets. We find that secular chaos in a planetesimal belt driven by large ($\gtrsim 10 M_{\oplus}$), distant ($\gtrsim 10 \, {\rm au}$) planets can sustain metal accretion onto a WD over Gyr time-scales. We constrain the total mass of planetesimals initially present within the chaotic zone by requiring that the predicted mass delivery rate to the Roche limit be consistent with the observed metal accretion rates of WDs with atmospheric pollution throughout the cooling sequence. Based on the occurrence of long-period exoplanets and exo-asteroid belts, we conclude that secular chaos can be a significant (perhaps dominant) channel for polluting solitary WDs. Secular chaos can also produce short-period planets and planetesimals around WDs in concert with various circularization mechanisms. We discuss prospects for detecting exoplanets driving secular chaos around WDs using direct imaging and microlensing. △ Less

Submitted 26 April, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 21 pages, 10 figures, 2 tables; accepted in MNRAS

Dynamical tidal Love numbers of rapidly rotating planets and stars

Authors: Janosz W. Dewberry, Dong Lai

Abstract: Tidal interactions play an important role in many astrophysical systems, but uncertainties regarding the tides of rapidly rotating, centrifugally distorted stars and gaseous planets remain. We have developed a precise method for computing the dynamical, non-dissipative tidal response of rotating planets and stars, based on summation over contributions from normal modes driven by the tidal potentia… ▽ More Tidal interactions play an important role in many astrophysical systems, but uncertainties regarding the tides of rapidly rotating, centrifugally distorted stars and gaseous planets remain. We have developed a precise method for computing the dynamical, non-dissipative tidal response of rotating planets and stars, based on summation over contributions from normal modes driven by the tidal potential. We calculate the normal modes of isentropic polytropes rotating at up to $\simeq90\%$ of their critical breakup rotation rates, and tabulate fits to mode frequencies and tidal overlap coefficients that can be used to compute the frequency-dependent, non-dissipative tidal response (via potential Love numbers $k_{\ell m}$). Although fundamental modes (f-modes) possess dominant tidal overlap coefficients at (nearly) all rotation rates, we find that the strong coupling of retrograde inertial modes (i-modes) to tesseral ($\ell>|m|$) components of the tidal potential produces resonances that may be relevant to gas giants like Jupiter and Saturn. The coupling of f-modes in rapid rotators to multiple components of both the driving tidal potential and the induced gravitational field also affect the tesseral response, leading to significant deviations from treatments of rotation that neglect centrifugal distortion and high-order corrections. For very rapid rotation rates ($\gtrsim 70\%$ of breakup), mixing between prograde f-modes and i-modes significantly enhances the sectoral ($\ell=|m|$) tidal overlap of the latter. The tidal response of very rapidly rotating, centrifugally distorted planets or stars can also be modified by resonant sectoral f-modes that are secularly unstable via the Chandrasekhar-Friedman-Schutz (CFS) mechanism. △ Less

Submitted 1 December, 2021; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 19 pages, 9 figures, 4 tables, accepted for publication in ApJ

Journal ref: ApJ. 925 (2022) 124

Dynamical Tides in Eccentric Binaries Containing Massive Main-Sequence Stars: Analytical Expressions

Authors: Yubo Su, Dong Lai

Abstract: Tidal evolution of eccentric binary systems containing at least one massive main-sequence (MS) star plays an important role in the formation scenarios of merging compact-object binaries. The dominant dissipation mechanism in such systems involves tidal excitation of outgoing internal gravity waves at the convective-radiative boundary and dissipation of the waves at the stellar envelope/surface. We… ▽ More Tidal evolution of eccentric binary systems containing at least one massive main-sequence (MS) star plays an important role in the formation scenarios of merging compact-object binaries. The dominant dissipation mechanism in such systems involves tidal excitation of outgoing internal gravity waves at the convective-radiative boundary and dissipation of the waves at the stellar envelope/surface. We have derived analytical expressions for the tidal torque and tidal energy transfer rate in such binaries for arbitrary orbital eccentricities and stellar rotation rates. These expressions can be used to study the spin and orbital evolution of eccentric binaries containing massive MS stars, such as the progenitors of merging neutron star binaries. Applying our results to the PSR J0045-7319 system, which has a massive B-star companion and an observed, rapidly decaying orbit, we find that for the standard radius of convective core based on non-rotating stellar models, the B-star must have a significant retrograde and differential rotation in order to explain the observed orbital decay rate. Alternatively, we suggest that the convective core may be larger as a result of rapid stellar rotation and/or mass transfer to the B-star in the recent past during the post-MS evolution of the pulsar progenitor. △ Less

Submitted 22 October, 2021; originally announced October 2021.

Comments: 10 pages, 8 figures, submitted to MNRAS

Dynamics of Colombo's Top: Tidal Dissipation and Resonance Capture, With Applications to Oblique Super-Earths, Ultra-Short-Period Planets and Inspiraling Hot Jupiters

Authors: Yubo Su, Dong Lai

Abstract: We present a comprehensive theoretical study on the spin evolution of a planet under the combined effects of tidal dissipation and gravitational perturbation from an external companion. Such a "spin + companion" system (called Colombo's top) appears in many [exo]planetary contexts. The competition between the tidal torque (which drives spin-orbit alignment and synchronization) and the gravitationa… ▽ More We present a comprehensive theoretical study on the spin evolution of a planet under the combined effects of tidal dissipation and gravitational perturbation from an external companion. Such a "spin + companion" system (called Colombo's top) appears in many [exo]planetary contexts. The competition between the tidal torque (which drives spin-orbit alignment and synchronization) and the gravitational torque from the companion (which drives orbital precession of the planet) gives rise to two possible spin equilibria ("Tidal Cassini Equilibria", tCE) that are stable and attracting: the "simple" tCE1, which typically has a low spin obliquity, and the "resonant" tCE2, which can have a significant obliquity. The latter arises from a spin-orbit resonance and can be broken when the tidal alignment torque is stronger than the precessional torque from the companion. We characterize the long-term evolution of the planetary spin (both magnitude and obliquity) for an arbitrary initial spin orientation, and develop a new theoretical method to analytically obtain the probability of resonance capture driven by tidal dissipation. Applying our general theoretical results to exoplanetary systems, we find that a super-Earth (SE) with an exterior companion can have a substantial probability of being trapped in the high-obliquity tCE2, assuming that SEs have a wide range of primordial obliquities. We also evaluate the recently proposed "obliquity tide" scenarios for the formation of ultra-short-period Earth-mass planets and for the orbital decay of hot Jupiter WASP-12b. We find in both cases that the probability of resonant capture into tCE2 is generally low and that such a high-obliquity state can be easily broken by the required orbital decay. △ Less

Submitted 2 November, 2021; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: 21 pages, 24 figures, accepted to MNRAS

The Impact of Stellar Clustering on the Observed Multiplicity of Super-Earth systems: Outside-in Cascade of Orbital Misalignments Initiated by Stellar Flybys

Authors: Laetitia Rodet, Dong Lai

Abstract: A recent study suggests that the observed multiplicity of super-Earth (SE) systems is correlated with stellar overdensities: field stars in high phase-space density environments have an excess of single-planet systems compared to stars in low density fields. This correlation is puzzling as stellar clustering is expected to influence mostly the outer part of planetary systems. Here we examine the p… ▽ More A recent study suggests that the observed multiplicity of super-Earth (SE) systems is correlated with stellar overdensities: field stars in high phase-space density environments have an excess of single-planet systems compared to stars in low density fields. This correlation is puzzling as stellar clustering is expected to influence mostly the outer part of planetary systems. Here we examine the possibility that stellar flybys indirectly excite the mutual inclinations of initially coplanar SEs, breaking their co-transiting geometry. We propose that flybys excite the inclinations of exterior substellar companions, which then propagate the perturbation to the inner SEs. Using analytical calculations of the secular coupling between SEs and companions, together with numerical simulations of stellar encounters, we estimate the expected number of "effective" flybys per planetary system that lead to the destruction of the SE co-transiting geometry. Our analytical results can be rescaled easily for various SE and companion properties (masses and semi-major axes) and stellar cluster parameters (density, velocity dispersion and lifetime). We show that for a given SE system, there exists an optimal companion architecture that leads to the maximum number of effective flybys; this results from the trade-off between the flyby cross section and the companion's impact on the inner system. Subject to uncertainties in the cluster parameters, we conclude that this mechanism is inefficient if the SE system has a single exterior companion, but may play an important role in "SE + two companions" systems that were born in dense stellar clusters. Whether this effect causes the observed correlation between planet multiplicity and stellar overdensities remains to be confirmed. △ Less

Submitted 9 November, 2021; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: 16 pages, 18 figures, accepted in MNRAS

Constraining Saturn's interior with ring seismology: effects of differential rotation and stable stratification

Authors: Janosz W. Dewberry, Christopher R. Mankovich, Jim Fuller, Dong Lai, Wenrui Xu

Abstract: Normal mode oscillations in Saturn excite density and bending waves in the C Ring, providing a valuable window into the planet's interior. Saturn's fundamental modes (f modes) excite the majority of the observed waves, while gravito-inertial modes (rotationally modified g modes) associated with stable stratification in the deep interior provide a compelling explanation for additional density waves… ▽ More Normal mode oscillations in Saturn excite density and bending waves in the C Ring, providing a valuable window into the planet's interior. Saturn's fundamental modes (f modes) excite the majority of the observed waves, while gravito-inertial modes (rotationally modified g modes) associated with stable stratification in the deep interior provide a compelling explanation for additional density waves with low azimuthal wavenumbers m. However, multiplets of density waves with nearly degenerate frequencies, including an m=3 triplet, still lack a definitive explanation. We investigate the effects of rapid and differential rotation on Saturn's oscillations, calculating normal modes for independently constrained interior models. We use a non-perturbative treatment of rotation that captures the full effects of the Coriolis and centrifugal forces, and consequently the mixing of sectoral f modes with g modes characterized by very different spherical harmonic degrees. Realistic profiles for differential rotation associated with Saturn's zonal winds can enhance these mode interactions, producing detectable oscillations with frequencies separated by less than 1%. Our calculations demonstrate that a three-mode interaction involving an f mode and two g modes can feasibly explain the finely split m=3 triplet, although the fine-tuning required to produce such an interaction generally worsens agreement with seismological constraints provided by m=2 density waves. Our calculations additionally demonstrate that sectoral f mode frequencies are measurably sensitive to differential rotation in Saturn's convective envelope. Finally, we find that including realistic equatorial antisymmetry in Saturn's differential rotation profile couples modes with even and odd equatorial parity, producing oscillations that could in principle excite both density and bending waves simultaneously. △ Less

Submitted 16 December, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 23 pages, 16 figures, published in PSJ. Replacement fixes typos identified while in press

Journal ref: Planet. Sci. J. 2 (2021) 198

Probing the Spins of Supermassive Black Holes with Gravitational Waves from Surrounding Compact Binaries

Authors: Bin Liu, Dong Lai

Abstract: Merging compact black-hole (BH) binaries are likely to exist in the nuclear star clusters around supermassive BHs (SMBHs), such as Sgr A$^\ast$. They may also form in the accretion disks of active galactic nuclei. Such compact binaries can emit gravitational waves (GWs) in the low-frequency band (0.001-1 Hz) that are detectable by several planned space-borne GW observatories. We show that the orbi… ▽ More Merging compact black-hole (BH) binaries are likely to exist in the nuclear star clusters around supermassive BHs (SMBHs), such as Sgr A$^\ast$. They may also form in the accretion disks of active galactic nuclei. Such compact binaries can emit gravitational waves (GWs) in the low-frequency band (0.001-1 Hz) that are detectable by several planned space-borne GW observatories. We show that the orbital axis of the compact binary may experience significant variation due to the frame-dragging effect associated with the spin of the SMBH. The dynamical behavior of the orbital axis can be understood analytically as a resonance phenomenon. We show that rate of change of the binary orbital axis encodes the information on the spin of the SMBH. Therefore detecting GWs from compact binaries around SMBHs, particularly the modulation of the waveform associated with the variation of the binary orbital axis, can provide a new probe on the spins of SMBHs. △ Less

Submitted 24 January, 2022; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: Published in ApJ

Spin Dynamics of Extrasolar Giant Planets in Planet-Planet Scattering

Authors: Yu-Cian Hong, Dong Lai, Jonathan I. Lunine, Philip D. Nicholson

Abstract: Planet-planet scattering best explains the eccentricity distribution of extrasolar giant planets. Past literature showed that the orbits of planets evolve due to planet-planet scattering. This work studies the spin evolution of planets in planet-planet scattering in 2-planet systems. Spin can evolve dramatically due to spin-orbit coupling made possible by the evolving spin and orbital precession d… ▽ More Planet-planet scattering best explains the eccentricity distribution of extrasolar giant planets. Past literature showed that the orbits of planets evolve due to planet-planet scattering. This work studies the spin evolution of planets in planet-planet scattering in 2-planet systems. Spin can evolve dramatically due to spin-orbit coupling made possible by the evolving spin and orbital precession during the planet-planet scattering phase. The main source of torque to planet spin is the stellar torque, and the total planet-plane torque contribution is negligible. As a consequence of the evolution of the spin, planets can end up with significant obliquity (the angle between a planet's own orbit normal and spin axis) like planets in our Solar System. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 7 pages, 9 figures

Jupiter's Dynamical Love Number

Authors: Dong Lai

Abstract: Recent observations by the {\it Juno} spacecraft have revealed that the tidal Love number $k_2$ of Jupiter is $4\%$ lower than the hydrostatic value. We present a simple calculation of the dynamical Love number of Jupiter that explains the observed "anomaly". The Love number is usually dominated by the response of the (rotation-modified) f-modes of the planet. Our method also allows for efficient… ▽ More Recent observations by the {\it Juno} spacecraft have revealed that the tidal Love number $k_2$ of Jupiter is $4\%$ lower than the hydrostatic value. We present a simple calculation of the dynamical Love number of Jupiter that explains the observed "anomaly". The Love number is usually dominated by the response of the (rotation-modified) f-modes of the planet. Our method also allows for efficient computation of high-order dynamical Love numbers. While the inertial-mode contributions to the Love numbers are negligible, a sufficiently strong stratification in a large region of the planet's interior would induce significant g-mode responses and influence the measured Love numbers. △ Less

Submitted 14 March, 2021; v1 submitted 10 March, 2021; originally announced March 2021.

Comments: 6 pages (updated Figs.3-4 from V1)

The Mass Ratio Distribution of Tertiary Induced Binary Black Hole Mergers

Authors: Yubo Su, Bin Liu, Dong Lai

Abstract: Many proposed scenarios for black hole (BH) mergers involve a tertiary companion that induces von Zeipel-Lidov-Kozai (ZLK) eccentricity cycles in the inner binary. An attractive feature of such mechanisms is the enhanced merger probability when the octupole-order effects, also known as the eccentric Kozai mechanism, are important. This can be the case when the tertiary is of comparable mass to the… ▽ More Many proposed scenarios for black hole (BH) mergers involve a tertiary companion that induces von Zeipel-Lidov-Kozai (ZLK) eccentricity cycles in the inner binary. An attractive feature of such mechanisms is the enhanced merger probability when the octupole-order effects, also known as the eccentric Kozai mechanism, are important. This can be the case when the tertiary is of comparable mass to the binary components. Since the octupole strength [$\propto (1-q)/(1+q)$] increases with decreasing binary mass ratio $q$, such ZLK-induced mergers favor binaries with smaller mass ratios. We use a combination of numerical and analytical approaches to fully characterize the octupole-enhanced binary BH mergers and provide analytical criteria for efficiently calculating the strength of this enhancement. We show that for hierarchical triples with semi-major axis ratio $a/a_{\rm out}\gtrsim 0.01$-$0.02$, the binary merger fraction can increase by a large factor (up to $\sim 20$) as $q$ decreases from unity to $0.2$. The resulting mass ratio distribution for merging binary BHs produced in this scenario is in tension with the observed distribution obtained by the LIGO/VIRGO collaboration, although significant uncertainties remain about the initial distribution of binary BH masses and mass ratios. △ Less

Submitted 14 June, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

Comments: 18 pages, 19 figures, accepted to MNRAS

On the Correlation between Hot Jupiters and Stellar Clustering: High-eccentricity Migration Induced by Stellar Flybys

Authors: Laetitia Rodet, Yubo Su, Dong Lai

Abstract: A recent observational study suggests that the occurrence of hot Jupiters (HJs) around solar-type stars is correlated with stellar clustering. We study a new scenario for HJ formation, called "Flyby Induced High-e Migration", that may help explain this correlation. In this scenario, stellar flybys excite the eccentricity and inclination of an outer companion (giant planet, brown dwarf, or low-mass… ▽ More A recent observational study suggests that the occurrence of hot Jupiters (HJs) around solar-type stars is correlated with stellar clustering. We study a new scenario for HJ formation, called "Flyby Induced High-e Migration", that may help explain this correlation. In this scenario, stellar flybys excite the eccentricity and inclination of an outer companion (giant planet, brown dwarf, or low-mass star) at large distance (10-300 au), which then triggers high-e migration of an inner cold Jupiter (at a few astronomical units) through the combined effects of von Zeipel-Lidov-Kozai (ZLK) eccentricity oscillation and tidal dissipation. Using semianalytical calculations of the effective ZLK inclination window, together with numerical simulations of stellar flybys, we obtain the analytic estimate for the HJ occurrence rate in this formation scenario. We find that this "flyby induced high-e migration" could account for a significant fraction of the observed HJ population, although the result depends on several uncertain parameters, including the density and lifetime of birth stellar clusters, and the occurrence rate of the "cold Jupiter + outer companion" systems. △ Less

Submitted 10 May, 2021; v1 submitted 15 February, 2021; originally announced February 2021.

Comments: 10 pages, 5 figures, accepted in ApJ

3D MHD Simulations of Accretion onto Stars with Tilted Magnetic and Rotational Axes

Authors: M. M. Romanova, A. V. Koldoba, G. V. Ustyugova, A. A. Blinova, D. Lai, R. V. E. Lovelace

Abstract: We present results of global three-dimensional (3D) magnetohydrodynamic (MHD) simulations of accretion onto magnetized stars where both the magnetic and rotational axes of the star are tilted about the rotational axis of the disc. We observed that initially the inner parts of the disc are warped, tilted, and process due to the magnetic interaction between the magnetosphere and the disc. Later, lar… ▽ More We present results of global three-dimensional (3D) magnetohydrodynamic (MHD) simulations of accretion onto magnetized stars where both the magnetic and rotational axes of the star are tilted about the rotational axis of the disc. We observed that initially the inner parts of the disc are warped, tilted, and process due to the magnetic interaction between the magnetosphere and the disc. Later, larger tilted discs form with the size increasing with the magnetic moment of the star. The normal vector to the discs are tilted at different angles, from 5-10 degrees up to 30-40 degrees. Small tilts may result from the winding of the magnetic field lines about the rotational axis of the star and the action of the magnetic force which tends to align the disc. Another possible explanation is the magnetic Bardeen-Petterson effect in which the disc settles in the equatorial plane of the star due to precessional and viscous torques in the disc. Tilted discs slowly precess with the time scale of the order of 50 Keplerian periods at the reference radius (approx. 3 stellar radii). Our results can be applied to different types of stars where signs of tilted discs and/or slow precession have been observed. △ Less

Submitted 4 July, 2021; v1 submitted 19 December, 2020; originally announced December 2020.

Comments: 14 pages, 13 figures, MNRAS, accepted. The new version of the paper is more complete compared with the earlier version

Excitation of Spin-Orbit Misalignments in Stellar Binaries with Circumbinary Disks: Application to DI Herculis

Authors: Kassandra R. Anderson, Dong Lai

Abstract: The large spin-orbit misalignments in the DI Herculis stellar binary system have resolved the decades-long puzzle of the anomalously slow apsidal precession rate, but raise new questions regarding the origin of the obliquities. This paper investigates obliquity evolution in stellar binaries hosting modestly-inclined circumbinary disks. As the disk and binary axes undergo mutual precession, each ob… ▽ More The large spin-orbit misalignments in the DI Herculis stellar binary system have resolved the decades-long puzzle of the anomalously slow apsidal precession rate, but raise new questions regarding the origin of the obliquities. This paper investigates obliquity evolution in stellar binaries hosting modestly-inclined circumbinary disks. As the disk and binary axes undergo mutual precession, each oblate star experiences a torque from its companion star, so that the spin and orbital axes undergo mutual precession. As the disk loses mass through a combination of winds and accretion, the system may be captured into a high-obliquity Cassini state (a spin-orbit resonance). The final obliquity depends on the details of the disk dispersal. We construct a simple disk model to emulate disk dispersal due to viscous accretion and photoevaporation, and identify the necessary disk properties for producing the observed obliquities in DI Herculis. The disk must be massive (at least $10 \%$ of the binary mass). If accretion onto the binary is suppressed, the observed high stellar obliquities are reproduced with a binary-disk inclination of $\sim 5^\circ - 10^\circ$, but if substantial accretion occurs, the inclination must be larger, $\sim 20^\circ - 30^\circ$. If moderate accretion occurs, initially the disk must lose its mass slowly, but eventually lose its remaining mass abruptly, analogous to the observed two-timescale behavior for disks around T-Tauri stars. The spin feedback on the binary orbit causes the binary-disk inclination to decay as the obliquity evolves, a feature that is absent from the standard Cassini state treatment. △ Less

Submitted 24 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 13 pages, 6 figures, accepted for publication in ApJ

Spin-Orbit Misalignments in Tertiary-Induced Black-Hole Binary Mergers: Theoretical Analysis

Authors: Yubo Su, Dong Lai, Bin Liu

Abstract: Black-hole (BH) binary mergers driven by gravitational perturbations of tertiary companions constitute an important class of dynamical formation channels for compact binaries detected by LIGO/VIRGO. Recent works have examined numerically the combined orbital and spin dynamics of BH binaries that undergo large Lidov-Kozai (LK) eccentricity oscillations induced by a highly inclined companion and mer… ▽ More Black-hole (BH) binary mergers driven by gravitational perturbations of tertiary companions constitute an important class of dynamical formation channels for compact binaries detected by LIGO/VIRGO. Recent works have examined numerically the combined orbital and spin dynamics of BH binaries that undergo large Lidov-Kozai (LK) eccentricity oscillations induced by a highly inclined companion and merge via gravitational wave radiation. However, the extreme eccentricity variations make such systems difficult to characterize analytically. In this paper, we develop an analytical formalism for understanding the spin dynamics of binary BHs undergoing LK-induced mergers. We show that, under certain conditions, the eccentricity oscillations of the binary can be averaged over to determine the long-term behavior of the BH spin in a smooth way. In particular, we demonstrate that the final spin-orbit misalignment angle $θ_{\rm sl}$ is often related to the binary's primordial spin orientation through an approximate adiabatic invariant. Our theory explains the "$90^\circ$ attractor" (as found in recent numerical studies) for the evolution of $θ_{\rm sl}$ when the initial BH spin is aligned with the orbital axis and the octupole LK effects are negligible -- such a "$90^\circ$ attractor" would lead to a small binary effective spin parameter $χ_{\rm eff}\sim 0$ even for large intrinsic BH spins. We calculate the deviation from adiabaticity in closed form as a function of the initial conditions. We also place accurate constraints on when this adiabatic invariant breaks down due to resonant spin-orbit interactions. We consider both stellar-mass and supermassive BH tertiary companions, and provide simple prescriptions for determining analytically the final spin-orbit misalignment angles of the merging BH binaries. △ Less

Submitted 24 February, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 21 pages, 17 figures, accepted by PRD

Journal ref: Phys. Rev. D 103, 063040 (2021)

Inclination Dynamics of Resonant Planets under the Influence of an Inclined External Companion

Authors: Laetitia Rodet, Dong Lai

Abstract: Recent observations suggest that a large fraction of Kepler super-Earth systems have external giant planet companions (cold Jupiters), which can shape the architecture of the inner planets, in particular their mutual inclinations. The dynamical perturbation from cold Jupiters may account for the population of misaligned planets in the Kepler data. The effectiveness of this mechanism can be hindere… ▽ More Recent observations suggest that a large fraction of Kepler super-Earth systems have external giant planet companions (cold Jupiters), which can shape the architecture of the inner planets, in particular their mutual inclinations. The dynamical perturbation from cold Jupiters may account for the population of misaligned planets in the Kepler data. The effectiveness of this mechanism can be hindered by a strong planet--planet coupling in the inner system. In this paper, we study how mean-motion resonances (MMRs) affect this coupling and the expected misalignment. We derive approximate analytical expressions for the mutual inclination excitations in the inner planet system induced by an inclined companion, for various period ratios and perturber properties. In most cases, the mutual inclination is proportional to a dimensionless parameter that characterizes the strength of the perturber relative to the coupling in the inner system. We show that the MMR strengthens the inner coupling, reducing the mutual inclination induced by the perturber by a factor of a few. We find that the resonance is resilient to the perturbation, and derive a criterion for the libration of the resonant angle. Our results have applications for constraining unseen planetary perturbers, and for understanding the architecture of multiplanet systems. △ Less

Submitted 5 March, 2021; v1 submitted 15 October, 2020; originally announced October 2020.

Comments: 17 pages, 14 figures, accepted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 502, Issue 3, April 2021, Pages 3746-3760

Enhanced Lidov-Kozai migration and the formation of the transiting giant planet WD1856+534b

Authors: Christopher E. O'Connor, Bin Liu, Dong Lai

Abstract: We investigate the possible origin of the transiting giant planet WD1856+534b, the first strong exoplanet candidate orbiting a white dwarf, through high-eccentricity migration (HEM) driven by the Lidov-Kozai (LK) effect. The host system's overall architecture is an hierarchical quadruple in the '2+2' configuration, owing to the presence of a tertiary companion system of two M-dwarfs. We show that… ▽ More We investigate the possible origin of the transiting giant planet WD1856+534b, the first strong exoplanet candidate orbiting a white dwarf, through high-eccentricity migration (HEM) driven by the Lidov-Kozai (LK) effect. The host system's overall architecture is an hierarchical quadruple in the '2+2' configuration, owing to the presence of a tertiary companion system of two M-dwarfs. We show that a secular inclination resonance in 2+2 systems can significantly broaden the LK window for extreme eccentricity excitation ($e \gtrsim 0.999$), allowing the giant planet to migrate for a wide range of initial orbital inclinations. Octupole effects can also contribute to the broadening of this 'extreme' LK window. By requiring that perturbations from the companion stars be able to overcome short-range forces and excite the planet's eccentricity to $e \simeq 1$, we obtain an absolute limit of $a_{1} \gtrsim 8 \, {\rm AU} \, (a_{3} / 1500 \, {\rm AU})^{6/7}$ for the planet's semi-major axis just before migration (where $a_{3}$ is the semi-major axis of the 'outer' orbit). We suggest that, to achieve a wide LK window through the 2+2 resonance, WD1856b likely migrated from $30 \, {\rm AU} \lesssim a_{1} \lesssim 60 \, {\rm AU}$, corresponding to $\sim 10$--$20 \, {\rm AU}$ during the host's main-sequence phase. We discuss possible difficulties of all flavours of HEM affecting the occurrence rate of short-period giant planets around white dwarfs. △ Less

Submitted 11 January, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 8 pages, 4 figures, 1 table; MNRAS 501 (2021) 507-514

Hierarchical Black-Hole Mergers in Multiple Systems: Constrain the Formation of GW190412, GW190814 and GW190521-like events

Authors: Bin Liu, Dong Lai

Abstract: The merging black-hole (BH) binaries GW190412, GW190814 and GW190521 from the third LIGO/VIRGO observing run exhibit some extraordinary properties, including highly asymmetric masses, significant spin, and component mass in the "mass gap". These features can be explained if one or both components of the binary are the remnants of previous mergers. In this paper, we explore hierarchical mergers in… ▽ More The merging black-hole (BH) binaries GW190412, GW190814 and GW190521 from the third LIGO/VIRGO observing run exhibit some extraordinary properties, including highly asymmetric masses, significant spin, and component mass in the "mass gap". These features can be explained if one or both components of the binary are the remnants of previous mergers. In this paper, we explore hierarchical mergers in multiple stellar systems, taking into account the natal kick and mass loss due to the supernova explosion (SN) on each component, as well as the merger kick received by the merger remnant. The binaries that have survived the SNe and kicks generally have too wide orbital separations to merge by themselves, but can merge with the aid of an external companion that gives rise to Lidov-Kozai oscillations. The BH binaries that consist of second-generation BHs can also be assembled in dense star clusters through binary interactions. We characterize the parameter space of these BH binaries by merger fractions in an analytical approach. Combining the distributions of the survived binaries, we further constrain the parameters of the external companion, using the analytically formulated tertiary perturbation strength. We find that to produce the three LIGO/VIRGO O3 events, the external companions must be at least a few hundreds $M_\odot$, and fall in the intermediate-mass BH and supermassive BH range. We suggest that GW190412, GW190814 and GW190521 could all be produced via hierarchical mergers in multiples, likely in a nuclear star cluster, with the final merger induced by a massive BH. △ Less

Submitted 20 January, 2021; v1 submitted 21 September, 2020; originally announced September 2020.

Comments: 19 pages, 20 figures, accepted for publication in MNRAS

Evolution of inspiralling neutron star binaries: effects of tidal interactions and orbital eccentricities

Authors: Jie-Shuang Wang, Dong Lai

Abstract: Neutron star (NS) binaries formed dynamically may have significant eccentricities while emitting gravitational waves (GWs) in the LIGO/VIRGO band. We study tidal effects in such eccentric inspiralling NS binaries using a set of hybrid equations. The NS is modelled as a compressible ellipsoid, which can deform nonlinearly due to tidal forces, while the orbit evolution is treated with the post-Newto… ▽ More Neutron star (NS) binaries formed dynamically may have significant eccentricities while emitting gravitational waves (GWs) in the LIGO/VIRGO band. We study tidal effects in such eccentric inspiralling NS binaries using a set of hybrid equations. The NS is modelled as a compressible ellipsoid, which can deform nonlinearly due to tidal forces, while the orbit evolution is treated with the post-Newtonian (PN) theory up to 2.5-PN order. We find that in general, the tidal interaction can accelerate the inspiral, and cause orbital frequency and phase shifts. For circular inspirals, our calculations reproduce previous linear result at large binary separations, but incorporate the dynamical response of the NS at small separations. For eccentric inspirals, the frequency and phase shifts oscillate considerably near pericenter passages, and the oscillating amplitudes increase with eccentricities. As a result, the GW phase is also significantly influenced by the tidal effect. At merger, the cumulative GW phase shift can reach more than 10 radians (for typical NS mass $1.4M_\odot$ and radius 11.6 km), much larger than the circular inspiral case. Although the event rate of eccentric NS mergers is likely low, the detection of such mergers could provide a useful constraint on the NS equation of state. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 9 pages, 8 figures, accepted to Phys. Rev. D

Strong Scatterings of Cold Jupiters and their Influence on Inner Low-mass Planet Systems: Theory and Simulations

Authors: Bonan Pu, Dong Lai

Abstract: Recent observations have indicated a strong connection between compact ($a \lesssim 0.5$ au) super-Earth and mini-Neptune systems and their outer ($a \gtrsim$ a few au) giant planet companions. We study the dynamical evolution of such inner systems subject to the gravitational effect of an unstable system of outer giant planets, focussing on systems whose end configurations feature only a single r… ▽ More Recent observations have indicated a strong connection between compact ($a \lesssim 0.5$ au) super-Earth and mini-Neptune systems and their outer ($a \gtrsim$ a few au) giant planet companions. We study the dynamical evolution of such inner systems subject to the gravitational effect of an unstable system of outer giant planets, focussing on systems whose end configurations feature only a single remaining outer giant. In contrast to similar studies which used on N-body simulations with specific (and limited) parameters or scenarios, we implement a novel hybrid algorithm which combines N-body simulations with secular dynamics with aims of obtaining analytical understanding and scaling relations. We find that the dynamical evolution of the inner planet system depends crucially on $N_{\mathrm{ej}}$, the number of mutual close encounters between the outer planets prior to eventual ejection/merger. When $N_{\mathrm{ej}}$ is small, the eventual evolution of the inner planets can be well described by secular dynamics. For larger values of $N_{\mathrm{ej}}$, the inner planets gain orbital inclination, inclination and eccentricity in a stochastic fashion analogous to Brownian motion. We develop a theoretical model, and compute scaling laws for the final orbital parameters of the inner system. We show that our model can account for the observed eccentric super-Earths/mini-Neptunes with inclined cold Jupiter companions, such as HAT-P-11, Gliese 777 and $π$ Men. △ Less

Submitted 2 September, 2021; v1 submitted 13 August, 2020; originally announced August 2020.

Comments: accepted for publication in MNRAS

Tidal Dissipation Impact on the Eccentric Onset of Common Envelope Phases in Massive Binary Star Systems

Authors: Michelle Vick, Morgan MacLeod, Dong Lai, Abraham Loeb

Abstract: Tidal dissipation due to turbulent viscosity in the convective regions of giant stars plays an important role in shaping the orbits of pre-common envelope systems. Such systems are possible sources of transients and close compact binary systems that will eventually merge and produce detectable gravitational wave signals. Most previous studies of the onset of common envelope episodes have focused o… ▽ More Tidal dissipation due to turbulent viscosity in the convective regions of giant stars plays an important role in shaping the orbits of pre-common envelope systems. Such systems are possible sources of transients and close compact binary systems that will eventually merge and produce detectable gravitational wave signals. Most previous studies of the onset of common envelope episodes have focused on circular orbits and synchronously rotating donor stars under the assumption that tidal dissipation can quickly spin up the primary and circularize the orbit before the binary reaches Roche-lobe overflow (RLO). We test this assumption by coupling numerical models of the post main sequence stellar evolution of massive-stars with the model for tidal dissipation in convective envelopes developed in Vick & Lai (2020) $-$ a tidal model that is accurate even for highly eccentric orbits with small pericentre distances. We find that, in many cases, tidal dissipation does not circularize the orbit before RLO. For a $10~M_\odot$ ($15~M_\odot$) primary star interacting with a $1.4~M_\odot$ companion, systems with pericentre distances within 3 AU (6 AU) when the primary leaves the main sequence will retain the initial orbital eccentricity when the primary grows to the Roche radius. Even in systems that tidally circularize before RLO, the donor star may be rotating subsynchronously at the onset of mass transfer. Our results demonstrate that some possible precursors to double neutron star systems are likely eccentric at the Roche radius. The effects of pre-common envelope eccentricity on the resulting compact binary merit further study. △ Less

Submitted 23 March, 2021; v1 submitted 12 August, 2020; originally announced August 2020.

Comments: 15 pages, 9 figures, Accepted to MNRAS

Giant Planet Scatterings and Collisions: Hydrodynamics, Merger-Ejection Branching Ratio, and Properties of the Remnants

Authors: Jiaru Li, Dong Lai, Kassandra R. Anderson, Bonan Pu

Abstract: Planetary systems with sufficiently small orbital spacings can experience planetary mergers and ejections. The branching ratio of mergers vs ejections depends sensitively on the treatment of planetary close encounters. Previous works have adopted a simple "sticky-sphere" prescription, whose validity is questionable. We apply both smoothed particle hydrodynamics and $N$-body integrations to investi… ▽ More Planetary systems with sufficiently small orbital spacings can experience planetary mergers and ejections. The branching ratio of mergers vs ejections depends sensitively on the treatment of planetary close encounters. Previous works have adopted a simple "sticky-sphere" prescription, whose validity is questionable. We apply both smoothed particle hydrodynamics and $N$-body integrations to investigate the fluid effects in close encounters between gas giants and the long-term evolution of closely-packed planetary systems. Focusing on parabolic encounters between Jupiter-like planets with $M_J$ and $2M_J$, we find that quick mergers occur when the impact parameter $r_p$ (the pericenter separation between the planets) is less than $2R_J$, and the merger conserved 97% of the initial mass. Strong tidal effects can affect the "binary-planet" orbit when $r_p$ is between $2R_J$ and $4R_J$. We quantify these effects using a set of fitting formulae that can be implemented in $N$-body codes. We run a suite of $N$-body simulations with and without the formulae for systems of two giant planets initially in unstable, nearly circular coplanar orbits. The fluid (tidal) effects significantly increase the branching ratio of planetary mergers relative to ejections by doubling the effective collision radius. While the fluid effects do not change the distributions of semi-major axis and eccentricity of each type of remnant planets (mergers vs surviving planets in ejections), the overall orbital properties of planet scattering remnants are strongly affected due to the change in the branching ratio. We also find that the merger products have broad distributions of spin magnitudes and obliquities. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Comments: 13 pages, 14 Figures, 4 Tables, Submitted to MNRAS

Planetary Spin and Obliquity from Mergers

Authors: Jiaru Li, Dong Lai

Abstract: In planetary systems with sufficiently small inter-planet spacing, close encounters can lead to planetary collisions/mergers or ejections. We study the spin property of the merger products of two giant planets in a statistical manner using numerical simulations and analytical modeling. Planetary collisions lead to rapidly rotating objects and a broad range of obliquities. We find that, under typic… ▽ More In planetary systems with sufficiently small inter-planet spacing, close encounters can lead to planetary collisions/mergers or ejections. We study the spin property of the merger products of two giant planets in a statistical manner using numerical simulations and analytical modeling. Planetary collisions lead to rapidly rotating objects and a broad range of obliquities. We find that, under typical conditions for two-planet scatterings, the distributions of spin magnitude $S$ and obliquity $Θ_{\rm SL}$ of the merger products have simple analytical forms: $f_{S} \propto S$ and $f_{\cosΘ_{\rm SL}} \propto (1-\cos^2Θ_{\rm SL})^{-1/2}$. Though parameter studies, we determine the regime of validity for the analytical distributions of spin and obliquity. Since planetary mergers is a major outcome of planet-planet scatterings, observational search for the spin/obliquity signatures of exoplanets would provide important constraints on the dynamical history of planetary systems. △ Less

Submitted 15 May, 2020; originally announced May 2020.