Issue 
A&A
Volume 604, August 2017



Article Number  A40  
Number of page(s)  17  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201730460  
Published online  01 August 2017 
Asymptotic g modes: Evidence for a rapid rotation of the solar core
^{1} Université Côte d’Azur, Observatoire Côte d’Azur, CNRS, Laboratoire Lagrange, CS 34229, Nice Cedex 4, France
email: eric.fossat@oca.eu
^{2} Institut d’Astrophysique Spatiale, Université ParisSud and CNRS (UMR 8617), Bâtiment 121, 91405 Orsay Cedex, France
^{3} Laboratoire AIM ParisSaclay, CEA/DRFCNRSUniv. Paris Diderot – IRFU/SAp, Centre de Saclay, 91191 GifsurYvette, France
^{4} LAB, 2 rue de l’Observatoire, BP89, 33271 Floirac Cedex, France
^{5} Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Spain
^{6} Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
^{7} DAP/IRFU/CEA, UMR AIM, University ParisSaclay, CE Saclay, 91191 GifsurYvette, France
^{8} Department of Physics and Astronomy Department, University of California at Los Angeles, 430 Portola Plaza, Box 951547, Los Angeles, CA 900951547, USA
^{9} LPHEA Laboratory, Oukaimeden Observatory, Cadi Ayyad University FSSM, BP 2390 Marrakech, Morocco
Received: 18 January 2017
Accepted: 24 May 2017
Context. Over the past 40 years, helioseismology has been enormously successful in the study of the solar interior. A shortcoming has been the lack of a convincing detection of the solar g modes, which are oscillations driven by gravity and are hidden in the deepest part of the solar body – its hydrogenburning core. The detection of g modes is expected to dramatically improve our ability to model this core, the rotational characteristics of which have, until now, remained unknown.
Aims. 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.
Methods. The GOLF instrument on board the SOHO space observatory has provided two decades of fulldisk helioseismic data. The search for g modes in GOLF measurements has been extremely difficult because of solar and instrumental noise. In the present study, the p modes of the GOLF signal are analyzed differently: 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.
Results. For oscillatory periods in the range between 9 and nearly 48 h, 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 unambiguously provide their asymptotic period equidistance and rotational splittings, in excellent agreement with the requirements of the asymptotic approximations. When the period equidistance has been measured, all of the individual frequencies of each mode can be determined. Previously, pmode helioseismology allowed the gmode period equidistance parameter P_{0} to be bracketed inside a narrow range, between approximately 34 and 35 min. Here, P_{0} is measured to be 34 min 01 s, with a 1 s uncertainty. The previously unknown gmode splittings have now been measured from a nonsynodic reference with very high accuracy, and they imply a mean weighted rotation of 1277 ± 10 nHz (9day period) of their kernels, resulting in a rapid rotation frequency of 1644 ± 23 nHz (period of one week) of the solar core itself, which is a factor 3.8 ± 0.1 faster than the rotation of the radiative envelope.
Conclusions. 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.
Key words: Sun: helioseismology / Sun: oscillations / Sun: rotation / Sun: interior
© ESO, 2017
1. Introduction
A conference was organized in 2013 in Tucson, Arizona, to celebrate 50 years of seismology of the sun and stars. Solar seismology became known in the early 1980s as helioseismology, and was almost totally devoted in this conference to the analysis of solar p modes, sometimes referred to as solar music, as only one presentation among 87 was devoted to g modes. It was a historical presentation (Appourchaux & Pallé 2013) that confirmed how difficult the quest for g modes has been in these 40 years.
The observational gmode quest began 40 years ago, when Severny et al. (1976) reported the detection of a 160min oscillation using a modified differential Babcock solar magnetograph and nine days of data. This period of 2 h and 40 min could not be attributed to a p mode in any reasonable solar model. In 1976, the field was still in its infancy, with the phenomena being referred to as solar oscillations and with no detections of individual p modes. Although the potential of p modes in probing the deep solar interior was still not clear, the potential of g modes was obvious, as these are mostly trapped inside the solar core. The 160min oscillation thus quickly became a hot topic, and at least ten years passed before it finally became clear that it was an artifact.
Following the first identification of global p modes by Claverie et al. (1979) and the first South Pole expedition, which allowed individual low spherical harmonic degree p modes to be revealed (Grec et al. 1980; 1983), helioseismology quickly grew on different fronts: on the ground with worldwide networks that were developed by the Birmingham Solar Oscillation Network (BiSON) in the UK, the International Research on the Interior of the Sun (IRIS) in France, and the Global Oscillation Network Group (GONG) in the USA In space the initial studies of ESA and NASA led to the selection of the SOHO space project in 1988. The SOHO payload included the Global Oscillations at Low Frequency (GOLF, Gabriel et al. 1995), the Michelson Doppler Imager for Solar Oscillation Investigation (MDI/SOI), and the Variability of Irradiance and Gravity Oscillations (VIRGO) instruments. As stated by Appourchaux & Pallé (2013), in these days, the detection of g modes was thought to only be a matter of time. In the end, it still required a further 20 years of effort. Some possible results were published during these 20 years using GOLF data. One l = 1 mixed pg mode was reported at 284.7 μHz (Gabriel et al. 2002), and several candidates were suggested in the vicinity of 220 μHz (TurckChièze et al. 2004). García et al. (2007) reported the possible collective detection of a gmode signature in a lower frequency range between 25 and 140 μHz, with a possibly rapid rotation of the solar core up to five times the rotation of the remaining radiative zone. Possible detections of five l = 1g modes were reported by García et al. (2011) in the frequency range between 60 and 140 μHz, with rotational splittings in the range 850 to 950 nHz consistent with a fast sidereal rotation (1760 to 1960 nHz) of the g modes themselves, measured using the rotationcorrected, maveraged spectrum technique from Salabert et al. (2009). An extensive review of different techniques of detection, as well as possible results, can be found in Appourchaux et al. (2010). So far, none of these possible results has been confirmed.
2. Looking differently and elsewhere
The search for g modes must overcome the difficulty that the modes are trapped in the deepest part of the solar body and are evanescent in the convective layers, so that their amplitudes are not easily detectable at the solar surface. The pmode amplitudes themselves decrease dramatically with decreasing frequency (left side of Fig. 1), falling below the increasing background noise level before reaching the frequency range around the fundamental pmode oscillation, where a few mixed modes (driven by both compressibility and buoyancy) are thought to exist. In this range, roughly speaking between 0.25 and 0.30 mHz, the lowest frequency p modes and highest frequency g modes are expected to exist with similar amplitudes at the solar surface, an amplitude of the order of 1 mm s^{1} with a large uncertainty on any prediction (see again Appourchaux et al. 2010, for a comprehensive review), and they have not been detected so far. The amplitude uncertainty remains very high in the gmode domain, especially at lower frequencies, where the background noise increases. It has never been demonstrated that the signaltonoise ratio (S/N) at these lower frequencies can improve in any way. This is the fundamental reason why purported detections of g modes in this region of the spectrum have always been doubted and have never been confirmed, despite numerous efforts.
Kennedy et al. (1993) suggested that the g modes might be detected as a source of modulation of pmode frequencies, which would manifest themselves as sidelobes in the pmode frequency profiles. However, as noted by the authors, convective noise is a problem. The sidelobes cannot be narrower than the pmode lines, and the p modes are more sensitive to turbulent motions near the surface than they are to the gmode motions in the deepest solar structure. Despite serious efforts in this direction, this approach has not been successful so far. Another approach was also suggested by Duvall (2004): modifications of timedistance helioseismology to enable travel time measurements for short time intervals. This has not been successful either.
Fig. 1 Power spectrum of 16.5 years of fulldisk velocity recorded by the GOLF instrument on board the SOHO spacecraft. 
Although based on the same idea of time distance measurement, our approach is nonetheless different. We seek to identify a differential parameter in the pmode frequency spectrum that can significantly decrease the sensitivity to this convective turbulence. Such a parameter should ideally exploit the pmode frequencies measured by the GOLF instrument, taking advantage of the fact that these fulldisk p modes are modes of the lowest degrees, that is, from 0 to 3, which propagate through the solar core.
The apparent ideal candidate would be the socalled small separation, that is, the frequency separation between p modes of degrees 0 and 2. This was immediately identified, after being first resolved by the South Pole data set in 1980 (Grec et al. 1980; 1983), as a parameter probing the solar core because it measures a difference of two kernels that are essentially identical in the upper layers and different in the deepest layers. Unfortunately, it requires more than two days of measurement to be resolved, so that only g modes with periods longer than these two days could be tracked. These very long period g modes are not accessible in practice, largely because of the high density of g modes per frequency bin.
The secondbest choice is then the socalled large separation, or rather half this large separation (between even 0−2 and odd 1−3 pairs of p modes), which defines the quasivertical lines in the socalled Echelle diagram. The large separation is linked to the propagation of the acoustic waves through the Sun and is proportional to the mean density of the star. Thus, any modulation of the solar structure produces a modulation of the large separation. If these modulations are produced by the action of g modes on the geometry or density of the deep solar core, they should present the same spectral properties as those expected for the g modes. We can also seek for gmode signatures in the large separation signal. The large separation offers an additional benefit of a lower sensitivity to surface effects when adequate care is taken to measure this separation. A detailed study of the Echelle diagram shows that the almost vertical lines have a general Sshape that contains the signature of the depth of the helium ionization layer. This layer, not far below the solar surface, is again a turbulent region, which we wish to avoid. We therefore have to limit the frequency range used for measuring this large separation to the straightest line part of this Echelle diagram, in practice, to the frequency range between 2.32 and 3.74 mHz, which is defined precisely to ensure that it contains an integer number of equidistances.
Even in this limited range, the frequencies plotted in the Echelle diagram still display small oscillations around the straight lines that contain the signature of the depth of the convective zone. This boundary layer, named the tachocline, is another dynamic region that we also wish to avoid here. It is then essential to use the broadest possible frequency range in measuring the large separation to reduce the sensitivity of its fluctuations to this tachocline region. The range 2.32−3.74 mHz appears to be the best compromise between the need for limiting the range to the most linear part and the need for the broadest possible range.
There are two more reasons to justify the choice of this frequency range limitation. Below 2.3 mHz, the individual spectra are not sufficiently sensitive, and the pmode peaks are lost in background noise with 8 h of integration (see below). At the other end, above 3.7 mHz, there is a sharp identified increase of the pmode parameter sensitivity to all solar surface effects (Salabert et al. 2002). Frequency and linewidth changes with the solar cycle are much larger than in the lower pmode frequency domain, which would introduce more undesired errors in the search for the subtle modulations produced in the solar core.
We therefore use this semilarge separation, measured between 2.32 and 3.74 mHz, as a tool for the gmode search. A pmode spectrum obtained from an 8h time series reveals the discrete structure of its frequency distribution, characterized by this pseudoperiodicity on the order of 67.5 μHz. This spacing corresponds to the inverse of the round trip time T of an acoustic wave traveling through the solar center. One interesting consequence is that the fulldisk velocity signal has a strong correlation peak at a lag T of about 4 h and 6 min (which is on the order of the inverse of 67.5 μHz). This was used to advantage by the IRIS network team (Fossat et al. 1999) in an optimal gapfilling technique. In the present context, we consider T as the optimal parameter to be monitored in the search for its possible modulations produced by a change of structure located inside this core. Because we cannot filter the frequency spectrum before sampling the time series T(t), we need to be careful about aliasing, which occurs when the spectrum is folded about the Nyquist frequency. An easy and efficient technique consists of sampling the signal twice faster than the integration time, which means a 4h sampling in the present case. The Nyquist frequency is then doubled from 1/16h = 17.36 μHz to 34.72 μHz. This avoids the spectrum folding, which would be around 17.36 μHz if we were sampling at 8h intervals. The search for gmode frequencies is then possible in the lowfrequency range that corresponds to periods longer than 8 h, even if the part of the spectrum approaching 34.72 μHz (inverse of 8 h) is progressively filtered and contains aliases of frequencies coming from beyond 34.72 μHz. It still provides access to potential gmode frequencies in this range.
We worked with a GOLF timeseries calibrated in velocities as explained in García et al. (2005), which starts at 0:00:30 (T.A.I.) on April 11, 1996. It is sampled at 80 s. For the 16.5 year interval with 8h segments sampled every 4 h, we start with 36 131 possible samples, but because of data gaps (especially the socalled SOHO vacation in 1998), we finally have 34 612 8h pmode spectra.
Fig. 2 From the same 16.5year data set as Fig. 1, this is an average of 34 612 spectra computed from 8h selections taken at intervals of 4 h. 
Fig. 3 Example of one of the 34612 GOLF power spectra, limited to the range 2.32−3.74 mHz and divided by the envelope of the mean spectrum of Fig. 2. 
The general process developed for obtaining the time series T(t) can be summarized as follows:

1.
Split the GOLF time series into 8h long subtime series at 4h intervals.

2.
Compute the power spectrum of each 8h long time series oversampled by filling zeroes up to 10^{6} s to obtain a convenient 1 μHz frequency bin (the mean of all these spectra is shown in Fig. 2).

3.
Divide each power spectrum by a Gaussian function (exp(−x^{2}/ 2σ^{2})) centered at 3.22 mHz with σ = 0.39 mHz, as an approximation of the envelope of the mean spectrum, and select the [2.32−3.74] frequency range. An example of such a bandlimited envelopecorrected spectrum is shown in Fig. 3.

4.
Compute the power spectrum of each of these windowed power spectra oversampled by filling zeroes up to 125 mHz to obtain a time bin of 8 s. These power spectra of power spectra result in a function that has time as its abscissa, with their main peak at this period T of approximately 4 h, as presented in Fig. 4, which shows one example of these power spectra, or in Fig. 5, which is the sum of 34 612 similar spectra.

5.
Measure the precise location of the centroid of the main peak around its mean value of 14807 s using a secondorder polynomial fit.
We now describe these steps in some more detail.
The nearly 67.5 μHz equidistance is measured by a Fourier analysis that requires two levels of oversampling. A simple fast Fourier transform (FFT) power spectrum of power spectrum would return the initial sampling time of 80 s, which is not fine enough for our purpose. The spectra shown in Figs. 2 and 3 are oversampled to a convenient bin of 1 μHz (step 2), and then the 34612 spectra of spectra (see Figs. 4 and 5) are also oversampled to a bin of 8 s, which was chosen to be smaller than the rms scatter of T (step 4). We show that these oversampling steps do not affect the statistical properties of the final product of this step of analysis, that is, the power spectrum of T (Fig. 7).
For the most precise estimation of the equidistance through this Fourier method, it is also necessary to ensure that the entire selected frequency range contributes and not just the few highest pmode peaks around 3.2 mHz. Each of these individual frequency spectra is then “flattened” (step 3) by dividing it by a Gaussian function (centered at 3.22 mHz with σ = 0.39 mHz) that approximates the envelope of the mean (16.5 year) spectrum (Fig. 2), thus allowing optimal use to be made of the resolution provided by the 1.42 mHz bandwidth. Figure 3 is an example of such a selected and flattened individual spectrum. Each of the 34 612 individual envelopecorrected 8h pmode spectra is then analyzed by Fourier transformation, from which its respective mean pmode separation is determined. This fifth step of the analysis is not made through a fit of the GOLF signal autocorrelation peak around 4 h because that would not permit the correction of the power envelope, which is an essential step of this analysis.
Fig. 4 Power spectrum of the spectrum shown in Fig. 3. The abscissa is in seconds and shows a main peak at a time T of about 14 800 s, which is the round trip travel time of an acoustic wave through the solar center. 
We note that our spectra such as shown in Fig. 4 are computed from spectra such as shown in Fig. 3, minus their mean value, and translated at abscissa zero. This means that the first point of Fig. 4 is at zero and suppresses the 5min modulation, which would remain inside the 4h peak, whose position can then be measured more precisely. These spectra are arbitrarily normalized at the height of the main peak, since each spectrum is no longer exactly an autocorrelation that would automatically be normalized at 1 on the zero lag. Only the precise position of the 4h peak centroid (our parameter T) is relevant in our analysis, and our general guideline is then to minimize the amplitude of the random fluctuations of this parameter T around its mean value.
Each value of T is estimated by a secondorder polynomial fit around the peak profile (step 5). The choice of the width over which this fit is computed is again guided by the goal of minimizing the random fluctuations of T. A priori it should be made at least on the true bin size defined by the inverse of the range from 2.32 to 3.74 mHz, that is, about 700 s. In practice, the polynomial fit made on a range of ± 800 s around 14 800 s minimizes the scatter of T. This is on the order of the full width of this peak and not at all the best fit on the peak profile, but it is the least noisy estimate of the peak centroid. The profile of this peak could be sensitive to the way the raw intensities are corrected and then calibrated into velocities, which is indeed the case for the amplitude of the peak. However, the centroid of the peak is our relevant parameter here, and there is no reason for it to be calibration dependent, either through physics assumptions or mathematical processes. This is indeed one of the reasons why we selected this parameter in our gmode search. The peak fit sometimes fails; a limit of 240 s was therefore set to the values of the fluctuations of T to eliminate these unrealistic poorfit results. Less than 1 percent of the values is rejected, so that the number of nonzero values of T is finally 34 261. The effects are all negligible for clock and timing. The SOHO on board timer is always maintained within 20 ms of TAI. We checked that the few shifts that occurred in the GOLF sampling during the operations had no signature in all of the processes we describe here. We also checked that neglecting the SOHOSun distance variation through the orbit has no influence on our results.
Fig. 6 16.5year time series of 36 131 values of the time T fluctuations defined in Fig. 4. The number of nonzero values is 34 261. 
Figure 6 shows the time series made up from 36 131 consecutive values of T_{i} fluctuations (where 1 ≤ i ≤ 36131, including gaps). This is the new “signal” to be analyzed in the hunt for g modes. It is interesting to know that the mean value of T, when it is measured as we did, is 14 807 s, or 4 h, 6 min, 47 s. This can be seen as the mean return travel time of a sound wave from the visible surface through a diameter. The measured rms value of the scatter of T around this mean value is 52 s. The distribution around the mean value is nearly Gaussian with this σ of 52 s, with slightly extended tails on both sides that correspond to the poor fits and were cut, as mentioned above, at ± 240 s.
We checked that our parameter T is remarkably insensitive to surface effects such as solar activity or orbital velocity residuals (again, it was carefully defined for that purpose). For instance, Fig. 6 just marginally shows a slight increase of noise amplitude with time, while the GOLF signal itself shows it more substantially, as well as showing some consequences of a few technical changes on the instrument through these years that are totally invisible here. T is also insensitive to the detailed characteristics of the calibration steps used to process the GOLF signal, which will in any case reveal the pmode spectrum in our selected range.
Figure 7 shows the resulting power spectrum of the 16.5year time series of T fluctuations. Although, as discussed above, this spectrum is limited to an extremely low frequency domain (below 34.7 μHz, i.e., the lower 0.5 percent of the bandwidth shown in Fig. 1), it is remarkably flat, whereas the GOLF signal spectrum itself is characterized by a very strong 1 /f^{2} dependence in this range (Fig. 1), where it would not make sense to look for g modes. In addition, Fig. 8 shows an enlargement on the extreme lowest frequency part that illustrates the absence of signature of either the solar cycle period or the orbital period of one year, again demonstrating the very low sensitivity of T to the surface velocity.
Fig. 8 Enlargement of the lowest frequency range of Fig. 7. The arrows indicate the frequencies of the solar cycle and of the oneyear orbit. 
We now use this power spectrum to determine lowfrequency asymptotic values of gmode parameters, taking the very low frequency range into account over which this asymptotic search is best possible. In this analysis, the relevant asymptotic gmode properties for a given degree are on one hand the equidistant period spacing between consecutive radial harmonics, and on the other hand, assuming a uniformly rotating core, a constant frequency splitting value for the triplet structure of the l = 1 dipolar modes and for the quintuplet structure of the l = 2 quadrupolar modes. We show below that only degrees 1 and 2 have indeed been detected.
3. Real GOLF data versus artificial data sets
At various steps of the analysis below, it will be important to check that the analysis itself is not responsible for creating the results that we will claim. For this reason, it is essential to compare the real GOLF data with artificial datasets presenting all the GOLF pmode amplitudes and noise characteristics and for which the temporal coverage is identical to the GOLF window function.
In order to create a simulated time series with the same spectral content and the same statistical behavior as the actual GOLF data, we use the fact that the p modes are damped and randomly excited so that the data can be regarded as a stationary random process. Therefore, the periodogram of the time series, calculated as the square modulus of the FFT of the time series, can be regarded as a particular realization of the power spectrum density (PSD) of the solar signal. The PSD is the asymptotic limit of the Fourier transform of the autocorrelation of the solar signal observed for an infinite duration. When the data are evenly spaced, the FFT calculated on the same number of points as the time series for a finite observation is producing independent values, bin by bin, for both the real an imaginary parts, which are also normally distributed, so that its square modulus follows a χ^{2} distribution. It is then possible to produce simulated time series y_{sim} with the same spectral contents and statistical properties by multiplying the square root of the PSD of the data by two random normal distributions to produce the real and imaginary parts of the Fourier transform of the time series, in the form ${\mathit{y}}_{\mathrm{sim}}\mathrm{=}\mathit{FF}{\mathit{T}}^{1}\mathrm{\left(}\mathit{X}\mathrm{\right)}\mathit{,}$(1)with $\begin{array}{ccc}& & \mathit{X}\mathrm{\left(}{\mathit{\nu}}_{\mathit{i}}\mathrm{\right)}\mathrm{=}\sqrt{\mathit{PSD}}\mathrm{(}{\mathit{N}}_{\mathrm{1}}\mathrm{+}\mathrm{j}{\mathit{N}}_{\mathrm{2}}\mathrm{)}\mathit{,}\mathrm{0}\mathit{<}{\mathit{\nu}}_{\mathit{i}}\\ & & \mathit{X}\mathrm{(}\mathrm{}{\mathit{\nu}}_{\mathit{i}}\mathrm{)}\mathrm{=}\mathit{X}\mathrm{(}{\mathit{\nu}}_{\mathit{i}}{\mathrm{)}}^{\mathit{\star}}\mathit{,}\end{array}$where N_{1} and N_{2} are two normally distributed random series with half the number of points of the GOLF full time series. This method has been extensively studied by Percival (1992) and has been used in helioseismic data analysis (see, for instance, FierryFraillon et al. 1998; Appourchaux et al. 1998).
In order to estimate the PSD of the GOLF data, we used the method developed by Welch (1967). We calculated 120 periodograms corresponding to successive 90day periods of the total 16.5year time series overlapped by 50% and averaged them. The resulting spectrum was convolved by a triangular 0.1 μHz window function to obtain a smoother estimate of the PSD. Figure 9 shows the square root of the PSD superimposed on the GOLF time series FFT modulus. The PSD estimate contains all the known properties of the pmodes spectrum, as its resolution is sufficient to resolve all the components of the split p modes. Therefore, any combination of split modes that could propagate through any subsequent step of our analysis would exist in the simulated data. In contrast, the possible modulation of the pmode separation by g modes is completely suppressed in the averaging. It should be noted that the simulated data not only have the same spectral contents for the p modes, including the presence of all modes up to l = 4, with correct profile widths, amplitudes and asymmetries, but also quantitatively reproduce the surrounding noise level.
Fig. 9 Estimate of the amplitude spectrum (square root of the density power spectrum) in blue, superposed on the modulus of the Fourier transform of the GOLF time series (in magenta). This estimate was used to calculate simulated time series. 
This process can be compared to the method described by Percival (1992). We added one specific step to take the GOLF data set specificity into account, however. The GOLF data duty cycle is close to 96% in this 16.5year time series. A few gaps exist, in particular at the beginning of the time series (the SOHO vacation). Some variations of the mean amplitude can also be seen during the 16 years of observations, corresponding to the change from one side of the sodium lines to the other, as well as a yearly variation linked to the orbital ellipticity. These features would result in an imperfect independence between the frequency bins in the FFT. This was taken into account in the simulated data. For that purpose, the time series y_{sim} were multiplied by the total GOLF window function and by an amplitude modulation calculated from the GOLF dataset with a 3month moving average. By doing so, we ensure that the simulated spectra have the same statistical behavior as the original time series. It is then possible to analyze the simulated spectra with each step of the procedure that was applied to the original GOLF data.
Ten such artificial GOLF time series were produced. They were used for the first step of our analysis detailed so far: the production of the time series of T(t). The rms values of the fluctuations of these ten time series lie between 49.5 and 51 s, all close to but slightly lower than the 52 s obtained for the GOLF data set.
Fig. 10 Autocorrelation function A(ν) of the frequency range 5.5−34.7 μHz of the power spectrum shown in Fig. 7. The first value displayed is bin number 7, i.e., 0.013 μHz. The vertical lines are at 210, 630, and 1260 nHz. 
4. Rotational splitting and solar core rotation
Very many modes of degrees 1 and 2 are available that have detectable amplitudes in the very low frequency range we are considering. The best first step in this case, as suggested many years ago (Fossat et al. 1988), is to look for the rotational splitting that is assumed to be approximately independent of the radial order for a given degree of the g modes in the asymptotic approximation. The obvious approach involves computing the autocorrelation of the power spectrum shown in Fig. 7. To maximize the statistical improvement, we used the broadest possible frequency range, just avoiding the very lowest part, where a slight increase in noise is visible, and where we also know that the gmode density per unit frequency (if this spectrum contains any gmode signature, of course) becomes too crowded to permit any possible use. We also used a version of this spectrum smoothed over 6 bins, which allows contributions from peaks with frequencies that deviate slightly from the strict asymptotic equality of all individual splittings. Figure 10 shows this autocorrelation computed over the frequency range between 5.5 μHz and the Nyquist frequency of 34.7 μHz. About 100 modes of degree 1 and 170 modes of degree 2 could exist in such a broad range.
The main peak of this autocorrelation, around 210 nHz, is highly significant for two reasons. First, when we use the standard deviation σ of the autocorrelation, which has a Gaussian distribution, the peak at 210 nHz is at 4.7σ above the mean level of the autocorrelation. With the 6bin smoothing and a Gaussian distribution of 2000 bins of this autocorrelation in the range displayed in Fig. 10, that is, around 300 independent values, this implies a low probability of lower than 10^{3} for this peak to be produced by a random fluctuation. Second, this peak remains the main peak in various attempts to reduce the statistics, either by exploring a narrower frequency range, or by reducing the duration of the 16.5year initial data set to about half. Figures 11 and 12 show the autocorrelation computed on the two independent frequency ranges, from 5.5 to 13.2 μHz in Fig. 11, and from 13.2 to 34.7 μHz in Fig. 12. This shows that this correlation peak is clearly present over a very broad frequency range. This is a first indication that many asymptotic gmode signatures are present in this broad range, and this is confirmed in the next section by means of the other asymptotic parameter, namely the period equidistance. Figure 13 shows the increase in 210 nHz peak significance with increasing the duration of the time series. The length of the time series can be adjusted by adding data points to either the beginning or the end, and this figure is the average of the two possibilities. This figure shows that the autocorrelation clearly becomes significant with about 8 to 9 years of data.
The ten sets of artificial data never produce an autocorrelation value above 4σ through the same analysis, and never produce a highest peak at 210 nHz. This excludes the idea that the peak could be an artifact produced by a subtle effect of the pmode splittings themselves. As an example of these autocorrelations from artificial data sets, Fig. 14 does not show any outstanding peak, and the highest peak, not at 210 nHz, stands at 3.3σ. As we mentioned, interpreting this as a gmode splitting signature assumes a constant splitting value of the asymptotic approximation, and this is probably true only within 3 percent in the highest part of our frequency range (Berthomieu & Provost 1991). It can be verified that the smoothing of the power spectrum improves the S/N of this correlation peak. The maximum benefit is obtained by a smoothing on 6 to 8 bins, which is consistent with this small departure from the asymptotic approximation. If the same autocorrelation is computed without the 6bin smoothing, the main peak at 210 nHz is spread across 6 bins, and we can note that bin number 2 at an abscissa of 1.92 nHz is at lower than 0.02, which proves that there is essentially no correlation between consecutive bins in the power spectrum.
Fig. 11 Autocorrelation function of the frequency range 5.5−13.2 μHz of the power spectrum shown in Fig. 7. 
Fig. 12 Autocorrelation function of the frequency range 13.2−34.7 μHz of the power spectrum shown in Fig. 7. 
Figure 15 shows the histogram of the ten autocorrelations of the same frequency range of the power spectra of artificial time series, all taken together and individually centered around their mean value (which is not exactly zero because of the 6bin smoothing and the slow variations of the background noise). The rms value is 0.0166 and the highest of these 20 000 values is at 3.8σ. This indeed corresponds to a probability of 1/20 000 in a normal gaussian distribution. The 210 nHz peak of Fig. 10 (where the rms value is also 0.0166, which confirms the internal consistency of the simulation) that stands at 4.7σ, has a probability of 2 × 10^{6} to be a random outcome. These numbers are consistent with the estimated probability of lower than 10^{3} mentioned above that an autocorrelation of 300 independent values contains one of these as high as this main peak, as 300 × 2 × 10^{6} = 0.6 × 10^{3}. Its position is shown by an arrow in Fig. 15.
The statistical significance of this correlation peak is then quite strong, and this test on artificial data demonstrates that it is not produced by a subtle effect of the pmode splittings. The subsequent analysis has proved that it can really be interpreted as the signature of gmode triplets.
Fig. 13 S/N of the 210 nHz peak of Fig. 10 as a function of the duration of the time series, from 3 to 16.5 years. 
Fig. 14 Autocorrelation function of the frequency range 5.5−34.7 μHz of a power spectrum resulting from an artificial data set. No peak is visible at 210 nHz, and the highest peak stands at 3.3σ to be compared to the 4.7σ of the peak obtained from the GOLF data set. 
Fig. 15 Histogram of ten figures similar to Fig. 10 obtained with artificial data, simulated as described in text. 
Before interpreting further Fig. 10, we have to mention that special care must be taken here because of the specific process we used for this gmode detection; our detection agent is the set of lowdegree p modes, located inside the Sun and not on Earth or at the Lagrange point. The gmode signal is then seen from a reference that rotates inside the radiative and convective zones. This will have crucial consequences for the rotation analysis.
When we assume that below the convection zone, the rotation rate depends only on radius r and not on the colatitude θ (Ω(r,θ) = Ω(r) for r ≤ r_{cz}), then we can define a mean rotation rate sensed by the asymptotic g modes by ${\mathrm{\Omega}}_{\mathrm{g}}\mathrm{=}\frac{{\mathrm{\int}}_{\mathrm{0}}^{{\mathit{r}}_{\mathrm{cz}}}\mathrm{\Omega}\mathrm{\left(}\mathit{r}\mathrm{\right)}\mathrm{(}{\frac{\mathit{N}\mathrm{\left(}\mathit{r}\mathrm{\right)}}{\mathit{r}}}^{\mathrm{)}}\mathrm{d}\mathit{r}}{{\mathrm{\int}}_{\mathrm{0}}^{{\mathit{r}}_{\mathrm{cz}}}\mathrm{(}{\frac{\mathit{N}\mathrm{\left(}\mathit{r}\mathrm{\right)}}{\mathit{r}}}^{\mathrm{)}}\mathrm{d}\mathit{r}}\mathrm{=}{\mathrm{\int}}_{\mathrm{0}}^{{\mathit{r}}_{\mathrm{cz}}}\mathit{K}\mathrm{\left(}\mathit{r}\mathrm{\right)}\mathrm{\Omega}\mathrm{\left(}\mathit{r}\mathrm{\right)}\mathrm{d}\mathit{r}\mathit{,}$(4)where r_{cz} is the base of the convection zone and N(r) is the BruntVäisälä frequency. The sidereal asymptotic splittings are then be given by Ω_{g}/ 2 and 5/6 Ω_{g} for l = 1 and l = 2, respectively (Berthomieu & Provost 1991). For lowfrequency g modes, below 33 μHz with  n  ≥ 20, the relative difference between the mode splittings and these asymptotic values is smaller than 3% (see Fig. 5 of Berthomieu & Provost 1991).
In this work, the rotational splittings are not seen from an inertial frame. Instead we search for them in the signal of the lowdegree p modes that propagate through the rotating Sun. Following Komm et al. (2003), we can define the mean rotation rate sensed by p modes as the solidbody rotation of the Sun given by ${\mathrm{\Omega}}_{\mathrm{p}}\mathrm{=}\frac{{\mathrm{\int}}_{{\mathit{r}}_{\mathrm{c}}}^{{\mathit{R}}_{\mathrm{\odot}}}\mathit{\rho}{\mathit{r}}^{\mathrm{4}}{\mathrm{\int}}_{\mathrm{0}}^{\mathit{\pi}}\frac{\mathrm{3}}{\mathrm{4}}{\mathrm{sin}}^{\mathrm{3}}\mathit{\theta}\mathrm{\Omega}\mathrm{\left(}\mathit{r,\theta}\mathrm{\right)}\mathrm{d}\mathit{r}\mathrm{d}\mathit{\theta}}{{\mathrm{\int}}_{{\mathit{r}}_{\mathrm{c}}}^{{\mathit{R}}_{\mathrm{\odot}}}\mathit{\rho}{\mathit{r}}^{\mathrm{4}}\mathrm{d}\mathit{r}}\mathit{,}$(5)where ρ is the density, R_{⊙} is the solar radius, and r_{c} ≃ 0.2 R_{⊙} is the boundary below which pmode splittings become rather insensitive to the rotation rate. From the inversion of MDI splittings, Komm et al. (2003) reported Ω_{p} = 434.5 ± 5.9 nHz. In analogy with the synodicsidereal correction that is necessary when the splittings are observed from a terrestrial perspective, the expected asymptotic gmode splitting components seen in the pmode signal are given for l = 1 and l = 2 by
$\begin{array}{ccc}\mathit{s}\mathrm{\left(}\mathrm{1}\mathit{,m}\mathrm{\right)}& \mathrm{=}& \\ \mathit{s}\mathrm{\left(}\mathrm{2}\mathit{,m}\mathrm{\right)}& \mathrm{=}& \end{array}$respectively. When we use only the lowdegree modes (those used in this gmode search), no inversion is possible, but a mean rotation rate Ω_{p} = 432.5 ± 1 nHz has been reported by Fossat et al. (2002). Taking these two estimates into account, we used a mean value of Ω_{p} = 433.5 ± 7 nHz for the rotation of our “instrumental reference” (i.e., the p modes that corotate with the Sun).
Because our instrumental reference is located inside the rotating Sun, Eqs. (6) and (7) may imply a surprising consequence. If the gmode kernel rotation Ω_{g} were on the same order as the pmode kernel rotation, close to 433 nHz (case of a nearly rigid rotation down to the solar center), then Eq. (6) would imply a negative value of the splitting s(1,1), which would indeed be equal to half the value of the pmode rotation. A negative splitting means that the positions of the two tesseral peaks corresponding to m = ± 1 are reversed. It also means that an observer located on the solar surface, with a constant internal rotation down to the solar center, would “see” the g modes rotating in the opposite direction, at half the global rotation rate. As in the case of the Foucault pendulum, this is a nonintuitive consequence of the Coriolis effect, which is itself due to the nearly horizontal motions of the g modes.
Our main peak at about 210 nHz, close to half of 433 nHz (although significantly different), could correspond to this situation if it is to be seen as a signature of the s(1,1) splitting. This would mean, incidentally, that the m = 0g modes are detected by this method. Helioseismologists are used to the fact that p modes with an odd value of (m + l) are not visible in fulldisk observations. This is an observational bias that arises because the visibility is limited to half the total solar surface. In the present case, the “observation” of g modes is made from inside the Sun itself, so that this bias disappears and our measurements are sensitive to all m components. The same benefit should then exist for the m = ± 1 of l = 2, which are also not directly accessible in fulldisk velocity pmode measurements, but are also expected to contribute here in the gmode detection.
However, the autocorrelation does not see the difference between a negative and a positive splitting. We are then confronted by an ambiguity between a possible nearly solid rotation producing a splitting s(1,1) = −210 nHz and a much faster gmode rotation of more than 1200 nHz, which would produce a splitting s(1,1) = + 210 nHz.
To solve this ambiguity, we need to also detect the signature of the l = 2g modes. Their Coriolis coefficient, 5/6, is quite different so that the same 210 nHz splitting of l = 1 is associated with very different values of the l = 2 splitting in the two possible solutions. The l = 2gmodes’ splitting signature is certainly seen with a lower S/N, as shown in Fig. 10, where no dominant series of peaks can be clearly assigned to them. However, they are potentially many, and they consist of quintuplets, so that the increased statistical efficiency of the autocorrelation may still offer a chance that some of the nondominant peaks of Fig. 10 contain their signature.
In the case of slow rotation, the single s(2,1) splitting should be negative, around −60 nHz, while in the second case, its positive value is expected around 630 nHz. No significant peak can be seen in Fig. 10 around 60 nHz, while the second highest peak of this figure indeed stands near 630 nHz. The third highest peak, although not highly significant by itself, is located near 1260 nHz, and could be the signature of s(2,2).
On the other hand, the 210 nHz high peak could mathematically also be the signature of any other splitting component such as s(1, ± 1) or s(2,1). Equations (6) and (7) must be satisfied, and this will help to reduce the ambiguity.
Fig. 16 Sum S_{cor} (Eq. (8)) of the three possible gmode splitting signatures in the autocorrelation, as a function of the gmode kernel rotation Ω_{g}. It identifies a fast rotation rate of 1277 ± 10 nHz. 
Fig. 17 S/N of the main peak of Fig. 16 as a function of the duration of the time series analyzed, from 3 to 16.5 years. 
The simplest way of checking the significance of these various possible signatures of gmode splittings is to plot the sum of the autocorrelation values at the three abscissae s(1,1), s(2,1), and s(2,2) as a function of Ω_{g}, ${\mathit{S}}_{\mathrm{cor}}\mathrm{=}\mathrm{A}\mathrm{\left(}\mathit{s}\mathrm{\right(}\mathrm{1}\mathit{,}\mathrm{1}\mathrm{\left)}\mathrm{\right)}\mathrm{+}\mathrm{A}\mathrm{\left(}\mathit{s}\mathrm{\right(}\mathrm{2}\mathit{,}\mathrm{1}\mathrm{\left)}\mathrm{\right)}\mathrm{+}\mathrm{A}\mathrm{\left(}\mathit{s}\mathrm{\right(}\mathrm{2}\mathit{,}\mathrm{2}\mathrm{\left)}\mathrm{\right)}\mathit{,}$(8)with s(l,m) defined by Eqs. (6) and (7), as abscissae of the ordinate A plotted in Fig. 10.
This is shown in Fig. 16, where the sharp highest peak clearly favors the rapid rotation solution. In this figure, Ω_{g} is sampled at 1 nHz intervals, and the three abscissae s(1,1), s(2,1), and s(2,2) are rounded to the integer value of the bin number obtained from Eqs. (6) and (7). An easy Monte Carlo test demonstrates the high significance of this peak: a random selection is made of three abscissae in Fig. 10, and their three ordinates are added. One million such selections can easily be made, and repeating this one hundred times provides on average one occurrence of such a peak at 5.2σ (with respect to a mean value of 0.04 and an rms scatter of 0.029 in each set of one million selections). This occurrence always corresponds to contributions of the 210 nHz peak plus two of the three or four highest other peaks, and only one of these solutions satisfies Eqs. (6) and (7). Figure 17 shows that the S/N of this gmode rotation peak at about 1280 nHz increases with increasing data set duration. This increase does not differ much from the increase of the 210 nHz peak alone, an indication that taking the assumed l = 2 contributions into account does not only add noise. This 1280 nHz peak is consistent with the addition of the three highest peaks in Fig. 10, within their respective uncertainties of a few nHz.
For the statistical significance of the l = 2 splitting contribution, it is possible to reproduce Fig. 16 by using separately only the contribution of the l = 1 splitting, or the contributions of the first two splitting values of l = 2. This is shown in Fig. 18 in blue for l = 1 and red for l = 2. The blue curve shows the ambiguity on the rotation rate created by Eq. (6), which is visible as a symmetry around 866 nHz (2 × 433), so that two dominant peaks that are equally high support either a slow or a fast rotation. The red curve, corresponding to l = 2 and Eq. (7), displays a symmetry around 520 nHz (433 × 6/5) and clearly helps selecting one unique solution that is precisely located on one of the other two peaks. They are slightly rescaled here for better visibility, and Fig. 16 shows the sum of these two curves. The rapid rotation of about 1277 ± 10 nHz could therefore almost be measured by the only contribution of the l = 2g modes, but the case becomes very compelling with the two contributions taken together. The error bar of 10 nHz around 1277 nHz takes the three contributions of the three correlation peaks into account. We can then conclude that the autocorrelation shown in Fig. 10 indeed contains the significant signatures of gmode splittings of degrees 1 and 2. This indicates a rapid rotation of the solar core, and we can move to the next step of the analysis, but we return to this important point of the rotation in the general discussion.
Fig. 18 Same as Fig. 16, but computed separately for the unique contribution of l = 1, A(s(1,1)) in blue, or for the contribution of two components of l = 2, A(s(2,1)) + A(s(2,2)) in red. This illustrates that the l = 2 contribution confirms one of the two possible rotation rates detected by l = 1 with ambiguity. 
It is also interesting to compare Fig. 18 with the same calculation made on the same artificial data set, shown in Fig. 19. It is clearly more random.
5. Equidistant period spacing. The case of l = 1
For a sufficiently low frequency, the secondorder asymptotic period is given by ${\mathit{P}}_{\mathit{n,l}}\mathrm{=}{\mathit{P}}_{\mathit{l}}\mathrm{\left(}\mathrm{\right}\mathit{n}\mathrm{}\mathrm{+}\mathit{l}\mathit{/}\mathrm{2}\mathrm{+}\mathit{\theta}\mathrm{)}\mathrm{+}\mathit{O}\mathrm{(}\mathrm{1}\mathit{/}{\mathit{P}}_{\mathit{n,l}}\mathrm{)}\mathit{,}$(9)where P_{l} is the equidistant period spacing, defined by ${\mathit{P}}_{\mathit{l}}\mathrm{=}{\mathit{P}}_{\mathrm{0}}\mathit{/}\sqrt{\mathrm{\left(}\mathit{l}\mathrm{\right(}\mathit{l}\mathrm{+}\mathrm{1}\mathrm{\left)}\mathrm{\right)}}\mathit{,}$(10)and θ is a phase factor that tends to a constant when the radial order increases (Berthomieu & Provost 1991; Provost et al. 2000). P_{0} is directly related to the denominator of the asymptotic gmode rotation kernel given by Eq. (4), ${\mathit{P}}_{\mathrm{0}}\mathrm{=}\frac{\mathrm{2}{\mathit{\pi}}^{\mathrm{2}}}{{\mathrm{\int}}_{\mathrm{0}}^{{\mathit{r}}_{\mathrm{cz}}}\frac{\mathit{N}\mathrm{\left(}\mathit{r}\mathrm{\right)}}{\mathit{r}}\mathrm{d}\mathit{r}}\mathrm{\xb7}$(11)Pmode helioseismology has bracketed the value of P_{0} to be around 34 min, with an uncertainty smaller than one minute. The two asymptotic values of P_{l} should then be approximately 24 min for l = 1 and 14 min for l = 2. We separately address the case of l = 1 , which provides a highly significant result. The case of l = 2 is expected to be more difficult.
In this section we assume that all of the g modes participating in our analysis indeed precisely follow this equidistant period spacing asymptotic approximation. We discuss this hypothesis in the general discussion section.
Given the expected equidistant period spacing of about 24 min, there must exist about 100 gmode dipole frequencies in the range explored to produce Fig. 7. We do not know how many of these really contribute to this figure because they are not individually detectable. We know from Figs. 11 and 12 that they are many, however. A priori we should explore the broadest possible frequency range to search for the equidistant period spacing. Another difficulty arises, however, which is the density of g modes in the lowest part of this range. Because they are equidistant in period, the consecutive modes approach each other with decreasing frequency. At 6 μHz, two consecutive frequencies are separated by about 51 nHz, while their s(1,1)−s(1,−1) splitting is 420 nHz. This means that more than eight consecutive g modes overlap. This density provides an interesting statistical benefit in the autocorrelation, while it might create confusion in this part of the analysis.
Provost & Berthomieu (1986) and ChristensenDalsgaard & Berthomieu (1991) estimated that the asymptotic approximations could be regarded as valid within about 10^{3}. In our frequency range this means a few bins. The 6bin smoothing we used for the autocorrelation has indeed proved to optimize the visibility of most of the following figures.
The search for equidistant period spacing makes use of the assumed known splitting to produce one model of many multiplets that cover part of or all the frequency range used to measure the rotation. This model has two free parameters: the equidistant period spacing P_{l}, and the position of this equidistant series, which is defined by the value P_{min,l} of the shortest period of the model.
A 2D exploration (P_{l} in abscissa and P_{min,l} in ordinate) can be computed by adding all of the individual contributions of the spectrum of T (Fig. 7) along the abscissae of this model. This can be written as ${\mathit{C}}_{\mathit{l}}\mathrm{\left(}{\mathit{P}}_{\mathit{l}}\mathit{,}{\mathit{P}}_{\mathrm{min}\mathit{,l}}\mathrm{\right)}\mathrm{=}\sum _{\mathit{i}\mathrm{=}\mathrm{0}}^{{\mathit{N}}_{\mathit{l}}\mathrm{}\mathrm{1}}\sum _{\mathit{m}\mathrm{=}\mathrm{}\mathit{l}}^{\mathit{l}}{\mathit{P}}_{\mathrm{S}}\mathrm{\left(}\mathrm{\right(}{\mathit{P}}_{\mathrm{min}\mathit{,l}}\mathrm{+}\mathit{i}{\mathit{P}}_{\mathit{l}}{\mathrm{)}}^{1}\mathrm{+}\mathit{m}\mathrm{\times}{\mathit{s}}_{\mathit{l}}\mathrm{)}\mathit{,}$(12)where P_{S}(ν) is the value of the power spectrum of T at the frequency ν. The frequencies are approximated to the integer number of bins, which is close enough with the 6bin smoothing used in all the next steps of the analysis. N_{l} and s_{l} = s(l,m) /m are respectively the number of modes and the single splitting value used in the model for degree l.
For l = 1, using s_{1} = 210 nHz, we started with the broadest possible frequency range of 5.5−34.7 μHz, and then attempted to reduce it by starting above 5.5 μHz and determining whether this would increase or decrease the S/N. The best result is obtained with a set of N_{1} = 76 modes, starting just above 7 μHz.
Figure 20 shows C_{1}(P_{1},P_{min,1}), which uses this model, including the three components of 76 assumed l = 1gmodes triplet ν_{n,1} and ν_{n,1} ± 210 nHz. This shows the resulting summation of these 228 values, displayed for a 120 s range of the equidistant period spacing centered on 24 min, 20 s (abscissa), and a 2000 s time window starting at 30 000 s (ordinate axis). The horizontal range is much wider than the estimated uncertainty in the previous knowledge of the period spacing, and the vertical range is wider than this spacing, so that the first two periods of the equidistant series can be seen. The two red spots horizontally define an equidistant period spacing of 1443 s with at worst an uncertainty of 1 s, and they are indeed vertically separated by 1443 s. The shortest period of the series is seen at about 30 360 s. Figure 21 shows the series of maxima of each line in Fig. 20. This demonstrates that there is no other maximum that can compare to the red spots in the broad range of the 2D figure that covers more than all the possibilities (4.7σ against 2.1σ for the second highest value). It also reveals, quantitatively more precisely than the color view, the first two gmode periods of the series, at 30 360 ± 15 s and 31 803 ± 15 s, respectively, which are indeed separated by the spacing of 1443 s. This is obtained by construction because, as the above formula shows, replacing P by P + P_{1} leads to a replacement of P_{S}((P + iP_{1})^{1}) by P_{S}((P + (i + 1)P_{1})^{1}) and thus to a shift of the sum by one term only. When we assume that the asymptotic regime has been reached in this frequency range and by comparing these values with the frequency table of Mathur et al. (2007), these peaks, corresponding to frequencies of 32.94 ± 0.015 μHz and 31.44 ± 0.015 μHz, can tentatively be identified with the signatures of dipolar g modes of radial orders n = −20 and n = −21, that is, P_{min,1} = P_{− 20,1} in Eq. (12). They might differ slightly from the real frequencies of these two modes if the strict asymptotic regime starts beyond the 20th radial harmonics. We return to this point in the general discussion.
Fig. 20 Sum C_{1}(P_{1},P_{min,1}) (Eq. (12)) of the 228 values of a model based on 76 gmode dipole triplets with the splitting value s_{1} = 210 nHz, applied to the power spectrum of Fig. 7, smoothed over 6 bins. The two free parameters of the model are the equidistant period spacing P_{1}, defining the abscissa, and the period P_{min,1} of the first selected mode, defining the ordinate axis. The first (lowest) sharp peak is located at 30 360 ± 15 s (FWHM) and covers the radial order range from −20 to −95, the second (upper part of the figure) covers the radial order range from −21 to −96. The equidistant period spacing is 1443 s, or 24 min 03 s, within less than 1 s. 
Fig. 21 Each value is the maximum of a row of Fig. 20, so that the abscissa is the ordinate of Fig. 20. The interval between these peaks is the l = 1 equidistant period spacing of 1443 s. The power scale is the mean power per individual gmode peak. 
Fig. 22 As in Fig. 20, the sum of the 228 values of a model based on 76 gmode dipole triplets and applied on an artificial data set, with an assumed splitting value of 267 nHz that corresponds to the highest correlation peak in Fig. 14. No significant equidistant period is shown, as demonstrated quantitatively by the next figure. 
Fig. 23 Each value of this figure is the maximum of a row of Fig. 22, so that its abscissa is the ordinate of Fig. 22. This figure relative to an artificial data set must be compared to Fig. 21 obtained with the real GOLF data. 
At this step it is useful to make a comparison with the same analysis made on an artificial data set. The equivalent of Fig. 20 with the artificial time series we used earlier is shown in Fig. 22, where no specific period equidistance is clearly visible. The equivalent of Fig. 21, shown in Fig. 23, confirms that it is not the analysis that produces the results. We note that Figs. 21 and 23, as well as the next figures of the series (Figs. 24, 25, 27, and 35) are all presented with the same power unit of s^{2}Hz^{1}, as the original spectrum of Fig. 7, so that the amplitudes of the displayed peaks are directly comparable.
Given the very high S/N of Fig. 21, the idea of independently carrying out the two steps of this analysis (rotation on one hand, equidistance on the other hand) can be attempted by assuming that the rotational splitting is unknown and by assuming a contribution of the central component of the split triplet (the signature of the simple splitting s(1,1) would not be visible in the autocorrelation otherwise). Then, the search for the equidistant period spacing can be made by the same analysis by reducing the model to a simple series of equidistant central periods, instead of the triplets. We have to expect a loss of sensitivity of about $\sqrt{\mathrm{3}}$, but it can still be successful because these peaks are at 4.7σ. The result is shown in Fig. 24. The S/N is indeed significantly reduced, but the series of equidistant periods starting at 30 360 s is still clearly identified from the central gmode frequencies alone. After the 210 nHz peak in the autocorrelation, this is a second and independent strong indication of a signature of asymptotic g modes of degree 1 in the power spectrum of Fig. 7.
We can conclude from this exercise that the central components of the g mode indeed clearly contribute to the previous figures, and we can try to proceed one step further and directly detect their mean contribution in the power spectrum of T(t) itself. We cannot see the triplets individually, but since we know their 76 periods, we know 76 individual frequencies and can compute the mean of 76 sections of Fig. 7 power spectrum, each one centered on one of these 76 mode frequencies and broader than the total splitting of s(1,1)−s(1,−1). The statistical benefit for each one of the three tesseral components is comparable to what was shown in Fig. 24, where we were searching for the best solution for efficiently adding 76 individual central components. Now we have this best solution (the positions of the 76 periods, equidistant by 1443 s), and we test whether their addition in the frequency domain will also reveal the two tesseral splitting components. The following formula was used to derive Fig. 25, over the limited range of 1.2 μHz: $\mathrm{M}\mathrm{\left(}\mathit{\nu}\mathrm{\right)}\mathrm{=}\sum _{\mathit{n}\mathrm{=}\mathrm{}\mathrm{20}}^{95}{\mathit{P}}_{\mathrm{S}}\mathrm{(}{\mathit{\nu}}_{\mathit{n,}\mathrm{1}}\mathrm{}\mathit{\nu}\mathrm{)}\mathit{,}$(13)where ν is the frequency and ν_{n,1} is the inverse of P_{n,1} = P_{− 20,1} + (  n  −20) × P_{1} (rounded to the integer value of the bin number).
Figure 25 exhibits the two split tesseral components with an S/N that can be compared to the S/N of Fig. 21, and with a splitting value that is indeed consistent with the 210 nHz value seen in the autocorrelation, and which can be here estimated at 209.5 ± 2 nHz. Figure 25 also shows that the central component does not only contribute, but contributes with an amplitude about 1.5 times larger than the other two components. This will have to be understood, and it explains, or justifies, the fact that the double splitting s(1,1)−s(1,−1) is not detected in the autocorrelation: with these relative amplitudes, its signature must be three times smaller, so that it can easily disappear in the background noise.
It is also interesting to note that the mean heights of these peaks can be translated into the modulation amplitude of the sound wave travel time through the solar diameter. The product of height × linewidth is the square of this amplitude. With a mean height of 10^{7} s^{2} Hz^{1} and a line width of about 10 nHz, this results in a modulation amplitude of about $\sqrt{{\mathrm{10}}^{1}}\mathrm{\simeq}\mathrm{0.3}$ s, or as a fraction of the travel time of 4 h, a relative amplitude modulation of 2 × 10^{5}.
Fig. 24 Sum of the 76 values based on the same model based on 76 g modes; this model is simplified to the 76 central gmode frequencies, thus ignoring knowledge of their splitting. 
Fig. 25 Mean profile M(ν) of 76 equidistant (in period) dipole modes whose central frequencies are defined by Fig. 20, shown for a 1.2 μHz bandwidth (Eq. (13)). The triplet splitting of this profile, even when contaminated by the background noise and many overlapping neighbors, is nevertheless directly measured on the modes themselves, at 209.5 ± 2 nHz. 
6. Equidistant period spacing. The case of l = 2
The same general analysis can be applied to the search for the quadrupole modes. This is of course significantly more difficult because of the smaller detected amplitudes. In the autocorrelation, the main peak of the l = 1 splitting signature is the sum of (at most) 200 products of individual peaks corresponding to m = 0 and m = ± 1. In the same frequency range, there are about 170 frequencies of quadrupolar modes, composed of quintuplets. The first signature of their splitting in the autocorrelation is then the sum, at most, of 4 × 170 = 680 similar products of individual peaks corresponding to neighbor tesseral orders. Despite this encouraging ratio of 680 versus 200, the mean signature in the autocorrelation is about twice smaller (relative heights of the peaks at 630 and 210 nHz), which indicates a mean amplitude, individual mode by individual mode, smaller by at least a factor 4. An important additional difficulty is the very dense overlap of these modes, especially in the lowest frequency range. They are closer to each other than the dipole modes by a factor $\sqrt{\mathrm{3}}$, while their total splitting width is broader by a large factor: 2520/420 = 6. This makes it difficult to explore the very lowest part of the frequency range (where each quintuplet is contaminated by the contribution of more than 80 neighbors) and thus reduces the statistical benefit.
On the other hand, we have in hand the knowledge of the splitting and even of the equidistance period spacing, which is firmly constrained by the asymptotic equation Eq. (10), which implies a value close to ${\mathit{P}}_{\mathrm{2}}\mathrm{=}{\mathit{P}}_{\mathrm{1}}\mathit{/}\sqrt{\mathrm{3}}\mathrm{\simeq}\mathrm{833}$ s. We can then attempt the search for equidistant period spacing in a limited range near the expected value of 833 s.
In the case of Fig. 20 for l = 1, we used a gmode model of 76 triplets covering the frequency range from 7 to 33 μHz. This does not work with l = 2 because of the large overlap near 7 μHz. In an attempt to find a compromise between the need for a good S/N and the excess of confusion created by the overlap of the modes, we tried here the frequency range starting at 8 μHz. Figure 26 is then similar to Fig. 20. It shows C_{2}(P_{2},P_{min,2}) (Eq. (12) for l = 2) for a model of N_{2} = 112l = 2gmodes quintuplets and a single splitting value of s_{2} = 630 nHz. This quantity is displayed over a range of possible equidistance values P_{2} between 825 and 840 s, and a range of 1500 s for the first value P_{min,2} of the 112 periods, starting at 29 000 s. We can indeed see a possible signature close to the expected equidistance of 833 s, with the visibility of the first two periods separated by about 833 s.
Figure 27 is then a vertical cut of Fig. 26 at 832.8 s, extended to the first three periods. It confirms the equidistant period spacing of 832.8 s with an S/N of 3.65 and also shows the beginning of the decrease in efficiency of the method when pushing it too far toward very long periods. We show more reasons to accept this l = 2 detection in the general discussion.
When we again compare this with the same theoretical frequency table of Mathur et al. (2007), the first mode of our model can tentatively be identified with the radial order n = −33 (i.e., P_{min,2} = P_{− 33,2} in Eq. (12)) at its period of 29 380 ± 30 s, and thus a frequency of 34.04 ± 0.03 μHz. We return to the error bars in the general discussion.
Once again, the same analysis can be made on the artificial data set, and Fig. 28 confirms that it does not show any signature of period equidistance.
Fig. 26 Sum C_{2}(P_{2},P_{min,2}) (Eq. (12)) of the 560 components of 112 l = 2 quintuplets, with a single splitting value of s_{2} = 630 nHz, and the two free parameters are the period spacing P_{2} in the abscissa, in the range 825−840 s, and the first central period of the series P_{min,2} in the ordinate, as in Fig. 20. The equidistance is seen at 832.8 s. 
Fig. 27 Vertical cut of Fig. 26 at 832.8 s, extended to a third period. It shows a good S/N and clearly locates the first period of the series at 29 380 ± 30 s. This peak is the mean of the 560 components of the 112 l = 2g modes in the range of radial orders from −33 to −144. The other two peaks cover the ranges (−34,−145) and (−35,−146). 
Fig. 28 Same as Fig. 26 for the artificial data set. This cut is naturally random as no l = 2 splitting was identified with the simulated data set. Consequently, no equidistant period spacing can be identified. 
7. Discussion
The reliability of important results such as the solar gmode detection and the solar core rotation measurement requires a strong statistical and physical case. The statistical significance of all results presented so far is comfortable, and we have also shown that the correlation peak at 210 nHz can be obtained from independent data sets in two nonoverlapping frequency ranges. We discuss the solar core rotation first, and then show that independent results can also be obtained for the period spacing equidistances.
7.1. Rotation of the solar core
An interesting additional exercise can be carried out here to reinforce, were it needed, the reliability of this measurement. It consists of considering the problem in the opposite direction. We assume that the autocorrelation does contain the signatures of the three gmode splittings s(1,1), s(2,1), and s(2,2) and that we do not know their frequencies, they just have to follow, all together, Eqs. (6) and (7). We do not know the value of Ω_{g}, and this is the most interesting part of this exercise. Then we take the pmode rotation Ω_{p} as the unknown parameter of this exercise. We scan a broad range of Ω_{p} possibilities, and for each possibility we scan the values of Ω_{g} that cover the range of Fig. 16, and we select the highest value of the sum Eq. (8), and we finally plot this highest sum as a function of Ω_{p}. We still do not need to know Ω_{g}. The goal is to determine the best value of Ω_{p} that is consistent with Eqs. (6) and (7). In other words, we use the only assumption of the existence of g modes as a tool for measuring the rotation rate sensed by the p modes. Figure 29 shows that the g modes indeed confirm the 433 nHz rotation rate of the p modes, just a little less accurately than their direct surface measurements. We already know (from Figs. 16 and 18) that this Ω_{p} = 433 nHz corresponds to Ω_{g} = 1277 nHz. We show here that taking Eqs. (6) and (7) together, applied to the autocorrelation of Fig. 10, without any a priori knowledge of Ω_{p} or Ω_{g}, provides an optimal solution with the pair Ω_{p} = 433 nHz and Ω_{g} = 1277 nHz. As 433 nHz is already well known to be the mean global pmode rotation, this confirms the strong confidence in the measurement of the gmode mean rotation Ω_{g} = 1277 nHz. The uncertainty of this rotation rate is due to an uncertainty in both the measured gmode splittings and the 433 nHz pmode rotation. They are both very small, and the total does not exceed ± 10 nHz, that is, the uncertainty is smaller than 1 percent.
The mean rotation rate sensed by the asymptotic gmodes Ω_{g} = 1277 ± 10 nHz is a weighted average below the convection zone (r ≤ r_{cz}). The weighting function is given by the kernel K(r) (Fig. 30). However, the pmode splittings have already given a good estimate of the mean rotation rate in the radiative interior r_{c} ≤ r ≤ r_{cz}. In this zone, all results published so far are compatible with a radiative zone that rotates rigidly at a mean rate of 433 ± 10 nHz. The area of the kernel K(r) above r_{c} = 0.2 R_{⊙} is 30 ± 1%, which leads to an estimate of the mean value of the rotation rate below r_{c} of 1644 ± 23 nHz (or oneweek period), that is, a mean rotation of the solar core 3.8 ± 0.1 times faster than the rotation of the mean radiative zone.
A rapid rotation of the solar core has been suggested by García et al. (2007), who also gave an approximate estimate of the P_{1} value with a precision much lower than that of our analysis: rotation between three and five times faster than that of the radiative zone, and P_{1} between 22 and 26 min (obtained in a higher frequency range not yet very close to the asymptotic period equidistance). Our more accurate results do not contradict the P_{1} estimate, but they cannot be reconciled with the more precise values of the rotation given later by García et al. (2011), which were based on the possible identification of five dipolar g modes: the claimed high splittings of these modes, although suffering large scatter, would imply a gmode rotation of more than 1800 nHz, to be compared to our result of 1277 ± 10 nHz.
Fig. 29 Assumed unknown pmode rotation rate tested against the optimum visibility of the three gmode splitting values that are assumed to exist, but are unknown so far. The wellknown 433 nHz pmode rotation rate is confirmed by the g modes, only a little less accurately than when measured directly on the p modes themselves. 
7.2. Equidistant period spacing of the asymptotic gmode
The broad frequency range across which this parameter has been measured from the dipole l = 1g modes implies an interesting accuracy on its value, as 76 equidistant modes have been exploited. The degree of freedom for targeting the last mode from the first mode is extremely small. From the results displayed in Figs. 20 and 21, the equidistance can be estimated to be P_{1} = 1443.1 ± 0.5 s.
Figure 21, produced by the sum of 76 equidistant (in period) symmetrical triplets, does not permit much doubt about this interpretation as a sequence of l = 1g modes. However, two different tests can be made to reinforce this assessment. The very high S/N once again permits us to divide either the period range or the original time series itself into two independent parts that can then provide totally independent results. Figures 31 and 32 illustrate these two possibilities.
Figure 31 was obtained in the same way as Fig. 21, with the two independent time series of the first and second halves of the original 16.5 years. A factor of 0.7 was applied on the second half to display the two curves with a similar scale, as the background noise is about 30 percent higher in the second time series analysis. Despite an obvious reduction of the S/N in comparison with Fig. 21, the first two periods of the series are still clearly visible and stand precisely at the same values of 30 360 and 31 804 s. We recall that these two analyses are made on completely independent data sets.
Fig. 31 Same as Fig. 21, computed on the same sequence of 76 layers of the power spectrum, equidistant in periods. This sum is computed separately here for the independent first half (red) and second half (blue) of the 16.5year time series. The power scale of this and the next three figures is arbitrary, as they contain two curves that are rescaled to facilitate comparison. 
Fig. 32 Following Figs. 21 and 31, this illustration uses two independent sequences of 38 intervals that are equidistant in periods, the first half (shorter periods, radial order range from −20 to −57 for the first of the four peaks) plotted in blue, and the second half (longer periods, radial order range from −58 to −95 for the first peak, rescaled by a factor 0.8) shown in red, both computed on the complete time series of 16.5 years. The first four periods of the sequences are shown. 
In Fig. 32, the complete 16.5year time series is used, but the sequence of 76 equidistant periods is separated into two. We used a first range that contains the first 38 periods with the same starting value of 30 000 s, and another range that also contains 38 more periods, starting at 84 834 s, which is 38 periods beyond the first range. This figure illustrates the first four periods of two series of 38 equidistant periods that are completely independent. Not only is the 1443 s periodicity clearly identified in these two independent parts, but the second periodicity is seen as an exact extension of the first. This gives confidence in the high accuracy of the period equidistance measurement, and it shows that many g modes are indeed contributing in a very broad range of periods. We note that the second series (red) shows peaks that look significantly broader. This is due to the display in period, while the power spectrum itself is obtained in frequency. At lower frequencies, the same frequency bin contains a broader interval of periods.
The same two separations into independent time series and independent frequency ranges can be analyzed with the l = 2g modes. We have to expect noisier results. The visibility is indeed not as good as it is with l = 1, but as shown by Figs. 33 and 34, these separations clearly confirm that the l = 2 modes are really present during the 16 years and also in the broad range of frequencies or periods. Consequently, this detection also validates the result for the rotation of the solar core.
Fig. 33 Figure to be compared to Fig. 27 for the l = 2g modes, computed for the same sequence of 112 layers of the power spectrum, equidistant in periods. This sum is computed here separately for the independent first half (red) and second half (blue) of the 16.5year time series. The first three periods of each sequence are shown. 
Fig. 34 As in Fig. 32 for l = 1, this illustration for l = 2 uses two independent sequences of 56 intervals equidistant in period, the first half (shorter periods, ranging from n = −33 to −88) is plotted in blue and the second half (longer periods, ranging from n = −89 to −144, slightly rescaled) is shown in red, both computed on the complete time series of 16.5 years. The first three periods of each sequence are shown. 
One last point has to be addressed in connection with this asymptotic gmode equidistant period spacing: a possible departure from the asymptotic regime. All of the results shown so far assume a perfect equidistance in period, and their impressive sharpness supports this assumption. However, a small departure from this equidistance cannot be excluded and would slightly modify both the value of the real equidistance and the precise values of individual periods and frequencies. Provost & Berthomieu (1986) cautioned that at such accuracy of an equidistance measured within less than 1 s, we cannot pretend to have reached the asymptote as early as at n = −20. The difference between the periods and the linear asymptote should decrease following an hyperbolic trend of the radial order n, and the difference between consecutive periods should become constant within less than 1 s well beyond n = −20. This was tested by taking a nonzero value of α in Eq. (14), derived from Eq. (9), thereby introducing the hyperbolic decreasing term in the definition of our 76period model, ${\mathit{P}}_{\mathit{n,}\mathrm{1}}\mathrm{=}{\mathit{P}}_{\mathrm{min}\mathit{,}\mathrm{1}}\mathrm{+}{\mathit{P}}_{\mathrm{1}}\mathrm{\left(}\mathrm{\right}\mathit{n}\mathrm{}\mathrm{}\mathrm{20}\mathrm{)}\mathrm{+}\frac{\mathit{\alpha}}{\mathrm{\left}\mathit{n}\mathrm{\right}}\mathrm{\xb7}$(14)We then scanned a broad range of values of the parameter α to look for a possible improvement of the S/N of Fig. 21. The result is that, indeed, a small but significant improvement can be obtained with a positive value of α. The contrast of Fig. 21 increases by 6 percent with α = 1350 s, and decreases irregularly on both sides of this value. This value of α = 1350 s corresponds to an asymptotic equidistance of 1443.6 s and to a departure from this equidistance of about 3 s at n = −20. It also implies a slight change in the period of the first identified g mode of radial order n = −20. A very precise adjustment can be made by optimizing the cumulated visibility of the mean triplet profile of these dipolar g modes. When we recompute the sum of the 76 spectral bands centered on the 76 frequencies defined by our model including this additional term with α = 1350 s, and then the average of this profile and its leftright flipped image, we have a sensitive tool for targeting the precise position of the central frequencies (if it is not precisely located, the central peak is not at zero in this plot and the average of the profile and its leftright image is clearly not optimal). Incidentally, optimizing the visibility of the central peak simultaneously optimizes the combined visibility of the split sidelobes, which are the average of the two m = ± 1 components. This means that these two split components are precisely symmetrical in frequency around the central peak. Figure 35 shows this mean symmetrized l = 1 triplet profile, averaged over the 76 modes that are not quite equidistant in period. The splitting is clearly visible at its value of 209.5 nHz. This fine analysis artificially equalizes the amplitudes of the m = ± 1 tesseral components, but it has the additional benefit of improving the visibility of the triplet by a factor $\sqrt{\mathrm{2}}$ thanks to its precise leftright symmetry, and it brackets the value of the n = −20 period at 30 385 ± 3 s with α = 1350 s. However, this period and the equidistance are not independent, and both slightly depend on α.
Fig. 35 Mean symmetrized profile of the 76 dipolar g modes, following the model that includes a slight departure from the asymptotic equidistance. The comparison with Fig. 25 shows that taking this slight departure into account improves the visibility. 
We conclude this discussion by setting the asymptotic equidistance at 1443.6 s with a conservative uncertainty of about 0.7 s. The period of the radial order −20 g mode, which is not yet on the asymptotic series, is then 30385 ± 30 s. Its distance from the asymptote is 1350/20 = 67.5 s.
In the case of l = 2, the first accessible g modes correspond to radial orders around 35. They are consequently closer to the asymptotic approximation, and no departure from this approximation can be detected by a similar analysis. The best choice for a similar α parameter is zero. The asymptotic equidistance remains at 832.8 s with an uncertainty on the same order: 0.7 s. Assuming a perfect equidistance gives an uncertainty smaller than 0.7 s, but we have to take into account the small, but not measured, departure from the exact equidistance of the detected modes. For the same reason, even if we can measure the first period within a few seconds, this apparent accuracy assumes that everything is perfectly asymptotic. Taking the existence of a small but not detected slight departure from this linear approximation into account, the true uncertainty on the first period of our model can again be estimated to 30 s. This is the number that we have mentioned above.
Taken together, the measured asymptotic equidistances of the modes of degrees 1 and 2 provide a value of the asymptotic parameter P_{0} (which is usually given in min rather than in seconds) P_{0} = 34 mn 01 ± 1 s.
7.3. About the gmode amplitudes
An interesting question arises concerning the amplitudes of the detected modes. This amplitude is not directly measured as a velocity of the waves’ motions, it is indirectly detected by its effects on the sound waves’ travel time integrated through the volume of the gmode kernels. It is then expressed in seconds, to be taken as a fraction of this travel time. A careful calibration of the power spectrum shown in Fig. 7 and of all ensuing byproducts such as the mean profile shown in Fig. 25 leads to a mean amplitude of 0.3 s for each individual component of each contributing dipolar g mode. This is a small fraction, on the order of 2 × 10^{5}, of the 4h travel time. It is less for the l = 2 components, close to 0.1 s. These modulations occur only inside the gmode kernels, however, and this means mostly inside the solar core itself where the sound waves travel much faster than in the radiative and convective envelopes. Its relative amplitude with respect to the travel time through the solar core is then rather on the order of 2−3 × 10^{4}, which is not so small, considering the many existing g modes that are revealed by this analysis.
The meaning of this amplitude in terms of the frequency modulation rate on the p modes themselves remains to be investigated. It seems clear, however, that the search for g modes in the side lobes of the p modes could not give a similar result, either because the side lobes would have a tiny amplitude, much lower than the noise level between the modes, or because the nonlinear effects would spread the gmode signature across many peaks, as noted by Lou (2001), and would only contribute to the noise.
When this is considered together with the rapid rotation of the deep solar layers implied by the values of gmode splitting identified here, a new and very interesting challenge is raised for the study of solar core dynamics, as well as for physics itself. Previous studies dedicated to the individual search of some gravity modes have shown the potential of the Sun to place constraints on particle physics, including the search for dark matter in our Galactic environment (TurckChièze & Lopes 2012; TurckChièze et al. 2012).
8. Results and conclusions
An original method has been applied to the search for solar g modes, based on an analysis of the round trip time of acoustic waves traveling through the solar diameter. This travel time is affected by temporal modulations, which are presumably caused by the presence of g modes that shake the structure of the core region. These modulations provide the signatures of many solar g modes with periods ranging between 0.4 and 2 days. This very low frequency range makes it possible to take statistical advantage of the asymptotic approximations, and thus to obtain extremely accurate results. A time series of 16.5 years of GOLF data was used for this analysis, and the results can also be obtained by separating this time series into two independent halves. This provides an extremely high level of confidence. The quantitative results include the asymptotic equidistance of gmode periods, the periods themselves, and the rotation of the solar core.
The mean rotation rate sensed by the asymptotic g modes, Ω_{g} = 1277 ± 10 nHz, is a weighted average below the convection zone (r ≤ r_{cz}). This leads to a mean value of the rotation rate below r_{c}, of 1644 ±23 nHz (oneweek period), that is, a mean rotation of the solar core that is 3.8 ±0.1 times faster than the mean radiative zone rotation.
This rapid rotation nevertheless remains difficult to explain by models describing a pure angular momentum evolution without adding new dynamical processes such as internal magnetic breaking, which could have appeared when the Sun was young (TurckChièze et al. 2010). Future more detailed analyses of the inversion will undoubtedly shed light on this intriguing aspect of solar dynamics. Clarifying the transition between a nearly solid rotation of the deep radiative envelope and a rapid rotation of the core just beneath this envelope will be an interesting challenge.
Besides the rotation rate of the solar core, the other quantitative parameters obtained by this analysis are the asymptotic period equidistances:

for l = l, this is P_{1} = 1443.6 s = 24 min 03.6 s ± 0.7 s, and

for l = 2, this is P_{2} = 832.8 s = 13 min 52.8 s ± 0.7 s.
These two values taken together provide the parameter in Eq. (11):

P_{0} = 2041 s = 34 min 01 s ± 1 s.
We have also proposed a tentative identification of the shortest detected gmode periods:

for l = 1, this is P_{− 20,1} = 30385 ± 30 s, ν_{− 20,1} = 32.91 ± 0.03 μHz (very likely). The next 75 l = 1 mode periods can be obtained by P_{n,1} = P_{− 20,1} + (  n  −20)(P_{1}−67.5/  n ) for −95 ≤ n ≤ −21.

for l = 2, this is P_{− 33,2} = 29 380 ± 30 s, ν_{− 33,2} = 34.04 ± 0.03 μHz (likely). The next 111 l = 2 modes periods can be obtained by P_{n,2} = P_{− 33,2} + P_{2}(  n  −33) for −144 ≤ n ≤ −34.
Our results in this very delicate and long quest have been obtained from two independent time series, as well as two independent frequency ranges. This lends them great credence. They open exciting new questions about the gmode amplitudes, their excitation mechanism in such a broad and lowfrequency range, and the solar rotation evolution. They provide strong new constraints for a dramatically improved description of the solar core.
Acknowledgments
The GOLF instrument on board SOHO is a cooperative effort of scientists, engineers, and technicians, to whom we are indebted. For more details, see Gabriel et al. (1995) or http://www.ias.upsud.fr/golf/. GOLF data are available at MEDOC data and operations centre (CNES/CNRS/Univ. ParisSud), http://medoc.ias.upsud.fr/. They are available as well on the GOLF web site. The data file used in this analysis has been calibrated by R. García. The high quality of the GOLF data is due equally to the outstanding performance of the SOHO platform and overall system. SOHO is a project of international collaboration between ESA and NASA. We would like to acknowledge the support received continuously during more than three decades from CNES. D.S. acknowledges the financial support from the CNES GOLF grant and the Observatoire de la Côte d’Azur for support during his stays. R.K.U. acknowledges support from NASA for his participation in this project and thanks John Bahcall for enthusiastic encouragement for the gmode search. We also wish to acknowledge R. M. Bonnet for his leadership in making the SOHO mission and GOLF a reality and P. Delache, who stimulated and motivated many of us in his broadly acknowledged role of gmode ambassador. We wish to thank the referee for the initial suggestions that definitely helped to improve the general structure and to make the presentation more convincing.
References
 Appourchaux, T., & Pallé, P. 2013, The History of the mode Quest, in Fifty Years of Seismology of the Sun and Stars, eds. K. Jain, S.C. Tripathy, F. Hill, J. W. Leibacher, & A. A. Pevtsov (San Francisco: ASP), ASP Conf. Ser., 478, 125 [Google Scholar]
 Appourchaux, T., Gizon, L., & RabelloSoares, M.C. 1998, A&ARv, 132, 107 [Google Scholar]
 Appourchaux, T., Belkacem, K., Broomhall, A.M., et al. 2010, A&ARv, 18, 197. [NASA ADS] [CrossRef] [Google Scholar]
 Berthomieu, G., & Provost, J. 1991, Sol. Phys., 133, 127 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J., & Berthomieu, G. 1991, Theory of solar oscillations in Solar Interior and atmosphere, eds. A. N. Cox, W. C. Livingston, & M. S. Matthews (Tucson: University of Arizona Press), Space Sci. Ser., 401 [Google Scholar]
 Claverie, A., Isaak, G., McLeod, C., van der Raay, H., & Roca Cortés, T. 1979, Nature, 282, 591 [NASA ADS] [CrossRef] [Google Scholar]
 Duvall, T. L. Jr. 2004, A New Method to Search for Solar Gravity Mode Oscillations, in SOHO 14/ GONG 2004 Helio and Asteroseismology: Towards a Golden Future, ed. D. Danesy (Noordwijk: ESA Publications Division), ESA SP, 559, 412 [Google Scholar]
 Fierry Fraillon, D., Gelly, B., Schmider, F. X., et al. 1998, A&A, 333, 362 [NASA ADS] [Google Scholar]
 Fossat, E., Grec., G., Gavryusev, V., & Gavryuseva, E. 1988, Searching for solar modes, Tests of a statistical method, in Seismology of the Sun and SunLike Stars, ed. E. J. Rolfe (Noordwijk: ESA Publications Division), ESA SP, 286, 393 [Google Scholar]
 Fossat, E., Kholikov, S., Gelly, B., et al. 1999, A&A, 343, 608 [NASA ADS] [Google Scholar]
 Fossat, E., Salabert, D., Cacciani, et al. 2002, Eleven years of IRIS frequencies and splittings, in Local and global helioseismology: the present and future, ed. H. SawayaLacoste (Noordwijk: ESA Publications Division), ESA SP, 517, 139 [Google Scholar]
 Gabriel, A. H., Grec, G., Charra, J., et al. 1995, Sol. Phys., 162, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Gabriel, A. H., Baudin, F., Boumier, P., et al. 2002, A&A, 390, 1119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García, R. A., TurckChièze, S., Boumier, P., et al. 2005, A&A, 442, 385 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García, R. A., TurckChièze, S., JiménezReyes, S., et al. 2007, Science, 316, 1591 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 García, R. A., Salabert, D., Ballot, J., et al. 2011, J. Phys. Conf. Ser., 271, 012046 [NASA ADS] [CrossRef] [Google Scholar]
 Grec, G., Fossat, E., & Pomerantz, M. 1980, Nature, 288, 541 [NASA ADS] [CrossRef] [Google Scholar]
 Grec, G., Fossat, E., & Pomerantz, M. 1983, Sol. Phys., 82, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Kennedy, J. R., Jefferies, S. M., & Hill, F. 1993, Solar GMode Signatures in PMode Signals, in GONG 1992, Seismic Investigation of the Sun and Stars, ed. T. M. Brown (San Francisco: ASP), ASP Conf. Ser., 42, 273 [Google Scholar]
 Komm, R., Howe, R., Durney, B. R., & Hill, F. 2003, ApJ, 586, 650 [NASA ADS] [CrossRef] [Google Scholar]
 Lou, Y. Q. 2001, ApJ, 556, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Mathur, S., TurckChièze, S., Couvidat, S., & García, R. A. 2007, ApJ, 668, 594 [NASA ADS] [CrossRef] [Google Scholar]
 Percival, D. B. 1992, Comp. Sci. Stat., 24, 534 [Google Scholar]
 Provost, J., & Berthomieu, G. 1986, A&A, 165, 218 [NASA ADS] [Google Scholar]
 Provost, J., Berthomieu, G., & Morel, P. 2000, A&A, 353, 775 [NASA ADS] [Google Scholar]
 Salabert, D., JiménezReyes, S. J., Fossat, E., Gelly, B., & Schmider, F. X. 2002, Variability of mode parameters in 11 years of IRIS++ data, in Proc. SOHO 11 Symp., From Solar Min to Max: Half a Solar Cycle with SOHO, ed. A. Wilson (Noordwijk: ESA Publications Division), ESA SP, 508, 95 [Google Scholar]
 Salabert, D., Leibacher, J., Appourchaux, T., & Hill, F. 2009, ApJ, 696, 653 [NASA ADS] [CrossRef] [Google Scholar]
 Severny, A. B., Kotov, V. A., & Tsap, T. T. 1976, Nature, 259, 87 [NASA ADS] [CrossRef] [Google Scholar]
 TurckChièze, S., & Lopes, I. 2012, RA&A, 12, 8 [Google Scholar]
 TurckChièze, S., García, R. A., Couvidat, S., et al. 2004, ApJ, 604, 455 [NASA ADS] [CrossRef] [Google Scholar]
 TurckChièze, S., Palacios, A., Marques, J. P., & Nghiem, P. A. P. 2010, ApJ, 715, 1539 [NASA ADS] [CrossRef] [Google Scholar]
 TurckChièze, S., García, R. A., Lopes, I., et al. 2012, ApJ, 746, L12 [NASA ADS] [CrossRef] [Google Scholar]
 Welch, P. D. 1967, IEEE Transactions on Audio and Electroacoustics, 15, 70 [CrossRef] [Google Scholar]
