• ### Improving 1D Stellar Models with 3D Atmospheres(1610.07323)

Oct. 24, 2016 astro-ph.SR
Stellar evolution codes play a major role in present-day astrophysics, yet they share common issues. In this work we seek to remedy some of those by the use of results from realistic and highly detailed 3D hydrodynamical simulations of stellar atmospheres. We have implemented a new temperature stratification extracted directly from the 3D simulations into the Garching Stellar Evolution Code to replace the simplified atmosphere normally used. Secondly, we have implemented the use of a variable mixing-length parameter, which changes as a function of the stellar surface gravity and temperature -- also derived from the 3D simulations. Furthermore, to make our models consistent, we have calculated new opacity tables to match the atmospheric simulations. Here, we present the modified code and initial results on stellar evolution using it.
• ### Surface-effect corrections for the solar model(1606.01030)

June 3, 2016 astro-ph.SR
Solar p-mode oscillations exhibit a systematic offset towards higher frequencies due to shortcomings in the 1D stellar structure models, especially, the lack of turbulent pressure in the superadiabatic layers just below the optical surface, arising from the convective velocity field. We study the influence of the turbulent expansion, chemical composition, and magnetic fields on the stratification in the upper layers of the solar models in comparison with solar observations. Furthermore, we test alternative <3D> averages for improved results on the oscillation frequencies. We appended temporally and spatially averaged <3D> stratifications to 1D models to compute adiabatic oscillation frequencies that we then tested against observations. We also developed depth-dependent corrections for the solar 1D model, for which we expanded the geometrical depth to match the pressure stratification of the solar <3D> model, and we reduced the density that is caused by the turbulent pressure. We obtain the same results with our <3D> models as have been reported previously. Our depth-dependent corrected 1D models match the observations to almost a similar extent as the <3D> model. We find that correcting for the expansion of the geometrical depth and the reducing of the density are both equally necessary. Interestingly, the influence of the adiabatic exponent Gam1 is less pronounced than anticipated. The turbulent elevation directly from the <3D> model does not match the observations properly. Considering different reference depth scales for the <3D> averaging leads to very similar frequencies. Solar models with high metal abundances in their initial chemical composition match the low-frequency part much better. We find a linear relation between the p-mode frequency shift and the vertical magnetic field strength with dvnl = 26.21Bz microHz/kG, which is able to render the solar activity cycles correctly.
• ### On the relevance of bubbles and potential flows for stellar convection(1601.05811)

Jan. 21, 2016 astro-ph.SR
Recently Pasetto et al. have proposed a new method to derive a convection theory appropriate for the implementation in stellar evolution codes. Their approach is based on the simple physical picture of spherical bubbles moving within a potential flow in dynamically unstable regions, and a detailed computation of the bubble dynamics. Based on this approach the authors derive a new theory of convection which is claimed to be parameter free, non-local and time-dependent. This is a very strong claim, as such a theory is the holy grail of stellar physics. Unfortunately we have identified several distinct problems in the derivation which ultimately render their theory inapplicable to any physical regime. In addition we show that the framework of spherical bubbles in potential flows is unable to capture the essence of stellar convection, even when equations are derived correctly.
• ### The Stagger-grid: A grid of 3D stellar atmosphere models - III. The relation to mixing-length convection theory(1403.1062)

