BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Bill Paxton

Bill Paxton contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.HEastro-ph.GAastro-ph.IMastro-phastro-ph.COhep-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
24works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

24 published item(s)

preprint2020arXiv

A Massive Star's Dying Breaths: Pulsating Red Supergiants and Their Resulting Type IIP Supernovae

Massive stars undergo fundamental-mode and first-overtone radial pulsations with periods of 100-1000 days as Red Supergiants (RSGs). At large amplitudes, these pulsations substantially modify the outer envelope's density structure encountered by the outgoing shock wave from the eventual core collapse of these $M>9M_\odot$ stars. Using Modules for Experiments in Stellar Astrophysics (MESA), we model the effects of fundamental-mode and first-overtone pulsations in the RSG envelopes, and the resulting Type IIP supernovae (SNe) using MESA+STELLA. We find that, in the case of fundamental mode pulsations, SN plateau observables such as the luminosity at day 50, $L_{50}$, time-integrated shock energy $ET$, and plateau duration $t_{\rm p}$ are consistent with radial scalings derived considering explosions of non-pulsating stars. Namely, most of the effect of the pulsation is consistent with the behavior expected for a star of a different size at the time of explosion. However, in the case of overtone pulsations, the Lagrangian displacement is not monotonic. Therefore, in such cases, excessively bright or faint SN emission at different times reflects the underdense or overdense structure of the emitting region near the SN photosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Carbon Shell or Core Ignitions in White Dwarfs Accreting from Helium Stars

White dwarfs accreting from helium stars can stably burn at the accreted rate and avoid the challenge of mass loss associated with unstable Helium burning that is a concern for many Type Ia supernovae scenarios. We study binaries with helium stars of mass $1.25 M_\odot\le M_{\rm{He}} \le 1.8 M_\odot$, which have lost their hydrogen rich envelopes in an earlier common envelope event and now orbit with periods ($P_{\rm orb}$) of several hours with non-rotating $0.84$ and $1.0 M_\odot$ C/O WDs. The helium stars fill their Roche lobes (RLs) after exhaustion of central helium and donate helium on their thermal timescales (${\sim}10^5$yr). As shown by others, these mass transfer rates coincide with the steady helium burning range for WDs, and grow the WD core up to near the Chandrasekhar mass ($M_{\rm Ch}$) and a core carbon ignition. We show here, however, that many of these scenarios lead to an ignition of hot carbon ashes near the outer edge of the WD and an inward going carbon flame that does not cause an explosive outcome. For $P_{\rm orb} = 3$ hours, $1.0 M_\odot$ C/O WDs with donor masses $M_{\rm He}\gtrsim1.8 M_\odot$ experience a shell carbon ignition, while $M_{\rm He}\lesssim1.3 M_\odot$ will fall below the steady helium burning range and undergo helium flashes before reaching core C ignition. Those with $1.3 M_\odot \lesssim M_{\rm He} \lesssim 1.7 M_\odot$ will experience a core C ignition. We also calculate the retention fraction of accreted helium when the accretion rate leads to recurrent weak helium flashes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

AM Canum Venaticorum Progenitors with Helium Star Donors and the resultant Explosions

