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