Many stars do not live alone, but instead have one or more stellar companions. Observations show that these binaries, triples and higher-order multiples are common. Whereas the evolution of single stars and binaries have been studied extensively, the same is not true for the evolution of stellar triples. To fill this gap in our general understanding of stellar lives, we aim to systematically explore the long-term evolution of triples and to map out the most common evolutionary pathways that triples go through. We quantitatively study how triples evolve, which processes are most relevant, and how this differs from binary evolution. We simulate the evolution of several large populations of triples with a population synthesis approach. We make use of the triple evolution code TRES to simulate the evolution of each triple in a consistent way; including three-body dynamics (based on the secular approach), stellar evolution and their mutual influences. We simulate the evolution of the system up until mass transfer starts, the system becomes dynamically unstable, or a Hubble time has passed. We find that stellar interactions are common in triples. Compared to a binary population, we find that the fraction of systems that can undergo mass transfer is about 2 to 3 times larger in triples. Moreover, whereas in binaries the orbits typically reach circularisation before Roche-lobe overflow, this is not true anymore in triples. In our simulations, about 40% of systems retain an eccentric orbit. Additionally, we discuss various channels of triple evolution in detail such as those where the secondary or the tertiary is the first star to initiate a mass transfer event.

By performing Monte Carlo simulations of the evolution of binary primordial black hole (PBH) systems, we estimate the effect of distant encounters with single PBHs upon the coalescence time and merger rate of binary PBHs. We find that, for models where PBHs compose a large fraction of dark matter, $\fPBH\sim 1$, the expected fractional change in coalescence time is negligible, of order $10^{-6}$ for most binaries. For models with significantly lower PBH abundances, $\fPBH\ll 1$, we find that the average change in binary lifetime due to encounters can be as large as $\mathcal{O}(10^{-2})$, with a small number of binaries experiencing an order unity change in lifetime. In the absence of encounters, we also compare the use of an analytic approximation for the coalescence time to numerically evolving the binary system, finding that the analytic approximation results in an order $10\%$ error in the coalescence time. However, when these effects are taken into consideration, there is a negligible change to the calculated merger rate, placing previous constraints on the PBH abundance arising from observed gravitational wave signals from merging binary black holes on a more secure footing.

The gravitational wave source GW190412 is a binary black hole (BBH) merger with three unique properties: i) its mass ratio is about 0.28, the lowest found so far, ii) it has a relatively high positive effective spin parameter \chi_eff=0.25, and iii) it is observed to be precessing due to in-plane projected spin of the binary with an in-plane precession parameter \chi_p=0.3. The two main formation channels of BBH formation fail to account for GW190412: field formation scenarios cannot explain the observed precession unless by invoking large natal kicks, and dynamical assembly in dense stellar systems is inefficient in producing such low mass-ratio BBH mergers. Here, we investigate whether "double mergers" in wide hierarchical quadruple systems in the "3+1" configuration could explain the unique properties of GW190412. In this scenario, a compact object quadruple system experiences two mergers: first, two compact objects in the innermost orbit merge due to secular chaotic evolution. At a later time, the merged compact object coalesces with another compact object due to secular Lidov-Kozai oscillations. We find that our scenario is consistent with GW190412. In particular, we find a preferential projected spin around \chi_p=0.2. However, the likelihood of a double merger is small and the formation efficiency of these systems is uncertain. If GW190412 originated from a double merger in a 3+1 quadruple, we find a strong constraint that the first merger likely occurred between roughly equal-mass BHs in the innermost orbit, since the recoil velocity from unequal-mass BHs would otherwise have disrupted the system.

Statistics of hierarchical systems containing three or more stars are continuously improving. The Multiple Star Catalog (MSC) is currently the most comprehensive catalogue of multiple-star systems and contains component masses, orbital periods, and additional information. The systems in the MSC are interesting for several reasons, including the long-term dynamical evolution of few-body systems. Although the secular evolution of triples and quadruples has been explored before, a systematic study of the systems in the MSC including also quintuples and sextuples has not been carried out. Here, we explore the main-sequence (MS) evolution of stars from the MSC based on approximately 2x10^5 secular dynamical integrations. We estimate statistical probabilities for strong interactions during the MS such as tidal evolution and mass transfer, and the onset of dynamical instability. Depending on the assumed model for the unknown orbital elements, we find that the fraction of noninteracting systems is largest for triples (~0.9), and decreases to ~0.6-0.8 for sextuples. The fraction of strong interactions increases from ~0.1 to ~0.2 from triples to sextuples, and the fraction of dynamically unstable systems increases from ~0.001 to ~0.1-0.2. The larger fractions of strong interactions and dynamical instability in systems with increasing multiplicity can be attributed to increasingly complex secular evolution in these systems. Our results indicate that a significant fraction of high-multiplicity systems interact or become dynamically unstable already during the MS, with an increasing importance as the number of stars increases.

The SecularMultiple code, presented in two previous papers of this series, integrates the long-term dynamical evolution of multiple systems with any number of bodies and hierarchical structure, provided that the system is composed of nested binaries. In the formalism underlying SecularMultiple, we previously averaged over all orbits in the system. This approximation significantly speeds up numerical integration of the equations of motion, making large population synthesis studies possible. However, the orbit averaging approximation can break down when the secular evolution timescale of the system is comparable to or shorter than any of the orbital periods in the system. Here, we present an update to SecularMultiple in which we incorporate hybrid integration techniques, and orbit-averaging corrections. With this update, the user can specify which orbits should be integrated directly (without averaging), or assuming averaged orbits. For orbits that are integrated directly, we implemented two integration techniques, one which is based on the regularised Kustaanheimo-Stiefel equations of motion in element form. We also implemented analytical orbit-averaging corrections for pairwise interactions to quadrupole order. The updates presented here provide more flexibility for integrating the long-term dynamical evolution of hierarchical multiple systems. By effectively combining direct integration and orbit averaging the long-term evolution can be accurately computed, but with significantly lower computational cost compared to existing direct N-body codes. We give a number of examples in which the new features are beneficial. Our updated code, which is written in C++ supplemented with a user-friendly interface in Python, is freely available.

