Site Map
Follow us on Facebook Follow us on Twitter YouTube channel
Centro de Astrofísica da Universidade do Porto

Inter-comparison of the g-, f- and p-modes calculated using different oscillation codes for a given stellar model

A. Moya Bedón, J. Christensen-Dalsgaard, S. Charpinet, Y. Lebreton, A. Miglio, J. Montalbán, M. J. P. F. G. Monteiro, J. Provost, I. W. Roxburgh, R. Scuflaire, J.-C. Suárez, M. D. Suran


In order to make asteroseismology a powerful tool to explore stellar interiors, different numerical codes should give the same oscillation frequencies for the same input physics. Any differences found when comparing the numerical values of the eigenfrequencies will be an important piece of information regarding the numerical structure of the code. The ESTA group was created to analyze the non-physical sources of these differences. The work presented in this report is a part of Task 2 of the ESTA group. Basically the 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 first step in this comparison is presented here. 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 possi ble 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 number of mesh points and their distribution are crucial for interpreting the results. 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. This influence depends on the set of the eigenfunctions used for the solution of the differential equations. 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. This effect can provide differences for a given equilibrium model substantially larger than those resulting from the use of different codes or numerical techniques; the actual differences between the values of G used by the different codes account for much of the frequency differences found here.

Astrophysics and Space Science
Volume 316, Page 231
April 2008

>> PDF>> ADS>> DOI

Institute of Astrophysics and Space Sciences

Institute of Astrophysics and Space Sciences (IA) is a new but long anticipated research infrastructure with a national dimension. It embodies a bold but feasible vision for the development of Astronomy, Astrophysics and Space Sciences in Portugal, taking full advantage and fully realizing the potential created by the national membership of the European Space Agency (ESA) and the European Southern Observatory (ESO). IA resulted from the merging the two most prominent research units in the field in Portugal: the Centre for Astrophysics of the University of Porto (CAUP) and the Center for Astronomy and Astrophysics of the University of Lisbon (CAAUL). It currently hosts more than two-thirds of all active researchers working in Space Sciences in Portugal, and is responsible for an even greater fraction of the national productivity in international ISI journals in the area of Space Sciences. This is the scientific area with the highest relative impact factor (1.65 times above the international average) and the field with the highest average number of citations per article for Portugal.

Proceed on CAUP's website|Go to IA website