We present the identification of very low frequency g modes in the asymptotic regime and two important parameters that have long been waited for: the core rotation rate, and the asymptotic equidistant period spacing of these g modes. The GOLF instrument on board the SOHO space observatory has provided two decades of full-disk helioseismic data. In the present study, we search for possible collective frequency modulations that are produced by periodic changes in the deep solar structure. Such modulations provide access to only very low frequency g modes, thus allowing statistical methods to take advantage of their asymptotic properties. For oscillatory periods in the range between 9 and nearly 48 hours, almost 100 g modes of spherical harmonic degree 1 and more than 100 g modes of degree 2 are predicted. They are not observed individually, but when combined, they unambiguouslyprovide their asymptotic period equidistance and rotational splittings, in excellent agreement with the requirements of the asymptotic approximations. Previously, p-mode helioseismology allowed the g-mode period equidistance parameter $P_0$ to be bracketed inside a narrow range, between approximately 34 and 35 minutes. Here, $P_0$ is measured to be 34 min 01 s, with a 1 s uncertainty. The previously unknown g-mode splittings have now been measured from a non-synodic reference with very high accuracy, and they imply a mean weighted rotation of 1277 $\pm$ 10 nHz (9-day period) of their kernels, resulting in a rapid rotation frequency of 1644 $\pm$ 23 nHz (period of one week) of the solar core itself, which is a factor 3.8 $\pm$ 0.1 faster than the rotation of the radiative envelope. The g modes are known to be the keys to a better understanding of the structure and dynamics of the solar core. Their detection with these precise parameters will certainly stimulate a new era of research in this field.

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.

With the space-borne missions CoRoT and Kepler, a large amount of asteroseismic data is now available. So-called global oscillation parameters are inferred to characterize the large sets of stars, to perform ensemble asteroseismology, and to derive scaling relations. The mean large separation is such a key parameter. It is therefore crucial to measure it with the highest accuracy. As the conditions of measurement of the large separation do not coincide with its theoretical definition, we revisit the asymptotic expressions used for analysing the observed oscillation spectra. Then, we examine the consequence of the difference between the observed and asymptotic values of the mean large separation. The analysis is focused on radial modes. We use series of radial-mode frequencies to compare the asymptotic and observational values of the large separation. We propose a simple formulation to correct the observed value of the large separation and then derive its asymptotic counterpart. We prove that, apart from glitches due to stellar structure discontinuities, the asymptotic expansion is valid from main-sequence stars to red giants. Our model shows that the asymptotic offset is close to 1/4, as in the theoretical development. High-quality solar-like oscillation spectra derived from precise photometric measurements are definitely better described with the second-order asymptotic expansion. The second-order term is responsible for the curvature observed in the \'echelle diagrams used for analysing the oscillation spectra and this curvature is responsible for the difference between the observed and asymptotic values of the large separation. Taking it into account yields a revision of the scaling relations providing more accurate asteroseismic estimates of the stellar mass and radius.

We review some major open issues in the current modelling of low and intermediate mass, main sequence stars based on seismological studies. The solar case was discussed in a companion paper, here several issues specific to other stars than the Sun are illustrated with a few stars observed with CoRoT and expectations from Kepler data.

We have analyzed data from a multi-site campaign to observe oscillations in the F5 star Procyon. The data consist of high-precision velocities that we obtained over more than three weeks with eleven telescopes. A new method for adjusting the data weights allows us to suppress the sidelobes in the power spectrum. Stacking the power spectrum in a so-called echelle diagram reveals two clear ridges that we identify with even and odd values of the angular degree (l=0 and 2, and l=1 and 3, respectively). We interpret a strong, narrow peak at 446 muHz that lies close to the l=1 ridge as a mode with mixed character. We show that the frequencies of the ridge centroids and their separations are useful diagnostics for asteroseismology. In particular, variations in the large separation appear to indicate a glitch in the sound-speed profile at an acoustic depth of about 1000 s. We list frequencies for 55 modes extracted from the data spanning 20 radial orders, a range comparable to the best solar data, which will provide valuable constraints for theoretical models. A preliminary comparison with published models shows that the offset between observed and calculated frequencies for the radial modes is very different for Procyon than for the Sun and other cool stars. We find the mean lifetime of the modes in Procyon to be 1.29 +0.55/-0.49 days, which is significantly shorter than the 2-4 days seen in the Sun.

