Issue 
A&A
Volume 629, September 2019



Article Number  A26  
Number of page(s)  16  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201935434  
Published online  29 August 2019 
Signature of solar g modes in firstorder pmode frequency shifts
^{1}
MaxPlanckInstitut für Sonnensystemforschung, JustusvonLiebigWeg 3, 37077 Göttingen, Germany
email: boening@mps.mpg.de
^{2}
MaxPlanckInstitut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
^{3}
Institut für Astrophysik, GeorgAugustUniversität Göttingen, FriedrichHundPlatz 1, 37077 Göttingen, Germany
^{4}
Center for Space Science, NYUAD Institute, New York University Abu Dhabi, PO Box 129188, Abu Dhabi, UAE
Received:
8
March
2019
Accepted:
27
June
2019
Context. Solar gravity modes (g modes) are buoyancy waves that are trapped in the solar radiative zone and have been very difficult to detect at the surface. Solar g modes would complement solar pressure modes (p modes) in probing the central regions of the Sun, for example the rotation rate of the core.
Aims. A detection of g modes using changes in the large frequency separation of p modes has recently been reported. However, it is unclear how p and g modes interact. The aim of this study is to evaluate to what extent g modes can perturb the frequencies of p modes.
Methods. We computed the firstorder perturbation to global pmode frequencies due to a flow field and perturbations to solar structure (e.g. density and sound speed) caused by a g mode. We focused on longperiod g modes and assumed that the gmode perturbations are constant in time. The surface amplitude of g modes is assumed to be 1 mm s^{−1}, which is close to the observational limit set by Doppler observations.
Results. Gravity modes do perturb pmode frequencies to first order if the harmonic degree of the g mode is even and if its azimuthal order is zero. The effect is extremely small. For dipole and quadrupole p modes, all frequency shifts are smaller than 0.1 nHz, or 2 × 10^{−8} in relative numbers. This is because the relative perturbation to solar structure quantities caused by a g mode of realistic amplitude is of the order of 10^{−6}–10^{−5}. Additionally, we find that structural changes dominate over advection. Surprisingly, the interaction of g and p modes takes place to a large part near the surface, where p modes spend most of their propagation times and g modes generate the largest relative changes to solar structure. This is due to the steep density stratification, which compensates the evanescent behaviour of g modes in the convection zone.
Conclusions. It appears to be impossible to detect g modes solely through their signature in pmode frequency shifts. Whether g modes leave a detectable signature in pmode travel times under a given observational setup remains an open question.
Key words: Sun: helioseismology / Sun: interior / Sun: oscillations / waves
© V. G. A. Böning et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Open Access funding provided by Max Planck Society.
1. Introduction
Using observations of pressure modes (p modes), helioseismology has been very successful in determining solar interior structure and differential rotation in the convection zone and in most of the radiative zone (see, e.g. Howe 2009; Aerts et al. 2010; Basu 2016). However, a number of outstanding questions about the deep solar interior remain, such as the rotation of the solar core.
Solar internal gravity modes (g modes), if detected, would provide complementary information to the p modes. This is because g modes have most of their kinetic energy in the central regions of the Sun. In addition, they have a smaller radial wavelength in this region (e.g. Berthomieu & Provost 1990, 1991), which would give a higher spatial resolution than the p modes.
The quest for detecting solar g modes has, therefore, attracted considerable attention in the past with several reported detections (see Appourchaux et al. 2010; Appourchaux & Pallé 2013, for reviews). Unfortunately, none of the reported detections have been confirmed so far.
In line with Kennedy et al. (1993) and Duvall et al. (2004), Fossat et al. (2017) proposed a new method for detecting g modes in the firstarrival travel time of lowdegree p modes using diskintegrated observations. Under the simplest interpretation, this pmode travel time is inversely proportional to the large frequency separation of p modes, which is the frequency difference between p modes with identical harmonic degrees whose radial orders differ by one. The pmode travel time may be affected by a gmode oscillation through changes in sound speed, for example. Fossat et al. (2017) measured the pmode travel time using eighthour intervals of Global Oscillations at Low Frequency (GOLF; Gabriel et al. 1995) data and a sampling rate of four hours. The idea is that such a time series can sample changes in solar structure and dynamics due to highorder g modes that have periods longer than eight hours. Using this method, Fossat et al. (2017) and Fossat & Schmider (2018) reported periodic signals in the power spectrum of the data that were interpreted in terms of the rotational splittings of the g modes, implying a solar core rotating nearly four times faster than the radiative zone. Subsequently, Schunker et al. (2018) reproduced the signal of Fossat et al. (2017), but warned that it is very sensitive to a number of choices and parameters in the analysis procedure. Very recently, Scherrer & Gough (2019) also showed that the statistical significance of the signal depends on the helioseismic instrument used (from SOHO, SDO, or from the ground).
The motivation for the present study is to understand to what extent g modes can affect pmode frequencies and whether and how such an effect would be measurable. We do not address the question of the signature of g modes in the observed pmode travel times, which will require a separate study in order to include the details of the measurement procedure.
One may first ask the general question whether p and g modes should interact at all. In a linear wave equation, the superposition principle holds. It ensures that any two wave fields can be added and still fulfil the wave equation. As a consequence, the presence of one wave does not change the properties of a second wave at all. The situation is different, however, when nonlinear effects are taken into account. Nonlinearities in the wave equation include nonlinearities in the Navier–Stokes equations. For example, this is the case for advection, temporal changes in the background (sound speed, density, pressure), dissipation, and nonadiabatic effects (e.g. Chaps. 2 and 3 in Hamilton & Blackstock 1998). In a nonlinear wave equation, two solutions to the wave equation do not, in general, add up to a third solution because the waves interact. This effect may in principle be exploited. A successful example of detecting the presence of waves in the Sun using other waves is given by the recent detection of solar equatorial Rossby waves, which were detected because they advect p modes (Löptien et al. 2018; Liang et al. 2019; Hanasoge & Mandal 2019).
In classical acoustics, the nonlinear scattering of sound by sound has been observed experimentally (e.g. Thuras et al. 1935) and modelled theoretically (e.g. Westervelt 1957, 1963; Hamilton & Blackstock 1998, and references therein). For example, Westervelt (1957) obtained an approximative nonlinear acoustic wave equation by inserting a solution to the linear wave equation into the full nonlinear Navier–Stokes equations. In addition, the nonlinearity introduced by a second order term in the equation of state was taken into account. In the presence of sound waves of two different frequencies, the generation of waves with frequencies equal to the sum and difference of the original frequencies is observed (e.g. Thuras et al. 1935).
The nonlinear coupling of stellar oscillation eigenmodes has been studied by several authors. Perdang & Blacher (1982) used a Hamiltonian formalism and showed that a nonlinear coupling can lead to chaotic motions with an approximate periodicity. Dziembowski (1982) studied nonlinear interactions of two and three modes and found that twomode interaction has an amplitudelimiting effect on stellar pulsations (see also Kumar & Goldreich 1989). Based on Dziembowski (1982), Wentzel (1987) showed that two p modes can generate a g mode by a nonlinear interaction. Wentzel (1987) argued that the evanescent behaviour of the g modes in the convection zone is overcompensated by the increase of pmode sensitivity towards the surface. As a consequence, the coupling would take place predominantly in the convection zone. The nonlinear coupling of p and g modes was also studied for tidally interacting neutron star binaries by Weinberg et al. (2013), where the coupling is found to alter the phase of the emitted gravitational wave and thereby its detection signal. Dziembowski (1982) considered modes that are in a resonance condition, that is one frequency is the sum or difference of the other two frequencies. In Weinberg et al. (2013), the mode frequencies are nonresonant, but the radial wave numbers of the p and g modes are similar in a region in the stellar interior, which means that the radial wave numbers fulfil a resonance condition.
A possible signature of g modes in the pmode frequency signal has been studied previously by Kennedy et al. (1993), Lou (2001), and Scherrer & Gough (2019). In these studies, the gmode period was assumed to be much longer than the pmode period, and hence the perturbation to solar structure caused by a g mode was assumed to be timeindependent on shorter temporal scales. Kennedy et al. (1993) used perturbation theory on the asymptotic pmode dispersion relation and found that g modes can produce a frequency modulation of the p modes. This would create sidelobes near each pmode frequency, which was also found by Lou (2001). Using degenerate perturbation theory, Scherrer & Gough (2019) gave an expression for the firstorder perturbation of p modes by g modes and estimated the order of magnitude of the signal. They found that the signal reported by Fossat & Schmider (2018) would be due to relative frequency shifts of the order of 10^{−5}, which would correspond to a gmode amplitude of 50 cm s^{−1} and which is much higher than the observational upper limit on gmode amplitudes (for a review, see Appourchaux et al. 2010). Both Kennedy et al. (1993) and Scherrer & Gough (2019) found that only g modes with even harmonic degrees can couple to p modes to leading order.
As the nonlinearity in the wave equation is certainly a small effect, the interaction of p modes with g modes can be studied using perturbation theory. The simplest approach is to model the effect of a timeindependent perturbation on mode frequencies (Kennedy et al. 1993; Scherrer & Gough 2019). In addition, it is possible to consider the change in eigenfunctions (e.g. Lavely & Ritzwoller 1992), changes in mode amplitude ratios (e.g. Schad et al. 2011, 2013), or the coupling of oscillations at different frequencies (e.g. Woodard 2007, 2016; Hanasoge et al. 2017).
In this work, we expand on the approach taken by Kennedy et al. (1993) and Scherrer & Gough (2019), and study the firstorder perturbation of pmode frequencies by a frozen g mode, that is we assume that the perturbation caused by the g mode is constant in time. Our aim is to obtain quantitative estimates of the induced frequency shifts, and to understand in a more detailed manner where the interaction is taking place and which of the perturbative effects of the g mode is causing the interaction.
The effect of aspherical perturbations on pmode frequencies has been studied by Gough (1993), on which Kennedy et al. (1993) and Scherrer & Gough (2019) are based. The aforementioned work was done under the assumption that the perturbation to solar structure preserves hydrostatic equilibrium. In this paper, we include the effect of a departure from hydrostatic equilibrium as it is the case for a g mode. In addition, we provide detailed quantitative estimates of the induced firstorder frequency shift and we provide a detailed analysis of the resulting selection rules. To do that, we will include a number of effects such as changes to solar structure through changes in sound speed, pressure, density, gravitational potential, and advection due to the flow field from a g mode.
In addition, we need numerical solar p and g mode eigenfunctions which can be computed using stellar oscillation codes like ADIPLS (ChristensenDalsgaard 2008) or GYRE (Townsend & Teitler 2013). Such a computation is based on a stellar model. In this work, we use model S (ChristensenDalsgaard et al. 1996) and a result from the MESA stellar evolution code (Paxton et al. 2011).
This paper is structured as follows. We first introduce the oscillation equations (Sect. 2) and the changes in structure and dynamics due to a g mode (Sect. 3). In Sect. 4 we setup the perturbation theory and we derive the equations for the firstorder frequency shift of a p mode in the presence of a g mode, including selection rules. We then present numerical results in Sect. 5 and conclusions in Sect. 6.
2. Linearised equations of stellar oscillations
To model the effect of solar g on p modes, we briefly review the equations of stellar oscillations following Aerts et al. (2010, Chap. 3) and Basu (2016). Solar g and p modes are characterised by an eigenfunction ξ(r) which describes the displacement of a parcel of gas during an oscillation, depending on the location r in the solar interior. The displacement eigenfunction is a solution to the linear wave equation
where ω is the cyclic eigenfrequency of the mode and the wave operator ℒ can be written
Here, ρ, p, and g denote density, pressure, and gravitational acceleration, and u_{0} is a background flow, which we assume to be equal to differential rotation (e.g. Schou et al. 1998),
Quantities with a zero subscript are values in hydrostatic equilibrium of the Sun. Quantities with a one as subscript denote Eulerian firstorder perturbations to hydrostatic equilibrium and their explicit dependence on the displacement is given in Appendix A.
It can be shown (LyndenBell & Ostriker 1967) that the wave operator ℒ is selfadjoint under the scalar product
for a suitable boundary condition, for example a vacuum boundary with ∇ ⋅ ξ = 0. As a consequence, two eigenfunctions ξ and ξ′ with eigenvalues ω^{2} ≠ ω′^{2} are orthogonal with regard to the above scalar product. For comparison purposes, we use four different boundary conditions in this work, see Sect. 3 and Appendix G.
In hydrostatic equilibrium and in absence of a background flow (u_{0} = 0), the Sun is assumed to be a spherically symmetric star. In this case, the eigenfunctions can be written in terms of spherical harmonics,
where Ψ_{lm} is a vector spherical harmonic (Barrera et al. 1985), the quantities are unit vectors in the directions of r, θ, ϕ, the functions ξ_{r}, ξ_{h} are radial and horizontal displacement eigenfunctions, and the quantum numbers l, m, n denote harmonic degree, azimuthal order, and radial order, respectively. We adopt the usual convention that the radial order n for a p mode is positive and for a g mode is negative (e.g. Aerts et al. 2010, Chap. 3) and normalise the eigenfunctions so that
The linearity of the wave equation implies that any mode, for example a p mode, is not affected to first order by the restoring force of a different mode, for example a g mode. By first order, we mean here the linearization of the Navier–Stokes equations that yields the linear wave equation. However, as g modes do perturb the hydrostatic equilibrium of the Sun, they may have some affect on p modes at a higher order, if nonlinear effects are taken into account.
3. Perturbation to solar structure and flows due to g modes
In order to determine the effect of a g mode on pmode frequencies, we need to determine the perturbation of a g mode to solar structure. To do so, we computed eigenfunctions for g and p modes using the codes GYRE (Townsend & Teitler 2013; Townsend et al. 2018) and ADIPLS (ChristensenDalsgaard 2008). The computations are based on three nonrotating solar equilibrium models, namely solar model S (ChristensenDalsgaard et al. 1996), a version of model S where the density was smoothed slightly, and a solar model obtained from MESA (Paxton et al. 2011), see Appendix G. Furthermore, we tested the effect of several surface boundary conditions like zero Lagrangian pressure perturbation (vacuum boundary) at the surface and two conditions for an isothermal atmosphere (ChristensenDalsgaard 2008; Unno et al. 1989), as well as the atmospheric boundary condition from Dziembowski (1971). The Eulerian perturbations to pressure, density, and gravitational potential for a particular eigenmode are a direct output from the GYRE code while they also can be obtained from the displacement eigenfunctions directly (GYRE and ADIPLS).
As an example, we show eigenfunctions for l_{g} = 2 g modes with n_{g} = −10 and n_{g} = −36 in Figs. 1 and 2. The eigenfunctions were computed using the slightly smoothed version of model S, the GYRE code, and the boundary condition of an isothermal atmosphere from ChristensenDalsgaard (2008). The panels in the left columns of Figs. 1 and 2 show that the kinetic energy density of the g modes is concentrated in the radiative zone. The g modes are trapped in this region, the eigenfunctions are oscillatory and the absolute changes to solar structure variables caused by the g mode, i.e. Eulerian perturbations to density, pressure, and squared sound speed c^{2}, reach their maximum amplitude. In the convection zone, the amplitude of the kinetic energy density is lower and decays. The mode is evanescent in this region. This is well known (e.g. ChristensenDalsgaard 2003).
Fig. 1. Eigenfunctions and Eulerian perturbations to solar structure for a typical g mode with l_{g} = 2, n_{g} = −10, ν_{g} = 103.2 μHz (compare to Fig. 5.10 in ChristensenDalsgaard 2003). The mode was normalised to a surface amplitude of 1 mm s^{−1}, which is at the order of observational upper limits (e.g. Appourchaux et al. 2010). Displayed are (from top to bottom) radial and horizontal eigenfunctions, and the Eulerian perturbations to density, pressure, and squared sound speed. Left column: kinetic energy density and absolute changes to solar structure caused by the mode are centred near the solar core. Right column: relative changes to solar structure are very small, and in the case of density, pressure, and sound speed, reach the largest values near the surface. The exact behaviour of the relative Eulerian changes to solar structure at the surface depends on the choice of boundary condition used in the eigenfunction computation. 

