Showing 1–50 of 69 results for author: Arnett, D

3D simulations of convective shell Neon-burning in a massive star

Authors: C. Georgy, F. Rizzuti, R. Hirschi, V. Varma, W. D. Arnett, C. Meakin, M. Mocak, A. StJ. Murphy, T. Rauscher

Abstract: The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulen… ▽ More The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same timescale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 17 pages, 20 figures, accepted for publication in MNRAS

3D stellar evolution: hydrodynamic simulations of a complete burning phase in a massive star

Authors: F. Rizzuti, R. Hirschi, W. D. Arnett, C. Georgy, C. Meakin, A. StJ. Murphy, T. Rauscher, V. Varma

Abstract: Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multi-dimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an… ▽ More Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multi-dimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M$_\odot$ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution, thus revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions and compact remnants. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 12 pages, 15 figures. Accepted for publication in MNRAS

Realistic 3D hydrodynamics simulations find significant turbulent entrainment in massive stars

Authors: F. Rizzuti, R. Hirschi, C. Georgy, W. D. Arnett, C. Meakin, A. StJ. Murphy

Abstract: Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multi-dimensional processes over short timescal… ▽ More Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multi-dimensional processes over short timescales, thanks to the recent advances in computing resources. Among these processes, turbulent entrainment leading to mixing across convective boundaries is one of the least understood and most impactful. Here we present the results from a set of hydrodynamics simulations of the neon-burning shell in a massive star, and interpret them in the framework of the turbulent entrainment law from geophysics. Our simulations differ from previous studies in their unprecedented degree of realism in reproducing the stellar environment. Importantly, the strong entrainment found in the simulations highlights the major flaws of the current implementation of convective boundary mixing in 1D stellar models. This study therefore calls for major revisions of how convective boundaries are modelled in 1D, and in particular the implementation of entrainment in these models. This will have important implications for supernova theory, nucleosynthesis, neutron stars and black holes physics. △ Less

Submitted 16 August, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 7 pages, 6 figures, accepted for publication in MNRAS

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Authors: R. Andrassy, J. Higl, H. Mao, M. Mocák, D. G. Vlaykov, W. D. Arnett, I. Baraffe, S. W. Campbell, T. Constantino, P. V. F. Edelmann, T. Goffrey, T. Guillet, F. Herwig, R. Hirschi, L. Horst, G. Leidi, C. Meakin, J. Pratt, F. Rizzuti, F. K. Roepke, P. Woodward

Abstract: Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility o… ▽ More Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study. △ Less

Submitted 26 January, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 20 pages, 19 figures, new Fig. B.2, new clarifications, electronic materials updated, accepted for publication in A&A

Journal ref: A&A 659, A193 (2022)

Convective core entrainment in 1D main sequence stellar models

Authors: L. J. A. Scott, R. Hirschi, C. Georgy, W. D. Arnett, C. Meakin, E. A. Kaiser, S. Ekström, N. Yusof

Abstract: 3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M$_{\odot}$ at solar metallicity ($Z=0.014$) and compare them to previous generations of mod… ▽ More 3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M$_{\odot}$ at solar metallicity ($Z=0.014$) and compare them to previous generations of models and observations on the main sequence. The boundary penetrability, quantified by the bulk Richardson number, $Ri_{\mathrm{B}}$, varies with mass and to a smaller extent with time. The variation of $Ri_{\mathrm{B}}$ with mass is due to the mass dependence of typical convective velocities in the core and hence the luminosity of the star. The chemical gradient above the convective core dominates the variation of $Ri_{\mathrm{B}}$ with time. An entrainment law method can therefore explain the apparent mass dependence of convective boundary mixing through $Ri_{\mathrm{B}}$. New models including entrainment can better reproduce the mass dependence of the main sequence width using entrainment law parameters $A \sim 2 \times 10^{-4}$ and $n=1$. We compare these empirically constrained values to the results of 3D hydrodynamics simulations and discuss implications. △ Less

Submitted 11 March, 2021; v1 submitted 10 March, 2021; originally announced March 2021.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

Relative Importance of Convective Uncertainties in Massive Stars

Authors: Etienne A. Kaiser, Raphael Hirschi, W. David Arnett, Cyril Georgy, Laura J. A. Scott, Andrea Cristini

Abstract: In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called `overshoot', namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary th… ▽ More In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called `overshoot', namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses $15$, $20$ and $25\,\rm{M}_\odot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main-sequence with more CBM, which is more in agreement with observations. Furthermore during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves red-wards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the $^{12}$C to $^{16}$O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduce and we find differences of up to $70\%$ for the core masses and the total mass of the star. △ Less

Submitted 17 July, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: 24 pages, 13 figures, 3 tables, accepted for publication in MNRAS. Stellar tracks and profiles are available at http://www.astro.keele.ac.uk/shyne/datasets/convective-boundary-mixing-uncertainties

Catching Element Formation In The Act

Authors: Chris L. Fryer, Frank Timmes, Aimee L. Hungerford, Aaron Couture, Fred Adams, Wako Aoki, Almudena Arcones, David Arnett, Katie Auchettl, Melina Avila, Carles Badenes, Eddie Baron, Andreas Bauswein, John Beacom, Jeff Blackmon, Stephane Blondin, Peter Bloser, Steve Boggs, Alan Boss, Terri Brandt, Eduardo Bravo, Ed Brown, Peter Brown, Steve Bruenn. Carl Budtz-Jorgensen, Eric Burns , et al. (194 additional authors not shown)

Abstract: Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray… ▽ More Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 14 pages including 3 figures

Report number: LA-UR-18-29748

Dependence of Convective Boundary Mixing on Boundary Properties and Turbulence Strength

Authors: A. Cristini, R. Hirschi, C. Meakin, D. Arnett, C. Georgy, I. Walkington

Abstract: Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with th… ▽ More Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with the PROMPI code. We performed a mean field analysis of the simulations within the Reynolds-averaged Navier-Stokes framework. Both the vertical RMS velocity within the convective region and the bulk Richardson number of the boundaries are found to scale with the driving luminosity as expected from theory. The positions of the convective boundaries were estimated through the composition profiles across them, and the strength of convective boundary mixing was determined by analysing the boundaries within the framework of the entrainment law. We find that the entrainment is approximately inversely proportional to the bulk Richardson number. Although the entrainment law does not encompass all the processes occurring at boundaries, our results support the use of the entrainment law to describe convective boundary mixing in 1D models, at least for the advanced phases. The next steps and challenges ahead are also discussed. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

3D Simulations and MLT: II. Onsager's Ideal Turbulence

Authors: W. David Arnett, Raphael Hirschi, Simon W. Campbell, Miroslav Mocák, Cyril Georgy, Casey Meakin, Andrea Cristini, Laura J. A. Scott, Etienne A. Kaiser, Maxime Viallet

Abstract: We simulate stellar convection at high Reynolds number (Re$\lesssim$7000) with causal time stepping but no explicit viscosity. We use the 3D Euler equations with shock capturing (Colella & Woodward 1984). Anomalous dissipation of turbulent kinetic energy occurs as an emergent feature of advection ("Onsager damping"), caused by the moderate shocks which terminate the turbulent kinetic energy spectr… ▽ More We simulate stellar convection at high Reynolds number (Re$\lesssim$7000) with causal time stepping but no explicit viscosity. We use the 3D Euler equations with shock capturing (Colella & Woodward 1984). Anomalous dissipation of turbulent kinetic energy occurs as an emergent feature of advection ("Onsager damping"), caused by the moderate shocks which terminate the turbulent kinetic energy spectrum; see also (Perry 2021). In strongly stratified stellar convection the asymptotic limit for the global damping length of turbulent kinetic energy is $\ell_d \sim \langle u^3 \rangle /\langle ε\rangle$. This "dissipative anomaly" (Onsager 1949) fixes the value of the "mixing length parameter", $α= \ell_{\rm MLT}/H_P =\overline{\langleΓ_1\rangle}$, which is $\sim\, 5/3$ for complete ionization. The estimate is numerically robust, agrees to within 10% with estimates from stellar evolution with constant $α$. For weak stratification $\ell_d$ shrinks to the depth of a thin convective region. Our flows are filamentary, produce surfaces of separation at boundary layers, resolve the energy-containing eddies, and develop a turbulent cascade down to the grid scale which agrees with the $4096^3$ direct numerical simulation of Kaneda (2003). The cascade converges quickly, and satisfies a power-law velocity spectrum similar to Kolmogorov (1941). Our flows exhibit intermittency, anisotropy, and interactions between coherent structures, features missing from K41 theory. We derive a dissipation rate from Reynolds stresses which agrees with (i) our flows, (ii) experiment (Warhaft 2002), and (iii) high Re simulations of the Navier-Stokes equations (Iyer, et al. 2018). △ Less