Binary-binary interactions are important in a number of astrophysical contexts including dense stellar systems such as globular clusters. Although less frequent than binary-single encounters, binary-binary interactions lead to a much richer range of possibilities such as the formation of stable triple systems. Here, we focus on the regime of distant binary-binary encounters, i.e., two binaries approaching each other on an unbound orbit with a periapsis distance Q much larger than the internal binary separations. This `secular' regime gives rise to changes in the orbital eccentricities and orientations, which we study using analytic considerations and numerical integrations. We show that `direct' interactions between the three orbits only occur starting at a high expansion order of the Hamiltonian (hexadecupole order), and that the backreaction of the outer orbit on the inner two orbits at lower expansion orders is weak. Therefore, to good approximation, one can obtain the changes of each orbit by using previously-known analytic results for binary-single interactions, and replacing the mass of the third body with the total mass of the companion binary. Nevertheless, we find some dependence of the `binarity' of the companion binary, and derive explicit analytic expressions for the secular changes that are consistent with numerical integrations. In particular, the eccentricity and inclination changes of orbit 1 due to orbit 2 scale as \eps_{SA1} (a_2/Q)^2 [m_3 m_4/(m_3+m_4)^2], where \eps_{SA1} is the approximate quadrupole-order change, and a_2 and (m_3,m_4) are the companion binary orbital semimajor axis and component masses, respectively. Our results are implemented in several Python scripts that are freely available.

We investigate secular evolution in hierarchical quadruple systems as a formation channel of mass-gap black holes (with masses of about $3-5$ $M_{\odot}$) in systems that will eventually lead to binary black hole mergers detectable by ground-based gravitational wave detectors (LIGO/Virgo). We show that in a 3+1 hierarchical system, two episodes of induced mergers would first cause two neutron stars to merge and form a mass-gap black hole, which will subsequently merge with another (more massive) black hole through a second induced merger. We demonstrate that such systems are stable to flybys, and their formation would predict a high mass ratio and eccentric merger of a mass-gap black hole with a more massive black hole companion. Such a formation channel may explain observed gravitational wave events such as the recently-discovered LIGO/Virgo events S190814bv and S190924h.

The isolated binary evolution model for merging neutron stars (NSs) involves processes such as mass transfer, common-envelope evolution, and natal kicks, all of which are poorly understood. Also, the predicted NS-NS merger rates are typically lower than the rates inferred from the LIGO GW170817 event. Here, we investigate merger rates of NS and black hole (BH)-NS binaries in hierarchical triple-star systems. In such systems, the tertiary can induce Lidov-Kozai (LK) oscillations in the inner binary, accelerating its coalescence, and potentially enhancing compact object merger rates. However, since compact objects originate from massive stars, the prior evolution should also be taken into account. Natal kicks, in particular, could significantly reduce the rates by unbinding the tertiary before it can affect the inner binary through LK evolution. We carry out simulations of massive triples taking into account stellar evolution starting from the main sequence, secular and tidal evolution, and the effects of supernovae. For large NS birth kicks (\sigma_k = 265 km/s), we find that the triple NS-NS merger rate (several hundred Gpc^{-3} yr^{-1}) is lower by a factor of ~2-3 than the binary rate, but for no kicks (\sigma_k = 0 km/s), the triple rate (several thousand Gpc^{-3} yr^{-1}) is comparable to the binary rate. Our results indicate that a significant fraction of NS-NS mergers could originate from triples if a substantial portion of the NS population is born with low kick velocities, as indicated by other work. However, uncertainties and open questions remain because of our simplifying assumption of dynamical decoupling after inner binary interaction has been triggered.

