BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Adrian S. Hamers

Adrian S. Hamers appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPastro-ph.HEastro-ph.GAastro-ph.COastro-ph.IMgr-qcMachine Learningnlin.CDphysics.comp-ph

Catalog footprint

What is connected

22works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

22 published item(s)

preprint2022arXiv

A statistical view on stable and unstable Roche lobe overflow of a tertiary star onto the inner binary in triple systems

In compact stellar triple systems, an evolved tertiary star can overflow its Roche lobe around the inner binary. Subsequently, the tertiary star can transfer mass to the inner binary in a stable manner, or Roche lobe overflow (RLOF) can be unstable and lead to common-envelope (CE) evolution. In the latter case, the inner binary enters the extended envelope of the tertiary star and spirals in towards the donor's core, potentially leading to mergers or ejections. Although studied in detail for individual systems, a comprehensive statistical view on the various outcomes of triple RLOF is lacking. Here, we carry out 10^5 population synthesis simulations of tight triples, self-consistently taking into account stellar evolution, binary interactions, and gravitational dynamics. Also included are prescriptions for the long-term evolution of stable triple mass transfer, and triple CE evolution. Although simple and ignoring hydrodynamic effects, these prescriptions allow for a qualitative statistical study. We find that triple RLOF occurs in ~0.06% of triples, with ~64% leading to stable mass transfer, and ~36% to triple CE evolution. Triple CE is most often (~76%) followed by one or multiple mergers in short succession, most likely an inner binary merger of two main-sequence stars. Other outcomes of triple CE are a binary+single system (~23%, most of which not involving exchange interactions), and a stable triple (~1%). We also estimate the rate of Type Ia supernovae involving white dwarf mergers following triple RLOF, but find only a negligible contribution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Algebraic and machine learning approach to hierarchical triple-star stability

We present two approaches to determine the dynamical stability of a hierarchical triple-star system. The first is an improvement on the Mardling-Aarseth stability formula from 2001, where we introduce a dependence on inner orbital eccentricity and improve the dependence on mutual orbital inclination. The second involves a machine learning approach, where we use a multilayer perceptron (MLP) to classify triple-star systems as `stable' and `unstable'. To achieve this, we generate a large training data set of 10^6 hierarchical triples using the N-body code MSTAR. Both our approaches perform better than previous stability criteria, with the MLP model performing the best. The improved stability formula and the machine learning model have overall classification accuracies of 93 % and 95 % respectively. Our MLP model, which accurately predicts the stability of any hierarchical triple-star system within the parameter ranges studied with almost no computation required, is publicly available on Github in the form of an easy-to-use Python script.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How important is secular evolution for black hole and neutron star mergers in 2+2 and 3+1 quadruple-star systems?

Mergers of black holes (BHs) and neutron stars (NSs) result in the emission of gravitational waves that can be detected by LIGO. In this paper, we look at 2+2 and 3+1 quadruple-star systems, which are common among massive stars, the progenitors of BHs and NSs. We carry out a detailed population synthesis of quadruple systems using the MSE code, which seamlessly takes into consideration stellar evolution, binary and tertiary interactions, $N$-body dynamics, and secular evolution. We find that, although secular evolution plays a role in compact object (BH and NS) mergers, (70--85) \% (depending on the model assumptions) of the mergers are solely due to common envelope (CE) evolution. Significant eccentricities in the LIGO band (higher than 0.01) are only obtained with zero supernova (SNe) kicks and are directly linked to the role of secular evolution. A similar outlier effect is seen in the $χ_{\mathrm{eff}}$ distribution, with negative values obtained only with zero SNe kicks. When kicks are taken into account, there are no systems that evolve into a quadruple consisting of four compact objects. For our fiducial model, we estimate the merger rates (in units of $\Gpcyr$) in 2+2 quadruples (3+1 quadruples) to be 10.8 $\pm$ 0.9 (2.9 $\pm$ 0.5), 5.7 $\pm$ 0.6 (1.4 $\pm$ 0.4) and 0.6 $\pm$ 0.2 (0.7 $\pm$ 0.3) for BH-BH, BH-NS and NS-NS mergers respectively. The BH-BH merger rates represent a significant fraction of the current LIGO rates, whereas the other merger rates fall short of LIGO estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Equilibrium Tide: An Updated Prescription for Population Synthesis Codes

