BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Norbert Langer

Norbert Langer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEastro-ph.GAastro-phastro-ph.IMastro-ph.COastro-ph.EPcond-mat.quant-gasgr-qcphysics.flu-dyn

Catalog footprint

What is connected

74works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

74 published item(s)

preprint2025arXiv

A new sample of massive B-type contact binary candidates from the OGLE survey of the Magellanic Clouds

Massive contact binaries (CBs) are key to understanding close-binary evolution and stellar mergers, yet their study has been limited by the scarcity of observed systems, particularly of B-type binaries expected to dominate this class. We bridge this gap by mining a large sample of massive CB candidates from the OGLE-IV database, increasing their known numbers in the Magellanic Clouds by nearly an order of magnitude. Using main-sequence colour-magnitude limits, an observationally informed period-luminosity-colour relation for CBs, and a high morph-parameter cut ($c\geq0.7$), we identified 68 O- and B-type binaries that exhibit smooth, sinusoidal light curves with nearly equal eclipse depths. We then isolated a bona fide sample of 37 CB candidates (28 in the LMC and 9 in the SMC) that match theoretical colour-magnitude and period distributions derived from an extensive grid of MESA binary models. The bona fide sample, dominated by B-type systems with $P\approx0.6-1$ d, agrees with the predicted population and may contain many $q\approx1$ binaries, as expected from models showing mass equalization preceding temperature equalization during nuclear-timescale contact. Synthetic PHOEBE light curves of contact and near-contact phases of MESA models reveal a degeneracy between these configurations, suggesting possible misidentifications among these systems. Spectroscopic follow-up is required to test these predictions and refine the evolutionary framework of massive CBs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stellar mergers as the origin of the blue main-sequence band in young star clusters

Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that they follow a nearly flat mass function, which supports their unusual formation path. Our results imply that the cluster stars gain their mass in two different ways, by disk accretion leading to rapid rotation, contributing to the red MS, or by binary merger leading to slow rotation and populating the blue MS. We also derive the approximate merger time of the individual stars of the blue MS component, and find a strong early peak in the merger rate, with a lower level merger activity prevailing for tens of Myr. This supports recent binary formation models, and explains new velocity dispersion measurements for members of young star clusters. Our findings shed new light on the origin of the bi-modal mass, spin, and magnetic field distributions of main-sequence stars.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stripped-Envelope Stars in Different Metallicity Environments I. Evolutionary Phases, Classification and Populations

Massive stars that become stripped of their hydrogen envelope through binary interaction or winds can be observed either as Wolf-Rayet stars, if they have optically thick winds, or as transparent-wind stripped-envelope stars. We approximate their evolution through evolutionary models of single helium stars, and compute detailed model grids in the initial mass range 1.5 to 70 M$_{\odot}$ for metallicities between 0.01 and 0.04, from core helium ignition until core collapse. Throughout their lifetime, some stellar models expose the ashes of helium burning. We propose that models that have nitrogen-rich envelopes are candidate WN stars, while models with a carbon-rich surface are candidate WC stars during core helium burning, and WO stars afterwards. We measure metallicity dependance of the total lifetime of our models and the duration of their evolutionary phases. We propose an analytic estimate of the wind optical depth to distinguish models of Wolf-Rayet stars from transparent-wind stripped-envelope stars, and find that the luminosity ranges at which WN, WC and WO type stars can exist is a strong function of metallicity. We find that all carbon-rich models produced in our grids have optically thick winds and match the luminosity distribution of observed populations. We construct population models and predict the numbers of transparent-wind stripped-envelope stars and Wolf-Rayet stars, and derive their number ratios at different metallicities. We find that as metallicity increases, the number of transparent-wind stripped-envelope stars decreases and the number of Wolf-Rayet stars increases. At high metallicities WC and WO type stars become more common. We apply our population models to nearby galaxies, and find that populations are more sensitive to the transition luminosity between Wolf-Rayet stars and transparent-wind helium stars than to the metallicity dependent mass loss rates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Nature of Unseen Companions in Massive Single-Line Spectroscopic Binaries

Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The R136 star cluster dissected with Hubble Space Telescope/STIS. III. The most massive stars and their clumped winds

Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krtička & Kubát (2018) matches best the observed mass-loss rates.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bonn Optimized Stellar Tracks (BoOST). Simulated Populations of Massive and Very Massive Stars for Astrophysical Applications

Massive and very massive stars can play important roles in stellar populations by ejecting strong stellar winds and exploding in energetic phenomena. It is therefore imperative that their behavior be properly accounted for in synthetic model populations. We present nine grids of stellar evolutionary model sequences, together with finely resolved interpolated sequences and synthetic populations, of stars with 9-500 Msun and with metallicities ranging from Galactic metallicity down to 1/250 Zsun. The stellar models were computed with the Bonn evolutionary code with consistent physical ingredients, and covering core hydrogen- and core helium-burning phases. The interpolation and population synthesis were performed with our newly developed routine 'synStars'. Eight of the grids represent slowly rotating massive stars with a normal or classical evolutionary path, while one grid represents fast-rotating, chemically homogeneously evolving models. The grids contain data on stellar wind properties such as estimated wind velocity and kinetic energy of the wind, as well as common stellar parameters such as mass, radius, surface temperature, luminosity, mass-loss rate, and surface abundances of 34 isotopes. We also provide estimates of the helium and carbon-oxygen core mass for calculating the mass of stellar remnants. The Bonn Optimized Stellar Tracks (BoOST) project is published as simple tables that include stellar models, interpolated tracks, and synthetic populations. Covering the broadest mass and metallicity range of any published massive star evolutionary model sets to date, BoOST is ideal for further scientific applications such as star formation studies in both low- and high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Explodability fluctuations of massive stellar cores enable asymmetric compact object mergers such as GW190814

The first three observing runs with Advanced LIGO and Virgo have resulted in the detection of binary black hole mergers (BBH) with highly unequal mass components, which are difficult to reconcile with standard formation paradigms. The most representative of these is GW190814, a highly asymmetric merger between a 23 M$_{\odot}$ black hole and a 2.6 M$_{\odot}$ compact object. Here, we explore recent results suggesting that a sizeable fraction of stars with pre-collapse carbon-oxygen core masses above 10 M$_{\odot}$, and extending up to at least 30 M$_{\odot}$, may produce objects inside the so-called lower mass gap that bridges the division between massive pulsars and BHs in Galactic X-ray binaries. We demonstrate that such an explosion landscape would naturally cause a fraction of massive binaries to produce GW190814-like systems instead of symmetric-mass BBHs. We present examples of specific evolutionary channels leading to the formation of GW190814 and GW200210, a 24+2.8 M$_{\odot}$ merger discovered during the O3b observing run. We estimate the merger-rate density of these events in our scenario to be $\mathcal{O}$(5%) of the total BBH merger rate. Finally, we discuss the broader implications of this formation channel for compact object populations, and its possible relevance to less asymmetric merger events such as GW200105 and GW200115

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of close binary evolution on the main-sequence morphology of young star clusters

Star clusters are the building blocks of galaxies. They are composed of stars of nearly equal age and chemical composition, allowing us to use them as chronometers and as testbeds for gauging stellar evolution. It has become clear recently that massive stars are formed preferentially in close binaries, in which mass transfer will drastically change the evolution of the stars. This is expected to leave a significant imprint in the distribution of cluster stars in the Hertzsprung-Russell diagram. Our results, based on a dense model grid of more than 50,000 detailed binary-evolution calculations, indeed show several distinct, coeval main-sequence (MS) components, most notably an extended MS turnoff region, and a group of near-critical rotating stars that is spread over a large luminosity range on the red side of the classical MS. We comprehensively demonstrate the time evolution of the features in an animation, and we derive analytic expressions to describe these features. We find quantitative agreement with results based on recent photometric and spectroscopic observations. We conclude that while other factors may also be at play, binary evolution has a major impact on the MS morphology of young star clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Internal circulation in tidally locked massive binary stars -- Consequences for double black hole formation

Steady-state circulation currents are predicted in tidally deformed binary stars, which are believed to be progenitors of double black-hole merger events. This work aims to quantitatively characterise the steady-state circulation currents in components of a tidally locked binary system and to explore the effects of such currents on numerical models. Previous results describing the circulation in a single rotating star and a binary star are used to deduce a new prescription for the internal circulation in tidally locked binaries. We explore the effect of this prescription numerically for binary systems with primary masses between 25 and 100 solar masses. When comparing circulation velocities in the radial direction for the single rotating star and binary star, it is found that the average circulation velocity in the binary star may be described as an enhancement to the circulation velocity in a single rotating star. This velocity enhancement is a simple function depending on the masses of the binary components and amounts to a factor of approximately two when the components have equal masses. It is found that the ehancement causes the formation of double helium stars through efficient mixing to occur for systems with higher initial orbital periods, lower primary masses and lower mass ratios, compared to the standard circulation scenario. Taking into account appropriate distributions for primary mass, initial period and mass ratio, models with enhanced mixing predict 2.4 times more double helium stars being produced in the parameter space than models without. We conclude that the effects of companion-induced circulation have strong implications for the formation of close binary black holes. Not only do the predicted detection rates increase but double black-hole systems with mass ratios as low as 0.8 may be formed when companion-induced circulation is taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling of Magneto-Rotational Stellar Evolution I. Method and first applications