Open with DEXTER 
Fig. 2. Same as Fig. 1, but for a g mode with n_{g} = −36, ν_{g} = 31.2 μHz. See also Figs. B.4 and B.5. 

Open with DEXTER 
On the other hand, the panels in the right columns of Figs. 1 and 2 show that relative changes to density, pressure, and sound speed are very small, namely at the order of 10^{−6} − 10^{−5} for a g mode normalised to a surface amplitude of 1 mm s^{−1}. This is close to the upper limit for the gmode amplitude from observations (e.g. Appourchaux et al. 2010, and references therein). Even more interestingly, the relative changes to solar structure caused by the g mode reach local maxima close to the solar surface and their amplitude is larger than in the core, at least for g modes of lower radial order. For modes with high radial order, the relative changes to solar structure are larger but of similar order of magnitude in the core compared to the surface. This concentration of relative changes to solar structure near the surface is due to the sharp decrease of density and pressure in this region and takes place despite the concentration of kinetic energy near the core and despite the evanescent character of the g modes in the convection zone. To guarantee the correctness of the eigenfunction computation, we have performed a number of tests, see Appendix B. These tests include a comparison of eigenfunctions to asymptotic formulae, see Figs. B.4 and B.5.
The question is now, how p modes are sensitive to the perturbation to solar structure from a g mode. It is well known that p modes are trapped between their lower turning point, which mostly depends on their horizontal surface phase speed , and their upper turning point close to the photosphere (e.g. ChristensenDalsgaard 2003). The kinetic energy of p modes is concentrated in this region and it is increasing towards the surface. In addition, p modes spend most of their lifetime close to the surface due to the strong increase of the sound speed in the solar interior.
As a consequence, the sensitivity of p modes to relative changes in solar structure increases towards the surface (e.g. Figs. 10 and 11 in Basu 2016). At the same time, relative changes to solar structure due to g modes reach considerable amplitudes near the surface as we have illustrated in Figs. 1 and 2.
We thus conclude that the coupling of a g mode with a p mode takes place to a large extent near the solar surface, if there is any such coupling present in the Sun. The hypothesis that the increase of pmode sensitivity in the convection zone would overcompensate the decrease in amplitude of the g mode was already put forward by Wentzel (1987).
4. pmode frequency perturbation due to a frozen g mode
In the following, we derive equations for the firstorder perturbation to a pmode eigenfrequency in the presence of a g mode. To do that, we assume that the period of the g mode is much longer than that of the p mode. Furthermore, we assume that the gmode period is long compared to the observational time scale. To simplify the analysis, we therefore assume that the perturbation to solar structure caused by a g mode is constant in time (see also Kennedy et al. 1993; Scherrer & Gough 2019). In this picture, the timedependence of the g mode may be added after the frequency shift has been calculated.
The perturbed wave operator for changes to solar structure has been derived before (see, e.g. Basu 2016, and references therein). However, these derivations were made for a perturbation to solar structure that preserves hydrostatic equilibrium. As the perturbation in our case is due to a gmode oscillation, this assumption is no longer valid. We therefore rederive the perturbed wave operator for a perturbation to solar structure that does not necessarily conserve hydrostatic equilibrium, see Appendix C.
The perturbed wave operator ℒ + δℒ acts on a pmode with eigenfunction ξ_{p} and frequency ω_{p}, where the index p stands for a p mode with
On the other hand, the perturbation is caused by a g mode with eigenfunction ξ_{g} and frequency ω_{g} with index
Consequently, we can express changes to all structure quantities, namely density, pressure, gravitational acceleration, and squared soundspeed, as well as the flow field associated to the gmode oscillation in terms of the gmode eigenfunction or Eulerian perturbations (see also Eqs. (A.1)–(A.3)),
where Γ_{1, 0}(r) is the first adiabatic exponent of the equilibrium solar model (ChristensenDalsgaard 2008; Aerts et al. 2010). In Eq. (13), we have assumed that the Lagrangian perturbation to the first adiabatic exponent during the gmode oscillation is zero. The adiabatic exponent is a material property and reflects the chemical composition and the ionization state of the plasma (see, e.g. Kippenhahn & Weigert 1990, Chap. 14). During the oscillation, the chemical composition does not change when moving with the plasma. The degree of ionization, however, is a function of the gas pressure and temperature of the plasma. As temperature and pressure do change during the oscillation, the degree of ionization can change as well. This may have a small effect on the adiabatic exponent in the ionization zones of the Sun, which we neglect in this analysis.
4.1. Frequency shifts in the rotating Sun
The pmode eigenfunctions of the wave equation (Eq. (1)) for the rotating Sun can be expressed in terms of the solutions of the eigenfunctions in the nonrotating case, ξ_{p}, to first order in quasidegenerate perturbation theory (Lavely & Ritzwoller 1992), according to
where we note that in Eq. (16) only terms of the same m as in ξ_{p} enter. The perturbation to the eigenfunction, , has a much smaller amplitude than the unperturbed eigenfunction if the harmonic degree l_{p} of the unperturbed mode is small (compare to Woodard 1989). Similarly, the perturbed eigenfrequencies can be written, to first order, as (e.g. Ritzwoller & Lavely 1991)
We assume that rotation does not couple p modes and g modes since their frequencies differ largely.
For g modes, the same expansion as in Eqs. (15)–(18) is possible, and we note that the perturbation to the eigenfunction, , has a much smaller amplitude than the eigenfunction of the nonrotating model also for g modes in the case of small harmonic degrees as considered in this paper.
Given the eigenfunctions and mode frequencies for the rotating solar model, we can now evaluate the frequency perturbation of the p mode due to a g mode using nondegenerate perturbation theory. The use of nondegenerate perturbation theory will be justified later because the frequency shift from the g modes is much smaller than any frequency difference between the eigenfrequencies of the rotating eigenmodes. In this case, the frequency shift from a g mode is simply given by
where indicates that the perturbation is caused by a g mode in a rotating Sun, . Using this and (see Eqs. (6) and (15)), Eq. (19) can be written to leading order as
Here, all terms involving or are of higher than first order and have therefore been omitted. Hence, we evaluate δ(ω_{p}^{2}) and the frequency shift as 2π δν_{p} = δω_{p} = δ(ω_{p}^{2})/(2ω_{p}) in the following.
The scalar product in Eq. (20) leads to a number of integrals over one gmode and two pmode eigenfunctions, which are given in detail in Appendix D for the terms involved in δℒ. Similar integrals have been found by other authors for the coupling of g and p modes (e.g. Dziembowski 1982; Kennedy et al. 1993; Scherrer & Gough 2019), as well as for flows and aspherical perturbations to solar structure preserving hydrostatic equilibrium that are decomposed in spherical harmonics (e.g. Lavely & Ritzwoller 1992; Gough 1993). We note here that we omitted one term (δℒ_{4}) for simplicity. This term models the change of the Eulerian perturbation to the gravitational potential during the p mode oscillation, or in other words changes to Φ_{1, p}, in the presence of the g mode. This term is expected to have a small contribution to the wave equation.
4.2. Selection rules
Our derivation of the firstorder frequency shift gives rise to selection rules for mode coupling, which are derived in Appendix E. The horizontal integrals over three spherical harmonics appearing in Appendix D can be evaluated using Gaunt’s formula (e.g., Dahlen & Tromp 1998, see also Hanasoge et al. 2017; Kiefer et al. 2017; Fournier et al. 2018), see Appendices E and F.
We find that a firstorder frequency shift of a p mode by a g mode is possible only if all of the following conditions are satisfied.

