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

2021 ◽  
Vol 502 (3) ◽  
pp. 4479-4512
Author(s):  
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 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.

2004 ◽  
Vol 194 ◽  
pp. 243-243
Author(s):  
Stephen Justham ◽  
Ulrich Kolb

Studies of binary evolution have, until recently, neglected thermal timescale mass transfer (TTMT). Recent work has suggested that this previously poorly studied area is crucial in the understanding of systems across the compact binary spectrum. We use the state-of-the-art binary population synthesis code BiSEPS (Willems and Kolb, 2002, MNRAS 337 1004-1016). However, the present treatment of TTMT is incomplete due to the nonlinear behaviour of stars in their departure from gravothermal ‘equilibrium’. Here we show work that should update the ultrafast stellar evolution algorithms within BiSEPS to make it the first pseudo-analytic code that can follow TTMT properly. We have generated fits to a set of over 300 Case B TTMT sequences with a range of intermediate-mass donors. These fits produce very good first approximations to both HR diagrams and mass-transfer rates (see figures 1 and 2), which we later hope to improve and extend. They are already a significant improvement over the previous fits.


Author(s):  
J. J. Eldridge ◽  
N. -Y. Guo ◽  
N. Rodrigues ◽  
E. R. Stanway ◽  
L. Xiao

Abstract We use the results of a supernova light-curve population synthesis to predict the range of possible supernova light curves arising from a population of single-star progenitors that lead to type IIP supernovae. We calculate multiple models varying the initial mass, explosion energy, nickel mass and nickel mixing and then compare these to type IIP supernovae with detailed light curve data and pre-explosion imaging progenitor constraints. Where a good fit is obtained to observations, we are able to achieve initial progenitor and nickel mass estimates from the supernova lightcurve that are comparable in precision to those obtained from progenitor imaging. For 2 of the 11 IIP supernovae considered our fits are poor, indicating that more progenitor models should be included in our synthesis or that our assumptions, regarding factors such as stellar mass loss rates or the rapid final stages of stellar evolution, may need to be revisited in certain cases. Using the results of our analysis we are able to show that most of the type IIP supernovae have an explosion energy of the order of log(Eexp/ergs) = 50.52 ± 0.10 and that both the amount of nickel in the supernovae and the amount of mixing may have a dependence on initial progenitor mass.


2003 ◽  
Vol 214 ◽  
pp. 109-112 ◽  
Author(s):  
Zhanwen Han ◽  
Philipp Podsiadlowski

Using binary population synthesis (BPS), we studied the birthrates of SNe Ia for two progenitor models – the single degenerate model and the double degenerate model. We find that the birthrates from both models are within a factor of a few comparable to those inferred observationally. For each model, we investigate different star-formation histories (single star burst or constant star formation), different metallicities (Z = 0.02, 0.004, 0.001), different parameters for the BPS model (mass transfer efficiency during stable mass transfer, common-envelope ejection parameters) and obtain the evolution of birthrates with time.


1998 ◽  
Vol 11 (1) ◽  
pp. 430-432
Author(s):  
Ted Von Hippel

The study of cluster white dwarfs (WDs) has been invigorated recently bythe Hubble Space Telescope (HST). Recent WD studies have been motivated by the new and independent cluster distance (Renzini et al. 1996), age (von Hippel et al. 1995; Richer et al. 1997), and stellar evolution (Koester & Reimers 1996) information that cluster WDs can provide. An important byproduct of these studies has been an estimate of the WD mass contribution in open and globular clusters. The cluster WD mass fraction is of importance for understanding the dynamical state and history of star clusters. It also bears an important connection to the WD mass fractions of the Galactic disk and halo. Current evidence indicates that the open clusters (e.g. von Hippel et al. 1996; Reid this volume) have essentially the same luminosity function (LF) as the solar neighborhood population. The case for the halo is less clear, despite the number of very good globular cluster LFs down to nearly 0.1 solar masses (e.g. Cool et al. 1996; Piotto, this volume), as the field halo LF is poorly known. For most clusters dynamical evolution should cause evaporation of the lowest mass members, biasing clusters to have flatter present-day mass functions (PDMFs) than the disk and halo field populations. Dynamical evolution should also allow cluster WDs to escape, though not in the same numbers as the much lower mass main sequence stars. The detailed connection between cluster PDMFs and the field IMF awaits elucidation from observations and the new combined N-body and stellar evolution models (Tout, this volume). Nevertheless, the WD mass fraction of clusters already provides an estimate for the WD mass fraction of the disk and halo field populations. A literature search to collect cluster WDs and a simple interpretive model follow. This is a work in progress and the full details of the literature search and the model will be published elsewhere.