Submitted 23 April, 2021; v1 submitted 10 October, 2018; originally announced October 2018.

Comments: 31 pages, 6 figures. The time evolution and the fly-through movies may be found at \url{http://www.astro.keele.ac.uk/shyne/321D/convection-and-convective-boundary-mixing/visualisations}

3D Simulations and MLT: I. Renzini's Critique

Authors: W. David Arnett, Casey Meakin, Raphael Hirschi, Andrea Cristini, Cyril Georgy, Simon Campbell, Laura Scott, Etienne Kaiser

Abstract: Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES)… ▽ More Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES) automatically connect large scale turbulence to a Kolmogorov cascade below the grid scale, allowing turbulent boundary layers to remove singularities that appear in the theory. Interactions between coherent structures give multi-modal behavior, driving intermittency and fluctuations. Reynolds averaging (RA) allows us to abstract the essential features of this dynamical behavior of boundaries which are appropriate to stellar evolution, and consider how they relate static boundary conditions (Richardson, Schwarzschild or Ledoux). We clarify several questions concerning when and why MLT works, and does not work, using both analytical theory and 3D high resolution numerical simulations. The composition gradients and boundary layer structure which are produced by our simulations suggest a self-consistent approach to boundary layers, removing the need for ad hoc procedures for 'convective overshooting' and `semi-convection'. In a companion paper we quantify the adequacy of our numerical resolution, determine of the length scale of dissipation (the `mixing length') without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and extend MLT to deal with strong stratification. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Relative Importance of Convective Uncertainties

Authors: Etienne A. Kaiser, Raphael Hirschi, W. David Arnett, Andrea Cristini, Cyril Georgy, Laura J. A. Scott

Abstract: Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of… ▽ More Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of extra mixing, and their impact on the early evolution of massive stars. With increasing convective boundary mixing (CBM), we find a convergence of the two different boundary locations, a decreasing blue to red super giant ratio and a reduced importance of semiconvection. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 3 pages, 1 figure; to appear in Springer Proceedings in Physics (Proc. of Itl. Conf. "Nuclei in the Cosmos XV", LNGS Assergi, Italy, June 2018)

Turbulent Mixing and Nuclear Burning in Stellar Interiors

Authors: Miroslav Mocák, Casey Meakin, Simon Wattana Campbell, David Arnett

Abstract: The turbulent burning of nuclei is a common phenomenon in the evolution of stars. Here we examine a challenging case: the merging of the neon and oxygen burning shells in a 23 M$_{\odot}$ star. A previously unknown quasi-steady state is established by the interplay between mixing, turbulent transport, and nuclear burning. The resulting stellar structure has two burning shells within a single conve… ▽ More The turbulent burning of nuclei is a common phenomenon in the evolution of stars. Here we examine a challenging case: the merging of the neon and oxygen burning shells in a 23 M$_{\odot}$ star. A previously unknown quasi-steady state is established by the interplay between mixing, turbulent transport, and nuclear burning. The resulting stellar structure has two burning shells within a single convection zone. We find that the new neon burning layer covers an extended region of the convection zone, with the burning peak occurring substantially below where the Damköhler number first becomes equal to unity. These characteristics differ from those predicted by 1D stellar evolution models of similar ingestion events. We develop the mean-field turbulence equations that govern compositional evolution, and use them to interpret our data set. An important byproduct is a means to quantify sub-grid-scale effects intrinsic to the numerical hydrodynamic scheme. For implicit large eddy simulations, the analysis method is particularly powerful because it can reveal where and how simulated flows are modified by resolution, and provide straightforward physical interpretations of the effects of dissipation or induced transport. Focusing on the mean-field composition variance equations for our analysis, we recover a Kolmogorov rate of turbulent dissipation without it being imposed, in agreement with previous results which used the turbulent kinetic energy equation. △ Less

Submitted 20 August, 2018; originally announced August 2018.

Comments: 17 pages, 8 figures, Submitted to MNRAS

Progenitors of Core-Collapse Supernovae

Authors: Raphael Hirschi, David Arnett, Andrea Cristini, Cyril Georgy, Casey Meakin, Ian Walkington

Abstract: Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss som… ▽ More Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss some of the major uncertainties in the modelling of massive stars, with a particular emphasis on the treatment of convection in 1D stellar evolution codes. Finally, we present new 3D hydrodynamic simulations of convection in carbon burning and list key points to take from 3D hydrodynamic studies for the development of new prescriptions for convective boundary mixing in 1D stellar evolution codes. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 10 pages, 5 figures, "SN 1987A, 30 years later", Proceedings IAU Symposium No. 331, 2017; A. Marcowith, G. Dubner, A. Ray, A. Bykov, & M. Renaud, eds

The First 3D Simulations of Carbon Burning in a Massive Star

Authors: Andrea Cristini, Casey Meakin, Raphael Hirschi, David Arnett, Cyril Georgy, Maxime Viallet

Abstract: We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs… ▽ More We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs via turbulent entrainment as in the case of oxygen burning. The expansion of the boundary into the surrounding stable region and the entrainment rate are smaller at the bottom boundary because it is stiffer than the upper boundary. The results of this and similar studies call for improved CBM prescriptions in 1D stellar evolution models. △ Less

Submitted 19 December, 2017; originally announced December 2017.

Comments: 5 pages, 3 figures. Published in IAUS 329 on 28/07/17

Journal ref: Proceedings of the International Astronomical Union, 12(S329), 237-241 (2017)

A. G. W. Cameron 1925-2005, Biographical Memoir, National Academy of Sciences

Authors: David Arnett

Abstract: Alastair Graham Walker Cameron was an astrophysicist and planetary scientist of broad interests and exceptional originality. A founder of the field of nuclear astrophysics, he developed the theoretical understanding of the chemical elements’ origins and made pioneering connections between the abundances of elements in meteorites to advance the theory that the Moon originated from a giant impact… ▽ More Alastair Graham Walker Cameron was an astrophysicist and planetary scientist of broad interests and exceptional originality. A founder of the field of nuclear astrophysics, he developed the theoretical understanding of the chemical elements’ origins and made pioneering connections between the abundances of elements in meteorites to advance the theory that the Moon originated from a giant impact with the young Earth by an object at least the size of Mars. Cameron was an early and persistent exploiter of computer technology in the theoretical study of complex astronomical systems—including nuclear reactions in supernovae, the structure of neutron stars, and planetary collisions. △ Less

Submitted 17 August, 2017; originally announced August 2017.

Comments: 12 pages, I picture

The Sirius System and its Astrophysical Puzzles: Hubble Space Telescope and Ground-Based Astrometry

Authors: Howard E. Bond, Gail H. Schaefer, Ronald L. Gilliland, Jay B. Holberg, Brian D. Mason, Irving W. Lindenblad, Miranda Seitz-McLeese, W. David Arnett, Pierre Demarque, Federico Spada, Patrick A. Young, Martin A. Barstow, Matthew R. Burleigh, Donald Gudehus

Abstract: Sirius, the seventh-nearest stellar system, is a visual binary containing the metallic-line A1 V star Sirius A, brightest star in the sky, orbited in a 50.13-year period by Sirius B, the brightest and nearest white dwarf (WD). Using images obtained over nearly two decades with the Hubble Space Telescope (HST), along with photographic observations covering almost 20 years, and nearly 2300 historica… ▽ More Sirius, the seventh-nearest stellar system, is a visual binary containing the metallic-line A1 V star Sirius A, brightest star in the sky, orbited in a 50.13-year period by Sirius B, the brightest and nearest white dwarf (WD). Using images obtained over nearly two decades with the Hubble Space Telescope (HST), along with photographic observations covering almost 20 years, and nearly 2300 historical measurements dating back to the 19th century, we determine precise orbital elements for the visual binary. Combined with the parallax and the motion of the A component, these elements yield dynamical masses of 2.063+/-0.023 Msun and 1.018+/-0.011 Msun for Sirius A and B, respectively. Our precise HST astrometry rules out third bodies orbiting either star in the system, down to masses of ~15-25 Mjup. The location of Sirius B in the H-R diagram is in excellent agreement with theoretical cooling tracks for WDs of its dynamical mass, and implies a cooling age of ~126 Myr. The position of Sirius B in the mass-radius plane is also consistent with WD theory, assuming a carbon-oxygen core. Including the pre-WD evolutionary timescale of the assumed progenitor, the total age of Sirius B is about 228+/-10 Myr. We calculated evolutionary tracks for stars with the dynamical mass of Sirius A, using two independent codes. We find it necessary to assume a slightly sub-solar metallicity, of about 0.85 Zsun, to fit its location in the luminosity-radius plane. The age of Sirius A based on these models is about 237-247 Myr, with uncertainties of +/-15 Myr, consistent with that of the WD companion. We discuss astrophysical puzzles presented by the Sirius system, including the probability that the two stars must have interacted in the past, even though there is no direct evidence for this, and the orbital eccentricity remains high. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Comments: Accepted by Astrophysical Journal

