Evolutionary Models Of Massive Stars Research Articles

Exploring extreme brightness variations in blue supergiant MACHO 80.7443.1718: Evidence for companion-driven enhanced mass loss

Context.The evolution of massive stars is dominated by interactions within binary and multiple systems. In order to accurately model this evolution, it is necessary to investigate all possible forms of an interaction in binary systems that may affect the evolution of the components. One of the “laboratories” plausible for this kind of investigation is the massive eccentric binary system MACHO 80.7443.1718 (ExtEV), which exhibits an exceptionally large amplitude of light variability close to the periastron passage of its 32.8-day orbit.Aims.We examine whether the light variability of ExtEV can be explained by a wind-wind collision (WWC) binary system model. We also critically review other models proposed to explain the light curve of ExtEV.Methods.We conducted an analysis of (i) the broadband multicolor photometry of ExtEV spanning a wide range of wavelengths from the ultraviolet to near-infrared, (ii) the time-series space photometry from the Transiting Exoplanet Survey Satellite (TESS), (iii) ground-based JohnsonUBVphotometry, and (iv) time-series high-resolution spectroscopy. To derive the parameters of the primary component of the system, we fit the spectral energy distribution (SED) and calculated evolutionary models of massive stars that included mass loss. Using radial-velocity data, we determined the spectroscopic parameters of the system. We also fit an analytical model of light variations to the TESS light curve of ExtEV.Results.The ExtEV system exhibits an infrared excess, indicating an increased mass-loss rate. The system does not match the characteristics of B[e] stars, however. We rule out the possibility of the presence of a Keplerian disk around the primary component. We also argue that the scenario with periodic Roche-lobe overflow at periastron may not be consistent with the observations of ExtEV. Analysis of the SED suggests that the primary component has a radius of about 30R⊙and a luminosity of ∼6.6 × 105L⊙. With the analysis of the radial-velocity data, we refine the orbital parameters of ExtEV and find evidence for the presence of a tertiary component in the system. Using evolutionary models we demonstrate that the primary component’s mass is between 25 and 45M⊙. We successfully reproduced the light curve of ExtEV with our analytical model, showing that the dominant processes shaping its light curve can be attributed to the atmospheric eclipse and light scattered in the WWC cone. We also estimate the primary’s mass loss rate due to stellar wind for 4.5 × 10−5M⊙yr−1.Conclusions.ExtEV is most likely not an extreme eccentric ellipsoidal variable, but rather an exceptional WWC binary system. The mass loss rate we derived exceeds theoretical predictions by up to two orders of magnitude. This implies that the wind in the system is likely enhanced by tidal interactions, rotation, and possibly also tidally excited oscillations. Therefore, ExtEV represents a rare evolutionary phase of a binary system that may help to understand the role of companion-driven enhanced mass loss in the evolution of massive binary systems.

A Multiwavelength Classification and Study of Red Supergiant Candidates in NGC 6946

We have combined resolved stellar photometry from Hubble Space Telescope (HST), Spitzer, and Gaia to identify red supergiant (RSG) candidates in NGC 6946, based on their colors, proper motions, visual morphologies, and spectral energy distributions. We start with a large sample of 17,865 RSG candidates based solely on HST near-infrared photometry. We then chose a small sample of 385 of these candidates with Spitzer matches for a more detailed study. Using evolutionary models and isochrones, we isolate a space where RSGs would be found in our photometry catalogs. We then visually inspect each candidate and compare them to Gaia catalogs to identify and remove foreground stars. As a result, we classify 95 potential RSGs, with 40 of these being in our highest-quality sample. We fit the photometry of the populations of stars in the regions surrounding the RSGs to infer their ages. Placing our best candidate RSG stars into three age bins between 1 and 30 Myr, we find 27.5% of the candidates falling between 1–10 Myr, 37.5% between 10–20 Myr, and 35% between 20–30 Myr. A comparison of our results to the models of massive star evolution shows some agreement between model luminosities and the luminosities of our candidates for each age. Three of our candidates appear significantly more consistent with binary models than single-star evolution models.

The effects of surface fossil magnetic fields on massive star evolution: IV. Grids of models at Solar, LMC, and SMC metallicities