By applying the theory of slowly rotating stars to the Sun, the solar quadrupole and octopole moments J2 and J4 were computed using a solar model obtained from CESAM stellar evolution code (Morel, 1997) combined with a recent model of solar differential rotation deduced from helioseismology (Corbard et al., 2002). This model takes into account a near-surface radial gradient of rotation which was inferred and quantified from MDI f-mode observations by Corbard and Thompson (2002). The effect of this observational near-surface gradient on the theoretical values of the surface parameters J2, J4 is investigated. The results show that the octopole moment J4 is much more sensitive than the quadrupole moment J2 to the subsurface radial gradient of rotation.

Solar gravity modes (or g modes) -- oscillations of the solar interior for which buoyancy acts as the restoring force -- have the potential to provide unprecedented inference on the structure and dynamics of the solar core, inference that is not possible with the well observed acoustic modes (or p modes). The high amplitude of the g-mode eigenfunctions in the core and the evanesence of the modes in the convection zone make the modes particularly sensitive to the physical and dynamical conditions in the core. Owing to the existence of the convection zone, the g modes have very low amplitudes at photospheric levels, which makes the modes extremely hard to detect. In this paper, we review the current state of play regarding attempts to detect g modes. We review the theory of g modes, including theoretical estimation of the g-mode frequencies, amplitudes and damping rates. Then we go on to discuss the techniques that have been used to try to detect g modes. We review results in the literature, and finish by looking to the future, and the potential advances that can be made -- from both data and data-analysis perspectives -- to give unambiguous detections of individual g modes. The review ends by concluding that, at the time of writing, there is indeed a consensus amongst the authors that there is currently no undisputed detection of solar g modes.

We have carried out a multi-site campaign to measure oscillations in the F5 star Procyon A. We obtained high-precision velocity observations over more than three weeks with eleven telescopes, with almost continuous coverage for the central ten days. This represents the most extensive campaign so far organized on any solar-type oscillator. We describe in detail the methods we used for processing and combining the data. These involved calculating weights for the velocity time series from the measurement uncertainties and adjusting them in order to minimize the noise level of the combined data. The time series of velocities for Procyon shows the clear signature of oscillations, with a plateau of excess power that is centred at 0.9 mHz and is broader than has been seen for other stars. The mean amplitude of the radial modes is 38.1 +/- 1.3 cm/s (2.0 times solar), which is consistent with previous detections from the ground and by the WIRE spacecraft, and also with the upper limit set by the MOST spacecraft. The variation of the amplitude during the observing campaign allows us to estimate the mode lifetime to be 1.5 d (+1.9/-0.8 d). We also find a slow variation in the radial velocity of Procyon, with good agreement between different telescopes. These variations are remarkably similar to those seen in the Sun, and we interpret them as being due to rotational modulation from active regions on the stellar surface. The variations appear to have a period of about 10 days, which presumably equals the stellar rotation period or, perhaps, half of it. The amount of power in these slow variations indicates that the fractional area of Procyon covered by active regions is slightly higher than for the Sun.

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.

Since the beginning of this century we have attended a blooming of the gravity-mode research thanks to the unprecedented quality of the data available, either from space with SoHO, or from the ground-based networks as BiSON or GONG. From the first upper limit of the gravity-mode amplitudes fixed at 10 mm/s at 200 microHz given by Appourchaux et al. (2000), on one hand, a peak was supposed to be a component of the l=1, n=1 mixed mode (Garcia et al. 2001a, b; Gabriel et al. 2002) and, on the other hand, a couple of patterns --multiplets-- were attributed to gravity modes (Turck-Chieze et al. 2004; Mathur et al. 2007). One of these patterns, found around 220 microHz, could be labeled as the l=2, n =-3 g mode, which is expected to be the one with the highest surface amplitude (Cox and Guzik 2004). Finally, in 2007, Garcia et al. were able to measure the fingertips of the dipole gravity modes looking for their asymptotic properties. In the present paper we present an update of the recent developments on this subject with special attention to the 220 microHz region, the dipole asymptotic properties and the impact of the incoming g-mode observations on the knowledge of the solar structure and rotation profile.

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.

Procyon A is a bright F5IV star in a binary system. Although the distance, mass and angular diameter of this star are all known with high precision, the exact evolutionary state is still unclear. Evolutionary tracks with different ages and different mass fractions of hydrogen in the core pass, within the errors, through the observed position of Procyon A in the Hertzsprung-Russell diagram. For more than 15 years several different groups have studied the solar-like oscillations in Procyon A to determine its evolutionary state. Although several studies independently detected power excess in the periodogram, there is no agreement on the actual oscillation frequencies yet. This is probably due to either insufficient high-quality data (i.e., aliasing) or due to intrinsic properties of the star (i.e., short mode lifetimes). Now a spectroscopic multi-site campaign using 10 telescopes world-wide (minimizing aliasing effects) with a total time span of nearly 4 weeks (increase the frequency resolution) is performed to identify frequencies in this star and finally determine its properties and evolutionary state.

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%