Synergies between Asteroseismology and Three-dimensional Simulations of Stellar Turbulence

Authors: W. David Arnett, E. Moravveji

Abstract: Turbulent mixing of chemical elements by convection has fundamental effects on the evolution of stars. The standard algorithm at present, mixing-length theory (MLT), is intrinsically local, and must be supplemented by extensions with adjustable parameters. As a step toward reducing this arbitrariness, we compare asteroseismically inferred internal structures of two Kepler slowly pulsating B stars… ▽ More Turbulent mixing of chemical elements by convection has fundamental effects on the evolution of stars. The standard algorithm at present, mixing-length theory (MLT), is intrinsically local, and must be supplemented by extensions with adjustable parameters. As a step toward reducing this arbitrariness, we compare asteroseismically inferred internal structures of two Kepler slowly pulsating B stars (SPB's; $M\sim 3.25 M_\odot$) to predictions of 321D turbulence theory, based upon well-resolved, truly turbulent three-dimensional simulations (Arnett , et al. 2015, Christini, et al. 2016) which include boundary physics absent from MLT. We find promising agreement between the steepness and shapes of the theoretically-predicted composition profile outside the convective region in 3D simulations and in asteroseismically constrained composition profiles in the best 1D models of the two SPBs. The structure and motion of the boundary layer, and the generation of waves, are discussed. △ Less

Submitted 26 January, 2017; originally announced January 2017.

Comments: 5 pages, 2 figures, accepted by ApJ Letters 1/26/2016

Asteroseismic Constraints on the Models of Hot B Subdwarfs: Convective Helium-Burning Cores

Authors: J. -T. Schindler, Elizabeth M. Green, W. David Arnett

Abstract: Asteroseismology of non-radial pulsations in Hot B Subdwarfs (sdB stars) offers a unique view into the interior of core-helium-burning stars. Ground-based and space-borne high precision light curves allow for the analysis of pressure and gravity mode pulsations to probe the structure of sdB stars deep into the convective core. As such asteroseismological analysis provides an excellent opportunity… ▽ More Asteroseismology of non-radial pulsations in Hot B Subdwarfs (sdB stars) offers a unique view into the interior of core-helium-burning stars. Ground-based and space-borne high precision light curves allow for the analysis of pressure and gravity mode pulsations to probe the structure of sdB stars deep into the convective core. As such asteroseismological analysis provides an excellent opportunity to test our understanding of stellar evolution. In light of the newest constraints from asteroseismology of sdB and red clump stars, standard approaches of convective mixing in 1D stellar evolution models are called into question. The problem lies in the current treatment of overshooting and the entrainment at the convective boundary. Unfortunately no consistent algorithm of convective mixing exists to solve the problem, introducing uncertainties to the estimates of stellar ages. Three dimensional simulations of stellar convection show the natural development of an overshooting region and a boundary layer. In search for a consistent prescription of convection in one dimensional stellar evolution models, guidance from three dimensional simulations and asteroseismological results is indispensable. △ Less

Submitted 7 December, 2016; originally announced December 2016.

Comments: 10 pages, 3 figures, proceedings of the Joint TASC2 - KASC9 Workshop - SPACEINN - HELAS8 Conference "Seismology of the Sun and the Distant Stars 2016". To be published by the EPJ Web of Conferences

Pre-nebular light curves of type I supernovae

Authors: W. David Arnett, Christopher L. Fryer, Thomas Matheson

Abstract: We compare analytic predictions of supernova light curves with recent high quality data from SN2011fe (Ia), from KSN2011b (Ia), and the Palomar Transient Factory (PTF) and the La Silla-QUEST variability survey (LSQ) (Ia). Because of the steady, fast cadence of observations, KSN2011b provides unique new information on SNe Ia: the smoothness of the light curve, which is consistent with significant… ▽ More We compare analytic predictions of supernova light curves with recent high quality data from SN2011fe (Ia), from KSN2011b (Ia), and the Palomar Transient Factory (PTF) and the La Silla-QUEST variability survey (LSQ) (Ia). Because of the steady, fast cadence of observations, KSN2011b provides unique new information on SNe Ia: the smoothness of the light curve, which is consistent with significant large-scale mixing during the explosion, possibly due to 3D effects (e.g., Rayleigh-Taylor instabilities), and provides support for a slowly-varying leakage (mean opacity). For a more complex light curve (SN2008D, SNIb), we separate the luminosity due to multiple causes and indicate the possibility of a radioactive plume. The early rise in luminosity is shown to be affected by the opacity (leakage rate) for thermal and non-thermal radiation. A general derivation of Arnett's rule again shows that it depends upon {\em all} processes heating the plasma, not just radioactive ones, so that SNe Ia will differ from SNe Ibc if the latter have multiple heating processes. △ Less

Submitted 13 August, 2017; v1 submitted 26 November, 2016; originally announced November 2016.

Comments: 18 pages, 7 figures; accepted ApJ 2017; expanded comparative discussion of types Ia, b, and c; succinct presentation of light curve algorithm to help users catch implementation errors

3D Hydrodynamic Simulations of Carbon Burning in Massive Stars

Authors: Andrea Cristini, Casey Meakin, Raphael Hirschi, David Arnett, Cyril Georgy, Maxime Viallet

Abstract: We present the first detailed three-dimensional (3D) hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon burning shell within a 15$\,\textrm{M}_\odot$ one-dimensional stellar evolution model. We consider models with $128^3$, $256^3$, $512^3$ and $1024^3$ zones. The turbulent flow proper… ▽ More We present the first detailed three-dimensional (3D) hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon burning shell within a 15$\,\textrm{M}_\odot$ one-dimensional stellar evolution model. We consider models with $128^3$, $256^3$, $512^3$ and $1024^3$ zones. The turbulent flow properties of these carbon burning simulations are very similar to the oxygen burning case. We performed a mean field analysis of the kinetic energy budgets within the Reynolds-averaged Navier-Stokes framework. For the upper convective boundary region, we find that the numerical dissipation is insensitive to resolution for linear mesh resolutions above 512 grid points. For the stiffer, more stratified lower boundary, our highest resolution model still shows signs of decreasing sub-grid dissipation suggesting it is not yet numerically converged. We find that the widths of the upper and lower boundaries are roughly 30% and 10% of the local pressure scale heights, respectively. The shape of the boundaries is significantly different from those used in stellar evolution models. As in past oxygen-shell burning simulations, we observe entrainment at both boundaries in our carbon-shell burning simulations. In the large Péclet number regime found in the advanced phases, the entrainment rate is roughly inversely proportional to the bulk Richardson number, Ri$_{\rm B}$ ($\propto $Ri${\rm_B}^{-α}$, $0.5\lesssim α\lesssim 1.0$). We thus suggest the use of Ri$_{\rm B}$ as a means to take into account the results of 3D hydrodynamics simulations in new 1D prescriptions of convective boundary mixing. △ Less

Submitted 18 June, 2017; v1 submitted 17 October, 2016; originally announced October 2016.

Comments: 26 pages, 15 figures, accepted for publication in MNRAS, movie available at the following URL: http://www.astro.keele.ac.uk/shyne/321D/convection-and-convective- boundary-mixing/visualisations/very-high-resolution-movie-of-the-c-shell/view

Key Issues Review: Numerical studies of turbulence in stars

Authors: W. David Arnett, Casey Meakin

Abstract: The numerical simulation of turbulence in stars has led to a rich set of possibilities regarding stellar pulsations, asteroseismology, thermonuclear yields, and formation of neutron stars and black holes. The breaking of symmetry by turbulent flow grows in amplitude as collapse is approached, which insures that the conditions at the onset of collapse are not spherical. This lack of spherical symme… ▽ More The numerical simulation of turbulence in stars has led to a rich set of possibilities regarding stellar pulsations, asteroseismology, thermonuclear yields, and formation of neutron stars and black holes. The breaking of symmetry by turbulent flow grows in amplitude as collapse is approached, which insures that the conditions at the onset of collapse are not spherical. This lack of spherical symmetry has important implications for the mechanism of explosion and ejected nucleosynthesis products. Numerical resolution of several different types of three--dimensional (3D) stellar simulations are compared; it is suggested that core collapse simulations may be under-resolved. New physical effects which appear in 3D are summarized. Connections between simulations of progenitor explosion and observations of supernova remnants (SNR) are discussed. Present treatment of boundaries, for mixing regions during He--burning, requires revision. △ Less

Submitted 17 March, 2016; originally announced March 2016.

Comments: 8 pages, 1 figure, 1 table, submitted to Reports on Progress in Physics

Convective Properties of Rotating Two-Dimensional Core-Collapse Supernova Progenitors

Authors: Emmanouil Chatzopoulos, Sean M. Couch, W. David Arnett, F. X. Timmes

Abstract: We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow… ▽ More We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in one-dimensional (1D) stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions. △ Less

Submitted 25 January, 2016; v1 submitted 21 January, 2016; originally announced January 2016.

Comments: 10 pages, 9 figures

Linking 1D Evolutionary to 3D Hydrodynamical Simulations of Massive Stars

Authors: Andréa Cristini, Casey Meakin, Raphael Hirschi, David Arnett, Cyril Georgy, Maxime Viallet

Abstract: Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective… ▽ More Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective boundaries in stellar models has been studied for decades but still lacks a definitive solution. This paper presents initial results of a study on convective boundary mixing (CBM) in massive stars. The 'stiffness' of a convective boundary can be quantified using the bulk Richardson number ($\textrm{Ri}_B$), the ratio of the potential energy for restoration of the boundary to the kinetic energy of turbulent eddies. A 'stiff' boundary ($\textrm{Ri}_B \sim 10^4$) will suppress CBM, whereas in the opposite case a 'soft' boundary ($\textrm{Ri}_B \sim 10$) will be more susceptible to CBM. One of the key results obtained so far is that lower convective boundaries (closer to the centre) of nuclear burning shells are 'stiffer' than the corresponding upper boundaries, implying limited CBM at lower shell boundaries. This is in agreement with 3D hydrodynamic simulations carried out by Meakin and Arnett [The Astrophysical Journal 667:448-475, 2007]. This result also has implications for new CBM prescriptions in massive stars as well as for nuclear burning flame front propagation in Super-Asymptotic Giant Branch stars and also the onset of novae. △ Less

Submitted 7 January, 2016; originally announced January 2016.

Comments: Accepted for publication (12/12/15) in the Physica Scripta focus issue on Turbulent Mixing and Beyond

Towards 21st Century Stellar Models: Star Clusters, Supercomputing, and Asteroseismology

Authors: S. W. Campbell, T. N. Constantino, V. D'Orazi, C. Meakin, D. Stello, J. Christensen-Dalsgaard, C. Kuehn, G. M. De Silva, W. D. Arnett, J. C. Lattanzio, B. T. MacLean

Abstract: Stellar models provide a vital basis for many aspects of astronomy and astrophysics. Recent advances in observational astronomy -- through asteroseismology, precision photometry, high-resolution spectroscopy, and large-scale surveys -- are placing stellar models under greater quantitative scrutiny than ever. The model limitations are being exposed and the next generation of stellar models is neede… ▽ More Stellar models provide a vital basis for many aspects of astronomy and astrophysics. Recent advances in observational astronomy -- through asteroseismology, precision photometry, high-resolution spectroscopy, and large-scale surveys -- are placing stellar models under greater quantitative scrutiny than ever. The model limitations are being exposed and the next generation of stellar models is needed as soon as possible. The current uncertainties in the models propagate to the later phases of stellar evolution, hindering our understanding of stellar populations and chemical evolution. Here we give a brief overview of the evolution, importance, and substantial uncertainties of core helium burning stars in particular and then briefly discuss a range of methods, both theoretical and observational, that we are using to advance the modelling. △ Less

Submitted 15 December, 2015; originally announced December 2015.

Comments: To appear in Astronomische Nachrichten, special issue "Reconstruction the Milky Way's History: Spectroscopic surveys, Asteroseismology and Chemo-dynamical models", Guest Editors C. Chiappini, J. Montalbán, and M. Steffen, AN 2016 (in press)

Toward a consistent use of overshooting parametrizations in 1D stellar evolution codes

Authors: Maxime Viallet, Casey Meakin, Vincent Prat, David Arnett

Abstract: Several parametrizations for overshooting in 1D stellar evolution calculations coexist in the literature. These parametrizations are used somewhat arbitrarily in stellar evolution codes, based on what works best for a given problem, or even for historical reasons related to the development of each code. We bring attention to the fact that these different parametrizations correspond to different ph… ▽ More Several parametrizations for overshooting in 1D stellar evolution calculations coexist in the literature. These parametrizations are used somewhat arbitrarily in stellar evolution codes, based on what works best for a given problem, or even for historical reasons related to the development of each code. We bring attention to the fact that these different parametrizations correspond to different physical regimes of overshooting, depending whether the effects of radiation are dominant, marginal, or negligible. Our analysis is based on previously published theoretical results, as well as multidimensional hydrodynamical simulations of stellar convection where the interaction between the convective region and a stably-stratified region is observed. Although the underlying hydrodynamical processes are the same, the outcome of the overshooting process is profoundly affected by radiative effects. Using a simple picture of the scales involved in the overshooting process, we show how three regimes are obtained, depending on the importance of radiative effects. These three regimes correspond to the different behaviors observed in hydrodynamical simulations so far, and to the three types of parametrizations used in 1D codes. We suggest that the existing parametrizations for overshooting should coexist in 1D stellar evolution codes, and should be applied consistently at convective boundaries depending on the local physical conditions. △ Less

Submitted 26 June, 2015; v1 submitted 9 June, 2015; originally announced June 2015.

Comments: 5 pages, 2 figures, to appear in A&A as a regular paper. Last version: language editing + typos in Eq. (6) & (9) corrected

Journal ref: A&A 580, A61 (2015)

The Three Dimensional Evolution to Core Collapse of a Massive Star

Authors: Sean M. Couch, Emmanouil Chatzopoulos, W. David Arnett, F. X. Timmes

Abstract: We present the first three dimensional (3D) simulation of the final minutes of iron core growth in a massive star, up to and including the point of core gravitational instability and collapse. We self-consistently capture the development of strong convection driven by violent Si burning in the shell surrounding the iron core. This convective burning builds the iron core to its critical (Chandrasek… ▽ More We present the first three dimensional (3D) simulation of the final minutes of iron core growth in a massive star, up to and including the point of core gravitational instability and collapse. We self-consistently capture the development of strong convection driven by violent Si burning in the shell surrounding the iron core. This convective burning builds the iron core to its critical (Chandrasekhar) mass and collapse ensues, driven by electron capture and photodisintegration. The non-spherical structure and motion (turbulent fluctuations) generated by 3D convection is substantial at the point of collapse. We examine the impact of such physically-realistic 3D initial conditions on the core-collapse supernova mechanism using 3D simulations including multispecies neutrino leakage. We conclude that non-spherical progenitor structure should not be ignored, and has a significant and favorable impact on the likelihood for neutrino-driven explosions. △ Less

Submitted 22 June, 2015; v1 submitted 7 March, 2015; originally announced March 2015.

Comments: 7 pages, 5 figures, accepted for publication in ApJ Letters. Movies may be viewed at http://flash.uchicago.edu/~smc/progen3d

Beyond Mixing-length Theory: a step toward 321D

Authors: W. David Arnett, Casey Meakin, Maxime Viallet, Simon W. Campbell, John Lattanzio, Miroslav Mocák

Abstract: We examine the physical basis for algorithms to replace mixing-length theory (MLT) in stellar evolutionary computations. Our 321D procedure is based on numerical solutions of the Navier-Stokes equations. These implicit large eddy simulations (ILES) are three-dimensional (3D), time-dependent, and turbulent, including the Kolmogorov cascade. We use the Reynolds-averaged Navier-Stokes (RANS) formulat… ▽ More We examine the physical basis for algorithms to replace mixing-length theory (MLT) in stellar evolutionary computations. Our 321D procedure is based on numerical solutions of the Navier-Stokes equations. These implicit large eddy simulations (ILES) are three-dimensional (3D), time-dependent, and turbulent, including the Kolmogorov cascade. We use the Reynolds-averaged Navier-Stokes (RANS) formulation to make concise the 3D simulation data, and use the 3D simulations to give closure for the RANS equations. We further analyze this data set with a simple analytical model, which is non-local and time-dependent, and which contains both MLT and the Lorenz convective roll as particular subsets of solutions. A characteristic length (the damping length) again emerges in the simulations; it is determined by an observed balance between (1) the large-scale driving, and (2) small-scale damping. The nature of mixing and convective boundaries is analyzed, including dynamic, thermal and compositional effects, and compared to a simple model. We find that (1) braking regions (boundary layers in which mixing occurs) automatically appear {\it beyond} the edges of convection as defined by the Schwarzschild criterion, (2) dynamic (non-local) terms imply a non-zero turbulent kinetic energy flux (unlike MLT), (3) the effects of composition gradients on flow can be comparable to thermal effects, and (4) convective boundaries in neutrino-cooled stages differ in nature from those in photon-cooled stages (different Péclet numbers). The algorithms are based upon ILES solutions to the Navier-Stokes equations, so that, unlike MLT, they do not require any calibration to astronomical systems in order to predict stellar properties. Implications for solar abundances, helioseismology, asteroseismology, nucleosynthesis yields, supernova progenitors and core collapse are indicated. △ Less

Submitted 9 June, 2015; v1 submitted 1 March, 2015; originally announced March 2015.

Comments: 22 pages, 4 figures, 2 tables; significantly re-written, critique of Pasetto, et al. model added, accepted for publication by ApJ

Exploring stellar evolution models of sdB stars using MESA

Authors: Jan-Torge Schindler, Elizabeth M. Green, W. David Arnett

Abstract: Stellar evolution calculations have had great success reproducing the observed atmospheric properties of different classes of stars. Recent detections of g-mode pulsations in evolved He burning stars allow a rare comparison of their internal structure with stellar models. Asteroseismology of subdwarf B stars suggests convective cores of $0.22-0.28\,M_\odot$, $\gtrsim 45\,\%$ of the total stellar m… ▽ More Stellar evolution calculations have had great success reproducing the observed atmospheric properties of different classes of stars. Recent detections of g-mode pulsations in evolved He burning stars allow a rare comparison of their internal structure with stellar models. Asteroseismology of subdwarf B stars suggests convective cores of $0.22-0.28\,M_\odot$, $\gtrsim 45\,\%$ of the total stellar mass. Previous studies found significantly smaller convective core masses ($\lesssim 0.19\,M_\odot$) at a comparable evolutionary stage. We evolved stellar models with MESA (Modules for Experiments in Stellar Astrophysics) to explore how well the interior structure inferred from asteroseismology can be reproduced by standard algorithms. Our qualitative evolutionary paths, position in the $\log g-T_{\rm{eff}}$ diagram and model timescales are consistent with previous results. SdB masses from our full evolutionary sequences fall within the range of the empirical sdB mass distribution, but are nearly always lower than the median. Using standard MLT with atomic diffusion we find convective core masses of $\sim 0.17-0.18 M_\odot$, averaged over the entire sdB lifetime. We can increase the convective core sizes to be as large as those inferred from asteroseismology, but only for extreme values of the overshoot parameter (overshoot gives numerically unstable and physically unrealistic behavior at the boundary). High resolution three-dimensional (3D) simulations of turbulent convection in stars suggest that the Schwarzschild criterion for convective mixing sytematically {\it underestimates the actual extent of mixing because a boundary layer forms.} Accounting for this would decrease the errors in both sdB total and convective core masses. △ Less

Submitted 26 May, 2015; v1 submitted 29 October, 2014; originally announced October 2014.

Comments: This paper has been updated to a standard case with OPALII opacities. 13 pages, 6 figures. Submitted to ApJ, Accepted May 2015

Linking 1D Stellar Evolution to 3D Hydrodynamical Simulations

Authors: Andrea Cristini, Raphael Hirschi, Cyril Georgy, Casey Meakin, David Arnett, Maxime Viallet

Abstract: In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness'. Models for a 15Mo star were computed. We f… ▽ More In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness'. Models for a 15Mo star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are `stiffer' according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae. △ Less

Submitted 28 October, 2014; originally announced October 2014.

Comments: 2 pages, 1 figure. To appear in proceedings of the IAU Symposium 307: New Windows on Massive Stars: Asteroseismology, Interferometry and Spectropolarimetry

3D and Some Other Things Missing from the Theory of Massive Star Evolution

Authors: W. David Arnett

Abstract: This is a sketch of a 321D approximation for stellar convection which is nonlocal, and thus has nonzero fluxes of KE (to be published in more detail elsewhere). Boundary conditions are discussed in a fluid dynamics context (i.e., predictions for overshoot, semiconvection and entrainment are analyzed). We plan to add this as an option to MESA. Inclusion of KE fluxes seems to help resolve the solar… ▽ More This is a sketch of a 321D approximation for stellar convection which is nonlocal, and thus has nonzero fluxes of KE (to be published in more detail elsewhere). Boundary conditions are discussed in a fluid dynamics context (i.e., predictions for overshoot, semiconvection and entrainment are analyzed). We plan to add this as an option to MESA. Inclusion of KE fluxes seems to help resolve the solar abundance problem (Asplund 2009). Smaller cores may ease the explosion problems with core collapse supernova simulations. △ Less

Submitted 1 August, 2014; originally announced August 2014.

Comments: 11 pages, 2 figures, IAU S307

Compressible Hydrodynamic Mean-Field Equations in Spherical Geometry and their Application to Turbulent Stellar Convection Data

Authors: Miroslav Mocák, Casey Meakin, Maxime Viallet, David Arnett

Abstract: We present a statistical analysis of turbulent convection in stars within our Reynolds-Averaged Navier Stokes (RANS) framework in spherical geometry which we derived from first principles. The primary results reported in this document include: (1) an extensive set of mean-field equations for compressible, multi-species hydrodynamics, and (2) corresponding mean-field data computed from various simu… ▽ More We present a statistical analysis of turbulent convection in stars within our Reynolds-Averaged Navier Stokes (RANS) framework in spherical geometry which we derived from first principles. The primary results reported in this document include: (1) an extensive set of mean-field equations for compressible, multi-species hydrodynamics, and (2) corresponding mean-field data computed from various simulation models. Some supplementary scale analysis data is also presented. The simulation data which is presented includes: (1) shell convection during oxygen burning in a 23 solar mass supernova progenitor, (2) envelope convection in a 5 solar mass red giant, (3) shell convection during the helium flash, and (4) a hydrogen injection flash in a 1.25 solar mass star. These simulations have been partially described previously in Meakin [2006], Meakin and Arnett [2007a,b, 2010], Arnett et al. [2009, 2010], Viallet et al. [2011, 2013a,b] and Mocak et al. [2009, 2011]. New data is also included in this document with several new domain and resolution configurations as well as some variations in the physical model such as convection zone depth and driving source term. The long term goal of this work is to aid in the development of more sophisticated models for treating hydrodynamic phenomena (e.g., turbulent convection) in the field of stellar evolution by providing a direct link between 3D simulation data and the mean fields which are modeled by 1D stellar evolution codes. As such, this data can be used to test previously proposed turbulence models found in the literature and sometimes used in stellar modeling. This data can also serve to test basic physical principles for model building and inspire new prescriptions for use in 1D evolution codes. △ Less

Submitted 21 January, 2014; originally announced January 2014.

Comments: 176 Pages, 115 Figures

Chaos and Turbulent Nucleosynthesis Prior to a Supernova Explosion

Authors: W. David Arnett, Casey Meakin, Maxime Viallet

Abstract: Three-dimensional (3D), time dependent numerical simulations, of flow of matter in stars, now have sufficient resolution to be fully turbulent. The late stages of the evolution of massive stars, leading up to core collapse to a neutron star (or black hole), and often to supernova explosion and nucleosynthesis, are strongly convective because of vigorous neutrino cooling and nuclear heating. Unlike… ▽ More Three-dimensional (3D), time dependent numerical simulations, of flow of matter in stars, now have sufficient resolution to be fully turbulent. The late stages of the evolution of massive stars, leading up to core collapse to a neutron star (or black hole), and often to supernova explosion and nucleosynthesis, are strongly convective because of vigorous neutrino cooling and nuclear heating. Unlike models based on current stellar evolutionary practice, these simulations show a chaotic dynamics characteristic of highly turbulent flow. Theoretical analysis of this flow, both in the Reynolds-averaged Navier-Stokes (RANS) framework and by simple dynamic models, show an encouraging consistency with the numerical results. It may now be possible to develop physically realistic and robust procedures for convection and mixing which (unlike 3D numerical simulation) may be applied throughout the long life times of stars. In addition, a new picture of the presupernova stages is emerging which is more dynamic and interesting (i.e., predictive of new and newly observed phenomena) than our previous one. △ Less

Submitted 28 January, 2014; v1 submitted 11 December, 2013; originally announced December 2013.

Comments: 11 pages, 2 figures, Submitted to AIP Advances: Stardust, added figures and modest rewriting

Preparing for an Explosion: Hydrodynamic Instabilities and Turbulence in Presupernovae

Authors: Nathan Smith, David Arnett

Abstract: Both observations and direct numerical simulations are discordant with predictions of conventional stellar evolution codes for the latest stages of a massive star's life prior to core collapse. We suggest that the problem lies in the treatment of turbulent convection in these codes, which ignores finite amplitude fluctuations in velocity and temperature, and their nonlinear interaction with nuclea… ▽ More Both observations and direct numerical simulations are discordant with predictions of conventional stellar evolution codes for the latest stages of a massive star's life prior to core collapse. We suggest that the problem lies in the treatment of turbulent convection in these codes, which ignores finite amplitude fluctuations in velocity and temperature, and their nonlinear interaction with nuclear burning. The hydrodynamic instabilities that may arise prompt us to discuss a number of far-reaching implications for the fates of massive stars. In particular, we explore connections to enhanced presupernova mass loss, unsteady nuclear burning and consequent eruptions, swelling of the stellar radius that may trigger violent interactions with a companion star, and potential modifications to the core structure that could dramatically impact calculations of the core-collapse mechanism itself. These modifications may be of fundamental importance to the interpretation of measured isotopic anomalies in meteorites, changing the predictions of both mixing and detailed nucleosynthesis, and of young core-collapse supernova remnants. They may also make possible the development of an early warning system for the detection of impending core collapse. △ Less

Submitted 5 February, 2014; v1 submitted 18 July, 2013; originally announced July 2013.

Comments: 13 pages, 1 figure, rewritten in response to referee, accepted ApJ

Journal ref: ApJ, 785:82, 2014

A panchromatic view of the restless SN2009ip reveals the explosive ejection of a massive star envelope

Authors: R. Margutti, D. Milisavljevic, A. M. Soderberg, R. Chornock, B. A. Zauderer, K. Murase, C. Guidorzi, N. E. Sanders, P. Kuin, C. Fransson, E. M. Levesque, P. Chandra, E. Berger, F. B. Bianco, P. J. Brown, P. Challis, E. Chatzopoulos, C. C. Cheung, C. Choi, L. Chomiuk, N. Chugai, C. Contreras, M. R. Drout, R. Fesen, R. J. Foley , et al. (32 additional authors not shown)

Abstract: The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and X… ▽ More The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified. △ Less

Submitted 31 May, 2013; originally announced June 2013.

Comments: Visit https://www.cfa.harvard.edu/~rmargutt/ for a 5-minute presentation highlighting the major results and to hear a sonification of the SN2009ip explosion

The Age and Stellar Parameters of the Procyon Binary System

Authors: James Liebert, Gilles Fontaine, Patrick A. Young, Kurtis A. Williams, David Arnett

Abstract: The Procyon AB binary system (orbital period 40.838 years, a newly-refined determination), is near and bright enough that the component radii, effective temperatures, and luminosities are very well determined, although more than one possible solution to the masses has limited the claimed accuracy. Preliminary mass determinations for each component are available from HST imaging, supported by groun… ▽ More The Procyon AB binary system (orbital period 40.838 years, a newly-refined determination), is near and bright enough that the component radii, effective temperatures, and luminosities are very well determined, although more than one possible solution to the masses has limited the claimed accuracy. Preliminary mass determinations for each component are available from HST imaging, supported by ground-based astrometry and an excellent Hipparcos parallax; we use these for our preferred solution for the binary system. Other values for the masses are also considered. We have employed the TYCHO stellar evolution code to match the radius and luminosity of the F5 IV-V primary star to determine the system's most likely age as 1.87 +/- 0.13 Gyr. Since prior studies of Procyon A found its abundance indistinguishable from solar, the solar composition of Asplund, Grevesse & Sauval (Z=0.014) is assumed for the HR Diagram flitting. An unsuccessful attempt to fit using the older solar abundance scale of Grevesse & Sauval (Z=0.019) is also reported. For Procyon B, eleven new sequences for the cooling of non-DA white dwarfs have been calculated, to investigate the dependence of the cooling age on (1) the mass, (2) the core composition, (3) the helium layer mass, and (4) heavy-element opacities in the helium envelope. Our calculations indicate a cooling age of 1.19+/-0.11 Gyr, which implies that the progenitor mass of Procyon B was 2.59(+0.44,-0.26) Msun. In a plot of initial vs final mass of white dwarfs in astrometric binaries or star clusters (all with age determinations), the Procyon B final mass lies several sigma below a straight line fit. △ Less

Submitted 2 May, 2013; originally announced May 2013.

Comments: Accepted for the Astrophysical Journal, in press

Turbulent convection in stellar interiors. III. Mean-field analysis and stratification effects

Authors: Maxime Viallet, Casey Meakin, David Arnett, Miroslav Mocak

Abstract: We present 3D implicit large eddy simulations (ILES) of the turbulent convection in the envelope of a 5 Msun red giant star and in the oxygen-burning shell of a 23 Msun supernova progenitor. The numerical models are analyzed in the framework of 1D Reynolds-Averaged Navier-Stokes (RANS) equations. The effects of pressure fluctuations are more important in the red giant model, owing to larger strati… ▽ More We present 3D implicit large eddy simulations (ILES) of the turbulent convection in the envelope of a 5 Msun red giant star and in the oxygen-burning shell of a 23 Msun supernova progenitor. The numerical models are analyzed in the framework of 1D Reynolds-Averaged Navier-Stokes (RANS) equations. The effects of pressure fluctuations are more important in the red giant model, owing to larger stratification of the convective zone. We show how this impacts different terms in the mean-field equations. We clarify the driving sources of kinetic energy, and show that the rate of turbulent dissipation is comparable to the convective luminosity. Although our flows have low Mach number and are nearly adiabatic, our analysis is general and can be applied to photospheric convection as well. The robustness of our analysis of turbulent convection is supported by the insensitivity of the mean-field balances to linear mesh resolution. We find robust results for the turbulent convection zone and the stable layers in the oxygen-burning shell model, and robust results everywhere in the red giant model, but the mean fields are not well converged in the narrow boundary regions (which contain steep gradients) in the oxygen-burning shell model. This last result illustrates the importance of unresolved physics at the convective boundary, which governs the mixing there. △ Less

Submitted 13 March, 2013; v1 submitted 27 December, 2012; originally announced December 2012.

Comments: 26 pages, 20 figures, Accepted for publication in ApJ

Hydrodynamic Studies of the Evolution of Recurrent, Symbiotic, and Dwarf Novae: The White Dwarf Components are Growing in Mass

Authors: S. Starrfield, F. X. Timmes, C. Iliadis, W. R. Hix, W. D. Arnett, C. Meakin, W. M. Sparks

Abstract: Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely lum… ▽ More Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely luminous and it seems likely that quiescent nuclear burning is ongoing on the accreting WD. A fundamental question is the secular evolution of the WD. Do the repeated outbursts or quiescent burning in these accreting systems cause the WD to gain or lose mass? If it is gaining mass, can it eventually reach the Chandrasekhar Limit and become a supernova (a SN Ia if it can hide the hydrogen and helium in the system)? In order to better understand these systems, we have begun a new study of the evolution of Thermonuclear Runaways (TNRs) in the accreted envelopes of WDs using a variety of initial WD masses, luminosities and mass accretion rates. We use our 1-D hydro code, NOVA, which includes the new convective algorithm of Arnett, Meakin and Young, the Hix and Thielemann nuclear reaction solver, the Iliadis reaction rate library, the Timmes equation of state, and the OPAL opacities. We assume a solar composition (Lodders abundance distribution) and do not allow any mixing of accreted material with core material. This assumption strongly influences our results. We report here (1) that the WD grows in mass for all simulations so that canonical `steady burning' does not occur, and (2) that only a small fraction of the accreted matter is ejected in some (but not all) simulations. We also find that the accreting systems, before thermonuclear runaway, are too cool to be seen in X-ray searches for SN Ia progenitors. △ Less

Submitted 26 November, 2012; originally announced November 2012.

Comments: 11 pages 10 figures in the 2011 Asiago Symbiotic Workshop

Journal ref: Baltic Astronomy, Vol 21, pages 76 - 87, 2012

The Ring Nebula Around the Blue Supergiant SBW1: Pre-Explosion Snapshot of a SN 1987A Twin

Authors: Nathan Smith, W. David Arnett, John Bally, Adam Ginsburg, Alexei V. Filippenko

Abstract: SBW1 is a B supergiant surrounded by a ring nebula that is a nearby twin of SN 1987A. We present images and spectra of SBW1 obtained with HST, Spitzer, and Gemini South. HST images of SBW1 do not exhibit long Rayleigh-Taylor fingers, which are presumed to cause the hotspots in the SN1987A ring, but instead show a geometrically thin clumpy ring. The radial mass distribution and size scales of inhom… ▽ More SBW1 is a B supergiant surrounded by a ring nebula that is a nearby twin of SN 1987A. We present images and spectra of SBW1 obtained with HST, Spitzer, and Gemini South. HST images of SBW1 do not exhibit long Rayleigh-Taylor fingers, which are presumed to cause the hotspots in the SN1987A ring, but instead show a geometrically thin clumpy ring. The radial mass distribution and size scales of inhomogeneities in SBW1's ring closely resemble those in the SN1987A ring, but the more complete disk expected to reside at the base of the RT fingers is absent. This structure may explain why portions of the SN1987A ring between the hotspots have not yet brightened after more than 15 years. The model we suggest does not require a fast wind colliding with a previous red supergiant wind. More surprisingly, images of SBW1 also reveal diffuse emission filling the interior of the ring in H-alpha and thermal-IR emission; 190K dust dominates the 8-20 micron luminosity (but contains only 1e-5Msun). Cooler (85K dust resides in the equatorial ring itself (dust mass of 5e-3Msun). Diffuse emission extends inward to 1 arcsecond from the central star, where a paucity of emission suggests an inner hole excavated by the wind. We propose that diffuse emission inside the ring arises from an ionized flow of material photoevaporated from the dense ring, and it prevents the supergiant wind from advancing in the equator. This inner emission could correspond to a structure hypothesized to reside around SN1987A that was never directly detected. We suggest that photoionization can play an important dynamical role in shaping the ring nebula, and we speculate that this might help explain the origin of the polar rings around SN1987A. abridged. △ Less

Submitted 24 November, 2012; v1 submitted 19 November, 2012; originally announced November 2012.

Comments: 20 pages, 11 figures, accepted in MNRAS

Theoretical Studies of Accretion of Matter onto White Dwarfs and the Single Degenerate Scenario for Supernovae of Type Ia

Authors: S. Starrfield, C. Iliadis, F. X. Timmes, W. R. Hix, W. D. Arnett, C. Meakin, W. M. Sparks

Abstract: We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and s… ▽ More We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown. △ Less

Submitted 22 October, 2012; originally announced October 2012.

Comments: 25 Pages, Bulletin of the Astronomical Society of India (BASI) in press

Evidence for Type Ia Supernova Diversity from Ultraviolet Observations with the Hubble Space Telescope

Authors: Xiaofeng Wang, Lifan Wang, Alexei V. Filippenko, Eddie Baron, Markus Kromer, Dennis Jack, Tianmeng Zhang, Greg Aldering, Pierre Antilogus, David Arnett, Dietrich Baade, Brian J. Barris, Stefano Benetti, Patrice Bouchet, Adam S. Burrows, Ramon Canal, Enrico Cappellaro, Raymond Carlberg, Elisa di Carlo, Peter Challis, Arlin Crotts, John I. Danziger, Massimo Della Valle, Michael Fink, Ryan J. Foley , et al. (71 additional authors not shown)

Abstract: We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample.… ▽ More We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects. △ Less

Submitted 6 February, 2012; v1 submitted 26 October, 2011; originally announced October 2011.

Comments: 17 pages, 13 figures, accepted by ApJ

Journal ref: Astrophys.J.749:126-142,2012

Towards Realistic Progenitors of Core-Collapse Supernovae

Authors: W. David Arnett, Casey Meakin

Abstract: Two-dimensional (2D) hydrodynamical simulations of progenitor evolution of a 23 solar mass star, close to core collapse (about 1 hour, in 1D), with simultaneously active C, Ne, O, and Si burning shells, are presented and contrasted to existing 1D models (which are forced to be quasi-static). Pronounced asymmetries, and strong dynamical interactions between shells are seen in 2D. Although instigate… ▽ More Two-dimensional (2D) hydrodynamical simulations of progenitor evolution of a 23 solar mass star, close to core collapse (about 1 hour, in 1D), with simultaneously active C, Ne, O, and Si burning shells, are presented and contrasted to existing 1D models (which are forced to be quasi-static). Pronounced asymmetries, and strong dynamical interactions between shells are seen in 2D. Although instigated by turbulence, the dynamic behavior proceeds to sufficiently large amplitudes that it couples to the nuclear burning. Dramatic growth of low order modes is seen, as well as large deviations from spherical symmetry in the burning shells. The vigorous dynamics is more violent than that seen in earlier burning stages in the 3D simulations of a single cell in the oxygen burning shell, or in 2D simulations not including an active Si shell. Linear perturbative analysis does not capture the chaotic behavior of turbulence (e.g., strange attractors such as that discovered by Lorenz), and therefore badly underestimates the vigor of the instability. The limitations of 1D and 2D models are discussed in detail. The 2D models, although flawed geometrically, represent a more realistic treatment of the relevant dynamics than existing 1D models, and present a dramatically different view of the stages of evolution prior to collapse. Implications for interpretation of SN1987A, abundances in young supernova remnants, pre-collapse outbursts, progenitor structure, neutron star kicks, and fallback are outlined. While 2D simulations provide new qualitative insight, fully 3D simulations are needed for a quantitative understanding of this stage of stellar evolution. The necessary properties of such simulations are delineated. △ Less

Submitted 28 January, 2011; originally announced January 2011.

Comments: 26 pages, 1 table, 4 figures

Turbulent Cells in Stars: I. Fluctuations in Kinetic Energy and Luminosity

Authors: W. David Arnett, Casey Meakin

Abstract: Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett, 2007b) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a general instability of the convective cell, requiring only a localized source of heating or cooling. Such fluctuations are shown to be suppressed in simulations of stellar evolution which use mixing-len… ▽ More Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett, 2007b) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a general instability of the convective cell, requiring only a localized source of heating or cooling. Such fluctuations are shown to be suppressed in simulations of stellar evolution which use mixing-length theory (MLT). Quantitatively similar behavior occurs in the model of a convective roll (cell) of Lorenz (1963), which is known to have a strange attractor that gives rise to chaotic fluctuations in time of velocity and, as we show, luminosity. Study of simulations suggests that the behavior of a Lorenz convective roll may resemble that of a cell in convective flow. We examine some implications of this simplest approximation, and suggest paths for improvement. Using the Lorenz model as representative of a convective cell, a multiple-cell model of a convective layer gives total luminosity fluctuations which are suggestive of irregular variables (red giants and supergiants (Schwarzschild 1975)), and of the long secondary period feature in semi-regular AGB variables (Stothers 2010, Wood, Olivier and Kawaler 2004). This "tau-mechanism" is a new source for stellar variability, which is inherently non-linear (unseen in linear stability analysis), and one closely related to intermittency in turbulence. It was already implicit in the 3D global simulations of Woodward, Porter and Jacobs (2003). This fluctuating behavior is seen in extended 2D simulations of CNeOSi burning shells (Arnett and Meakin 2011b), and may cause instability which leads to eruptions in progenitors of core collapse supernovae PRIOR to collapse. △ Less