Oct. 6, 2014 astro-ph.SR
We investigate the relation between 1D atmosphere models that rely on the mixing length theory and models based on full 3D radiative hydrodynamic (RHD) calculations to describe convection in the envelopes of late-type stars. The adiabatic entropy value of the deep convection zone, s_bot, and the entropy jump, {\Delta}s, determined from the 3D RHD models, are matched with the mixing length parameter, {\alpha}_MLT, from 1D hydrostatic atmosphere models with identical microphysics (opacities and equation-of-state). We also derive the mass mixing length, {\alpha}_m, and the vertical correlation length of the vertical velocity, C[v_z,v_z], directly from the 3D hydrodynamical simulations of stellar subsurface convection. The calibrated mixing length parameter for the Sun is {\alpha}_MLT (s_bot) = 1.98. For different stellar parameters, {\alpha}_MLT varies systematically in the range of 1.7 - 2.4. In particular, {\alpha}_MLT decreases towards higher effective temperature, lower surface gravity and higher metallicity. We find equivalent results for {\alpha}_MLT ({\Delta}s). Also, we find a tight correlation between the mixing length parameter and the inverse entropy jump. We derive an analytical expression from the hydrodynamic mean field equations that motivates the relation to the mass mixing length, {\alpha}_m, and find that it exhibits qualitatively a similar variation with stellar parameter (between 1.6 and 2.4) with a solar value of {\alpha}_m = 1.83. The vertical correlation length scaled with the pressure scale height yields for the Sun 1.71, but displays only a small systematic variation with stellar parameters, the correlation length slightly increasing with Teff. We derive mixing length parameters for various stellar parameters that can be used to replace a constant value. Within any convective envelope, {\alpha}_m and related quantities vary a lot.
• ### Neutrino and axion bounds from the globular cluster M5 (NGC 5904)(1311.1669)

Nov. 7, 2013 hep-ph, astro-ph.CO, astro-ph.SR
The red-giant branch (RGB) in globular clusters is extended to larger brightness if the degenerate helium core loses too much energy in "dark channels." Based on a large set of archival observations, we provide high-precision photometry for the Galactic globular cluster M5 (NGC 5904), allowing for a detailed comparison between the observed tip of the RGB with predictions based on contemporary stellar evolution theory. In particular, we derive 95% confidence limits of $g_{ae}<4.3\times10^{-13}$ on the axion-electron coupling and $\mu_\nu<4.5\times10^{-12}\,\mu_{\rm B}$ (Bohr magneton $\mu_{\rm B}=e/2m_e$) on a neutrino dipole moment, based on a detailed analysis of statistical and systematic uncertainties. The cluster distance is the single largest source of uncertainty and can be improved in the future.
• ### Particle-physics constraints from the globular cluster M5: Neutrino Dipole Moments(1308.4627)

Stellar evolution is modified if energy is lost in a "dark channel" similar to neutrino emission. Comparing modified stellar evolution sequences with observations provides some of the most restrictive limits on axions and other hypothetical low-mass particles and on non-standard neutrino properties. In particular, a putative neutrino magnetic dipole moment mu_nu enhances the plasmon decay process, postpones helium ignition in low-mass stars, and therefore extends the red-giant branch (RGB) in globular clusters (GCs). The brightness of the tip of the RGB (TRGB) remains the most sensitive probe for mu_nu and we revisit this argument from a modern perspective. Based on a large set of archival observations, we provide high-precision photometry for the Galactic GC M5 (NGC5904) and carefully determine its TRGB position. On the theoretical side, we add the extra plasmon decay rate brought about by mu_nu to the Princeton-Goddard-PUC stellar evolution code. Different sources of uncertainty are critically examined. The main source of systematic uncertainty is the bolometric correction and the main statistical uncertainty derives from the distance modulus based on main-sequence fitting. (Other measures of distance, e.g., the brightness of RR Lyrae stars, are influenced by the energy loss that we wish to constrain.) The statistical uncertainty of the TRGB position relative to the brightest RGB star is less important because the RGB is well populated. We infer an absolute I-band brightness of M_I=-4.17+/-0.13 mag for the TRGB compared with the theoretical prediction of -3.99+/-0.07 mag, in reasonable agreement with each other. A significant brightness increase caused by neutrino dipole moments is constrained such that mu_nu<2.6x10^-12mu_B(68% CL), where mu_B is the Bohr magneton, and mu_nu<4.5x10^-12 mu_B(95% CL). In these results, statistical and systematic errors have been combined in quadrature.
• ### S-process in low-mass extremely metal-poor stars(1308.2224)