While magnetic fields have long been considered to be important for the evolution of magnetic non-degenerate stars and compact stars, it has become clear in recent years that actually all of the stars are deeply affected. This is particularly true regarding their internal angular momentum distribution, but magnetic fields may also influence internal mixing processes and even the fate of the star. We propose a new framework for stellar evolution simulations, in which the interplay between magnetic field, rotation, mass loss, and changes in the stellar density and temperature distributions are treated self-consistently. For average large-scale stellar magnetic fields which are symmetric to the axis of rotation of the star, we derive 1D evolution equations for the toroidal and poloidal components from the mean-field MHD equation by applying Alfven's theorem, and a conservative form of the angular momentum transfer due to the Lorentz force is formulated. We implement our formalism into a numerical stellar evolution code and simulate the magneto-rotational evolution of 1.5 M$_\odot$ stars. The Lorentz force aided by the $Ω$ effect imposes torsional Alfven waves propagating through the magnetized medium, leading to near-rigid rotation within the Alfven timescale. Our models with different initial spins and B-fields can reproduce the main observed properties of Ap/Bp stars. Calculations continued to the red-giant regime show a pronounced core-envelope coupling, which reproduces the core and surface rotation periods determined by asteroseismic observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Type Ia supernovae from non-accreting progenitors

Type Ia supernovae (SNe Ia) are manifestations of stars deficient of hydrogen and helium disrupting in a thermonuclear runaway. While explosions of carbon-oxygen white dwarfs are thought to account for the majority of events, part of the observed diversity may be due to varied progenitor channels. We demonstrate that helium stars with masses between $\sim$1.8 and 2.5 M$_{\odot}$ may evolve into highly degenerate, near-Chandrasekhar mass cores with helium-free envelopes that subsequently ignite carbon and oxygen explosively at densities $\sim(1.8-5.9)\times 10^{9}$g cm$^{-3}$. This happens either due to core growth from shell burning (when the core has a hybrid CO/NeO composition), or following ignition of residual carbon triggered by exothermic electron captures on $^{24}$Mg (for a NeOMg-dominated composition). We argue that the resulting thermonuclear runaways is likely to prevent core collapse, leading to the complete disruption of the star. The available nuclear energy at the onset of explosive oxygen burning suffices to create ejecta with a kinetic energy of $\sim$10$^{51}$ erg, as in typical SNe Ia. Conversely, if these runaways result in partial disruptions, the corresponding transients would resemble SN Iax events similar to SN 2002cx. If helium stars in this mass range indeed explode as SNe Ia, then the frequency of events would be comparable to the observed SN Ib/c rates, thereby sufficing to account for the majority of SNe Ia in star-forming galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Single Star Path to Be Stars

Be stars are rapidly rotating B main sequence stars, which show line emission due to an outflowing disc. By studying the evolution of rotating single star models, we can assess their contribution to the observed Be star populations. We identify the main effects which are responsible for single stars to approach critical rotation as functions of initial mass and metallicity, and predict the properties of populations of rotating single stars. We perform population synthesis with single star models of initial masses ranging between 3 and 30 solar masses, initial equatorial rotation velocities between 0 and 600 kms$^{-1}$ at compositions representing the Milky Way, Large and Small Magellanic Clouds. These models include efficient core-envelope coupling mediated by internal magnetic fields and correspond to the maximum efficiency of Be star production. We predict Be star fractions and the positions of fast rotating stars in the colour-magnitude diagram. We identify stellar wind mass-loss and the convective core mass fraction as the key parameters which determine the time dependance of the stellar rotation rates. Using empirical distributions of initial rotational velocities,our single star models can reproduce the trends observed in Be star fractions with mass and metallicity. However,they fail to produce a significant number of stars rotating very close to critical. We also find that rapidly rotating Be stars in the Magellanic Clouds should have significant surface nitrogen enrichments, which may be in conflict with abundance determinations of Be stars. Single star evolution may explain the high number of Be stars if 70 to 80% of critical rotationwould be sufficient to produce the Be phenomenon. However even in this case, the unexplained presence of many Be stars far below the cluster turn-off indicates the importance of the binary channel for Be star production.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A new route towards merging massive black holes

Recent advances in gravitational-wave astronomy make the direct detection of gravitational waves from the merger of two stellar-mass compact objects a realistic prospect. Evolutionary scenarios towards mergers of double compact objects generally invoke common-envelope evolution which is poorly understood, leading to large uncertainties in merger rates. We explore the alternative scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary which remain fully mixed due to their tidally induced high spin. We use the public stellar-evolution code MESA to systematically study this channel by means of detailed simulations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass ratios close to one, in two mass ranges, ~25...60msun and >~ 130msun, with pair instability supernovae (PISNe) being produced in-between. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for it. We map the initial parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, Kerr parameters and predict event rates. We typically find that for Z~<Z_sun/10, there is one BH+BH merger for ~1000 core-collapse supernovae. The advanced LIGO (aLIGO) detection rate is more uncertain and depends on the metallicity evolution. Deriving upper and lower limits from a local and a global approximation for the metallicity distribution of massive stars, we estimate aLIGO detection rates (at design limit) of ~19-550 yr^(-1) for BH+BH mergers below the PISN gap and of ~2.1-370 yr^(-1) above the PISN gap. Even with conservative assumptions, we find that aLIGO should soon detect BH+BH mergers from the MOB scenario and that these could be the dominant source for aLIGO detections.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Diagnostic of the unstable envelopes of Wolf-Rayet stars

The envelopes of stars near the Eddington limit are prone to various instabilities. A high Eddington factor in connection with the Fe opacity peak leads to convective instability, and a corresponding envelope inflation may induce pulsational instability. Here, we investigate the occurrence and consequences of both instabilities in models of Wolf-Rayet stars. We determine the convective velocities in the sub-surface convective zones to estimate the amplitude of the turbulent velocity at the base of the wind that potentially leads to the formation of small-scale wind structures, as observed in several WR stars. We also investigate the effect of mass loss on the pulsations of our models. We approximated solar metallicity WR stars by models of mass-losing helium stars, and we characterized the properties of convection in the envelope adopting the standard MLT. Our results show the occurrence of sub-surface convective regions in all studied models. Small surface velocity amplitudes are predicted for models with masses below 10Msun. For models with M>10Msun, the surface velocity amplitudes are of the order of 10km/s. Moreover we find the occurrence of pulsations for stars in the mass range 9-14Msun, while mass loss appears to stabilize the more massive WR stars. We confront our results with observationally derived line variabilities of 17 WN stars. The data suggest variability to occur for stars above 10Msun, which is increasing linearly with mass above this value, in agreement with our results. We further find some of our models to be unstable to radial pulsations, and predict local magnetic fields of the order of hundreds of Gauss in WR stars more massive than 10Msun. Our study relates the surface velocity fluctuations induced by sub-surface convection to the formation of clumping in the inner part of the wind. From this mechanism, we expect a stronger variability in more massive WR stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inferring supernova IIb/Ib/Ic ejecta properties from light curves and spectra: Correlations from radiative-transfer models

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875Å at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772Å line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$μ$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Metallicity dependence of turbulent pressure and macroturbulence in stellar envelopes

Macroturbulence, introduced as a fudge to reproduce the width and shape of stellar absorption lines, reflects gas motions in stellar atmospheres. While in cool stars, it is thought to be caused by convection zones immediately beneath the stellar surface, the origin of macroturbulence in hot stars is still under discussion. Recent works established a correlation between the turbulent-to-total pressure ratio inside the envelope of stellar models and the macroturbulent velocities observed in corresponding Galactic stars. To probe this connection further, we evaluated the turbulent pressure that arises in the envelope convective zones of stellar models in the mass range 1-125 Msun based on the mixing-length theory and computed for metallicities of the Large and Small Magellanic Cloud. We find that the turbulent pressure contributions in models with these metallicities located in the hot high-luminosity part of the Hertzsprung-Russel (HR) diagram is lower than in similar models with solar metallicity, whereas the turbulent pressure in low-metallicity models populating the cool part of the HR-diagram is not reduced. Based on our models, we find that the currently available observations of hot massive stars in the Magellanic Clouds appear to support a connection between macroturbulence and the turbulent pressure in stellar envelopes. Multidimensional simulations of sub-surface convection zones and a larger number of high-quality observations are necessary to test this idea more rigorously.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Models of low-mass helium white dwarfs including gravitational settling, thermal and chemical diffusion, and rotational mixing

A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study, for the first time, the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rejuvenation of stellar mergers and the origin of magnetic fields in massive stars

