Marc H. Pinsonneault

Solar Models: Current Epoch and Time Dependences, Neutrinos, and Helioseismological Properties

We calculate accurate solar models and report the detailed time dependences of important solar quantities. We use helioseismology to constrain the luminosity evolution of the Sun and report the discovery of semiconvection in evolved solar models that include diffusion. In addition, we compare the computed sound speeds with the results of p-mode observations by BiSON, GOLF, GONG, LOWL, and MDI instruments. We contrast the neutrino predictions from a set of eight standard-like solar models and four deviant (or deficient) solar models with the results of solar neutrino experiments. For solar neutrino and helioseismological applications, we present present-epoch numerical tabulations of characteristics of the standard solar model as a function of solar radius, including the principal physical and composition variables, sound speeds, neutrino fluxes, and functions needed for calculating solar neutrino oscillations.

Solar models with helium and heavy element diffusion

Helium and heavy-element diffusion are both included in precise calculations of solar models. In addition, improvements in the input data for solar interior models are described for nuclear reaction rates, the solar luminosity, the solar age, heavy-element abundances, radiative opacities, helium and metal diffusion rates, and neutrino interaction cross sections. The effects on the neutrino fluxes of each change in the input physics are evaluated separately by constructing a series of solar models with one additional improvement added at each stage. The effective 1 σ uncertainties in the individual input quantities are estimated and used to evaluate the uncertainties in the calculated neutrino fluxes and the calculated event rates for solar neutrino experiments. The calculated neutrino event rates, including all of the improvements, are 9.3-1.4+1.2 SNU for the 37Cl experiment and 137-7+8 SNU for the 71Ga experiments. The calculated flux of 7Be neutrinos is 5.1 (1.00-0.07+0.06)×10^9 cm^-2 s^-1 and the flux of 8B neutrinos is 6.6(1.00-0.17+0.14)×10^6 cm^-2 s^-1. The primordial helium abundance found for this model is Y=0.278. The present-day surface abundance of the model is Ys=0.247, in agreement with the helioseismological measurement of Ys=0.242±0.003 determined by Hernandez and Christensen-Dalsgaard (1994). The computed depth of the convective zone is R=0.712R⊙, in agreement with the observed value determined from p-mode oscillation data of R=0.713±0.003R⊙ found by Christensen-Dalsgaard et al. (1991). Although the present results increase the predicted event rate in the four operating solar neutrino experiments by almost 1 σ (theoretical uncertainty), they only slightly increase the difficulty of explaining the existing experiments with standard physics (i.e., by assuming that nothing happens to the neutrinos after they are created in the center of the sun). For an extreme model in which all diffusion (helium and heavy-element diffusion) is neglected, the event rates are 7.0-1.0+0.9 SNU for the 37Cl experiment and 126-6+6 SNU for the 71Ga experiments, while the 7Be and 8B neutrino fluxes are, respectively, 4.5(1.00-0.07+0.06)×10^9 cm^-2 s^-1 and 4.9(1.00-0.17+0.14)×10^6 cm^-2 s^-1. For the no-diffusion model, the computed value of the depth of the convective zone is R=0.726R⊙, which disagrees with the observed helioseismological value. The calculated surface abundance of helium, Ys=0.268, is also in disagreement with the p-mode measurement. The authors conclude that helioseismology provides strong evidence for element diffusion and therefore for the somewhat larger solar neutrino event rates calculated in this paper.

Evolutionary models of the rotating sun

This paper reviews current work on the evolution of a differentially rotating solar model. Although we discuss global features of the evolution with rotation in general terms, the specific models described are those computed with the new Yale Rotating Evolution Code (YREC). YREC uses the Kippenhahn and Thomas (1970, KT) formalism as implemented by Endal and Sofia (1976), although the numerical formulation of our code is totally new. Particular calculations that we describe include the effects of different initial total angular momentum, the consequences of varying the properties and magnitude of angular momentum losses by wind torquing, and the consequences of specific composition and angular momentum redistribution mechanisms. This paper is a progress report which points out the complexity of the problem, and the need for a broad-based observational program to solve it. Because the final solution is not yet in hand, we outline the steps that, in our estimation, need to be undertaken in order to make progress.

Revised stellar properties of Kepler targets for the quarter 1-16 transit detection run

We present revised properties for 196,468 stars observed by the NASA Kepler mission and used in the analysis of Quarter 1-16 (Q1-Q16) data to detect and characterize transiting planets. The catalog is based on a compilation of literature values for atmospheric properties (temperature, surface gravity, and metallicity) derived from different observational techniques (photometry, spectroscopy, asteroseismology, and exoplanet transits), which were then homogeneously fitted to a grid of Dartmouth stellar isochrones. We use broadband photometry and asteroseismology to characterize 11,532 Kepler targets which were previously unclassified in the Kepler Input Catalog (KIC). We report the detection of oscillations in 2762 of these targets, classifying them as giant stars and increasing the number of known oscillating giant stars observed by Kepler by ~20% to a total of ~15,500 stars. Typical uncertainties in derived radii and masses are ~40% and ~20%, respectively, for stars with photometric constraints only, and 5%-15% and ~10% for stars based on spectroscopy and/or asteroseismology, although these uncertainties vary strongly with spectral type and luminosity class. A comparison with the Q1-Q12 catalog shows a systematic decrease in radii of M dwarfs, while radii for K dwarfs decrease or increase depending on the Q1-Q12 provenance (KIC or Yonsei-Yale isochrones). Radii of F-G dwarfs are on average unchanged, with the exception of newly identified giants. The Q1-Q16 star properties catalog is a first step toward an improved characterization of all Kepler targets to support planet-occurrence studies.

