A Bayesian approach to the modelling of α Cen A

Michael Bazot, S. Bourguignon, J. Christensen-Dalsgaard

Research output: Contribution to journalArticle

Abstract

Determining the physical characteristics of a star is an inverse problem consisting of estimating the parameters of models for the stellar structure and evolution, and knowing certain observable quantities. We use a Bayesian approach to solve this problem for α Cen A, which allows us to incorporate prior information on the parameters to be estimated, in order to better constrain the problem. Our strategy is based on the use of a Markov chain Monte Carlo (MCMC) algorithm to estimate the posterior probability densities of the stellar parameters: mass, age, initial chemical composition, etc. We use the stellar evolutionary code astec to model the star. To constrain this model both seismic and non-seismic observations were considered. Several different strategies were tested to fit these values, using either two free parameters or five free parameters in astec. We are thus able to show evidence that MCMC methods become efficient with respect to more classical grid-based strategies when the number of parameters increases. The results of our MCMC algorithm allow us to derive estimates for the stellar parameters and robust uncertainties thanks to the statistical analysis of the posterior probability densities. We are also able to compute odds for the presence of a convective core in α Cen A. When using core-sensitive seismic observational constraints, these can rise above ∼40 per cent. The comparison of results to previous studies also indicates that these seismic constraints are of critical importance for our knowledge of the structure of this star.

Original languageEnglish (US)
Pages (from-to)1847-1866
Number of pages20
JournalMonthly Notices of the Royal Astronomical Society
Volume427
Issue number3
DOIs
StatePublished - Dec 11 2012

Fingerprint

Markov chains
modeling
Markov chain
stars
stellar structure
stellar evolution
estimates
inverse problem
parameter
statistical analysis
Monte Carlo method
chemical composition
estimating
grids

Keywords

  • Methods: numerical
  • Methods: statistical
  • Stars: evolution
  • Stars: fundamental parameters
  • Stars: individual: α Cen A
  • Stars: oscillations

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

A Bayesian approach to the modelling of α Cen A. / Bazot, Michael; Bourguignon, S.; Christensen-Dalsgaard, J.

In: Monthly Notices of the Royal Astronomical Society, Vol. 427, No. 3, 11.12.2012, p. 1847-1866.

Research output: Contribution to journalArticle

Bazot, Michael ; Bourguignon, S. ; Christensen-Dalsgaard, J. / A Bayesian approach to the modelling of α Cen A. In: Monthly Notices of the Royal Astronomical Society. 2012 ; Vol. 427, No. 3. pp. 1847-1866.
@article{e26e26fdb63f4f889a1e1d041ae92bb0,
title = "A Bayesian approach to the modelling of α Cen A",
abstract = "Determining the physical characteristics of a star is an inverse problem consisting of estimating the parameters of models for the stellar structure and evolution, and knowing certain observable quantities. We use a Bayesian approach to solve this problem for α Cen A, which allows us to incorporate prior information on the parameters to be estimated, in order to better constrain the problem. Our strategy is based on the use of a Markov chain Monte Carlo (MCMC) algorithm to estimate the posterior probability densities of the stellar parameters: mass, age, initial chemical composition, etc. We use the stellar evolutionary code astec to model the star. To constrain this model both seismic and non-seismic observations were considered. Several different strategies were tested to fit these values, using either two free parameters or five free parameters in astec. We are thus able to show evidence that MCMC methods become efficient with respect to more classical grid-based strategies when the number of parameters increases. The results of our MCMC algorithm allow us to derive estimates for the stellar parameters and robust uncertainties thanks to the statistical analysis of the posterior probability densities. We are also able to compute odds for the presence of a convective core in α Cen A. When using core-sensitive seismic observational constraints, these can rise above ∼40 per cent. The comparison of results to previous studies also indicates that these seismic constraints are of critical importance for our knowledge of the structure of this star.",
keywords = "Methods: numerical, Methods: statistical, Stars: evolution, Stars: fundamental parameters, Stars: individual: α Cen A, Stars: oscillations",
author = "Michael Bazot and S. Bourguignon and J. Christensen-Dalsgaard",
year = "2012",
month = "12",
day = "11",
doi = "10.1111/j.1365-2966.2012.21818.x",
language = "English (US)",
volume = "427",
pages = "1847--1866",
journal = "Monthly Notices of the Royal Astronomical Society",
issn = "0035-8711",
publisher = "Oxford University Press",
number = "3",

}

TY - JOUR

T1 - A Bayesian approach to the modelling of α Cen A

AU - Bazot, Michael

AU - Bourguignon, S.

AU - Christensen-Dalsgaard, J.

PY - 2012/12/11

Y1 - 2012/12/11

N2 - Determining the physical characteristics of a star is an inverse problem consisting of estimating the parameters of models for the stellar structure and evolution, and knowing certain observable quantities. We use a Bayesian approach to solve this problem for α Cen A, which allows us to incorporate prior information on the parameters to be estimated, in order to better constrain the problem. Our strategy is based on the use of a Markov chain Monte Carlo (MCMC) algorithm to estimate the posterior probability densities of the stellar parameters: mass, age, initial chemical composition, etc. We use the stellar evolutionary code astec to model the star. To constrain this model both seismic and non-seismic observations were considered. Several different strategies were tested to fit these values, using either two free parameters or five free parameters in astec. We are thus able to show evidence that MCMC methods become efficient with respect to more classical grid-based strategies when the number of parameters increases. The results of our MCMC algorithm allow us to derive estimates for the stellar parameters and robust uncertainties thanks to the statistical analysis of the posterior probability densities. We are also able to compute odds for the presence of a convective core in α Cen A. When using core-sensitive seismic observational constraints, these can rise above ∼40 per cent. The comparison of results to previous studies also indicates that these seismic constraints are of critical importance for our knowledge of the structure of this star.

AB - Determining the physical characteristics of a star is an inverse problem consisting of estimating the parameters of models for the stellar structure and evolution, and knowing certain observable quantities. We use a Bayesian approach to solve this problem for α Cen A, which allows us to incorporate prior information on the parameters to be estimated, in order to better constrain the problem. Our strategy is based on the use of a Markov chain Monte Carlo (MCMC) algorithm to estimate the posterior probability densities of the stellar parameters: mass, age, initial chemical composition, etc. We use the stellar evolutionary code astec to model the star. To constrain this model both seismic and non-seismic observations were considered. Several different strategies were tested to fit these values, using either two free parameters or five free parameters in astec. We are thus able to show evidence that MCMC methods become efficient with respect to more classical grid-based strategies when the number of parameters increases. The results of our MCMC algorithm allow us to derive estimates for the stellar parameters and robust uncertainties thanks to the statistical analysis of the posterior probability densities. We are also able to compute odds for the presence of a convective core in α Cen A. When using core-sensitive seismic observational constraints, these can rise above ∼40 per cent. The comparison of results to previous studies also indicates that these seismic constraints are of critical importance for our knowledge of the structure of this star.

KW - Methods: numerical

KW - Methods: statistical

KW - Stars: evolution

KW - Stars: fundamental parameters

KW - Stars: individual: α Cen A

KW - Stars: oscillations

UR - http://www.scopus.com/inward/record.url?scp=84869756372&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84869756372&partnerID=8YFLogxK

U2 - 10.1111/j.1365-2966.2012.21818.x

DO - 10.1111/j.1365-2966.2012.21818.x

M3 - Article

AN - SCOPUS:84869756372

VL - 427

SP - 1847

EP - 1866

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

SN - 0035-8711

IS - 3

ER -