ABSTRACT Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the mesa software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3–60 M⊙), surface equatorial magnetic field strength (0–50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly rotating, nitrogen-enriched (‘Group 2’) stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of (i) revisiting the derived stellar parameters of known magnetic stars, and (ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes.

3D hydrodynamics simulations of massive stars with the PROMPI code

Evolutionary models of massive stars are affected by uncertainties related to three-dimensional processes like convection. However, it is not computationally possible to simulate at the same time both the entire lifetime of a star and its fluid motions. 3D hydrodynamics simulations can be employed to study the fluid behaviour in very high detail for a short time-scale, which can be used to improve 1D models. The code PROMPI is a multi-dimensional hydrodynamics software developed specifically for simulations of the late phases in massive star interiors. In these proceedings, we present results from a new PROMPI simulation of the neon-burning shell in a 15 M© star, focusing on the interplay between convective turbulence and nuclear reactions. In particular, we show how it is possible to use the chemical composition to study the time evolution of the convective zone, which is subjected to turbulent entrainment.

Weakening the wind with ULLYSES: Examining the Bi-Stability Jump

AbstractRadiation-driven mass-loss is an important, but still highly debated, driver for the evolution of massive stars. Current massive star evolution models rely on the theoretical prediction that low luminosity massive stars experience a sudden increase in mass loss below a stellar effective temperature of about 20 000 K. However, novel radiation-driven mass-loss rate predictions show no such bi-stability jump, which effects the post main-sequence evolution of massive stars. The ULLYSES data set provides a unique opportunity to investigate the theoretical bi-stability jump dichotomy and may help to assess the existence of the bi-stability jump in massive star winds. By utilising UV spectra from ULLYSES combined with X-shooter optical data we obtain empirical mass-loss rate constraints, that are no longer degenerate to the effects of wind clumping, and derive novel empirical constraints on the mass-loss behavior across the temperature range of the bi-stability jump. Current preliminary results do not show a clear presence of a bi-stability jump.

Pulsational properties of ten new slowly pulsating B stars

Context. Slowly pulsating B (SPB) stars are upper main-sequence multi-periodic pulsators that show non-radial g-mode oscillations driven by the κ mechanism acting on the iron bump. These multi-periodic pulsators have great asteroseismic potential and can be employed for the calibration of stellar structure and evolution models of massive stars. Aims. We collected a sample of ten hitherto unidentified SPB stars with the aim of describing their pulsational properties and identifying pulsational modes. Methods. Photometric time series data from various surveys were collected and analyzed using diverse frequency search algorithms. We calculated astrophysical parameters and investigated the location of our sample stars in the log Teff vs. log L/L⊙ diagram. Current pulsational models were calculated and used for the identification of pulsational modes in our sample stars. An extensive grid of stellar models along with their g-mode eigenfrequencies was calculated and subsequently cross-matched with the observed pulsational frequencies. The best-fit models were then used in an attempt to constrain stellar parameters such as mass, age, metallicity, and convective overshoot. Results. We present detected frequencies, corresponding g-mode identifications, and the masses and ages of the stellar models producing the best frequency cross-matches. We partially succeeded in constraining stellar parameters, in particular concerning mass and age. Where applicable, rotation periods have been derived from the spacing of triplet component frequencies. No evolved SPB stars are present in our sample. We identify two candidate high-metallicity objects (HD 86424 and HD 163285), one young SPB star (HD 36999), and two candidate young SPB stars (HD 61712 and HD 61076). Conclusions. We demonstrate the feasibility of using ground-based observations to perform basic asteroseismological analyses of SPB stars. Our results significantly enlarge the sample of known SPB stars with reliable pulsational mode identifications, which provides important input parameters for modeling attempts aiming to investigate the internal processes at work in upper main-sequence stars.

New predictions for radiation-driven, steady-state mass-loss and wind-momentum from hot, massive stars

Context. Radiation-driven mass loss plays a key role in the life cycles of massive stars. However, basic predictions of such mass loss still suffer from significant quantitative uncertainties. Aims. We develop new radiation-driven, steady-state wind models for massive stars with hot surfaces, suitable for quantitative predictions of global parameters like mass-loss and wind-momentum rates. Methods. The simulations presented here are based on a self-consistent, iterative grid solution to the spherically symmetric, steady-state equation of motion, using full non-local thermodynamic equilibrium radiative transfer solutions in the co-moving frame to derive the radiative acceleration. We do not rely on any distribution functions or parametrization for computation of the line force responsible for the wind driving. The models start deep in the subsonic and optically thick atmosphere and extend up to a large radius at which the terminal wind speed has been reached. Results. In this first paper, we present models representing two prototypical O-stars in the Galaxy, one with a higher stellar mass M*∕M⊙ = 59 and luminosity log10L*∕L⊙ = 5.87 (spectroscopically an early O supergiant) and one with a lower M*∕M⊙ = 27 and log10L*∕L⊙ = 5.1 (a late O dwarf). For these simulations, basic predictions for global mass-loss rates, velocity laws, and wind momentum are given, and the influence from additional parameters like wind clumping and microturbulent speeds is discussed. A key result is that although our mass-loss rates agree rather well with alternative models using co-moving frame radiative transfer, they are significantly lower than those predicted by the mass-loss recipes normally included in models of massive-star evolution. Conclusions. Our results support previous suggestions that Galactic O-star mass-loss rates may be overestimated in present-day stellar evolution models, and that new rates might therefore be needed. Indeed, future papers in this series will incorporate our new models into such simulations of stellar evolution, extending the very first simulations presented here toward larger grids covering a range of metallicities, B supergiants across the bistability jump, and possibly also Wolf-Rayet stars.

No evidence for intermediate-mass black holes in the globular clusters ω Cen and NGC 6624

ABSTRACT We compare the results of a large grid of N-body simulations with the surface brightness and velocity dispersion profiles of the globular clusters ω Cen and NGC 6624. Our models include clusters with varying stellar-mass black hole retention fractions and varying masses of a central intermediate-mass black hole (IMBH). We find that an $\sim 45\, 000$ M⊙ IMBH, whose presence has been suggested based on the measured velocity dispersion profile of ω Cen, predicts the existence of about 20 fast-moving, m > 0.5 M⊙, main-sequence stars with a (1D) velocity v > 60 km s−1 in the central 20 arcsec of ω Cen. However, no such star is present in the HST/ACS proper motion catalogue of Bellini et al. (2017), strongly ruling out the presence of a massive IMBH in the core of ω Cen. Instead, we find that all available data can be fitted by a model that contains 4.6 per cent of the mass of ω Cen in a centrally concentrated cluster of stellar-mass black holes. We show that this mass fraction in stellar-mass BHs is compatible with the predictions of stellar evolution models of massive stars. We also compare our grid of N-body simulations with NGC 6624, a cluster recently claimed to harbour a 20 000 M⊙ black hole based on timing observations of millisecond pulsars. However, we find that models with MIMBH > 1000 M⊙ IMBHs are incompatible with the observed velocity dispersion and surface brightness profile of NGC 6624, ruling out the presence of a massive IMBH in this cluster. Models without an IMBH provide again an excellent fit to NGC 6624.

Constraining mixing in massive stars in the Small Magellanic Cloud

Context. The evolution of massive stars is strongly influenced by internal mixing processes such as semiconvection, convective core overshooting, and rotationally induced mixing. None of these processes are currently well constrained. Aims. We investigate models for massive stars in the Small Magellanic Cloud (SMC), for which stellar-wind mass loss is less important than for their metal-rich counterparts. We aim to constrain the various mixing efficiencies by comparing model results to observations. Methods. For this purpose, we use the stellar-evolution code MESA to compute more than 60 grids of detailed evolutionary models for stars with initial masses of 9…100 M⊙, assuming different combinations of mixing efficiencies of the various processes in each grid. Our models evolve through core hydrogen and helium burning, such that they can be compared with the massive main sequence and supergiant population of the SMC. Results. We find that for most of the combinations of the mixing efficiencies, models in a wide mass range spend core-helium burning either only as blue supergiants, or only as red supergiants. The latter case corresponds to models that maintain a shallow slope of the hydrogen/helium (H/He) gradient separating the core and the envelope of the models. Only a small part of the mixing parameter space leads to models that produce a significant number of blue and red supergiants, which are both in abundance in the SMC. Some of our grids also predict a cut-off in the number of red supergiants above log L/L⊙ = 5…5.5. Interestingly, these models contain steep H/He gradients, as is required to understand the hot, hydrogen-rich Wolf-Rayet stars in the SMC. We find that unless it is very fast, rotation has a limited effect on the H/He profiles in our models. Conclusions. While we use specific implementations of the considered mixing processes, they comprehensively probe the two first-order structural parameters, the core mass and the H/He gradient in the core-envelope interface. Our results imply that in massive stars, mixing during the main-sequence evolution leads to a moderate increase in the helium core masses, and also that the H/He gradients above the helium cores become very steep. Our model grids can be used to further refine the various mixing efficiencies with the help of future observational surveys of the massive stars in the SMC, and thereby help to considerably reduce the uncertainties in models of massive star evolution.

Massive star evolution: rotation, winds, and overshooting vectors in the mass-luminosity plane

Context. Massive star evolution is dominated by various physical effects, including mass loss, overshooting, and rotation, but the prescriptions of their effects are poorly constrained and even affect our understanding of the main sequence. Aims. We aim to constrain massive star evolution models using the unique test-bed eclipsing binary HD 166734 with new grids of MESA stellar evolution models, adopting calibrated prescriptions of overshooting, mass loss, and rotation. Methods. We introduce a novel tool, called the mass-luminosity plane or M−L plane, as an equivalent to the traditional HR diagram, utilising it to reproduce the test-bed binary HD 166734 with newly calibrated MESA stellar evolution models for single stars. Results. We can only reproduce the Galactic binary system with an enhanced amount of core overshooting (αov = 0.5), mass loss, and rotational mixing. We can utilise the gradient in the M−L plane to constrain the amount of mass loss to 0.5–1.5 times the standard prescription test-bed, and we can exclude extreme reduction or multiplication factors. The extent of the vectors in the M−L plane leads us to conclude that the amount of core overshooting is larger than is normally adopted in contemporary massive star evolution models. We furthermore conclude that rotational mixing is mandatory to obtain the correct nitrogen abundance ratios between the primary and secondary components (3:1) in our test-bed binary system. Conclusions. Our calibrated grid of models, alongside our new M−L plane approach, present the possibility of a widened main sequence due to an increased demand for core overshooting. The increased amount of core overshooting is not only needed to explain the extended main sequence, but the enhanced overshooting is also needed to explain the location of the upper-luminosity limit of the red supergiants. Finally, the increased amount of core overshooting has – via the compactness parameter – implications for supernova explodability.

Spectroscopic and physical parameters of Galactic O-type stars

Context.Massive stars play a key role in the evolution of galaxies and our Universe.Aims.Our goal is to compare observed and predicted properties of single Galactic O stars to identify and constrain uncertain physical parameters and processes in stellar evolution and atmosphere models.Methods.We used a sample of 53 objects of all luminosity classes and with spectral types from O3 to O9.7. For 30 of these, we determined the main photospheric and wind parameters, including projected rotational rates accounting for macroturbulence, and He and N surface abundances, using optical spectroscopy and applying the model atmosphere code FASTWIND. For the remaining objects, similar data from the literature, based on analyses by means of the CMFGEN code, were used instead. The properties of our sample were then compared to published predictions based on two grids of single massive star evolution models that include rotationally induced mixing.Results.Any of the considered model grids face problem in simultaneously reproducing the stellar masses, equatorial gravities, surface abundances, and rotation rates of our sample stars. The spectroscopic masses derived for objects below ~30M⊙tend to be smaller than the evolutionary ones, no matter which of the two grids have been used as a reference. While this result may indicate the need to improve the model atmosphere calculations (e.g. regarding the treatment of turbulent pressure), our analysis shows that the established mass problem cannot be fully explained in terms of inaccurate parameters obtained by quantitative spectroscopy or inadequate model values ofVroton the zero age main sequence. Within each luminosity class, we find a close correlation of N surface abundance and luminosity, and a stronger N enrichment in more massive and evolved O stars. Additionally, we also find a correlation of the surface nitrogen and helium abundances. The large number of nitrogen-enriched stars above ~30M⊙argues for rotationally induced mixing as the most likely explanation. However, none of the considered models can match the observed trends correctly, especially in the high mass regime.Conclusions.We confirm mass discrepancy for objects in the low mass O-star regime. We conclude that the rotationally induced mixing of helium to the stellar surface is too strong in some of the models. We also suggest that present inadequacies of the models to represent the N enrichment in more massive stars with relatively slow rotation might be related (among other issues) to problematic efficiencies of rotational mixing. We are left with a picture in which invoking binarity and magnetic fields is required to achieve a more complete agreement of the observed surface properties of a population of massive main-sequence stars with corresponding evolutionary models.

Confronting Models of Massive Star Evolution and Explosions with Remnant Mass Measurements

Abstract The mass distribution of compact objects provides a fossil record that can be studied to uncover information on the late stages of massive star evolution, the supernova explosion mechanism, and the dense matter equation of state. Observations of neutron star masses indicate a bimodal Gaussian distribution, while the observed black hole mass distribution decays exponentially for stellar-mass black holes. We use these observed distributions to directly confront the predictions of stellar evolution models and the neutrino-driven supernova simulations of Sukhbold et al. We find strong agreement between the black hole and low-mass neutron star distributions created by these simulations and the observations. We show that a large fraction of the stellar envelope must be ejected, either during the formation of stellar-mass black holes or prior to the implosion through tidal stripping due to a binary companion, in order to reproduce the observed black hole mass distribution. We also determine the origins of the bimodal peaks of the neutron star mass distribution, finding that the low-mass peak (centered at ∼1.4 M ⊙) originates from progenitors with M ZAMS ≈ 9–18 M ⊙. The simulations fail to reproduce the observed peak of high-mass neutron stars (centered at ∼1.8 M ⊙) and we explore several possible explanations. We argue that the close agreement between the observed and predicted black hole and low-mass neutron star mass distributions provides new, promising evidence that these stellar evolution and explosion models capture the majority of relevant stellar, nuclear, and explosion physics involved in the formation of compact objects.

Wolf–Rayet stars in the Small Magellanic Cloud as testbed for massive star evolution

Context. The majority of the Wolf–Rayet (WR) stars represent the stripped cores of evolved massive stars who lost most of their hydrogen envelope. Wind stripping in single stars is expected to be inefficient in producing WR stars in metal-poor environments such as the Small Magellanic Cloud (SMC). While binary interaction can also produce WR stars at low metallicity, it is puzzling that the fraction of WR binaries appears to be about 40%, independent of the metallicity.Aim. We aim to use the recently determined physical properties of the twelve known SMC WR stars to explore their possible formation channels through comparisons with stellar models.Methods. We used the MESA stellar evolution code to construct two grids of stellar models with SMC metallicity. One of these consists of models of rapidly rotating single stars, which evolve in part or completely chemically homogeneously. In a second grid, we analyzed core helium burning stellar models assuming constant hydrogen and helium gradients in their envelopes.Results. We find that chemically homogeneous evolution is not able to account for the majority of the WR stars in the SMC. However, in particular the apparently single WR star SMC AB12, and the double WR system SMC AB5 (HD 5980) appear consistent with this channel. We further find a dichotomy in the envelope hydrogen gradients required to explain the observed temperatures of the SMC WR stars. Shallow gradients are found for the WR stars with O star companions, while much steeper hydrogen gradients are required to understand the group of hot apparently single WR stars.Conclusions. The derived shallow hydrogen gradients in the WR component of the WR+O star binaries are consistent with predictions from binary models where mass transfer occurs early, in agreement with their binary properties. Since the hydrogen profiles in evolutionary models of massive stars become steeper with time after the main sequence, we conclude that most of the hot (Teff > 60 kK ) apparently single WR stars lost their envelope after a phase of strong expansion, e.g., as the result of common envelope evolution with a lower mass companion. The so far undetected companions, either main sequence stars or compact objects, are then expected to still be present. A corresponding search might identify the first immediate double black hole binary progenitor with masses as high as those detected in GW150914.

Modeling the early evolution of massive OB stars with an experimental wind routine

Stellar evolution models of massive stars are very sensitive to the adopted mass-loss scheme. The magnitude and evolution of mass-loss rates significantly affect the main sequence evolution, and the properties of post-main sequence objects, including their rotational velocities. Driven by potential discrepancies between theoretically predicted and observationally derived mass-loss rates in the OB star range, we particularly aim to investigate the response to mass-loss rates that are lower than currently adopted, in parallel with the mass-loss behavior at the "first" bi-stability jump. We perform 1D hydrodynamical model calculations of single $20 - 60 \, M_{\odot}$ Galactic ($Z = 0.014$) stars where the effects of stellar winds are already significant during the main sequence phase. We develop an experimental wind routine to examine the behavior and response of the models under the influence of different mass-loss rates. This observationally guided, simple and flexible wind routine is not a new mass-loss description but a useful tool based on the Wind-momentum Luminosity Relation and other scaling relations, and provides a meaningful base for various tests and comparisons. The main result of this study indicates a dichotomy when accounting for currently debated problems regarding mass-loss rates of hot massive stars. In a fully diffusive approach, and for commonly adopted initial rotational velocities, lower mass-loss rates than theoretically predicted require to invoke an additional source of angular momentum loss (either due to bi-stability braking, or yet unidentified) to brake down surface rotational velocities. [...]

The evolution of magnetic hot massive stars: Implementation of the quantitative influence of surface magnetic fields in modern models of stellar evolution

AbstractLarge-scale dipolar surface magnetic fields have been detected in a fraction of OB stars, however only few stellar evolution models of massive stars have considered the impact of these fossil fields. We are performing 1D hydrodynamical model calculations taking into account evolutionary consequences of the magnetospheric-wind interactions in a simplified parametric way. Two effects are considered: i) the global mass-loss rates are reduced due to mass-loss quenching, and ii) the surface angular momentum loss is enhanced due to magnetic braking. As a result of the magnetic mass-loss quenching, the mass of magnetic massive stars remains close to their initial masses. Thus magnetic massive stars - even at Galactic metallicity - have the potential to be progenitors of ‘heavy’ stellar mass black holes. Similarly, at Galactic metallicity, the formation of pair instability supernovae is plausible with a magnetic progenitor.