Standard solar model

A set of solar models have been constructed, each based on a single modification to the physics of a reference solar model. In addition, a model combining several of the improvements has been calculated to provide a best solar model. Improvements were made to the nuclear reaction rates, the equation of state, the opacities, and the treatment of the atmosphere. The impact on both the structure and the frequencies of the low-l p-modes of the model to these improvements are discussed. It is found that the combined solar model, which is based on the best physics available (and does not contain any ad hoc assumptions), reproduces the observed oscillation spectrum (for low-l) within the errors associated with the uncertainties in the model physics (primarily opacities).

How uncertain are solar neutrino predictions

Abstract Solar neutrino fluxes and sound speeds are calculated using a systematic reevaluation of nuclear fusion rates. The largest uncertainties are identified and their effects on the solar neutrino fluxes are estimated.

Ensemble asteroseismology of solar-type stars with the NASA Kepler mission.

Measurements of 500 Sun-like stars show that their properties differ from those predicted by stellar population models. In addition to its search for extrasolar planets, the NASA Kepler mission provides exquisite data on stellar oscillations. We report the detections of oscillations in 500 solar-type stars in the Kepler field of view, an ensemble that is large enough to allow statistical studies of intrinsic stellar properties (such as mass, radius, and age) and to test theories of stellar evolution. We find that the distribution of observed masses of these stars shows intriguing differences to predictions from models of synthetic stellar populations in the Galaxy.

Helioseismological Implications of Recent Solar Abundance Determinations

We show that standard solar models are in good agreement with the helioseismologically determined sound speed and density as a function of solar radius, the depth of the convective zone, and the surface helium abundance, as long as those models do not incorporate the most recent heavy-element abundance determinations. However, sophisticated new analyses of the solar atmosphere infer lower abundances of the lighter metals (like C, N, O, Ne, and Ar) than the previously widely used surface abundances. We show that solar models that include the lower heavy-element abundances disagree with the solar profiles of sound speed and density as well as the depth of the convective zone and the helium abundance. The disagreements for models with the new abundances range from factors of several to many times the quoted uncertainties in the helioseismological measurements. The disagreements are at temperatures that are too low to affect significantly solar neutrino emission. If errors in the calculated OPAL opacities are solely responsible for the disagreements, then the corrections in the opacity must extend from 2 × 106 K (R = 0.7 R☉) to 5 × 106 K (R = 0.4 R☉), with opacity increases of the order of 10%.

The Problem of Hipparcos Distances to Open Clusters. I. Constraints from Multicolor Main-Sequence Fitting

Parallax data from the Hipparcos mission allow the direct distance to open clusters to be compared with the distance inferred from main-sequence (MS) fitting. There are surprising differences between the two distance measurements, indicating either the need for changes in the cluster compositions or reddening, underlying problems with the technique of MS fitting, or systematic errors in the Hipparcos parallaxes at the 1 mas level. We examine the different possibilities, focusing on MS fitting in both metallicity-sensitive B-V and metallicity-insensitive V-I for five well-studied systems (the Hyades, Pleiades, α Per, Praesepe, and Coma Ber). The Hipparcos distances to the Hyades and α Per are within 1 σ of the MS-fitting distance in B-V and V-I, while the Hipparcos distances to Coma Ber and the Pleiades are in disagreement with the MS-fitting distance at more than the 3 σ level. There are two Hipparcos measurements of the distance to Praesepe; one is in good agreement with the MS-fitting distance and the other disagrees at the 2 σ level. The distance estimates from the different colors are in conflict with one another for Coma but in agreement for the Pleiades. Changes in the relative cluster metal abundances, age related effects, helium, and reddening are shown to be unlikely to explain the puzzling behavior of the Pleiades. We present evidence for spatially dependent systematic errors at the 1 mas level in the parallaxes of Pleiades stars. The implications of this result are discussed.

SEISMIC EVIDENCE FOR A RAPIDLY ROTATING CORE IN A LOWER-GIANT-BRANCH STAR OBSERVED WITH KEPLER

Rotation is expected to have an important influence on the structure and the evolution of stars. However, the mechanisms of angular momentum transport in stars remain theoretically uncertain and very complex to take into account in stellar models. To achieve a better understanding of these processes, we desperately need observational constraints on the internal rotation of stars, which until very recently was restricted to the Sun. In this paper, we report the detection of mixed modes—i.e., modes that behave both as g modes in the core and as p modes in the envelope—in the spectrum of the early red giant KIC 7341231, which was observed during one year with the Kepler spacecraft. By performing an analysis of the oscillation spectrum of the star, we show that its non-radial modes are clearly split by stellar rotation and we are able to determine precisely the rotational splittings of 18 modes. We then find a stellar model that reproduces very well the observed atmospheric and seismic properties of the star. We use this model to perform inversions of the internal rotation profile of the star, which enables us to show that the core of the star is rotating at least five times faster than the envelope. This will shed new light on the processes of transport of angular momentum in stars. In particular, this result can be used to place constraints on the angular momentum coupling between the core and the envelope of early red giants, which could help us discriminate between the theories that have been proposed over the last few decades.