All Figures
Fig. 1 Power spectrum of 16.5 years of fulldisk velocity recorded by the GOLF instrument on board the SOHO spacecraft. 

In the text 
Fig. 2 From the same 16.5year data set as Fig. 1, this is an average of 34 612 spectra computed from 8h selections taken at intervals of 4 h. 

In the text 
Fig. 3 Example of one of the 34612 GOLF power spectra, limited to the range 2.32−3.74 mHz and divided by the envelope of the mean spectrum of Fig. 2. 

In the text 
Fig. 4 Power spectrum of the spectrum shown in Fig. 3. The abscissa is in seconds and shows a main peak at a time T of about 14 800 s, which is the round trip travel time of an acoustic wave through the solar center. 

In the text 
Fig. 5 Sum of 34 612 power spectra similar to the spectrum shown in Fig. 4. 

In the text 
Fig. 6 16.5year time series of 36 131 values of the time T fluctuations defined in Fig. 4. The number of nonzero values is 34 261. 

In the text 
Fig. 7 Power spectrum P_{S}(ν) of the time series shown in Fig. 6. 

In the text 
Fig. 8 Enlargement of the lowest frequency range of Fig. 7. The arrows indicate the frequencies of the solar cycle and of the oneyear orbit. 

In the text 
Fig. 9 Estimate of the amplitude spectrum (square root of the density power spectrum) in blue, superposed on the modulus of the Fourier transform of the GOLF time series (in magenta). This estimate was used to calculate simulated time series. 

