BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Adrian S. Hamers

Adrian S. Hamers contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRastro-ph.EPastro-ph.COastro-ph.IMgr-qcMachine Learningnlin.CDphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

12 published item(s)

preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Algebraic and machine learning approach to hierarchical triple-star stability

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Equilibrium Tide: An Updated Prescription for Population Synthesis Codes

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Binary-binary scattering in the secular limit

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access