Submitted 30 July, 2011; v1 submitted 8 December, 2010; originally announced December 2010.

Comments: 30 pages, 13 figures

Time-dependent Turbulence in Stars

Authors: W. David Arnett, Casey Meakin

Abstract: Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett 2007) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a global instability in the convective region, which has been suppressed in calculations of stellar evolution which use mixing-length theory (MLT). Quantitatively similar behavior occurs in the model of a c… ▽ More Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett 2007) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a global instability in the convective region, which has been suppressed in calculations of stellar evolution which use mixing-length theory (MLT). Quantitatively similar behavior occurs in the model of a convective roll (cell) of Lorenz (1963), which is known to have a strange attractor that gives rise to random fluctuations in time.An extension of the Lorenz model, which includes Kolmogorov damping and nuclear burning, is shown to exhibit bursty, recurrent fluctuations like those seen in the 3D simulations. A simple model of a convective layer (composed of multiple Lorenz cells) gives luminosity fluctuations which are suggestive of irregular variables (red giants and supergiants, Schwarzschild 1975). Apparent inconsistencies between Arnett, Meakin, and Young (2009) and Nordlund, Stein, and Asplund (2009) on the nature of convective driving have been resolved, and are discussed. △ Less

Submitted 27 October, 2010; originally announced October 2010.