The chemical composition of Galactic ring nebulae around massive stars

We present deep spectra of ring nebulae associated with Wolf-Rayet (WR) and O-type stars: NGC 6888, G2.4+1.4, RCW 58, S 308, NGC 7635 and RCW 52. The data have been taken with the 10m Gran Telescopio Canarias and the 6.5m Clay Telescope. We extract spectra of several apertures in some of the objects. We derive C$^{++}$ and O$^{++}$ abundances from faint recombination lines in NGC 6888 and NGC 7635, permitting to derive their C/H and C/O ratios and estimate the abundance discrepancy factor (ADF) of O$^{++}$. The ADFs are larger than the typical ones of normal HII regions but similar to those found in the ionised gas of star-forming dwarf galaxies. We find that chemical abundances are rather homogeneous in the nebulae where we have spectra of several apertures: NGC 6888, NGC 7635 and G2.4+1.4. We obtain very high values of electron temperature in a peripheral zone of NGC 6888, finding that shock excitation can reproduce its spectral properties. We find that all the objects associated with WR stars show N enrichment. Some of them also show He enrichment and O deficiency as well as a lower Ne/O than expected, this may indicate the strong action of the ON and NeNa cycles. We have compared the chemical composition of NGC 6888, G2.4+1.4, RCW 58 and S 308 with the nucleosynthesis predicted by stellar evolution models of massive stars. We find that non-rotational models of stars of initial masses between 25 and 40 solar masses seem to reproduce the observed abundance ratios of most of the nebulae.

