450 research outputs found

    Numerical Tests of Rotational Mixing in Massive Stars with the new Population Synthesis Code BONNFIRES

    Full text link
    We use our new population synthesis code BONNFIRES to test how surface abundances predicted by rotating stellar models depend on the numerical treatment of rotational mixing, such as spatial resolution, temporal resolution and computation of mean molecular weight gradients. We find that even with identical numerical prescriptions for calculating the rotational mixing coefficients in the diffusion equation, different timesteps lead to a deviation of the coefficients and hence surface abundances. We find the surface abundances vary by 10-100% between the model sequences with short timestep of 0.001Myr to model sequences with longer timesteps. Model sequences with stronger surface nitrogen enrichment also have longer main-sequence lifetimes because more hydrogen is mixed to the burning cores. The deviations in main-sequence lifetimes can be as large as 20%. Mathematically speaking, no numerical scheme can give a perfect solution unless infinitesimally small timesteps are used. However, we find that the surface abundances eventually converge within 10% between modelling sequences with sufficiently small timesteps below 0.1Myr. The efficiency of rotational mixing depends on the implemented numerical scheme and critically on the computation of the mean molecular weight gradient. A smoothing function for the mean molecular weight gradient results in stronger rotational mixing. If the discretization scheme or the computational recipe for calculating the mean molecular weight gradient is altered, re-calibration of mixing parameters may be required to fit observations. If we are to properly understand the fundamental physics of rotation in stars, it is crucial that we minimize the uncertainty introduced into stellar evolution models when numerically approximating rotational mixing processes.Comment: 8 pages, 6 figures, accepted by A&

    The occurrence of classical Cepheids in binary systems

    Full text link
    Classical Cepheids, like binary stars, are laboratories for stellar evolution and Cepheids in binary systems are especially powerful ones. About one-third of Galactic Cepheids are known to have companions and Cepheids in eclipsing binary systems have recently been discovered in the Large Magellanic Cloud. However, there are no known Galactic binary Cepheids with orbital periods less than one year. We compute population synthesis models of binary Cepheids to compare to the observed period and eccentricity distributions of Galactic Cepheids as well as to the number of observed eclipsing binary Cepheids in the LMC. We find that our population synthesis models are consistent with observed binary properties of Cepheids. Furthermore, we show that binary interaction on the red giant branch prevents some red giant stars from becoming classical Cepheids. Such interactions suggest that the binary fraction of Cepheids should be significantly less than that of their main-sequence progenitors, and that almost all binary Cepheids have orbital periods longer than one year. If the Galactic Cepheid spectroscopic binary fraction is about 35%, then the spectroscopic binary fraction of their intermediate mass main sequence progenitors is about 40-45%.Comment: 7 pages, 3 figures, resubmitted to A&

    BONNSAI: a Bayesian tool for comparing stars with stellar evolution models

    Get PDF
    Powerful telescopes equipped with multi-fibre or integral field spectrographs combined with detailed models of stellar atmospheres and automated fitting techniques allow for the analysis of large number of stars. These datasets contain a wealth of information that require new analysis techniques to bridge the gap between observations and stellar evolution models. To that end, we develop BONNSAI (BONN Stellar Astrophysics Interface), a Bayesian statistical method, that is capable of comparing all available observables simultaneously to stellar models while taking observed uncertainties and prior knowledge such as initial mass functions and distributions of stellar rotational velocities into account. BONNSAI can be used to (1) determine probability distributions of fundamental stellar parameters such as initial masses and stellar ages from complex datasets, (2) predict stellar parameters that were not yet observationally determined and (3) test stellar models to further advance our understanding of stellar evolution. An important aspect of BONNSAI is that it singles out stars that cannot be reproduced by stellar models through Ο‡2\chi^{2} hypothesis tests and posterior predictive checks. BONNSAI can be used with any set of stellar models and currently supports massive main-sequence single star models of Milky Way and Large and Small Magellanic Cloud composition. We apply our new method to mock stars to demonstrate its functionality and capabilities. In a first application, we use BONNSAI to test the stellar models of Brott et al. (2011a) by comparing the stellar ages inferred for the primary and secondary stars of eclipsing Milky Way binaries. Ages are determined from dynamical masses and radii that are known to better than 3%. We find that the stellar models reproduce the Milky Way binaries well. BONNSAI is available through a web-interface at http://www.astro.uni-bonn.de/stars/bonnsai.Comment: Accepted for publication in A&A; 15 pages, 10 figures, 4 tables; BONNSAI is available through a web-interface at http://www.astro.uni-bonn.de/stars/bonnsa

    The HD5980 multiple system: Masses and evolutionary status

    Get PDF
    New spectroscopic observations of the LBV/WR multiple system HD5980 in the Small Magellanic Cloud are used to address the question of the masses and evolutionary status of the two very luminous stars in the 19.3d eclipsing binary system. Two distinct components of the N V 4944 A line are detected in emission and their radial velocity variations are used to derive masses of 61 and 66 Mo, under the assumption that binary interaction effects on this atomic transition are negligible. We propose that this binary system is the product of quasi-chemically homogeneous evolution with little or no mass transfer. Thus, both of these binary stars may be candidates for gamma-ray burst progenitors or even pair instability supernovae. Analysis of the photospheric absorption lines belonging to the third-light object in the system confirm that it consists of an O-type star in a 96.56d eccentric orbit (e=0.82) around an unseen companion. The 5:1 period ratio and high eccentricities of the two binaries suggest that they may constitute a hierarchical quadruple system.Comment: 27 pages, 8 tables, 15 figures; accepted A

    The Origin of B-Type Runaway Stars: Non-LTE Abundances as a Diagnostic

    Get PDF
    There are two accepted mechanisms to explain the origin of runaway OB-type stars: the Binary Supernova Scenario (BSS), and the Cluster Ejection Scenario (CES). In the former, a supernova explosion within a close binary ejects the secondary star, while in the latter close multi-body interactions in a dense cluster cause one or more of the stars to be ejected from the region at high velocity. Both mechanisms have the potential to affect the surface composition of the runaway star. TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and carbon, nitrogen, magnesium and silicon abundances for a sample of B-type runaways. These same analytical tools were used by Hunter et al. (2009) for their analysis of 50 B-type open cluster Galactic stars (i.e. non-runaways). Effective temperatures were deduced using the silicon-ionization balance technique, surface gravities from Balmer line profiles and microturbulent velocities derived using the Si spectrum. The runaways show no obvious abundance anomalies when compared with stars in the open clusters. The runaways do show a spread in composition which almost certainly reflects the Galactic abundance gradient and a range in the birthplaces of the runaways in the Galactic disk. Since the observed Galactic abundance gradients of C, N, Mg and Si are of a similar magnitude, the abundance ratios (e.g., N/Mg) are, as obtained, essentially uniform across the sample

    The Art of Modeling Stellar Mergers and the Case of the B[e] Supergiant R4 in the Small Magellanic Cloud

    Get PDF
    Most massive stars exchange mass with a companion, leading to evolution which is altered drastically from that expected of stars in isolation. Such systems are the result of unusual binary evolution pathways and, as such, may be used to place stringent constraints on the physics of these interactions. We use the R4 system's B[e] supergiant, which has been postulated to be the product of a binary stellar merger, to guide our understanding of such outcomes by comparing observations of R4 to the results of simulations of mergers performed with the 3d hydrodynamics code FLASH. Our approach tailors the simulation initial conditions to the observed properties of R4 and implements realistic stellar profiles generated by the 1d stellar evolution code MESA onto the 3d grid, resolving the merger inspiral to within 0.02 RβŠ™0.02\, R_{\odot}. We then map the merger remnant into MESA to track its evolution on the HR diagram over a period of 10410^4 years. This generates models for a B[e] supergiant with stellar properties, age, and nebula structure in qualitative agreement with that of the R4 system. Our calculations provide concrete evidence to support the idea that R4 was originally a member of a triple system in which the inner binary merged after its most massive member evolved off the main sequence, producing a new object that is of similar mass yet significantly more luminous than the A supergiant companion. The potential applications of the code framework presented in this paper are wide ranging and can be used to generate models of a variety of merger stellar remnants
    • …
    corecore