Comments: 8 pages, 2 figures, IAU Symposium 271 "Astrophysical Dynamics: From Galaxies to Stars", Nice, FR, 2010

Presupernova Structure of Massive Stars

Authors: Casey A. Meakin, Tuguldur Sukhbold, David Arnett

Abstract: Issues concerning the structure and evolution of core collapse progenitor stars are discussed with an emphasis on interior evolution. We describe a program designed to investigate the transport and mixing processes associated with stellar turbulence, arguably the greatest source of uncertainty in progenitor structure, besides mass loss, at the time of core collapse. An effort to use precision obse… ▽ More Issues concerning the structure and evolution of core collapse progenitor stars are discussed with an emphasis on interior evolution. We describe a program designed to investigate the transport and mixing processes associated with stellar turbulence, arguably the greatest source of uncertainty in progenitor structure, besides mass loss, at the time of core collapse. An effort to use precision observations of stellar parameters to constrain theoretical modeling is also described. △ Less

Submitted 2 June, 2010; originally announced June 2010.

Comments: Proceedings for invited talk at High Energy Density Laboratory Astrophysics conference, Caltech, March 2010. Special issue of Astrophysics and Space Science, submitted for peer review: 7 pages, 3 figures

Some Properties of the Kinetic Energy Flux and Dissipation in Turbulent Stellar Convection Zones