Linking 1D evolutionary to 3D hydrodynamical simulations of massive stars

Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective boundaries in stellar models has been studied for decades but still lacks a definitive solution. This paper presents initial results of a study on convective boundary mixing (CBM) in massive stars. The ‘stiffness’ of a convective boundary can be quantified using the bulk Richardson number (), the ratio of the potential energy for restoration of the boundary to the kinetic energy of turbulent eddies. A ‘stiff’ boundary () will suppress CBM, whereas in the opposite case a ‘soft’ boundary () will be more susceptible to CBM. One of the key results obtained so far is that lower convective boundaries (closer to the centre) of nuclear burning shells are ‘stiffer’ than the corresponding upper boundaries, implying limited CBM at lower shell boundaries. This is in agreement with 3D hydrodynamic simulations carried out by Meakin and Arnett (2007 Astrophys. J. 667 448–75). This result also has implications for new CBM prescriptions in massive stars as well as for nuclear burning flame front propagation in super-asymptotic giant branch stars and also the onset of novae.

Period change of massive binaries from combined photometric and spectroscopic data in Cygnus OB2

Context. Mass loss is an important property in evolution models of massive stars. As up to 90% of the massive stars have a visual or spectroscopic companion and many of them exhibit mass exchange, mass-loss rates can be acquired through the period study of massive binaries. Aims. Using our own photometric observations as well as archival data, we look for variations in orbital periods of seven massive eclipsing binary systems in the Cygnus OB2 association and estimate their mass-loss rates and stellar parameters. Methods. We use a Bayesian parameter estimation method to simultaneously fit the period and period change to all available data and a stellar modelling tool to model the binary parameters from photometric and radial-velocity data. Results. Four out of the seven selected binaries show non-zero period change values at two-sigma confidence level. We also report for the first time the eclipsing nature of a star MT059.

