Theoretical astrophysics and plasma physics at RPC

Abstract:

Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.

Abstract:

Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system's chaotic nature. We show that an ensemble of converged solutions for the planar chaotic three-body problem obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at a fixed computational cost and up to 100 million times faster than the numerical integrator. In addition, we demonstrate the importance of training an ANN using converged solutions from an arbitrary precise integrator, relative to solutions computed by a conventional fixed precision integrator, which can introduce errors in the training data, due to numerical round-off and time discretization, that are learned by the ANN. Our results provide evidence that, for computationally challenging regions of phase space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black hole binary systems or the origin of the core collapse in dense star clusters.