(I.)
The harmonic degree l_{g} of the g mode is even.

(II.)
The azimuthal order m_{g} of the g mode is zero.

(III.)
The harmonic degree l_{g} of the g mode is not larger than twice the harmonic degree of the p mode (l_{g} ≤ 2l_{p}).
Kennedy et al. (1993) and Scherrer & Gough (2019) have found previously that only g modes of even harmonic degree produce a pmode frequency shift at first order, see also Gough (1993) for perturbations to solar structure that preserve hydrostatic equilibrium. Although Scherrer & Gough (2019) do not note a selection rule regarding m_{g}, it can be derived using their framework of degenerate perturbation theory, where both the g mode and rotation where treated as perturbations to the solar model.
This can be seen by considering the coupling matrix that leads to the perturbed pmode eigenfunctions and perturbed pmode frequencies at first order. If m_{g} ≠ 0, the diagonal of that matrix is equal to the frequency perturbation from rotation. The coupling due to g modes appears at the offdiagonal matrix elements where m_{p} − m_{p′} = m_{g}. In the notation of Scherrer & Gough (2019), this corresponds to μ − μ′=m, see Eqs. (25)–(27) in their work. All other matrix elements are equal to zero. As the offdiagonal matrix elements due to the g mode are much smaller than the diagonal matrix elements due to rotation, the offdiagonal matrix elements only lead to a second order perturbation in the eigenfrequency and the firstorder frequency shift from the g mode is zero.
Somewhat similarly, the selection rules III can be obtained using the framework of quasidegenerate perturbation theory outlined in Lavely & Ritzwoller (1992, see also Schad 2013).
Here, we have modelled the firstorder frequency shift in the p modes but we did not model the change in pmode travel time as measured by Fossat et al. (2017). Given the above selection rules, global pmode frequency shifts can be ruled out as a cause for the observations of Fossat et al. (2017). However, in our eyes, it is not straightforward to make a direct connection between frequency shifts and the measurement of the roundtrip travel time by Fossat et al. (2017). For example, it is in principle possible that perturbations to the pmode eigenfunctions cause a signature of g modes in this measurement.
5. Numerical results
In the following, we numerically evaluate the analytical expression for the firstorder perturbation to the p mode eigenfrequency that we derived in the preceding section. Here, we present results for g modes with harmonic degree l_{g} = 2 and p modes with l_{p} = 2. Results for l_{p} = 1 are very similar. The g mode eigenfunctions are normalised to a surface amplitude of 1 mm s^{−1}, which is at the order of observational upper limits for the amplitude of solar g modes (e.g. Appourchaux et al. 2010). The frequency shifts computed in this section can therefore either be seen as approximative upper limits, or they can be rescaled to different gmode surface amplitudes.
Inspired by the analysis of Fossat et al. (2017), we assume that the observational time scale is eight hours and therefore apply our results to g modes with periods longer than that, which means frequencies ν_{g} < 34.7 μHz. Moreover, we obtain results also for loworder g modes with frequencies above 34.7 μHz for information purposes, even though our assumptions of a clear separation of observational time scale and gmode period may break down in this range. We do not consider results obtained for n_{g} < −200 (ν_{g} < 5.67 μHz) in the following. We checked however, that this limitation does not effect the major conclusions presented here. In addition to the pmode range considered by Fossat et al. (2017), we obtain results for p modes with frequencies ν_{p} < 5 mHz. The results presented in this section were obtained from eigenmodes computed with GYRE and solar model S (ChristensenDalsgaard et al. 1996), and we assumed an isothermal atmosphere as surface boundary condition (ChristensenDalsgaard 2008), see also Appendix G.
5.1. The effect is tiny in magnitude
As a first result, we find that the magnitude of the frequency shift of a p mode caused by a frozen g mode is extremely small, see Fig. 3. All frequency shifts δν_{p} computed are found to be below 0.1 nHz. The relative frequency shifts δν_{p}/ν_{p} are all below 2 × 10^{−8}. These numbers are in very broad agreement with the estimates presented in Scherrer & Gough (2019), who reported an order of magnitude of 2 × 10^{−8} when converted to the gmode amplitude assumed in this work.
Fig. 3. Firstorder perturbations to the pmode frequency due to a gmode oscillation assumed to be frozen over the observation period. Results are shown for g modes and p modes with harmonic degrees l_{g} = l_{p} = 2 and azimuthal orders m_{g} = m_{p} = 0. Left: frequency shift as a function of radial order for mode frequencies ν_{p} < 5 mHz and ν_{g} > 5.67 μHz. The range of modes studied by Fossat et al. (2017), Fossat & Schmider (2018), and Schunker et al. (2018) is inside the dashed rectangle, extending outwards of the plotted region at n_{g} > 200. Right: cuts through the results shown in the dashed region in the left panel for a selection of gmode radial orders. 