Authors: Casey Meakin, David Arnett

Abstract: We investigate simulated turbulent flow within thermally driven stellar convection zones. Different driving sources are studied, including cooling at the top of the convectively unstable region, as occurs in surface convection zones; and heating at the base by nuclear burning. The transport of enthalpy and kinetic energy, and the distribution of turbulent kinetic energy dissipation are studied.… ▽ More We investigate simulated turbulent flow within thermally driven stellar convection zones. Different driving sources are studied, including cooling at the top of the convectively unstable region, as occurs in surface convection zones; and heating at the base by nuclear burning. The transport of enthalpy and kinetic energy, and the distribution of turbulent kinetic energy dissipation are studied. We emphasize the importance of global constraints on shaping the quasi-steady flow characteristics, and present an analysis of turbulent convection which is posed as a boundary value problem that can be easily incorporated into standard stellar evolution codes for deep, efficient convection. Direct comparison is made between the theoretical analysis and the simulated flow and very good agreement is found. Some common assumptions traditionally used to treat quasi-steady turbulent flow in stellar models are briefly discussed. The importance and proper treatment of convective boundaries are indicated. △ Less

Submitted 23 January, 2010; originally announced January 2010.

Comments: Refereed Conference Proceeding, HELAS Workshop on 'Synergies between solar and stellar modelling', Rome, June 2009, to be published in Astrophys. Space Sci