In the text 
Fig. 10 Autocorrelation function A(ν) of the frequency range 5.5−34.7 μHz of the power spectrum shown in Fig. 7. The first value displayed is bin number 7, i.e., 0.013 μHz. The vertical lines are at 210, 630, and 1260 nHz. 

In the text 
Fig. 11 Autocorrelation function of the frequency range 5.5−13.2 μHz of the power spectrum shown in Fig. 7. 

In the text 
Fig. 12 Autocorrelation function of the frequency range 13.2−34.7 μHz of the power spectrum shown in Fig. 7. 

In the text 
Fig. 13 S/N of the 210 nHz peak of Fig. 10 as a function of the duration of the time series, from 3 to 16.5 years. 

In the text 
Fig. 14 Autocorrelation function of the frequency range 5.5−34.7 μHz of a power spectrum resulting from an artificial data set. No peak is visible at 210 nHz, and the highest peak stands at 3.3σ to be compared to the 4.7σ of the peak obtained from the GOLF data set. 

In the text 
Fig. 15 Histogram of ten figures similar to Fig. 10 obtained with artificial data, simulated as described in text. 

In the text 
Fig. 16 Sum S_{cor} (Eq. (8)) of the three possible gmode splitting signatures in the autocorrelation, as a function of the gmode kernel rotation Ω_{g}. It identifies a fast rotation rate of 1277 ± 10 nHz. 