Aug. 9, 2013 astro-ph.SR
Extremely metal-poor low-mass stars experience an ingestion of protons into the helium-rich layer during the core He-flash, resulting in the production of neutrons through the reactions ^{12}C(p,\gamma)^{13}N(\beta)^{13}C(\alpha,n)^{16}O. This is a potential site for the production of s-process elements in extremely metal-poor stars not occurring in more metal-rich counterparts. Observationally, the signatures of s-process elements in the two most iron deficient stars observed to date, HE1327-2326 & HE0107-5240, still await for an explanation. We investigate the possibility that low-mass EMP stars could be the source of s-process elements observed in extremely iron deficient stars, either as a result of self-enrichment or in a binary scenario as the consequence of a mass transfer episode. We present evolutionary and post-processing s-process calculations of a 1Msun stellar model with metallicities Z=0, 10^{-8} and 10^{-7}. We assess the sensitivity of nucleosynthesis results to uncertainties in the input physics of the stellar models, particularly regarding the details of convective mixing during the core He-flash. Our models provide the possibility of explaining the C, O, Sr, and Ba abundance for the star HE0107-5240 as the result of mass-transfer from a low-mass EMP star. The drawback of our model is that if mass would be transferred before the primary star enters the AGB phase, nitrogen would be overproduced and the ^{12}C/^{13}C abundance ratio would be underproduced in comparison to the observed values. Our results show that low-mass EMP stars cannot be ruled out as the companion stars that might have polluted HE1327-2326 & HE0107-5240 and produced the observed s-process pattern. However, more detailed studies of the core He-flash and the proton ingestion episode are needed to determine the robustness of our predictions.
• ### Helium burning in moderate-mass stars(1108.4634)

Aug. 24, 2011 astro-ph.SR
The evolution of low- and intermediate mass stars at the onset and during core helium burning is reviewed. Particular emphasis is laid on structural differences, which may allow to identify a star's nature and evolutionary phase in spite of the fact that it is found in a region of the Hertzsprung-Russell-Diagram objects from both mass ranges may populate. Seismic diagnostics which are sensitive to the temperature and density profile at the border of the helium core and outside of it may be the most promising tool.
• ### The age of the Milky Way halo stars from the Sloan Digital Sky Survey(1105.2022)

July 21, 2011 astro-ph.GA
We determined the age of the stellar content of the Galactic halo by considering main-sequence turn-off stars. From the large number of halo stars provided by Sloan Digital Sky Survey, we could accurately detect the turn-off as a function of metallicity of the youngest dominant population, which was done by looking at the hottest (bluest) stars of a population. Using the turn-off temperature of a population of a given metallicity, we looked for the isochrones with that turn-off temperature and metallicity and found no age gradient as a function of metallicity. This would mean that this dominating population of the Galactic halo formed rapidly, probably during the collapse of the proto-Galactic gas. Moreover, we could find a significant number of stars with hotter temperatures than the turn-off, which might be blue horizontal branch (BHB) stars, blue stragglers, or main sequence stars that are younger than the dominant population and were probably formed in external galaxies and accreted later on to our Milky Way. Motivated by the current debate about the efficiency of gravitational settling (atomic diffusion) in the interior of old solar-type stars, we used isochrones with and without settling to determine the ages. When ignoring diffusion in the isochrones we obtained ages of 14-16 Gyr. This result is a strong argument against inhibited diffusion in old halo field stars, since it results in a serious conflict with the age of the Universe of 13.7 Gyr. The age obtained including diffusion in the isochrones was 10-12 Gyr, which agrees with the absolute age of the old globular clusters in the inner halo.
• ### Stellar atmosphere parameters with MAx, a MAssive compression of x^2 for spectral fitting(1001.1223)

May 4, 2010 astro-ph.GA, astro-ph.IM
MAx is a new tool to estimate parameters from stellar spectra. It is based on the maximum likelihood method, with the likelihood compressed in a way that the information stored in the spectral fluxes is conserved. The compressed data are given by the size of the number of parameters, rather than by the number of flux points. The optimum speed-up reached by the compression is the ratio of the data set to the number of parameters. The method has been tested on a sample of low-resolution spectra from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey for the estimate of metallicity, effective temperature and surface gravity, with accuracies of 0.24 dex, 130K and 0.5 dex, respectively. Our stellar parameters and those recovered by the SEGUE Stellar Parameter Pipeline agree reasonably well. A small sample of high-resolution VLT-UVES spectra were also used to test the method and the results have been compared to a more classical approach. The speed and multi-resolution capability of MAx combined with its performance compared with other methods indicates that it will be a useful tool for the analysis of upcoming spectral surveys.
• ### On using the CMD morphology of M67 to test solar abundances(1004.3308)

