450 research outputs found
Numerical Tests of Rotational Mixing in Massive Stars with the new Population Synthesis Code BONNFIRES
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
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
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
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
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
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
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 . We then map the
merger remnant into MESA to track its evolution on the HR diagram over a period
of 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
- β¦