Turbulent Mixing in Stars: Theoretical Hurdles

Authors: W. David Arnett, Casey Meakin

Abstract: A program is outlined, and first results described, in which fully three-dimensional, time dependent simulations of hydrodynamic turbulence are used as a basis for theoretical investigation of the physics of turbulence in stars. The inadequacy of the treatment of turbulent convection as a diffusive process is discussed. A generalization to rotation and magnetohydrodynamics is indicated, as are c… ▽ More A program is outlined, and first results described, in which fully three-dimensional, time dependent simulations of hydrodynamic turbulence are used as a basis for theoretical investigation of the physics of turbulence in stars. The inadequacy of the treatment of turbulent convection as a diffusive process is discussed. A generalization to rotation and magnetohydrodynamics is indicated, as are connection to simulations of 3D stellar atmospheres. △ Less

Submitted 15 December, 2009; originally announced December 2009.

Comments: 5 pages, 1 figure, IAU Symposium 265, 2009

LSST Science Book, Version 2.0

Authors: LSST Science Collaboration, Paul A. Abell, Julius Allison, Scott F. Anderson, John R. Andrew, J. Roger P. Angel, Lee Armus, David Arnett, S. J. Asztalos, Tim S. Axelrod, Stephen Bailey, D. R. Ballantyne, Justin R. Bankert, Wayne A. Barkhouse, Jeffrey D. Barr, L. Felipe Barrientos, Aaron J. Barth, James G. Bartlett, Andrew C. Becker, Jacek Becla, Timothy C. Beers, Joseph P. Bernstein, Rahul Biswas, Michael R. Blanton, Joshua S. Bloom , et al. (223 additional authors not shown)