THE TRACE OF THE CNO CYCLE IN THE RING NEBULA NGC 6888

We present new results on the chemical composition of the Galactic ring nebula NGC6888 surrounding the WN6(h) star WR136. The data are based on deep spectroscopical observations taken with the High Dispersion Spectrograph at the 8.2m Subaru Telescope. The spectra cover the optical range from 3700 to 7400 A. The effect of the CNO cycle is well identified in the abundances of He, N, and O, while elements not involved in the synthesis such as Ar, S, and Fe present values consistent with the solar vicinity and the ambient gas. The major achievement of this work is the first detection of the faint CII 4267 recombination line in a Wolf-Rayet nebula. This allows to estimate the C abundance in NGC6888 and therefore investigate for the first time the trace of the CNO cycle in a ring nebula around a Wolf-Rayet star. Although the detection of the CII line has a low signal-to-noise ratio, the C abundance seems to be higher than the predictions of recent stellar evolution models of massive stars. The Ne abundance also show a puzzling pattern with an abundance of about 0.5 dex lower than the solar vicinity, which may be related to the action of the NeNa cycle. Attending to the constraints imposed by the dynamical timescale and the He/H and N/O ratios of the nebula, the comparison with stellar evolution models indicates that the initial mass of the stellar progenitor of NGC6888 is between 25 Msun and 40 Msun.