Open with DEXTER 
Generally, we find the frequency shifts to be greatest for large pmode and large gmode frequencies. The study by Fossat et al. (2017) considered a range of g mode frequencies, which does not cover the g modes with the largest frequencies, see Fig. 3 (right panel and dashed region in left panel). For these modes, the p mode frequency shifts are smaller by about one order of magnitude and do not exceed 0.01 nHz, or relative shifts of 10^{−9}.
The observational upper limit on the gmode amplitude depends on mode frequency, see for instance Fig. 15 in Appourchaux et al. (2010). While it is around 1 mm s^{−1} for the g modes with the lowest radial orders, it is at the order of 1 cm s^{−1} for modes with moderately higher radial orders, for example ν_{g} ≈ 35 μHz, n_{g} ≈ −33. This is where the analysis of Fossat et al. (2017) is sensitive to. Even if we rescale the frequency shifts for all modes with ν_{g} < 35 μHz to an amplitude of 1 cm s^{−1}, the maximum relative frequency shift among all modes does not change. This can also be seen by inspecting the relevant region in Fig. 3.
5.2. Contributions from density, pressure, sound speed, and flow
The question naturally arises, which of the effects of the g mode, the induced changes to density, pressure, sound speed, gravitional acceleration, or its flow field, is the dominant contributor to the frequency shift. In order to answer this question, we sum the contributions from the individual terms derived in Appendix D for each of these quantities. The results are shown in Fig. 4. In addition, we also show the contribution of an example term that usually is of large amplitude (δℒ_{3−2}, see Eq. (D.21)), and which is due to changes in density. It can be seen that this term is usually about one order of magnitude larger than the total density contribution.
Fig. 4. Contribution to the frequency shift of p modes by changes to solar structure variables ρ, p, c^{2}, g, and the flow field u caused by a g mode in absolute values. The total frequency shift δν_{p} and the contribution from an example term (δℒ_{3−2}, see Eq. (D.21)) are also given for comparison. The azimuthal order of the modes is m_{p} = m_{g} = 0. 

Open with DEXTER 
Similarly, the total contribution from relevant structural changes, that is from density, pressure, and soundspeed, is again about an order of magnitude smaller than their individual contribution. This is because these terms tend to cancel each other to a large extent. The contribution from changes to gravity is negligible.
Even though the contribution of the flow is about two orders of magnitude smaller than from individual density terms, it still has a relevant contribution to the total frequency shift of order 10%. In some cases, this number is larger, if the total frequency shift is close to zero.
In general, however, the effect of structural changes dominates over the effect of a flow caused by the g mode.
5.3. Regions of interaction
We further investigate in which region in the Sun the interaction predominantly occurs. To do that, we rewrite the total integral leading to the frequency shift from Eq. (20) as one radial integral. In addition, we decompose this integrand in the contribution from changes to solar structure and from the flow, see Fig. 5. We find that the total integrand takes the largest values very close to the surface and that it is nearly identical to the integrand from the structure terms. This is in agreement with our finding from Sect. 3 that the relative changes to solar structure from a g mode are concentrated near the surface. However, the integrands change sign near the surface, and the total contribution from the nearsurface layers is thus smaller than the magnitude of the integrands suggest.
Fig. 5. Radial integrands for computing the frequency shift of a p mode with l_{p} = 2, m_{p} = 0, n_{p} = 20 caused by a g mode with l_{g} = 2, m_{g} = 0, n_{g} = −36. When integrating over radius, the total integrand yields the total frequency shift δν_{p} (both panels). The contributions to the total integrand resulting from the structural changes and the flow are also shown (right panel). The spike at the surface boundary is due to numerical difficulties in computing second order derivatives at the boundary. Therefore, the three outermost radial grid points were discarded. 

