Search

Helium-burning stars, in particular Cepheids, are especially difficult to model, as the choice of free parameters can greatly impact the shape of the blue loops - the part of the evolutionary track at which instability strip is crossed. Contemporary one-dimensional stellar evolution codes, like MESA (Modules for Experiments in Stellar Astrophysics), come with a large number of free parameters that allow to model the physical processes in stellar interiors under many assumptions. The uncertainties that arise from this freedom are rarely discussed in the literature despite their impact on the evolution of the model. We calculate a grid of evolutionary models with MESA, varying several controls, like solar mixture of heavy elements, mixing length theory prescription, nuclear reaction rates, scheme to determine convective boundaries, atmosphere model, temporal and spatial resolution, and quantify their impact on age and location of the evolutionary track on the HR diagram from the main sequence till the end of core-helium burning. Our investigation was conducted for a full range of masses and metallicities expected for classical Cepheids. The uncertainties are significant, especially during core-helium burning, reaching or exceeding the observational uncertainties of log Teff and log L for detached eclipsing binary systems. For models below 9 solar masses, thin convective shells develop and evolve erratically, not allowing the models to converge. A careful inspection of Kippenhahn diagrams and convergence study is advised for a given mass and metallicity, to assess how severe this problem is and to what extent it may affect the evolution.

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

This first VELOCE data release comprises 18,225 high-precision RV measurements of 258 bona fide classical Cepheids on both hemispheres collected mainly between 2010 and 2022, alongside 1161 additional observations of 164 other stars. The median per-observation RV uncertainty is 0.037 km/s, and some reach 0.002 km/s. Non-variable standard stars characterize RV zero-point stability and provide a base for future cross-calibrations. We determined zero-point differences between VELOCE and 31 literature data sets using template fitting and measured linear period changes of 146 Cepheids. Seventy six spectroscopic binary Cepheids and 14 candidates are identified using VELOCE data alone and are investigated in detail in a companion paper (VELOCE II). Several new insights into Cepheid pulsations were obtained, including: a) the most detailed description of the Hertzsprung progression by RVs; b) the identification of double-peaked bumps in the RV curve; c) clear evidence that virtually all Cepheids feature spectroscopic variability signals that lead to modulated RV variability. We identified 36 such stars, of which 4 also exhibit orbital motion. Linear radius variations depend strongly on pulsation period and a steep increase in slope of the $\Delta$R/p versus logP-relation is found near 10d, challenging the existence of a tight relation between Baade-Wesselink projection factors and pulsation periods. We investigated the accuracy of RV time series measurements, v$_\gamma$, and RV amplitudes published in Gaia's DR3 and determined an average offset of 0.65 \pm 0.11 km/s relative to VELOCE. We recommend adopting a single set of template correlation parameters for distinct classes of large-amplitude variable stars to avoid systematic offsets in v$_\gamma$ among stars belonging to the same class. Peak-to-peak amplitudes of Gaia RVs exhibit significant (16%) dispersion compared to VELOCE. [abridged], Comment: A&A in press, 46 pages, 35 figures, 20 tables. Some data only available via the CDS at publication. VELOCE DR1 data will be made public in FITS format via zenodo.org at DOI: 10.5281/zenodo.10793507 upon publication of the paper in A&A

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron degenerate ignition events. We strengthen MESA's implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator split nuclear burning mode. We close by discussing major updates to MESA's software infrastructure that enhance source code development and community engagement., Comment: 50 pages, 29 figures; Accepted to ApJS

We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models without $^{22}$Ne are fit to the observed periods and other properties. The root-mean-square residuals to the $\simeq$ 150$-$400 s low-order g-mode periods are typically in the range of $\sigma_{rms}$ $\lesssim$ 0.3 s, for a fit precision of $\sigma_{rms}/ P$ $\lesssim$ 0.3 %. We find average relative period shifts of $\Delta P/P$ $\simeq$ $\pm$ 0.5 % for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of $\Delta P/P$ depending on the specific g-mode, abundance of $^{22}$Ne, effective temperature, and mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when $^{22}$Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of $^{22}$Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including $^{22}$Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model fitting processes., Comment: 23 pages, 19 figures total. Accepted to the ApJ. Presented results at 2021 AAS January meeting

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development., Comment: 57 pages, 57 figures; Accepted to ApJS

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction., Comment: 7 pages, 5 figures; Accepted for publication in ApJ Letters

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact., Comment: 64 pages, 61 figures; Accepted to AAS Journals

