Researcher profile

Dong Lai

Dong Lai contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Dynamical Instability in Multi-Orbiter Systems with Gas Friction

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.

preprint2022arXiv

Evolution of Stellar Orbits Around Merging Massive Black-Hole Binary

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.

preprint2022arXiv

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

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 &#34;closeness&#34; and the property of the resulting binary BH mergers. For a typical &#34;SMBH + 2 BHs&#34; 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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

Secular oscillations in multi-planet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This &#34;secular chaos&#34; 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&#39; 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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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&#39;s overall architecture is an hierarchical quadruple in the &#39;2+2&#39; 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 &#39;extreme&#39; LK window. By requiring that perturbations from the companion stars be able to overcome short-range forces and excite the planet&#39;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&#39;s semi-major axis just before migration (where $a_{3}$ is the semi-major axis of the &#39;outer&#39; 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&#39;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.

preprint2021arXiv

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

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 &#34;mass gap&#34;. 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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Fast radio burst counterparts and their implications for the central engine

While the radiation mechanism of fast radio bursts (FRBs) is unknown, coherent curvature radiation and synchrotron maser are promising candidates. We find that both radiation mechanisms work for a neutron star (NS) central engine with $B\gtrsim 10^{12}$ G, while for the synchrotron maser, the central engine can also be an accreting black hole (BH) with $B\gtrsim 10^{12}$ G and a white dwarf (WD) with $B\sim 10^8-10^9$ G. We study the electromagnetic counterparts associated with such central engines, i.e., nebulae for repeating FRBs and afterglows for non-repeating FRBs. In general, the energy spectrum and flux density of the counterpart depend strongly on its size and total injected energy. We apply the calculation to the nebula of FRB 121102 and find that the persistent radio counterpart requires the average energy injection rate into the nebula to be between $2.7\times10^{39}~{\rm erg/s}$ and $1.5\times10^{44}~{\rm erg/s}$, and the minimum injected energy be $6.0\times10^{47}~{\rm erg}$ in around $7$ yr. Consequently, we find that for FRB 121102 and its nebula: (1) WD and accretion BH central engines are disfavored; (2) a rotation-powered NS central engine works when $1.2\times10^{12}~{\rm G}\lesssim B\lesssim 7.8\times10^{14}~{\rm G}$ with initial period $P<180$ ms, but the radio emission must be more efficient than that in typical giant pulses of radio pulsars; and (3) a magnetic-powered NS central engine works when its internal magnetic field $B\gtrsim 10^{16}$ G. We also find that the radio-emitting electrons in the nebula could produce a significant rotation measure (RM), but cannot account for the entire observed RM of FRB 121102.

preprint2020arXiv

Merging Compact Binaries Near a Rotating Supermassive Black Hole: Eccentricity Excitation due to Apsidal Precession Resonance

We study the dynamics of merging compact binaries near a rotating supermassive black hole (SMBH) in a hierarchical triple configuration. We include various general relativistic effects that couple the inner orbit, the outer orbit and the spin of the SMBH. During the orbital decay due to gravitational radiation, the inner binary can encounter an &#34;apsidal precession resonance&#34; and experience eccentricity excitation. This resonance occurs when the apsidal precession rate of the inner binary matches that of the outer binary, with the precessions driven by both Newtonian interactions and various post-Newtonian effects. The eccentricity excitation requires the outer orbit to have a finite eccentricity, and is most effective for triples with small mutual inclinations, in contrast to the well-studied Lidov-Kozai effect. The resonance and the associated eccentricity growth may occur while the binary emits gravitational waves in the low-frequency band, and may be detectable by future space-based gravitational wave detectors.

preprint2020arXiv

Periodic Fast Radio Bursts with Neutron Star Free/Radiative Precession

The CHIME/FRB collaboration recently reported the detection of a 16 day periodicity in the arrival times of radio bursts from FRB 180916.J0158+65. We study the possibility that the observed periodicity arises from free precession of a magnetized neutron star, and put constraints on different components of the star&#39;s magnetic fields. Using a simple geometric model, where radio bursts are emitted from a rotating neutron star magnetosphere, we show that the emission pattern as a function of time can match that observed from FRB 180916.J0158+65.

preprint2020arXiv

Planetary Spin and Obliquity from Mergers

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.

preprint2020arXiv