We present an updated prescription for the equilibrium tides suitable for population synthesis codes. A grid of 1D evolutionary models was created and the viscous time-scale was calculated for each detailed model. A metallicity dependent power-law relation was fitted to both the convective cores and convective envelopes of the models. The prescription was implemented into the population synthesis code BSE and predicts an 16.5% reduction in the overall number of merges, with those involving main-sequence stars most affected. The new prescription also reduces the overall supernova rate by 3.6% with individual channels being differently affected. The single degenerate Ia supernova occurrence is reduced by 12.8%. The merging of two Carbon Oxygen white dwarfs to cause a Ia supernova occurs 16% less frequently. The number of sub-synchronously rotating stars in close binaries is substantially increased with our prescription, as is the number of non-circularized systems at the start of common-envelope evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Chaos in self-gravitating many-body systems: Lyapunov time dependence of $N$ and the influence of general relativity

In self-gravitating $N$-body systems, small perturbations introduced at the start, or infinitesimal errors that are produced by the numerical integrator or are due to limited precision in the computer, grow exponentially with time. For Newton's gravity, we confirm earlier results that for relatively homogeneous systems, this rate of growth per crossing time increases with $N$ up to $N \sim 30$, but that for larger systems, the growth rate has a weaker scaling with $N$. For concentrated systems, however, the rate of exponential growth continues to scale with $N$. In relativistic self-gravitating systems, the rate of growth is almost independent of $N$. This effect, however, is only noticeable when the system's mean velocity approaches the speed of light to within three orders of magnitude. The chaotic behavior of systems with more than a dozen bodies for the usually adopted approximation of only solving the pairwise interactions in the Einstein-Infeld-Hoffmann equation of motion is qualitatively different than when the interaction terms (or cross terms) are taken into account. This result provides a strong motivation for follow-up studies on the microscopic effect of general relativity on orbital chaos, and on the influence of higher-order cross-terms in the Taylor-series expansion of the Einstein-Infeld-Hoffmann equations of motion.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multiple Stellar Evolution: a population synthesis algorithm to model the stellar, binary, and dynamical evolution of multiple-star systems

In recent years, observations have shown that multiple-star systems such as hierarchical triple and quadruple-star systems are common, especially among massive stars. They are potential sources of interesting astrophysical phenomena such as compact object mergers, leading to supernovae, and gravitational wave events. However, many uncertainties remain in their often complex evolution. Here, we present the population synthesis code Multiple Stellar Evolution (MSE), designed to rapidly model the stellar, binary, and dynamical evolution of multiple-star systems. MSE includes a number of new features not present in previous population synthesis codes: (1) an arbitrary number of stars, as long as the initial system is hierarchical, (2) dynamic switching between secular and direct N-body integration for efficient computation of the gravitational dynamics, (3) treatment of mass transfer in eccentric orbits, which occurs commonly in multiple-star systems, (4) a simple treatment of tidal, common-envelope, and mass transfer evolution in which the accretor is a binary instead of a single star, (5) taking into account planets within the stellar system, and (6) including gravitational perturbations from passing field stars. MSE, written primarily in the C++ language, will be made publicly available and has few prerequisites; a convenient Python interface is provided. We give a detailed description of MSE and illustrate how to use the code in practice. We demonstrate its operation in a number of examples.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A census of main-sequence interactions in the Multiple Star Catalog

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Binary-binary scattering in the secular limit

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Secular dynamics of hierarchical multiple systems composed of nested binaries, with an arbitrary number of bodies and arbitrary hierarchical structure. III. Suborbital effects: hybrid integration techniques and orbit-averaging corrections

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Was GW190412 born from a hierarchical 3+1 quadruple configuration?

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 χ_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 χ_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 χ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Formation and merging of Mass Gap Black Holes in Gravitational Wave Merger Events from Wide Hierarchical Quadruple Systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A triple origin for the lack of tight coplanar circumbinary planets around short-period binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical constraints on the origin of hot and warm Jupiters with close friends

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mass transfer between debris discs during close stellar encounters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Secular dynamics of hierarchical multiple systems composed of nested binaries, with an arbitrary number of bodies and arbitrary hierarchical structure. First applications to multiplanet and multistar systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Secular dynamics of hierarchical quadruple systems: the case of a triple system orbited by a fourth body

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Secular dynamics of multiplanet systems: implications for the formation of hot and warm Jupiters via high-eccentricity migration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

White dwarf pollution by planets in stellar binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probing the formation of planetesimals in the Galactic Centre using Sgr A* flares

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic dynamics of stars near a supermassive black hole

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Population synthesis of triple systems in the context of mergers of carbon-oxygen white dwarfs

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.

0 citations0 reviewsNo code linkedOpen access