Open with DEXTER 
To investigate this further, we computed the contributions to the total integral from the convection zone and from the radiative zone, which are shown in Fig. 6. It can be seen that the contribution from the radiative zone is about an order of magnitude smaller than from the convection zone (n_{g} = −36 in Fig. 6). A similar conclusion was reached by Wentzel (1987) by the simple argument that the evanescent behaviour of the g modes in the convection zone is overcompensated by the increase in sensitivity of the p modes. However, the contribution from the radiative zone can reach a comparable magnitude as from the convection zone for higherorder modes, for example n_{g} ≤ −100, see Fig. 6.
Fig. 6. Contribution to the frequency shift from the convection zone and from the radiative zone as a function of pmode frequency for several typical g modes considered in this study. Results are shown for l_{p} = l_{g} = 2, m_{p} = m_{g} = 0. 

Open with DEXTER 
Following ChristensenDalsgaard & Thompson (1997), one may ask whether the above conclusions stay the same if we use Lagrangian instead of Eulerian perturbations to solar structure. These authors found that contributions to the frequency shift from, e.g, Eulerian changes to sound speed and density nearly cancel as it is the case in our study. If one uses Lagrangian perturbations, such a near cancellation does not occur as strongly in their computation.
We therefore recomputed our main results using Lagrangian perturbations to solar structure for a selected number of mode pairs as test cases. As expected, the resulting total frequency shifts agree with the ones from the Eulerian perturbations (successful for low radial order gmodes n_{g}≲40).
We further considered the example of quadrupolar modes with n_{g} = − 36 and n_{p} = 20 in more detail. For this pair, the radial integrand from the Lagrangian variables has the same order of magnitude as the Eulerian one and a qualitatively similar behaviour. Both assume maximal values in the same region close to the surface. Their actual functional dependence on depth is different, which is consistent with ChristensenDalsgaard & Thompson (1997). The contributions from the radiative zone and the convection zone are of comparable amplitude in the Lagrangian case, but of opposite sign. These contributions have a larger amplitude as in the Eulerian case, which leads to additional nearcancellation. This is found for all pairs of modes we considered. Thus, Lagrangian perturbations do not seem to us as a better choice for studying this problem.
For both Eulerian and Lagrangian variables, we find in general that the contribution from nearsurface layers is of similar magnitude as the contribution from the convection zone. We thus conclude that the interaction of g and p modes is taking place to a large part near the surface, whether solar structure perturbations are computed in Eulerian or Lagrangian variables.
5.4. Dependence on solar model and surface boundary conditions
Due to the large magnitude of the integrands near the surface, one may of course ask whether and how the numerical results depend on the surface boundary condition applied in the eigenfunction computation.
We therefore computed eigenfunctions using a number of different boundary conditions, which are implemented in the GYRE code and which are detailed in Appendix G.
In addition, we checked whether our results depend on the solar model used to compute the eigenfunctions. The second derivative of the equilibrium density enters in our results, although it is not quite smooth in the original version of model S. This leads to spikes in some of the integrands near the core and wiggles near the surface, see for example the right panel in Fig. 5. We therefore also computed eigenfunctions using a version of model S, where the density was slightly smoothed, see details in Appendix G. In order to get a rough idea about the magnitude of possible systematic errors, we also obtained results for a 1 M_{⊙} stellar model obtained using MESA (Paxton et al. 2011).
We find that the numerical results depend to some extent on the boundary condition and solar model used. Qualitative results, however, are unchanged, like the overall magnitude of the effect, selection rules, the interaction taking place to a large part near the surface, that the effects of density, pressure and sound speed nearly cancel, and that the changes to solar structure have a larger effect than the flow.
For all boundary conditions and models, the general pattern and the order of magnitude is the same as in Fig. 3. For some boundary conditions, models and mode combinations, however, parts of the pattern visible in the right panel in Fig. 3 are shifted by about 0.5−1.0 × 10^{−12} Hz.
There are a number of good reasons why the observed sensitivity to solar model and surface boundary conditions is not surprising. Firstly, we found that the interaction of p and g modes takes place to a large part near the surface. This layer has thus to be accurately modelled to guarantee an accurate answer. Secondly, p modes are sensitive to relative perturbations to solar structure. Close to the surface, the relative Eulerian perturbations from a g mode depend on the surface boundary conditions used in the eigenfunction computation. Thirdly, we found that the effects of pressure, density, and sound speed almost cancel, with the final frequency shift being nearly two orders of magnitude smaller than the largest term involved in the sum. Substracting terms of similar magnitude may result in an amplification of systematic errors.
As a consequence, all effects have to be modelled in a very accurate way if quantitative predictions from this method are to be used.
Finally, it is unknown whether convection has effects on g modes that are not included in the eigenfunction computation such as by turbulence. For example, if the amplitude of the gmode eigenfunction was significantly reduced at the surface compared to the core, the contribution of the radiative zone to the frequency shifts would be enhanced. Attenuation of p and g modes has not been included in our model.
6. Conclusion
In this paper, we have developed a model for the firstorder frequency shift of solar p modes caused by solar g modes in the framework of nondegenerate timeindependent perturbation theory. To do so, the timedependent perturbation to solar structure caused by a g mode is assumed to be constant at the observational time scale and the degeneracy in mode frequency is assumed to be lifted by rotation.
Our study was motivated by a recent study by Fossat et al. (2017), who proposed to detect g modes by measuring their impact on the wave travel time across the Sun, which is inversely proportional to the large frequency separation of pmode frequencies.
Firstly, we find that indeed g modes can perturb global pmode eigenfrequencies. However, the effect is extremely small. For a gmode surface amplitude of A_{g} = 1 mm s^{−1}, all frequency shifts computed with our model are smaller than 0.1 nHz, or in relative numbers
where results for p modes with l_{p} = 1, 2 and g modes with l_{g} = 2 are taken into account. Measuring such a small effect would require an extremely long observation time. Our results are in broad agreement with the estimates obtained by Scherrer & Gough (2019).
Secondly, we find that the interaction of g and p modes takes place to a large part in nearsurface layers. This is due to two reasons. On the one hand, p modes are most sensitive to relative changes to solar structure close to the surface. On the other hand, relative changes caused by g modes are relatively large in this region, too. This is the case despite the large amplitude of g modes in the solar core, and caused by the sharp decrease of pressure and density near the surface. Furthermore, the effect from structural changes dominates over the effect from the flow caused by the g mode.
Thirdly, we find that a firstorder frequency shift of a p mode by a g mode is only possible if all of the following selection rules are satisfied. The harmonic degree l_{g} of the g mode has to be even, the azimuthal order m_{g} of the g mode has to be zero, and the harmonic degree l_{g} of the g mode has to be smaller than or equal to twice the harmonic degree of the p mode (l_{g} ≤ 2l_{p}). The condition that the harmonic degree of the g mode has to be even has been previously obtained by Kennedy et al. (1993) and Scherrer & Gough (2019).
We caution that we only modelled the firstorder frequency shift in the p modes but not the change in pmode travel time as measured by Fossat et al. (2017). For example, we did not model the effect of gmodes on the pmode eigenfunctions, which may result in different selection rules.
While challenging, the exploration of various methods for detecting g modes certainly remains a worthwhile endeavour.
Acknowledgments
We thank Jesper Schou, Aaron Birch, and Damien Fournier for scientific discussions, Earl Bellinger and Warrick Ball for help with installing and running GYRE, and Richard Townsend for answering questions on eigenfunction computations via the GYRE forum. We thank the anonymous referee for a number of useful suggestions that improved the manuscript. We used the tomso python package written by Warrick Ball for reading ADIPLS eigenfunctions and solar models. This work was supported in part by the Max Planck Society through a grant on PLATO Science. The computational resources were provided by the German Data Center for SDO through a grant from the German Aerospace Center (DLR). We acknowledge partial support from the European Research Council Synergy Grant WHOLE SUN #810218. HH acknowledges the support of the Erasmus Mundus Joint Masters Programme, AstroMundus. Parts of this work were submitted as a MSc. thesis by HH.
References
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Dordrecht: Springer) [CrossRef] [Google Scholar]
 Appourchaux, T., & Pallé, P. L. 2013, in Fifty Years of Seismology of the Sun and Stars, eds. K. Jain, S. C. Tripathy, F. Hill, J. W. Leibacher, & A. A. Pevtsov, ASP Conf. Ser., 478, 125 [NASA ADS] [Google Scholar]
 Appourchaux, T., Belkacem, K., Broomhall, A.M., et al. 2010, A&ARv, 18, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Barrera, R. G., Estevez, G. A., & Giraldo, J. 1985, Eur. J. Phys., 6, 287 [CrossRef] [Google Scholar]
 Basu, S. 2016, Liv. Rev. Sol. Phys., 13, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Berthomieu, G., & Provost, J. 1990, A&A, 227, 563 [NASA ADS] [Google Scholar]
 Berthomieu, G., & Provost, J. 1991, Sol. Phys., 133, 127 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J. 2003, Lecture Notes on Stellar Oscillations, 5th edn. Available online at http://usersphys.au.dk/jcd/oscilnotes/ on August 2, 2017 [Google Scholar]
 ChristensenDalsgaard, J. 2008, A&AS, 316, 113 [Google Scholar]
 ChristensenDalsgaard, J., & Thompson, M. J. 1997, MNRAS, 284, 527 [NASA ADS] [Google Scholar]
 ChristensenDalsgaard, J., Dappen, W., Ajukov, S. V., et al. 1996, Science, 272, 1286 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Dahlen, F., & Tromp, J. 1998, Theoretical Global Seismology (Princeton: Princeton University Press) [Google Scholar]
 Duvall, Jr., T. L. 2004, in SOHO 14 Helio and Asteroseismology: Towards a Golden Future, ed. D. Danesy, ESA Spec. Publ., 559, 412 [NASA ADS] [Google Scholar]
 Dziembowski, W. A. 1971, Acta Astron., 21, 289 [NASA ADS] [Google Scholar]
 Dziembowski, W. A. 1982, Acta Astron., 32, 147 [NASA ADS] [Google Scholar]
 Fossat, E., & Schmider, F. X. 2018, A&A, 612, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fossat, E., Boumier, P., Corbard, T., et al. 2017, A&A, 604, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fournier, D., Hanson, C. S., Gizon, L., & Barucq, H. 2018, A&A, 616, A156 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gabriel, A. H., Grec, G., Charra, J., et al. 1995, Sol. Phys., 162, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O. 1993, in Astrophysical Fluid Dynamics  Les Houches 1987, eds. J. P. Zahn, & J. ZinnJustin, 399 [Google Scholar]
 Hamilton, M. F., & Blackstock, D. T. 1998, Nonlinear Acoustics (San Diego: Academic press) [Google Scholar]
 Hanasoge, S., & Mandal, K. 2019, ApJ, 871, L32 [NASA ADS] [CrossRef] [Google Scholar]
 Hanasoge, S. M., Woodard, M., Antia, H. M., Gizon, L., & Sreenivasan, K. R. 2017, MNRAS, 470, 1404 [NASA ADS] [CrossRef] [Google Scholar]
 Howe, R. 2009, Liv. Rev. Sol. Phys., 6, 1 [Google Scholar]
 Kennedy, J. R., Jefferies, S. M., & Hill, F. 1993, in GONG 1992. Seismic Investigation of the Sun and Stars, ed. T. M. Brown, ASP Conf. Ser., 42, 273 [NASA ADS] [Google Scholar]
 Kiefer, R., Schad, A., & Roth, M. 2017, ApJ, 846, 162 [NASA ADS] [CrossRef] [Google Scholar]
 Kippenhahn, R., & Weigert, A. 1990, Stellar Structure and Evolution (Berlin, Heidelberg: SpringerVerlag) [CrossRef] [Google Scholar]
 Kumar, P., & Goldreich, P. 1989, ApJ, 342, 558 [NASA ADS] [CrossRef] [Google Scholar]
 Lavely, E. M., & Ritzwoller, M. H. 1992, Phil. Trans. R. Soc. London Ser. A, 339, 431 [NASA ADS] [CrossRef] [Google Scholar]
 Liang, Z.C., Gizon, L., Birch, A. C., & Duvall, T. L. 2019, A&A, 626, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Löptien, B., Gizon, L., Birch, A. C., et al. 2018, Nat. Astron., 2, 568 [NASA ADS] [CrossRef] [Google Scholar]
 Lou, Y.Q. 2001, ApJ, 556, L121 [NASA ADS] [CrossRef] [Google Scholar]
 LyndenBell, D., & Ostriker, J. P. 1967, MNRAS, 136, 293 [NASA ADS] [CrossRef] [Google Scholar]
 Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Perdang, J., & Blacher, S. 1982, A&A, 112, 35 [NASA ADS] [Google Scholar]
 Provost, J., Berthomieu, G., & Morel, P. 2000, A&A, 353, 775 [NASA ADS] [Google Scholar]
 Ritzwoller, M. H., & Lavely, E. M. 1991, ApJ, 369, 557 [NASA ADS] [CrossRef] [Google Scholar]
 Schad, A. 2013, PhD Thesis, Universität Freiburg im Breisgau, Germany [Google Scholar]
 Schad, A., Timmer, J., & Roth, M. 2011, ApJ, 734, 97 [NASA ADS] [CrossRef] [Google Scholar]
 Schad, A., Timmer, J., & Roth, M. 2013, ApJ, 778, L38 [NASA ADS] [CrossRef] [Google Scholar]
 Scherrer, P. H., & Gough, D. O. 2019, ApJ, 877, 42 [NASA ADS] [CrossRef] [Google Scholar]
 Schou, J., Antia, H. M., Basu, S., et al. 1998, ApJ, 505, 390 [NASA ADS] [CrossRef] [Google Scholar]
 Schunker, H., Schou, J., Gaulme, P., & Gizon, L. 2018, Sol. Phys., 293, 95 [NASA ADS] [CrossRef] [Google Scholar]
 Thuras, A. L., Jenkins, R. T., & O’Neil, H. T. 1935, J. Acoust. Soc. Am., 6, 173 [NASA ADS] [CrossRef] [Google Scholar]
 Townsend, R. H. D., & Teitler, S. A. 2013, MNRAS, 435, 3406 [NASA ADS] [CrossRef] [Google Scholar]
 Townsend, R. H. D., Goldstein, J., & Zweibel, E. G. 2018, MNRAS, 475, 879 [NASA ADS] [CrossRef] [Google Scholar]
 Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial Oscillations of Stars, 2nd edn. (Tokyo: University of Tokyo Press) [Google Scholar]
 Weinberg, N. N., Arras, P., & Burkart, J. 2013, ApJ, 769, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Wentzel, D. G. 1987, ApJ, 319, 966 [NASA ADS] [CrossRef] [Google Scholar]
 Westervelt, P. J. 1957, J. Acoust. Soc. Am., 29, 199 [NASA ADS] [CrossRef] [Google Scholar]
 Westervelt, P. J. 1963, J. Acoust. Soc. Am., 35, 535 [NASA ADS] [CrossRef] [Google Scholar]
 Woodard, M. F. 1989, ApJ, 347, 1176 [NASA ADS] [CrossRef] [Google Scholar]
 Woodard, M. F. 2007, ApJ, 668, 1189 [NASA ADS] [CrossRef] [Google Scholar]
 Woodard, M. F. 2016, MNRAS, 460, 3292 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Wave operator