We explore the outcome of mass transfer via Roche lobe overflow (RLOF) of $M_{\rm He}\lesssim0.51 M_\odot$ pure helium burning stars in close binaries with white dwarfs (WDs). The evolution is driven by the loss of angular momentum through gravitational wave radiation (GWR), and both stars are modeled using Modules for Experiments in Stellar Astrophysics (MESA). The donors have masses of $M_{\rm He}=0.35, 0.4, \&\ 0.51M_\odot$ and accrete onto WDs of mass $M_{\rm WD}$ from $0.6M_\odot$ to $1.26M_\odot$. The initial orbital periods ($P_{\rm{orb}}$) span 20 to 80 minutes. For all cases, the accretion rate onto the WD is below the stable helium burning range, leading to accumulation of helium followed by unstable ignition. The mass of the convective core in the donors is small enough so that the WD accretes enough helium-rich matter to undergo a thermonuclear runaway in the helium shell before any carbon-oxygen enriched matter is transferred. The mass of the accumulated helium shell depends on $M_{\rm WD}$ and the accretion rate. We show that for $M_{\rm He}\gtrsim0.4 M_\odot$ and $M_{\rm WD}\gtrsim0.8 M_\odot$, the first flash is likely vigorous enough to trigger a detonation in the helium layer. These thermonuclear runaways may be observed as either faint and fast .Ia SNe, or, if the carbon in the core is also detonated, Type Ia SNe. Those that survive the first flash and eject mass will have a temporary increase in orbital separation, but GWR drives the donor back into contact, resuming mass transfer and triggering several subsequent weaker flashes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dark Stars: Improved Models and First Pulsation Results

We use the stellar evolution code MESA to study dark stars. Dark stars (DSs), which are powered by dark matter (DM) self-annihilation rather than by nuclear fusion, may be the first stars to form in the Universe. We compute stellar models for accreting DSs with masses up to 10^6 M_{sun}. The heating due to DM annihilation is self-consistently included, assuming extended adiabatic contraction of DM within the minihalos in which DSs form. We find remarkably good overall agreement with previous models, which assumed polytropic interiors. There are some differences in the details, with positive implications for observability. We found that, in the mass range of 10^4 -10^5 M_{sun}, our DSs are hotter by a factor of 1.5 than those in Freese et al.(2010), are smaller in radius by a factor of 0.6, denser by a factor of 3 - 4, and more luminous by a factor of 2. Our models also confirm previous results, according to which supermassive DSs are very well approximated by (n=3)-polytropes. We also perform a first study of dark star pulsations. Our DS models have pulsation modes with timescales ranging from less than a day to more than two years in their rest frames, at z ~ 15, depending on DM particle mass and overtone number. Such pulsations may someday be used to identify bright, cool objects uniquely as DSs; if properly calibrated, they might, in principle, also supply novel standard candles for cosmological studies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Super-Eddington Stellar Winds Driven by Near-Surface Energy Deposition