In the text 
Fig. 17 S/N of the main peak of Fig. 16 as a function of the duration of the time series analyzed, from 3 to 16.5 years. 

In the text 
Fig. 18 Same as Fig. 16, but computed separately for the unique contribution of l = 1, A(s(1,1)) in blue, or for the contribution of two components of l = 2, A(s(2,1)) + A(s(2,2)) in red. This illustrates that the l = 2 contribution confirms one of the two possible rotation rates detected by l = 1 with ambiguity. 

In the text 
Fig. 19 Same as Fig. 18 for an artificial data set. 

In the text 
Fig. 20 Sum C_{1}(P_{1},P_{min,1}) (Eq. (12)) of the 228 values of a model based on 76 gmode dipole triplets with the splitting value s_{1} = 210 nHz, applied to the power spectrum of Fig. 7, smoothed over 6 bins. The two free parameters of the model are the equidistant period spacing P_{1}, defining the abscissa, and the period P_{min,1} of the first selected mode, defining the ordinate axis. The first (lowest) sharp peak is located at 30 360 ± 15 s (FWHM) and covers the radial order range from −20 to −95, the second (upper part of the figure) covers the radial order range from −21 to −96. The equidistant period spacing is 1443 s, or 24 min 03 s, within less than 1 s. 

In the text 
Fig. 21 Each value is the maximum of a row of Fig. 20, so that the abscissa is the ordinate of Fig. 20. The interval between these peaks is the l = 1 equidistant period spacing of 1443 s. The power scale is the mean power per individual gmode peak. 