Approximately 10% of massive OBA main-sequence (MS) and pre-MS stars harbour strong, large-scale magnetic fields. At the same time there is a dearth of magnetic stars in close binaries. A process generating strong magnetic fields only in some stars must be responsible with the merging of pre-MS and MS stars being suggested as one such channel. Stars emerging from the coalescence of two MS stars are rejuvenated, appearing younger than they are. They can therefore be identified by comparison with reference clocks. Here we predict the rejuvenation of MS merger products over a wide range of masses and binary configurations calibrated to smoothed-particle-hydrodynamical merger models. We find that the rejuvenation is of the order of the nuclear timescale and is strongest in the lowest-mass mergers and the most evolved binary progenitors with the largest mass ratios. These predictions allow us to put constraints on the binary progenitors of merger products. We show that the magnetic stars HR 2949 and $τ$ Sco are younger than the potential binary companion HR 2948 and the Upper-Sco association, respectively, making them promising merger candidates. We find that the age discrepancies and the potential binary progenitors of both are consistent with them being rejuvenated merger products, implying that their magnetic fields may originate from this channel. Searching for age discrepancies in magnetic stars is therefore a powerful way to explore which fraction of magnetic stars may have obtained their strong magnetic field in MS mergers and to improve our understanding of magnetism in massive stars and their remnants.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cold gas in hot star clusters: the wind from the red supergiant W26 in Westerlund 1

The massive red supergiant (RSG) W26 in Westerlund 1 is one of a growing number of RSGs shown to have winds that are ionized from the outside in. The fate of this dense wind material is important for models of second generation star formation in massive star clusters. Mackey et al. (2014) showed that external photoionization can stall the wind of RSGs and accumulate mass in a dense static shell. We use 1D R-HD simulations of an externally photoionized wind to predict the Halpha and [NII] emission arising from photoionized winds both with and without a dense shell. We analyse spectra of the Halpha and [NII] emission in the environment around W26 and compare them with predicted synthetic emission. Simulations of slow winds that are decelerated into a dense shell show strongly limb-brightened line emission, with line radial velocities that are independent of the wind speed. Faster winds (>22 km/s) do not form a dense shell, have less limb-brightening, and the line radial velocity is a good tracer of the wind speed. The brightness of the [NII] and Halpha lines as a function of distance from W26 agrees reasonably well with observations when only the line flux is considered. The radial velocity disagrees, however: the brightest observed emission is blueshifted by ~25 km/s relative to the radial velocity of the star, whereas a spherically symmetric wind has the brightest emission at zero radial velocity. Our results show that the bright nebula surrounding W26 must be asymmetric; we suggest it is confined by external ram pressure from the wind of the nearby supergiant W9. We obtain a lower limit on the nitrogen abundance within the nebula of 2.35 times solar. The line ratio strongly favours photoionization over shock ionization, and so even if the observed nebula is pressure confined there should still be an ionization front and a photoionization-confined shell closer to the star.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detecting stellar-wind bubbles through infrared arcs in HII regions

Mid-infrared arcs of dust emission are often seen near ionizing stars within HII regions. A possible explanations for these arcs is that they could show the outer edges of asymmetric stellar wind bubbles. We use two-dimensional, radiation-hydrodynamics simulations of wind bubbles within HII regions around individual stars to predict the infrared emission properties of the dust within the HII region. We assume that dust and gas are dynamically well-coupled and that dust properties (composition, size distribution) are the same in the HII region as outside it, and that the wind bubble contains no dust. We post-process the simulations to make synthetic intensity maps at infrared wavebands using the TORUS code. We find that the outer edge of a wind bubble emits brightly at 24um through starlight absorbed by dust grains and re-radiated thermally in the infrared. This produces a bright arc of emission for slowly moving stars that have asymmetric wind bubbles, even for cases where there is no bow shock or any corresponding feature in tracers of gas emission. The 24um intensity decreases exponentially from the arc with increasing distance from the star because the dust temperature decreases with distance. The size distribution and composition of the dust grains has quantitative but not qualitative effects on our results. Despite the simplifications of our model, we find good qualitative agreement with observations of the HII region RCW120, and can provide physical explanations for any quantitative differences. Our model produces an infrared arc with the same shape and size as the arc around CD -38 11636 in RCW120, and with comparable brightness. This suggests that infrared arcs around O stars in HII regions may be revealing the extent of stellar wind bubbles, although we have not excluded other explanations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Envelope inflation in Wolf-Rayet stars and extended supernova shock breakout signals

Massive, luminous stars reaching the Eddington limit in their interiors develop very dilute, extended envelopes. This effect is called envelope inflation. If the progenitors of Type Ib/c supernovae, which are believed to be Wolf-Rayet (WR) stars, have inflated envelopes then the shock breakout signals diffuse in them and can extend their rise times significantly. We show that our inflated, hydrogen-free, WR stellar models with a radius of ~Rsun can have shock breakout signals longer than ~60 s. The puzzlingly long shock breakout signal observed in the Type Ib SN 2008D can be explained by an inflated progenitor envelope, and more such events might argue in favour of existence of inflated envelopes in general.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Extended supernova shock breakout signals from inflated stellar envelopes

Stars close to the Eddington luminosity can have large low-density inflated envelopes. We show that the rise times of shock breakout signals from supernovae can be extended significantly if supernova progenitors have an inflated stellar envelope. If the shock breakout occurs in such inflated envelopes, the shock breakout signals diffuse in them, and their rise time can be significantly extended. Then, the rise times of the shock breakout signals are dominated by the diffusion time in the inflated envelope rather than the light-crossing time of the progenitors. We show that our inflated Wolf-Rayet star models whose radii are of the order of the solar radius can have shock breakout signals which are longer than ~100 sec. The existence of inflated envelopes in Wolf-Rayet supernova progenitors may be related to the mysterious long shock breakout signal observed in Type Ib SN 2008D. Extended shock breakout signals may provide evidence for the existence of inflated stellar envelopes and can be used to constrain the physical properties of these enigmatic structures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Massive main sequence stars evolving at the Eddington limit

The evolution of massive stars even on the main sequence is not yet well understood. Due to the steep mass-luminosity relation, massive main sequence stars become very luminous. This brings their envelopes very close to the Eddington limit. We are analysing stellar evolutionary models in which the Eddington limit is reached and exceeded, and explore the rich diversity of physical phenomena which take place in their envelopes, and investigate their observational consequences. We use the grids of detailed stellar models by Brott et al. (2011) and Koehler et al. (2015), to investigate the envelope properties of core hydrogen burning massive stars. We find that at the stellar surface, the Eddington limit is almost never reached, even for stars up to 500 Msun. When an appropriate Eddington limit is defined locally in the stellar envelope, most stars more massive than 40 Msun actually exceed this limit, in particular in the partial ionization zones of iron, helium or hydrogen. While most models adjust their structure such that the local Eddington limit is exceeded at most by a few per cent, our most extreme models do so by a factor of more than seven. We find that the local violation of the Eddington limit has severe consequences for the envelope structure, as it leads to envelope inflation, convection, density inversions and possibly to pulsations. We find that all models with luminosities higher than 4*10^5 Lsun i.e. stars above 40 Msun show inflation, with a radius increase of up to a factor of about 40. We find that the hot edge of the S Dor variability region coincides with a line beyond which our models are inflated by more than a factor of two, indicating a possible connection between S Dor variability and inflation. Furthermore, our coolest models show highly inflated envelopes with masses of up to several solar masses, and appear to be candidates to produce major LBV eruptions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Observational consequences of turbulent pressure in the envelopes of massive stars

The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on the mixing length theory, and search for observable consequences. We find that the turbulent pressure fraction can amount up to ~5% in OB supergiants, and to ~30% in cooler supergiants. The resulting structural changes are, however, not significantly affecting the evolutionary tracks compared to previous calculations. Instead, a comparison of macroturbulent velocities derived from high quality spectra of OB stars with the turbulent pressure fraction obtained in corresponding stellar models reveals a strong correlation of these two quantities. We discuss a possible physical connection, and conclude that turbulent pressure fluctuations may drive high-order oscillations, which - as conjectured earlier - manifest themselves as macroturbulence in the photospheres of hot luminous stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Progenitors of ultra-stripped supernovae

The explosion of ultra-stripped stars in close binaries may explain new discoveries of weak and fast optical transients. We have demonstrated that helium star companions to neutron stars (NSs) may evolve into naked metal cores as low as ~1.5 Msun, barely above the Chandrasekhar mass limit, by the time they explode. Here we present a new systematic investigation of the progenitor evolution leading to such ultra-stripped supernovae (SNe), in some cases yielding pre-SN envelopes of less than 0.01 Msun. We discuss the nature of these SNe (electron-capture vs iron core-collapse) and their observational light-curve properties. Ultra-stripped SNe are highly relevant for binary pulsars, as well as gravitational wave detection of merging NSs by LIGO/VIRGO, since these events are expected to produce mainly low-kick NSs in the mass range 1.10-1.80 Msun.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radiative-transfer models for supernovae IIb/Ib/Ic from binary-star progenitors

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relating turbulent pressure and macroturbulence across the HR diagram with a possible link to γ-Dor stars

A significant fraction of the envelope of low- and intermediate-mass stars is unstable to convection, leading to sub-surface turbulent motion. Here, we consider and include the effects of turbulence pressure in our stellar evolution calculations. In search of an observational signature, we compare the fractional contribution of turbulent pressure to the observed macroturbulent velocities in stars at different evolutionary stages. We find a strong correlation between the two quantities, similar to what was previously found for massive OB stars. We therefore argue that turbulent pressure fluctuations of finite amplitude may excite high-order, high-angular degree stellar oscillations, which manifest themselves at the surface an additional broadening of the spectral lines, i.e., macroturbulence, across most of the HR diagram. When considering the locations in the HR diagram where we expect high-order oscillations to be excited by stochastic turbulent pressure fluctuations, we find a close match with the observational $γ$ Doradus instability strip, which indeed contains high-order, non-radial pulsators. We suggest that turbulent pressure fluctuations on a percentual level may contribute to the $γ$ Dor phenomenon, calling for more detailed theoretical modelling in this direction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Revealing the binary origin of Type Ic superluminous supernovae through nebular hydrogen emission