The Eulerian perturbations to the solar structure variables are (e.g. ChristensenDalsgaard 2003; Basu 2016):
As a consequence, the wave operator ℒ in Eq. (2) can be written as
Appendix B: Checking the eigenfunction computation
To check the correctness of the eigenfunction computations, we have performed a number of tests.
Firstly, we reproduced plots of eigenfunctions in the literature, namely Berthomieu & Provost (1990, Figs. 6, 7a, 7b), Provost et al. (2000, Fig. 1), and ChristensenDalsgaard (2003, Figs. 5.10 and 7.4) with the results from our eigenfunction computation, see Figs. B.1–B.3. The results look the same as in the literature.
Fig. B.1. Kinetic energy density (ρr^{2}(ξ_{r}(r)^{2} + l(l + 1)ξ_{h}(r)^{2}) [g cm]) as a function of fractional solar radius r/R_{⊙} for modes with degree l_{p} = l_{g} = 1. The radial orders and frequencies of the modes are given in each panel. Compare to Fig. 1 from Provost et al. (2000). We note that the mode frequencies are slightly different compared to Provost et al. (2000) due to slight differences in the solar model. 

Open with DEXTER 
Fig. B.2. Scaled radial gmode eigenfunction. Compare to Fig. 7.4 in ChristensenDalsgaard (2003). 