In the text 
Fig. 22 As in Fig. 20, the sum of the 228 values of a model based on 76 gmode dipole triplets and applied on an artificial data set, with an assumed splitting value of 267 nHz that corresponds to the highest correlation peak in Fig. 14. No significant equidistant period is shown, as demonstrated quantitatively by the next figure. 

In the text 
Fig. 23 Each value of this figure is the maximum of a row of Fig. 22, so that its abscissa is the ordinate of Fig. 22. This figure relative to an artificial data set must be compared to Fig. 21 obtained with the real GOLF data. 

In the text 
Fig. 24 Sum of the 76 values based on the same model based on 76 g modes; this model is simplified to the 76 central gmode frequencies, thus ignoring knowledge of their splitting. 

In the text 
Fig. 25 Mean profile M(ν) of 76 equidistant (in period) dipole modes whose central frequencies are defined by Fig. 20, shown for a 1.2 μHz bandwidth (Eq. (13)). The triplet splitting of this profile, even when contaminated by the background noise and many overlapping neighbors, is nevertheless directly measured on the modes themselves, at 209.5 ± 2 nHz. 

In the text 
Fig. 26 Sum C_{2}(P_{2},P_{min,2}) (Eq. (12)) of the 560 components of 112 l = 2 quintuplets, with a single splitting value of s_{2} = 630 nHz, and the two free parameters are the period spacing P_{2} in the abscissa, in the range 825−840 s, and the first central period of the series P_{min,2} in the ordinate, as in Fig. 20. The equidistance is seen at 832.8 s. 