We propose that nebular H-alpha emission as detected in the Type Ic superluminous supernova iPTF13ehe stems from matter which is stripped from a companion star when the supernova ejecta collide with it. The temporal evolution, the line broadening, and the overall blueshift of the emission are consistent with this interpretation. We scale the nebular H-alpha luminosity predicted for Type Ia supernovae in single-degenerate systems to derive the stripped mass required to explain the H-alpha luminosity of iPTF13ehe. We find a stripped mass of 0.1 - 0.9 solar masses, assuming that the supernova luminosity is powered by radioactivity or magnetar spin down. Because a central heating source is required to excite the H-alpha emission, an interaction-powered model is not favored for iPTF13ehe if the H-alpha emission is from stripped matter. We derive a companion mass of more than 20 solar masses and a binary separation of less than about 20 companion radii based on the stripping efficiency during the collision, indicating that the supernova progenitor and the companion formed a massive close binary system. If Type Ic superluminous supernovae generally occur in massive close binary systems, the early brightening observed previously in several Type Ic superluminous supernovae may also be due to the collision with a close companion. Observations of nebular hydrogen emission in future Type Ic superluminous supernovae will enable us to test this interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The strange evolution of the Large Magellanic Cloud Cepheid OGLE-LMC-CEP1812

Classical Cepheids are key probes of both stellar astrophysics and cosmology as standard candles and pulsating variable stars. It is important to understand Cepheids in unprecedented detail in preparation for upcoming GAIA, JWST and extremely-large telescope observations. Cepheid eclipsing binary stars are ideal tools for achieving this goal, however there are currently only three known systems. One of those systems, OGLE-LMC-CEP1812, raises new questions about the evolution of classical Cepheids because of an apparent age discrepancy between the Cepheid and its red giant companion. We show that the Cepheid component is actually the product of a stellar merger of two main sequence stars that has since evolved across the Hertzsprung gap of the HR diagram. This post-merger product appears younger than the companion, hence the apparent age discrepancy is resolved. We discuss this idea and consequences for understanding Cepheid evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ultra-stripped supernovae: progenitors and fate

The explosion of ultra-stripped stars in close binaries can lead to ejecta masses < 0.1 M_sun and may explain some of the recent discoveries of weak and fast optical transients. In Tauris et al. (2013), it was demonstrated that helium star companions to neutron stars (NSs) may experience mass transfer and evolve into naked ~1.5 M_sun metal cores, barely above the Chandrasekhar mass limit. Here we present a systematic investigation of the progenitor evolution leading to ultra-stripped supernovae (SNe). In particular, we examine the binary parameter space leading to electron-capture (EC SNe) and iron core-collapse SNe (Fe CCSNe), respectively, and determine the amount of helium ejected with applications to their observational classification as Type Ib or Type Ic. We mainly evolve systems where the SN progenitors are helium star donors of initial mass M_He = 2.5 - 3.5 M_sun in tight binaries with orbital periods of P_orb = 0.06 - 2.0 days, and hosting an accreting NS, but we also discuss the evolution of wider systems and of both more massive and lighter - as well as single - helium stars. In some cases we are able to follow the evolution until the onset of silicon burning, just a few days prior to the SN explosion. We find that ultra-stripped SNe are possible for both EC SNe and Fe CCSNe, and that the amount of helium ejected is correlated with P_orb - the tightest systems even having donors being stripped down to envelopes of less than 0.01 M_sun. We estimate the rise time of ultra-stripped SNe to be in the range 12 hr - 8 days, and light curve decay times between 1 and 50 days. Ultra-stripped SNe may produce NSs in the mass range 1.10 - 1.80 M_sun and are highly relevant for LIGO/VIRGO since most (possibly all) merging double NS systems have evolved through this phase. Finally, we discuss the low-velocity kicks which might be imparted on these resulting NSs at birth. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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 $χ^{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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dust Production Factories in the Early Universe: Formation of Carbon Grains in Red-supergiant Winds of Very Massive Population III Stars

We investigate the formation of dust in a stellar wind during the red-supergiant (RSG) phase of a very massive Population III star with the zero-age main sequence mass of 500 M_sun. We show that, in a carbon-rich wind with a constant velocity, carbon grains can form with a lognormal-like size distribution, and that all of the carbon available for dust formation finally condense into dust for wide ranges of the mass-loss rate ((0.1-3)x10^{-3} M_sun yr^{-1}) and wind velocity (1-100 km s^{-1}). We also find that the acceleration of the wind driven by newly formed dust suppresses the grain growth but still allows more than half of gas-phase carbon to be finally locked up in dust grains. These results indicate that at most 1.7 M_sun of carbon grains can form in total during the RSG phase of 500 M_sun Population III stars. Such a high dust yield could place very massive primordial stars as important sources of dust at the very early epoch of the universe if the initial mass function of Population III stars was top-heavy. We also briefly discuss a new formation scenario of carbon-rich ultra-metal-poor stars considering the feedback from very massive Population III stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effects of stellar evolution and ionizing radiation on the environments of massive stars

We discuss two important effects for the astrospheres of runaway stars: the propagation of ionizing photons far beyond the astropause, and the rapid evolution of massive stars (and their winds) near the end of their lives. Hot stars emit ionizing photons with associated photoheating that has a significant dynamical effect on their surroundings. 3D simulations show that HII regions around runaway O stars drive expanding conical shells and leave underdense wakes in the medium they pass through. For late O stars this feedback to the interstellar medium is more important than that from stellar winds. Late in life, O stars evolve to cool red supergiants more rapidly than their environment can react, producing transient circumstellar structures such as double bow shocks. This provides an explanation for the bow shock and linear bar-shaped structure observed around Betelgeuse.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electron-capture supernovae exploding within their progenitor wind

The most massive stars on the asymptotic giant branch (AGB), so called super-AGB stars, are thought to produce supernovae (SNe) triggered by electron captures in their degenerate O+Ne+Mg cores. Super-AGB stars are expected to have slow winds with high mass-loss rates, so their wind density is high. The explosions of super-AGB stars are therefore presumed to occur in this dense wind. We provide the first synthetic light curves (LCs) for such events by exploding realistic electron-capture supernova (ecSN) progenitors within their super-AGB winds. We find that the early LC, i.e. before the recombination wave reaches the bottom of the H-rich envelope of SN ejecta (the plateau phase), is not affected by the dense wind. However, after the plateau phase, the luminosity remains higher when the super-AGB wind is taken into account. We compare our results to the historical LC of SN 1054, the progenitor of the Crab Nebula, and show that the explosion of an ecSN within an ordinary super-AGB wind can explain the LC features. We conclude that SN 1054 could have been a Type IIn SN without any extra extreme mass loss which was previously suggested to be necessary to account for its early high luminosity. We also show that our LCs match Type IIn SNe with an early plateau phase (`Type IIn-P') and suggest that they are ecSNe within super-AGB winds. Although some ecSNe can be bright in the optical spectral range due to the large progenitor radius, their X-ray luminosity from the interaction does not necessarily get as bright as other Type IIn SNe whose optical luminosities are also powered by the interaction. Thus, we suggest that optically-bright X-ray-faint Type IIn SNe can emerge from ecSNe. Optically-faint Type IIn SNe, such as SN 2008S, can also originate from ecSNe if their H-rich envelope masses are small. Some of them can be observed as `Type IIn-b' SNe due to the small H-rich envelope mass.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Explosion and nucleosynthesis of low redshift pair instability supernovae

Both recent observations and stellar evolution models suggest that pair-instability supernovae (PISNe) could occur in the local Universe, at metallicities below Z_Sun/3. Previous PISN models were mostly produced at very low metallicities in the context of the early Universe. We present new PISNe models at a metallicity of Z=0.001, which are relevant for the local Universe. We take the self-consistent stellar evolutionary models of pair-instability progenitors with initial masses of 150 and 250 solar masses at metallicity of Z=0.001 by Langer et al. (2007) and follow the evolution of these models through the supernova explosions, using a hydrodynamics stellar evolution code with an extensive nuclear network including 200 isotopes. Both models explode as PISNe without leaving a compact stellar remnant. Our models produce a nucleosynthetic pattern that is generally similar to that of Population III PISN models, which is mainly characterized by the production of large amounts of alpha-elements and a strong deficiency of the odd-charged elements. However, the odd-even effect in our models is significantly weaker than that found in Population III models. The comparison with the nucleosynthetic yields from core-collapse supernovae at a similar metallicity (Z=0.002) indicates that PISNe could have strongly influenced the chemical evolution below Z=0.002, assuming a standard initial mass function. The odd-even effect is predicted to be most prominent for the intermediate mass elements between silicon and calcium. With future observations of chemical abundances in Population II stars, our result can be used to constrain the number of PISNe that occurred during the past evolution of our Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interacting supernovae from photoionization-confined shells around red supergiant stars

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pulsations of red supergiant pair-instability supernova progenitors leading to extreme mass loss

Recent stellar evolution models show consistently that very massive metal-free stars evolve into red supergiants shortly before they explode. We argue that the envelopes of these stars, which will form pair-instability supernovae, become pulsationally unstable and that this will lead to extreme mass-loss rates despite the tiny metal content of the envelopes. We investigate the pulsational properties of such models and derive pulsationally induced mass-loss rates, which take the damping effects of the mass loss on the pulsations selfconsistently into account. We find that the pulsations may induce mass-loss rates of ~ 1e-4 - 1e-2 Msun/yr shortly before the explosions, which may create a dense circumstellar medium. Our results show that very massive stars with dense circumstellar media may stem from a wider initial mass range than pulsational-pair instability supernovae. The extreme mass loss will cease when so much of the hydrogen-rich envelope is lost that the star becomes more compact and stops pulsating. The helium core of these stars therefore remains unaffected, and their fate as pair-instability supernovae remains unaltered. The existence of dense circumstellar media around metal-free pair-instability supernovae can make them brighter and bluer, and they may be easier to detect at high redshifts than previously expected. We argue that the mass-loss enhancement in pair-instability supernova progenitors can naturally explain some observational properties of superluminous supernovae: the energetic explosions of stars within hydrogen-rich dense circumstellar media with little 56Ni production and the lack of a hydrogen-rich envelope in pair-instability supernova candidates with large 56Ni production.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sher 25: pulsating but apparently alone

The blue supergiant Sher25 is surrounded by an asymmetric, hourglass-shaped circumstellar nebula, which shows similarities to the triple-ring structure seen around SN1987A. From optical spectroscopy over six consecutive nights, we detect periodic radial velocity variations in the stellar spectrum of Sher25 with a peak-to-peak amplitude of ~12 km/s on a timescale of about 6 days, confirming the tentative detec-tion of similar variations by Hendry et al. From consideration of the amplitude and timescale of the signal, coupled with observed line profile variations, we propose that the physical origin of these variations is related to pulsations in the stellar atmosphere, rejecting the previous hypothesis of a massive, short-period binary companion. The radial velocities of two other blue supergiants with similar bipolar nebulae, SBW1 and HD 168625, were also monitored over the course of six nights, but these did not display any significant radial velocity variations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Circumstellar Medium of Massive Stars in Motion

The circumstellar medium around massive stars is strongly impacted by stellar winds, radiation, and explosions. We use numerical simulations of these interactions to constrain the current properties and evolutionary history of various stars by comparison with observed circumstellar structures. Two- and three-dimensional simulations of bow shocks around red supergiant stars have shown that Betelgeuse has probably only recently evolved from a blue supergiant to a red supergiant, and hence its bow shock is very young and has not yet reached a steady state. We have also for the first time investigated the magnetohydrodynamics of the photoionised H II region around the nearby runaway O star Zeta Oph. Finally, we have calculated a grid of models of bow shocks around main sequence and evolved massive stars that has general application to many observed bow shocks, and which forms the basis of future work to model the explosions of these stars into their pre-shaped circumstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The formation of low-mass helium white dwarfs orbiting pulsars: Evolution of low-mass X-ray binaries below the bifurcation period

Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs) which have been spun up to high rotation rates via accretion of matter from a companion star in a low-mass X-ray binary (LMXB). However, many details of this recycling scenario remain to be understood. Here we investigate binary evolution in close LMXBs to study the formation of radio MSPs with low-mass helium white dwarf companions (He WDs) in tight binaries with orbital periods P_orb = 2-9 hr. In particular, we examine: i) if such observed systems can be reproduced from theoretical modelling using standard prescriptions of orbital angular momentum losses (i.e. with respect to the nature and the strength of magnetic braking), ii) if our computations of the Roche-lobe detachments can match the observed orbital periods, and iii) if the correlation between WD mass and orbital period (M_WD, P_orb) is valid for systems with P_orb < 2 days. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in ~ 400 close LMXB systems with different initial values of donor star mass, NS mass, orbital period and the so-called gamma-index of magnetic braking. Subsequently, we followed the orbital and the interior evolution of the detached low-mass (proto) He WDs, including stages with residual shell hydrogen burning. We find that a severe fine-tuning is necessary to reproduce the observed MSPs in tight binaries with He WD companions of mass < 0.20 M_sun, which suggests that something needs to be modified or is missing in the standard input physics of LMXB modelling. We demonstrate that the theoretically calculated (M_WD, P_orb)-relation is in general also valid for systems with P_orb < 2 days, although with a large scatter in He WD masses between 0.15-0.20 M_sun. The results of the thermal evolution of the (proto) He WDs are reported in a follow-up paper (Paper II).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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%.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The timescale of low-mass proto-helium white dwarf evolution