Tidal Evolution of Eccentric Binaries Driven by Convective Turbulent Viscosity

Tidal dissipation due to convective turbulent viscosity shapes the evolution of a variety of astrophysical binaries. For example, this type of dissipation determines the rate of orbital circularization in a binary with a post-main sequence star that is evolving toward a common envelope phase. Viscous dissipation can also influence binaries with solar-type stars, or stars with a close-in giant planet. In general, the effective viscosity in a convective stellar envelope depends on the tidal forcing frequency $ω_{\rm tide}$; when $ω_{\rm tide}$ is larger than the turnover frequency of convective eddies, the viscosity is reduced. Previous work has focused on binaries in nearly circular orbits. However, for eccentric orbits, the tidal potential has many forcing frequencies. In this paper, we develop a formalism for computing tidal dissipation that captures the effect of frequency-dependent turbulent viscosity and is valid for arbitrary binary eccentricities. We also present an alternative simpler formulation that is suitable for very high eccentricities. We apply our formalisms to a giant branch (GB) star model and a solar-type star model. We find that a range of pseudosynchronous rotation rates are possible for both stellar models, and the pseudosynchronous rate can differ from the prediction of the commonly-used weak tidal friction theory by up to a factor of a few. We also find that tidal decay and circularization due to turbulent viscosity can be a few orders of magnitude faster than predicted by weak tidal friction in GB stars on eccentric, small pericentre orbits, but is suppressed by a few orders of magnitude in solar-type stars due to viscosity reduction.

preprint2019arXiv

Circumbinary Accretion from Finite and Infinite Disks

We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite&#39;&#39;, steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite&#39;&#39; disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${\sim}0.2Ω_{\rm{b}}$ for $q_{\rm{b}}>0.5$ and $Ω_{\rm{b}}$ for $q_{\rm{b}}<0.5$, ($Ω_{\rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue&#39;&#39; $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{\rm b}^2Ω_{\rm{b}}$ (with $a_{\rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{\rm{b}}\gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite&#39;&#39; disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.

preprint2019arXiv

In-Situ Scattering of Warm Jupiters and Implications for Dynamical Histories

Many warm Jupiters (WJs) have substantial eccentricities, which are linked to their formation and migration histories. This paper explores eccentricity excitation of WJs due to planet-planet scattering, beginning with 3-4 planets in unstable orbits, with the innermost planet placed in the range (0.1 - 1)AU. Such a setup is consistent with either in-situ formation or arrival at sub-AU orbits due to disk migration. Most previous N-body experiments have focused on &#34;cold&#34; Jupiters at several AU, where scattering results in planet ejections, efficiently exciting the eccentricities of surviving planets. In contrast, scattering at sub-AU distances results in a mixture of collisions and ejections, and the final eccentricities of surviving planets are unclear. We conduct scattering experiments for a range of planet masses and initial spacings, including the effect of general relativistic apsidal precession, and systematically catalogue the scattering outcomes and properties of surviving planets. A comparable number of one-planet and two-planet systems are produced. Two-planet systems arise exclusively through planet-planet collisions, and tend to have low eccentricities/mutual inclinations and compact configurations. One-planet systems arise through a combination of ejections and collisions, resulting in higher eccentricities. The observed eccentricity distribution of solitary WJs (lacking detection of a giant planet companion) is consistent with roughly 60% of the systems having undergone in-situ scattering, and the remaining experiencing a quiescent history.

preprint2019arXiv

Pulsed Disc Accretion Driven by Hot Jupiters

We present 2D hydrodynamical simulations of hot Jupiters orbiting near the inner edge of protoplanetary discs. We systemically explore how the accretion rate at the inner disc edge is regulated by a giant planet of different mass, orbital separation and eccentricity. We find that a massive (with planet-to-star mass ratio $\gtrsim 0.003$) eccentric ($e_p\gtrsim 0.1$) planet drives a pulsed accretion at the inner edge of the disc, modulated at one or two times the planet&#39;s orbital frequency. The amplitude of accretion variability generally increases with the planet mass and eccentricity, although some non-monotonic dependences are also possible. Applying our simulation results to the T Tauri system CI Tau, where a young hot Jupiter candidate has been detected, we show that the observed luminosity variability in this system can be explained by pulsed accretion driven by an eccentric giant planet.