In the text 
Fig. 27 Vertical cut of Fig. 26 at 832.8 s, extended to a third period. It shows a good S/N and clearly locates the first period of the series at 29 380 ± 30 s. This peak is the mean of the 560 components of the 112 l = 2g modes in the range of radial orders from −33 to −144. The other two peaks cover the ranges (−34,−145) and (−35,−146). 

In the text 
Fig. 28 Same as Fig. 26 for the artificial data set. This cut is naturally random as no l = 2 splitting was identified with the simulated data set. Consequently, no equidistant period spacing can be identified. 

In the text 
Fig. 29 Assumed unknown pmode rotation rate tested against the optimum visibility of the three gmode splitting values that are assumed to exist, but are unknown so far. The wellknown 433 nHz pmode rotation rate is confirmed by the g modes, only a little less accurately than when measured directly on the p modes themselves. 

In the text 
Fig. 30 Asymptotic gmode rotational kernel K(r) (Eq. (4)). 

In the text 
Fig. 31 Same as Fig. 21, computed on the same sequence of 76 layers of the power spectrum, equidistant in periods. This sum is computed separately here for the independent first half (red) and second half (blue) of the 16.5year time series. The power scale of this and the next three figures is arbitrary, as they contain two curves that are rescaled to facilitate comparison. 

In the text 
Fig. 32 Following Figs. 21 and 31, this illustration uses two independent sequences of 38 intervals that are equidistant in periods, the first half (shorter periods, radial order range from −20 to −57 for the first of the four peaks) plotted in blue, and the second half (longer periods, radial order range from −58 to −95 for the first peak, rescaled by a factor 0.8) shown in red, both computed on the complete time series of 16.5 years. The first four periods of the sequences are shown. 

In the text 
Fig. 33 Figure to be compared to Fig. 27 for the l = 2g modes, computed for the same sequence of 112 layers of the power spectrum, equidistant in periods. This sum is computed here separately for the independent first half (red) and second half (blue) of the 16.5year time series. The first three periods of each sequence are shown. 

In the text 
Fig. 34 As in Fig. 32 for l = 1, this illustration for l = 2 uses two independent sequences of 56 intervals equidistant in period, the first half (shorter periods, ranging from n = −33 to −88) is plotted in blue and the second half (longer periods, ranging from n = −89 to −144, slightly rescaled) is shown in red, both computed on the complete time series of 16.5 years. The first three periods of each sequence are shown. 

In the text 
Fig. 35 Mean symmetrized profile of the 76 dipolar g modes, following the model that includes a slight departure from the asymptotic equidistance. The comparison with Fig. 25 shows that taking this slight departure into account improves the visibility. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.