A large number of low-mass (< 0.20 M_sun) helium white dwarfs (He WDs) have recently been discovered. The majority of these are orbiting another WD or a millisecond pulsar (MSP) in a close binary system; a few examples are found to show pulsations or to have a main-sequence star companion. There appear to be discrepancies between the current theoretical modelling of such low-mass He WDs and a number of key observed cases, indicating that their formation scenario remains to be fully understood. Here we investigate the formation of detached proto-He WDs in close-orbit low-mass X-ray binaries (LMXBs). Our prime focus is to examine the thermal evolution and the contraction phase towards the WD cooling track and investigate how this evolution depends on the WD mass. Our calculations are then compared to the most recent observational data. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in a large number of close-orbit LMXBs. Subsequently, we followed the evolution of the detached low-mass proto-He WDs, including stages with residual shell hydrogen burning and vigorous flashes caused by unstable CNO burning. We find that the time between Roche-lobe detachment until the low-mass proto-He WD reaches the WD cooling track is typically Delta_t_proto = 0.5 - 2 Gyr, depending systematically on the WD mass and therefore on its luminosity. The lowest WD mass for developing shell flashes is ~0.21 M_sun for progenitor stars of mass M2 <= 1.5 M_sun (and ~0.18 M_sun for M2 = 1.6 M_sun). The long timescale of low-mass proto-He WD evolution can explain a number of recent observations, including some MSP systems hosting He WD companions with very low surface gravities and high effective temperatures. We find no evidence for Delta_t_proto to depend on the occurrence of flashes and thus question the suggested dichotomy in thermal evolution of proto-WDs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Wind bubbles within H II regions around slowly moving stars

Interstellar bubbles around O stars are driven by a combination of the star's wind and ionizing radiation output. The wind contribution is uncertain because the boundary between the wind and interstellar medium is difficult to observe. Mid-infrared observations (e.g., of the H II region RCW 120) show arcs of dust emission around O stars, contained well within the H II region bubble. These arcs could indicate the edge of an asymmetric stellar wind bubble, distorted by density gradients and/or stellar motion. We present two-dimensional, radiation-hydrodynamics simulations investigating the evolution of wind bubbles and H II regions around massive stars moving through a dense (n=3000 cm^{-3}), uniform medium with velocities ranging from 4 to 16 km/s. The H II region morphology is strongly affected by stellar motion, as expected, but the wind bubble is also very aspherical from birth, even for the lowest space velocity considered. Wind bubbles do not fill their H II regions (we find filling factors of 10-20%), at least for a main sequence star with mass M~30 Msun. Furthermore, even for supersonic velocities the wind bow shock does not significantly trap the ionization front. X-ray emission from the wind bubble is soft, faint, and comes mainly from the turbulent mixing layer between the wind bubble and the H II region. The wind bubble radiates <1 per cent of its energy in X-rays; it loses most of its energy by turbulent mixing with cooler photoionized gas. Comparison of the simulations with the H II region RCW 120 shows that its dynamical age is <=0.4 Myr and that stellar motion <=4 km/s is allowed, implying that the ionizing source is unlikely to be a runaway star but more likely formed in situ. The region's youth, and apparent isolation from other O or B stars, makes it very interesting for studies of massive star formation and of initial mass functions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

3D Smoothed Particle Hydrodynamics Models of Betelgeuse's Bow Shock

Betelgeuse, the bright red supergiant (RSG) in Orion, is a runaway star. Its supersonic motion through the interstellar medium has resulted in the formation of a bow shock, a cometary structure pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. We show that the bow shock morphology depends substantially on the growth timescale for Rayleigh-Taylor versus Kelvin-Helmholtz instabilities. We discuss our models in light of the recent Herschel, GALEX and VLA observations. If the mass in the bow shock shell is low (~few x 0.001 Msun), as seems to be implied by the AKARI and Herschel observations, then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular, smooth bow shock shell is consistent with this conclusion. We further discuss the implications of our results, in particular, the possibility that Betelgeuse may have only recently entered the RSG phase.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Massive Pulsar in a Compact Relativistic Binary

Many physically motivated extensions to general relativity (GR) predict significant deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 +/- 0.04 solar mass pulsar in a 2.46-hr orbit with a 0.172 +/- 0.003 solar mass white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of H II regions around exiled O stars