In the present work we show robust indications of the existence of g modes in the Sun using 10 years of GOLF data. The present analysis is based on the exploitation of the collective properties of the predicted low-frequency (25 to 140 microHz) g modes: their asymptotic nature, which implies a quasi equidistant separation of their periods for a given angular degree (l). The Power Spectrum (PS) of the Power Spectrum Density (PSD), reveals a significant structure indicating the presence of features (peaks) in the PSD with near equidistant periods corresponding to l=1 modes in the range n=-4 to n=-26. The study of its statistical significance of this feature was fully undertaken and complemented with Monte Carlo simulations. This structure has a confidence level better than 99.86% not to be due to pure noise. Furthermore, a detailed study of this structure suggests that the gravity modes have a much more complex structure than the one initially expected (line-widths, magnetic splittings...). Compared to the latest solar models, the obtained results tend to favor a solar core rotating significantly faster than the rest of the radiative zone. In the framework of the Phoebus group, we have also applied the same methodology to other helioseismology instruments on board SoHO and ground based networks.

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.

The solar magnetism is no more considered as a purely superficial phenomenon. The SoHO community has shown that the length of the solar cycle depends on the transition region between radiation and convection. Nevertheless, the internal solar (stellar) magnetism stays poorly known. Starting in 2008, the American instrument HMI/SDO and the European microsatellite PICARD will enrich our view of the Sun-Earth relationship. Thus obtaining a complete MHD solar picture is a clear objective for the next decades and it requires complementary observations of the dynamics of the radiative zone. For that ambitious goal, space prototypes are being developed to improve gravity mode detection. The Sun is unique to progress on the topology of deep internal magnetic fields and to understand the complex mechanisms which provoke photospheric and coronal magnetic changes and possible longer cycles important for human life. We propose the following roadmap in Europe to contribute to this "impressive" revolution in Astronomy and in our Sun-Earth relationship: SoHO (1995-2007), PICARD (2008-2010), DynaMICS (2009-2017) in parallel to SDO (2008-2017) then a world-class mission located at the L1 orbit or above the solar poles.

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.

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.

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.

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.

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.

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.

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.

Recent studies have established that peaks in solar oscillation power spectra are not Lorentzian in shape, but have a distinct asymmetry. Fitting a symmetric Lorentzian profile to the peaks therefore produces a shift in frequency of the modes. Accurate determination of low-frequency modes is essential to infer the structure of the solar core by inversion of the mode frequencies. In this paper we investigate how the changes in frequencies of low-degree modes obtained by fitting symmetric and asymmetric peak profiles change the inferred properties of the solar core. We use data obtained by the Global Oscillations at Low Frequencies (GOLF) project on board the SoHO spacecraft. Two different solar models and inversion procedures are used to invert the data to determine the sound speed in the solar core. We find that for a given set of modes no significant difference in the inferred sound-speed results from taking asymmetry into account when fitting the low-degree modes.

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.

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}.

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.

Inversions of rotational splittings have shown that there exists at the base of the solar convection zone a region called the tachocline in which high radial gradients of the rotation rate occur. The usual linear regularization methods tend to smooth out any high gradients in the solution, and may not be appropriate for the study of this zone. In this paper we use, in the helioseismic context of rotation inversions, regularization methods that have been developed for edge-preserving regularization in computed imaging. It is shown from Monte-Carlo simulations that this approach can lead directly to results similar to those reached by linear inversions which however required some assumptions on the shape of the transition in order to be deconvolved. The application of this method to LOWL data leads to a very thin tachocline. From the discussions on the parameters entering the inversion and the Monte-Carlo simulations, our conclusion is that the tachocline width is very likely below 0.05R_sun which lowers our previous estimate of 0.05+/- 0.03R_sun obtained from the same dataset (Corbard et al. 1998).

Inversions of rotational splittings have shown that the surface layers and the so-called solar tachocline at the base of the convection zone are regions in which high radial gradients of the rotation rate occur. The usual regularization methods tend to smooth out every high gradients in the solution and may not be appropriate for the study of a zone like the tachocline. In this paper we use nonlinear regularization methods that are developed for edge-preserving regularization in computed imaging (e.g. Blanc-F\'eraud et al. 1995) and we apply them in the helioseismic context of rotational inversions.

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)

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.

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.