Normalized to: Morel, P.
[1]
oai:arXiv.org:1806.10533 [pdf] - 1759330
Impacts of radiative accelerations on solar-like oscillating
main-sequence stars
Submitted: 2018-06-27
Chemical element transport processes are among the crucial physical processes
needed for precise stellar modelling. Atomic diffusion by gravitational
settling nowadays is usually taken into account, and is essential for
helioseismic studies. On the other hand, radiative accelerations are rarely
accounted for, act differently on the various chemical elements, and can
strongly counteract gravity in some stellar mass domains. In this study we aim
at determining whether radiative accelerations impact the structure of
solar-like oscillating main-sequence stars observed by asteroseismic space
missions. We implemented the calculation of radiative accelerations in the
CESTAM code using the Single-Valued Parameter method. We built and compared
several grids of stellar models including gravitational settling, but some with
and others without radiative accelerations. We found that radiative
accelerations may not be neglected for stellar masses larger than
1.1~M$_{\odot}$ at solar metallicity. The difference in age due to their
inclusion in models can reach 9\% for the more massive stars of our grids. We
estimated that the percentage of the PLATO core program stars whose modelling
would require radiative accelerations ranges between 33 and 58\% depending on
the precision of the seismic data. We conclude that, in the context of Kepler,
TESS, and PLATO missions, which provide (or will provide) high quality seismic
data, radiative accelerations can have a significant effect when inferring the
properties of solar-like oscillators properly. This is particularly important
for age inferences. However, the net effect for each individual star results
from the competition between atomic diffusion including radiative accelerations
and other internal transport processes. This will be investigated in a
forthcoming companion paper.
[2]
oai:arXiv.org:1709.03054 [pdf] - 1588159
A semi-empirical model for magnetic braking of solar-type stars
Submitted: 2017-09-10
We develop new angular momentum evolution models for stars with masses of
$0.5$ to $1.6~\rm M_\odot$ and from the pre-main-sequence (\rm PMS) through the
end of their main-sequence (\rm MS) lifetime. The parametric models include
magnetic braking based on numerical simulations of magnetised stellar winds,
mass loss rate prescription, core-envelope decoupling as well as disk locking
phenomena.
We have also accounted for recent developments in modelling dramatically
weakened magnetic braking in stars more evolved than the Sun. We fit the free
parameters in our model by comparing model predictions to rotational
distributions of a number of stellar clusters as well as individual field
stars. Our model reasonably successfully reproduces the rotational behaviour of
stars during the \rm PMS phase to the zero-age main-sequence (\rm ZAMS) spin
up, sudden \rm ZAMS spin down, and convergence of the rotation rates
afterwards. We find that including core-envelope decoupling improves our models
especially for low-mass stars at younger ages. In addition, by accounting for
the almost complete suppression of magnetic braking at slow spin periods, we
provide better fits to observations of stellar rotations compared to previous
models.
[3]
oai:arXiv.org:1708.04937 [pdf] - 1587170
Equation of state SAHA-S meets stellar evolution code CESAM2k
Submitted: 2017-08-16
We present an example of an interpolation code of the SAHA-S equation of
state that has been adapted for use in the stellar evolution code CESAM2k. The
aim is to provide the necessary data and numerical procedures for its
implementation in a stellar code. A technical problem is the discrepancy
between the sets of thermodynamic quantities provided by the SAHA-S equation of
state and those necessary in the CESAM2k computations. Moreover, the
independent variables in a practical equation of state (like SAHA-S) are
temperature and density, whereas for modelling calculations the variables
temperature and pressure are preferable. Specifically for the CESAM2k code,
some additional quantities and their derivatives must be provided. To provide
the bridge between the equation of state and stellar modelling, we prepare
auxiliary tables of the quantities that are demanded in CESAM2k. Then we use
cubic spline interpolation to provide both smoothness and a good approximation
of the necessary derivatives. Using the B-form of spline representation
provides us with an efficient algorithm for three-dimensional interpolation.
The table of B-spline coefficients provided can be directly used during stellar
model calculations together with the module of cubic spline interpolation. This
implementation of the SAHA-S equation of state in the CESAM2k stellar structure
and evolution code has been tested on a solar model evolved to the present. A
comparison with other equations of state is briefly discussed. The choice of a
regular net of mesh points for specific primary quantities in the SAHA-S
equation of state, together with accurate and consistently smooth tabulated
values, provides an effective algorithm of interpolation in modelling
calculations. The proposed module of interpolation procedures can be easily
adopted in other evolution codes.
[4]
oai:arXiv.org:1612.01331 [pdf] - 1533005
Evolution of lithium abundance in the Sun and solar twins
Submitted: 2016-12-05
Evolution of the 7Li abundance in the convection zone of the Sun during
different stages of its life time is considered to explain its low photospheric
value in comparison with that of the solar system meteorites. Lithium is
intensively and transiently burned in the early stages of evolution (pre-main
sequence, pMS) when the radiative core arises, and then the Li abundance only
slowly decreases during the main sequence (MS). We study the rates of lithium
burning during these two stages. In a model of the Sun, computed ignoring pMS
and without extra-convective mixing (overshooting) at the base of the
convection zone, the lithium abundance does not decrease significantly during
the MS life time of 4.6 Gyr. Analysis of helioseismic inversions together with
post-model computations of chemical composition indicates the presence of the
overshooting region and restricts its thickness. It is estimated to be
approximately half of the local pressure scale height 0.5Hp) which corresponds
to 3.8 % of the solar radius. Introducing this extra region does not noticeably
deplete lithium during the MS stage. In contrast, at the pMS stage, an
overshooting region with a value of approximately 0.18Hp is enough to produce
the observed lithium depletion. If we conclude that the dominant lithium
burning takes place during the pMS stage, the dispersion of the lithium
abundance in solar twins is explained by different physical conditions,
primarily during the early stage of evolution before the MS.
[5]
oai:arXiv.org:1410.4780 [pdf] - 936228
Benchmark stars for Gaia: fundamental properties of the Population II
star HD140283 from interferometric, spectroscopic and photometric data
Creevey, O.;
Thévenin, F.;
Berio, P.;
Heiter, U.;
von Braun, K.;
Mourard, D.;
Bigot, L.;
Boyajian, T. S.;
Kervella, P.;
Morel, P.;
Pichon, B.;
Chiavassa, A.;
Nardetto, N.;
Perraut, K.;
Meilland, A.;
Alister, H. A. Mc;
Brummelaar, T. A. ten;
Sturmann, C. Farrington. J.;
Sturmann, L.;
Turner, N.
Submitted: 2014-10-17, last modified: 2014-12-15
We determined the fundamental properties of HD 140283 by obtaining new
interferometric and spectroscopic measurements and combining them with
photometry from the literature. The interferometric measurements were obtained
using the visible interferometer VEGA on the CHARA array and we determined a 1D
limb-darkened angular diameter of 0.353 +/- 0.013 milliarcseconds. Using
photometry from the literature we derived the bolometric flux with two
solutions: a zero-reddening one of Fbol = 3.890 +/- 0.066 1E-8 erg/s/cm2 and a
solution with a maximum of Av = 0.1 mag, Fbol= 4.220 +/- 0.067 1E-8 erg/s/cm2.
The interferometric Teff is thus 5534 +/- 103 K or 5647 +/- 105 K and its
radius is R = 2.21 +/- 0.08 Rsol. Spectroscopic measurements of HD140283 were
obtained using HARPS, NARVAL, and UVES and a 1D LTE analysis of H-alpha line
wings yields Teff(Halpha) = 5626 +/- 75 K. Using fine-tuned stellar models
including diffusion of elements we then determined the mass M and age t of
HD140283. Once the metallicity has been fixed, the age of the star depends on
M, initial helium abundance Yi and mixing-length parameter alpha, only two of
which are independent. We need to adjust alpha to much lower values than the
solar one (~2) in order to fit the observations, and if Av = 0.0 mag then 0.5 <
alpha < 1. We give an equation to estimate t from M, Yi (alpha) and Av.
Establishing a reference alpha = 1.00 and adopting Yi = 0.245 we derive a mass
and age of HD140283: M = 0.780 +/- 0.010 Msol and t = 13.7 +/- 0.7 Gyr (Av =
0.0) or M = 0.805 +/- 0.010 Msol and t = 12.2 +/- 0.6 Gyr (Av=0.1 mag). Our
stellar models yield an initial metallicity of [Z/X]i = -1.70 and logg = 3.65
+/- 0.03. Asteroseismic observations are critical for overcoming limitations in
our results.
[6]
oai:arXiv.org:1409.7477 [pdf] - 1222258
Evolution of Exoplanets and their Parent Stars
Submitted: 2014-09-26
Studying exoplanets with their parent stars is crucial to understand their
population, formation and history. We review some of the key questions
regarding their evolution with particular emphasis on giant gaseous exoplanets
orbiting close to solar-type stars. For masses above that of Saturn, transiting
exoplanets have large radii indicative of the presence of a massive
hydrogen-helium envelope. Theoretical models show that this envelope
progressively cools and contracts with a rate of energy loss inversely
proportional to the planetary age. The combined measurement of planetary mass,
radius and a constraint on the (stellar) age enables a global determination of
the amount of heavy elements present in the planet interior. The comparison
with stellar metallicity shows a correlation between the two, indicating that
accretion played a crucial role in the formation of planets. The dynamical
evolution of exoplanets also depends on the properties of the central star. We
show that the lack of massive giant planets and brown dwarfs in close orbit
around G-dwarfs and their presence around F-dwarfs are probably tied to the
different properties of dissipation in the stellar interiors. Both the
evolution and the composition of stars and planets are intimately linked.
[7]
oai:arXiv.org:1309.7179 [pdf] - 1179540
Models of solar surface dynamics: impact on eigenfrequencies and radius
Submitted: 2013-09-27
We study the effects of different descriptions of the solar surface
convection on the eigenfrequencies of p-modes. 1-D evolution calculations of
the whole Sun and 3-D hydrodynamic and magnetohydrodynamic simulations of the
current surface are performed. These calculations rely on realistic physics.
Averaged stratifications of the 3-D simulations are introduced in the 1-D solar
evolution or in the structure models. The eigenfrequencies obtained are
compared to those of 1-D models relying on the usual phenomenologies of
convection and to observations of the MDI instrument aboard SoHO. We also
investigate how the magnetic activity could change the eigenfrequencies and the
solar radius, assuming that, 3 Mm below the surface, the upgoing plasma advects
a 1.2 kG horizontal field.
All models and observed eigenfrequencies are fairly close below 3 mHz. Above
3 mHz the eigenfrequencies of the phenomenological convection models are above
the observed eigenfrequencies. The frequencies of the models based on the 3-D
simulations are slightly below the observed frequencies. Their maximum
deviation is ~ 3 \mu Hz at 3 mHz but drops below 1 \mu Hz at 4 mHz. Replacing
the hydrodynamic by the magnetohydrodynamic simulation increases the
eigenfrequencies. The shift is negligible below 2.2 mHz and then increases
linearly with frequency to reach ~ 1.7 \mu Hz at 4 mHz. The impact of the
simulated activity is a 14 milliarcsecond shrinking of the solar layers near
the optical depth unity.
[8]
oai:arXiv.org:1309.2009 [pdf] - 1380630
Differential asteroseismic study of seismic twins observed by CoRoT;
Comparison of HD 175272 with HD 181420
Ozel, N.;
Mosser, B.;
Dupret, M. A.;
Bruntt, H.;
Barban, C.;
Deheuvels, S.;
García, R. A.;
Michel, E.;
Samadi, R.;
Baudin, F.;
Mathur, S.;
Régulo, C.;
Auvergne, M.;
Morel, P.;
Pichon, B.
Submitted: 2013-09-08
The CoRoT short asteroseismic runs give us the opportunity to observe a large
variety of late-type stars through their solar-like oscillations. We report the
observation and modeling of the F5V star HD 175272. Our aim is to define a
method for extracting as much information as possible from a noisy oscillation
spectrum. We followed a differential approach that consists of using a
well-known star as a reference to characterize another star. We used classical
tools such as the envelope autocorrelation function to derive the global
seismic parameters of the star. We compared HD 175272 with HD 181420 through a
linear approach, because they appear to be asteroseismic twins. The comparison
with the reference star enables us to substantially enhance the scientific
output for HD 175272. First, we determined its global characteristics through a
detailed seismic analysis of HD 181420. Second, with our differential approach,
we measured the difference of mass, radius and age between HD 175272 and HD
181420. We have developed a general method able to derive asteroseismic
constraints on a star even in case of low-quality data. %This method is based
on the comparison to a star with common seismic and classical properties.
Seismic data allow accurate measurements of radii and masses differences
between the two stars. This method can be applied to stars with interesting
properties but low signal-to-noise ratio oscillation spectrum, such as stars
hosting an exoplanet or members of a binary system.
[9]
oai:arXiv.org:1211.1271 [pdf] - 1157563
Seismic diagnostics for transport of angular momentum in stars 1.
Rotational splittings from the PMS to the RGB
Marques, J. P.;
Goupil, M. J.;
Lebreton, Y.;
Talon, S.;
Palacios, A.;
Belkacem, K.;
Ouazzani, R. -M.;
Mosser, B.;
Moya, A.;
Morel, P.;
Pichon, B.;
Mathis, S.;
Zahn, J. -P.;
Turck-Chièze, S.;
Nghiem, P A. P.
Submitted: 2012-11-06
Rotational splittings are currently measured for several main sequence stars
and a large number of red giants with the space mission Kepler. This will
provide stringent constraints on rotation profiles. Our aim is to obtain
seismic constraints on the internal transport and surface loss of angular
momentum of oscillating solar-like stars. To this end, we study the evolution
of rotational splittings from the pre-main sequence to the red-giant branch for
stochastically excited oscillation modes. We modified the evolutionary code
CESAM2K to take rotationally induced transport in radiative zones into account.
Linear rotational splittings were computed for a sequence of $1.3 M_{\odot}$
models. Rotation profiles were derived from our evolutionary models and
eigenfunctions from linear adiabatic oscillation calculations. We find that
transport by meridional circulation and shear turbulence yields far too high a
core rotation rate for red-giant models compared with recent seismic
observations. We discuss several uncertainties in the physical description of
stars that could have an impact on the rotation profiles. For instance, we find
that the Goldreich-Schubert-Fricke instability does not extract enough angular
momentum from the core to account for the discrepancy. In contrast, an increase
of the horizontal turbulent viscosity by 2 orders of magnitude is able to
significantly decrease the central rotation rate on the red-giant branch. Our
results indicate that it is possible that the prescription for the horizontal
turbulent viscosity largely underestimates its actual value or else a mechanism
not included in current stellar models of low mass stars is needed to slow down
the rotation in the radiative core of red-giant stars.
[10]
oai:arXiv.org:1207.5954 [pdf] - 541890
Fundamental properties of the Population II fiducial stars HD 122563 and
Gmb 1830 from CHARA interferometric observations
Creevey, O. L.;
Thévenin, F.;
Boyajian, T. S.;
Kervella, P.;
Chiavassa, A.;
Bigot, L.;
Mérand, A.;
Heiter, U.;
Morel, P.;
Pichon, B.;
Alister, H. A. Mc;
Brummelaar, T. A. ten;
Collet, R.;
van Belle, G. T.;
Foresto, V. Coudé du;
Farrington, C.;
Goldfinger, P. J.;
Sturmann, J.;
Sturmann, L.;
Turner, N.
Submitted: 2012-07-25
We have determined the angular diameters of two metal-poor stars, HD 122563
and Gmb 1830, using CHARA and Palomar Testbed Interferometer observations. For
the giant star HD 122563, we derive an angular diameter theta_3D = 0.940 +-
0.011 milliarcseconds (mas) using limb-darkening from 3D convection simulations
and for the dwarf star Gmb 1830 (HD 103095) we obtain a 1D limb-darkened
angular diameter theta_1D = 0.679 +- 0.007 mas. Coupling the angular diameters
with photometry yields effective temperatures with precisions better than 55 K
(Teff = 4598 +- 41 K and 4818 +- 54 K --- for the giant and the dwarf star,
respectively). Including their distances results in very well-determined
luminosities and radii (L = 230 +- 6 L_sun, R = 23.9 +- 1.9 R_sun and L = 0.213
+- 0.002 L_sun, R = 0.664 +- 0.015 R_sun, respectively). We used the CESAM2k
stellar structure and evolution code in order to produce models that fit the
observational data. We found values of the mixing-length parameter alpha (which
describes 1D convection) that depend on the mass of the star. The masses were
determined from the models with precisions of <3% and with the well-measured
radii excellent constraints on the surface gravity are obtained (log g = 1.60
+- 0.04, 4.59 +- 0.02, respectively). The very small errors on both log g and
Teff provide stringent constraints for spectroscopic analyses given the
sensitivity of abundances to both of these values. The precise determination of
Teff for the two stars brings into question the photometric scales for
metal-poor stars.
[11]
oai:arXiv.org:1110.0985 [pdf] - 421811
The diameter of the CoRoT target HD 49933. Combining the 3D limb
darkening, asteroseismology, and interferometry
Bigot, L.;
Mourard, D.;
Berio, P.;
Thévenin, F.;
Ligi, R.;
Tallon-Bosc, I.;
Chesneau, O.;
Delaa, O.;
Nardetto, N.;
Perraut, K.;
Stee, Ph.;
Boyajian, T.;
Morel, P.;
Pichon, B.;
Kervella, P.;
Schmider, F. X.;
McAlister, H.;
Brummelaar, T. Ten;
Ridgway, S. T.;
Sturmann, J.;
Sturmann, L.;
Turner, N.;
Farrington, C.;
Goldfinger, P. J.
Submitted: 2011-10-05
Context. The interpretation of stellar pulsations in terms of internal
structure depends on the knowledge of the fundamental stellar parameters.
Long-base interferometers permit us to determine very accurate stellar radii,
which are independent constraints for stellar models that help us to locate the
star in the HR diagram. Aims: Using a direct interferometric determination of
the angular diameter and advanced three-dimensional (3D) modeling, we derive
the radius of the CoRoT target HD 49933 and reduce the global stellar parameter
space compatible with seismic data. Methods: The VEGA/CHARA
spectro-interferometer is used to measure the angular diameter of the star. A
3D radiative hydrodynamical simulation of the surface is performed to compute
the limb darkening and derive a reliable diameter from visibility curves. The
other fundamental stellar parameters (mass, age, and Teff) are found by fitting
the large and small p-mode frequency separations using a stellar evolution
model that includes microscopic diffusion. Results: We obtain a limb-darkened
angular diameter of {\theta}LD = 0.445 \pm 0.012 mas. With the Hipparcos
parallax, we obtain a radius of R = 1.42 \pm 0.04 Rsun. The corresponding
stellar evolution model that fits both large and small frequency separations
has a mass of 1.20 \pm 0.08 Msun and an age of 2.7 Gy. The atmospheric
parameters are Teff = 6640 \pm 100 K, log g = 4.21 \pm 0.14, and [Fe/H] =
-0.38.
[12]
oai:arXiv.org:1002.3461 [pdf] - 1025258
Survival of a convective core in low-mass solar-like pulsator HD 203608
Submitted: 2010-02-18
A 5-night asteroseismic observation of the F8V star HD 203608 was conducted
in August 2006 with HARPS, followed by an analysis of the data, and a
preliminary modeling of the star (Mosser et al. 2008). The stellar parameters
were significantly constrained, but the behavior of one of the seismic indexes
(the small spacing d01) could not be fitted with the observed one, even with
the best considered models. We study the possibility of improving the agreement
between models and observations by changing the physical properties of the
inner parts of the star (to which d01 is sensitive). We show that, in spite of
its low mass, it is possible to produce models of HD 203608 with a convective
core. No such model was considered in the preliminary modeling. In practice, we
obtain these models here by assuming some extra mixing at the edge of the early
convective core. We optimize the model parameters using the Levenberg-Marquardt
algorithm. The agreement between the new best model with a convective core and
the observations is much better than for the models without. All the
observational parameters are fitted within 1-sigma observational error bars.
This is the first observational evidence of a convective core in an old and
low-mass star such as HD 203608. In standard models of low-mass stars, the core
withdraws shortly after the ZAMS. The survival of the core until the present
age of HD 203608 provides very strong constraints on the size of the mixed zone
associated to the convective core. Using overshooting as a proxy to model the
processes of transport at the edge of the core, we find that to reproduce both
global and seismic observations, we must have alpha_ov = 0.17 +/- 0.03 Hp for
HD 203608. We revisit the process of the extension of the core lifetime due to
overshooting in the particular case of HD 203608.
[13]
oai:arXiv.org:0907.1448 [pdf] - 26141
Ledoux's convection criterion in evolution and asteroseismology of
massive stars
Submitted: 2009-07-09
Saio et al. (2006) have shown that the presence of an intermediate convective
zone (ICZ) in post-main sequence models could prevent the propagation of
g-modes in the radiative interior and hence avoid the corresponding radiative
damping. The development of such a convective region highly depends on the
structure of the star in the mu-gradient region surrounding the convective core
during the main sequence phase. In particular,the development of this ICZ
depends on physical processes such as mass loss, overshooting (Chiosi & Maeder
1986, Chiosi et al. 1992, see also Godart et al., these proceedings) and
convective instability criterion (Schwarzschild's or Ledoux's criteria). In
this paper we study the consequences of adopting the Ledoux's criterion on the
evolution of the convective regions in massive stars (15 and 20 Msun), and on
the pulsation spectrum of these new B-type variables (also called SPBsg).
[14]
oai:arXiv.org:0804.2144 [pdf] - 11736
Thorough analysis of input physics in CESAM and CLES codes
Submitted: 2008-04-14
This contribution is not about the quality of the agreement between stellar
models computed by CESAM and CLES codes, but more interesting, on what
ESTA-Task~1 run has taught us about these codes and about the input physics
they use. We also quantify the effects of different implementations of the same
physics on the seismic properties of the stellar models, that in fact is the
main aim of ESTA experiments.
[15]
oai:arXiv.org:0801.2019 [pdf] - 8954
CESAM: a free code for stellar evolution calculations
Submitted: 2008-01-14
The Cesam code is a consistent set of programs and routines which perform
calculations of 1D quasi-hydrostatic stellar evolution including microscopic
diffusion of chemical species and diffusion of angular momentum. The solution
of the quasi-static equilibrium is performed by a collocation method based on
piecewise polynomials approximations projected on a B-spline basis; that allows
stable and robust calculations, and the exact restitution of the solution, not
only at grid points, even for the discontinuous variables. Other advantages are
the monitoring by only one parameter of the accuracy and its improvement by
super-convergence. An automatic mesh refinement has been designed for adjusting
the localisations of grid points according to the changes of unknowns. For
standard models, the evolution of the chemical composition is solved by stiffly
stable schemes of orders up to four; in the convection zones mixing and
evolution of chemical are simultaneous. The solution of the diffusion equation
employs the Galerkin finite elements scheme; the mixing of chemicals is then
performed by a strong turbulent diffusion. A precise restoration of the
atmosphere is allowed for.
[16]
oai:arXiv.org:0704.2294 [pdf] - 335826
Sensitivity of low degree oscillations to the change in solar abundances
Submitted: 2007-04-18, last modified: 2007-04-22
Context. The most recent determination of the solar chemical composition,
using a time-dependent, 3D hydrodynamical model of the solar atmosphere,
exhibits a significant decrease of C, N, O abundances compared to their
previous values. Solar models that use these new abundances are not consistent
with helioseismological determinations of the sound speed profile, the surface
helium abundance and the convection zone depth. Aims. We investigate the effect
of changes of solar abundances on low degree p-mode and g-mode characteristics
which are strong constraints of the solar core. We consider particularly the
increase of neon abundance in the new solar mixture in order to reduce the
discrepancy between models using new abundances and helioseismology. Methods.
The observational determinations of solar frequencies from the GOLF instrument
are used to test solar models computed with different chemical compositions. We
consider in particular the normalized small frequency spacings in the low
degree p-mode frequency range. Results. Low-degree small frequency spacings are
very sensitive to changes in the heavy-element abundances, notably neon. We
show that by considering all the seismic constraints, including the small
frequency spacings, a rather large increase of neon abundance by about (0.5 +/-
0.05)dex can be a good solution to the discrepancy between solar models that
use new abundances and low degree helioseismology, subject to adjusting
slightly the solar age and the highest abundances. We also show that the change
in solar abundances, notably neon, considerably affects g-mode frequencies,
with relative frequency differences between the old and the new models higher
than 1.5%
[17]
oai:arXiv.org:astro-ph/0703426 [pdf] - 90235
Microscopic Diffusion in Stellar Evolution Codes: First Comparison
results of ESTA-Task~3
Lebreton, Y.;
Montalban, J.;
Christensen-Dalsgaard, J.;
Theado, S.;
Hui-Bon-Hoa, A.;
Monteiro, M. J. P. F. G.;
Degl'Innocenti, S.;
Marconi, M.;
Morel, P.;
Moroni, P. G. Prada;
Weiss, A.
Submitted: 2007-03-15
We present recent work undertaken by the Evolution and Seismic Tools Activity
(ESTA) team of the CoRoT Seismology Working Group. The new ESTA-Task 3 aims at
testing, comparing and optimising stellar evolution codes which include
microscopic diffusion of the chemical elements resulting from pressure,
temperature and concentration gradients. The results already obtained are
globally satisfactory, but some differences between the different numerical
tools appear that require further investigations.
[18]
oai:arXiv.org:astro-ph/0605685 [pdf] - 82372
Report on the CoRoT Evolution and Seismic Tools Activity
Monteiro, M. K. J. P. F. G.;
Lebreton, Y.;
Montalban, J.;
Christensen-Dalsgaard, J.;
Castro, M.;
Degl'Innocenti, S.;
Moya, A.;
Roxburgh, I. W.;
Scuflaire, R.;
Baglin, A.;
Cunha, M. S.;
Eggenberger, P.;
Fernandes, J.;
Goupil, M. J.;
Hui-Bon-Hoa, A.;
Marconi, M.;
Marques, J. P.;
Michel, E.;
Miglio, A.;
Morel, P.;
Pichon, B.;
Moroni, P. G. Prada;
Provost, J.;
Ruoppo, A.;
Suarez, J. -C.;
Suran, M.;
Teixeira, T. C.
Submitted: 2006-05-27
We present the work undertaken by the Evolution and Seismic Tools Activity
(ESTA) team of the CoRoT Seismology Working Group. We have focused on two main
tasks: Task 1 - now finished - has aimed at testing, comparing and optimising
seven stellar evolution codes which will be used to model the internal
structure and evolution of the CoRoT target stars. Task 2, still underway, aims
at testing, comparing and optimising different seismic codes used to calculate
the oscillations of models for different types of stars. The results already
obtained are quite satisfactory, showing minor differences between the
different numerical tools provided the same assumptions on the physical
parameters are made. This work gives us confidence on the numerical tools that
will be available to interpret the future CoRoT seismic data.
[19]
oai:arXiv.org:astro-ph/0511602 [pdf] - 77960
The effects of moderately fast shellular rotation on adiabatic
oscillations
Submitted: 2005-11-21, last modified: 2005-11-22
We investigate adiabatic oscillations for delta Scuti star models, taking
into account a moderate rotation velocity ~100 \km/s. The resulting oscillation
frequencies include corrections for rotation up to second order in the rotation
rate including those of near degeneracy. Effects of either a uniform rotation
or a rotation profile assuming local angular momentum conservation of the form
Omega=Omega(r) on oscillation frequencies are compared. As expected, important
differences (around 3 microHz) are obtained in the $g$ and mixed mode regions.
For higher frequency p modes, differences range between 1 microHz and 3
microHz. Such differences are likely to be detectable with future space
missions such as COROT, where precisions in frequency around 0.5 microHz will
be reachable.
[20]
oai:arXiv.org:astro-ph/0501420 [pdf] - 1456403
VLTI/VINCI diameter constraints on the evolutionary status of delta Eri,
xi Hya, eta Boo
Submitted: 2005-01-19
Using VLTI/VINCI angular diameter measurements, we constrain the evolutionary
status of three asteroseismic targets: the stars $\delta$ Eri, $\xi$ Hya,
$\eta$ Boo. Our predictions of the mean large frequency spacing of these stars
are in agreement with published observational estimations. Looking without
success for a companion of $\delta$ Eri we doubt on its classification as an RS
CVn star.
[21]
oai:arXiv.org:astro-ph/0312068 [pdf] - 61281
Modeling of Alpha Cen and Procyon using VLTI observations
Submitted: 2003-12-02, last modified: 2003-12-11
We present a novel approach to model the nearby stars Alpha Cen A & B and
Procyon A using asteroseismic and interferometric constraints. Using the VINCI
instrument installed at the VLT Interferometer (VLTI), the angular diameters of
the Alpha Centauri system were measured with a relative precision of 0.2% and
0.6%, respectively. From these values, we derive linear radii of R[A] = 1.224
+/- 0.003 R_sun and R[B] = 0.863 +/- 0.005 R_sun. These radii are in excellent
agreement with the models of Thevenin et al. (2002), that use asteroseismic
frequencies as constraints (Bouchy & Carrier 2001; Bouchy & Carrier 2002). With
the same instrument, we also measured the angular diameter of Procyon A. Using
the Hipparcos parallax, we obtain a linear radius of 2.048 +/- 0.025 R_sun. We
use this result together with spectroscopic and photometric constraints to
model this star with the CESAM code. We also computed the adiabatic oscillation
spectrum of our model of Procyon A, giving a mean large frequency separation of
Delta nu_0 = 54.8 uHz, in agreement with the seismic observations by Martic et
al. (2001). Our model favours a mass around 1.4 M_sun for Procyon A.
[22]
oai:arXiv.org:astro-ph/0309784 [pdf] - 59611
VINCI / VLTI observations of Main Sequence stars
Kervella, P.;
Thevenin, F.;
Morel, P.;
Provost, J.;
Berthomieu, G.;
Segransan, D.;
Queloz, D.;
Borde, P.;
Di Folco, E.;
Forveille, T.
Submitted: 2003-09-29
Main Sequence (MS) stars are by far the most numerous class in the Universe.
They are often somewhat neglected as they are relatively quiet objects (but
exceptions exist), though they bear testimony of the past and future of our
Sun. An important characteristic of the MS stars, particularly the solar-type
ones, is that they host the large majority of the known extrasolar planets.
Moreover, at the bottom of the MS, the red M dwarfs pave the way to
understanding the physics of brown dwarfs and giant planets. We have measured
very precise angular diameters from recent VINCI/VLTI interferometric
observations of a number of MS stars in the K band, with spectral types between
A1V and M5.5V. They already cover a wide range of effective temperatures and
radii. Combined with precise Hipparcos parallaxes, photometry, spectroscopy as
well as the asteroseismic information available for some of these stars, the
angular diameters put strong constraints on the detailed models of these stars,
and therefore on the physical processes at play.
[23]
oai:arXiv.org:astro-ph/0309148 [pdf] - 1468558
The diameter and evolutionary state of Procyon A
Submitted: 2003-09-04
We report the angular diameter measurement obtained with the VINCI/VLTI
instrument on the nearby star Procyon A (alpha CMi A, F5IV-V), at a relative
precision of +/- 0.9%. We obtain a uniform disk angular diameter in the K band
of UD = 5.376 +/- 0.047 mas and a limb darkened value of LD = 5.448 +/- 0.053
mas. Together with the Hipparcos parallax, this gives a linear diameter of
2.048 +/- 0.025 D_sun. We use this result in combination with spectroscopic,
photometric and asteroseismic constraints to model this star with the CESAM
code. One set of modeling parameters that reproduces the observations within
their error bars are an age of 2314 Myr, an initial helium mass fraction Y_i =
0.301 and an initial mass ratio of heavy elements to hydrogen Z/X_i = 0.0314.
We also computed the adiabatic oscillation spectrum of our model of Procyon A,
giving a mean large frequency separation of Delta nu_0 = 54.7 microHz. This
value is in agreement with the seismic observations by Martic et al. (1999,
2001). The interferometric diameter and the asteroseismic large frequency
spacing together suggest a mass closer to 1.4 M_sun than to 1.5 M_sun. We
conclude that Procyon is currently ending its life on the main sequence, as its
luminosity class indicates.
[24]
oai:arXiv.org:astro-ph/0306604 [pdf] - 57685
The interferometric diameter and internal structure of Sirius A
Submitted: 2003-06-28
The interferometric observations of dwarf stars in the solar neighbourhood,
combined with Hipparcos parallaxes provide very precise values of their linear
diameters. In this paper, we report the direct measurement of the angular
diameter of the bright star Sirius A with the VINCI/VLTI instrument. We obtain
a uniform disk angular diameter of UD = 5.936 +/- 0.016 mas in the K band and a
limb darkened value of LD = 6.039 +/- 0.019 mas. In combination with the
Hipparcos parallax of 379.22 +/- 1.58 mas, this translates into a linear
diameter of 1.711 +/- 0.013 Dsun. Using the VINCI/VLTI interferometric diameter
and the published properties of Sirius A, we derive internal structure models
corresponding to ages between 200 and 250 +/- 12 Myr. This range is defined
mainly by the hypothesis on the mass of the star, the overshoot and the
metallicity.
[25]
oai:arXiv.org:astro-ph/0303634 [pdf] - 363025
The diameters of Alpha Centauri A and B - A comparison of the
asteroseismic and VINCI/VLTI views
Submitted: 2003-03-28, last modified: 2003-03-29
We compare the first direct angular diameter measurements obtained on our
closest stellar neighbour, Alpha Centauri, to recent model diameters
constrained by asteroseismic observations. Using the VINCI instrument installed
at ESO's VLT Interferometer (VLTI), the angular diameters of the two main
components of the system, Alpha Cen A and B, were measured with a relative
precision of 0.2% and 0.6%, respectively. Particular care has been taken in the
calibration of these measurements, considering that VINCI is estimating the
fringe visibility using a broadband K filter. We obtain uniform disk angular
diameters for Alpha Cen A and B of UD[A] = 8.314 +/- 0.016 mas and UD[B] =
5.856 +/- 0.027 mas, and limb darkened angular diameters of LD[A] = 8.511 +/-
0.020 mas and LD[B] = 6.001 +/- 0.034 mas. Combining these values with the
parallax from Soderhjelm (1999), we derive linear diameters of D[A] = 1.224 +/-
0.003 Dsun and D[B] = 0.863 +/- 0.005 Dsun. These values are compatible with
the masses published by Thevenin et al.(2002) for both stars.
[26]
oai:arXiv.org:astro-ph/0206283 [pdf] - 1468419
Asteroseismology and calibration of alpha Cen binary system
Submitted: 2002-06-17
Using the oscillation frequencies of alpha Cen A recently discovered by
Bouchy & Carrier, the available astrometric, photometric and spectroscopic
data, we tried to improve the calibration of the visual binary system alpha
Cen. With the revisited masses of Pourbaix et al. (2002) we do not succeed to
obtain a solution satisfying all the seismic observational constraints.
Relaxing the constraints on the masses, we have found an age t_alpha
Cen=4850+-500 Myr, an initial helium mass fraction Y_i = 0.300+-0.008, and an
initial metallicity (Z/X)_i=0.0459+-0.0019, with M_A=1.100+-0.006M_o and
M_B=0.907+-0.006M_o for alpha Cen A&B.
[27]
oai:arXiv.org:astro-ph/0205434 [pdf] - 49508
Atomic diffusion in star models of type earlier than G
Submitted: 2002-05-24
We find that the radiative diffusivity associated with the radiative
viscosity generates a mixing that diminishes the efficiency of the
gravitational settling in the external layers of stellar models corresponding
to types earlier than G. The surface abundances of chemical species predicted
by the models are successfully compared with the abundances determined in
members of the Hyades open cluster. Our modeling depends on an efficiency
parameter, which is evaluated to a value close to unity, that we calibrate in
this study.
[28]
oai:arXiv.org:astro-ph/0110004 [pdf] - 45062
The Zeta Herculis binary system revisited. Calibration and seismology
Submitted: 2001-09-29
We have revisited the calibration of the visual binary system Zeta Herculis
with the goal to give the seismological properties of the G0 IV sub-giant Zeta
Her A. We have used the most recent physical and observational data. For the
age we have obtained 3387 Myr, for the masses respectively 1.45 and 0.98 solar
mass, for the initial helium mass fraction 0.243, for the initial mass ratio of
heavy elements to hydrogen 0.0269 and for the mixing-length parameters
respectively 0.92 and 0.90 using the Canuto & Mazitelli (1991, 1992) convection
theory. Our results do not exclude that Zeta Her A is itself a binary
sub-system; the mass of the hypothetical unseen companion would be smaller than
0.05 solar mass. The adiabatic oscillation spectrum of Zeta Her A is found to
be a complicated superposition of acoustic and gravity modes; some of them have
a dual character. This greatly complicates the classification of the non-radial
modes. The echelle diagram used by the observers to extract the frequencies
will work for ell=0, 2, 3. The large difference is found to be of the order of
42 mu Hz, in agreement with the Martic et al. (2001) seismic observations.
[29]
oai:arXiv.org:astro-ph/0010543 [pdf] - 38890
Shear turbulence beneath the solar tachocline
Submitted: 2000-10-26
Helioseismic inversions of the Sun's internal angular velocity profile show
that the rotation changes from differential in latitude in the convection zone
to almost uniform in the radiative region below. The transition occurs in a
thin layer, the tachocline, which is the seat of strong shear in the vertical
direction. In this Note we examine whether this rotation profile can lead to
shear turbulence at the top of the radiation zone. By using the standard solar
model, we show that such turbulence can be generated only in a narrow region
0.695R_sun < r < 0.713R_sun at the equator and even in narrower layers at
higher latitudes. We conclude that the turbulence generated by this vertical
shear is unlikely to play a significant role in the transport of matter and
angular momentum, and that other mechanisms must be invoked to achieve this.
[30]
oai:arXiv.org:astro-ph/0010180 [pdf] - 38527
Calibrations of alpha Cen A & B
Submitted: 2000-10-10
Detailed evolutionary models of the visual binary alpha Centauri, including
pre main-sequence evolution, have been performed using the masses recently
determined by Pourbaix et al. (1999). Models have been constructed using the
CEFF equation of state, OPAL opacities, NACRE thermonuclear reaction rates and
microscopic diffusion. A chi^2-minimization is performed to derive the most
reliable set of modeling parameters wp={t_alpha Cen, Yi, [Fe/H]i, alpha A,
alpha B}, where t alpha Cen is the age of the system, Yi the initial helium
content, [Fe/H]i the initial metallicity and, alpha A and alpha B the
convection parameters of the two components. Using the basic Bohm-Vitense
(1958) mixing-length theory of convection, we derive wp BV={2710 Myr,
0.284,0.257, 1.53, 1.57}. We obtain a noticeably smaller age than estimated
previously, in agreement with Pourbaix et al. (1999), mainly because of the
larger masses. If convective core overshoot is considered we get wp ov={3530
Myr, 0.279,0.264,1.64,1.66}. The use of Canuto & Mazitelli (1991, 1992)
convection theory leads to the set wp CM={4086 Myr, 0.271, 0.264, 0.964,
0.986}. Using the observational constraints adopted by Guenther & Demarque
(2000), and the basic mixing-length theory, we obtain wp GD={5640 Myr, 0.300,
0.296, 1.86, 1.97} and surface lithium depletions close to their observed
values.
A seismological analysis of our calibrated models has been performed. The
determination of large and small spacings between the frequencies of acoustic
oscillations from seismic observations would help to discriminate between the
models of alpha Cen computed with different masses and to confirm or rules out
the new determination of masses.
[31]
oai:arXiv.org:astro-ph/9912307 [pdf] - 1469871
A calibration of iota Pegasi system
Submitted: 1999-12-15
Recent observations provide determinations of individual masses, chemical
composition and metallicity of the components of the spectroscopic and
interferometric binary \iota Peg (Boden et al. 1998). Using updated physics, to
calibrate the system, we have computed using the stellar evolutionary code
CESAM, evolutionary sequences of stellar models with the masses of \iota Peg A
1.326 Msol and iota Peg B 0.819 Msol (Boden et al. loc.cit) and with
different values of the mixing-length parameter \alpha, the helium Y and the
heavy element Z initial mass fraction with the constraint of the observed
metallicity. Adopting effective temperatures and luminosities, as derived from
observations with the bolometric corrections, and the empirical scale of
temperatures of Alonso et al. (1995, 1996), we find \alpha_A=1.46,
\alpha_B=1.36, Y=0.278, Z=0.017. The evolution time, including pre-main
sequence, is found within 40 My < t_ev < 0.5Gy. The calibrated models of \iota
Peg. A and B are non homogeneous zero age main sequence stars with the
evolutionary time t_ev=56 My. Due to the large uncertainties of their
determinations, the values derived for the mixing-length parameters are smaller
than the solar one but however marginally compatible with it. Our results ought
to be improved as soon as a more accurate value of the magnitude difference in
the V filter will be available. Detailed spectroscopic analysis for both
components looks practicable, so it is urgently needed.
[32]
oai:arXiv.org:astro-ph/9910388 [pdf] - 108942
About the time of evolution of a Solar Model
Submitted: 1999-10-21
The evolution of a solar model is initialized with homogeneous models of
either, pre-main sequence (P-models) or, zero-age main sequence (Z-models). The
zero-age of a solar model is conventionally referenced as the time where the
nuclear reactions just begin to dominate gravitation as the primary source of
energy. Fixing the physics, we found that the structure of P- and Z-models
computed with the same physics are almost similar soon after the exhaustion of
their convective core. This similarity gives a connection between the age of
the Sun t_sun and the time t_cal elapsed in the calculation of calibrated solar
models. We found that a Z-model calibrated with t_cal=t_sun and a P-model
calibrated with t_cal=t_sun+25My, are indistinguishable at the relative
accuracy level of a few 10^{-4}.
[33]
oai:arXiv.org:astro-ph/9907381 [pdf] - 107606
Solar Models and NACRE thermonuclear reaction rates
Submitted: 1999-07-27
Using the most recent updated physics, calibrated solar models have been
computed with the new thermonuclear reaction rates of NACRE, the recently
available European compilation. Comparisons with models computed with the
reaction rates of Caughlan & Fowler (\cite{cf88}) and of Adelberger et al.
(\cite{a98}) are made for global structure, expected neutrinos fluxes, chemical
composition and sound speed profiles, helioseismological properties of p-modes
and g-modes.
[34]
oai:arXiv.org:astro-ph/9902160 [pdf] - 105179
Microscopic diffusion and subdwarfs
Submitted: 1999-02-10
The recent distance determinations by HIPPARCOS have allowed to locate very
precisely a reasonable set of subdwarfs in the HR diagram. A detailed
comparison with evolutionary models is now possible and it has been recently
claimed that the observed lower part of the main-sequence is cooler than the
computed one. In this paper we show that microscopic diffusion could be the
physical explanation of such a systematic shift. We emphazise the fact that, as
diffusion alters surface abundances, the initial metal content of the stellar
material is higher than the present one. A "calibration" procedure is necessary
when comparing observed and theoretical HR diagram, to take into account this
effect.
[35]
oai:arXiv.org:astro-ph/9902068 [pdf] - 105087
Survival of 6li and 7li in Metal-Poor Stars
Submitted: 1999-02-04
The relationship between the respective depletion factors for 6Li and 7Li in
metal-poor stars is studied for two simple models of depletion by nuclear
burning by mixing with deeper layers, below the bottom of the convective zone.
The depletion of 6Li is usually much larger than the depletion of 7Li because
it occurs during the pre-mainsequence already.
[36]
oai:arXiv.org:astro-ph/9806272 [pdf] - 101865
Standard Solar Models in the Light of New Helioseismic Constraints: I
The Solar Core
Submitted: 1998-06-19, last modified: 1998-06-22
In this paper, we examine a new updated solar model taking advantage of
recent reexamination of the nuclear reactions rates and of the microscopic
diffusion of helium and heavy elements. Our best model fits the helioseismic
data reasonably well: the base of the convective zone is at r=0.715 Ro, the
photospheric helium in mass fraction is 0.243 and the sound speed square
difference between the Sun and the model is < 1%. This model leads to a
reestimate of the neutrinos fluxes : 7.18 SNU for the chlorine experiment,
127.2 SNU for the gallium detector and 4.82 10^6/cm2/s for the B8 neutrino
flux. Acoustic and gravity mode predictions are also estimated. We present 3
models which respect the helioseismic constraints and fit the seismic
observations better by taking advantage of the known physical uncertainties:
nuclear reaction rates, CNO abundances and microscopic diffusion. We conclude
that the standard model in side its inherent uncertainties, is robust in the
light of the present acoustic mode detection and that the neutrino fluxes
prediction is not well constrained.
[37]
oai:arXiv.org:astro-ph/9805091 [pdf] - 101312
Atmospheric structure and acoustic cut-off frequency of roAp stars
Submitted: 1998-05-07
Some of the rapidly oscillating (CP2) stars, have frequencies which are
larger than the theoretical acoustic cut-off frequency. As the cut-off
frequency depends on the T(tau) relation in the atmosphere, we have computed
models and adiabatic frequencies for pulsating Ap stars with T(tau) laws based
on Kurucz model atmospheres and on Hopf's purely radiative relation.
The frequency-dependent treatment of radiative transfer as well as an
improved calculation of the radiative pressure in Kurucz model atmospheres
increase the theoretical acoustic cut-off frequency by about 200 microHz, which
is closer to the observations.
For alpha Cir we find models with Kurucz atmospheres which have indeed a
cut-off frequency beyond the largest observed frequency and which are well
within the Teff - L error box. For HD 24712 only models which are hotter by
about 100 K and less luminous by nearly 10% than what is actually the most
probable value would have an acoustic cut-off frequency large enough.
One may thus speculate that the old controversy about a mismatch between
observed largest frequencies and theoretical cut-off frequencies of roAp star
models is resolved. However, the observational errors for the astrophysical
fundamental parameters have to be reduced further and the model atmospheres
refined.
Further details can be found in Audard et al. (1997)
[38]
oai:arXiv.org:astro-ph/9712126 [pdf] - 99615
The acoustic cut-off frequency of roAp stars
Submitted: 1997-12-09
Some of the rapidly oscillating (roAp) stars, have frequencies which are
larger than the theoretical acoustic cut-off frequency. We have computed models
and adiabatic frequencies for pulsating Ap stars with T(tau) laws based on
Kurucz model atmospheres and on the Hopf's purely radiative relation. We
compare the values of the cut-off frequency derived from expressions for the
potential from Vorontsov & Zarkhov (1989), from Gough (1986), and from the
approximation of an isothermal atmosphere. Kurucz model atmospheres increase
the theoretical acoustic cut-off frequency for roAp stars by about 200 mu$Hz,
which is closer to the observations. We compare our computations with
observations of the 3 roAp stars HD 24712, alphaCir, and HD 134214. For alpha
Cir we find models with Kurucz atmospheres which have indeed a cut-off
frequency beyond the largest observed frequency. For HD 24712 only models which
are hotter by about 100 K and less luminous by nearly 10% than what is actually
the most probable value would have an acoustic cut-off frequency large enough.
One may thus speculate that the old controversy about a mismatch between
observed largest frequencies and theoretical cut-off frequencies of roAp star
models is resolved. However, the model atmospheres have to be refined by
investigating NLTE effects, among others, before a definite conclusion can be
drawn.
[39]
oai:arXiv.org:astro-ph/9710319 [pdf] - 99078
Inferring the equatorial solar tachocline from frequency splittings
Submitted: 1997-10-28, last modified: 1997-10-30
Helioseismic inversions, carried out for several years on various
ground-based and spatial observations, have shown that the solar rotation rate
presents two principal regimes: a quasi-rigid rotation in the radiative
interior and a latitude-dependent rotation in the whole convection zone. The
thin layer, named solar tachocline, between these two regimes is difficult to
infer through inverse techniques because of the ill-posed nature of the problem
that requires regularization techniques which, in their global form, tend to
smooth out any high gradient in the solution. Thus, most of the previous
attempts to study the rotation profile of the solar tachocline have been
carried out through forward modeling. In this work we show that some
appropriate inverse techniques can also be used and we compare the ability of
three 1D inverse techniques combined with two automatic strategies for the
choice of the regularization parameter, to infer the solar tachocline profile
in the equatorial plane.
Our work, applied on LOWL (LOWL is an abbreviation for low degree denoted by
L) two years dataset, argue in favor of a very sharp (0.05+/-0.03R_sun)
transition zone located at 0.695+/-0.005R_sun which is in good agreement with
the previous forward analysis carried out on Global Oscillations Network Group
(GONG), Big Bear Solar Observatory (BBSO) and LOWL datasets.
[40]
oai:arXiv.org:astro-ph/9705251 [pdf] - 97513
Updated Solar Models
Submitted: 1997-05-30
Solar models computed with mass loss, microscopic diffusion of helium and
heavy elements and with updated physics have been evolved from the pre-main
sequence to present days; they are compared to the observational constraints
including lithium depletion and to the Solar Reference Model (SRM) of Basu et
al. (1996), derived by inversion. Microscopic diffusion significantly improves
the agreement with the observed solar frequencies and agree with the SRM within
+- 0.2% for the sound velocity and +- 1% for the density, but slightly worsens
the neutrino problem. Neither microscopic diffusion nor overshooting explain
the observed lithium depletion consistently with helioseismological
constraints, while a mass loss process does it. Models computed with OPAL
equation of state and opacities are in better agreement with the seimic sound
speed. To reach the level of precision of helioseismological observations the
accuracy of solar models still needs to be improved by one magnitude; any such
improvement will necessitate equation of state and opacity data taking into
account of detailed changes in the mixture.
[41]
oai:arXiv.org:astro-ph/9704069 [pdf] - 97031
A Reestimate of the Protosolar 2H/1H ratio from 3He/4He solar wind
measurements
Submitted: 1997-04-09
We reanalyze the inference of the protosolar abundance of deuterium made by
Geiss (1993) from measurements of 3He/4He in the solar wind. We use an
evolutionary solar model with microscopic diffusion, constrained to fit the
present age, radius and luminosity, as well as the observed ratio of heavy
elements to hydrogen. The protosolar 2H/1H is obtained from the best fit of
3He/4He in the solar wind. Taking for the protosolar 3He/4He the value measured
in Jupiter by the Galileo probe (Niemann et al. 1996), we derive the protosolar
2H/1H=(3.01 +- 0.17) 10^{-5} ratio. Compared to the present interstellar medium
value (Linsky et al. 1993), this result is compatible with models of the
chemical evolution of the Galaxy in the solar neighborhood; it is also
marginally compatible with the Jovian 2H/1H=(5 +- 2) 10^{-5} ratio measured by
Galileo.