On the seismic age and heavy-element abundance of the Sun

We estimate the main-sequence age and heavy-element abundance of the Sun by means of an asteroseismic calibration of theoretical solar models using only low-degree acoustic modes from the BiSON. The method can therefore be applied also to other solar-type stars, such as those observed by the NASA satellite Kepler and the planned ground-based Danish-led SONG network. The age, 4.60+/-0.04 Gy, obtained with this new seismic method, is similar to, although somewhat greater than, today's commonly adopted values, and the surface heavy-element abundance by mass, Zs=0.0142+/-0.0005, lies between the values quoted recently by Asplund et al. (2009) and by Caffau et al. (2009). We stress that our best-fitting model is not a seismic model, but a theoretically evolved model of the Sun constructed with `standard' physics and calibrated against helioseismic data.Comment: 16 pages, 11 figures, 5 tables, accepted for publication in MNRA

Progress report on solar age calibration

We report on an ongoing investigation into a seismic calibration of solar models designed for estimating the main-sequence age and a measure of the chemical abundances of the Sun. Only modes of low degree are employed, so that with appropriate modification the procedure could be applied to other stars. We have found that, as has been anticipated, a separation of the contributions to the seismic frequencies arising from the relatively smooth, glitch-free, background structure of the star and from glitches produced by helium ionization and the abrupt gradient change at the base of the convection zone renders the procedure more robust than earlier calibrations that fitted only raw frequencies to glitch-free asymptotics. As in the past, we use asymptotic analysis to design seismic signatures that are, to the best of our ability, contaminated as little as possible by those uncertain properties of the star that are not directly associated with age and chemical composition. The calibration itself, however, employs only numerically computed eigenfrequencies. It is based on a linear perturbation from a reference model. Two reference models have been used, one somewhat younger, the other somewhat older than the Sun. The two calibrations, which use BiSON data, are more-or-less consistent, and yield a main-sequence age $t_\odot=4.68\pm0.02$Gy, coupled with a formal initial heavy-element abundance $Z=0.0169\pm0.0005$. The error analysis has not yet been completed, so the estimated precision must be taken with a pinch of salt.Comment: 8 pages, 3 figures, in L. Deng, K.L. Chan, C. Chiosi, eds, The Art of Modelling Stars in the 21st Century, Proc. IAU Symp. No. 252, invited contributed pape

ASTEC -- the Aarhus STellar Evolution Code

The Aarhus code is the result of a long development, starting in 1974, and still ongoing. A novel feature is the integration of the computation of adiabatic oscillations for specified models as part of the code. It offers substantial flexibility in terms of microphysics and has been carefully tested for the computation of solar models. However, considerable development is still required in the treatment of nuclear reactions, diffusion and convective mixing.Comment: Astrophys. Space Sci, in the pres

The internal rotation profile of the B-type star KIC10526294 from frequency inversion of its dipole gravity modes and statistical model comparison

The internal angular momentum distribution of a star is key to determine its evolution. Fortunately, the stellar internal rotation can be probed through studies of rotationally-split non-radial oscillation modes. In particular, detection of non-radial gravity modes (g modes) in massive young stars has become feasible recently thanks to the Kepler space mission. Our aim is to derive the internal rotation profile of the Kepler B8V star KIC 10526294 through asteroseismology. We interpret the observed rotational splittings of its dipole g modes using four different approaches based on the best seismic models of the star and their rotational kernels. We show that these kernels can resolve differential rotation the radiative envelope if a smooth rotational profile is assumed and the observational errors are small. Based on Kepler data, we find that the rotation rate near the core-envelope boundary is well constrained to $163\pm89$ nHz. The seismic data are consistent with rigid rotation but a profile with counter-rotation within the envelope has a statistical advantage over constant rotation. Our study should be repeated for other massive stars with a variety of stellar parameters in order to deduce the physical conditions that determine the internal rotation profile of young massive stars, with the aim to improve the input physics of their models.Comment: 52 pages, 32 figures, accepted for publication in The Astrophysical Journa

High-precision abundances of elements in Kepler LEGACY stars. Verification of trends with stellar age