At least 25 per cent of massive stars are ejected from their parent cluster, becoming runaways or exiles, travelling with often-supersonic space velocities through the interstellar medium (ISM). Their overpressurised H II regions impart kinetic energy and momentum to the ISM, compress and/or evaporate dense clouds, and can constrain properties of both the star and the ISM. Here we present one-, two-, and (the first) three-dimensional simulations of the H II region around a massive star moving supersonically through a uniform, magnetised ISM, with properties appropriate for the nearby O star Zeta Oph. The H II region leaves an expanding overdense shell behind the star and, inside this, an underdense wake that should be filled with hot gas from the shocked stellar wind. The gas column density in the shell is strongly influenced by the ISM magnetic field strength and orientation. H-alpha emission maps show the H II region remains roughly circular, although the star is displaced somewhat from the centre of emission. For our model parameters, the kinetic energy feedback from the H II region is comparable to the mechanical luminosity of the stellar wind, and the momentum feedback rate is >100X larger than that from the wind and about 10X larger than the total momentum input rate available from radiation pressure. Compared to the star's eventual supernova explosion, the kinetic energy feedback from the H II region over the star's main sequence lifetime is >100X less, but the momentum feedback is up to 4X larger. H II region dynamics are found to have only a small effect on the ISM conditions that a bow shock close to the star would encounter.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution towards and beyond accretion-induced collapse of massive white dwarfs and formation of millisecond pulsars

Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs), formed via type Ib/c core-collapse supernovae (SNe), which have been spun up to high rotation rates via accretion from a companion star in a low-mass X-ray binary (LMXB). In an alternative formation channel, NSs are produced via the accretion-induced collapse (AIC) of a massive white dwarf (WD) in a close binary. Here we investigate binary evolution leading to AIC and examine if NSs formed in this way can subsequently be recycled to form MSPs and, if so, how they can observationally be distinguished from pulsars formed via the standard core-collapse SN channel in terms of their masses, spins, orbital periods and space velocities. Numerical calculations with a detailed stellar evolution code were used for the first time to study the combined pre- and post-AIC evolution of close binaries. We investigated the mass transfer onto a massive WD in 240 systems with three different types of non-degenerate donor stars: main-sequence stars, red giants, and helium stars. When the WD is able to accrete sufficient mass (depending on the mass-transfer rate and the duration of the accretion phase) we assumed it collapses to form a NS and we studied the dynamical effects of this implosion on the binary orbit. Subsequently, we followed the mass-transfer epoch which resumes once the donor star refills its Roche lobe and calculated the continued LMXB evolution until the end. We demonstrate that the final properties of these MSPs are, in general, remarkably similar to those of MSPs formed via the standard core-collapse SN channel. However, the resultant MSPs created via the AIC channel preferentially form in certain orbital period intervals. Finally, we discuss the link between AIC and young NSs in globular clusters. Our calculations are also applicable to progenitor binaries of SNe Ia under certain conditions. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Massive Binary Stars and Self-Enrichment of Globular Clusters

Globular clusters contain many stars with surface abundance patterns indicating contributions from hydrogen burning products, as seen in the anti-correlated elemental abundances of e.g. sodium and oxygen, and magnesium and aluminium. Multiple generations of stars can explain this phenomenon, with the second generation forming from a mixture of pristine gas and ejecta from the first generation. We show that massive binary stars may be a source of much of the material that makes this second generation of stars. Mass transfer in binaries is often non-conservative and the ejected matter moves slowly enough that it can remain inside a globular cluster and remain available for subsequent star formation. Recent studies show that there are more short-period massive binaries than previously thought, hence also more stars that interact and eject nuclear-processed material.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical models for the circumstellar medium around Betelgeuse

The nearby red supergiant (RSG) Betelgeuse has a complex circumstellar medium out to at least 0.5 parsecs from its surface, shaped by its mass-loss history within the past 0.1 Myr, its environment, and its motion through the interstellar medium (ISM). In principle its mass-loss history can be constrained by comparing hydrodynamic models with observations. Observations and numerical simulations indicate that Betelgeuse has a very young bow shock, hence the star may have only recently become a RSG. To test this possibility we calculated a stellar evolution model for a single star with properties consistent with Betelgeuse. We incorporated the resulting evolving stellar wind into 2D hydrodynamic simulations to model a runaway blue supergiant (BSG) undergoing the transition to a RSG near the end of its life. The collapsing BSG wind bubble induces a bow shock-shaped inner shell which at least superficially resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We investigate whether this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The ages of young star clusters, massive blue stragglers and the upper mass limit of stars: analysing age dependent stellar mass functions

Massive stars rapidly change their masses through strong stellar winds and mass transfer in binary systems. We show that such mass changes leave characteristic signatures in stellar mass functions of young star clusters which can be used to infer their ages and to identify products of binary evolution. We model the observed present day mass functions of the young Galactic Arches and Quintuplet star clusters using our rapid binary evolution code. We find that shaping of the mass function by stellar wind mass loss allows us to determine the cluster ages to 3.5$\pm$0.7 Myr and 4.8$\pm$1.1 Myr, respectively. Exploiting the effects of binary mass exchange on the cluster mass function, we find that the most massive stars in both clusters are rejuvenated products of binary mass transfer, i.e. the massive counterpart of classical blue straggler stars. This resolves the problem of an apparent age spread among the most luminous stars exceeding the expected duration of star formation in these clusters. We perform Monte Carlo simulations to probe stochastic sampling, which support the idea of the most massive stars being rejuvenated binary products. We find that the most massive star is expected to be a binary product after 1.0$\pm$0.7 Myr in Arches and after 1.7$\pm$1.0 Myr in Quintuplet. Today, the most massive 9$\pm$3 stars in Arches and 8$\pm$3 in Quintuplet are expected to be such objects. Our findings have strong implications for the stellar upper mass limit and solve the discrepancy between the claimed 150 $\mathrm{M}_\odot$ limit and observations of fours stars with initial masses of 165-320 $\mathrm{M}_\odot$ in R136 and of SN 2007bi, which is thought to be a pair-instability supernova from an initial 250 $\mathrm{M}_\odot$ star. Using the stellar population of R136, we revise the upper mass limit to values in the range 200-500 $\mathrm{M}_\odot$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Classical Cepheids Require Enhanced Mass Loss

Measurements of rates of period change of Classical Cepheids probe stellar physics and evolution. Additionally, better understanding of Cepheid structure and evolution provides greater insight into their use as standard candles and tools for measuring the Hubble constant. Our recent study of the period change of the nearest Cepheid, Polaris, suggested that it is undergoing enhanced mass loss when compared to canonical stellar evolution model predictions. In this work, we expand the analysis to rates of period change measured for about 200 Galactic Cepheids and compare them to population synthesis models of Cepheids including convective core overshooting and enhanced mass loss. Rates of period change predicted from stellar evolution models without mass loss do not agree with observed rates whereas including enhanced mass loss yields predicted rates in better agreement with observations. This is the first evidence that enhanced mass loss as suggested previously for Polaris and delta Cephei must be a ubiquitous property of Classical Cepheids.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Double bow shocks around young, runaway red supergiants: application to Betelgeuse

A significant fraction of massive stars are moving supersonically through the interstellar medium (ISM), either due to disruption of a binary system or ejection from their parent star cluster. The interaction of their wind with the ISM produces a bow shock. In late evolutionary stages these stars may undergo rapid transitions from red to blue and vice versa on the Hertzsprung-Russell diagram, with accompanying rapid changes to their stellar winds and bow shocks. Recent 3D simulations of the bow shock produced by the nearby runaway red supergiant (RSG) Betelgeuse, under the assumption of a constant wind, indicate that the bow shock is very young (<30000 years old), hence Betelgeuse may have only recently become a RSG. To test this possibility, we have calculated stellar evolution models for single stars which match the observed properties of Betelgeuse in the RSG phase. The resulting evolving stellar wind is incorporated into 2D hydrodynamic simulations in which we model a runaway blue supergiant (BSG) as it undergoes the transition to a RSG near the end of its life. We find that the collapsing BSG wind bubble induces a bow shock-shaped inner shell around the RSG wind that resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We suggest that this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Formation of millisecond pulsars with CO white dwarf companions - II. Accretion, spin-up, true ages and comparison to MSPs with He white dwarf companions

Millisecond pulsars (MSPs) are mainly characterised by their spin periods, B-fields and masses - quantities which are largely affected by previous interactions with a companion star in a binary system. In this paper, we investigate the formation mechanism of MSPs by considering the pulsar recycling process in both intermediate-mass X-ray binaries (IMXBs) and low-mass X-ray binaries (LMXBs). The IMXBs mainly lead to the formation of binary MSPs with a massive carbon-oxygen (CO) or an oxygen-neon-magnesium white dwarf (ONeMg WD) companion, whereas the LMXBs form recycled pulsars with a helium white dwarf (He WD) companion. We discuss the accretion physics leading to the spin-up line in the PPdot-diagram and demonstrate that such a line cannot be uniquely defined. We derive a simple expression for the amount of accreted mass needed for any given pulsar to achieve its equilibrium spin and apply this to explain the observed differences of the spin distributions of recycled pulsars with different types of companions. From numerical calculations we present further evidence for significant loss of rotational energy in accreting X-ray MSPs in LMXBs during the Roche-lobe decoupling phase (Tauris 2012) and demonstrate that the same effect is negligible in IMXBs. We examine the recycling of pulsars with CO WD companions via Case BB Roche-lobe overflow (RLO) of naked helium stars in post common envelope binaries. We find that such pulsars typically accrete of the order 0.002-0.007 M_sun which is just about sufficient to explain their observed spin periods. We introduce isochrones of radio MSPs in the PPdot-diagram to follow their spin evolution and discuss their true ages from comparison with observations. Finally, we apply our results of the spin-up process to the massive pulsar J1614-2230 (Paper I) and put new constraints on the birth masses of a number of recycled pulsars. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The period change of the Cepheid Polaris suggests enhanced mass loss