Dense stellar systems such as globular clusters are believed to harbor merging binary black holes (BHs). The evolution of such binaries is driven by interactions with other stars, most notably, binary-single interactions. Traditionally, so-called `strong' interactions are believed to be the driving force in this evolution. However, we recently showed that more distant, i.e., "weak" or "secular" encounters, can have important implications for the properties of merging BH binaries in globular clusters. This motivates more detailed understanding of the effects of secular encounters on a binary. In another previous paper, we analytically calculated expressions for the changes of the eccentricity and angular-momentum vectors taking into account second-order (SO) perturbation theory, and showed that, for highly eccentric binaries, the new expressions give rise to behavior that is not captured by first-order (FO) theory. Here, we extend our previous work to third order (TO) perturbation theory. We also include terms up to and including octupole order. The latter are nonzero for binaries with unequal component masses. In addition, we consider the effects of post-Newtonian terms, and we determine the steady-state distribution due to the cumulative effect of secular encounters by computing the associated angular-momentum diffusion coefficients, and applying the Fokker-Planck equation. Together with our previous work, the results in this paper provide a framework for incorporating the effects of distant encounters on binaries in models of cluster evolution, such as Monte Carlo codes.

Type Ia supernovae (SNe Ia) could arise from mergers of carbon-oxygen white dwarfs (WDs) triggered by Lidov-Kozai (LK) oscillations in hierarchical triple-star systems. However, predicted merger rates are several orders of magnitude lower than the observed SNe Ia rate. The low predicted rates can be attributed in part to the fact that many potential WD-merger progenitor systems, with high mutual orbital inclination, merge or interact before the WD stage. Recently, evidence was found for the existence of natal kicks imparted on WDs with a typical magnitude of 0.75 km/s. In triples, kicks change the mutual inclination and in general increase the outer orbit eccentricity, bringing the triple into an active LK regime at late stages and avoiding the issue of pre-WD merger or interaction. Stars passing by the triple can result in similar effects. However, both processes can also disrupt the triple. In this paper, we quantitatively investigate the impact of WD kicks and flybys on the rate of WD mergers using detailed simulations. We find that WD kicks and flybys combine to increase the predicted WD merger rates by a factor of ~2.5, resulting in a time-integrated rate of ~1.1e-4 Msun^-1. Despite the significant boost, the predicted rates are still more than one order of magnitude below the observed rate of ~1e-3 Msun^-1. However, many systematic uncertainties still remain in our calculations, in particular the potential contributions from tighter triples, dynamically unstable systems, unbound systems due to WD kicks, and quadruple systems.

Binary-single interactions play a crucial role in the evolution of dense stellar systems such as globular clusters. In addition, they are believed to drive black hole (BH) binary mergers in these systems. A subset of binary-single interactions are secular encounters, for which the third body approaches the binary on a relatively wide orbit, and such that it is justified to average the equations of motion over the binary's orbital phase. Previous works used first-order perturbation theory to compute the effects of such secular encounters on the binary. However, this approach can break down for highly eccentric binaries, which are important for BH binary mergers and gravitational wave sources. Here, we present an analytic computation using second-order perturbation techniques, valid to the quadrupole-order approximation. In our calculation, we take into account the instantaneous back-reaction of the binary to the third body, and compute corrections to previous first-order results. Using singly-averaged and direct 3-body integrations, we demonstrate the validity of our expressions. In particular, we show that the eccentricity change for highly eccentric binaries can reach a plateau, associated with a large inclination change, and can even reverse sign. These effects are not captured by previous first-order results. We provide a simple script to conveniently evaluate our analytic expressions, including routines for numerical integration and verification.

We present the first systematic study on how distant weak interactions impact the dynamical evolution of merging binary black holes (BBHs) in dense stellar clusters. Recent studies indicate that dense clusters are likely to significantly contribute to the rate of merging BBHs observable through gravitational waves (GWs), and that many of these mergers will appear with notable eccentricities measurable in the LISA and LIGO sensitivity bands. This is highly interesting, as eccentricity can be used to distinguish between different astrophysical merger channels. However, all of these recent studies are based on various Monte Carlo (MC) techniques that only include strong interactions for the dynamical evolution of BBHs, whereas any binary generally undergoes orders-of-magnitude more weak interactions than strong. It is well known that weak interactions primarily lead to a change in the binary's eccentricity, which for BBHs implies that weak interactions can change their GW inspiral time and thereby their merger probability. With this motivation, we perform MC simulations of BBHs evolving in dense clusters under the influence of both weak and strong interactions. We find that including weak interactions leads to a notable increase in the number of BBHs that merge inside their cluster, which correspondingly leads to a higher number of eccentric LISA sources. These preliminary results illustrate the importance of including weak interactions for accurately modeling how BBHs merge in clusters, and how to link their emitted GW signals to their astrophysical environment.

Most studies of mass transfer in binary systems assume circular orbits at the onset of Roche lobe overflow. However, there are theoretical and observational indications that mass transfer could occur in eccentric orbits. In particular, eccentricity could be produced via sudden mass loss and velocity kicks during supernova explosions, or Lidov-Kozai (LK) oscillations in hierarchical triple systems, or, more generally, secular evolution in multiple-star systems. However, current analytic models of eccentric mass transfer are faced with the problem that they are only well defined in the limit of very high eccentricities, and break down for less eccentric and circular orbits. This provides a major obstacle to implementing such models in binary and higher-order population synthesis codes, which are useful tools for studying the long-term evolution of a large number of systems. Here, we present a new analytic model to describe the secular orbital evolution of binaries undergoing conservative mass transfer. The main improvement of our model is that the mass transfer rate is a smoothly varying function of orbital phase, rather than a delta function centered at periapsis. Consequently, our model is in principle valid for any eccentricity, thereby overcoming the main limitation of previous works. We implement our model in an easy-to-use and publicly available code that can be used as a basis for implementations of our model into population synthesis codes. We investigate the implications of our model in a number of applications with circular and eccentric binaries, and triples undergoing LK oscillations.

The (yet-to-be confirmed) discovery of a Neptune-sized moon around the ~3.2 Jupiter-mass planet in Kepler 1625 puts interesting constraints on the formation of the system. In particular, the relatively wide orbit of the moon around the planet, at ~40 planetary radii, is hard to reconcile with planet formation theories. We demonstrate that the observed characteristics of the system can be explained from the tidal capture of a secondary planet in the young system. After a quick phase of tidal circularization, the lunar orbit, initially much tighter than 40 planetary radii, subsequently gradually widened due to tidal synchronization of the spin of the planet with the orbit, resulting in a synchronous planet-moon system. Interestingly, in our scenario the captured object was originally a Neptune-like planet, turned into a moon by its capture.

The distribution of the inner orbital periods of solar-type main-sequence (MS) triple star systems is known to be peaked at a few days, and this has been attributed to tidal evolution combined with eccentricity excitation due to Lidov-Kozai oscillations. Solar-type MS quadruple star systems also show peaks in their inner orbital period distributions at a few days. Here, we investigate the natural question whether tidal evolution combined with secular evolution can explain the observed inner orbital period distributions in quadruple stars. We carry out population synthesis simulations of solar-type MS quadruple star systems in both the 2+2 (two binaries orbiting each other's barycentre) and 3+1 (triple orbited by a fourth star) configurations. We take into account secular gravitational and tidal evolution, and the effects of passing stars. We assume that no short-period systems are formed initially, and truncate the initial orbital period distributions below 10 d accordingly. We find that, due to secular and tidal evolution, the inner orbital period distributions develop tails at short periods. Although qualitatively consistent with the observations, we find that our simulated orbital period distributions only quantitatively agree with the observations for the 3+1 systems. The observed 2+2 systems, on the other hand, show an enhancement of systems around 10 d, which is not reproduced in the simulations. This suggests that the inner orbital periods of 2+2 systems are not predominantly shaped by tidal and secular evolution, but by other processes, most likely occurring during the stellar formation and early evolution.

Binaries within the sphere of influence of a massive black hole (MBH) in galactic nuclei are susceptible to the Lidov-Kozai (LK) mechanism, which can drive orbits to high eccentricities and trigger strong interactions within the binary such as the emission of gravitational waves (GWs), and mergers of compact objects. These events are potential sources for GW detectors such as Advanced LIGO and VIRGO. The LK mechanism is only effective if the binary is highly inclined with respect to its orbit around the MBH (within a few degrees of 90 deg), implying low rates. However, close to an MBH, torques from the stellar cluster give rise to the process of vector resonant relaxation (VRR). VRR can bring a low-inclination binary into an `active' LK regime in which high eccentricities and strong interactions are triggered in the binary. Here, we study the coupled LK-VRR dynamics, with implications for LIGO and VIRGO GW sources. We carry out Monte Carlo simulations and find that the merger fraction enhancement due to LK-VRR dynamics is up to a factor of ~10 for the lower end of assumed MBH masses (M_MBH = 1e4 MSun), and decreases sharply with increasing M_MBH. We find that, even in our most optimistic scenario, the baseline BH-BH merger rate is small, and the enhancement by LK-VRR coupling is not large enough to increase the rate to well above the LIGO/VIRGO lower limit, 12 Gpc^{-3} yr^{-1}. For the Galactic Center, the LK-VRR-enhanced rate is ~100 times lower than the LIGO/VIRGO limit, and for M_MBH = 1e4 MSun, the rate barely reaches 12 Gpc^{-3} yr^{-1}.

