-
In-depth analysis of solar models with high-metallicity abundances and updated opacity tables
Authors:
G. Buldgen,
A. Noels,
R. Scuflaire,
A. M. Amarsi,
N. Grevesse,
P. Eggenberger,
J. Colgan,
C. J. Fontes,
V. A. Baturin,
A. V. Oreshina,
S. V. Ayukov,
P. Hakel,
D. P. Kilcrease
Abstract:
Due to the high quality constraints available for the Sun, we can carry out combined analyses using neutrino, spectroscopic and helioseismic observations. Such studies lay the ground for future improvements of key physical components of solar and stellar models, such as the equation of state, radiative opacities or prescriptions for macroscopic transport processes of chemicals which are then used…
▽ More
Due to the high quality constraints available for the Sun, we can carry out combined analyses using neutrino, spectroscopic and helioseismic observations. Such studies lay the ground for future improvements of key physical components of solar and stellar models, such as the equation of state, radiative opacities or prescriptions for macroscopic transport processes of chemicals which are then used to study other stars in the Universe. We study the existing degeneracies in solar models using the recent high-metallicity spectroscopic abundances by comparing them to helioseismic and neutrino data and discuss how their properties are impacted by changes in various physical ingredients. We carry out a detailed study of solar models computed with a high-metallicity composition from the literature based on averaged-3D models that was claimed to solve the solar problem. The properties of the solar models are significantly affected by using the recent OPLIB opacities and the inclusion of macroscopic transport. The properties of the standard solar models computed using the OPAL opacities are similar to those using the OP opacities. We show that a modifying the temperature gradient just below the base of the convective zone is required to erase the discrepancies in solar models, particularly in the presence of macroscopic mixing. This can be simulated by a local increase of opacity of a few percent. We conclude that the existing degeneracies and issues in solar modelling are not erased by an increase in the solar metallicity in contradiction to was suggested in recent papers. Therefore, standard solar models cannot be used as an argument for a high metallicity composition. While further work is required to improve solar models, we note that direct helioseismic inversions indicate a low metallicity in the convective envelope, in agreement with spectroscopic analyses based on full 3D models.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
Helioseismic determination of the solar metal mass fraction
Authors:
G. Buldgen,
A. Noels,
V. A. Baturin,
A. V. Oreshina,
S. V. Ayukov,
R. Scuflaire,
A. M. Amarsi,
N. Grevesse
Abstract:
Context. The metal mass fraction of the Sun Z is a key constraint in solar modelling, but its value is still under debate. The standard solar chemical composition of the late 2000s have the ratio of metals to hydrogen Z/X = 0.0181, with a small increase to 0.0187 in 2021, as inferred from 3D non-LTE spectroscopy. However, more recent work on a horizontally and temporally averaged <3D> model claim…
▽ More
Context. The metal mass fraction of the Sun Z is a key constraint in solar modelling, but its value is still under debate. The standard solar chemical composition of the late 2000s have the ratio of metals to hydrogen Z/X = 0.0181, with a small increase to 0.0187 in 2021, as inferred from 3D non-LTE spectroscopy. However, more recent work on a horizontally and temporally averaged <3D> model claim Z/X = 0.0225, consistent with the high values of twenty-five years ago based on 1D LTE spectroscopy. Aims. We aim to determine a precise and robust value of the solar metal mass fraction from helioseismic inversions, thus providing independent constraints from spectroscopic methods. Methods. We devise a detailed seismic reconstruction technique of the solar envelope, combining multiple inversions and equations of state to accurately and precisely determine the metal mass fraction value. Results. We show that a low value of the solar metal mass fraction corresponding to Z/X = 0.0187 is favoured by helioseismic constraints and that a higher metal mass fraction corresponding to Z/X = 0.0225 are strongly rejected by helioseismic data. Conclusions. We conclude that direct measurement of the metal mass fraction in the solar envelope favours a low metallicity, in line with the 3D non-LTE spectroscopic determination of 2021. A high metal mass fraction as measured using a <3D> model in 2022 is disfavoured by helioseismology for all modern equations of state used to model the solar convective envelope.
△ Less
Submitted 25 August, 2023;
originally announced August 2023.
-
Higher metal abundances do not solve the solar problem
Authors:
G. Buldgen,
P. Eggenberger,
A. Noels,
R. Scuflaire,
A. M. Amarsi,
N. Grevesse,
S. Salmon
Abstract:
Context. The Sun acts as a cornerstone of stellar physics. Thanks to spectroscopic, helioseismic and neutrino flux observations, we can use the Sun as a laboratory of fundamental physics in extreme conditions. The conclusions we draw are then used to inform and calibrate evolutionary models of all other stars in the Universe. However, solar models are in tension with helioseismic constraints. The…
▽ More
Context. The Sun acts as a cornerstone of stellar physics. Thanks to spectroscopic, helioseismic and neutrino flux observations, we can use the Sun as a laboratory of fundamental physics in extreme conditions. The conclusions we draw are then used to inform and calibrate evolutionary models of all other stars in the Universe. However, solar models are in tension with helioseismic constraints. The debate on the ``solar problem'' has hitherto led to numerous publications discussing potential issues with solar models and abundances. Aims. Using the recently suggested high-metallicity abundances for the Sun, we investigate whether standard solar models, as well as models with macroscopic transport reproducing the solar surface lithium abundances and analyze their properties in terms of helioseismic and neutrino flux observations. Methods. We compute solar evolutionary models and combine spectroscopic and helioseismic constraints as well as neutrino fluxes to investigate the impact of macroscopic transport on these measurements. Results. When high-metallicity solar models are calibrated to reproduce the measured solar lithium depletion, tensions arise with respect to helioseismology and neutrino fluxes. This is yet another demonstration that the solar problem is also linked to the physical prescriptions of solar evolutionary models and not to chemical composition alone. Conclusions. A revision of the physical ingredients of solar models is needed in order to improve our understanding of stellar structure and evolution. The solar problem is not limited to the photospheric abundances if the depletion of light elements is considered. In addition, tighter constraints on the solar beryllium abundance will play a key role in the improvement of solar models.
△ Less
Submitted 13 December, 2022;
originally announced December 2022.
-
Asteroseismology of evolved stars with EGGMiMoSA I. Theoretical mixed-mode patterns from the subgiant to the RGB phase
Authors:
M. Farnir,
C. Pinçon,
M-A. Dupret,
A. Noels,
R. Scuflaire
Abstract:
This study is the first of a series of papers that provide a technique to analyse the mixed-modes frequency spectra and characterise the structure of stars on the subgiant and red-giant branches. We define seismic indicators, relevant of the stellar structure and study their evolution on a grid of models. The proposed method, EGGMiMoSA, relies on the asymptotic description of mixed modes, defines…
▽ More
This study is the first of a series of papers that provide a technique to analyse the mixed-modes frequency spectra and characterise the structure of stars on the subgiant and red-giant branches. We define seismic indicators, relevant of the stellar structure and study their evolution on a grid of models. The proposed method, EGGMiMoSA, relies on the asymptotic description of mixed modes, defines initial guesses for the parameters, and uses a Levenberg-Marquardt technique to adjust the mixed-modes pattern efficiently. We follow the evolution of the mixed-modes parameters along a grid of models from the subgiant phase to the RGB bump and extend past works. We show the impact of the mass and composition on their evolution. The evolution of the period spacing $Δπ_1$, pressure offset $ε_p$, gravity offset $ε_g$, and coupling factor $q$ as a function of $Δν$ is little affected by the chemical composition and it follows two different regimes depending on the evolutionary stage. On the subgiant branch, the models display a moderate core-envelope density contrast. The evolution of $Δπ_1$, $ε_p$, $ε_g$, and $q$ thus significantly changes with the mass. Also, we demonstrate that, at fixed Z/X and with proper measurements of $Δπ_1$ and $Δν$, we may unambiguously constrain the mass, radius and age of a subgiant star. Conversely, on the red-giant branch, the core-envelope density contrast becomes very large. Consequently, the evolution of $ε_p$, $ε_g$ and $q$ as a function of $Δν$ becomes independent of the mass. This is also true for $Δπ_1$ in stars with masses $\lesssim 1.8M_\odot$ because of core electron degeneracy. This degeneracy is lifted for higher masses, again allowing for a precise measurement of the age. Overall, our computations qualitatively agree with past observed and theoretical studies.
△ Less
Submitted 16 July, 2021;
originally announced July 2021.
-
Standard solar models: a perspective from updated solar neutrino fluxes and the gravity-mode period spacing
Authors:
Sébastien Salmon,
Gaël Buldgen,
Arlette Noels,
Patrick Eggenberger,
Richard Scuflaire,
Georges Meynet
Abstract:
Context: The Sun is by far a privileged target for testing stellar models with unique precision. A recent concern appeared with the progress in the solar surface abundances derivation that has led to a decrease of the solar metallicity. While the ancient high-metallicity models were in fair agreement with other solar observational indicators, it is no longer the case for low-metallicity models. Re…
▽ More
Context: The Sun is by far a privileged target for testing stellar models with unique precision. A recent concern appeared with the progress in the solar surface abundances derivation that has led to a decrease of the solar metallicity. While the ancient high-metallicity models were in fair agreement with other solar observational indicators, it is no longer the case for low-metallicity models. Recent collection of data are however promising to shed a new light on it. For instance, the Borexino collaboration released in 2020 the first-ever complete estimate of neutrinos emitted in the CNO cycle. It has reaffirmed the role of the neutrino constraints in the solar modelling process and its associated issues. In parallel, newly claimed detection of solar gravity modes of oscillations offers another opportunity of probing the stratification in the Sun's central layers. Aims: We propose to combine the diagnoses from neutrinos and helioseismology, both from pressure and gravity modes, for assessing the predictions of solar models. We compare in detail the different physical prescriptions currently at disposal for stellar model computations. Results: The CNO neutrino flux confirms a preference for high-metallicity models. Nevertheless, we found that mild modification of the nuclear screening factors can re-match low-metallicity model predictions to observed fluxes, although it does not restore the agreement with the helioseismic frequency ratios. Neither the high-metallicity or low-metallicity models are able to reproduce the gravity-mode period spacing. The disagreement is huge, more than 100$σ$ to the reported value. Reversely, the family of standard models narrows the expected range of the Sun's period spacing: between $\sim$2150 to $\sim$2190~s. Moreover, we show this indicator can constrain the chemical mixture, opacity, and to a lower extent nuclear reactions in solar models.
△ Less
Submitted 3 May, 2021;
originally announced May 2021.
-
On the impact of the structural surface effect on global stellar properties and asteroseismic analyses
Authors:
Andreas Christ Sølvsten Jørgensen,
Josefina Montalbán,
George C. Angelou,
Andrea Miglio,
Achim Weiss,
Richard Scuflaire,
Arlette Noels,
Jakob Rørsted Mosumgaard,
Víctor Silva Aguirre
Abstract:
In a series of papers, we have recently demonstrated that it is possible to construct stellar structure models that robustly mimic the stratification of multi-dimensional radiative magneto-hydrodynamic simulations at every time-step of the computed evolution. The resulting models offer a more realistic depiction of the near-surface layers of stars with convective envelopes than parameterizations,…
▽ More
In a series of papers, we have recently demonstrated that it is possible to construct stellar structure models that robustly mimic the stratification of multi-dimensional radiative magneto-hydrodynamic simulations at every time-step of the computed evolution. The resulting models offer a more realistic depiction of the near-surface layers of stars with convective envelopes than parameterizations, such as mixing length theory, do. In this paper, we explore how this model improvement impacts on seismic and non-seismic properties of stellar models across the Hertzsprung-Russell diagram. We show that the improved description of the outer boundary layers alters the predicted global stellar properties at different evolutionary stages. In a hare and hound exercise, we show that this plays a key role for asteroseismic analyses, as it, for instance, often shifts the inferred stellar age estimates by more than 10 per cent. Improper boundary conditions may thus introduce systematic errors that exceed the required accuracy of the PLATO space mission. Moreover, we discuss different approximations for how to compute stellar oscillation frequencies. We demonstrate that the so-called gas $Γ_1$ approximation performs reasonably well for all main-sequence stars. Using a Monte Carlo approach, we show that the model frequencies of our hybrid solar models are consistent with observations within the uncertainties of the global solar parameters when using the so-called reduced $Γ_1$ approximation.
△ Less
Submitted 23 September, 2020;
originally announced September 2020.
-
Seismic Solar Models from Ledoux discriminant inversions
Authors:
G. Buldgen,
P. Eggenberger,
V. A. Baturin,
T. Corbard,
J. Christensen-Dalsgaard,
S. J. A. J. Salmon,
A. Noels,
A. V. Oreshina,
R. Scuflaire
Abstract:
The Sun constitutes an excellent laboratory of fundamental physics. With the advent of helioseismology, we were able to probe its internal layers with unprecedented precision. However, the current state of solar modelling is still stained by tedious issues. One of these problems is related to the disagreement between models computed with recent photospheric abundances and helioseismic constraints.…
▽ More
The Sun constitutes an excellent laboratory of fundamental physics. With the advent of helioseismology, we were able to probe its internal layers with unprecedented precision. However, the current state of solar modelling is still stained by tedious issues. One of these problems is related to the disagreement between models computed with recent photospheric abundances and helioseismic constraints. We use solar evolutionary models as initial conditions for reintegrations of their structure using Ledoux discriminant inversions. The resulting models are defined as seismic solar models, satisfying the equations of hydrostatic equilibrium. They will allow us to better constrain the internal structure of the Sun and provide complementary information to that of evolutionary models. These seismic models were computed using various reference models with different equations of state, abundances and opacity tables. We check the robustness of our approach by confirming the good agreement of our seismic models in terms of sound speed, density and entropy proxy inversions as well as frequency-separation ratios of low-degree pressure modes. Our method allows us to determine with an excellent accuracy the Ledoux discriminant profile of the Sun and compute full profiles of this quantity. Our models show an agreement with seismic data of ~0.1% in sound speed, density and entropy proxy as well as with the observed frequency-separation ratios. They surpass all standard and non-standard evolutionary models including ad-hoc changes aiming at reproducing helioseismic constraints. The obtained seismic Ledoux discriminant profile as well as the consistent structure obtained from our procedure paves the way for renewed attempts at constraining the solar modelling problem and the missing physical processes acting in the solar interior by breaking free from the hypotheses of evolutionary models.
△ Less
Submitted 20 July, 2020;
originally announced July 2020.
-
Chronologically dating the early assembly of the Milky Way
Authors:
Josefina Montalbán,
John Ted Mackereth,
Andrea Miglio,
Fiorenzo Vincenzo,
Cristina Chiappini,
Gaël Buldgen,
Benoît Mosser,
Arlette Noels,
Richard Scuflaire,
Mathieu Vrard,
Emma Willett,
Guy R. Davies,
Oliver Hall,
Martin Bo Nielsen,
Saniya Khan,
Ben M. Rendle,
Walter E. van Rossem,
Jason W. Ferguson,
William J. Chaplin
Abstract:
The standard cosmological model ($Λ$-CDM) predicts that galaxies are built through hierarchical assembly on cosmological timescales$^{1,2}$. The Milky Way, like other disc galaxies, underwent violent mergers and accretion of small satellite galaxies in its early history. Thanks to Gaia-DR2$^3$ and spectroscopic surveys$^4$, the stellar remnants of such mergers have been identified$^{5-7}$. The chr…
▽ More
The standard cosmological model ($Λ$-CDM) predicts that galaxies are built through hierarchical assembly on cosmological timescales$^{1,2}$. The Milky Way, like other disc galaxies, underwent violent mergers and accretion of small satellite galaxies in its early history. Thanks to Gaia-DR2$^3$ and spectroscopic surveys$^4$, the stellar remnants of such mergers have been identified$^{5-7}$. The chronological dating of such events is crucial to uncover the formation and evolution of the Galaxy at high redshift, but it has so far been challenging owing to difficulties in obtaining precise ages for these oldest stars. Here we combine asteroseismology -- the study of stellar oscillations -- with kinematics and chemical abundances, to estimate precise stellar ages ($\sim$ 11%) for a sample of stars observed by the $\mathit{Kepler}$ space mission$^8$. Crucially, this sample includes not only some of the oldest stars that were formed inside the Galaxy, but also stars formed externally and subsequently accreted onto the Milky Way. Leveraging this resolution in age, we provide compelling evidence in favour of models in which the Galaxy had already formed a substantial population of its stars (which now reside mainly in its thick disc) before the in-fall of the satellite galaxy Gaia-Enceladus/Sausage$^{5,6}$ around 10 billions years ago
△ Less
Submitted 17 May, 2021; v1 submitted 2 June, 2020;
originally announced June 2020.
-
Investigating surface correction relations for RGB stars
Authors:
Andreas Christ Sølvsten Jørgensen,
Josefina Montalbán,
Andrea Miglio,
Ben M. Rendle,
Guy R. Davies,
Gaël Buldgen,
Richard Scuflaire,
Arlette Noels,
Patrick Gaulme,
Rafael A. García
Abstract:
State-of-the-art stellar structure and evolution codes fail to adequately describe turbulent convection. For stars with convective envelopes, such as red giants, this leads to an incomplete depiction of the surface layers. As a result, the predicted stellar oscillation frequencies are haunted by systematic errors, the so-called surface effect. Different empirically and theoretically motivated corr…
▽ More
State-of-the-art stellar structure and evolution codes fail to adequately describe turbulent convection. For stars with convective envelopes, such as red giants, this leads to an incomplete depiction of the surface layers. As a result, the predicted stellar oscillation frequencies are haunted by systematic errors, the so-called surface effect. Different empirically and theoretically motivated correction relations have been proposed to deal with this issue. In this paper, we compare the performance of these surface correction relations for red giant branch stars. For this purpose, we apply the different surface correction relations in asteroseismic analyses of eclipsing binaries and open clusters. In accordance with previous studies of main-sequence stars, we find that the use of different surface correction relations biases the derived global stellar properties, including stellar age, mass, and distance estimates. We furthermore demonstrate that the different relations lead to the same systematic errors for two different open clusters. Our results overall discourage from the use of surface correction relations that rely on reference stars to calibrate free parameters. Due to the demonstrated systematic biasing of the results, the use of appropriate surface correction relations is imperative to any asteroseismic analysis of red giants. Accurate mass, age, and distance estimates for red giants are fundamental when addressing questions that deal with the chemo-dynamical evolution of the Milky Way galaxy. In this way, our results also have implications for fields, such as galactic archaeology, that draw on findings from stellar physics.
△ Less
Submitted 6 June, 2020; v1 submitted 28 April, 2020;
originally announced April 2020.
-
Evolutionary models for ultracool dwarfs
Authors:
Catarina S. Fernandes,
Valerie Van Grootel,
Sebastien J. A. J. Salmon,
Bernhard Aringer,
Adam J. Burgasser,
Richard Scuflaire,
Pierre Brassard,
Gilles Fontaine
Abstract:
Ultracool dwarfs have emerged as key targets for searches of transiting exoplanets. Precise estimates of the host parameters (including mass, age, and radius) are fundamental to constrain the physical properties of orbiting exoplanets. We have extended our evolutionary code CLES (Code Liégeois d'Evolution Stellaire) to the ultracool dwarf regime. We include relevant equations of state for H, He, a…
▽ More
Ultracool dwarfs have emerged as key targets for searches of transiting exoplanets. Precise estimates of the host parameters (including mass, age, and radius) are fundamental to constrain the physical properties of orbiting exoplanets. We have extended our evolutionary code CLES (Code Liégeois d'Evolution Stellaire) to the ultracool dwarf regime. We include relevant equations of state for H, He, as well as C and O elements to cover the temperature-density regime of ultracool dwarf interiors. For various metallicities, we couple the interior models to two sets of model atmospheres as surface boundary conditions. We show that including C and O in the EOS has a significant effect close the H-burning limit mass. The typical systematic error associated with uncertainties in input physics in evolutionary models is $\sim 0.0005 M_\odot$. We test model results against observations for objects whose parameters have been determined from independent techniques. We are able to reproduce dynamical mass measurements of LSPM J1314+1320AB within $1σ$ with the condition of varying the metallicity (determined from calibrations) up to $2.5σ$. For GJ 65AB, a $2σ$ agreement is obtained between individual masses from differential astrometry and those from evolutionary models. We provide tables of ultracool dwarf models for various masses and metallicities that can be used as reference when estimating parameters for ultracool objects.
△ Less
Submitted 13 June, 2019;
originally announced June 2019.
-
AIMS - A new tool for stellar parameter determinations using asteroseismic constraints
Authors:
B. M. Rendle,
G. Buldgen,
A. Miglio,
D. Reese,
A. Noels,
G. R. Davies,
T. L. Campante,
W. J. Chaplin,
M. N. Lund,
J. S. Kuszlewicz,
L. J. A. Scott,
R. Scuflaire,
W. H. Ball,
J. Smetana,
B. Nsamba
Abstract:
A key aspect in the determination of stellar properties is the comparison of observational constraints with predictions from stellar models. Asteroseismic Inference on a Massive Scale (AIMS) is an open source code that uses Bayesian statistics and a Markov Chain Monte Carlo approach to find a representative set of models that reproduce a given set of classical and asteroseismic constraints. These…
▽ More
A key aspect in the determination of stellar properties is the comparison of observational constraints with predictions from stellar models. Asteroseismic Inference on a Massive Scale (AIMS) is an open source code that uses Bayesian statistics and a Markov Chain Monte Carlo approach to find a representative set of models that reproduce a given set of classical and asteroseismic constraints. These models are obtained by interpolation on a pre-calculated grid, thereby increasing computational efficiency. We test the accuracy of the different operational modes within AIMS for grids of stellar models computed with the Liège stellar evolution code (main sequence and red giants) and compare the results to those from another asteroseismic analysis pipeline, PARAM. Moreover, using artificial inputs generated from models within the grid (assuming the models to be correct), we focus on the impact on the precision of the code when considering different combinations of observational constraints (individual mode frequencies, period spacings, parallaxes, photospheric constraints,...). Our tests show the absolute limitations of precision on parameter inferences using synthetic data with AIMS, and the consistency of the code with expected parameter uncertainty distributions. Interpolation testing highlights the significance of the underlying physics to the analysis performance of AIMS and provides caution as to the upper limits in parameter step size. All tests demonstrate the flexibility and capability of AIMS as an analysis tool and its potential to perform accurate ensemble analysis with current and future asteroseismic data yields.
△ Less
Submitted 9 January, 2019;
originally announced January 2019.
-
Combining multiple structural inversions to constrain the Solar modelling problem
Authors:
G. Buldgen,
S. J. A. J. Salmon,
A. Noels,
R. Scuflaire,
J. Montalban,
V. A. Baturin,
P. Eggenberger,
V. K. Gryaznov,
I. L. Iosilevskiy,
G. Meynet,
W. J. Chaplin,
A. Miglio,
A. V. Oreshina,
O. Richard,
A. N. Starostin
Abstract:
The Sun is the most studied of all stars. It is a reference for all other observed stars and a laboratory of fundamental physics helping us understand processes occuring in conditions irreproducible on Earth. However, our understanding of the Sun is currently stained by the solar modelling problem which can stem from various causes, such as the opacities, the equation of state and the mixing of ch…
▽ More
The Sun is the most studied of all stars. It is a reference for all other observed stars and a laboratory of fundamental physics helping us understand processes occuring in conditions irreproducible on Earth. However, our understanding of the Sun is currently stained by the solar modelling problem which can stem from various causes, such as the opacities, the equation of state and the mixing of chemical elements. We combine inversions of sound speed, an entropy proxy and the Ledoux discriminant with constraints such as the position of the base of the convective zone and the photospheric helium abundance. We test various combinations of standard ingredients for solar modelling such as abundance tables, equation of state, formalism for convection and diffusion and opacity tables and study the diagnostic potential of the inversions to constrain ad-hoc modifications of the opacity profile and additional mixing below the convective envelope. Combining inversions provides stringent constraints on the modifications on the models, far beyond what is achievable only from sound speed inversions. We constrain the form and amplitude of the opacity increase required and show that a 15% increase at log T=6.35 provides a significant improvement but that a more global increase of the opacity, within the uncertainties of current tables, coupled with an additional mixing at the bottom of the convective zone gives the best agreement for low metallicity models. We show that high metallicity models do not satisfy all the inversion results. We conclude that the solar problem likely occurs from various small sources, as ingredients such as the equation of state or the formalism of convection can induce small but significant changes and that using phase shift analyses combined with our approach is the next step for a better understanding of the inaccuracies of solar models just below the convective envelope.
△ Less
Submitted 24 September, 2018;
originally announced September 2018.
-
Mean density inversions for red giants and red clump stars
Authors:
Gaël Buldgen,
B. Rendle,
T. Sonoi,
G. R. Davies,
A. Miglio,
S. J. A. J. Salmon,
D. R. Reese,
D. Bossini,
P. Eggenberger,
A. Noels,
R. Scuflaire
Abstract:
Since the CoRoT and Kepler missions, the availability of high quality seismic spectra for red giants has made them the standard clocks and rulers for Galactic Archeology. With the expected excellent data from the TESS and PLATO missions, red giants will again play a key role in Galactic studies and stellar physics, thanks to the precise masses and radii determined by asteroseismology. The determin…
▽ More
Since the CoRoT and Kepler missions, the availability of high quality seismic spectra for red giants has made them the standard clocks and rulers for Galactic Archeology. With the expected excellent data from the TESS and PLATO missions, red giants will again play a key role in Galactic studies and stellar physics, thanks to the precise masses and radii determined by asteroseismology. The determination of these quantities is often based on so-called scaling laws, which have been used extensively for main-sequence stars. We show how the SOLA inversion technique can provide robust determinations of the mean density of red giants within 1 per cent of the real value, using only radial oscillations. Combined with radii determinations from Gaia of around 2 per cent precision, this approach provides robust, less model-dependent masses with an error lower than 10 per cent. It will improve age determinations, helping to accurately dissect the Galactic structure and history. We present results on artificial data of standard models, models including an extended atmosphere from averaged 3D simulations and non-adiabatic frequency calculations to test surface effects, and on eclipsing binaries. We show that the inversions provide very robust mean density estimates, using at best seismic information. However, we also show that a distinction between red-giant branch and red-clump stars is required to determine a reliable estimate of the mean density. The stability of the inversion enables an implementation in automated pipelines, making it suitable for large samples of stars.
△ Less
Submitted 25 August, 2018;
originally announced August 2018.
-
Stellar parameters for TRAPPIST-1
Authors:
Valerie Van Grootel,
Catarina S. Fernandes,
Michaël Gillon,
Emmanuel Jehin,
Jean Manfroid,
Richard Scuflaire,
Adam J. Burgasser,
Artem Burdanov,
Laetitia Delrez,
Brice-Olivier Demory,
Julien de Wit,
Didier Queloz,
Amaury H. M. J. Triaud
Abstract:
TRAPPIST-1 is an ultracool dwarf star transited by seven Earth-sized planets, for which thorough characterization of atmospheric properties, surface conditions encompassing habitability and internal compositions is possible with current and next generation telescopes. Accurate modeling of the star is essential to achieve this goal. We aim to obtain updated stellar parameters for TRAPPIST- 1 based…
▽ More
TRAPPIST-1 is an ultracool dwarf star transited by seven Earth-sized planets, for which thorough characterization of atmospheric properties, surface conditions encompassing habitability and internal compositions is possible with current and next generation telescopes. Accurate modeling of the star is essential to achieve this goal. We aim to obtain updated stellar parameters for TRAPPIST- 1 based on new measurements and evolutionary models, compared to those used in discovery studies. We present a new measurement for the parallax of TRAPPIST-1, 82.4 $\pm$ 0.8 mas, based on 188 epochs of observations with the TRAPPIST and Liverpool Telescopes from 2013 to 2016. This revised parallax yields an updated luminosity of $L_*=(5.22\pm0.19)\times 10^{-4} L_{\odot}$, very close to the previous estimate but almost twice more precise. We next present an updated estimate for TRAPPIST-1 stellar mass, based on two approaches: mass from stellar evolution modeling, and empirical mass derived from dynamical masses of equivalently classified ultracool dwarfs in astrometric binaries. We combine them through a Monte-Carlo approach to derive a semi-empirical estimate for the mass of TRAPPIST-1. We also derive estimate for the radius by combining this mass with stellar density inferred from transits, as well as estimate for the effective temperature from our revised luminosity and radius. Our final results are $M_*=0.089 \pm 0.006 M_{\odot}$, $R_* = 0.121 \pm 0.003 R_{\odot}$, and $T_{\rm eff} =$ 2516 $\pm$ 41 K. Considering the degree to which TRAPPIST-1 system will be scrutinized in coming years, these revised and more precise stellar parameters should be considered when assessing the properties of TRAPPIST-1 planets.
△ Less
Submitted 5 December, 2017;
originally announced December 2017.
-
Determining the metallicity of the solar envelope using seismic inversion techniques
Authors:
Gaël Buldgen,
S. J. A. J. Salmon,
A. Noels,
R. Scuflaire,
M. A. Dupret,
D. R. Reese
Abstract:
The solar metallicity issue is a long-lasting problem of astrophysics, impacting multi- ple fields and still subject to debate and uncertainties. While spectroscopy has mostly been used to determine the solar heavy elements abundance, helioseismologists at- tempted providing a seismic determination of the metallicity in the solar convective enveloppe. However, the puzzle remains since two independ…
▽ More
The solar metallicity issue is a long-lasting problem of astrophysics, impacting multi- ple fields and still subject to debate and uncertainties. While spectroscopy has mostly been used to determine the solar heavy elements abundance, helioseismologists at- tempted providing a seismic determination of the metallicity in the solar convective enveloppe. However, the puzzle remains since two independent groups prodived two radically different values for this crucial astrophysical parameter. We aim at provid- ing an independent seismic measurement of the solar metallicity in the convective enveloppe. Our main goal is to help provide new information to break the current stalemate amongst seismic determinations of the solar heavy element abundance. We start by presenting the kernels, the inversion technique and the target function of the inversion we have developed. We then test our approach in multiple hare-and-hounds exercises to assess its reliability and accuracy. We then apply our technique to solar data using calibrated solar models and determine an interval of seismic measurements for the solar metallicity. We show that our inversion can indeed be used to estimate the solar metallicity thanks to our hare-and-hounds exercises. However, we also show that further dependencies in the physical ingredients of solar models lead to a low accuracy. Nevertheless, using various physical ingredients for our solar models, we determine metallicity values between 0.008 and 0.014.
△ Less
Submitted 1 September, 2017;
originally announced September 2017.
-
Inversions of the Ledoux discriminant: a closer look at the tachocline
Authors:
Gaël Buldgen,
S. J. A. J. Salmon,
M. Godart,
A. Noels,
R. Scuflaire,
M. A. Dupret,
D. R. Reese,
J. Colgan,
C. J. Fontes,
P. Eggenberger,
P. Hakel,
D. P. Kilcrease,
O. Richard
Abstract:
Modelling the base of the solar convective envelope is a tedious problem. Since the first rotation inversions, solar modellers are confronted with the fact that a region of very limited extent has an enormous physical impact on the Sun. Indeed, it is the transition region from differential to solid body rotation, the tachocline, which furthermore is influenced by turbulence and is also supposed to…
▽ More
Modelling the base of the solar convective envelope is a tedious problem. Since the first rotation inversions, solar modellers are confronted with the fact that a region of very limited extent has an enormous physical impact on the Sun. Indeed, it is the transition region from differential to solid body rotation, the tachocline, which furthermore is influenced by turbulence and is also supposed to be the seat of the solar magnetic dynamo. Moreover, solar models show significant disagreement with the sound speed profile in this region. In this paper, we show how helioseismology can provide further constraints on this region by carrying out an inversion of the Ledoux discriminant. We compare these inversions for Standard Solar Models built using various opacity tables and chemical abundances and discuss the origins of the discrepancies between Solar Models and the Sun.
△ Less
Submitted 1 September, 2017;
originally announced September 2017.
-
Seismic inversion of the solar entropy: A case for improving the Standard Solar Model
Authors:
G. Buldgen,
S. J. A. J. Salmon,
A. Noels,
R. Scuflaire,
D. R. Reese,
M-A. Dupret,
J. Colgan,
C. J. Fontes,
P. Eggenberger,
P. Hakel,
D. P. Kilcrease,
S. Turck-Chièze
Abstract:
The Sun is the most constrained and well-studied of all stars. As a consequence, the physical ingredients entering solar models are used as a reference to study all other stars observed in the Universe. However, our understanding of the solar structure is still imperfect, as illustrated by the current debate on the heavy element abundances in the Sun. We wish to provide additional information on t…
▽ More
The Sun is the most constrained and well-studied of all stars. As a consequence, the physical ingredients entering solar models are used as a reference to study all other stars observed in the Universe. However, our understanding of the solar structure is still imperfect, as illustrated by the current debate on the heavy element abundances in the Sun. We wish to provide additional information on the solar structure by carrying out structural inversions of a new physical quantity, a proxy of the entropy of the solar plasma which properties are very sensitive to the temperature gradient below the convective zone. We use new structural kernels to carry out direct inversions of an entropy proxy of the solar plasma and compare the solar structure to various standard solar models built using various opacity tables and chemical abundances. We also link our results to classical tests commonly found in the literature. Our analysis allows us to probe more efficiently the uncertain regions of the solar models, just below the convective zone, paving the way for new in-depth analyses of the Sun taking into account additional physical uncertainties of solar models beyond the specific question of chemical abundances.
△ Less
Submitted 17 July, 2017;
originally announced July 2017.
-
Prospects for asteroseismic inference on the envelope helium abundance in red giant stars
Authors:
A. -M. Broomhall,
A. Miglio,
J. Montalban,
P. Eggenberger,
W. J. Chaplin,
Y. Elsworth,
R. Scuflaire,
P. Ventura,
G. A. Verner
Abstract:
Regions of rapid variation in the internal structure of a star are often referred to as acoustic glitches since they create a characteristic periodic signature in the frequencies of p modes. Here we examine the localized disturbance arising from the helium second ionization zone in red giant branch and clump stars. More specifically, we determine how accurately and precisely the parameters of the…
▽ More
Regions of rapid variation in the internal structure of a star are often referred to as acoustic glitches since they create a characteristic periodic signature in the frequencies of p modes. Here we examine the localized disturbance arising from the helium second ionization zone in red giant branch and clump stars. More specifically, we determine how accurately and precisely the parameters of the ionization zone can be obtained from the oscillation frequencies of stellar models. We use models produced by three different generation codes that not only cover a wide range of stages of evolution along the red giant phase but also incorporate different initial helium abundances. We discuss the conditions under which such fits robustly and accurately determine the acoustic radius of the second ionization zone of helium. The determined radii of the ionization zones as inferred from the mode frequencies were found to be coincident with the local maximum in the first adiabatic exponent described by the models, which is associated with the outer edge of the second ionization zone of helium. Finally, we consider whether this method can be used to distinguish stars with different helium abundances. Although a definite trend in the amplitude of the signal is observed any distinction would be difficult unless the stars come from populations with vastly different helium abundances or the uncertainties associated with the fitted parameters can be reduced. However, application of our methodology could be useful for distinguishing between different populations of red giant stars in globular clusters, where distinct populations with very different helium abundances have been observed.
△ Less
Submitted 27 March, 2014;
originally announced March 2014.
-
Transit confirmation and improved stellar and planet parameters for the super-Earth HD 97658 b and its host star
Authors:
V. Van Grootel,
M. Gillon,
D. Valencia,
N. Madhusudhan,
D. Dragomir,
A. R. Howe,
A. S. Burrows,
B. -O. Demory,
D. Deming,
D. Ehrenreich,
C. Lovis,
M. Mayor,
F. Pepe,
D. Queloz,
R. Scuflaire,
S. Seager,
D. Segransan,
S. Udry
Abstract:
Super-Earths transiting nearby bright stars are key objects that simultaneously allow for accurate measurements of both their mass and radius, providing essential constraints on their internal composition. We present here the confirmation, based on Spitzer transit observations, that the super-Earth HD 97658 b transits its host star. HD 97658 is a low-mass ($M_*=0.77\pm0.05\,M_{\odot}$) K1 dwarf, a…
▽ More
Super-Earths transiting nearby bright stars are key objects that simultaneously allow for accurate measurements of both their mass and radius, providing essential constraints on their internal composition. We present here the confirmation, based on Spitzer transit observations, that the super-Earth HD 97658 b transits its host star. HD 97658 is a low-mass ($M_*=0.77\pm0.05\,M_{\odot}$) K1 dwarf, as determined from the Hipparcos parallax and stellar evolution modeling. To constrain the planet parameters, we carry out Bayesian global analyses of Keck-HIRES radial velocities, and MOST and Spitzer photometry. HD 97658 b is a massive ($M_P=7.55^{+0.83}_{-0.79} M_{\oplus}$) and large ($R_{P} = 2.247^{+0.098}_{-0.095} R_{\oplus}$ at 4.5 $μ$m) super-Earth. We investigate the possible internal compositions for HD 97658 b. Our results indicate a large rocky component, by at least 60% by mass, and very little H-He components, at most 2% by mass. We also discuss how future asteroseismic observations can improve the knowledge of the HD 97658 system, in particular by constraining its age. Orbiting a bright host star, HD 97658 b will be a key target for coming space missions TESS, CHEOPS, PLATO, and also JWST, to characterize thoroughly its structure and atmosphere.
△ Less
Submitted 26 February, 2014; v1 submitted 24 February, 2014;
originally announced February 2014.
-
Testing convective-core overshooting using period spacings of dipole modes in red giants
Authors:
J. Montalban,
A. Miglio,
A. Noels,
M. -A. Dupret,
R. Scuflaire,
P. Ventura
Abstract:
Uncertainties on central mixing in main sequence (MS) and core He-burning (He-B) phases affect key predictions of stellar evolution such as late evolutionary phases, chemical enrichment, ages etc. We propose a test of the extension of extra-mixing in two relevant evolutionary phases based on period spacing Delta_P of solar-like oscillating giants. From stellar models and their corresponding adiaba…
▽ More
Uncertainties on central mixing in main sequence (MS) and core He-burning (He-B) phases affect key predictions of stellar evolution such as late evolutionary phases, chemical enrichment, ages etc. We propose a test of the extension of extra-mixing in two relevant evolutionary phases based on period spacing Delta_P of solar-like oscillating giants. From stellar models and their corresponding adiabatic frequencies (respectively computed with ATON and LOSC codes) we provide the first predictions of the observable Delta_P for stars in the red giant branch (RGB) and in the red clump (RC). We find: i) a clear correlation between Delta_P and the mass of the helium core (M_He); the latter in intermediate-mass stars depends on the MS overshooting, hence it can be used to set constraints on extra mixing during MS when coupled with chemical composition; ii) a linear dependence of the average value of the asymptotic period spacing (<Delta_P>_a) during the He-B phase on the size of the convective core. A first comparison with the inferred asymptotic period spacing for Kepler RC stars suggests the need for extra mixing also during this phase, as evinced from other observational facts.
△ Less
Submitted 13 February, 2013;
originally announced February 2013.
-
Successful Asteroseismology for a Better Characterization of the Exoplanet HAT-P-7b
Authors:
M. Oshagh,
A. Grigahcène,
O. Benomar,
M. -A. Dupret,
M. J. P. F. G. Monteiro,
R. Scuflaire,
N. C. Santos
Abstract:
It is well known that asteroseismology is the unique technique permitting the study of the internal structure of pulsating stars using their pulsational frequencies, which is per se very important. It acquires an additional value when the star turns out to be a planet host. In this case, the asteroseismic study output may be a very important input for the study of the planetary system. With this i…
▽ More
It is well known that asteroseismology is the unique technique permitting the study of the internal structure of pulsating stars using their pulsational frequencies, which is per se very important. It acquires an additional value when the star turns out to be a planet host. In this case, the asteroseismic study output may be a very important input for the study of the planetary system. With this in mind, we use the large time-span of the Kepler public data obtained for the star system HAT-P-7, first to perform an asteroseismic study of the pulsating star using Time-Dependent-Convection (TDC) models. Secondly, we make a revision of the planet properties in the light of the asteroseismic study.
△ Less
Submitted 2 January, 2013;
originally announced January 2013.
-
Towards precise asteroseismology of solar-like stars
Authors:
A. Grigahcène,
M. -A. Dupret,
S. G. Sousa,
M. J. P. F. G. Monteiro,
R. Garrido,
R. Scuflaire,
M. Gabriel
Abstract:
Adiabatic modeling of solar-like oscillations cannot exceed a certain level of precision for fitting individual frequencies. This is known as the problem of near surface effects on the mode physics. We present a theoretical study which addresses the problem of frequency precision in non-adiabatic models using a time-dependent convection treatment. We find that the number of acceptable model soluti…
▽ More
Adiabatic modeling of solar-like oscillations cannot exceed a certain level of precision for fitting individual frequencies. This is known as the problem of near surface effects on the mode physics. We present a theoretical study which addresses the problem of frequency precision in non-adiabatic models using a time-dependent convection treatment. We find that the number of acceptable model solutions is significantly reduced and more precise constraints can be imposed on the models. Results obtained for a specific star (beta Hydri) lead to very good agreement with both global and local seismic observables. This indicates that the accuracy of model fitting to seismic data is greatly improved when a more complete description of the interaction between convection and pulsation is taken into account.
△ Less
Submitted 27 December, 2011;
originally announced December 2011.
-
Adiabatic solar-like oscillations in red giant stars
Authors:
J. Montalban,
A. Miglio,
A. Noels,
R. Scuflaire,
P. Ventura,
F. D'Antona
Abstract:
Since the detection of non-radial solar-like oscillation modes in red giants with the CoRoT satellite, the interest in the asteroseismic properties of red giants and the link with their global properties and internal structure is substantially increasing. Moreover, more and more precise data are being collected with the space-based telescopes CoRoT and Kepler. In this paper we present a survey of…
▽ More
Since the detection of non-radial solar-like oscillation modes in red giants with the CoRoT satellite, the interest in the asteroseismic properties of red giants and the link with their global properties and internal structure is substantially increasing. Moreover, more and more precise data are being collected with the space-based telescopes CoRoT and Kepler. In this paper we present a survey of the most relevant theoretical and observational results obtained up to now concerning the potential of solar-like oscillations in red giants.
△ Less
Submitted 24 August, 2011;
originally announced August 2011.
-
Inference from adiabatic analysis of solar-like oscillations in Red giants
Authors:
Josefina Montalban,
Andrea Miglio,
Arlette Noels,
Richard Scuflaire,
Paolo Ventura
Abstract:
The clear detection with CoRoT and KEPLER of radial and non-radial solar-like oscillations in many red giants paves the way to seismic inferences on the structure of such stars. We present an overview of the properties of the adiabatic frequencies and frequency separations of radial and non-radial oscillation modes, highlighting how their detection allows a deeper insight into the properties of th…
▽ More
The clear detection with CoRoT and KEPLER of radial and non-radial solar-like oscillations in many red giants paves the way to seismic inferences on the structure of such stars. We present an overview of the properties of the adiabatic frequencies and frequency separations of radial and non-radial oscillation modes, highlighting how their detection allows a deeper insight into the properties of the internal structure of red giants. In our study we consider models of red giants in different evolutionary stages, as well as of different masses and chemical composition. We describe how the large and small separations computed with radial modes and with non-radial modes mostly trapped in the envelope depend on the stellar global parameters and evolutionary state, and we compare our theoretical predictions and first KEPLER data.Finally, we find that the properties of dipole modes constitute a promising seismic diagnostic of the evolutionary state of red-giant stars.
△ Less
Submitted 20 December, 2010;
originally announced December 2010.
-
Seismic diagnostics of red giants: first comparison with stellar models
Authors:
J. Montalbán,
A. Miglio,
A. Noels,
R. Scuflaire,
P. Ventura
Abstract:
The clear detection with CoRoT and KEPLER of radial and non-radial solar-like oscillations in many red giants paves the way to seismic inferences on the structure of such stars. We present an overview of the properties of the adiabatic frequencies and frequency separations of radial and non-radial oscillation modes for an extended grid of models. We highlight how their detection allows a deeper in…
▽ More
The clear detection with CoRoT and KEPLER of radial and non-radial solar-like oscillations in many red giants paves the way to seismic inferences on the structure of such stars. We present an overview of the properties of the adiabatic frequencies and frequency separations of radial and non-radial oscillation modes for an extended grid of models. We highlight how their detection allows a deeper insight into the internal structure and evolutionary state of red giants. In particular, we find that the properties of dipole modes constitute a promising seismic diagnostic tool of the evolutionary state of red-giant stars. We compare our theoretical predictions with the first 34 days of KEPLER data and predict the frequency diagram expected for red giants in the COROT exofield in the galactic center direction.
△ Less
Submitted 9 September, 2010;
originally announced September 2010.
-
Evidence for a sharp structure variation inside a red-giant star
Authors:
A. Miglio,
J. Montalbán,
F. Carrier,
J. De Ridder,
B. Mosser,
P. Eggenberger,
R. Scuflaire,
P. Ventura,
F. D'Antona,
A. Noels,
A. Baglin
Abstract:
The availability of precisely determined frequencies of radial and non-radial oscillation modes in red giants is finally paving the way for detailed studies of the internal structure of these stars. We look for the seismic signature of regions of sharp structure variation in the internal structure of the CoRoT target HR7349. We analyse the frequency dependence of the large frequency separation and…
▽ More
The availability of precisely determined frequencies of radial and non-radial oscillation modes in red giants is finally paving the way for detailed studies of the internal structure of these stars. We look for the seismic signature of regions of sharp structure variation in the internal structure of the CoRoT target HR7349. We analyse the frequency dependence of the large frequency separation and second frequency differences, as well as the behaviour of the large frequency separation obtained with the envelope auto-correlation function. We find evidence for a periodic component in the oscillation frequencies, i.e. the seismic signature of a sharp structure variation in HR7349. In a comparison with stellar models we interpret this feature as caused by a local depression of the sound speed that occurs in the helium second-ionization region. Using solely seismic constraints this allows us to estimate the mass (M=1.2^{+0.6}_{-0.4} Msun) and radius (R=12.2^{+2.1}_{-1.8} Rsun) of HR7349, which agrees with the location of the star in an HR diagram.
△ Less
Submitted 6 September, 2010;
originally announced September 2010.
-
Solar-like oscillations in a massive star
Authors:
K. Belkacem,
R. Samadi,
M. -J. Goupil,
L. Lefevre,
F. Baudin,
S. Deheuvels,
M. -A. Dupret,
T. Appourchaux,
R. Scuflaire,
M. Auvergne,
C. Catala,
E. Michel,
A. Miglio,
J. Montalban,
A. Thoul,
S. Talon,
A. Baglin,
A. Noels
Abstract:
Seismology of stars provides insight into the physical mechanisms taking place in their interior, with modes of oscillation probing different layers. Low-amplitude acoustic oscillations excited by turbulent convection were detected four decades ago in the Sun and more recently in low-mass main-sequence stars. Using data gathered by the Convection Rotation and Planetary Transits mission, we repor…
▽ More
Seismology of stars provides insight into the physical mechanisms taking place in their interior, with modes of oscillation probing different layers. Low-amplitude acoustic oscillations excited by turbulent convection were detected four decades ago in the Sun and more recently in low-mass main-sequence stars. Using data gathered by the Convection Rotation and Planetary Transits mission, we report here on the detection of solar-like oscillations in a massive star, V1449 Aql, which is a known large-amplitude (b Cephei) pulsator.
△ Less
Submitted 20 June, 2009;
originally announced June 2009.
-
A seismic approach to testing different formation channels of subdwarf B stars
Authors:
Haili Hu,
M. -A. Dupret,
C. Aerts,
G. Nelemans,
S. D. Kawaler,
A. Miglio,
J. Montalban,
R. Scuflaire
Abstract:
There are many unknowns in the formation of subdwarf B stars. Different formation channels are considered to be possible and to lead to a variety of helium-burning subdwarfs. All seismic models to date, however, assume that a subdwarf B star is a post-helium-flash-core surrounded by a thin inert layer of hydrogen. We examine an alternative formation channel, in which the subdwarf B star originat…
▽ More
There are many unknowns in the formation of subdwarf B stars. Different formation channels are considered to be possible and to lead to a variety of helium-burning subdwarfs. All seismic models to date, however, assume that a subdwarf B star is a post-helium-flash-core surrounded by a thin inert layer of hydrogen. We examine an alternative formation channel, in which the subdwarf B star originates from a massive (>~2 Msun) red giant with a non-degenerate helium-core. Although these subdwarfs may evolve through the same region of the log g-Teff diagram as the canonical post-flash subdwarfs, their interior structure is rather different. We examine how this difference affects their pulsation modes and whether it can be observed.
Using detailed stellar evolution calculations we construct subdwarf B models from both formation channels. The iron accumulation in the driving region due to diffusion, which causes the excitation of the modes, is approximated by a Gaussian function. The pulsation modes and frequencies are calculated with a non-adiabatic pulsation code. A detailed comparison of two subdwarf B models from different channels, but with the same log g and Teff, shows that their mode excitation is different. The excited frequencies are lower for the post-flash than for the post-non-degenerate subdwarf B star. This is mainly due to the differing chemical composition of the stellar envelope. A more general comparison between two grids of models shows that the excited frequencies of most post-non-degenerate subdwarfs cannot be well-matched with the frequencies of post-flash subdwarfs. In the rare event that an acceptable seismic match is found, additional information, such as mode identification and log g and Teff determinations, allows us to distinguish between the two formation channels.
△ Less
Submitted 22 August, 2008;
originally announced August 2008.
-
The CoRoT Evolution and Seismic Tools Activity: Goals and Tasks
Authors:
Y. Lebreton,
M. J. P. F. G. Monteiro,
J. Montalban,
A. Moya,
A. Baglin,
J. Christensen-Dalsgaard,
M. -J. Goupil,
E. Michel,
J. Provost,
I. W. Roxburgh,
R. Scuflaire,
the ESTA Team
Abstract:
The forthcoming data expected from space missions such as CoRoT require the capacity of the available tools to provide accurate models whose numerical precision is well above the expected observational errors. In order to secure that these tools meet the specifications, a team has been established to test and, when necessary, to improve the codes available in the community. The CoRoT evolution a…
▽ More
The forthcoming data expected from space missions such as CoRoT require the capacity of the available tools to provide accurate models whose numerical precision is well above the expected observational errors. In order to secure that these tools meet the specifications, a team has been established to test and, when necessary, to improve the codes available in the community. The CoRoT evolution and seismic tool activity (ESTA) has been set up with this mission.
Several groups have been involved. The present paper describes the motivation and the organisation of this activity, providing the context and the basis for the presentation of the results that have been achieved so far. This is not a finished task as future even better data will continue to demand more precise and complete tools for asteroseismology.
△ Less
Submitted 7 May, 2008;
originally announced May 2008.
-
Thorough analysis of input physics in CESAM and CLES codes
Authors:
Josefina Montalban,
Yveline Lebreton,
Andrea Miglio,
Richard Scuflaire,
Pierre Morel,
Arlette Noels
Abstract:
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 e…
▽ More
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.
△ Less
Submitted 14 April, 2008;
originally announced April 2008.
-
Grids of Stellar Models and Frequencies with CLES + LOSC
Authors:
J. Montalban,
A. Miglio,
A. Noels,
R. Scuflaire
Abstract:
We present a grid of stellar models, obtained with the CLES evolution code, following the specification of ESTA-Task1, and the corresponfing seismic properties, computed with the LOSC code. We provide a complete description of the corresponding files that will be available on the ESTA web-pages.
We present a grid of stellar models, obtained with the CLES evolution code, following the specification of ESTA-Task1, and the corresponfing seismic properties, computed with the LOSC code. We provide a complete description of the corresponding files that will be available on the ESTA web-pages.
△ Less
Submitted 22 February, 2008;
originally announced February 2008.
-
The Liege Oscillation Code
Authors:
R. Scuflaire,
J. Montalban,
S. Theado,
P. -O. Bourge,
A. Miglio,
M. Godart,
A. Thoul,
A. Noels
Abstract:
The Liege Oscillation code can be used as a stand-alone program or as a library of subroutines that the user calls from a Fortran main program of his own to compute radial and non-radial adiabatic oscillations of stellar models. We describe the variables and the equations used by the program and the methods used to solve them. A brief account is given of the use and the output of the program.
The Liege Oscillation code can be used as a stand-alone program or as a library of subroutines that the user calls from a Fortran main program of his own to compute radial and non-radial adiabatic oscillations of stellar models. We describe the variables and the equations used by the program and the methods used to solve them. A brief account is given of the use and the output of the program.
△ Less
Submitted 20 December, 2007;
originally announced December 2007.
-
CLES, Code Liegeois d'Evolution Stellaire
Authors:
R. Scuflaire,
S. Theado,
J. Montalban,
A. Miglio,
P. -O. Bourge,
M. Godart,
A. Thoul,
A. Noels
Abstract:
Cles is an evolution code recently developed to produce stellar models meeting the specific requirements of studies in asteroseismology. It offers the users a lot of choices in the input physics they want in their models and its versatility allows them to tailor the code to their needs and implement easily new features. We describe the features implemented in the current version of the code and…
▽ More
Cles is an evolution code recently developed to produce stellar models meeting the specific requirements of studies in asteroseismology. It offers the users a lot of choices in the input physics they want in their models and its versatility allows them to tailor the code to their needs and implement easily new features. We describe the features implemented in the current version of the code and the techniques used to solve the equations of stellar structure and evolution. A brief account is given of the use of the program and of a solar calibration realized with it.
△ Less
Submitted 20 December, 2007;
originally announced December 2007.
-
Inter-comparison of the g-, f- and p-modes calculated using different oscillation codes for a given stellar model
Authors:
A. Moya,
J. Christensen-Dalsgaard,
S. Charpinet,
Y. Lebreton,
A. Miglio,
J. Montalban,
M. J. P. F. G. Monteiro,
J. Provost,
I. W. Roxburgh,
R. Scuflaire,
J. C. Suarez,
M. Suran
Abstract:
In order to make astroseismology a powerful tool to explore stellar interiors, different numerical codes should give the same oscillation frequencies for the same input physics. This work is devoted to test, compare and, if needed, optimize the seismic codes used to calculate the eigenfrequencies to be finally compared with observations. The oscillation codes of nine research groups in the field…
▽ More
In order to make astroseismology a powerful tool to explore stellar interiors, different numerical codes should give the same oscillation frequencies for the same input physics. This work is devoted to test, compare and, if needed, optimize the seismic codes used to calculate the eigenfrequencies to be finally compared with observations. The oscillation codes of nine research groups in the field have been used in this study. The same physics has been imposed for all the codes in order to isolate the non-physical dependence of any possible difference. Two equilibrium models with different grids, 2172 and 4042 mesh points, have been used, and the latter model includes an explicit modelling of semiconvection just outside the convective core. Comparing the results for these two models illustrates the effect of the number of mesh points and their distribution in particularly critical parts of the model, such as the steep composition gradient outside the convective core. A comprehensive study of the frequency differences found for the different codes is given as well. These differences are mainly due to the use of different numerical integration schemes. The use of a second-order integration scheme plus a Richardson extrapolation provides similar results to a fourth-order integration scheme. The proper numerical description of the Brunt-Vaisala frequency in the equilibrium model is also critical for some modes. An unexpected result of this study is the high sensitivity of the frequency differences to the inconsistent use of values of the gravitational constant (G) in the oscillation codes, within the range of the experimentally determined ones, which differ from the value used to compute the equilibrium model.
△ Less
Submitted 16 November, 2007;
originally announced November 2007.
-
An asteroseismic study of the beta Cephei star theta Ophiuchi: constraints on global stellar parameters and core overshooting
Authors:
M. Briquet,
T. Morel,
A. Thoul,
R. Scuflaire,
A. Miglio,
J. Montalban,
M. -A. Dupret,
C. Aerts
Abstract:
We present a seismic study of the beta Cephei star theta Ophiuchi. Our analysis is based on the observation of one radial mode, one rotationally split l = 1 triplet and three components of a rotationally split l = 2 quintuplet for which the m-values were well identified by spectroscopy. We identify the radial mode as fundamental, the triplet as p_1 and the quintuplet as g_1. Our NLTE abundance a…
▽ More
We present a seismic study of the beta Cephei star theta Ophiuchi. Our analysis is based on the observation of one radial mode, one rotationally split l = 1 triplet and three components of a rotationally split l = 2 quintuplet for which the m-values were well identified by spectroscopy. We identify the radial mode as fundamental, the triplet as p_1 and the quintuplet as g_1. Our NLTE abundance analysis results in a metallicity and CNO abundances in full agreement with the most recent updated solar values. With X \in [0.71,0.7211] and Z \in [0.009,0.015], and using the Asplund et al. (2005) mixture but with a Ne abundance about 0.3 dex larger (Cunha et al. 2006), the matching of the three independent modes, enables us to deduce constrained ranges for the mass (M = 8.2 +/- 0.3 Msun) and central hydrogen abundance (X_c = 0.38 +/- 0.02) of theta Oph and to prove the occurrence of core overshooting (alpha_ov = 0.44 +/- 0.07). We also derive an equatorial rotation velocity of 29 +/- 7 km/s. Moreover, we show that the observed non-equidistance of the l=1 triplet can be reproduced by second order effects of rotation. Finally, we show that the observed rotational splitting of two modes cannot rule out a rigid rotation model.
△ Less
Submitted 22 June, 2007;
originally announced June 2007.
-
Report on the CoRoT Evolution and Seismic Tools Activity
Authors:
M. K. J. P. F. G. Monteiro,
Y. Lebreton,
J. Montalban,
J. Christensen-Dalsgaard,
M. Castro,
S. Degl'Innocenti,
A. Moya,
I. W. Roxburgh,
R. Scuflaire,
A. Baglin,
M. S. Cunha,
P. Eggenberger,
J. Fernandes,
M. J. Goupil,
A. Hui-Bon-Hoa,
M. Marconi,
J. P. Marques,
E. Michel,
A. Miglio,
P. Morel,
B. Pichon,
P. G. Prada Moroni,
J. Provost,
A. Ruoppo,
J. -C. Suarez
, et al. (2 additional authors not shown)
Abstract:
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, com…
▽ More
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.
△ Less
Submitted 27 May, 2006;
originally announced May 2006.
-
Analysis of MERCATOR data Part I: variable B stars
Authors:
P. De Cat,
M. Briquet,
C. Aerts,
K. Goossens,
S. Saesen,
J. Cuypers,
K. Yakut,
R. Scuflaire,
M. -A. Dupret,
many observers
Abstract:
We re-classified 31 variable B stars which were observed more than 50 times in the Geneva photometric system with the P7 photometer attached to the MERCATOR telescope (La Palma) during its first 3 years of scientific observations. HD89688 is a possible beta Cephei/slowly pulsating B star hybrid and the main mode of the COROT target HD180642 shows non-linear effects. The Maia candidates are re-cl…
▽ More
We re-classified 31 variable B stars which were observed more than 50 times in the Geneva photometric system with the P7 photometer attached to the MERCATOR telescope (La Palma) during its first 3 years of scientific observations. HD89688 is a possible beta Cephei/slowly pulsating B star hybrid and the main mode of the COROT target HD180642 shows non-linear effects. The Maia candidates are re-classified as either ellipsoidal variables or spotted stars. Although the mode identification is still ongoing, all the well-identified modes so far have a degree l = 0, 1 or 2.
△ Less
Submitted 8 November, 2005;
originally announced November 2005.
-
A new seismic analysis of Alpha Centauri
Authors:
A. Thoul,
R. Scuflaire,
A. Noels,
B. Vatovez,
M. Briquet,
M. -A. Dupret,
J. Montalban
Abstract:
Models of alpha Cen A & B have been computed using the masses determined by Pourbaix et al. (2002) and the data derived from the spectroscopic analysis of Neuforge and Magain (1997). The seismological data obtained by Bouchy and Carrier (2001, 2002) do help improve our knowledge of the evolutionary status of the system. All the constraints are satisfied with a model which gives an age of about 6…
▽ More
Models of alpha Cen A & B have been computed using the masses determined by Pourbaix et al. (2002) and the data derived from the spectroscopic analysis of Neuforge and Magain (1997). The seismological data obtained by Bouchy and Carrier (2001, 2002) do help improve our knowledge of the evolutionary status of the system. All the constraints are satisfied with a model which gives an age of about 6 Gyr for the binary.
△ Less
Submitted 20 March, 2003;
originally announced March 2003.
-
A photometric mode identification method, including an improved non-adiabatic treatment of the atmosphere
Authors:
M. -A. Dupret,
J. De Ridder,
P. De Cat,
C. Aerts,
R. Scuflaire,
A. Noels,
A. Thoul
Abstract:
We present an improved version of the method of photometric mode identification of Heynderickx et al. (1994). Our new version is based on the inclusion of precise non-adiabatic eigenfunctions determined in the outer stellar atmosphere according to the formalism recently proposed by Dupret et al.(2002). Our improved photometric mode identification technique is therefore no longer dependent on ad…
▽ More
We present an improved version of the method of photometric mode identification of Heynderickx et al. (1994). Our new version is based on the inclusion of precise non-adiabatic eigenfunctions determined in the outer stellar atmosphere according to the formalism recently proposed by Dupret et al.(2002). Our improved photometric mode identification technique is therefore no longer dependent on ad hoc parameters for the non-adiabatic effects. It contains the complete physical conditions of the outer atmosphere of the star, provided that rotation does not play a key role. We apply our improved method to the two slowly pulsating B stars HD 74560 and HD 138764 and to the beta Cephei star EN (16) Lac. Besides identifying the degree l of the pulsating stars, our method is also a tool for improving the knowledge of stellar interiors and atmospheres, by imposing constraints on parameters such as the metallicity and the mixing-length parameter alpha (a procedure we label non-adiabatic asteroseismology).
△ Less
Submitted 18 November, 2002;
originally announced November 2002.