We present a new model for computing the spectral evolution of stellar populations at ages between 100,000 yr and 20 Gyr at a resolution of 3 A across the whole wavelength range from 3200 to 9500 A for a wide range of metallicities. These predictions are based on a newly available library of observed stellar spectra. We also compute the spectral evolution across a larger wavelength range, from 91 A to 160 micron, at lower resolution. The model incorporates recent progress in stellar evolution theory and an observationally motivated prescription for thermally-pulsing stars on the asymptotic giant branch. The latter is supported by observations of surface brightness fluctuations in nearby stellar populations. We show that this model reproduces well the observed optical and near-infrared colour-magnitude diagrams of Galactic star clusters of various ages and metallicities. Stochastic fluctuations in the numbers of stars in different evolutionary phases can account for the full range of observed integrated colours of star clusters in the Magellanic Clouds. The model reproduces in detail typical galaxy spectra from the Early Data Release (EDR) of the Sloan Digital Sky Survey (SDSS). We exemplify how this type of spectral fit can constrain physical parameters such as the star formation history, metallicity and dust content of galaxies. Our model is the first to enable accurate studies of absorption-line strengths in galaxies containing stars over the full range of ages. Using the highest-quality spectra of the SDSS EDR, we show that this model can reproduce simultaneously the observed strengths of those Lick indices that do not depend strongly on element abundance ratios [abridged].

https://hal.archives-ouvertes.fr/hal-00012948/document

Stellar population synthesis at the resolution of 2003

Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source, robust, efficient, thread-safe libraries for a wide range of applications in computational stellar astrophysics. A one-dimensional stellar evolution module, MESAstar, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very low mass to massive stars, including advanced evolutionary phases. MESAstar solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. State-of-the-art modules provide equation of state, opacity, nuclear reaction rates, element diffusion data, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own explicitly defined public interface to facilitate independent development. Several detailed examples indicate the extensive verification and testing that is continuously performed and demonstrate the wide range of capabilities that MESA possesses. These examples include evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets to very old ages; the complete evolutionary track of a 1 M ☉ star from the pre-main sequence (PMS) to a cooling white dwarf; the solar sound speed profile; the evolution of intermediate-mass stars through the He-core burning phase and thermal pulses on the He-shell burning asymptotic giant branch phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; the complete evolutionary tracks of massive stars from the PMS to the onset of core collapse; mass transfer from stars undergoing Roche lobe overflow; and the evolution of helium accretion onto a neutron star. MESA can be downloaded from the project Web site (http://mesa.sourceforge.net/).

https://iopscience.iop.org/article/10.1088/0067-0049/192/1/3/pdf

Modules for Experiments in Stellar Astrophysics (MESA)

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

/pdf/modules-for-experiments-in-stellar-astrophysics-mesa-planets-29j324a9bq.pdf

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

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

https://iopscience.iop.org/article/10.1088/0067-0049/220/1/15/pdf

Modules for experiments in stellar astrophysics (mesa): binaries, pulsations, and explosions

There is almost universal agreement among astronomers that most of the mass in the Universe and most of the mass in the Galactic halo is dark. Many lines of reasoning suggest that the dark matter consists of some new, as yet undiscovered, weakly-interacting massive particle (WIMP). There is now a vast experimental effort being surmounted to detect WIMPS in the halo. The most promising techniques involve direct detection in low-background laboratory detectors and indirect detection through observation of energetic neutrinos from annihilation of WIMPs that have accumulated in the Sun and/or the Earth. Of the many WIMP candidates, perhaps the best motivated and certainly the most theoretically developed is the neutralino, the lightest superpartner in many supersymmetric theories. We review the minimal supersymmetric extension of the Standard Model and discuss prospects for detection of neutralino dark matter. We review in detail how to calculate the cosmological abundance of the neutralino and the event rates for both direct- and indirect-detection schemes, and we discuss astrophysical and laboratory constraints on supersymmetric models. We isolate and clarify the uncertainties from particle physics, nuclear physics, and astrophysics that enter at each step in the calculation. We briefly review other related dark-matter candidates and detection techniques.

Supersymmetric Dark Matter

In the light of recent significant progress on both the observational and theoretical fronts, we review the status of white dwarf stars as cosmochronometers. These objects represent the end products of stellar evolution for the vast majority of stars and, as such, can be used to constrain the ages of various populations of evolved stars in the Galaxy. For example, the oldest white dwarfs in the solar neighborhood (the remnants of the very first generation of intermediate-mass stars in the Galactic disk) are still visible and can be used, in conjunction with cooling theory, to estimate the age of the disk. More recent observations suggest the tantalizing possibility that a population of very old white dwarfs inhabits the Galactic halo. Such a population may contribute significantly to baryonic "dark" matter in the Milky Way and may be used to obtain an independent estimate of the age of the halo. In addition, white dwarf cosmochronology is likely to play a very significant role in the coming era of giant 8-10 m telescopes when faint white dwarf populations should be routinely discovered and studied in open and globular clusters.

/pdf/the-potential-of-white-dwarf-cosmochronology-348kctcvos.pdf

The Potential of White Dwarf Cosmochronology

An age of 9.3 + or - 2.0 Gyr is derived for the Galactic disk on the basis of comparisons between the sudden drop in the observed luminosity distribution and theoretical evolutionary white dwarf models and allowance for a mean prewhite-dwarf lifetime of 0.3 Gyr. To obtain the age of the universe, the time between the big bang and the first appearance of stars in the Galactic disk is added. The age of the universe is estimated to be 10.3 + or - 2.2 Gyr.

/pdf/an-independent-method-for-determining-the-age-of-the-3c8lahsc17.pdf

An independent method for determining the age of the universe

On presente une methode approximee de calcul des coefficients de diffusion pour des plasmas denses caracteristiques des enveloppes de naines blanches

Diffusion coefficients for stellar plasmas

We present new calculations that strongly reinforce the idea—originally proposed by Charpinet et al.—that pulsation modes are driven through an opacity bump due to a local enhancement of the iron abundance in the envelopes of sdB stars. Our improved models incorporate nonuniform iron abundance distributions obtained through the condition of diffusive equilibrium between gravitational settling and radiative levitation. They also include special Rosseland opacity tables that take into account the large variations of the iron abundance about the cosmic value that are predicted by equilibrium radiative levitation theory. For representative models with M = 0.48 M☉ and log g = 5.8, we find strong instabilities for low-order radial and nonradial (p and f) pulsation modes in the range 36,500 K Teff 29,000 K. The four pulsating sdB stars currently known all have effective temperatures in that range. In addition, one of our models with Teff = 34,000 K has a band of unstable modes with periods in the range 116-195 s, in excellent agreement with those of the known pulsators. We therefore claim that our proposed iron bump mechanism provides a natural explanation for the instabilities found in the newly discovered class of pulsating sdB stars.

Gilles Fontaine

Papers

The Potential of White Dwarf Cosmochronology

An independent method for determining the age of the universe

Diffusion coefficients for stellar plasmas

A Driving Mechanism for the Newly Discovered Class of Pulsating Subdwarf B Stars

Optical and Ultraviolet Analyses of ZZ Ceti Stars and Study of the Atmospheric Convective Efficiency in DA White Dwarfs