HARPS-N spectra with S/N > 250 and MARCS model atmospheres were used to derive abundances of C, O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Zn, and Y in ten stars from the Kepler LEGACY sample (including the binary pair 16 Cyg A and B) selected to have metallicities in the range -0.15 < [Fe/H] < +0.15 and ages between 1 and 7 Gyr. Stellar gravities were obtained from seismic data and effective temperatures were determined by comparing non-LTE iron abundances derived from FeI and FeII lines. Available non-LTE corrections were also applied when deriving abundances of the other elements. The results support the [X/Fe]-age relations previously found for solar twins. [Mg/Fe], [Al/Fe], and [Zn/Fe] decrease by ~0.1 dex over the lifetime of the Galactic thin disk due to delayed contribution of iron from Type Ia supernovae relative to prompt production of Mg, Al, and Zn in Type II supernovae. [Y/Mg] and [Y/Al], on the other hand, increase by ~0.3 dex, which can be explained by an increasing contribution of s-process elements from low-mass AGB stars as time goes on. The trends of [C/Fe] and [O/Fe] are more complicated due to variations of the ratio between refractory and volatile elements among stars of similar age. Two stars with about the same age as the Sun show very different trends of [X/H] as a function of elemental condensation temperature Tc and for 16 Cyg, the two components have an abundance difference, which increases with Tc. These anomalies may be connected to planet-star interactions.Comment: 13 pages with 7 figures. Accepted for publication in A&

Prospects for Measuring Differential Rotation in White Dwarfs Through Asteroseismology

We examine the potential of asteroseismology for exploring the internal rotation of white dwarf stars. Data from global observing campaigns have revealed a wealth of frequencies, some of which show the signature of rotational splitting. Tools developed for helioseismology to use many solar p-mode frequencies for inversion of the rotation rate with depth are adapted to the case of more limited numbers of modes of low degree. We find that the small number of available modes in white dwarfs, coupled with the similarity between the rotational-splitting kernels of the modes, renders direct inversion unstable. Accordingly, we adopt what we consider to be plausible functional forms for the differential rotation profile; this is sufficiently restrictive to enable us to carry out a useful calibration. We show examples of this technique for PG 1159 stars and pulsating DB white dwarfs. Published frequency splittings for white dwarfs are currently not accurate enough for meaningful inversions; reanalysis of existing data can provide splittings of sufficient accuracy when the frequencies of individual peaks are extracted via least-squares fitting or multipeak decompositions. We find that when mode trapping is evident in the period spacing of g modes, the measured splittings can constrain dOmega/dr.Comment: 26 pages, 20 postscript figures. Accepted for publication in The Astrophysical Journa

Excitation of solar-like oscillations across the HR diagram

We extend semi-analytical computations of excitation rates for solar oscillation modes to those of other solar-like oscillating stars to compare them with recent observations. Numerical 3D simulations of surface convective zones of several solar-type oscillating stars are used to characterize the turbulent spectra as well as to constrain the convective velocities and turbulent entropy fluctuations in the uppermost part of the convective zone of such stars. These constraints, coupled with a theoretical model for stochastic excitation, provide the rate 'P' at which energy is injected into the p-modes by turbulent convection. These energy rates are compared with those derived directly from the 3D simulations. The excitation rates obtained from the 3D simulations are systematically lower than those computed from the semi-analytical excitation model. We find that Pmax, the excitation rate maximum, scales as (L/M)^s where s is the slope of the power law and L and M are the mass and luminosity of the 1D stellar model built consistently with the associated 3D simulation. The slope is found to depend significantly on the adopted form of the eddy time-correlation ; using a Lorentzian form results in s=2.6, whereas a Gaussian one gives s=3.1. Finally, values of Vmax, the maximum in the mode velocity, are estimated from the computed power laws for Pmax and we find that Vmax increases as (L/M)^sv. Comparisons with the currently available ground-based observations show that the computations assuming a Lorentzian eddy time-correlation yield a slope, sv, closer to the observed one than the slope obtained when assuming a Gaussian. We show that the spatial resolution of the 3D simulations must be high enough to obtain accurate computed energy rates.Comment: 14 pages ; 7 figures ; accepted for publication in Astrophysics & Astronom

Variation in the frequency separations with activity and impact on stellar parameter determination

Frequency separations used to infer global properties of stars through asteroseismology can change depending on the strength and at what epoch of the stellar cycle the p-mode frequencies are measured. In the Sun these variations have been seen, even though the Sun is a low-activity star. In this paper, we discuss these variations and their impact on the determination of the stellar parameters (radius, mass and age) for the Sun. Using the data from maximum and minimum activity, we fitted an age for the Sun that differs on average by 0.2 Gyr: slightly older during minimum activity. The fitted radius is also lower by about 0.5% for the solar effective temperature during minimum.Comment: to be published in JPCS to be published in JPC