April 19, 2010 astro-ph.SR
The open cluster M67 has solar metallicity and an age of about 4Gyr. The turn-off mass is close to the minimum mass for which solar metallicity stars develop a convective core during main sequence evolution as a result of the development of hydrogen burning through the CNO-cycle. The morphology of the color-magnitude-diagram (CMD) of M67 around the turn-off shows a clear hook-like feature, direct sign that stars close to the turn-off have convective cores. VandenBerg et al. investigated the possibility of using the morphology of the M67 turn-off to put constraints on the solar metallicity, particularly CNO elements, for which solar abundances have been revised downwards by more than 30% over the last few years. Here, we extend their work filling in the gaps in their analysis. To this aim, we compute isochrones appropriate for M67 using new (low metallicity) and old (high metallicity) solar abundances and study whether the characteristic turn-off in the CMD of M67 can be reproduced or not. We also study the importance of other constitutive physics on determining the presence of such a hook, particularly element diffusion, overshooting and nuclear reaction rates. We find that using the new solar abundance determinations, with low CNO abundances, makes it more difficult to reproduce the characteristic CMD of M67. This result is in agreement with results by VandenBerg et al. However, changes in the constitutive physics of the models, particularly overshooting, can influence and alter this result to the extent that isochrones constructed with models using low CNO solar abundances can also reproduce the turn-off morphology in M67. We conclude that only if all factors affecting the turn-off morphology are completely under control (and this is not the case), M67 could be used to put constraints on solar abundances.
• ### Detection of solar-like oscillations from Kepler photometry of the open cluster NGC 6819(1001.0026)

Dec. 30, 2009 astro-ph.SR
Asteroseismology of stars in clusters has been a long-sought goal because the assumption of a common age, distance and initial chemical composition allows strong tests of the theory of stellar evolution. We report results from the first 34 days of science data from the Kepler Mission for the open cluster NGC 6819 -- one of four clusters in the field of view. We obtain the first clear detections of solar-like oscillations in the cluster red giants and are able to measure the large frequency separation and the frequency of maximum oscillation power. We find that the asteroseismic parameters allow us to test cluster-membership of the stars, and even with the limited seismic data in hand, we can already identify four possible non-members despite their having a better than 80% membership probability from radial velocity measurements. We are also able to determine the oscillation amplitudes for stars that span about two orders of magnitude in luminosity and find good agreement with the prediction that oscillation amplitudes scale as the luminosity to the power of 0.7. These early results demonstrate the unique potential of asteroseismology of the stellar clusters observed by Kepler.
• ### New Asymptotic Giant Branch models for a range of metallicities(0903.2155)

Aug. 5, 2009 astro-ph.SR
We present a new grid of stellar model calculations for stars on the Asymptotic Giant Branch between 1.0 and 6.0 M_sun. Our grid consists of 5 chemical mixtures between Z=0.0005 and Z=0.04, with both solar-like and $\alpha$-element enhanced metal ratios. We treat consistently the carbon-enhancement of the stellar envelopes by using opacity tables with varying C/O-ratio and by employing theoretical mass loss rates for carbon stars. The low temperature opacities have been calculated specifically for this project. For oxygen stars we use an empirical mass loss formalism. The third dredge-up is naturally obtained by including convective overshooting. Our models reach effective temperatures in agreement with earlier synthetic models, which included approximative carbon-enriched molecular opacities and show good agreement with empirically determined carbon-star lifetimes. A fraction of the models could be followed into the post-AGB phase, for which we provide models in a mass range supplementing previous post-AGB calculations. Our grid constitutes the most extensive set of AGB-models, calculated with the latest physical input data and treating carbon-enhancement due to the third dredge-up most consistently.
• ### Semi-empirical white dwarf initial-final mass relationships: a thorough analysis of systematic uncertainties due to stellar evolution models(0807.3567)