PULSATING B-TYPE STARS IN THE OPEN CLUSTER NGC 884: FREQUENCIES, MODE IDENTIFICATION, AND ASTEROSEISMOLOGY

Recent progress in the seismic interpretation of field beta Cep stars has resulted in improvements of the physics in the stellar structure and evolution models of massive stars. Further asteroseismic constraints can be obtained from studying ensembles of stars in a young open cluster, which all have similar age, distance and chemical composition. We present an observational asteroseismology study based on the discovery of numerous multi-periodic and mono-periodic B-stars in the open cluster NGC 884. We describe a thorough investigation of the pulsational properties of all B-type stars in the cluster. Overall, our detailed frequency analysis resulted in 115 detected frequencies in 65 stars. We found 36 mono-periodic, 16 bi-periodic, 10 tri-periodic, and 2 quadru-periodic stars and one star with 9 independent frequencies. We also derived the amplitudes and phases of all detected frequencies in the U, B, V and I filter, if available. We achieved unambiguous identifications of the mode degree for twelve of the detected frequencies in nine of the pulsators. Imposing the identified degrees and measured frequencies of the radial, dipole and quadrupole modes of five pulsators led to a seismic cluster age estimate of log(age/yr) =7.12-7.28 from a comparison with stellar models. Our study is a proof-of-concept for and illustrates the current status of ensemble asteroseismology of a young open cluster.