Open with DEXTER 
Fig. B.3. Scaled radial gmode eigenfunction. Compare to panel a in Fig. 5.10 from ChristensenDalsgaard (2003). 

Open with DEXTER 
Secondly, we compared our results for gmode eigenfunctions to asymptotic relations, both for the trapping region (Eqs. (7.129) and (7.131) in ChristensenDalsgaard 2003) and the evanescent region, which mostly coincides with the convection zone (Eqs. (7.137) and (7.142) in ChristensenDalsgaard 2003). We find that the asymptotic relations qualitatively well reproduce the displacement eigenfunctions, and in most cases, a reasonable approximation to the displacement eigenfunction can be obtained by glueing together eigenfunctions obtained from the asymptotic relations for the trapping and evanescent regions near the bottom of the convection zone, see Fig. B.4. From the asymptotic displacement eigenfunctions, asymptotic expressions for the structure eigenfunctions can be obtained using Eqs. (A.1), (A.2), and (D.7), which agree well with the full structure eigenfunctions near the surface, see Fig. B.5.
Fig. B.4. Comparison of radial (top) and horizontal (bottom) gmode eigenfunctions (computed with GYRE) with asymptotic formulae. For the convection zone (CZ) and radiative zone (RZ), different approximations were used, see Appendix B. The same mode as in Fig. 2 is shown. 

Open with DEXTER 
Thirdly, we checked for modes of low order that eigenfunctions obtained using GYRE and ADIPLS agree closely provided the same solar model and the same boundary conditions are used.
Fig. B.5. Comparison of gmode solar structure eigenfunctions (computed with GYRE) with asymptotic formulae for the convection zone (CZ). The same mode as in Fig. 2 is shown. See Appendix B for details. 

Open with DEXTER 
Appendix C: The perturbed wave operator without hydrostatic equilibrium
In the following, we compute the perturbation to the wave operator δℒ in the presence of perturbations to solar structure, δρ_{0}, δp_{0}, δ(c_{0}^{2}),δg_{0}, which potentially do not fulfil the hydrostatic equilibrium condition δ∇p_{0} = δ(ρ_{0}g_{0}). Omitting the zero subscript in perturbed quantities, for example δρ = δρ_{0}, and perturbing Eqs. (A.4)–(A.8), we obtain
In addition, in the presence of a perturbation to the flow field δu(r), this flow gives rise to a fifth perturbation term (e.g. ChristensenDalsgaard 2003)
such that .
Appendix D: Perturbation to pmode frequencies
We will now assume that the perturbation to solar structure is given by a g mode as outlined in Eqs. (9)–(14). Given these perturbations, we will determine the perturbation to pmode frequencies using Eq. (20) and Eqs. (C.1)–(C.5).
For brevity in notation, we use
D.1. Some useful relations
To ease the necessary derivations, we assume that in addition to the eigenfunctions ξ_{r, g}, ξ_{h, g}, the Eulerian perturbations p_{1, g}, ρ_{1, g}, Φ_{1, g} are given and we separate their radial and horizontal dependence (e.g. Aerts et al. 2010),
Similarly, we obtain using Eqs. (9)–(13) and Eqs. (D.3)–(D.4),
In addition, we need
and an expression for the divergence of the displacement eigenfunction,
where the function Div_{ξp}(r) is defined as the term in the large brackets in Eq. (D.9).
D.2. Solar structure terms
We first separately look at the perturbations to solar structure, which means contributions to δ(ω_{p}^{2}) arising from δℒ_{1} till δℒ_{3}. We omit one term (δℒ_{4}) for simplicity. This term is expected to have a small contribution to the wave equation. It models the change in Eulerian perturbation to the gravitational potential Φ_{1, p} due to the presence of the g mode. To simplify our derivations, we write δℒ_{i − j} for the jth term in δℒ_{i} in Eqs. (C.1)–(C.3) and we use the relations given in Appendix D.1.
We start by the first term in δℒ_{1} (Eq. (C.1)),
The frequency shift arising from this term can then be evaluated using Eq. (20) and Eq. (5) applied to the pmode eigenfunction,
Here, dΩ = sin θ d θ dϕ is the surface element on the unit sphere.
For the following terms, we proceed in a similar way and obtain
Using similar intermediate steps, we obtain the results for the following terms,
The next term requires a bit more work. We start with
To evaluate the last term in ⟨ξ_{p}, δℒ_{2−2}[ξ_{p}]⟩ arising from Eq. (D.17), we first note that ∇(Ψ_{l′m′} ⋅ Ψ_{lm}) is purely horizontal. The horizontal gradient applied to a function f(θ, ϕ) that is independent of radius, satisfies r∇f = ∇_{h}f. It can then be shown using Gauss’s divergence theorem and the property that
Using this, we can now evaluate
The remaining terms can again be evaluated in a similar way as in Eqs. (D.10)–(D.11),
D.3. Flow terms
We now compute the terms due to the flow field caused by the gmode that are involved in δℒ_{5}. The perturbation to the solar model due to flows is given by Eq. (14), so that, from Eq. (C.5),
The computation of ⟨ξ_{p}, δℒ_{5}[ξ_{p}]⟩ has essentially been done by Lavely & Ritzwoller (1992). This is because a g mode produces a poloidal flow field which can be expanded in vector spherical harmonics as (Lavely & Ritzwoller 1992, Eqs. (1)–(3), (C14))
From Eqs. (C31) till (C33) in Lavely & Ritzwoller (1992), we then obtain,
where the Wigner 3j symbols are given for instance in Dahlen & Tromp (1998), and where
Appendix E: Horizontal integrals and selection rules
We will now evaluate the horizontal integrals appearing in Appendix D and identify selection rules that must be satisfied for a firstorder frequency shift to be present.
We first start by observing that the integration over ϕ, which is part of all horizontal integrals, can always be written as
This integral vanishes when m′=m_{g} ≠ 0. We therefore conclude that a firstorder frequency shift is only present if the g mode has azimuthal degree m_{g} = 0.
We now start with the simplest among the horizontal integrals, which can be evaluated using Gaunt’s formula, see Eq. (F.1),
From the selection rules of the Wigner 3j symbols, we can deduct selection rules for this angular integral, which are summarised in Sect. 4.2. As we will see, the remaining angular integrals either satisfy the same selection rules or evaluate to zero.
Next we consider the angular integral . Using generalised spherical harmonics , for which a definition is given in Appendix F, the vector spherical harmonic Ψ_{lm} can be written as, see Eq. (C.139) in Dahlen & Tromp (1998),
where the vectors are defined as, see Eq. (C.49) in Dahlen & Tromp (1998),
These functions satisfy the relations and . Therefore, we can evaluate the horizontal integral as follows,
Now, Eq. (E.5) can be written in terms of Wigner 3j symbols using Gaunt’s formula (see Appendix F),
Since reversing the sign of the lower row in the Wigner 3j symbol gives a phase factor,
the integral (E.6) vanishes if l′=l_{g} is odd. Consequently, Eq. (E.6) results in the same selection rules as Eq. (E.2).
Next, we consider the second angular integral appearing in Eq. (D.16). We now use the relations
and (see Eq. (C.50) in Dahlen & Tromp 1998), which satisfy , and . We then obtain,
which again follows the same selection rules as the above angular integrals.
Similarly, we obtain
which again has the same selection rules.
We now turn to the selection rules from the flow terms derived in Appendix D.3. We begin by noting two common factors in all terms involved in the integral from the flow term, which therefore can be written (see Eqs. (D.28)–(D.33))
As a consequence, the contribution from the gmode induced flow has the same selection rules as the perturbation to solar structure.
Appendix F: Gaunt’s formula
The surface integral of the product of three generalised spherical harmonics can be written in terms of Wigner 3j symbols using the following formula (see, e.g. Dahlen & Tromp 1998, Eqs. (C.188) and (C.198)), which holds if N_{1} = N_{2} + N_{3},
The definition of generalised spherical harmonics follows slightly different conventions in the literature. We here follow the definition in Dahlen & Tromp (1998, Eqs. (C.90), (C.117), (C.112)),
We note that Y_{lm}^{0} = Y_{lm}.
The Wigner 3j symbols follow certain selection rules. The second Wigner 3j symbol in Eq. (F.1) is zero unless both of the following conditions are satisfied,
Additionally, this Wigner 3j symbols evaluates to zero, if
Appendix G: Solar models and surface boundary conditions
The numerical results presented in Sect. 5 depend to some extent on the solar model and surface boundary condition used. This dependence is summarised in Sect. 5.4. In the following, we give some details on solar model and surface boundary conditions that we used.
G.1. Solar models used
We use the following solar models, (i) model S (ChristensenDalsgaard et al. 1996), (ii) model S, where the density was slightly smoothed in radial direction (see details outlined in Appendix G.2), and (iii) a stellar model obtained by evolving a 1 M_{⊙} star from premain sequence to an age of 4.6 Gy using MESA (Paxton et al. 2011). The stellar model from MESA was obtained using the sample input files in the folder 1M_pre_ms_to_wd in the MESA test suite. We note that this stellar model was not exactly calibrated to solar values. Our main goal in using this model is to get an impression on how sensitive the numerical computation is to nonnegligible differences in solar structure.
G.2. Details on the smoothed version of model S
In order to evaluate the formulas for the frequency shift numerically, it is necessary to compute second derivatives of solar model quantities such as pressure (e.g. in Eq. (D.16)) and density (e.g. for the first derivative of the Eulerian density perturbation ρ_{1, g} in Eq. (D.22)). As the second derivative of the density in model S is not very smooth at some locations in the model interior, we compute a version of Model S that is smoothed locally, thereby keeping the changes in the overall density structure as small as possible.
The density is smoothed radially with a Gaussian window. The width of the Gaussian varies with the radial location in the solar model with a maximum value of σ_{max} = 4.87 Mm in the core and σ_{min} = 20 km at the outer boundary. Between these values, we vary the width smoothly using a cosine bell,
Furthermore, when applying the Gaussian window near the surface boundary, we extend the density with the values it would have in an isothermal atmosphere. This is needed in order to guarantee that the second derivative does not have a kink at the boundary.
G.3. Surface boundary conditions used
We compute eigenfunctions using the surface boundary conditions implemented in the GYRE code (Townsend & Teitler 2013). These are a surface boundary matched to an isothermal atmosphere as implemented in ADIPLS (ChristensenDalsgaard 2008), a surface boundary matched to an isothermal atmosphere in the formulation of Unno et al. (1989), the atmospheric boundary condition described in Dziembowski (1971), and a vacuum outer boundary condition (equivalent to a zero Lagrangian pressure perturbation, ∇ ⋅ ξ = 0).
All Figures
Fig. 1. Eigenfunctions and Eulerian perturbations to solar structure for a typical g mode with l_{g} = 2, n_{g} = −10, ν_{g} = 103.2 μHz (compare to Fig. 5.10 in ChristensenDalsgaard 2003). The mode was normalised to a surface amplitude of 1 mm s^{−1}, which is at the order of observational upper limits (e.g. Appourchaux et al. 2010). Displayed are (from top to bottom) radial and horizontal eigenfunctions, and the Eulerian perturbations to density, pressure, and squared sound speed. Left column: kinetic energy density and absolute changes to solar structure caused by the mode are centred near the solar core. Right column: relative changes to solar structure are very small, and in the case of density, pressure, and sound speed, reach the largest values near the surface. The exact behaviour of the relative Eulerian changes to solar structure at the surface depends on the choice of boundary condition used in the eigenfunction computation. 