Polaris is one of the most observed stars in the night sky, with recorded observations spanning more than 200 years. From these observations, one can study the real-time evolution of Polaris via the secular rate of change of the pulsation period. However, the measurements of the rate of period change do not agree with predictions from state-of-the-art stellar evolution models. We show that this may imply that Polaris is currently losing mass at a rate of $\dot{M} \approx 10^{-6} M_\odot$ yr$^{-1}$ based on the difference between modeled and observed rates of period change, consistent with pulsation-enhanced Cepheid mass loss. A relation between the rate of period change and mass loss has important implications for understanding stellar evolution and pulsation, and provides insight into the current Cepheid mass discrepancy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Core-collapse explosions of Wolf-Rayet stars and the connection to type IIb/Ib/Ic supernovae

We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Formation of millisecond pulsars with CO white dwarf companions - I. PSR J1614-2230: Evidence for a neutron star born massive

The recent discovery of a 2 M_sun binary millisecond pulsar (Demorest et al. 2010) has not only important consequences for the equation-of-state of nuclear matter at high densities but also raises the interesting question if the neutron star PSR J1614-2230 was born massive. The answer is vital for understanding neutron star formation in core collapse supernovae. Furthermore, this system raises interesting issues about the nature of the progenitor binary and how it evolved during its mass exchanging X-ray phase. In this paper we discuss the progenitor evolution of PSR J1614-2230. We have performed detailed stellar evolution modelling of intermediate-mass X-ray binaries undergoing Case A Roche-lobe overflow (RLO) and applied an analytic parameterization for calculating the outcome of either a common envelope evolution or the highly super-Eddington isotropic re-emission mode. We find two viable possibilities for the formation of the PSR J1614-2230 system: either it contained a 2.2-2.6 M_sun giant donor star and evolved through a common envelope and spiral-in phase or, more likely, it descended from a close binary system with a 4.0-5.0 M_sun main sequence donor star via Case A RLO. We conclude that the neutron star must have been born with a mass of ~1.95 M_sun or 1.7+-0.15 M_sun, respectively - which significantly exceeds neutron star birth masses in previously discovered radio pulsar systems. Based on the expected neutron star birth masses from considerations of stellar evolution and explosion models, we find that the progenitor star of PSR J1614-2230 is likely to have been more massive than 20 M_sun.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Is there a mass discrepancy in the Cepheid binary OGLE-LMC-CEP0227?

Context. The Cepheid mass discrepancy, the difference between masses predicted from stellar evolution and stellar pulsation calculations, is a challenge for the understanding of stellar astrophysics. Recent models of the eclipsing binary Cepheid OGLE-LMC-CEP-0227 have suggested that the discrepancy may be resolved. Aims. We explore for what physical parameters do stellar evolution models agree with the measured properties of OGLE-LMC-CEP0227 and compare to canonical stellar evolution models assuming no convective core overshooting. Methods. We construct state-of-the-art stellar evolution models for varying mass, metallicity, and convective core overshooting and compare the stellar evolution predictions with the observed properties. Results. The observed mass, effective temperature, and radius of the two stars in the binary system are well fit by numerous combinations of physical parameters, suggesting a Cepheid mass discrepancy of 10-20% relative to canonical stellar evolution models. Conclusions. The properties of the observed binary Cepheid suggest that the Cepheid mass discrepancy is still a challenge and requires more specific observations, such as the rate of period change, to better constrain and understand the necessary physics for stellar evolution models to resolve the discrepancy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the circumstellar medium of massive stars and how it may appear in GRB observations

Massive stars lose a large fraction of their original mass over the course of their evolution. These stellar winds shape the surrounding medium according to parameters that are the result of the characteristics of the stars, varying over time as the stars evolve, leading to both permanent and temporary features that can be used to constrain the evolution of the progenitor star. Because long Gamma-Ray Bursts (GRBs) are thought to originate from massive stars, the characteristics of the circumstellar medium (CSM) should be observable in the signal of GRBs. This can occur directly, as the characteristics of the GRB-jet are changed by the medium it collides with, and indirectly because the GRB can only be observed through the extended circumstellar bubble that surrounds each massive star. We use computer simulations to describe the circumstellar features that can be found in the vicinity of massive stars and discuss if, and how, they may appear in GRB observations. Specifically, we make hydrodynamical models of the circumstellar environment of a rapidly rotating, chemically near-homogeneous star, which is represents a GRB progenitor candidate model. The simulations show that the star creates a large scale bubble of shocked wind material, which sweeps up the interstellar medium in an expanding shell. Within this bubble, temporary circumstellar shells, clumps and voids are created as a result of changes in the stellar wind. Most of these temporary features have disappeared by the time the star reaches the end of its life, leaving a highly turbulent circumstellar bubble behind. Placing the same star in a high density environment simplifies the evolution of the CSM as the more confined bubble prohibits the formation of some of the temporary structures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rotating Massive Main-Sequence Stars I: Grids of Evolutionary Models and Isochrones

We present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. Each grid covers masses ranging from 5 to 60 Msun and initial rotation rates between 0 and about 600 km/s. To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective overshooting by using the observed drop in rotation rates for stars with surface gravities log g < 3.2 to determine the width of the main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in the Large Magellanic cloud. We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the most massive fast rotators. We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that the consideration of surface abundances can resolve this degeneracy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rotating Massive Main-Sequence Stars II: Simulating a Population of LMC early B-type Stars as a Test of Rotational Mixing

Rotational mixing in massive stars is a widely applied concept, with far reaching consequences for stellar evolution. Nitrogen surface abundances for a large and homogeneous sample of massive B-type stars in the LMC were obtained by the VLT-FLAMES Survey of Massive Stars. This sample is the first covering a broad range of projected stellar rotational velocities, with a large enough sample of high quality data to allow for a statistically significant analysis. We use the sample to provide the first rigorous test of the theory of rotational mixing in massive stars. We calculated a grid of stellar evolution models, using the FLAMES sample to calibrate some of the uncertain mixing processes. We developed a new population-synthesis code, which uses this grid to simulate a large population of stars with masses, ages and rotational velocity distributions consistent with those from the FLAMES sample. The synthesized population is then filtered by the selection effects in the observed sample, to enable a direct comparison between the empirical results and theoretical predictions. Our simulations reproduce the fraction of stars without significant nitrogen enrichment. The predicted number of rapid rotators with enhanced nitrogen is about twice as large as found observationally. Furthermore, a group of stars consisting of slowly rotating, nitrogen-enriched objects and another consisting of rapidly rotating un-enriched objects can not be reproduced by our single-star population synthesis. Additional physical processes appear to be required to understand the population of massive main-sequence stars from the FLAMES sample.We discuss the possible role of binary stars and magnetic fields in the interpretation of our results. We find that the population of slowly rotating nitrogen-enriched stars is unlikely produced via mass transfer and subsequent tidal spin-down in close binary systems

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Runaway stars as progenitors of supernovae and gamma-ray bursts

When a core collapse supernova occurs in a binary system, the surviving star as well as the compact remnant emerging from the SN, may reach a substantial space velocity. With binary population synthesis modelling at solar and one fifth of solar metallicity, we predict the velocities of such runaway stars or binaries. We compile predictions for runaway OB stars, red supergiants and Wolf-Rayet stars. For those stars or binaries which undergo a second stellar explosion we compute their further evolution and the distance travelled until a Type II or Type Ibc SN or a long or short gamma-ray burst occurs. We find our predicted population of OB runaway stars broadly matches the observed population of stars but, to match the fastest observed WR runaway stars, we require that black holes receive an asymmetric kick upon formation. We find that at solar metallicity Type Ic SN progenitors travel shorter distances than the progenitors of other SN types because they are typically more massive and thus have shorter lifetimes. Those of Type IIP SN can fly farthest about 48 pc on average at solar metallicity. In considering the consequences of assuming that the progenitors of long GRBs are spun-up secondary stars that experience quasi-homogeneous evolution, we find that such evolution has a dramatic effect on the population of runaway WR stars and that some 30 per cent of GRBs could occur a hundred parsecs or more from their initial positions. We also consider mergers of double compact object binaries consisting of neutron stars and/or black holes. We find the most common type of visible mergers are neutron star--black hole mergers that are roughly ten times more common than neutron star--neutron star mergers. We also find that there may be a population of low-velocity neutron stars that are ejected from a binary rather than by their own natal kick.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stellar Evolution with Pulsation-Driven Mass Loss: The Case of LMC Cepheids

The Cepheid mass discrepancy, the difference between mass estimates from stellar pulsation and stellar evolution models, is a long standing challenge for the understanding of stellar astrophysics. We discuss the current state of the mass discrepancy and test the role of pulsation-driven mass loss in state-of-the-art stellar evolution calculations of Large Magellanic Cepheids in resolving it. We find that Cepheid mass loss is a significant contributor to the mass discrepancy, but it is not clear if the metallicity dependence of Cepheid mass loss is consistent with the measured metallicity dependence of the mass discrepancy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Cepheid mass discrepancy and pulsation-driven mass loss

