Showing 1–50 of 52 results for author: Hamers, A S

Quadruple-star systems are not always nested triples: a machine learning approach to dynamical stability

Authors: Pavan Vynatheya, Rosemary A. Mardling, Adrian S. Hamers

Abstract: The dynamical stability of quadruple-star systems has traditionally been treated as a problem involving two `nested' triples which constitute a quadruple. In this novel study, we employed a machine learning algorithm, the multi-layer perceptron (MLP), to directly classify 2+2 and 3+1 quadruples based on their stability (or long-term boundedness). The training data sets for the classification, comp… ▽ More The dynamical stability of quadruple-star systems has traditionally been treated as a problem involving two `nested' triples which constitute a quadruple. In this novel study, we employed a machine learning algorithm, the multi-layer perceptron (MLP), to directly classify 2+2 and 3+1 quadruples based on their stability (or long-term boundedness). The training data sets for the classification, comprised of $5\times10^5$ quadruples each, were integrated using the highly accurate direct $N$-body code MSTAR. We also carried out a limited parameter space study of zero-inclination systems to directly compare quadruples to triples. We found that both our quadruple MLP models perform better than a `nested' triple MLP approach, which is especially significant for 3+1 quadruples. The classification accuracies for the 2+2 MLP and 3+1 MLP models are 94% and 93% respectively, while the scores for the `nested' triple approach are 88% and 66% respectively. This is a crucial implication for quadruple population synthesis studies. Our MLP models, which are very simple and almost instantaneous to implement, are available on GitHub, along with Python3 scripts to access them. △ Less

Submitted 7 August, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: Accepted for publication by MNRAS

Forming Hot Subluminous Stars from Hierarchical Triples -- I. The Role of an Outer Tertiary on Formation Channels

Authors: Holly P. Preece, Adrian S. Hamers, Tiara Battich, Abinaya S. Rajamuthukumar

Abstract: We present evolutionary pathways for creating hot subdwarf OB (sdOB) stars from hierarchical triple configurations. We use the population synthesis code Multiple Stellar Evolution (MSE) to follow the stellar, binary, and gravitational dynamical evolution of triple-star systems. To ascertain the effect of the outer tertiary, we also consider the evolution of the inner binary with the tertiary compo… ▽ More We present evolutionary pathways for creating hot subdwarf OB (sdOB) stars from hierarchical triple configurations. We use the population synthesis code Multiple Stellar Evolution (MSE) to follow the stellar, binary, and gravitational dynamical evolution of triple-star systems. To ascertain the effect of the outer tertiary, we also consider the evolution of the inner binary with the tertiary component removed. We find we are able to create sdOB stars in single, binary and triple configurations. We also demonstrate that it is possible to form sdOBs in systems which undergo triple common envelope evolution, when the tertiary star undergoes unstable mass transfer onto the inner binary. We are unable to create single or wide sdOB systems without involving a merger earlier in the evolution. The triples can produce sdOBs in binaries with wide, non-interacting companions through binary interactions, which is impossible in isolated binaries. Owing to the closeness of the inner binary in hierarchical triples the formation channels associated with stable mass transfer are suppressed when compared to the isolated binary case. △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: 10 Pages, 10 Figures, Accepted for publication

Algebraic and machine learning approach to hierarchical triple-star stability

Authors: Pavan Vynatheya, Adrian S. Hamers, Rosemary A. Mardling, Earl P. Bellinger

Abstract: 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 cla… ▽ More 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. △ Less

Submitted 6 September, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: Accepted for publication by MNRAS

The Equilibrium Tide: An Updated Prescription for Population Synthesis Codes

Authors: Holly P. Preece, Adrian S. Hamers, Patrick G. Neunteufel, Adam L. Schafer, Christopher A. Tout

Abstract: 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… ▽ More 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. △ Less

Submitted 23 June, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 22 pages, 17 Figures, Accepted for publication 13/05/22

A spectroscopic quadruple as a possible progenitor of sub-Chandrasekhar Type Ia supernovae

Authors: T. Merle, A. S. Hamers, S. Van Eck, A. Jorissen, M. Van der Swaelmen, K. Pollard, R. Smiljanic, D. Pourbaix, T. Zwitter, G. Traven, G. Gilmore, S. Randich, A. Gonneau, A. Hourihane, G. Sacco, C. C. Worley

Abstract: Binaries have received much attention as possible progenitors of Type Ia supernova (SNIa) explosions, but long-term gravitational effects in tight triple or quadruple systems could also play a key role in producing SNIa. Here we report on the properties of a spectroscopic quadruple (SB4) found within a star cluster: the 2+2 hierarchical system HD 74438. Its membership in the open cluster IC 2391 m… ▽ More Binaries have received much attention as possible progenitors of Type Ia supernova (SNIa) explosions, but long-term gravitational effects in tight triple or quadruple systems could also play a key role in producing SNIa. Here we report on the properties of a spectroscopic quadruple (SB4) found within a star cluster: the 2+2 hierarchical system HD 74438. Its membership in the open cluster IC 2391 makes it the youngest (43 My) SB4 discovered so far and among the quadruple systems with the shortest outer orbital period. The eccentricity of the 6 y outer period is 0.46 and the two inner orbits, with periods of 20.5 d and 4.4 d, and eccentricities of 0.36 and 0.15, are not coplanar. Using an innovative combination of ground-based high resolution spectroscopy and Gaia/Hipparcos astrometry, we show that this system is undergoing secular interaction that likely pumped the eccentricity of one of the inner orbits higher than expected for the spectral types of its components. We compute the future evolution of HD 74438 and show that this system is an excellent candidate progenitor of sub-Chandrasekhar SNIa through white dwarf (WD) mergers. Taking into account the contribution of this specific type of SNIa better accounts for the chemical evolution of iron-peak elements in the Galaxy than considering only near Chandrasekhar-mass SNIa. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: accepted for publication in Nature Astronomy

Properties and applications of a predicted population of runaway He-sdO/B stars ejected from single degenerate He-donor SNe

Authors: P. Neunteufel, H. Preece, M. Kruckow, S. Geier, A. S. Hamers, S. Justham, Ph. Podsiadlowski

Abstract: This study builds on previous works, producing the most extensive prediction of the properties of such a hypothetical population to date, taking into account both Chandrasekhar and non-Chandrasekhar mass events. These results are then used to define criteria for membership of this population and characterise putative subpopulations This study contains 6x10^6 individual ejection trajectories out… ▽ More This study builds on previous works, producing the most extensive prediction of the properties of such a hypothetical population to date, taking into account both Chandrasekhar and non-Chandrasekhar mass events. These results are then used to define criteria for membership of this population and characterise putative subpopulations This study contains 6x10^6 individual ejection trajectories out of the Galactic plane calculated with the stellar kinematics framework SHyRT, which are analysed with regard to their bulk observational properties. These are then put into context with the only previously identified population member US\,708 and applied to a number of other possible candidate objects. We find that two additional previously observed objects possess properties to warrant a designation as candidate objects. Characterisation of these object with respect to the predicted population finds all of them to be extreme in at least one astrometric observable. We find that current observations support a Galactic SN rate on the order of ~3x10^-7/yr to ~2x10^-6/yr, three orders of magnitude below the inferred Galactic SN Ia rate and two orders of magnitude below the formation rate of predicted He-donor progenitors. The number of currently observed population members suggests that the He-donor scenariois not a dominant contributor to the number of observed SNe Ia. However, we find that, even at the low event rate suggested, the majority of possibly detectable population members is still undetected. The extreme nature of current population members suggests that a still larger number of objects has simply evaded detection up to this point, hinting at a higher contribution than is currently supported by observation. - abridged - △ Less

Submitted 28 March, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 26 pages, 23 figures, 5 tables - accepted in A&A

Journal ref: A&A 663, A91 (2022)

An improved numerical fit to the peak harmonic gravitational wave frequency emitted by an eccentric binary

Authors: Adrian S. Hamers

Abstract: I present a numerical fit to the peak harmonic gravitational wave frequency emitted by an eccentric binary system in the post-Newtonian approximation. This fit significantly improves upon a previous commonly-used fit in population synthesis studies, in particular for eccentricities $\lesssim 0.8$. I present a numerical fit to the peak harmonic gravitational wave frequency emitted by an eccentric binary system in the post-Newtonian approximation. This fit significantly improves upon a previous commonly-used fit in population synthesis studies, in particular for eccentricities $\lesssim 0.8$. △ Less

Submitted 30 November, 2021; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: Research Notes of the American Astronomical Society (RNAAS). 3 pages, 1 figure

Journal ref: 2021, Research Notes of the American Astronomical Society (RNAAS), Volume 5, Number 11 275

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

Authors: Pavan Vynatheya, Adrian S. Hamers

Abstract: 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 evo… ▽ More 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. △ Less

Submitted 10 January, 2022; v1 submitted 27 October, 2021; originally announced October 2021.

Comments: Accepted for publication by ApJ (with minor revisions)

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

Authors: Adrian S. Hamers, Hila Glanz, Patrick Neunteufel

Abstract: 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 towar… ▽ More 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. △ Less

Submitted 5 January, 2022; v1 submitted 30 September, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJS. 23 pages, 13 figures

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

Authors: Simon F. Portegies Zwart, Tjarda C. N. Boekholt, Emiel Por, Adrian S. Hamers, Steve L. W. McMillan

Abstract: 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$, bu… ▽ More 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. △ Less

Submitted 17 November, 2021; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: A&A in press

Journal ref: A&A 659, A86 (2022)

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

Authors: Adrian S. Hamers, Hagai B. Perets, Todd A. Thompson, Patrick Neunteufel

Abstract: 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 co… ▽ More 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. △ Less

Submitted 2 December, 2021; v1 submitted 28 July, 2021; originally announced July 2021.

Comments: Accepted for publication in ApJ. 22 pages, 16 figures

First and second-generation black hole and neutron star mergers in 2+2 quadruples: population statistics

Authors: Adrian S. Hamers, Giacomo Fragione, Patrick Neunteufel, Bence Kocsis

Abstract: Recent detections of gravitational waves from mergers of neutron stars (NSs) and black holes (BHs) in the low and high-end mass gap regimes pose a puzzle to standard stellar and binary evolution theory. Mass-gap mergers may originate from successive mergers in hierarchical systems such as quadruples. Here, we consider repeated mergers of NSs and BHs in stellar 2+2 quadruple systems, in which secul… ▽ More Recent detections of gravitational waves from mergers of neutron stars (NSs) and black holes (BHs) in the low and high-end mass gap regimes pose a puzzle to standard stellar and binary evolution theory. Mass-gap mergers may originate from successive mergers in hierarchical systems such as quadruples. Here, we consider repeated mergers of NSs and BHs in stellar 2+2 quadruple systems, in which secular evolution can accelerate the merger of one of the inner binaries. Subsequently, the merger remnant may interact with the companion binary, yielding a second-generation merger. We model the initial stellar and binary evolution of the inner binaries as isolated systems. In the case of successful compact object formation, we subsequently follow the secular dynamical evolution of the quadruple system. When a merger occurs, we take into account merger recoil, and model subsequent evolution using direct N-body integration. With different assumptions on the initial properties, we find that the majority of first-generation mergers are not much affected by secular evolution, with their observational properties mostly consistent with isolated binaries. A small subset shows imprints of secular evolution through residual eccentricity in the LIGO band, and retrograde spin-orbit orientations. Second-generation mergers are ~10^7 times less common than first-generation mergers, and can be strongly affected by scattering (i.e., three-body interactions) induced by the first-generation merger. In particular, scattering can account for mergers within the low-end mass gap, although not the high-end mass gap. Also, in a few cases, scattering could explain highly eccentric LIGO sources and negative effective spin parameters. △ Less

Submitted 6 August, 2021; v1 submitted 5 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS (updated references). 38 pages, 91 figures. Includes supplementary material appended to the main manuscript

Predicted spatial and velocity distributions of ejected companion stars of helium accretion-induced thermonuclear supernovae

Authors: P. Neunteufel, M. Kruckow, S. Geier, A. S. Hamers

Abstract: Thermonuclear supernovae (SNe), a subset of which are the highly important SNe Type\,Ia, remain one of the more poorly understood phenomena known to modern astrophysics. In recent years, the single degenerate helium (He) donor channel, where a white dwarf star (WD) accretes He-rich matter from a hydrogen-depleted companion, has emerged as a promising candidate progenitor scenario for these events.… ▽ More Thermonuclear supernovae (SNe), a subset of which are the highly important SNe Type\,Ia, remain one of the more poorly understood phenomena known to modern astrophysics. In recent years, the single degenerate helium (He) donor channel, where a white dwarf star (WD) accretes He-rich matter from a hydrogen-depleted companion, has emerged as a promising candidate progenitor scenario for these events. An unresolved question in this scenario is the fate of the companion star, which would be evident as a runaway hot subdwarf (He sdO/B) in the aftermath of the SN event. Previous studies have shown that the kinematic properties of an ejected companion provide an opportunity to closer examination of the properties of an SN progenitor system. However, with the number of observed objects not matching predictions by theory, the viability of this mechanism is called into question. In this study, we first synthesize a population of companion stars ejected by the aforementioned mechanism, taking into account predicted ejection velocities, inferred population density in the Galactic (Gal.) mass distribution and subsequent kinematics in the Gal. potential. We then discuss the astrometric properties of this population. We present $10^{6}$ individual ejection trajectories, numerically computed with a newly developed, lightweight simulation framework. A peak in the density distribution for close objects is expected in the direction of the Gal. center. If the entire considered mass range is realized, the radial velocity distribution should show a peak at 500\kms. If only close US\,708 analogues are considered, there should be a peak at ($\sim750-850$)\kms. We show that the puzzling lack of confirmed surviving companion stars of thermonuclear SNe, though possibly an observation-related selection effect, may indicate a selection against high mass donors in the SD He donor channel. (-abridged-) △ Less

Submitted 20 January, 2021; v1 submitted 30 November, 2020; originally announced November 2020.

Comments: Accepted for publication by A&AL (15/1/2021). Main paper: 5 pages, Appendix: 5 pages, 9 figures

Journal ref: A&A 646, L8 (2021)

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

Authors: Adrian S. Hamers, Antti Rantala, Patrick Neunteufel, Holly Preece, Pavan Vynatheya

Abstract: 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 pre… ▽ More 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. △ Less

Submitted 18 February, 2021; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS (fixed typos and updated references). 35 pages, 8 figures. Code can be found at https://github.com/hamers/mse

Properties of von Zeipel-Lidov-Kozai oscillations in triple systems at the quadrupole order: relaxing the test particle approximation

Authors: Adrian S. Hamers

Abstract: Von Zeipel-Lidov-Kozai (ZLK) oscillations in hierarchical triple systems have important astrophysical implications such as triggering strong interactions and producing, e.g., Type Ia supernovae and gravitational wave sources. When considering analytic properties of ZLK oscillations at the lowest (quadrupole) expansion order, as well as complications due to higher-order terms, one usually assumes t… ▽ More Von Zeipel-Lidov-Kozai (ZLK) oscillations in hierarchical triple systems have important astrophysical implications such as triggering strong interactions and producing, e.g., Type Ia supernovae and gravitational wave sources. When considering analytic properties of ZLK oscillations at the lowest (quadrupole) expansion order, as well as complications due to higher-order terms, one usually assumes the test particle limit, in which one of the bodies in the inner binary is massless. Although this approximation holds well for, e.g., planetary systems, it is less accurate for systems with more comparable masses such as stellar triples. Whereas non-test-particle effects are usually taken into account in numerical simulations, a more analytic approach focusing on the differences between the test particle and general case (at quadrupole order) has, to our knowledge, not been presented. Here, we derive several analytic properties of secular oscillations in triples at the quadruple expansion order. The latter applies even to relatively compact triples, as long as the inner bodies are similar in mass such that octupole-order effects are suppressed. We consider general conditions for the character of the oscillations (circular versus librating), minimum and maximum eccentricities, and timescales, all as a function of $γ\equiv (1/2) L_1/G_2$, a ratio of inner-to-outer orbital angular momenta variables ($γ=0$ in the test particle limit). In particular, eccentricity oscillations are more effective at retrograde orientations for non-zero $γ$; assuming zero initial inner eccentricity, the maximum eccentricity peaks at $\cos(i_\mathrm{rel,0}) = -γ$, where $i_\mathrm{rel,0}$ is the initial relative inclination. We provide a Python script which can be used to quickly compute these properties. △ Less

Submitted 25 November, 2020; v1 submitted 6 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS (reference added in proof in Introduction). 16 pages, 17 figures. Code can be found at https://github.com/hamers/ZLK

WD 1856 b: a close giant planet around a white dwarf that could have survived a common-envelope phase

Authors: F. Lagos, M. R. Schreiber, M. Zorotovic, B. T. Gänsicke, M. P. Ronco, Adrian S. Hamers

Abstract: The discovery of a giant planet candidate orbiting the white dwarf WD 1856+534 with an orbital period of 1.4 d poses the questions of how the planet reached its current position. We here reconstruct the evolutionary history of the system assuming common envelope evolution as the main mechanism that brought the planet to its current position. We find that common envelope evolution can explain the p… ▽ More The discovery of a giant planet candidate orbiting the white dwarf WD 1856+534 with an orbital period of 1.4 d poses the questions of how the planet reached its current position. We here reconstruct the evolutionary history of the system assuming common envelope evolution as the main mechanism that brought the planet to its current position. We find that common envelope evolution can explain the present configuration if it was initiated when the host star was on the AGB, the separation of the planet at the onset of mass transfer was in the range 1.69-2.35 au, and if in addition to the orbital energy of the surviving planet either recombination energy stored in the envelope or another source of additional energy contributed to expelling the envelope. We also discuss the evolution of the planet prior to and following common envelope evolution. Finally, we find that if the system formed through common envelope evolution, its total age is in agreement with its membership to the Galactic thin disc. We therefore conclude that common envelope evolution is at least as likely as alternative formation scenarios previously suggested such as planet-planet scattering or Kozai-Lidov oscillations. △ Less

Submitted 27 November, 2020; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: 7 pages, 3 figures, 2 tables; accepted for publication in MNRAS

The ominous fate of exomoons around hot Jupiters in the high-eccentricity migration scenario

Authors: Alessandro A. Trani, Adrian S. Hamers, Aaron Geller, Mario Spera

Abstract: All the giant planets in the solar system host a large number of natural satellites. Moons in extrasolar systems are difficult to detect, but a Neptune-sized exomoon candidate has been recently found around a Jupiter-sized planet in the Kepler-1625bsystem. Due to their relative ease of detection, hot Jupiters (HJs), which reside in close orbits around their host stars with a period of a few days,… ▽ More All the giant planets in the solar system host a large number of natural satellites. Moons in extrasolar systems are difficult to detect, but a Neptune-sized exomoon candidate has been recently found around a Jupiter-sized planet in the Kepler-1625bsystem. Due to their relative ease of detection, hot Jupiters (HJs), which reside in close orbits around their host stars with a period of a few days, may be very good candidates to search for exomoons. It is still unknown whether the HJ population can host (or may have hosted) exomoons. One suggested formation channel for HJs is high-eccentricity migration induced by a stellar binary companion combined with tidal dissipation. Here, we investigate under which circumstances an exomoon can prevent or allow high-eccentricity migration of a HJ, and in the latter case, if the exomoon can survive the migration process. We use both semianalytic arguments, as well as direct N-body simulations including tidal interactions. Our results show that massive exomoons are efficient at preventing high-eccentricity migration. If an exomoon does instead allow for planetary migration, it is unlikely that the HJ formed can host exomoons since the moon will either spiral onto the planet or escape from it during the migration process. A few escaped exomoons can become stable planets after the Jupiter has migrated, or by tidally migrating themselves. The majority of the exomoons end up being ejected from the system or colliding with the primary star and the host planet. Such collisions might nonetheless leave observable features, such as a debris disc around the primary star or exorings around the close-in giant. △ Less

Submitted 25 October, 2020; v1 submitted 31 August, 2020; originally announced August 2020.

Comments: 12 Pages, 6 Figures, accepted for publication in MNRAS

The impact of distant fly-bys on the rate of binary primordial black hole mergers

Authors: Sam Young, Adrian S. Hamers

Abstract: 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, $f_\mathrm{PBH}\sim 1$, the expected fractional change in coalescence time is negligible, of order… ▽ More 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, $f_\mathrm{PBH}\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, $f_\mathrm{PBH}\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. △ Less

Submitted 27 October, 2020; v1 submitted 26 June, 2020; originally announced June 2020.

Comments: 26 pages, 7 figures. V2: Updated to match the published version

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

Authors: Adrian S. Hamers, Mohammadtaher Safarzadeh

Abstract: 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 channe… ▽ More 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. △ Less

Submitted 9 June, 2020; v1 submitted 6 May, 2020; originally announced May 2020.

Comments: Accepted for publication in ApJ. 14 pages, 12 figures

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

Authors: Adrian S. Hamers

Abstract: 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 numeric… ▽ More 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. △ Less

Submitted 17 April, 2020; originally announced April 2020.

Comments: Accepted for publication in MNRAS. 15 pages, 4 figures. Code can be found at https://github.com/hamers/secularmultiple

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

Authors: Adrian S. Hamers

Abstract: 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. Al… ▽ More 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. △ Less

Submitted 21 April, 2020; v1 submitted 20 February, 2020; originally announced February 2020.

Comments: Accepted for publication in MNRAS. 16 pages, 13 figures

Binary-binary scattering in the secular limit

Authors: Adrian S. Hamers, Johan Samsing

Abstract: 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 app… ▽ More 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. △ Less

Submitted 9 March, 2020; v1 submitted 12 February, 2020; originally announced February 2020.

Comments: Accepted for publication in MNRAS. 19 pages, 9 figures. Code can be found at https://github.com/hamers/flybys_bin

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

Authors: Mohammadtaher Safarzadeh, Adrian S. Hamers, Abraham Loeb, Edo Berger

Abstract: 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 star… ▽ More 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. △ Less

Submitted 1 December, 2019; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: Accepted for publication at ApJ Letters

Double neutron star mergers from hierarchical triple-star systems

Authors: Adrian S. Hamers, Todd A. Thompson

Abstract: 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-st… ▽ More 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 (σ_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 (σ_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. △ Less

Submitted 12 August, 2019; v1 submitted 18 July, 2019; originally announced July 2019.

Comments: Accepted for publication in ApJ. 16 pages, 8 figures

Analytic computation of the secular effects of encounters on a binary: third-order perturbation, octupole, and post-Newtonian terms; steady-state distribution

Authors: Adrian S. Hamers, Johan Samsing

Abstract: 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" en… ▽ More 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. △ Less

Submitted 22 July, 2019; v1 submitted 20 June, 2019; originally announced June 2019.

Comments: Accepted for publication in MNRAS. 17 pages, 8 figures. Code can be found at https://github.com/hamers/flybys

The Effect of Distant Encounters on Black Hole Binaries in Globular Clusters: Systematic Increase of In-cluster Mergers in the LISA band

Authors: Johan Samsing, Adrian S. Hamers, Jacob G. Tyles

Abstract: 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 measur… ▽ More 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. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Comments: 19 pages. 12 figures. comments very welcome

Journal ref: Phys. Rev. D 100, 043010 (2019)

The impact of white dwarf natal kicks and stellar flybys on the rates of Type Ia supernovae in triple-star systems

Authors: Adrian S. Hamers, Todd A. Thompson

Abstract: 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 o… ▽ More 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. △ Less

Submitted 11 July, 2019; v1 submitted 29 April, 2019; originally announced April 2019.

Comments: Accepted for publication in ApJ. 12 pages, 9 figures

Analytic computation of the secular effects of encounters on a binary: features arising from second-order perturbation theory

Authors: Adrian S. Hamers, Johan Samsing

Abstract: 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… ▽ More 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. △ Less

Submitted 28 June, 2019; v1 submitted 21 April, 2019; originally announced April 2019.

Comments: Accepted for publication in MNRAS. One reference updated. 18 pages, 7 figures. Code can be found at https://github.com/hamers/flybys

An analytic model for mass transfer in binaries with arbitrary eccentricity, with applications to triple-star systems

Authors: Adrian S. Hamers, Fani Dosopoulou

Abstract: 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,… ▽ More 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. △ Less

Submitted 2 February, 2019; v1 submitted 13 December, 2018; originally announced December 2018.

Comments: Minor changes to match ApJ publication. 16 pages, 10 figures. Code can be found at https://github.com/hamers/emt

Catching a planet: A tidal capture origin for the exomoon candidate Kepler 1625b I

Authors: Adrian S. Hamers, Simon F. Portegies Zwart

Abstract: 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… ▽ More 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. △ Less

Submitted 24 November, 2018; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: Accepted for publication in ApJL. 7 pages, 5 figures

Shrinking orbits in hierarchical quadruple star systems

Authors: Adrian S. Hamers

Abstract: 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… ▽ More 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. △ Less

Submitted 6 November, 2018; v1 submitted 22 October, 2018; originally announced October 2018.

Comments: Minor revision to match MNRAS publication. 15 pages, 13 figures

Stability of exomoons around the Kepler transiting circumbinary planets

Authors: Adrian S. Hamers, Maxwell X. Cai, Javier Roa, Nathan Leigh

Abstract: 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$-bo… ▽ More 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. △ Less

Submitted 2 August, 2018; v1 submitted 15 June, 2018; originally announced June 2018.

Comments: Accepted for publication by MNRAS. 12 pages, 7 figures

The impact of vector resonant relaxation on the evolution of binaries near a massive black hole: implications for gravitational wave sources

Authors: Adrian S. Hamers, Ben Bar-Or, Cristobal Petrovich, Fabio Antonini

Abstract: 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 VIR… ▽ More 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}. △ Less

Submitted 14 August, 2018; v1 submitted 25 May, 2018; originally announced May 2018.

Comments: Accepted for publication in ApJ. 22 pages, 20 figures

On the rates of type Ia supernovae originating from white dwarf collisions in quadruple star systems

Authors: Adrian S. Hamers

Abstract: 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… ▽ More 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. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: Accepted for publication in MNRAS. 18 pages, 13 figures

Secular dynamics of hierarchical multiple systems composed of nested binaries, with an arbitrary number of bodies and arbitrary hierarchical structure. II. External perturbations: flybys and supernovae

Authors: Adrian S. Hamers

Abstract: 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 separati… ▽ More 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. △ Less

Submitted 15 February, 2018; originally announced February 2018.

Comments: Accepted for publication in MNRAS. 24 pages, 7 figures

Hot Jupiters driven by high-eccentricity migration in globular clusters

Authors: Adrian S. Hamers, Scott Tremaine

Abstract: 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… ▽ More 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. △ Less

Submitted 29 November, 2017; v1 submitted 29 September, 2017; originally announced October 2017.

Comments: Accepted for publication in AJ. 19 pages, 15 figures

The rate of WD-WD head-on collisions in isolated triples is too low to explain standard type Ia supernovae

Authors: S. Toonen, H. B. Perets, A. S. Hamers

Abstract: 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-S… ▽ More 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. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: 13 pages, 4 figures, submitted to A&A

Journal ref: A&A 610, A22 (2018)

Secular chaotic dynamics in hierarchical quadruple systems, with applications to hot Jupiters in stellar binaries and triples

Authors: Adrian S. Hamers, Dong Lai

Abstract: 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 eccentricit… ▽ More 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. △ Less

Submitted 15 June, 2017; v1 submitted 5 May, 2017; originally announced May 2017.

Comments: Revised to match MNRAS publication. 16 pages, 10 figures

Relaxation near supermassive black holes driven by nuclear spiral arms: anisotropic hypervelocity stars, S-stars and tidal disruption events

Authors: Adrian S. Hamers, Hagai B. Perets

Abstract: 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), causi… ▽ More 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. △ Less

Submitted 10 July, 2017; v1 submitted 20 April, 2017; originally announced April 2017.

Comments: Accepted for publication in ApJ. 21 pages, 18 figures

Binary black hole mergers from field triples: properties, rates and the impact of stellar evolution

Authors: Fabio Antonini, Silvia Toonen, Adrian S. Hamers

Abstract: 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 evoluti… ▽ More 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. △ Less

Submitted 19 April, 2017; v1 submitted 20 March, 2017; originally announced March 2017.

Comments: Accepted for publication in ApJ. 10 pages, 6 figures

Hints for hidden planetary companions to hot Jupiters in stellar binaries

Authors: Adrian S. Hamers

Abstract: 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 bin… ▽ More 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. △ Less

Submitted 16 January, 2017; originally announced January 2017.

Comments: Accepted for publication in ApJL. 6 pages, 4 figures

On the formation of hot and warm Jupiters via secular high-eccentricity migration in stellar triples

Authors: Adrian S. Hamers

Abstract: 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 st… ▽ More 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. △ Less

Submitted 6 January, 2017; originally announced January 2017.

Comments: Accepted for publication in MNRAS. 14 pages, 11 figures

White dwarf pollution by planets in stellar binaries

Authors: Adrian S. Hamers, Simon F. Portegies Zwart

Abstract: 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 compa… ▽ More 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. △ Less

Submitted 5 July, 2016; originally announced July 2016.

Comments: Accepted for publication in MNRAS Letters. 5 pages, 5 figures

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

Authors: Adrian S. Hamers, Fabio Antonini, Yoram Lithwick, Hagai B. Perets, Simon F. Portegies Zwart

Abstract: 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 thre… ▽ More 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. △ Less

Submitted 15 December, 2016; v1 submitted 23 June, 2016; originally announced June 2016.

Comments: Revised to match MNRAS publication. 14 pages, 10 figures

Journal ref: 2017, Monthly Notices of the Royal Astronomical Society, Volume 464, Issue 1, p.688-701

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

Authors: Fabio Antonini, Adrian S. Hamers, Yoram Lithwick

Abstract: 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 place… ▽ More 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. △ Less

Submitted 18 August, 2016; v1 submitted 6 April, 2016; originally announced April 2016.

Comments: Accepted to AAS journals (AJ). 15 pages, 9 figures

Mass transfer between debris discs during close stellar encounters

Authors: Lucie Jilkova, Adrian S. Hamers, Michael Hammer, Simon Portegies Zwart

Abstract: 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… ▽ More 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. △ Less

Submitted 29 January, 2016; originally announced January 2016.

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

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

Authors: Adrian S. Hamers, Simon F. Portegies Zwart

Abstract: 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 (quadrup… ▽ More 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. △ Less

Submitted 9 June, 2016; v1 submitted 3 November, 2015; originally announced November 2015.

Comments: Revised to match MNRAS publication. 48 pages, 36 figures

Journal ref: 2016, Monthly Notices of the Royal Astronomical Society, Volume 459, Issue 3, p.2827-2874

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

Authors: Adrian S. Hamers, Hagai B. Perets, Simon F. Portegies Zwart

Abstract: 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… ▽ More 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. △ Less

Submitted 9 June, 2016; v1 submitted 5 June, 2015; originally announced June 2015.

Comments: Revised to match MNRAS publication. 24 pages, 22 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 455, Issue 3, p.3180-3200 (2016)

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

Authors: Adrian S. Hamers, Hagai B. Perets, Fabio Antonini, Simon F. Portegies Zwart

Abstract: 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 o… ▽ More 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. △ Less

Submitted 9 June, 2016; v1 submitted 9 December, 2014; originally announced December 2014.

Comments: Revised to match MNRAS publication. 27 pages, 16 figures

Journal ref: 2015, Monthly Notices of the Royal Astronomical Society, Volume 449, Issue 4, p.4221-4245

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

Authors: Adrian S. Hamers, Simon F. Portegies Zwart

Abstract: 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 planetes… ▽ More 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. △ Less

Submitted 21 October, 2014; originally announced October 2014.

Comments: Accepted for publication in MNRAS. 21 pages, 12 figures