Open with DEXTER  
In the text 
Fig. 2. Same as Fig. 1, but for a g mode with n_{g} = −36, ν_{g} = 31.2 μHz. See also Figs. B.4 and B.5. 

Open with DEXTER  
In the text 
Fig. 3. Firstorder perturbations to the pmode frequency due to a gmode oscillation assumed to be frozen over the observation period. Results are shown for g modes and p modes with harmonic degrees l_{g} = l_{p} = 2 and azimuthal orders m_{g} = m_{p} = 0. Left: frequency shift as a function of radial order for mode frequencies ν_{p} < 5 mHz and ν_{g} > 5.67 μHz. The range of modes studied by Fossat et al. (2017), Fossat & Schmider (2018), and Schunker et al. (2018) is inside the dashed rectangle, extending outwards of the plotted region at n_{g} > 200. Right: cuts through the results shown in the dashed region in the left panel for a selection of gmode radial orders. 

Open with DEXTER  
In the text 
Fig. 4. Contribution to the frequency shift of p modes by changes to solar structure variables ρ, p, c^{2}, g, and the flow field u caused by a g mode in absolute values. The total frequency shift δν_{p} and the contribution from an example term (δℒ_{3−2}, see Eq. (D.21)) are also given for comparison. The azimuthal order of the modes is m_{p} = m_{g} = 0. 

Open with DEXTER  
In the text 
Fig. 5. Radial integrands for computing the frequency shift of a p mode with l_{p} = 2, m_{p} = 0, n_{p} = 20 caused by a g mode with l_{g} = 2, m_{g} = 0, n_{g} = −36. When integrating over radius, the total integrand yields the total frequency shift δν_{p} (both panels). The contributions to the total integrand resulting from the structural changes and the flow are also shown (right panel). The spike at the surface boundary is due to numerical difficulties in computing second order derivatives at the boundary. Therefore, the three outermost radial grid points were discarded. 

Open with DEXTER  
In the text 
Fig. 6. Contribution to the frequency shift from the convection zone and from the radiative zone as a function of pmode frequency for several typical g modes considered in this study. Results are shown for l_{p} = l_{g} = 2, m_{p} = m_{g} = 0. 

Open with DEXTER  
In the text 
Fig. B.1. Kinetic energy density (ρr^{2}(ξ_{r}(r)^{2} + l(l + 1)ξ_{h}(r)^{2}) [g cm]) as a function of fractional solar radius r/R_{⊙} for modes with degree l_{p} = l_{g} = 1. The radial orders and frequencies of the modes are given in each panel. Compare to Fig. 1 from Provost et al. (2000). We note that the mode frequencies are slightly different compared to Provost et al. (2000) due to slight differences in the solar model. 

Open with DEXTER  
In the text 
Fig. B.2. Scaled radial gmode eigenfunction. Compare to Fig. 7.4 in ChristensenDalsgaard (2003). 

Open with DEXTER  
In the text 
Fig. B.3. Scaled radial gmode eigenfunction. Compare to panel a in Fig. 5.10 from ChristensenDalsgaard (2003). 

Open with DEXTER  
In the text 
Fig. B.4. Comparison of radial (top) and horizontal (bottom) gmode eigenfunctions (computed with GYRE) with asymptotic formulae. For the convection zone (CZ) and radiative zone (RZ), different approximations were used, see Appendix B. The same mode as in Fig. 2 is shown. 

Open with DEXTER  
In the text 
Fig. B.5. Comparison of gmode solar structure eigenfunctions (computed with GYRE) with asymptotic formulae for the convection zone (CZ). The same mode as in Fig. 2 is shown. See Appendix B for details. 

Open with DEXTER  
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.