Context. A longstanding challenge for understanding classical Cepheids is the Cepheid mass discrepancy, where theoretical mass estimates using stellar evolution and stellar pulsation calculations have been found to differ by approximately 10 - 20%. Aims. We study the role of pulsation-driven mass loss during the Cepheid stage of evolution as a possible solution to this mass discrepancy. Methods. We computed stellar evolution models with a Cepheid mass-loss prescription and various amounts of convective core overshooting. The contribution of mass loss towards the mass discrepancy is determined using these models, Results. Pulsation-driven mass loss is found to trap Cepheid evolution on the instability strip, allowing them to lose about 5 - 10% of their total mass when moderate convective core overshooting, an amount consistent with observations of other stars, is included in the stellar models. Conclusions. We find that the combination of moderate convective core overshooting and pulsation-driven mass loss can solve the Cepheid mass discrepancy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convective Core Overshoot & Mass Loss in Classical Cepheids: A Solution to the Mass Discrepancy?

We explore the role of mass loss and convective core overshoot in the evolution of Classical Cepheids. Stellar evolution models are computed with a recipe for pulsation-driven mass loss and it is found that mass loss alone is unable to account for the long-standing Cepheid mass discrepancy. However, the combination of mass loss and moderate convective core overshooting does provide a solution, bringing the amount of convective core overshooting in Cepheids closer to that found for other stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermohaline mixing in evolved low-mass stars

Thermohaline mixing has recently been proposed to occur in low-mass red giants, with large consequence for the chemical yields of low-mass stars. We investigate the role of thermohaline mixing during the evolution of stars between 1Msun and 3Msun, in comparison to other mixing processes acting in these stars. We use a stellar evolution code which includes rotational mixing, internal magnetic fields and thermohaline mixing. We confirm that during the red giant stage, thermohaline mixing has the potential to decrease the abundance of ^3He which is produced earlier on the main sequence. In our models we find that this process is working on the RGB only in stars with initial mass M \simle 1.5Msun. Moreover we report that thermohaline mixing is present also during core helium burning and beyond, and has the potential to change the surface abundances of AGB stars. While we find rotational and magnetic mixing to be negligible compared to the thermohaline mixing in the relevant layers, the interaction of thermohaline motions with the differential rotation may be essential to establish the time scale of thermohaline mixing in red giants. To explain the surface abundances observed at the bump in the luminosity function, the speed of the mixing process needs to be more than two orders of magnitude higher than in our models. However it is not clear if thermohaline mixing is the only physical process responsible for these surface abundance anomalies. Therefore, at this stage, it is not possible to calibrate the efficiency of thermohaline mixing against the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Turbulence and magnetic spots at the surface of hot massive stars

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred as the "iron convection zone" (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Type Ib/c supernovae in binary systems I. Evolution and properties of the progenitor stars

We investigate the evolution of Type Ib/c supernova (SN Ib/c) progenitors in close binary systems, using new evolutionary models that include the effects of rotation, with initial masses of 12 - 25 Msun for the primary components, and of single helium stars with initial masses of 2.8 - 20 Msun. We find that, despite the impact of tidal interaction on the rotation of primary stars, the amount of angular momentum retained in the core at the presupernova stage in different binary model sequences converge to a value similar to those found in previous single star models. This amount is large enough to produce millisecond pulsars, but too small to produce magnetars or long gamma-ray bursts. We employ the most up-to-date estimate for the Wolf-Rayet mass loss rate, and its implications for SN Ib/c progenitors are discussed in detail. In terms of stellar structure, SN Ib/c progenitors in binary systems are predicted to have a wide range of final masses even up to 7 Msun, with helium envelopes of 0.16 - 1.5 Msun. Our results indicate that, if the lack of helium lines in the spectra of SNe Ic were due to small amounts of helium, the distribution of both initial and final masses of SN Ic progenitors should be bimodal. Furthermore, we find that a thin hydrogen layer (0.001 - 0.01 Msun) is expected to be present in many SN Ib progenitors at the presupernova stage. We show that the presence of hydrogen, together with a rather thick helium envelope, can lead to a significant expansion of some SN Ib/c progenitors by the time of supernova explosion. This may have important consequences for the shock break-out and supernova light curve. We also argue that some SN progenitors with thin hydrogen layers produced via Case AB/B transfer might be related to Type IIb supernova progenitors with relatively small radii of about 10 Rsun.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Are the Models for Type Ia Supernova Progenitors Consistent with the Properties of Supernova Remnants?

We explore the relationship between the models for progenitor systems of Type Ia supernovae and the properties of the supernova remnants that evolve after the explosion. Most models for Type Ia progenitors in the single degenerate scenario predict substantial outflows during the presupernova evolution. Expanding on previous work, we estimate the imprint of these outflows on the structure of the circumstellar medium at the time of the supernova explosion, and the effect that this modified circumstellar medium has on the evolution of the ensuing supernova remnant. We compare our simulations with the observational properties of known Type Ia supernova remnants in the Galaxy (Kepler, Tycho, SN 1006), the Large Magellanic Cloud (0509-67.5, 0519-69.0, N103B), and M31 (SN 1885). We find that optically thick outflows from the white dwarf surface (sometimes known as accretion winds) with velocities above 200 km/s excavate large low-density cavities around the progenitors. Such large cavities are incompatible with the dynamics of the forward shock and the X-ray emission from the shocked ejecta in all the Type Ia remnants that we have examined.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Evolution of rapidly rotating metal-poor massive stars towards gamma-ray bursts

Recent models of rotating massive stars including magnetic fields prove it difficult for the cores of single stars to retain enough angular momentum to produce a collapsar and gamma-ray burst. At low metallicity, even very massive stars may retain a massive hydrogen envelope due to the weakness of the stellar winds, posing an additional obstacle to the collapsar model. Here, we consider the evolution of massive, magnetic stars where rapid rotation induces almost chemically homogeneous evolution. We find that in this case, the requirements of the collapsar model are rather easily fulfilled if the metallicity sufficiently small: 1) Rapidly rotating helium stars are formed without the need to remove the hydrogen envelope, avoiding mass-loss induced spin-down. 2) Angular momentum transport from the helium core to hydrogen envelope by magnetic torques is insignificant. We demonstrate this by calculating evolutionary models of massive stars with various metallicities, and derive an upper metallicity limit for this scenario based on currently proposed mass loss rates. Our models also suggest the existence of a lower CO-core mass limit of about 10 Msun -- which relates to an initial mass of only about 20 Msun within our scenario -- for GRB production. We argue that the relative importance of the considered GRB progenitor channel, compared to any channel related to binary stars, may increase with decreasing metallicity, and that this channel might be the major path to GRBs from first stars.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Chemical composition and origin of nebulae around Luminous Blue Variables

We use the analysis of the heavy element abundances (C, N, O, S) in circumstellar nebulae around Luminous Blue Variables to infer the evolutionary phase in which the material has been ejected. (1) We discuss the different effects that may have changed the gas composition of the nebula since it was ejected (2) We calculate the expected abundance changes at the stellar surface due to envelope convection in the red supergiant phase. If the observed LBV nebulae are ejected during the RSG phase, the abundances of the LBV nebulae require a significantly smaller amount of mass to be lost than assumed in evolutionary models. (3) We calculate the changes in the surface composition during the main sequence phase by rotation induced mixing. If the nebulae are ejected at the end of the MS-phase, the abundances in LBV nebulae are compatible with mixing times between 5 x 10^6 and 1 x 10^7 years. The existence of ON stars supports this scenario. (4) The predicted He/H ratio in the nebulae are significantly smaller than the current observed photospheric values of their central stars. Combining various arguments we show that the LBV nebulae are ejected during the blue SG phase and that the stars have not gone through a RSG phase. The chemical enhancements are due to rotation induced mixing, and the ejection is possibly triggered by near-critical rotation. During the ejection, the outflow was optically thick, which resulted in a large effective radius and a low effective temperature. This also explains the observed properties of LBV dust.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Pulsations in red supergiants with high L/M ratio -- Implications for the stellar and circumstellar structure of supernova progenitors

We investigate the pulsational properties of RSG models --- which we evolve from ZAMS masses in the range 10 to $20 \Msun$ --- by means of linear and non-linear calculations. We find period and growth rate of the dominant fundamental mode to increase with increasing luminosity-to-mass ratio $L/M$. Our models obtain relatively large $L/M$ values due to the inclusion of rotation in the evolutionary calculations; however, the largest values are obtained at and beyond central He-exhaustion due to major internal rearrangements of the nuclear burning regions. Our non-linear calculations as well as the behavior of the linear period and growth rate of the pulsations for periods approaching the Kelvin-Helmholtz time scale of the H-rich stellar envelope point towards the possibility of large amplitude pulsations. Such properties are similar to that found in AGB stars and suggest the possibility of a ``superwind'' to occur before the RSGs explode as supernovae. We conclude that changes in global stellar properties during the last few $10^4 \yr$ before core collapse may lead to drastic changes in the pulsational and wind properties of pre-supernova stars, with marked consequences for the immediate pre-supernova structure of the star and the circumstellar medium. We compare our results with observations of long-period OH/IR variables and discuss observational evidence for our scenario from observed supernova light curves, spectra and remnants.

0 citations0 reviewsNo code linkedOpen access