1985 ◽  
Vol 113 ◽  
pp. 139-160 ◽  
Author(s):  
Douglas C. Heggie

This review describes work on the evolution of a stellar system during the phase which starts at the end of core collapse. It begins with an account of the models of Hénon, Goodman, and Inagaki and Lynden-Bell, as well as evaporative models, and modifications to these models which are needed in the core. Next, these models are related to more detailed numerical calculations of gaseous models, Fokker-Planck models, N-body calculations, etc., and some problems for further work in these directions are outlined. The review concludes with a discussion of the relation between theoretical models and observations of the surface density profiles and statistics of actual globular clusters.


2015 ◽  
Vol 05 (03) ◽  
pp. 222-241 ◽  
Author(s):  
Seblu Humne Negu ◽  
Solomon Belay Tessema

2021 ◽  
Vol 922 (2) ◽  
pp. 174
Author(s):  
Kenny X. Van ◽  
Natalia Ivanova

Abstract We present a new method for constraining the mass transfer evolution of low-mass X-ray binaries (LMXBs)—a reverse population synthesis technique. This is done using the detailed 1D stellar evolution code MESA (Modules for Experiments in Stellar Astrophysics) to evolve a high-resolution grid of binary systems spanning a comprehensive range of initial donor masses and orbital periods. We use the recently developed convection and rotation-boosted (CARB) magnetic braking scheme. The CARB magnetic braking scheme is the only magnetic braking prescription capable of reproducing an entire sample of well-studied persistent LMXBs—those with mass ratios, periods, and mass transfer rates that have been observationally determined. Using the reverse population synthesis technique, where we follow any simulated system that successfully reproduces an observed LMXB backward, we have constrained possible progenitors for each observed well-studied persistent LMXB. We also determined that the minimum number of LMXB formations in the Milky Way is 1500 per Gyr if we exclude Cyg X-2. For Cyg X-2, the most likely formation rate is 9000 LMXB Gyr−1. The technique we describe can be applied to any observed LMXB with well-constrained mass ratio, period, and mass transfer rate. With the upcoming GAIA DR3 containing information on binary systems, this technique can be applied to the data release to search for progenitors of observed persistent LMXBs.


1996 ◽  
Vol 174 ◽  
pp. 375-376
Author(s):  
P.-Y. Longaretti ◽  
C. Lagoute

We have computed simplified globular cluster evolutionary tracks which take into account the effects of internal relaxation, of the cluster rotation, of the galactic tidal field, and, in a cruder way, of stellar evolution and of gravitational shocking. The objectives are first to quantify the influence of rotation in the dynamical evolution of globular clusters; and second, to investigate the evolution of globular cluster angular momentum and flattening (Lagoute and Longaretti 1995a, Longaretti and Lagoute 1995b,c).


1977 ◽  
Vol 42 ◽  
pp. 365-370
Author(s):  
Józef Smak

The mass loss from cataclysmic binaries seems an important and worth studying phenomenon for a number of reasons. It is probably enough to mention only two of them:(a) Whenever we can directly observe the ejected material, determine its amount and the rate of mass loss, as well as its chemical composition (this being the case of the expanding envelopes of novae), we are getting a good insight into the basic physical mechanisms responsible for the observed phenomena.(b) The mass loss (together with the mass transfer) and the loss of the orbital angular momentum are related directly to the dynamical evolution of a binary system and - indirectly - to the evolution of its components.


Sign in / Sign up

Export Citation Format

Share Document