Oct. 23, 2008 astro-ph
Using the most recent results about white dwarfs in 10 open clusters, we revisit semi-empirical estimates of the initial-final mass relation in star clusters, with emphasis on the use of stellar evolution models. We discuss the influence of these models on each step of the derivation. One intention of our work is to use consistent sets of calculations both for the isochrones and the white dwarf cooling tracks. The second one is to derive the range of systematic errors arising from stellar evolution theory. This is achieved by using different sources for the stellar models and by varying physical assumptions and input data. We find that systematic errors, including the determination of the cluster age, are dominating the initial mass values, while observational uncertainties influence the final mass primarily. After having determined the systematic errors, the initial-final mass relation allows us finally to draw conclusions about the physics of the stellar models, in particular about convective overshooting.
• ### CoRoT/ESTA-TASK 1 and TASK 3 comparison of the internal structure and seismic properties of representative stellar models: Comparisons between the ASTEC, CESAM, CLES, GARSTEC and STAROX codes(0801.0928)

Jan. 7, 2008 astro-ph
We compare stellar models produced by different stellar evolution codes for the CoRoT/ESTA project, comparing their global quantities, their physical structure, and their oscillation properties. We discuss the differences between models and identify the underlying reasons for these differences. The stellar models are representative of potential CoRoT targets. Overall we find very good agreement between the five different codes, but with some significant deviations. We find noticeable discrepancies (though still at the per cent level) that result from the handling of the equation of state, of the opacities and of the convective boundaries. The results of our work will be helpful in interpreting future asteroseismology results from CoRoT.
• ### Age-luminosity relations for low-mass metal-poor stars(astro-ph/9912310)

May 16, 2000 astro-ph
We present a grid of evolutionary calculations for metal-poor low-mass stars for a variety of initial helium and metal abundances. The intention is mainly to provide a database for deriving directly stellar ages of halo and globular cluster stars for which basic stellar parameters are known, but the tracks can also be used for isochrone or luminosity function construction, since they extend to the tip of the red giant branch. Fitting formulae for age-luminosity relations are provided as well. The uncertainties of the evolutionary ages due to inherent shortcomings in the models and due to the unclear effectiveness of diffusion are discussed. A first application to field single stars is presented.
• ### Evolution and surface abundances of red giants experiencing deep mixing(astro-ph/0002290)

Feb. 14, 2000 astro-ph
We have calculated the evolution of low metallicity red giant stars under the assumption of deep mixing between the convective envelope and the hydrogen burning shell. We find that the extent of the observed abundance anomalies, and in particular the universal O-Na anticorrelation, can be totally explained by mixing which does not lead to significant helium enrichment of the envelope. On the other hand, models with extremely deep mixing and strong helium enrichment predict anomalies of sodium and oxygen, which are much larger than the observed ones. This latter result depends solely on the nucleosynthesis inside the hydrogen burning shell, but not on the details of the mixing descriptions. These, however, influence the evolution of surface abundances with brightness, which we compare with the limited observational material available. Our models allow, nevertheless, to infer details on the depth and speed of the mixing process in several clusters. Models with strong helium enrichment evolve to high luminosities and show an increased mass loss. However, under peculiar assumptions, red giants reach very high luminosities even without extreme helium mixing. Due to the consequently increased mass loss, such models could be candidates for blue horizontal branch stars, and, at the same time, would be consistent with the observed abundance anomalies.
• ### Colour transformations for isochrones in the VI-plane(astro-ph/9904236)

April 19, 1999 astro-ph
Reliable temperature-colour transformations are a necessary ingredient for isochrones to be compared with observed colour-magnitude-diagrams of globular clusters. We show that both theoretical and empirical published transformations to (V-I) to a large extend exhibit significant differences between them. Based on these comparisons we argue for particular transformations for dwarfs and giants to be preferred. We then show that our selected combination of transformations results in fits of V-(V-I)-CMDs of a quality which is comparable to that of our earlier V-(B-V) isochrones for a wide range of cluster metallicities. The cluster parameters, such as reddening, are consistent with those derived in (B-V). Therefore, at least in the case of the fit with our own isochrones - based on the particular distance scale provided by our own horizontal branch models, and on the treatment of convection by the mixing-length theory having l/H_p calibrated on our solar model - the chosen transformations appear to lead to self-consistent (V-I) isochrones. Our isochrones are now well tested and self-consistent for B, V and I photometric data.
• ### Helioseismological constraint on solar axion emission(hep-ph/9807476)

Dec. 3, 1998 astro-ph, hep-ph