The Kepler mission has detected a number of transiting circumbinary planets (CBPs). Although currently not detected, exomoons could be orbiting some of these CBPs, and they might be suitable for harboring life. A necessary condition for the existence of such exomoons is their long-term dynamical stability. Here, we investigate the stability of exomoons around the Kepler CBPs using numerical $N$-body integrations. We determine regions of stability and obtain stability maps in the (a_m,i_pm) plane, where a_m is the initial exolunar semimajor axis with respect to the CBP, and i_pm is the initial inclination of the orbit of the exomoon around the planet with respect to the orbit of the planet around the stellar binary. Ignoring any dependence on i_pm, for most Kepler CBPs the stability regions are well described by the location of the 1:1 mean motion commensurability of the binary orbit with the orbit of the moon around the CBP. This is related to a destabilizing effect of the binary compared to the case if the binary were replaced by a single body, and which is borne out by corresponding 3-body integrations. For high inclinations, the evolution is dominated by Lidov-Kozai oscillations, which can bring moons in dynamically stable orbits to close proximity within the CBP, triggering strong interactions such as tidal evolution, tidal disruption, or direct collisions. This suggests that there is a dearth of highly-inclined exomoons around the Kepler CBPs, whereas coplanar exomoons are dynamically allowed.

We consider the evolution of stellar hierarchical quadruple systems in the 2+2 (two binaries orbiting each other's barycentre) and 3+1 (triple orbited by a fourth star) configurations. In our simulations, we take into account the effects of secular dynamical evolution, stellar evolution, tidal evolution and encounters with passing stars. We focus on type Ia supernovae (SNe Ia) driven by collisions of carbon-oxygen (CO) white dwarfs (WDs). Such collisions can arise from several channels: (1) collisions due to extremely high eccentricities induced by secular evolution, (2) collisions following a dynamical instability of the system, and (3) collisions driven by semisecular evolution. The systems considered here have initially wide inner orbits, with initial semilatus recti larger than 12 AU, implying no interaction if the orbits were isolated. However, taking into account dynamical evolution, we find that $\approx 0.4$ ($\approx 0.6$) of 2+2 (3+1) systems interact. In particular, Roche Lobe overflow can be triggered possibly in highly eccentric orbits, dynamical instability can ensue due to mass-loss-driven orbital expansion or secular evolution, or a semisecular regime can be entered. We compute the delay-time distributions (DTDs) of collision-induced SNe Ia, and find that they are flatter compared to the observed DTD. Moreover, our combined SNe Ia rates are $(3.7\pm0.7) \times 10^{-6} \,\mathrm{M}_\odot^{-1}$ and $(1.3\pm0.2) \times 10^{-6} \,\mathrm{M}_\odot^{-1}$ for 2+2 and 3+1 systems, respectively, three orders of magnitude lower compared to the observed rate, of order $10^{-3} \, \mathrm{M}_\odot^{-1}$. The low rates can be ascribed to interactions before the stars evolve to CO WDs. However, our results are lower limits given that we considered a subset of quadruple systems.

We extend the formalism of a previous paper to include the effects of flybys and instantaneous perturbations such as supernovae on the long-term secular evolution of hierarchical multiple systems with an arbitrary number of bodies and hierarchy, provided that the system is composed of nested binary orbits. To model secular encounters, we expand the Hamiltonian in terms of the ratio of the separation of the perturber with respect to the barycentre of the multiple system, to the separation of the widest orbit. Subsequently, we integrate over the perturber orbit numerically or analytically. We verify our method for secular encounters, and illustrate it with an example. Furthermore, we describe a method to compute instantaneous orbital changes to multiple systems, such as asymmetric supernovae and impulsive encounters. The secular code, with implementation of the extensions described in this paper, is publicly available within AMUSE, and we provide a number of simple example scripts to illustrate its usage for secular and impulsive encounters, and asymmetric supernovae. The extensions presented in this paper are a next step toward efficiently modeling the evolution of complex multiple systems embedded in star clusters.

Hot Jupiters (HJs) are short-period giant planets that are observed around ~ 1% of solar-type field stars. One possible formation scenario for HJs is high-eccentricity (high-e) migration, in which the planet forms at much larger radii, is excited to high eccentricity by some mechanism, and migrates to its current orbit due to tidal dissipation occurring near periapsis. We consider high-e migration in dense stellar systems such as the cores of globular clusters (GCs), in which encounters with passing stars can excite planets to the high eccentricities needed to initiate migration. We study this process via Monte Carlo simulations of encounters with a star+planet system including the effects of tidal dissipation, using an efficient regularized restricted three-body code. HJs are produced in our simulations over a significant range of the stellar number density n_*. Assuming the planet is initially on a low-eccentricity orbit with semimajor axis 1 au, for n_* < 1e3 pc^{-3} the encounter rate is too low to induce orbital migration, whereas for n_* > 1e6 pc^{-3} HJ formation is suppressed because the planet is more likely ejected from its host star, tidally disrupted, or transferred to a perturbing star. The fraction of planets that are converted to HJs peaks at ~ 2% for intermediate number densities of ~ 4e4 pc^{-3}. Warm Jupiters, giant planets with periods between 10 and 100 days, are produced in our simulations with an efficiency of up to ~ 0.5%. Our results suggest that HJs can form through high-e migration induced by stellar encounters in the centers of of dense GCs, but not in their outskirts where the densities are lower.

Type Ia supernovae (Ia-SNe) are thought to arise from the thermonuclear explosions of white dwarfs (WDs). The progenitors of such explosions are still highly debated; in particular the conditions leading to detonations in WDs are not well understood in most of the suggested progenitor models. Nevertheless, direct head-on collisions of two WDs were shown to give rise to detonations and produce Ia-SNe - like explosions, and were suggested as possible progenitors. The rates of such collisions in dense globular clusters are far below the observed rates of type Ia SNe, but it was suggested that quasi-secular evolution of hierarchical triples could produce a high rate of such collisions. Here we used detailed triple stellar evolution populations synthesis models coupled with dynamical secular evolution to calculate the rates of WD-WD collisions in triples and their properties. We explored a range of models with different realistic initial conditions and derived the expected SNe total mass, mass-ratio and delay time distributions for each of the models. We find that the SNe rate from WD-WD collisions is of the order of 0.1% of the observed Ia-SNe rate across all our models, and the delay-time distribution is almost uniform in time, and is inconsistent with observations. We conclude that SNe from WD-WD collisions in isolated triples can at most provide for a small fraction of Ia-SNe, and can not serve as the main progenitors of such explosions.

Nuclear spiral arms are small-scale transient spiral structures found in the centers of galaxies. Similarly to their galactic-scale counterparts, nuclear spiral arms can perturb the orbits of stars. In the case of the Galactic Center (GC), these perturbations can affect the orbits of stars and binaries in a region extending to several hundred parsecs around the supermassive black hole (MBH), causing diffusion in orbital energy and angular momentum. This diffusion process can drive stars and binaries to close approaches with the MBH, disrupting single stars in tidal disruption events (TDEs), or disrupting binaries, leaving a star tightly bound to the MBH, and an unbound star escaping the galaxy, i.e., a hypervelocity star (HVS). Here, we consider diffusion by nuclear spiral arms in galactic nuclei, specifying to the Milky Way GC. We determine nuclear spiral arm-driven diffusion rates using test-particle integrations, and compute disruption rates. Our TDE rates are up to 20% higher compared to relaxation by single stars. For binaries, the enhancement is up to a factor of ~100, and our rates are comparable to the observed numbers of HVSs and S-stars. Our scenario is complementary to relaxation driven by massive perturbers. In addition, our rates depend on the inclination of the binary with respect to the Galactic plane. Therefore, our scenario provides a novel potential source for the observed anisotropic distribution of HVSs. Nuclear spiral arms may also be important for accelerating the coalescence of binary MBHs, and for supplying nuclear star clusters with stars and gas.

Hierarchical quadruple systems arise naturally in stellar binaries and triples that harbour planets. Examples are hot Jupiters (HJs) in stellar triple systems, and planetary companions to HJs in stellar binaries. The secular dynamical evolution of these systems is generally complex, with secular chaotic motion possible in certain parameter regimes. The latter can lead to extremely high eccentricities and, therefore, strong interactions such as efficient tidal evolution. These interactions are believed to play an important role in the formation of HJs through high-eccentricity migration. Nevertheless, a deeper understanding of the secular dynamics of these systems is still lacking. Here we study in detail the secular dynamics of a special case of hierarchical quadruple systems in either the `2+2' or `3+1' configurations. We show how the equations of motion can be cast in a form representing a perturbed hierarchical three-body system, in which the outer orbital angular momentum vector is precessing steadily around a fixed axis. In this case, we show that eccentricity excitation can be significantly enhanced when the precession period is comparable to the Lidov-Kozai oscillation time-scale of the inner orbit. This arises from an induced large mutual inclination between the inner and outer orbits driven by the precession of the outer orbit, even if the initial mutual inclination is small. We present a simplified semi-analytic model that describes the latter phenomenon.

We consider the formation of binary black hole mergers through the evolution of field massive triple stars. In this scenario, favorable conditions for the inspiral of a black hole binary are initiated by its gravitational interaction with a distant companion, rather than by a common-envelope phase invoked in standard binary evolution models. We use a code that follows self-consistently the evolution of massive triple stars, combining the secular triple dynamics (Lidov-Kozai cycles) with stellar evolution. After a black hole triple is formed, its dynamical evolution is computed using either the orbit-averaged equations of motion, or a high-precision direct integrator for triples with weaker hierarchies for which the secular perturbation theory breaks down. Most black hole mergers in our models are produced in the latter non-secular dynamical regime. We derive the properties of the merging binaries and compute a black hole merger rate in the range (0.3- 1.3) Gpc^{-3}yr^{-1}, or up to ~2.5Gpc^{-3}yr^{-1} if the black hole orbital planes have initially random orientation. Finally, we show that black hole mergers from the triple channel have significantly higher eccentricities than those formed through the evolution of massive binaries or in dense star clusters. Measured eccentricities could therefore be used to uniquely identify binary mergers formed through the evolution of triple stars. While our results suggest up to ~10 detections per year with Advanced-LIGO, the high eccentricities could render the merging binaries harder to detect with planned space based interferometers such as LISA.

Searches for stellar companions to hot Jupiters (HJs) have revealed that planetary systems hosting a HJ are approximately three times more likely to have a stellar companion with a semimajor axis between 50 and 2000 AU, compared to field stars. This correlation suggests that HJ formation is affected by the stellar binary companion. A potential model is high-eccentricity migration, in which the binary companion induces high-eccentricity Lidov-Kozai (LK) oscillations in the proto-HJ orbit, triggering orbital migration driven by tides. A pitfall of this `binary-LK' model is that the observed stellar binaries hosting HJs are typically too wide to produce HJs in sufficient numbers, because of suppression by short-range forces. We propose a modification to the binary-LK model in which there is a second giant planet orbiting the proto-HJ at a semimajor axis of several tens of AU. Such companions are currently hidden to observations, but their presence could be manifested by a propagation of the perturbation of the stellar binary companion inwards to the proto-HJ, thereby overcoming the barrier imposed by short-range forces. Our model does not require the planetary companion orbit to be eccentric and/or inclined with respect to the proto-HJ, but its semimajor axis should lie in a specific range given the planetary mass and binary semimajor axis, and the inclination with respect to the binary should be near $40^\circ$ or $140^\circ$. Our prediction for planetary companions to HJs in stellar binaries should be testable by future observations.

Hot Jupiters (HJs) are Jupiter-like planets orbiting their host star in tight orbits of a few days. They are commonly believed not to have formed in situ, requiring inwards migration towards the host star. One of the proposed migration scenarios is secular high-eccentricity or high-$e$ migration, in which the orbit of the planet is perturbed to high eccentricity by secular processes, triggering strong tidal evolution and orbital migration. Previous theoretical studies have considered secular excitation in stellar binaries. Recently, a number of HJs have been observed in stellar triple systems. In the latter, the secular dynamics are much richer compared to stellar binaries, and HJs could potentially be formed more efficiently. Here, we investigate this possibility by modeling the secular dynamical and tidal evolution of planets in two hierarchical configurations in stellar triple systems. We find that the HJ formation efficiency is higher compared to stellar binaries, but only by at most a few tens of per cent. The orbital properties of the HJs formed in the simulations are very similar to HJs formed in stellar binaries, and similarly to studies of the latter we find no significant number of warm Jupiters. HJs are only formed in our simulations for triples with specific orbital configurations, and our constraints are approximately consistent with current observations. In future, this allows to rule out high-$e$ migration in stellar triples if a HJ is detected in a triple grossly violating these constraints.

Field stars are frequently formed in pairs, and many of these binaries are part of triples or even higher-order systems. Even though, the principles of single stellar evolution and binary evolution, have been accepted for a long time, the long-term evolution of stellar triples is poorly understood. The presence of a third star in an orbit around a binary system can significantly alter the evolution of those stars and the binary system. The rich dynamical behavior in three-body systems can give rise to Lidov-Kozai cycles, in which the eccentricity of the inner orbit and the inclination between the inner and outer orbit vary periodically. In turn, this can lead to an enhancement of tidal effects (tidal friction), gravitational-wave emission and stellar interactions such as mass transfer and collisions. The lack of a self-consistent treatment of triple evolution, including both three-body dynamics as well as stellar evolution, hinders the systematic study and general understanding of the long-term evolution of triple systems. In this paper, we aim to address some of these hiatus, by discussing the dominant physical processes of hierarchical triple evolution, and presenting heuristic recipes for these processes. To improve our understanding on hierarchical stellar triples, these descriptions are implemented in a public source code TrES which combines three-body dynamics (based on the secular approach) with stellar evolution and their mutual influences. Note that modeling through a phase of stable mass transfer in an eccentric orbit is currently not implemented in TrES , but can be implemented with the appropriate methodology at a later stage.

Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with three to five planets. In this variant of high-$e$ migration, the eccentricity of the orbit of the innermost planet is excited on secular time-scales, triggering orbital migration due to tidal dissipation. We use a secular dynamics code and carry out a population synthesis study. We find that HJs are only produced if the viscous time-scale is short ($\approx 0.014$ yr). In contrast, in up to $\approx 0.3$ of systems, the innermost planet is tidally disrupted. The orbital period distribution is peaked around 5 d, consistent with observations. The median HJ mass is $1\,M_\mathrm{J}$ with a maximum of $\approx 2 \, M_\mathrm{J}$, similar to observed HJs. Approximately 0.1 of the HJs have retrograde orbits with respect to the stellar spin. We do not find any warm Jupiters in our simulations, i.e. planets with semimajor axes between 0.1 and 1 AU.

Gas giants orbiting their host star within the ice line are thought to have migrated to their current locations from farther out. Here we consider the origin and dynamical evolution of observed Jupiters, focusing on hot and warm Jupiters with outer friends. We show that the majority of the observed Jupiter pairs (twenty out of twenty-four) will be dynamically unstable if the inner planet was placed at >~1AU distance from the stellar host. This finding is at odds with formation theories that invoke the migration of such planets from semi-major axes >~1AU due to secular dynamical processes (e.g., secular chaos, Lidov-Kozai oscillations) coupled with tidal dissipation. In fact, the results of N-body integrations show that the evolution of dynamically unstable systems does not lead to tidal migration but rather to planet ejections and collisions with the host star. This and other arguments lead us to suggest that most of the observed planets with a companion could not have been transported from further out through secular migration processes. More generally, by using a combination of numerical and analytic techniques we show that the high-e Lidov-Kozai migration scenario can only account for less than 10% of all gas giants observed between 0.1-1 AU. Simulations of multi-planet systems support this result. Our study indicates that rather than starting on highly eccentric orbits with orbital periods above one year, these "warm" Jupiters are more likely to have reached the region where they are observed today without having experienced significant tidal dissipation.

Approximately $0.2 \pm 0.2$ of white dwarfs (WDs) show signs of pollution by metals, which is likely due to the accretion of tidally disrupted planetary material. Models invoking planet-planet interactions after WD formation generally cannot explain pollution at cooling times of several Gyr. We consider a scenario in which a planet is perturbed by Lidov-Kozai oscillations induced by a binary companion and exacerbated by stellar mass loss, explaining pollution at long cooling times. Our computed accretion rates are consistent with observations assuming planetary masses between $\sim 0.01$ and $1\,M_\mathrm{Mars}$, although nongravitational effects may already be important for masses $\lesssim 0.3 \, M_\mathrm{Mars}$. The fraction of polluted WDs in our simulations, $\sim 0.05$, is consistent with observations of WDs with intermediate cooling times between $\sim 0.1$ and 1 Gyr. For cooling times $\lesssim 0.1$ Gyr and $\gtrsim 1$ Gyr, our scenario cannot explain the high observed pollution fractions of up to 0.7. Nevertheless, our results motivate searches for companions around polluted WDs.

We study the secular gravitational dynamics of quadruple systems consisting of a hierarchical triple system orbited by a fourth body. These systems can be decomposed into three binary systems with increasing semimajor axes, binaries A, B and C. The Hamiltonian of the system is expanded in ratios of the three binary separations, and orbit-averaged. Subsequently, we numerically solve the equations of motion. We study highly hierarchical systems that are well described by the lowest-order terms in the Hamiltonian. We find that the qualitative behaviour is determined by the ratio $\mathcal{R}_0$ of the initial Kozai-Lidov (KL) time-scales of the binary pairs AB and BC. If $\mathcal{R}_0\ll 1$, binaries AB remain coplanar if this is initially the case, and KL eccentricity oscillations in binary B are efficiently quenched. If $\mathcal{R}_0\gg 1$, binaries AB become inclined, even if initially coplanar. However, there are no induced KL eccentricity oscillations in binary A. Lastly, if $\mathcal{R}_0\sim 1$, complex KL eccentricity oscillations can occur in binary A that are coupled with the KL eccentricity oscillations in B. Even if binaries A and B are initially coplanar, the induced inclination can result in very high eccentricity oscillations in binary A. These extreme eccentricities could have significant implications for strong interactions such as tidal interactions, gravitational wave dissipation, and collisions and mergers of stars and compact objects. As an example, we apply our results to a planet+moon system orbiting a central star, which in turn is orbited by a distant and inclined stellar companion or planet, and to observed stellar quadruples.

Transiting circumbinary planets are more easily detected around short-period than long-period binaries, but none have yet been observed by {\it Kepler} orbiting binaries with periods shorter than seven days. In triple systems, secular Kozai-Lidov cycles and tidal friction (KLCTF) have been shown to reduce the inner orbital period from $\sim 10^4$ to a few days. Indeed, the majority of short-period binaries are observed to possess a third stellar companion. Using secular evolution analysis and population synthesis, we show that KLCTF makes it unlikely for circumbinary transiting planets to exist around short-period binaries. We find the following outcomes. (1) Sufficiently massive planets in tight and/or coplanar orbits around the inner binary can quench the KL evolution because they induce precession in the inner binary. The KLCTF process does not take place, preventing the formation of a short-period binary. (2) Secular evolution is not quenched and it drives the planetary orbit into a high eccentricity, giving rise to an unstable configuration, in which the planet is most likely ejected from the system. (3) Secular evolution is not quenched but the planet survives the KLCTF evolution. Its orbit is likely to be much wider than the currently observed inner binary orbit, and is likely to be eccentric and inclined with respect to the inner binary. These outcomes lead to two main conclusions: (1) it is unlikely to find a massive planet on a tight and coplanar orbit around a short-period binary, and (2) the properties of circumbinary planets in short-period binaries are constrained by secular evolution.

We present a method for studying the secular gravitational dynamics of hierarchical multiple systems consisting of nested binaries, which is valid for an arbitrary number of bodies and arbitrary hierarchical structure. We derive the Hamiltonian of the system and expand it in terms of the -- assumed to be -- small ratios $x_i$ of binary separations. At the lowest nontrivial expansion order (quadrupole order, second order in $x_i$), the Hamiltonian consists of terms which, individually, depend on binary pairs. At higher orders, in addition to terms depending on binary pairs, we also find terms which, individually, depend on more than two binaries. In general, at order $n$ in $x_i$, individual terms depend on at most $n-1$ binaries. We explicitly derive the Hamiltonian including all terms up and including third order in $x_i$ (octupole order), and including the binary pairwise terms up and including fifth order in $x_i$. These terms are orbit averaged, and we present a new algorithm for efficiently solving the equations of motion. This algorithm is highly suitable for studying the secular evolution of hierarchical systems with complex hierarchies, making long-term integrations of such systems feasible. We show that accurate results are obtained for multiplanet systems with semimajor axis ratios as large as $\approx 0.4$, provided that high-order terms are included. In addition to multiplanet systems with a single star, we apply our results to multistar systems with multiple planets.

We study mass transfers between debris discs during stellar encounters. We carried out numerical simulations of close flybys of two stars, one of which has a disc of planetesimals represented by test particles. We explored the parameter space of the encounters, varying the mass ratio of the two stars, their pericentre and eccentricity of the encounter, and its geometry. We find that particles are transferred to the other star from a restricted radial range in the disc and the limiting radii of this transfer region depend on the parameters of the encounter. We derive an approximate analytic description of the inner radius of the region. The efficiency of the mass transfer generally decreases with increasing encounter pericentre and increasing mass of the star initially possessing the disc. Depending on the parameters of the encounter, the transfer particles have a specific distributions in the space of orbital elements (semimajor axis, eccentricity, inclination, and argument of pericentre) around their new host star. The population of the transferred particles can be used to constrain the encounter through which it was delivered. We expect that many stars experienced transfer among their debris discs and planetary systems in their birth environment. This mechanism presents a formation channel for objects on wide orbits of arbitrary inclinations, typically having high eccentricity but possibly also close-to-circular (eccentricities of about 0.1). Depending on the geometry, such orbital elements can be distinct from those of the objects formed around the star.

Flares in X-ray and near infrared are observed above the quiescent emission of the supermassive black hole (SBH) in the Galactic Centre (GC) at a rate of approximately once per day. One proposed energy source for these flares is the tidal disruption of planetesimals with radius $\gtrsim 10$ km passing within $\sim$ 1 AU of the SBH. Very little is known about the formation and evolution of planetesimals in galactic nuclei such as the GC, making predictions for flaring event rates uncertain. We explore two scenarios for the formation of planetesimals in the GC: (1) in a large-scale cloud bound to the SBH, and (2) in debris discs around stars. We model their orbital evolution around the SBH using the Fokker-Planck equation and investigate the effect of gravitational interactions with various relevant perturbers. Our predicted flaring rate, $\approx 0.6 \, \mathrm{day^{-1}}$, is nearly independent of the distribution of perturbers. Moreover, it is insensitive to scenarios (1) or (2). The assumed number of planetesimals per star is consistent with debris discs around stars in the Solar neighbourhood. In scenario (1) this implies that the number of planetesimals formed in the large-scale cloud is strongly correlated with the number of stars, and this requires finetuning for our results to be consistent with the observed flaring rate. We favour the alternative explanation that planetesimals in the GC are formed in debris discs around stars, similar to the Solar neighbourhood.

General relativistic precession limits the ability of gravitational encounters to increase the eccentricity $e$ of orbits near a supermassive black hole (SBH). This "Schwarzschild barrier" (SB) has been shown to play an important role in the orbital evolution of stars like the galactic center S-stars. However, the evolution of orbits below the SB, $e>e_\mathrm{SB}$, is not well understood; the main current limitation is the computational complexity of detailed simulations. Here we present an $N$-body algorithm that allows us to efficiently integrate orbits of test stars around a SBH including general relativistic corrections to the equations of motion and interactions with a large ($\gtrsim 10^3$) number of field stars. We apply our algorithm to the S-stars and extract diffusion coefficients describing the evolution in angular momentum $L$. We identify three angular momentum regimes, in which the diffusion coefficients depend in functionally different ways on $L$. Regimes of lowest and highest $L$ are well-described in terms of non-resonant relaxation (NRR) and resonant relaxation (RR), respectively. In addition, we find a new regime of "anomalous relaxation" (AR). We present analytic expressions, in terms of physical parameters, that describe the diffusion coefficients in all three regimes, and propose a new, empirical criterion for the location of the SB in terms of the $L$-dependence of the diffusion coefficients. Subsequently we apply our results to obtain the steady-state distribution of angular momentum for orbits near a SBH.

Hierarchical triple systems are common among field stars yet their long-term evolution is poorly understood theoretically. In such systems Kozai cycles can be induced in the inner binary system during which the inner orbit eccentricity and the inclination between both binary orbits vary periodically. These cycles, combined with tidal friction and gravitational wave emission, can significantly affect the inner binary evolution. To investigate these effects quantitatively we perform a population synthesis study of triple systems and focus on evolutionary paths that lead to mergers of carbon-oxygen (CO) white dwarfs (WDs), which constitute an important candidate progenitor channel for type Ia supernovae (SNe Ia). We approach this problem by Monte Carlo sampling from observation-based distributions of systems within the primary mass range 1.0 - 6.5 M_Sun and inner orbit semi-major axes a_1 and eccentricities e_1 satisfying a_1 (1-e_1^2) > 12 AU, i.e. non-interacting in the absence of a tertiary component. We evolve these systems by means of a newly developed algorithm that couples secular triple dynamics with an existing binary population synthesis code. We find that the tertiary significantly alters the inner binary evolution in about 24% of all sampled systems. In particular, we find several channels leading to CO WD mergers. Amongst these is a novel channel in which a collision occurs in wide inner binary systems as a result of extremely high eccentricities induced by Kozai cycles. With assumptions on which CO WD mergers lead to a SN Ia explosion we estimate the SNe Ia delay time distribution resulting from triples and compare to a binary population synthesis study and to observations. Although we find that our triple rate is low, we have determined a lower limit of the triple-induced SNe Ia rate and further study is needed that includes triples with initially tighter inner orbits.