Abstract: A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south… ▽ More A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy. △ Less

Submitted 1 December, 2009; originally announced December 2009.

Comments: 596 pages. Also available at full resolution at http://www.lsst.org/lsst/scibook

Convection Theory and Sub-photospheric Stratification

Authors: W. David Arnett, Casey Meakin, Patrick A. Young

Abstract: As a step toward a complete theoretical integration of 3D compressible hydrodynamic simulations into stellar evolution, convection at the surface and sub-surface layers of the Sun is re-examined, from a restricted point of view, in the language of mixing-length theory (MLT) . Requiring that MLT use a hydrodynamically realistic dissipation length gives a new constraint on solar models. While the… ▽ More As a step toward a complete theoretical integration of 3D compressible hydrodynamic simulations into stellar evolution, convection at the surface and sub-surface layers of the Sun is re-examined, from a restricted point of view, in the language of mixing-length theory (MLT) . Requiring that MLT use a hydrodynamically realistic dissipation length gives a new constraint on solar models. While the stellar structure which results is similar to that obtained by YREC and Garching models, the theoretical picture differs. A new quantitative connection is made between macro-turbulence, micro-turbulence, and the convective velocity scale at the photosphere, which has finite values. The "geometric parameter" in MLT is found to correspond more reasonably with the size of the strong downward plumes which drive convection (Stein and Nordlund 1998), and thus has a physical interpretation even in MLT. Use of 3D simulations of both adiabatic convection and stellar atmospheres will allow the determination of the dissipation length and the geometric parameter (i.e., the entropy jump), with no astronomical calibration. A physically realistic treatment of convection in stellar evolution will require additional modifications beyond MLT, including effects of kinetic energy flux, entrainment (the most dramatic difference from MLT found by Meakin and Arnett 2007), rotation, and magnetic fields (Balbus 2009}. △ Less

Submitted 5 October, 2009; originally announced October 2009.

Comments: 16 pages, 5 figures, submitted to ApJ

Journal ref: Astrophys.J.710:1619-1626,2010

A faint type of supernova from a white dwarf with a helium-rich companion

Authors: H. B. Perets, A. Gal-Yam, P. Mazzali, D. Arnett, D. Kagan, A. V. Filippenko, W. Li, I. Arcavi, S. B. Cenko, D. B. Fox, D. C. Leonard, D. -S. Moon, D. J. Sand, A. M. Soderberg, R. J. Foley, M. Ganeshalingam, J. P. Anderson, P. A. James, E. O. Ofek, L. Bildsten, G. Nelemans, K. J. Shen, N. N. Weinberg, B. D. Metzger, A. L. Piro , et al. (3 additional authors not shown)

Abstract: Supernovae (SNe) are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their explosion. These are thought to appear as as type Ib/c and II SNe, and are associated with young stellar populations. A type Ia SN is thought to arise from the thermonuclear detonation of a white… ▽ More Supernovae (SNe) are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their explosion. These are thought to appear as as type Ib/c and II SNe, and are associated with young stellar populations. A type Ia SN is thought to arise from the thermonuclear detonation of a white dwarf star composed mainly of carbon and oxygen, whose mass approaches the Chandrasekhar limit. Such SNe are observed in both young and old stellar environments. Here we report our discovery of the faint type Ib SN 2005E in the halo of the nearby isolated galaxy, NGC 1032. The lack of any trace of recent star formation near the SN location (Fig. 1), and the very low derived ejected mass (~0.3 M_sun), argue strongly against a core-collapse origin for this event. Spectroscopic observations and the derived nucleosynthetic output show that the SN ejecta have high velocities and are dominated by helium-burning products, indicating that SN 2005E was neither a subluminous nor a regular SN Ia (Fig. 2). We have therefore found a new type of stellar explosion, arising from a low-mass, old stellar system, likely involving a binary with a primary white dwarf and a helium-rich secondary. The SN ejecta contain more calcium than observed in any known type of SN and likely additional large amounts of radioactive 44Ti. Such SNe may thus help resolve fundamental physical puzzles, extending from the composition of the primitive solar system and that of the oldest stars, to the Galactic production of positrons. △ Less

Submitted 24 May, 2010; v1 submitted 10 June, 2009; originally announced June 2009.

Comments: Revised to reflect published version in Nature, May 20th, 2010. Additional data and analysis are included

Journal ref: Nature 465:322-325,2010

The Reddening Towards Cassiopeia A's Supernova: Constraining the 56Ni Yield

Authors: Kristoffer A. Eriksen, David Arnett, Donald W. McCarthy, Patrick Young

Abstract: We present new reddening measurements towards the young supernova remnant Cassiopeia A, using two techniques not previously applied to this object. Our observations of the near-infrared [Fe II] 1.257 micron and 1.644 micron lines show the extinction to be highly variable across the remnant, increasing towards the west and the south, consistent with previous radio and X-ray observations. While th… ▽ More We present new reddening measurements towards the young supernova remnant Cassiopeia A, using two techniques not previously applied to this object. Our observations of the near-infrared [Fe II] 1.257 micron and 1.644 micron lines show the extinction to be highly variable across the remnant, increasing towards the west and the south, consistent with previous radio and X-ray observations. While the absolute value of AV as determined by the [Fe II] lines is uncertain due to conflicting calculations and observations of their intrinsic flux ratio, parts of the remnant without previous optical measurements show comparatively higher reddening. We find AV = 6.2 +/- 0.6 from the broadband shape of the infrared synchrotron emission of a knot within 13" of the expansion center. Given this reddening, the apparent faintness of the nascent supernova, and iron mass constraints from X-ray observations, we estimate an ejected mass of 56Ni of 0.058 - 0.16 MSun. Taken with gamma-ray observations of the 44Ti decay chain, this nickel mass is broadly consistent with the solar 44Ca/56Fe ratio. △ Less

Submitted 23 February, 2009; originally announced February 2009.

Comments: 8 pages, 6 figures, accepted for publication in The Astrophysical Journal

Journal ref: Astrophys.J.697:29-36,2009