We develop analytic and numerical models of the properties of super-Eddington stellar winds, motivated by phases in stellar evolution when super-Eddington energy deposition (via, e.g., unstable fusion, wave heating, or a binary companion) heats a region near the stellar surface. This appears to occur in luminous blue variables (LBVs), Type IIn supernovae progenitors, classical novae, and X-ray bursts. We show that when the wind kinetic power exceeds Eddington, the photons are trapped and behave like a fluid. Convection does not play a significant role in the wind energy transport. The wind properties depend on the ratio of a characteristic speed in the problem vc ~ (Edot G)^{1/5} (where Edot is the heating rate) to the stellar escape speed near the heating region vesc(r_h). For vc > vesc(r_h) the wind kinetic power at large radii Edot_w ~ Edot. For vc < vesc(r_h), most of the energy is used to unbind the wind material and thus Edot_w < Edot. Multidimensional hydrodynamic simulations without radiation diffusion using FLASH and one-dimensional hydrodynamic simulations with radiation diffusion using MESA are in good agreement with the analytic predictions. The photon luminosity from the wind is itself super-Eddington but in many cases the photon luminosity is likely dominated by `internal shocks&#39; in the wind. We discuss the application of our models to eruptive mass loss from massive stars and argue that the wind models described here can account for the broad properties of LBV outflows and the enhanced mass loss in the years prior to Type IIn core-collapse supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The primordial and evolutionary abundance variations in globular-cluster stars: a problem with two unknowns

We demonstrate that among the potential sources of the primordial abundance variations of the proton-capture elements in globular-cluster stars proposed so far, such as the hot-bottom burning in massive AGB stars and H burning in the convective cores of supermassive and fast-rotating massive MS stars, only the supermassive MS stars with M > 10,000 Msun can explain all the observed abundance correlations without any fine-tuning of model parameters. We use our assumed chemical composition for the pristine gas in M13 (NGC6205) and its mixtures with 50% and 90% of the material partially processed in H burning in the 60,000 Msun MS model star as the initial compositions for the normal, intermediate and extreme populations of low-mass stars in this globular cluster, as suggested by its O-Na anti-correlation. We evolve these stars from the zero-age MS to the RGB tip with the thermohaline and parametric prescriptions for the RGB extra mixing. We find that the 3He-driven thermohaline convection cannot explain the evolutionary decline of [C/Fe] in M13 RGB stars, which, on the other hand, is well reproduced with the universal values for the mixing depth and rate calibrated using the observed decrease of [C/Fe] with MV in the globular cluster NGC5466 that does not have the primordial abundance variations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular momentum transport within evolved low-mass stars

Asteroseismology of 1.0-2.0 Msun red giants by the Kepler satellite has enabled the first definitive measurements of interior rotation in both first ascent red giant branch (RGB) stars and those on the Helium burning clump. The inferred rotation rates are 10-30 days for the ~0.2Msun He degenerate cores on the RGB and 30-100 days for the He burning core in a clump star. Using the MESA code we calculate state-of-the-art stellar evolution models of low mass rotating stars from the zero-age main sequence to the cooling white dwarf (WD) stage. We include transport of angular momentum due to rotationally induced instabilities and circulations, as well as magnetic fields in radiative zones (generated by the Tayler-Spruit dynamo). We find that all models fail to predict core rotation as slow as observed on the RGB and during core He burning, implying that an unmodeled angular momentum transport process must be operating on the early RGB of low mass stars. Later evolution of the star from the He burning clump to the cooling WD phase appears to be at nearly constant core angular momentum. We also incorporate the adiabatic pulsation code, ADIPLS, to explicitly highlight this shortfall when applied to a specific Kepler asteroseismic target, KIC8366239. The MESA inlist adopted to calculate the models in this paper can be found at \url{https://authorea.com/1608/} (bottom of the document).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hybrid C-O-Ne White Dwarfs as Progenitors of Diverse SNe Ia

When carbon is ignited off-center in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the center. Whether the C flame will actually be able to reach the center depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the center. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a &#34;hybrid&#34; WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-center C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the center. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The dependence of the evolution of SN type Ia progenitors on the C burning rate uncertainty and parameters of convective boundary mixing

Evolution of a supernova type Ia progenitor requires formation of a CO white dwarf, which implies a dependence on the C burning rate (CBR). It can also be affected by the recently identified possibility of C flame quenching by convective boundary mixing. We present first results of our study of the combined effect of these two potential sources of uncertainty on the SN Ia progenitor evolution. We consider the possibility that the CBR is higher than its currently recommended value by as much as a factor of 1000 if unidentified resonances are important, or that it is significantly lower because of the hindrance effect. For stellar models that assume the Schwarzschild boundary for convection, the maximum initial mass for the formation of CO WDs increases from M_i ~ 5.5 Msun for the CBR factor of 1000 to M_i > 7.0 Msun for the CBR factor of 0.01. For C-flame quenching models, hybrid C-O-Ne WDs form for a range of initial mass of Delta M_i ~ 1 Msun, which increases a fraction of stars that form WDs capable of igniting C in a thermonuclear runaway. The most extreme case is found for the CBR factor of 0.1 that is supported by the hindrance model. This nuclear physics assumption, combined with C flame quenching, leads to the formation of a hybrid C-O-Ne WD with a mass of 1.3 Msun. Such WDs do not need to accrete much mass to reach the Chandrasekhar limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Advanced burning stages and fate of 8-10 Mo stars

The stellar mass range 8<M/Mo<12 corresponds to the most massive AGB stars and the most numerous massive stars. It is host to a variety of supernova progenitors and is therefore very important for galactic chemical evolution and stellar population studies. In this paper, we study the transition from super-AGB star to massive star and find that a propagating neon-oxygen burning shell is common to both the most massive electron capture supernova (EC-SN) progenitors and the lowest mass iron-core collapse supernova (FeCCSN) progenitors. Of the models that ignite neon burning off-center, the 9.5Mo model would evolve to an FeCCSN after the neon-burning shell propagates to the center, as in previous studies. The neon-burning shell in the 8.8Mo model, however, fails to reach the center as the URCA process and an extended (0.6 Mo) region of low Ye (0.48) in the outer part of the core begin to dominate the late evolution; the model evolves to an EC-SN. This is the first study to follow the most massive EC-SN progenitors to collapse, representing an evolutionary path to EC-SN in addition to that from SAGB stars undergoing thermal pulses. We also present models of an 8.75Mo super-AGB star through its entire thermal pulse phase until electron captures on 20Ne begin at its center and of a 12Mo star up to the iron core collapse. We discuss key uncertainties and how the different pathways to collapse affect the pre-supernova structure. Finally, we compare our results to the observed neutron star mass distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

GYRE: A New Open-Source Stellar Oscillation Code

We introduce GYRE, a new open-source stellar oscillation code which solves the adiabatic/non-adiabatic pulsation equations using a novel Magnus Multiple Shooting (MMS) numerical scheme. The code has a global error scaling of up to 6th order in the grid spacing, and can therefore achieve high accuracy with few grid points. It is moreover robust and efficiently makes use of multiple processor cores and/or nodes. We present an example calculation using GYRE, and discuss recent work to integrate GYRE into the asteroseismic optimization module of the MESA stellar evolution code.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hydrogen Burning on Accreting White Dwarfs: Stability, Recurrent Novae, and the Post-Novae Supersoft Source

We examine the properties of white dwarfs (WDs) accreting hydrogen-rich matter in and near the stable burning regime of accretion rates as modeled by time-dependent calculations done with Modules for Experiments in Stellar Astrophysics (MESA). We report the stability boundary for WDs of masses between 0.51 solar masses and 1.34 solar masses as found via time-dependent calculations. We also examine recurrent novae that are accreting at rates close to, but below, the stable burning limit and report their recurrence times and ignition masses. Our dense grid in accretion rates finds the expected minimum possible recurrence times as a function of the WD mass. This enables inferences to be made about the minimum WD mass possible to reach a specific recurrence time. We compare our computational models of post-outburst novae to the stably burning WDs and explicitly calculate the duration and effective temperature (Teff) of the post-novae WD in the supersoft phase. We agree with the measured turnoff time - Teff relation in M31 by Henze and collaborators, infer WD masses in the 1.0-1.3 solar masses range, and predict ejection masses consistent with those observed. We close by commenting on the importance of the hot helium layer generated by stable or unstable hydrogen burning for the short- and long-term evolution of accreting white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star&#39;s ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA&#39;s performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reproducing the observed abundances in RCB and HdC stars with post-double degenerate merger models - constraints on merger and post-merger simulations and physics processes

The R Coronae Borealis (RCB) stars are hydrogen-deficient, variable stars that are most likely the result of He-CO WD mergers. They display extremely low oxygen isotopic ratios, 16O/18O ~ 1 - 10, 12C/13C>=100, and enhancements up to 2.6dex in F and in s-process elements from Zn to La, compared to solar. These abundances provide stringent constraints on the physical processes during and after the double-degenerate merger. As shown before O-isotopic ratios observed in RCB stars cannot result from the dynamic double-degenerate merger phase, and we investigate now the role of the long-term 1D spherical post-merger evolution and nucleosynthesis based on realistic hydrodynamic merger progenitor models. We adopt a model for extra envelope mixing to represent processes driven by rotation originating in the dynamical merger. Comprehensive nucleosynthesis post-processing simulations for these stellar evolution models reproduce, for the first time, the full range of the observed abundances for almost all the elements measured in RCB stars: 16O/18O ratios between 9 and 15, C-isotopic ratios above 100, and ~1.4 - 2.35dex F enhancements, along with enrichments in s-process elements. The nucleosynthesis processes in our models constrain the length and temperature in the dynamic merger shell-of-fire feature as well as the envelope mixing in the post-merger phase. s-process elements originate either in the shell-of-fire merger feature or during the post-merger evolution, but the contribution from the AGB progenitors is negligible. The post-merger envelope mixing must eventually cease ~ 10^6yr after the dynamic merger phase, before the star enters the RCB phase.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The C-flame Quenching by Convective Boundary Mixing in Super-AGB Stars and the Formation of Hybrid C/O/Ne White Dwarfs and SN Progenitors

After off-center C ignition in the cores of super-AGB stars the C flame propagates all the way down to the center, trailing behind it the C-shell convective zone, and thus building a degenerate ONe core. This standard picture is obtained in stellar evolution simulations if the bottom C-shell convection boundary is assumed to be a discontinuity associated with a strict interpretation of the Schwarzschild condition for convective instability. However, this boundary is prone to additional mixing processes, such as thermohaline convection and convective boundary mixing. Using hydrodynamic simulations we show that, contrary to previous results, thermohaline mixing is too inefficient to interfere with the C-flame propagation. However, even a small amount of convective boundary mixing removes the physical conditions required for the C-flame propagation all the way to the center. This result holds even if we allow for some turbulent heat transport in the CBM region. As a result, super AGB stars build in their interiors hybrid C-O-Ne degenerate cores composed of a relatively large CO core (M_CO ~ 0.2 M_sun) surrounded by a thick ONe zone (M_ONe ~ 0.85 M_sun) with another thin CO layer above. If exposed by mass loss, these cores will become hybrid C-O-Ne white dwarfs. Otherwise, the ignition of C-rich material in the central core, surrounded by the thick ONe zone, may trigger a thermonuclear supernova explosion. The quenching of the C-flame may have implications for the ignition mechanism of SN Ia in the double-degenerate merger scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From the CMD of Omega Centauri and (super-)AGB stellar models to a Galactic plane passage gas purging chemical evolution scenario

[Abbreviated] We have investigated the color-magnitude diagram of Omega Centauri and find that the blue main sequence (bMS) can be reproduced only by models that have a of helium abundance in the range Y=0.35-$0.40. To explain the faint subgiant branch of the reddest stars (&#34;MS-a/RG-a&#34; sequence), isochrones for the observed metallicity ([Fe/H]\approx0.7) appear to require both a high age (~13Gyr) and enhanced CNO abundances ([CNO/Fe]\approx0.9$). Y~0.35 must also be assumed in order to counteract the effects of high CNO on turnoff colors, and thereby to obtain a good fit to the relatively blue turnoff of this stellar population. This suggest a short chemical evolution period of time (<1Gyr) for Omega Cen. Our intermediate-mass (super-)AGB models are able to reproduce the high helium abundances, along with [N/Fe]~2 and substantial O depletions if uncertainties in the treatment of convection are fully taken into account. These abundance features distinguish the bMS stars from the dominant [Fe/H] $\approx1.7$ population. The most massive super-AGB stellar models (M_zams>=6.8M_sun, M_He,core>=1.245M_sun) predict too large N-enhancements, which limits their role in contributing to the extreme populations. We show quantitatively that highly He- and N-enriched AGB ejecta have particularly efficient cooling properties. Based on these results and on the reconstruction of the orbit of Omega Cen with respect to the Milky Way we propose the galactic plane passage gas purging scenario for the chemical evolution of this cluster. Our model addresses the formation and properties of the bMS population (including their central location in the cluster). We follow our model descriptively through four passage events, which could explain not only some key properties of the bMS, but also of the MS-a/RGB-a and the s-enriched stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

MESA Models of Classical Nova Outbursts: The Multicycle Evolution and Effects of Convective Boundary Mixing

Novae are cataclysmic variables driven by accretion of H-rich material onto a white-dwarf (WD) star from its low-mass main-sequence binary companion. New time-domain observational capabilities, such as the Palomar Transient Factory and Pan-STARRS, have revealed a diversity of their behaviour that should be theoretically addressed. Nova outbursts depend sensitively on nuclear physics data, and more readily available nova simulations are needed in order to effectively prioritize experimental effort in nuclear astrophysics. In this paper we use the MESA stellar evolution code to construct multicycle nova evolution sequences with CO WD cores. We explore a range of WD masses and accretion rates as well as the effect of different cooling times before the onset of accretion. In addition, we study the dependence on the elemental abundance distribution of accreted material and convective boundary mixing at the core-envelope interface. Models with such convective boundary mixing display an enrichment of the accreted envelope with C and O from the underlying white dwarf that is commensurate with observations. We compare our results with the previous work and investigate a new scenario for novae with the 3He-triggered convection.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Response of Giant Stars To Dynamical-Timescale Mass Loss

We study the response of giant stars to mass loss. One-dimensional simulations of red and asymptotic giant branch stars with mass loss rates from $10^{-3}$ up to a few \msun/yr show in no case any significant radius increase. The largest radius increase of 0.2% was found in the case with the lowest mass loss rate. For dynamical-timescale mass loss rates, that may be encountered during a common envelope phase, the evolution is not adiabatic. The superadiabatic outer layer of the giant&#39;s envelope has a local thermal timescale comparable to the dynamical timescale. Therefore, this layer has enough time to readjust thermally. Moreover, the giant star is driven out of hydrostatic equilibrium and evolves dynamically. In these cases no increase of the stellar radius with respect to its initial value is found. If the mass loss rate is high enough, the superadiabaticity of the outer layer is lost progressively and a radiative zone forms due to a combination of thermal and dynamical readjustment. Conditions for unstable mass transfer based on adiabatic mass loss models that predict a significant radius increase, may need to be re-evaluated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Acoustic Signatures of the Helium Core Flash

All evolved stars with masses M <2 solar masses undergo an initiating off-center helium core flash in their 0.48 solar mass He core as they ascend the red giant branch (RGB). This off-center flash is the first of a few successive helium shell subflashes that remove the core electron degeneracy over 2 Myrs, converting the object into a He burning star. Though characterized by Thomas over 40 years ago, this core flash phase has yet to be observationally probed. Using the Modules for Experiments in Stellar Astrophysics (MESA) code, we show that red giant asteroseismology enabled by space-based photometry (i.e. Kepler and CoRoT) can probe these stars during the flash. The rapid (< 100,000 years) contraction of the red giant envelope after the initiating flash dramatically improves the coupling of the p-modes to the core g-modes, making the detection of l=1 mixed modes possible for these 2 Myrs. This duration implies that 1 in 35 stars near the red clump in the HR diagram will be in their core flash phase. During this time, the star has a g-mode period spacing of Delta P_g = 70-100 seconds, lower than the Delta P_g=250 seconds of He burning stars in the red clump, but higher than the RGB stars at the same luminosity. This places them in an underpopulated part of the large frequency spacing (Delta nu) vs. Delta P_g diagram that should ease their identification amongst the thousands of observed red giants.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modules for Experiments in Stellar Astrophysics (MESA)

Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source libraries for a wide range of applications in computational stellar astrophysics. A newly designed 1-D stellar evolution module, MESA star, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very-low mass to massive stars, including advanced evolutionary phases. MESA star solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. Independently usable modules provide equation of state, opacity, nuclear reaction rates, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own public interface. Examples include comparisons to other codes and show evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets; the complete evolution of a 1 Msun star from the pre-main sequence to a cooling white dwarf; the Solar sound speed profile; the evolution of intermediate mass stars through the thermal pulses on the He-shell burning AGB phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; evolutionary tracks of massive stars from the pre-main sequence to the onset of core collapse; stars undergoing Roche lobe overflow; and accretion onto a neutron star. Instructions for downloading and installing MESA can be found on the project web site (http://mesa.sourceforge.net/).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the alpha formalism for the common envelope interaction

The α-formalism is a common way to parametrize the common envelope interaction between a giant star and a more compact companion. The α parameter describes the fraction of orbital energy released by the companion that is available to eject the giant star&#39;s envelope. By using new, detailed stellar evolutionary calculations we derive a user-friendly prescription for the λ parameter and an improved approximation for the envelope binding energy, thus revising the α equation. We then determine α both from simulations and observations in a self consistent manner. By using our own stellar structure models as well as population considerations to reconstruct the primary&#39;s parameters at the time of the common envelope interaction, we gain a deeper understanding of the uncertainties. We find that systems with very low values of q (the ratio of the companion&#39;s mass to the mass of the primary at the time of the common envelope interaction) have higher values of α. A fit to the data suggests that lower mass companions are left at comparable or larger orbital separations to more massive companions. We conjecture that lower mass companions take longer than a stellar dynamical time to spiral in to the giant&#39;s core, and that this is key to allowing the giant to use its own thermal energy to help unbind its envelope. As a result, although systems with light companions might not have enough orbital energy to unbind the common envelope, they might stimulate a stellar reaction that results in the common envelope ejection.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Evolutionary implications of the new triple-alpha nuclear reaction rate for low mass stars

Context: Ogata et al. (2009; hereafter OKK) presented a theoretical determination of the triple-alpha nuclear reaction rate. Their rate differs from the NACRE rate by many orders of magnitude at temperatures relevant for low mass stars. Aims: We explore the evolutionary implications of adopting the OKK triple-alpha reaction rate in low mass stars and compare the results with those obtained using the NACRE rate. Methods: The triple-alpha reaction rates are compared by following the evolution of stellar models at 1 and 1.5 Msol with Z=0.0002 and Z=0.02. Results: Results show that the OKK rate has severe consequences for the late stages of stellar evolution in low mass stars. Most notable is the shortening--or disappearance--of the red giant phase. Conclusions: The OKK triple-alpha reaction rate is incompatible with observations of extended red giant branches and He burning stars in old stellar systems.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Ellipsoidal Oscillations Induced by Substellar Companions: A Prospect for the Kepler Mission

Hundreds of substellar companions to solar-type stars will be discovered with the Kepler satellite. Kepler&#39;s extreme photometric precision gives access to low-amplitude stellar variability contributed by a variety of physical processes. We discuss in detail the periodic flux modulations arising from the tidal force on the star due to a substellar companion. An analytic expression for the variability is derived in the equilibrium-tide approximation. We demonstrate analytically and through numerical solutions of the linear, nonadiabatic stellar oscillation equations that the equilibrium-tide formula works extremely well for stars of mass <1.4 Msun with thick surface convection zones. More massive stars with largely radiative envelopes do not conform to the equilibrium-tide approximation and can exhibit flux variations $\ga$10 times larger than naive estimates. Over the full range of stellar masses considered, we treat the oscillatory response of the convection zone by adapting a prescription that A. J. Brickhill developed for pulsating white dwarfs. Compared to other sources of periodic variability, the ellipsoidal lightcurve has a distinct dependence on time and system parameters. We suggest that ellipsoidal oscillations induced by giant planets may be detectable from as many as ~100 of the 10^5 Kepler target stars. (Abridged)

0 citations0 reviewsNo code linkedOpen access