Our goal is to perform in-depth ensemble asteroseismology of the young open cluster NGC2244 with the 2-wheel Kepler mission. While the nominal Kepler mission already implied a revolution in stellar physics for solar-type stars and red giants, it was not possible to perform asteroseismic studies of massive OB stars because such targets were carefully avoided in the FoV in order not to disturb the exoplanet hunting. Now is an excellent time to fill this hole in mission capacity and to focus on the metal factories of the Universe, for which stellar evolution theory is least adequate. Our white paper aims to remedy major shortcomings in the theory of stellar structure and evolution of the most massive stars by focusing on a large ensemble of stars in a carefully selected young open cluster. Cluster asteroseismology of very young stars such as those of NGC2244 has the major advantage that all cluster stars have similar age, distance and initial chemical composition, implying drastic restrictions for the stellar modeling compared to asteroseismology of single isolated stars with very different ages and metallicities. Our study requires long-term photometric measurements of stars with visual magnitude ranging from 6.5 to 15 in a large FoV with a precision better than 30 ppm for the brightest cluster members (magnitude below 9) up to 500 ppm for the fainter ones, which is well achievable with 2-Wheel Kepler, in combination with high-precision high-resolution spectroscopy and spectro-polarimetry of the brightest pulsating cluster members. These ground-based spectroscopic data will be assembled with the HERMES and CORALIE spectrographs (twin 1.2m Mercator and Euler telescopes, La Palma, Canary Islands and La Silla, Chile), as well as with the spectro-polarimetric NARVAL instrument (2m BLT at the Pic du Midi, French Pyrenees), to which we have guaranteed access., Comment: 10 pages, 3 figures, white paper submitted in response to the NASA call for community input for science investigations the Kepler 2-Wheel spacecraft

We present the CoRoT light curve of the bright B2.5V star HD 48977 observed during a short run of the mission in 2008, as well as a high-resolution spectrum gathered with the HERMES spectrograph at the Mercator telescope. We use several time series analysis tools to explore the nature of the variations present in the light curve. We perform a detailed analysis of the spectrum of the star to determine its fundamental parameters and its element abundances. We find a large number of high-order g-modes, and one rotationally induced frequency. We find stable low-amplitude frequencies in the p-mode regime as well. We conclude that HD 48977 is a new Slowly Pulsating B star with fundamental parameters found to be Teff = 20000 $\pm$ 1000 K and log(g)=4.2 $/pm$ 0.1. The element abundances are similar to those found for other B stars in the solar neighbourhood. HD 48977 was observed during a short run of the CoRoT satellite implying that the frequency precision is insufficient to perform asteroseismic modelling of the star. Nevertheless, we show that a longer time series of this star would be promising for such modelling. Our present study contributes to a detailed mapping of the instability strips of B stars in view of the dominance of g-mode pulsations in the star, several of which occur in the gravito-inertial regime., Comment: 7 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

Photoionization by the high-redshift ultraviolet radiation background heats low density gas before it falls into dark matter potential wells, and it eliminates the neutral hydrogen and singly ionized helium that dominate cooling of primordial gas at temperatures of $10^4-10^5$K. We investigate the influence of photoionization on galaxy formation using high-resolution simulations with a 1-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code. We find that the presence of a photoionizing background suppresses the formation of galaxies with circular velocities $v_{circ} \simlt 30\,$\kms\ and substantially reduces the mass of cooled baryons in systems with circular velocities up to $v_{circ} \sim 50\,$\kms. Above $v_{circ} \sim 75\,$\kms, photoionization has no significant effect. Photoionization exerts its influence primarily by heating gas before collapse; the elimination of line cooling processes is less important. We discuss the implications of these results for hierarchical theories of galaxy formation., Comment: 22 pages, uuencoded compressed postscript with figures. Submitted to ApJ. Latest version also available at http://cfata2.harvard.edu/thoul/papers.html

We investigate the early formation of bound objects with masses comparable to the cosmological Jeans mass (10^5 solar masses). We follow the growth of isolated spherically symmetric density peaks starting from the linear perturbative regime. The initial parameters correspond to density peaks of various widths and heights in a Cold Dark Matter cosmology. We use a one-dimensional spherical Lagrangian hydrodynamics code to follow the dynamical, thermal, and non-equilibrium chemical evolution of the gas. The system includes a collisionless dark matter component and a baryonic component composed of the nine species H, H^-, H^+, He, He^+, He^{++}, H_2, H_2^+, and e^-. All relevant chemical reactions between these species and their cooling mechanisms are included in the calculations. We find that radiative cooling by H_2 affects the collapse dynamics of the gas only after it has already virialized and become part of the bound object. Therefore, radiative cooling is unlikely to have triggered the initial collapse of perturbations at redshifts z>10. Nevertheless, objects with baryonic masses well below the linear-theory Jeans mass (<10^3 solar masses) collapse due to shell crossing by the dark matter. Such objects could be the progenitors of a primordial population of high-mass stars in the intergalactic medium., Comment: 40 pages, uuencoded compressed Postscript, 14 figures included as three separate files

We study the diffusion of helium and other heavy elements in the solar interior by solving exactly the set of flow equations developed by Burgers for a multi-component fluid, including the residual heat-flow terms. No approximation is made concerning the relative concentrations and no restriction is placed on the number of elements considered. We give improved diffusion velocities for hydrogen, helium, oxygen and iron, in the analytic form derived previously by Bahcall and Loeb. These expressions for the diffusion velocities are simple to program in stellar evolution codes and are expected to be accurate to $\sim 15\%$. Our complete treatment of element diffusion can be directly incorporated in a standard stellar evolution code by means of an exportable subroutine, but, for convenience, we also give simple analytical fits to our numerical results., Comment: TeX document, 25 pages, for hardcopy with figures contact BEST@IASSNS.BITNET. Institute for Advanced Study number AST 93/17