Giant star seismology

The internal properties of stars in the red-giant phase undergo significant changes on relatively short timescales. Long near-uninterrupted high-precision photometric timeseries observations from dedicated space missions such as CoRoT and Kepler have provided seismic inferences of the global and internal properties of a large number of evolved stars, including red giants. These inferences are confronted with predictions from theoretical models to improve our understanding of stellar structure and evolution. Our knowledge and understanding of red giants have indeed increased tremendously using these seismic inferences, and we anticipate that more information is still hidden in the data. Unraveling this will further improve our understanding of stellar evolution. This will also have significant impact on our knowledge of the Milky Way Galaxy as well as on exo-planet host stars. The latter is important for our understanding of the formation and structure of planetary systems.

opacity it is not trivial to study the internal structure of stars (Eddington, 1926). The only way to probe and study the internal stellar structure directly is through global stellar oscillations. Many stars across the Hertzsprung-Russell (HR) diagram oscillate (driven by different mechanisms), which allows for in-depth studies of stars and stellar evolution.

History
Asteroseismology -the study of global properties of stars and their internal structure through their global intrinsic oscillations -is already more than a century old. An early remarkable result was obtained for Cepheids revealing the period-luminosity relation (Leavitt and Pickering, 1912), which played an important role in measuring distances of galaxies and star clusters and ultimately the expansion of the universe. Observations of pulsators which display large (photometric) variations have been ongoing ever since. These stars include among others, Cepheids, high amplitude δ Scuti stars, RR Lyrae stars, white dwarfs and Mira's. An early result for which asteroseismic observations and stellar modelling were inconsistent was presented by Petersen (1973) on double-mode Cepheids. He showed that stellar masses inferred from the periods and period ratios were inconsistent with the location of the stars in the HR diagram as predicted by stellar evolution models. This problem motivated a revision of stellar opacity tables which led to very good agreement between the models and the observations (e.g., Moskalik et al, 1992). Despite the long-term research on such bright large amplitude oscillators there are still open questions concerning their oscillations (e.g., Blazkho effect, Blažko, 1907) and internal properties. For recent reviews on classical oscillators see e.g. Balona (2010); Handler (2013); Szabó et al (2015) and references therein.
Lower-amplitude oscillations were not detected for several more decades. In the 1960's oscillations were first discovered in the Sun (Leighton et al, 1962). Subsequent observations provided details of the solar interior and constraints for general stellar modelling beyond that previously possible Christensen-Dalsgaard, 2002). Seismology was pivotal in the solar neutrino problem (Bahcall, 1972;Trimble and Reines, 1973) and confirmed that the solution lies within particle physics (Elsworth et al, 1990a). Through sensitive neutrino detections, McDonald et al (2001) later showed that low energy neutrinos do indeed change flavour. The Nobel Prize 2015 for physics was awarded for this discovery. Unfortunately, our understanding of the Sun remains incomplete. The 'new solar abundances' (Asplund, 2005;Asplund et al, 2009) result in a solar structure that deviates much further from that obtained from the oscillations (Bahcall et al, 2005;Guzik, 2008;Antia, 2008, 2013) as compared with the 'old solar abundances' (Grevesse and Noels, 1993;Grevesse and Sauval, 1998). Additionally, the 11-yr solar activity cycle discovered by Schwabe (1843) and visible in the solar oscillations (Elsworth et al, 1990b;Howe et al, 2002;Hathaway, 2015) is not fully understood. This includes dynamo processes involving interaction between rotation and convection, restructuring the magnetic field between the toroidal and polar components, in a manner that is yet to be resolved, and giving rise to the surface manifestation as sunspots when flux tubes break through the solar surface (e.g., Charbonneau, 2014;Cameron and Schüssler, 2015). Additionally, the reasons for the long solar minimum between cycle 23 and 24 are still subject of discussion Basu, 2013;Jiang et al, 2015). Furthermore, oscillations sensitive to the activity suggest the presence of a 2-yr cycle in addition to the 11-yr cycle  which is yet to be understood.
The road to the detection of low-amplitude oscillations in stars other than the Sun was paved by radial-velocity measurements -much like the variations used to discover pulsations in the Sun. Ground-based spectroscopic surveys first identified excess oscillation power in Arcturus (Smith et al, 1987), Procyon , η Bootis (Kjeldsen et al, , 2003, ξ Hydrae (Frandsen et al, 2002), and µ Her (Bonanno et al, 2008). These seminal discoveries relied on single-site spectroscopic programs and inspired longer multi-site spectroscopic campaigns on nearby bright stars, such as α Cen A Butler et al, 2004) and B (Carrier and Bourban, 2003;Kjeldsen et al, 2005), β Hydrae , ν Indi , η Serpentis  and Procyon (Hekker et al, 2008;Arentoft et al, 2008;Bedding et al, 2010). Complementary photometric multi-site campaigns were also pursued (Stello et al, 2006a(Stello et al, , 2007 and although both strategies revealed detailed oscillations patterns, they were yet unable to fully constrain the internal structure of the stars (Miglio and Montalbán, 2005;Huber et al, 2011).
Space-based instruments contributed to the progression of the field with the Hubble Space Telescope (Gilliland, 2008) and the star-tracker of the WIRE (Wide field InfraRed Explorer) mission (Buzasi, 2002;Stello et al, 2008) both detecting lowamplitude oscillations in other stars. Soon after, the era of dedicated photometric space-based missions heralded a revolution for asteroseismology. Observations with the MOST mission (Microvariability and Oscillations of STars; Matthews et al, 2000) have contributed significantly to our understanding of classical pulsators such as δ Scuti stars (Casey et al, 2013), Slowly Pulsating B stars (SPB; Aerts et al, 2006;Jerzykiewicz et al, 2013) and rapidly oscillating Ap stars Gruberbauer et al, 2011), as well as other pulsators such as pre-main sequence stars  and Wolf-Rayet stars (David-Uraz et al, 2012). CoRoT (Convection, Rotation and planetary Transits; Baglin et al, 2006) and Kepler (Borucki et al, 2008) have contributed greatly to the detection of pulsations in, and understanding of, many different kinds of stars such as; massive stars Degroote et al, 2010;Kurtz et al, 2015), RR Lyrae stars (Kolenberg et al, 2010) and sub-dwarf B stars (Østensen et al, 2014). Additionally, many break-through results have been reported on low-amplitude oscillations in Sun-like stars, subgiants and red-giant stars (see recent reviews by Chaplin and Miglio (2013) and Hekker (2013) and references therein).

Giant star evolution
The stars that are of interest for this review are subgiant and red-giant stars, i.e. low-to intermediate-mass stars (0.8 -5 M ⊙ ) that have left the main sequence. These stars are of great interest as they undergo large internal structure changes on a relatively short timescale. Stars enter the subgiant phase when they have ceased central hydrogen burning and only burn hydrogen in a shell around the core. Subgiants decrease their surface temperature before becoming more luminous while ascending the red-giant branch (RGB). When a star evolves from the end of the main-sequence to the tip of the RGB its radius increases by one or two orders of magnitude. The luminosity peaks at the onset of helium-core burning, which occurs either in a rapid process (heliumflash) for low-mass stars with a degenerate core or more gently for higher-mass stars with a non-degenerate core. Following the onset of helium-core burning a star enters the horizontal branch during which the star returns to a quiescent phase with central helium and hydrogen shell burning. In this phase the luminosity and radius are significantly reduced compared with those at the tip of the RGB. This quiescent phase is also known for more metal-rich stars as red-clump (RC) or secondary-clump phase depending on the pre-helium ignition conditions in the core. When central helium burning ceases, hydrogen and helium shell burning are the remaining nuclear sources of energy. During this so-called asymptotic-giant-branch (AGB) phase the star again increases its luminosity and radius while the surface temperature decreases.

Observations of stellar oscillations
Oscillations can be determined from timeseries data of either intensity variations or radial-velocity (RV) variations. Intensity variations reflect the brightness variations of a star induced by stellar oscillations. RV variations reveal the outward and inward movement of the stellar surface due to stellar oscillations through a Doppler shift of the spectrum. Part of a timeseries of intensity variations of a red-giant star observed with the Kepler space telescope is shown in the top panel of Fig. 1.
There are significant differences between intensity and RV variation measurements in the sense of their sensitivity to other intrinsic stellar features. For example, granulation (the visible effect of convection at a star's surface) has higher amplitude in intensity variations than in RV variations, relative to the oscillations. Additionally, stars other than the Sun can be observed only in integrated light, causing cancellation effects that differ between intensity and RV observations. Since solar-like oscillations are mainly in the radial direction RV observations with the projection onto the line of sight have a reduced sensitivity to the oscillations near the limb, increasing the response to modes of slightly higher degree and hence the diagnostic potential, compared with intensity observations.
To extract oscillation frequencies the timeseries data are most commonly transformed to frequency space by a Fourier transform. The resulting Fourier power spectrum reveals the oscillation frequencies as sharp peaks. An example of a Fourier power spectrum of a red-giant star observed with the Kepler space telescope is shown in the bottom panel of Fig. 1.
In this review we will discuss the internal structure and structure changes of subgiants and red giants including AGB stars in more detail (Sect. 2). These stars oscillate with intrinsically damped oscillations stochastically excited by convection in the outer layers of the stars, i.e., some of the convective energy is converted into Fig. 1 A 20 day long subset of a Kepler photometric timeseries (or light curve) of a red-giant star (KIC 9145955, top panel) and the Fourier power spectrum of the same star using a 1060-days long dataset (bottom panel). Note that only the frequency range in which oscillations occur is shown here. energy of eigenmodes of the star. Consequently these oscillations allow for observational investigations of internal stellar structure. As such oscillations are present in the Sun, they are referred to as solar-like oscillations. Solar-like oscillations are expected to be present in all stars with turbulent outer layers. In Section 3 we discuss the diagnostics that can be obtained from timeseries data and Fourier power spectra. An overview of stellar pulsation theory is presented in Section 4. We highlight ground-breaking results from the past years in Section 5. Finally, we discuss some promising prospects of asteroseismology and stellar structure of giants in Section 6.

Giant star evolution
Here we provide an overview of the internal structures of low-to intermediate-mass stars (roughly 0.8 to 10 M ⊙ ) in their respective evolutionary stages. We discuss the low-mass stars and intermediate-mass stars separately. Low-mass stars are those stars that ignite helium in the core under degenerate conditions. This occurs in stars with masses between ∼0.48 M ⊙ and ∼2 M ⊙ . The lower limit is defined by the lower limit of the critical mass needed to ignite helium, while the upper limit depends on the chemical composition of the stars. If mass loss is small enough, and time long enough, all low-mass stars will go through the giant phases discussed below. At the current Fig. 2 Helium core mass at ignition vs. stellar mass for stellar models of solar metallicity computed with the MESA stellar evolution code (Paxton et al, 2011). age of the universe (∼14 Gyr) only stars with M > 0.8 M ⊙ will have reached these late stages of evolution, unless substantial mass transfer in a binary system has taken place. Intermediate-mass stars are stars that do not develop a degenerate core and have a more gentle onset of core-helium burning. These stars range in mass between ∼2 M ⊙ and 8-10 M ⊙ depending on metallicity (Kippenhahn et al, 2012). For stars with degenerate cores the helium-core mass at ignition is the same, while this decreases for more massive stars without degenerate cores (see Fig. 2). To illustrate the internal structure changes of low-and intermediate-mass stars we show the paths of a 1 M ⊙ and a 3 M ⊙ model in the Herzsprung-Russel diagram (Fig. 3) and so-called Kippenhahn diagrams of these models in Figs 4 and 5.
The description provided here is aimed to provide an insight in the many internal structure changes that a giant star undergoes. For a more complete picture and details of stellar evolution we refer the reader to Kippenhahn et al (2012).

End of Hydrogen core burning phase: end of main sequence
For stars with masses below 1-1.5 M ⊙ (depending on chemical composition) hydrogen (H) in the central regions fuse to helium via proton-proton (pp) chains under conditions of radiative energy transport. Conversely, stars with masses above this range develop a convective core on the main-sequence, with the CNO (Carbon-Nitrogen-Oxygen) catalytic reactions predominately responsible for the conversion of H into 4 He. CNO burning takes place at higher temperatures (T ) than proton-proton interactions. Additionally, the energy production (ǫ) of CNO burning is more temperaturedependent than that of proton-proton chains (ǫ CNO ∼ T 17 vs. ǫ pp ∼ T 4 ; Iliadis, 2007). The respective conditions in the core significantly impact the structure and evolution at the end of the main sequence and therefore the two regimes, i.e. stars with radiative cores and convective cores on the main sequence, are discussed separately.
Stars with radiative cores For stars with a radiative core the fusion of H to 4 He is predominantly enacted via the pp chains. The increase in the mean molecular weight (µ) resulting from hydrogen fusion affects hydrostatic equilibrium. To sustain pressure support to the overlying layers of the star the core must contract. This increases the temperature in the inner parts of the star, leading to an increase in the efficacy of energy transport and hence in the luminosity (∼0.7% increase in brightness every 100 Myr in a 1 M ⊙ star). This is matched by a corresponding increase in the energy generation.
The gradual depletion of H from the radiative core generates a smooth transition to a chemically inhomogeneous structure with hydrogen-burning in an extended shell around a growing inert degenerate helium core. Complete depletion of hydrogen in the centre marks the end of the main-sequence and the transition to the hydrogen shell burning phase.
Stars with convective cores There is a physical limit to which energy can be transported by radiation. The Schwarzschild criterion (Schwarzschild, 1906) states that if the temperature gradient inside a star is too steep, convection will take over as the primary means of energy transport. In addition to efficient energy transport, the associated bulk mass motions of convection ensure that the composition of any convective region is well mixed. In a more formal sense convection is activated once the radiative temperature gradient exceeds the adiabatic temperature gradient (Schwarzschild, 1906).
The Schwarzschild criterion defines the likely regimes in stellar interiors in which convection will develop. The first regime is where there is a large energy flux, the second where the stellar material is opaque to photons (such as when there is a large presence of molecules) and energy transfer by radiation is therefore inefficient. In stars with M 1.1 M ⊙ the central conditions are sufficient for the activation of the CNO cycle -reactions that are highly temperature dependent. A consequence for these stars is that the burning region becomes ever more centrally concentrated, and the large energy flux and steep temperature gradient drives a convective core.
An alternative to the Schwarzschild criterion is the condition derived by Ledoux (1947) which, in addition to the temperature gradients, takes into account the spatial variation of the mean molecular weight. In some stars, particularly those with M 2.25 M ⊙ , a composition gradient may develop outside the shrinking (Ledoux) convective core. Application of the Schwarzschild criterion would render this region convectively unstable, however according to Ledoux, the presence of a molecular weight gradient can have a stabilizing effect. The true behaviour of material in the stellar interior under these circumstances remains unclear (Gabriel et al, 2014). However, in stellar evolution codes some form of slow (slow compared to convection) 2) and the helium core burning phase for the 3 M ⊙ model (see Section 2.2.4). The models are computed using the MESA stellar evolution code (Paxton et al, 2011) with solar metallicity. The letters indicate different phases of evolution: A = zero-age main-sequence; B ′ = core hydrogen mass fraction ≈ 0.05, B = start of thick shell burning; C = maximum extent of thick shell (in mass); D = start of thin shell burning; E = maximum bump luminosity; F = minimum bump luminosity; G = tip of the red-giant branch; H = onset of helium-shell burning; I = end of helium-core burning, and correspond to the letters indicated in Figs 4 and 5. mixing is applied in these regions in order to match various observational constraints (see e.g. Lattanzio, 1983;Langer et al, 1985). This slow mixing process can also be referred to as semi-convection.
In stars below ∼ 2 M ⊙ with a convective core the core grows for most of the main-sequence lifetime. This is caused by the increasing contribution to the energy production from the highly temperature-sensitive CNO cycle, as oxygen is gradually converted to nitrogen in parts of the cycle. A further extension of the mixed core can be caused by material that approaches the boundary of stability with momentum and overshoots into the radiative layer. This process, convective-core overshoot, extends the burning region and brings in fresh H-rich fuel prolonging the main-sequence evolution. As with a radiative core the decrease in hydrogen abundance with nuclear burning and the resulting increase in µ causes an increase in luminosity. However, with a convective core this depletion takes place uniformly in an extended region. When hydrogen is nearly depleted in this region (X core ≈ 0.05, B ′ in Fig. 3) the star contracts to maintain the energy production. This contraction leads to an increase in effective temperature (T eff ) and luminosity (L) until hydrogen is completely depleted in the centre (the "hook" in the Hertzsprung-Russell diagram: B ′ -B in the 3M ⊙ track in Fig. 3). At this point the central burning and convection cease abruptly. This is the end of the main-sequence phase.

Hydrogen shell burning phase: subgiants and red-giant branch stars
Subgiants Without the central nuclear reactions the star must find an alternative way to generate energy. The star turns to another available source of energy in the release of gravitational potential. For stars with M 1.1 M ⊙ the central density is large enough that electron degeneracy dominates and provides significant pressure support. Therefore, low-mass stars can remain in thermal and hydrostatic equilibrium with a degenerate, isothermal core as they smoothly transition to hydrogen-shell burning. As a consequence the contraction phase and transition to giant is gradual and the timescale much longer compared with the higher-mass counterparts (see Section 2.2.2). Shell burning is initially in the form of an extended burning region outside the core. Slowly, the core mass increases as the ashes of hydrogen-shell burning are deposited on the He core. As the core is degenerate this corresponds to a reduction in radius, accompanied by an expansion of the envelope. The boundary where contraction changes to expansion is located near the hydrogen-burning shell. This behaviour is typical of a more general exchange in evolving stellar models between contraction and expansion at shell-burning sources, referred to as the mirror principle.
The shell-burning region, dominated by CNO cycling, is confined to increasingly narrow mass. The envelope expands and cools while the star evolves from the main sequence towards the Hayashi line. The Hayashi line is the locus in the Herzsprung Russell diagram of fully convective stars, where a star cannot decrease its temperature further (otherwise it cannot maintain hydrostatic equilibrium). Thus, upon approaching the Hayashi line further increase of the radius causes an increase in luminosity. A large convective region develops in the envelope due to the increased photospheric opacity at lower temperatures (contribution from H − ions). The star is now on the redgiant branch, just on the hot side of the Hayashi line with a large convective envelope on top of a small radiative core.

Red-giant branch stars
On the red-giant branch, H-burning ashes are deposited on the degenerate core as the shell burns and moves outwards. Due to the degenerate conditions an increase in mass results in contraction and heating of the inert helium core. This also heats the hydrogen-burning shell, which reacts to the higher temperatures by compressing the burning region and increasing the energy generation. The density contrast becomes so large that the envelope and core are essentially decoupled. Therefore, the hydrogen-burning in the shell is completely determined by the properties of the helium core and not by the envelope. Hence, the luminosity of the star is now related to the mass of the inert degenerate helium core and does no longer depend on the total mass of the star. Therefore, stars with the same core mass, but a spread in total mass, follow the same path in the Hertzsprung-Russell diagram.
First dredge-up The convective envelope penetrates deep into the star to the regions where the chemical composition has been altered by nuclear processes that took place during the main sequence. The processed material is then subsequently transported to the surface. This is the first "dredge-up" phase, i.e. chemical elements from deeper layers are dredged up towards the surface of the stars: for example the 12 C/ 13 C ratios are lowered. The convective region reaches a maximum depth in mass and recedes because of the advance of the hydrogen-shell, leaving behind a chemical (mean molecular weight) discontinuity.
Bump The hydrogen shell, in which burning takes place, moves gradually outwards (in mass) while the helium-core mass steadily increases. In a simplified picture the luminosity decreases when the hydrogen-burning shell reaches the chemical discontinuity left behind by the deepest extent of the convective envelope due to the decrease in the mean molecular weight at the chemical discontinuity, causing the luminosity to decrease again following L ∝ µ 7 M 7 core (Refsdal and Weigert, 1970). After this the core mass keeps increasing at constant µ causing a resumption of the increase in the luminosity.
In fact, the situation is more complex with the luminosity beginning to decrease prior to the shell burning through the discontinuity. As discussed by Christensen-Dalsgaard (2015) the reason is that the decrease in µ above the discontinuity starts affecting the hydrostatic structure, and hence the temperature, within and above the hydrogenburning shell before it reaches the discontinuity. This causes the decrease in the luminosity. This zig-zag in the evolution path is the so-called RGB-bump (see right inset in Fig. 3). The bump is visible for stars up to about 2.2 M ⊙ as an over-density of stars in stellar clusters at the bump luminosity. For more massive stars helium-core burn-ing starts before the hydrogen-burning shell reaches the composition discontinuity left behind by the first dredge-up. In these cases a bump-like structure is not present on the red-giant branch.
High-luminosity red-giant branch stars The process on the RGB continues till the core reaches a temperature of ∼10 8 K (at a core mass of ∼0.45 M ⊙ ) at which helium is ignited in a thermal run-away process. This is the so-called helium flash. We note here that additional mixing processes, such as thermohaline mixing (e.g. Eggleton et al, 2006Eggleton et al, , 2008Charbonnel and Zahn, 2007;Charbonnel and Lagarde, 2010;Angelou et al, 2011Angelou et al, , 2012, are necessary to include in models of stars ascending the red-giant branch to match the observed chemical compositions of these stars.

Onset of Helium burning: He-flash
At a temperature of ∼10 8 K in the inert helium core, helium fusion can be ignited in a triple-alpha process. 1 In the highly degenerate core the pressure does not depend on the temperature and hence there is no thermostatic control to expand and cool the core. The onset of (unstable) burning in these degenerate conditions results in a thermal runaway process creating for a very short time (of order a few hours!) an enormous overproduction of nuclear energy. This energy is absorbed by the expansion of non-degenerate layers outside the degenerate core and does not reach the stellar surface.
The onset of helium fusion takes place at the location of maximum temperature. The temperature is generally at its maximum in a concentric shell around the centre of the degenerate core due to gravo-thermal energy and neutrino losses. Stellar models predict that the first main helium flash is followed by a series of sub-flashes. Each subsequent sub-flash is located closer to the centre such that eventually the degeneracy in the centre is completely lifted and the star is back in equilibrium with helium burning in a convective core. As neither the energy of the flash, nor the energy of the subflashes reaches the stellar surface, the existence of sub-flashes in real stars is not observationally confirmed, i.e., they could be artefacts of stellar models.

Helium core burning phase: red clump
The star has now two sources of energy generation: hydrogen burning in a shell around the core producing helium, while in the core helium is consumed to produce carbon and oxygen. Due to the expansion and accompanying decrease in density and temperature of the hydrogen-burning shell after the helium-flash it generates less energy (this is however still the main source of energy generation in the star). Therefore, the luminosity decreases while the core expands and the envelope contracts (mirror principle, see Section 2.1.2). The star is back in equilibrium and settles in the red clump.
All stars that have gone through a helium flash have very similar core masses and hence very similar luminosities on the horizontal branch. Therefore, all stars that go through a He-flash end up at a very similar spot in the Herzsprung-Russell diagram, with only some dependence on T eff owing to their total masses (with lower masses being slightly hotter) and composition (stars with higher contents of heavier elements are cooler).
Properties of the convective core The helium core burning stars have a central convective core that becomes enriched in carbon and oxygen during helium burning. The opacity in the temperature-density regime present in the core is dominated by freefree transitions and increases with increasing carbon and oxygen abundance. This causes an increase in the radiative temperature gradient in the core leading to a discontinuity in the radiative temperature gradient at the boundary between the convective core and the radiative envelope. In models with no mixing beyond the convective core the nuclear burning gives rise to an increasing composition discontinuity at the edge of the core. Convective overshoot and/or semiconvection increases the size of the core and, depending on the implementation, may lead to a smooth composition profile or further discontinuities in composition (see Constantino et al, 2015, for an overview of these processes). As discussed in Sections 4.2.4 and 4.2.5 the detailed properties of the composition profile can have a strong effect on the behaviour of the oscillations of the star. Towards the end of the central helium burning sudden mixing in the models between the carbon-oxygen-rich core and the overlying helium-rich layers may occur at the edge of the core. This leads to an abrupt increase in the helium content of the core and a loop in the Herzsprung-Russell diagram. These are so-called breathing pulses. The occurrence and appearance of these breathing pulses in stellar models depend on the criterion used for convection and may only be an artefact of the way convection is included in models and may not be present in real stars.
After some time helium is exhausted in the convective core and the star will undergo some rapid evolution towards a shell-burning phase with burning taking place in a helium shell and hydrogen shell surrounding the core. This phase of evolution is the asymptotic giant branch (AGB).

Helium and Hydrogen shell burning phase: Asymptotic giant branch
In the low-mass regime, the asymptotic giant branch is characterized by an inert carbon-oxygen core surrounded by two burning shells of which the helium shell is thermally unstable. In this phase a star is again moving in the Herzsprung-Russell diagram towards the Hayashi line and at the same time increasing its luminosity and radius. In the early AGB (E-AGB) phase hydrogen is burning outwards and the temperature in this shell drops. Consequently, the hydrogen-shell supplies only a small fraction of the energy for some time. However, as the temperature in the hydrogen shell increases again in between thermal pulses (see below) burning is recovered and dominates the energy production.
Thermal pulses As the star ascends the asymptotic giant branch, the helium-burning shell narrows while providing most of the energy to the stellar surface. Eventually the helium-burning shell advances in mass towards the hydrogen-burning shell and their separation, the inter-shell region, becomes too narrow. This, along with the high temperature dependence of the helium-burning reactions, results in a thermal runaway and the onset of the thermally pulsing (TP) AGB phase.
The thermally pulsing AGB phase is characterized by long periods of quiescent hydrogen-shell burning, followed by instabilities of the helium-burning shell (Schwarzschild and Härm, 1965;Weigert, 1966). Each instability or 'helium-shell flash' grows in amplitude for the first 5-10 pulses before approaching a maximum in helium luminosity. As per Iben (1981) we describe the TP-AGB cycle by four distinct phases: -On phase: The sudden deposition of energy from the shell flash drives an intershell convection/burning zone cycling the products of the triple-alpha process into the region below the hydrogen shell. -Power down phase: Although the carbon-oxygen core is highly electron degenerate the shell-flash instability occurs in a non-degenerate region of the stellar interior. The energy generated by the flash helps drive expansion of the star and thereby extinguishing the hydrogen-burning shell. The flash is able to generate luminosities of the order 10 7 or 10 8 solar luminosities. However, due to the expansion this increase in luminosity is not manifested at the surface. -Third dredge-up 2 : As the star expands and cools, convection is able to penetrate beyond the hydrogen-burning shell into regions homogenised by the inter-shell convection zone. Hence the products of He burning are mixed into the envelope where they can be observed at the stellar surface. -The inter-pulse phase: Eventually the helium-burning luminosity drops below the surface luminosity and the outer regions can once more contract. The hydrogen shell can reignite where it provides most of the luminosity until once again the interior conditions arise for a successive helium-shell instability. The inter-pulse phase lasts for ≈10 4 years.
Third dredge-up Third dredge-up plays an important role in the chemical enrichment of the galaxy. The significant amount of carbon produced in these stars and brought to the stellar surface is in some cases able to raise the C/O abundance ratio to become larger than unity. In addition TP-AGB stars have been identified as a site of the s-process nucleosynthesis 3 , which is responsible for the production of half of the elements beyond iron. Through efficient mass-loss processes these elements are expelled into the interstellar medium.
Post-thermally-pulsing AGB phase At some point the envelope mass is insufficient to allow TP to continue. When the envelope mass drops below about a few (∼5) per cent of the total mass, the envelope contracts and shell burning extinguishes: the star becomes a white dwarf. This is the post-AGB phase. The number of thermal pulses the star experiences and the final white dwarf mass depend on the competition between mass loss and core growth in the AGB phase. For extensive reviews of the AGB stars we refer to Iben and Renzini (1983), Herwig (2005) and Karakas and Lattanzio (2014).

End of Hydrogen core burning phase: end of main sequence
Intermediate-mass stars have a convective core on the main sequence and show the same behaviour as low-mass stars with a convective core as described in section 2.1.1.

Hydrogen shell-burning phase: subgiants and red-giant branch stars
After a short resettling at the end of the main sequence, hydrogen-shell burning intensifies in the region around the core. As for the low-mass star this shell burning steepens the hydrogen profile at the edge of the He-core leading to a narrowing of the burning shell when the lower hydrogen tail has been consumed. The core contraction and envelope expansion below and above the shell burning layer respectively (mirror principle) increase the radius of the star.
Hertzsprung gap In intermediate-mass stars the core is non-degenerate following central hydrogen exhaustion. Schönberg and Chandrasekhar (1942) demonstrated that there is a maximum relative core mass that an isothermal, non-degenerate core can have whilst maintaining hydrostatic and thermal equilibrium and deriving all its energy from a nuclear burning shell source. Without electron degeneracy to supply additional pressure support, contraction on a Kelvin-Helmholtz timescale will develop when the core mass exceeds this Schönberg-Chandrasekhar limit, so that the star can maintain equilibrium between the energy it generates in the interior and that which it loses at the surface. Whether stars in the intermediate-mass range reach the Schönberg-Chandrasekhar limit depends on their hydrogen-exhausted core mass. This mass depends on the amount of overshoot the convective core experienced on the main sequence. Upon reaching the Schönberg-Chandrasekhar limit stars cross the Hertzsprung-Russell diagram rapidly and move onto the red-giant branch. The fast timescales involved in this phase of evolution leads to a dearth of intermediate mass observed in this region of the HR diagram: the so-called Hertzsprung gap. At the bottom of the red-giant branch, intermediate-mass stars develop a deep outer convective region to transport energy more efficiently and prevent the star from cooling beyond the Hayashi line. This deepening of the convective envelope causes a change in the surface chemical composition due to the first dredge-up (Section 2.1.2). At the same time the core continues to contract and heat. When the core temperature has increased to about 10 8 K, helium is ignited. We note that the bump (Section 2.1.2) is only present on the red-giant branch for stars with masses roughly below 2.2 M ⊙ . More massive stars ignite helium before the hydrogen-shell reaches the chemical discontinuity left behind by the first dredge-up.

Onset of Helium burning: non-degenerate ignition
In the non-degenerate regime, the luminosity at which helium ignites is a monotonically increasing function of the core mass. Because the pressure and temperature are related, the thermostatic feedback allows intermediate mass stars to gently ignite helium in their core.

Helium core burning phase: secondary clump
Stars in the secondary clump have two sources of burning. Firstly, helium burning that produces carbon and oxygen is present in the core. Like the CNO reactions, energy production via the triple-alpha process is highly temperature dependent. The reactions are concentrated towards the centre and give rise to a convective core. Burning proceeds quiescently with the core growing as a function of time. Secondly, hydrogen burning is taking place in a shell around the core. The latter provides most of the total energy output as a rather small release of nuclear energy is sufficient in the core to compensate for the energy loss from the core and prevent the core from contracting. As the core masses of these stars can be different when they ignite helium, these stars do not "clump" as the red-clump stars. Instead they form the secondary clump at higher luminosities and effective temperatures. Additionally these stars loop through the red-giant region during central helium burning (see left inset in the HR Diagram (Fig. 3) where the stars leave the Hayashi line to become hotter and move back to the Hayashi line). The temperature range that the loops cover increases with increasing stellar mass. The cause of these loops lies in the chemical composition profile in the central regions of stars that had a convective core on the main-sequence.

Helium and Hydrogen shell burning phase: Asymptotic giant branch
When helium burning terminates in the core the burning continues in two shells around the inert carbon-oxygen core. For stars with masses of 4-8 M ⊙ the hydrogen shell is at best barely active. This allows the convective envelope to penetrate down reaching layers through which the hydrogen shell has burned (similar as for low-mass stars on the red-giant branch). This so-called second dredge-up brings processed material that is generated by helium and hydrogen burning such as carbon, oxygen, nitrogen and helium to the surface. For all stars that do not ignite carbon in their core, the AGB proceeds in roughly the same fashion as for low-mass stars (see also Section 2.1.5). They experience episodic thermal pulses and efficient mass loss until they lose most of their envelope and become a white dwarf. However, their s-process nucleosythesis will differ greatly depending on stellar mass.

Why do stars become giants?
The description written above represents our current understanding of stellar evolution based on stellar models and observations. However, the reason why stars become red giants is actually not understood. The mirror phenomenon mentioned before seems to play an essential role in stars to become a red giant. However, it is not understood what physical mechanism(s) drive the mirror nor what other physical mechanisms are essential in a star to become a giant.
A number of studies have investigated the question 'Why do stars become giants?' proposing reasons related to central gravitational field (Hoppner and Weigert, 1973;Weiss, 1983), the effective equation of state (Eggleton and Cannon, 1991;Eggleton et al, 1998), gravothermal instability in the core (Iben, 1993), thermal instabilities in the stellar envelope (Renzini et al, 1992) and mean molecular weight gradient (Stancliffe et al, 2009). Some of these studies have been met with fierce opposition whilst others have devised conditions that were later shown to be necessary but not sufficient in all stars. Currently, it is clear that a strong gravitational field and a mean molecular weight gradient play an important role (Stancliffe et al, 2009) in stars for them to become a giant. For an extensive overview regarding the literature addressing 'Why do stars become giants?' we refer the reader to Sugimoto and Fujimoto (2000), while Faulkner (2005) provided a detailed analytical investigation inspired by the scientific legacy of Fred Hoyle.

Rotation
In the description of the internal structure of stars provided here, rotation has not been taken into account. However, it is plausible that all stars rotate as the clouds of gas and dust that they are formed from contain angular momentum. For slowly rotating stars it is generally assumed that second-order effects can be neglected and that the hydrostatic structure of the star is not affected. For low-mass stars on the main sequence and more evolved stars this will generally be the case. Intermediate-mass stars on and shortly after the main sequence may however be faster rotators before they slow down under the influence of magnetic braking. The rotation can have a strong impact on the thermal structure and radiative transfer, possibly inducing meridional flows. Additionally, the shape of the star may become aspherical, and mixing processes can alter. Work on including rotation in stellar structure models has been performed by for example Palacios et al (2006); Eggenberger et al ( , 2012. A critical issue is the treatment of the evolution of the internal angular velocity, including transport of angular momentum, which currently fails to reproduce the seismically inferred internal rotation (see Section 5.5). One of the main conclusions for red-giant stars so far is that rotationally induced mixing an meridional circulation do not provide enough mixing of chemicals to explain the abundance anomalies observed around the bump luminosity of globular clusters (see also Section 6).

Seismic diagnostics
Seismic diagnostics are by definition obtained from a signal that varies over time. In the context of oscillating stars, timeseries data (see top panel of Fig. 1) are most often taken from photometric fluxes or radial velocities (see Section 1). Some stellar parameters can be directly determined from the timeseries data. However, most seismic diagnostics are obtained from a Fourier transform (Figs 1, 6 and 7) of the timeseries data.
Important characteristics of the timeseries data are the total length or timespan of the data (T ) and the typical time sampling (δt). These translate in Fourier space in the frequency resolution δν = 1/T , and the highest frequency at which one can reliably obtain results, i.e. the Nyquist frequency ν Nyq = 1/(2δt). The Nyquist frequency is a hard limit for evenly sampled timeseries. However, astrophysical datasets are usually not exactly evenly sampled, which allows for measurements with higher frequencies, so-called super-Nyquist determinations (Murphy et al, 2013). It was shown by Eyer and Bartholdi (1999) that in cases of serious oversampling or undersampling the Nyquist frequency can be derived as ν Nyq = 1/(2p) with p being the greatest common divisor of all differences between consecutive observation times. In practise a realistic estimate of ν Nyq in the case of unevenly sampled data is to use the inverse of twice the median of all time differences between two consecutive measurements in the entire timeseries .
The power in Fourier transforms is commonly normalized using Parseval's theorem. This theorem states that the integral of the square of a function is equal to the integral of the square of its transform, i.e. the total power in the Fourier transform is equal to the total of the squared flux variations (for intensity) in the timeseries. Alternatively, Fourier transforms can be normalized using 'peak-scaling' in which the Fourier transform is normalized to recover the full sine-amplitude of an injected signal. The power can either be computed per bin, i.e. the frequency resolution, or per frequency unit. In the latter case it is the power density that is shown which has the advantage that its value does not depend on the frequency resolution.
The fact that integration times of the observations are not infinitely short causes apodization (η) of the power. This affects the power at all frequencies with the largest impact close to the Nyquist frequency. The apodization can be expressed as follows 4 : Gaps in the timeseries data impact on features in the power spectrum. This window function, i.e. the pattern of observations and gaps, can cause alias frequencies (or sidelobes) to occur at n(1/T gap ), with n an integer and T gap the typical time between gaps (for instance one day for ground-based single-site observations).
When the data quality between the observations varies significantly, such as can happen with ground-based spectroscopic (multi-site) campaigns, one can opt to compute a weighted Fourier transform. This weighting can be performed to optimize the noise level, but also to optimize the window function and reduce the sidelobes caused by gaps (Arentoft et al, 2009).
For a crash course on data analysis in asteroseismology including the statistics and uncertainties of timeseries data we refer the reader to Appourchaux (2014).

Variance and typical timescale in timeseries data
For stars with a convective outer layer the variations in the flux are dominated by granulation and oscillations. As both granulation, i.e. the visible pattern of convection, and oscillations depend on surface gravity it has been possible to calibrate relations between the variance in the flux and surface gravity Bastien et al, 2013) as well as between the typical timescales present in the timeseries data and surface gravity (Kallinger et al, 2016).

Background signal in Fourier spectrum
For stars with a convective outer layer a frequency-dependent background signal is present in the Fourier spectrum. This background consists of stellar intrinsic phenomena including activity features, such as spots and flares that are observable features of magnetic fields, rotation and granulation. In addition to intrinsic stellar background signal an observed power spectrum also includes white noise and instrumental effects. These instrumental effects can incorporate for instance degrading of CCDs but also incidental cosmic ray hits or telescope jitter. All these features together form a background on top of which the oscillation modes are visible as relatively narrow peaks. This background can be fitted with the following function (e.g. Harvey, 1985;Kallinger et al, 2014): where P ′ n is the white noise; a i , c i and b i are the rms amplitude, exponent of the i-th background component, and the frequency at which the power of the component is equal to half its value at zero frequency (the characteristic frequency), respectively. The factor η is the apodization defined in Eq. (1). The exponent c i provides a measure of the temporal correlation of the signal and determines the slope of the decay of the background in Fourier space (Mathur et al, 2011b). The number of functions (i) needed depends on the presence of activity, (super)granulation, or faculae (bright spots) on the star. Based on state-of the-art data, a two-component fit is in most cases necessary and sufficient (Karoff et al, 2013;Kallinger et al, 2014). An example of a background fit is shown in Fig. 6.

Oscillation signal
Oscillations establish themselves as a series of relatively narrow peaks on top of the background described in the previous subsection. The oscillations are confined to a limited range in frequency. Within this frequency range the individual modes of stochastic oscillations have a Lorentzian shape with a width that represents the lifetime of the mode and a height determined by the intrinsic visibility of the mode and geometrical effects. In the limit of infinite lifetime, i.e. an oscillation that appears coherent over the timespan of the observations, the mode is unresolved and takes the form of a sinc 2 function (see footnote 4) in power. Each oscillation mode is characterized by its quantum numbers: radial order n, related to the number of nodes in the radial direction (cf. Section 4.2.1); degree l, the number of nodal lines on the surface and azimuthal order m, the number of nodal lines crossing the stellar equator. The frequencies, width and amplitudes of the individual modes as well as the overall shape of, and patterns in, the oscillation power excess have valuable diagnostic power. Here, we first discuss the global features of the oscillation power excess, i.e. single measures that provide a diagnostic. We subsequently provide more details regarding the individual frequencies and the diagnostics that can be extracted from them.
In red giants all non-radial modes have a mixed character, with a gravity-mode behaviour (buoyancy is the restoring force) in the core and an acoustic behaviour (pressure is the restoring force) in the envelope. Observationally, the acoustic behaviour is most prominent, and we first discuss the related observed properties. Afterwards we consider the more profound aspects of the star that are revealed by the mixed nature of the modes.

Frequency of maximum oscillation power (ν max )
All solar-like oscillations in a star form together a bell-shaped power excess above the granulation and background signal centred around a specific frequency (see Fig. 7 where the amplitude of the power excess envelope is enhanced for visual purposes). This specific frequency is often referred to as frequency of maximum oscillation power or ν max . This frequency has been linked empirically to the acoustic cut-off frequency (Lamb, 1932, using the approximation for an isothermal atmosphere); here c is the adiabatic sound speed and H p is the pressure scale height. (See Eq. (25) and subsequent text for a theoretical explanation of the acoustic cut-off frequency.) It can be shown that ν max provides a direct measure of the surface gravity (g) when the effective temperature (T eff ) is known (e.g. Brown et al, 1991;: with M and R the stellar mass and radius, respectively. A theoretical basis for this relation has been investigated by Belkacem et al (2011) and a relation with the Mach number to the power three has been established, although this does not provide a full explanation.
The value of ν max can be estimated as the centroid of a Gaussian fit to the oscillation power excess . Alternatively, one can use the peak of the oscillation power excess in the smoothed power spectrum (Huber et al, 2009) or the first moment of the area under the smoothed power envelope (Hekker et al, 2010b). All these methods use slightly different, but equally valid, definitions of ν max and therefore can provide different values. Comparisons between values obtained with different methods show that this difference is generally within a few percent (Hekker et al, 2011Verner et al, 2011).
As stated above, for high-precision data such as the Kepler timeseries the main sources of the signal in the timeseries data are the granulation and the oscillations. The amplitudes of these signals are correlated with ν max , and therefore the frequency of maximum oscillation power can also be directly estimated from the variance in the timeseries .

Frequency pattern
Following asymptotic theory (Tassoul, 1980, and Section 4.1), acoustic oscillation modes (pressure as restoring force) of solar-like oscillators follow a distinct pattern: with ν cyclic oscillation frequency, ∆ν large frequency separation (Section 3.3.3), ǫ a phase term (Section 3.3.4) and d nl a small correction to the leading order asymptotics, which is zero for l = 0. Based on the asymptotic expression, Mosser et al (2011b) developed the universal pattern for red-giant spectra, according to which all parameters are assumed to be a function of ∆ν. This is equivalent to assuming that the underlying physics of the parameters varies as a function of the global stellar parameters. The universal pattern has the following form to describe pure acoustic modes (Mosser et al, 2011b(Mosser et al, , 2012a: where n max = ν max /∆ν. The phase term ǫ and non-radial correctiond l , withd 0 = 0, are described by scaling laws on the form A + B log ∆ν . The second order term (or curvature) in the asymptotic expression is represented by the parameter α = 0.015∆ν −0.32 (Mosser et al, 2012a).
The reason for the universal behaviour of red giants is not yet fully understood, as the asymptotic approximation is fundamentally related to the behaviour of the solution of the oscillation equations near the singularity at r = 0 (Gough, 1986b, 1993, andSection 4.2.3). In the case of red giants the conditions in the central region of the star are considerably different compared with the conditions in main-sequence stars (Section 4.2.3). Nevertheless, it seems that the asymptotic approximation and the universal pattern provide reasonable results for both models and observations. Only for the very luminous and expanded red giants close to the tip of the RGB the universality does not hold . For modes with low radial orders the dipole mode is located closer to the neighbouring quadrupole mode, providing a regular pattern of three modes together instead of a pattern of a cluster of 2 modes together (the 0-2 pair) alternating with one (dipole) mode.

Large frequency separation
The large frequency separation ∆ν is the separation between modes of the same degree and consecutive radial orders. ∆ν is proportional to the inverse of the acoustic diameter, i.e. the sound travel time across a stellar diameter. Furthermore, it can be shown that ∆ν is a direct probe of the mean density (ρ) of the star (Ulrich, 1986): with r the distance to the centre of the star, see also Eq. (31). The near-regular pattern of the large frequency separation can be measured in a global sense from the autocorrelation of a power spectrum (e.g., Huber et al, 2009), the power spectrum of a power spectrum (e.g., Hekker et al, 2010b), or the autocorrelation of the timeseries (EACF, Mosser and Appourchaux, 2009) with a cosine filter or Hanning function with a full-width at half maximum (FWHM) of the order of the FWHM of the power excess. The latter is mathematically equivalent to the power spectrum of the power spectrum. Additionally, Mosser et al (2011b) have used the universal pattern (sect. 3.3.2) in which the known patterns in the power spectrum are used to compute templates which are convolved with observed power spectra to determine ∆ν. Furthermore, the large frequency separation can be obtained from fits to individual frequencies (Kallinger et al, 2010), from pair-wise differences, or a linear fit of the frequencies versus their radial order. It has been shown that the different determinations are consistent within their uncertainties and definitions (e.g. Verner et al, 2011;Hekker et al, 2011Hekker et al, , 2012. Nevertheless, some biases depending on the number of radial orders that were used in the analysis have been identified . A convenient way to represent the power spectrum of solar-like oscillators is in anéchelle diagram (Grec et al, 1983), as shown for a subgiant star in Fig. 8. This is obtained by dividing the frequency spectrum into segments of length ∆ν and stacking the segments. According to Eq. (5) this should lead to roughly vertical sequences of points corresponding to different degrees.
We stress here that the actual separation between the frequencies of adjacent modes of the same degree varies as a function of frequency due to stellar internal structure properties, as also reflected by Eq. (6) and discussed further in Section 4.1. Structure variations that happen on scales that are comparable with or shorter than the oscillation wavelength, i.e. a glitch, cause a damped sinusoidal modulation in the frequencies (see for more details Sections 3.4.5 and 4.2.4). Structure changes that take place over longer scales cause a gradual change, or curvature, in the large frequency separation. Therefore the value of the large frequency separation may change depending on the frequency range that is taken into account.

Phase term (ǫ)
The asymptotic relation for acoustic modes (Eq. 5) also contains a phase term ǫ. This phase term is a dimensionless offset of the radial modes in anéchelle diagram (see Fig. 8). Its value is correlated with the determination of ∆ν. The value of ǫ can be determined from the universal pattern in which it is considered to be a function of ∆ν. Additionally, a (weighted 5 ) least-squares fit to the radial (l = 0) frequencies to simultaneously determine ∆ν and ǫ can be performed. These global methods average over the variation of ǫ with frequency. Additionally, a 'local' ǫ can be determined by only including the three central radial orders around ν max in the analysis .
The inferred value of ǫ is generally between 0.5 and 1.5, with a potential 'offset' of ±1. This 'offset' only reflects the observational limitations in that the radial order n cannot be measured independently and is not an offset of the actual value of ǫ.
The main diagnostic power of the global phase term lies in the mode identification of the different ridges in theéchelle diagram for stochastic oscillators with short mode lifetimes. In these cases the width of the oscillation signals in the Fourier space does not allow to resolve the small frequency separation (see next subsection) between the l = 0 and l = 2 modes and hence the l = 1 and l = 0, 2 ridge have very similar characteristics. In these cases ǫ, i.e. the location of the radial modes in anéchelle diagram, can be used to identify correctly the ridges .
The 'local' phase term can be used to distinguish between different evolutionary phases . This is caused by the fact that the differences in the core cause differences in the thermodynamic state of the envelope, which results in a different location of the second helium-ionisation zone for stars with an inert helium core compared with stars with helium-core burning. The location of the second helium-ionisation zone leaves a trace in small oscillatory deviations in the frequencies (see Sections 3.4.5 and 4.2.4 on glitches). This causes a difference in the 'local' phase term for stars in different evolutionary phases (Christensen-Dalsgaard et al, 2014).

Small frequency separation (odd and even degree modes)
A typical separation in frequency exists between odd-and even-degree modes. This is the so-called small frequency separation that can be approximated asymptotically, for main-sequence stars, by (e.g. Gough, 1986b, see also Eq. 32): This parameter is generally measured from frequency differences between observed individual frequencies. Eq. (8) is rather strongly weighted towards the central parts of the star. Hence, δν l l+2 (n) provides a measure of the helium content in the core of main-sequence stars, and with that δν l l+2 (n) is a diagnostic of stellar age (Christensen-Dalsgaard, 1984;Ulrich, 1986;Christensen-Dalsgaard, 1988). It was noted by Huber et al (2010) from early Kepler observations, that for more evolved stars (subgiants and giants) the small frequency separation between modes of degree 0 and 2 essentially scales as ∆ν. This is due to the fact that these stars have a highly concentrated cores such that the lower turning point of the pressure mode cavity lies outside the compact core. Therefore, this small separation does not provide a measure of the helium content in the core and hence is no longer an age diagnostic.

Small frequency separation (δν 01 )
In a Fourier spectrum dipole modes are located approximately mid-way between radial modes as per Eq. (5). The offset from the midpoint between the radial modes and the frequency of the dipole mode can be computed from a three-point difference 6 and indicated as δν 0 1 : In main-sequence stars δν 0 1 is known to be sensitive to the central physical conditions. For red giants this is not the case. It has been shown that δν 0 1 is correlated with the distance between the l = 1 turning point and the bottom of the convective envelope. The value of δν 0 1 takes small and negative values for stars ascending and descending the RGB where the turning points of acoustic l = 1 modes are well within the convective envelope. Stars in the He-core burning phase have a shallower convective envelope and the turning points of l = 1 modes are in the radiative region and δν 0 1 generally takes positive values . As shown in Section 4.2.5 (Figs 32 and 33) the evanescent region between the buoyancy and acoustic cavities is smaller in a red-clump star compared with a RGB star. This leads to a stronger coupling between modes in a red-clump star increasing the spread in the modes, which most likely leads to less regular behaviour of δν 0 1 in clump stars compared with RGB stars. Hence, the value of δν 0 1 and the regularity of the acoustic dipole spectrum can be used as a diagnostic to distinguish between different evolutionary phases . We note here that most of the work on δν 0 1 for red giants was performed before the discovery of the fact that all non-radial modes are to some extent mixed modes. Hence, this influence was not taken into account.

Period spacing (∆Π)
As mentioned above, non-radial modes in red giants all involve aspects of gravitymode behaviour. To analyse the relevant properties we note that an asymptotic approximation (see Section 4.1) reveals that gravity modes of a given degree appear with near constant separation in period, the so-called period spacing ∆Π, as reflected in the asymptotic behaviour of their periods Π n,l : with where N is the Bruntt-Väsälä frequency (Eq. 24), r 1 and r 2 the turning points of the gravity mode cavity and ǫ g is a phase accounting for the behaviour near the turning points. This defines the period spacing: Based on the regular nature, one can use similar techniques to determine the period spacing as used to determine ∆ν (power spectrum of the power spectrum, auto-correlations and differences between modes with consecutive order, see Section 3.3.3), but now in period space and centred on regions where dipole (or quadrupole) modes are expected. It is however important to account for the fact that in solarlike oscillators no pure gravity modes can be observed, but only modes that have mixed gravity-pressure nature (see also Section 3.4 and Section 4.2 for more details about mixed modes). Due to this mixed nature the regular spacing in period deviates from the asymptotically predicted value of 'pure' gravity oscillation modes in the vicinity of 'pure' pressure modes. The observed period spacing follows, in general, a predictable pattern for each acoustic radial order with an empirically determined Lorentzian shape (Mosser et al, 2011b(Mosser et al, , 2012c, and Section 4.2.3).
The underlying period spacing can be determined from the modulated observed period spacing using a periodéchelle diagram (see Fig. 9). In this diagram the period spacing can be obtained by aligning the g-dominated modes surrounded in a roughly symmetric way by the p-dominated modes Mosser et al, 2012c;Datta et al, 2015). The patterns in the periodéchelle diagram can also be described by an analytical expression based on the coupling between the gravity and acoustic cavity as derived by Mosser et al (2012c) (see also Section 4.2.3): where q is the coupling strength. This is set to q = 0 for no coupling and q = 1 indicates maximum coupling. The period spacing and coupling strength can be determined in an iterative manner, when assuming ǫ g is a fixed value (often zero, i.e. the pattern is assumed symmetric). Note that Buysschaert et al (2016) have employed Eq. (13) while leaving ǫ g also as a free parameter. They find that although ǫ g remains ill-defined, its determination improves the determination of the period spacing. An alternative way to determine period spacing and to constrain the evanescent region between the p and g cavities (and thus coupling strength q) is through the inertia ratio of dipole and radial modes (see Section 4.2). Benomar et al (2014) estimated mode inertias observationally from the measurements of mixed mode characteristics for the first time. The mode inertia ratio could develop to be a diagnostic that can potentially provide strong constraints on the stellar structure.
Recently, Mosser et al (2015); Vrard et al (2016) have proposed a way to stretch the period Fourier spectrum to remove the modulation due to the coupling between the pressure and gravity modes and obtain the regular underlying period spacing in a direct manner.
The period spacing provides a strong diagnostic on the central regions of the star. It can be used to distinguish between stars in different evolutionary phases, most notably between red giants with an inert He-core and red giants with He-core burning Mosser et al, 2011aMosser et al, , 2014. Furthermore, small deviations from the regular period spacing pattern can reveal localised stellar structure changes such as the chemical discontinuity due to the first dredge-up , and Section 4.2.4).

Individual oscillation modes
Individual oscillations modes in stochastic oscillations are characterized by their mode frequency ν central , line width Γ and height H. A resolved oscillation mode, i.e. a mode with a lifetime that is at least ∼10 times shorter than the timespan of the timeseries data (Hekker et al, 2010a) can be modelled by a Lorentzian profile in the power spectrum P(ν): The mode lifetime t damp is directly related to the mode line width through Γ = 1/(πt damp ). The mode height and width are highly correlated and relate through the root-mean-square flux amplitude A = √ πHΓ/2 (e.g., Chaplin et al, 2005) which is the area underneath the profile. Hence the amplitude is a more robust parameter. The amplitude of the modes contains information on the excitation and damping of the oscillations (see e.g. Houdek, 2012;Samadi et al, 2012, and Section 4.3). However, the visibility of the modes is a combination of the intrinsic amplitude and geometrical effects, i.e., cancellation of some of the signal due to the fact that only integrated light from the visible part of stars can be observed. In case the oscillations are not resolved (lifetime longer than the timespan of the timeseries data) they are often approximated with a sinc function (see footnote 4) in which case the amplitude is A = √ 2Hδν, where δν is the frequency resolution. In intermediate cases where the modes are partly resolved a mixture of the resolved and unresolved description has to be applied (see Section 4.3.2).

Mixed modes
Radial modes are always pure acoustic modes with pressure as the restoring force. Non-radial modes in red giants, however, always have a mixed nature, i.e. are mixed modes, for which buoyancy acts as restoring force in the deep interior of the star and pressure acts as restoring force in the outer layers of the star. In other words a mixed mode is a single mode with different behaviour in the different regions. Frequencies of mixed modes are shifted compared with pure acoustic or gravity modes by an amount depending on the coupling strength between the two (gravity and pressure) oscillation cavities (e.g. Deheuvels and Michel, 2010;Hekker and Mazumdar, 2014, and Section 4.2).
Mixed modes (and hence period spacings (Section 3.3.7)) are mostly observed in dipole (l = 1) modes as for these modes the coupling between the pressure and gravity mode cavity is stronger (narrower evanescent zone between the cavities, see Fig. 12) and also because the period spacings are larger due to the dependence on √ l(l + 1) (see Eq. 10), and thus better resolved.
In case there are only a few mixed modes present in for instance subgiants, Michel (2010, 2011) showed that l = 1 avoided crossings in subgiant stars involve more than two modes and induce a characteristic distortion in the l = 1 ridge in theéchelle diagram (see Fig. 16). This can be used to constrain stellar models. Deheuvels and Michel (2011) have done so by matching the observed large frequency separation and frequency of the avoided crossing with values obtained from models. This results in a precise age estimate given the mass and physics of the models.
Following the analysis by Deheuvels and Michel (2011), Bedding (2012) introduced a powerful tool for analysing these mixed modes. He showed that by replicating theéchelle diagram horizontally the full structure of the avoided crossings can be displayed. Benomar et al (2012) subsequently used this to develop a method to fit the avoided crossings and determine the minimum separation between the two branches as a measure of the coupling strength (see Section 4.2.2). This provides a useful diagnostics of the stellar mass for subgiant stars.

Mode identification
The degree l of the observed individual frequencies can be obtained from the known pattern of the stochastic oscillation (Section 3.3.2) through anéchelle diagram (Fig. 8) where the ridges of radial and quadrupole modes are close together with the dipole ridge appearing at about a 0.5∆ν offset. In case the mode lifetimes are very short the radial modes can be so broad that they overlap the quadrupole modes. In that case the phase term (ǫ) can be useful to distinguish the odd and even ridge (Section 3.3.4). The (acoustic) radial order n p can be estimated from Eq. (5), i.e., from the ratio of the frequency of the mode over the large frequency separation bearing in mind that n is an integer and that the phase term ǫ takes values between 0.5 and 1.5. Note that for mixed modes the total radial order n consists of the nodes in the acoustic cavity n p as defined here and the nodes in the buoyancy cavity. The buoyancy radial order n g is indicated by definition with a negative number and can take large values (see also Section 4.2.1 and Fig. 18).
Mode identification can also be performed from spectroscopic data. Red-or blueshifted parts of the surface of a star leave traces in the shape of a spectral line profile. Over the course of the pulsation the blue-and red-shifted parts change and hence the line-profile shape changes. Therefore, the amplitude and phase of the line-profile variations at a particular frequency are fundamentally different for radial and nonradial oscillation modes. This technique has mainly been developed for stars with coherent oscillations (Zima et al, 2004;Zima, 2008), but also proved useful to reveal for the first time the presence of non-radial oscillations in red giants (Hekker et al, 2006;Hekker and Aerts, 2010).

Surface effect
Due to incomplete modelling of the convective outer layers of stars and the strong non-adiabatic behaviour of the oscillations in the superficial layers there is an offset between modelled and observed frequencies, the so-called surface effect. This offset is a function of frequency, but independent of degree at least for the Sun. Using the solar offset Kjeldsen et al (2008) proposed a widely used power-law correction that can be applied to other stars. An alternative procedure is to directly scale the solar offset on a suitable frequency scale (e.g., Christensen-Dalsgaard, 2012b). Additionally, Aerts et al (2010) have shown that due to their larger mode mass in the inner regions, mixed modes in subgiants are less sensitive to the incorrect modelling of the non-adiabatic outer layers, and hence the solar calibrated offset needs to be adapted. Ball and Gizon (2014), based on Gough (1990), have developed a correction based on mode inertia. So far this method has only been tested for observations of the Sun and solar-like stars as well as for models of different mass and metallicity covering a significant portion of the HR diagram from the main-sequence to red giants for models with T eff < 6500 K (Schmitt and Basu, 2015). In their work Schmitt and Basu (2015) concluded that the two-term model proposaed by Ball and Gizon (2014) works much better than other models across a large portion of the HR diagram, including the red giants.
For main-sequence stars it is possible to eliminate the influence of the nearsurface region by using frequency-separation ratios such as which are essentially independent of the near-surface problems (e.g. Roxburgh and Vorontsov, 2003;Roxburgh, 2005). As for the small frequency separations (δν 02 , Section 3.3.5) the usefulness of the frequency-separation ratios for red-giant stars is limited due to the universality of the frequency patterns (e.g. Huber et al, 2010). Additionally, in giants all non-radial modes are mixed modes, which have different sensitivities to the surface effect.

Rotational splitting
Rotation splits the non-radial modes into 2l + 1 singlets of different azimuthal orders m. These incorporate modes travelling with the rotational direction (prograde modes) and in the opposite direction (retrograde modes) in addition to the original mode unperturbed by rotation, i.e., in cyclic frequency (see also Eq. (123) given in angular frequency): In cases of slow rotation, which is in general the case for subgiants and red-giant stars, the assumption of symmetric splittings is often valid. The relative heights or visibility of the different azimuthal orders in a multiplet are indicative of the inclination angle with respect to the rotation axis at which we view the system (Gizon and Solanki, 2003). When viewing the system pole on, only the m = 0 mode is visible, while for an equator-on system all modes with even l − m are visible. A rotationally split mode can be fitted with a set of Lorentzian functions: where Ψ l m (i) is the visibility of the mode which depends on the inclination angle (i) and δν n l m is the rotational splitting. Here we assumed that all modes in a multiplet are excited to the same intrinsic average height, as may be reasonable for stochastically excited modes observed for a long time compared with the lifetime of the modes.
To detect the average rotational splitting one can also use the fact that the rotational splitting is approximately symmetric in slow rotators and apply for instance the EACF (autocorrelation of the timeseries, Mosser et al, 2011b) with a very narrow Hanning filter in the range of a non-radial rotationally split modes. It has been shown that the rotational splitting of modes with mixed character (i.e. all non-radial modes in red-giant stars) is a function of the mixed character (Beck et al, 2012;Mosser et al, 2012b;Goupil et al, 2013, and Section 4.4). Hence the rotational splitting of mixed modes is a diagnostic to probe the radial internal rotation profile.

Glitches
Information on specific transition regions in a star, such as the boundaries of convective zones or ionization zones of helium or hydrogen, can be obtained from the fact that, at such boundaries, the properties of the star change on a scale substantially smaller than the local wavelength of the oscillations (e.g. Vorontsov, 1988;Gough, 1990, and Section 4.2.4). These sharp features (also called glitches) cause oscillatory variations in the frequencies with respect to the pattern described in Eq. (5), or the corresponding pattern satisfied by mixed modes (cf. Eq. 13). The period of the variation depends on the location of the feature, while the amplitude depends on the detailed properties of the feature. For acoustic modes the oscillatory behaviour due to glitches can be measured in frequency differences (∆ν) but more often in second differences, and can be described by a damped oscillator: in the case of the glitch due to the He II ionisation zone, where c 0 indicates the amplitude of the oscillation and c 2 a characteristic width of the e-folding time of the damped oscillator. The acoustic depth is indicated by τ gl and φ gl is a constant that accounts for the phase. Eq. (18) was derived by Houdek and Gough (2007) and applied in various forms by e.g. Miglio et al (2010); Mazumdar et al (2014); Broomhall et al (2014); Verma et al (2014). In addition to the intrinsic limits of the sharpness of the glitch compared with the local wavelength of the oscillation mode (see Section 4.2.4), the observational data also provide natural limits. The depth in the star at which a glitch can be measured depends on the frequency range covered by the oscillations, i.e. the largest period in Eq. (18) that can be measured, while the minimum period is defined by the resolution of the measured frequencies. From these limitations we find that the bottom of the convection zone is located too deep in red giants to be measured. Furthermore, the helium I and hydrogen ionisation zones are located close to the surface, which make them very challenging to measure. Cunha et al (2015) and Section 4.2.4 provided a detailed discussion of the effect of buoyancy glitches on mixed modes, although no solid observational results on this have been presented so far.
Note that the diagnostic use of glitches can be done completely independent of stellar models.

Scaling relations and grid-based modelling
The scaling relations Eqs (4) and (7) can be used to obtain the mean density and surface gravity (and from that stellar mass and radius) in a direct manner, i.e. so-called direct method. These scaling relations are exceptionally good given that these relations do not account for metallicity differences, nor do they account for any knowledge we have about stellar evolution. To take account of this knowledge, it is also possible to compare the observables {∆ν, ν max , T eff , [Fe/H], π, L} or a subset thereof with a grid of models. In this so-called grid-based modelling one does account for knowledge of stellar structure and evolution, as well as metallicity. However, in case the scaling relations are used to determine ∆ν and ν max from the models suitable reference values (with uncertainties) to which one scales have to be adopted. Alternatively, one can compute individual frequencies for the models and derive ∆ν from that.
Reference values Both the direct method and grid-based modelling are based on the scaling relations (Eqs 4 and 7) which assume that the scaling is valid in a consistent way between the reference and the observed star. Often reference values based on the Sun are used; however, this may not be correct for stars with different properties, such as a different metallicity or rotation rate, or stars in different evolution phases, as their stellar internal structures are different. This has indeed been confirmed by White et al (2011) for stars with different metallicities along the mainsequence and just beyond that. Mosser et al (2013b) proposed new reference values based on empirically derived relations using results obtained using the universal pattern (Mosser et al, 2011b). It is currently unclear whether these newly derived reference values are also valid when other methods are used to derive ∆ν and ν max . For instance, Hekker et al (2013) failed to confirm the empirical relations quantitatively using stellar models. Furthermore, Miglio et al (2012) showed that the difference in internal temperature structure (hence sound speed) between RGB and RC stars has a significant impact on the scaling relations. Therefore, they applied a correction for red-clump stars in the open clusters NGC 6719 and NGC 6819 based on the masses of the RGB stars in the respective clusters and theoretical models. Subsequently, Christensen-Dalsgaard et al (2014) showed that the difference in the variation of the phase term ǫ (see Section 3.3.4) with frequency between RC and RGB stars is related to differences in the thermodynamic state of the convection zone. This supports the findings of Miglio et al (2012) that RC and RGB stars have internal structures that are significantly different, which calls for corrections to the scaling relations. Recently, Guggenberger et al (2016) proposed a new reference function for the ∆ν scaling relation (Eq. 7) that accounts for metallicity differences and is applicable for stars on the main sequence up till past the RGB bump on the red-giant branch. Additionally, Sharma et al (2016) and Serrenelli et al. (in preparation) devised methods to apply a correction between ∆ν scaling and ∆ν from frequencies in model calculations.

'Boutique' modelling
Oscillation frequencies can be compared with stellar models to infer the internal stellar structure. This is often done on a star by star basis, hence 'boutique' or 'detailed' modelling. This modelling can be done in both a forward and an inverse approach.
Forward modelling In forward modelling the observed frequencies and stellar parameters (log g, T eff , [Fe/H]) are matched with stellar models. This is mostly done using χ 2 minimisation either in a direct manner or through Singular Value Decomposition (SVD) (e.g. Brown et al, 1994). In this procedure the surface term is accounted for by one of the prescriptions mentioned above Ball and Gizon, 2014), or by using the frequency-separation ratios (Roxburgh and Vorontsov, 2003;Roxburgh, 2005). Detailed descriptions of the different methods currently applied can be found in the appendix of Ballard et al (2014); Silva Aguirre et al (2015).
Inverse modelling An inverse problem is a general framework that is used to convert observed measurements into information about a physical object or system. For stars the individual frequencies can be used to obtain information about the internal structure of stars. An inverse problem is however by nature ill-posed, and many frequencies probing the star to different depths are needed to obtain meaningful results. Most commonly inversions have been used to study the stellar rotation at different depths in the stars (e.g. Deheuvels et al, 2012Deheuvels et al, , 2014Di Mauro et al, 2016). It is however in principle also possible to perform structure inversions. These have been very powerful in determining the internal structure of the Sun (e.g., Gough et al, 1996, and references therein). For other stars the lack of observed frequencies probing the star to different depth has precluded detailed structure inversions. However, inversion techniques have been used in asteroseismic analyses to constrain specific properties of stars (e.g. Reese et al, 2012;Buldgen et al, 2016).
For a complete overview of stellar inversions we refer to Basu (2014Basu ( , 2016

Uncertainties
A result for any of the diagnostics mentioned above can only be meaningful when it is accompanied by an uncertainty. Uncertainties can be obtained using different statistics. The two main streams are a frequentist approach and a Bayesian approach.
In a frequentist approach one assigns probability to measurable events that can be measured an infinite number of times where the number of parameters remains fixed. Maximum likelihood estimation is a typical example of a frequentist approach. In a Bayesian approach the data are taken as fixed and are described with parameters that are unknown and probabilistic. Among other attributes, this generally means that a Bayesian analysis provides a more realistic uncertainty. See Appourchaux (2014) for detailed descriptions and formulations of both frequentist and Bayesian uncertainty calculations.

Theory of stellar pulsations
The general theory of stellar pulsations has been presented in considerable detail by, for example, Unno et al (1989) and Aerts et al (2010). However, oscillations of evolved stars present special properties which are important for the understanding of the observations. Thus, here we provide some background which is useful in the interpretation of the observations of the oscillations of such stars, relating the frequencies and other aspects of the oscillations to the properties of the stars. We are dealing with low-amplitude oscillations, which can be regarded as small perturbations to the equilibrium structure. Formally, these can be described using linearized perturbation analysis of the general equations of hydrodynamics. An important result concerns the geometrical properties of the modes of spherically symmetric stars. For the modes that are relevant the properties can be described by spherical harmonics Y m l (θ, φ) = P m l (cos θ)e imφ as functions of co-latitude θ and longitude φ; here P m l is a Legendre function, characterized by the degree l and the azimuthal order m (see Sect 3.3 for a definition). The time dependence of a mode is conveniently written as e −iωt where ω is the angular frequency, which is in general complex. The displacement vector can be written, as a function of position r and time t, where Re denotes the real part, ξ r and ξ h are the radial-and horizontal-displacement amplitude functions that depend only on the distance r to the centre, and a r , a θ and a φ are unit vectors in a spherical polar coordinate system. Other oscillating variables, such as the pressure perturbation, vary as the real part of Y m l (θ, φ)e −iωt . We can separate the frequency into real and imaginary parts as ω = ω r +iω i . Then the dependence of the oscillations on longitude φ and time t is essentially Unless m = 0 this describes a wave running in the φ direction, growing or decaying with time depending on whether ω i is positive or negative. In much of the following we consider adiabatic oscillations where processes causing excitation or damping are neglected. Then the frequency ω is real and we ignore the distinction between ω and ω r . We note that observed oscillations are typically discussed in terms of the cyclic frequency ν = ω r /2π, as done in the previous sections. However, for the theoretical analysis it is more convenient to use the angular frequency ω.
After this separation of variables we are left with differential equations that depend just on r. Combined with suitable boundary conditions, this is a relatively straightforward numerical problem, which determines the frequencies ω as eigenvalues. However, the physical treatment of the near-surface layers still suffers from substantial uncertainties, particularly when the mode energetics is taken into account (see Section 3.4.3).
To evaluate the diagnostic potential of solar-like oscillations in giant stars and interpret the inferences that are made, an understanding of the properties of the oscillations is required. Moreover, these properties are fascinating in their own right. A full utilization of the observed data requires detailed comparison of the observed frequencies, and other properties, with computations for stellar models. However, a great deal of insight as well as powerful diagnostics can be obtained from asymptotic analyses of the oscillations, to which we turn next.

Asymptotic theory
We first concentrate on the oscillation frequencies and overall properties of the eigenfunctions, assuming the oscillations to be adiabatic, and analyse the behaviour of the modes and their relation to stellar structure in terms of their asymptotic properties. Solar-like oscillations are generally of high radial order. Hence it is common to ignore the perturbation to the gravitational potential, in the so-called Cowling approximation (Cowling, 1941). Dziembowski (2012) pointed out, however, that this approximation is questionable for dipolar modes (with l = 1). Here the perturbation to the gravitational potential gives rise to a slowly varying component of the solution which may have a significant effect on the properties of the modes. The effects of this on the analysis discussed below deserves further investigation.
The Cowling approximation reduces the equations of adiabatic oscillations to a second-order system, greatly simplifying the analysis. The equations are often expressed on the form for suitable choices of the dependent variable X and the function K. The choice of X and K may depend on the specific properties that are being investigated. A convenient formulation was derived by Gough (Deubner and Gough, 1984;Gough, 1993) based on an analysis by Lamb (1932). Here X = c 2 ρ 1/2 div δ δ δr, where c is the adiabatic sound speed and ρ is density. The corresponding approximation to K is which is determined by three characteristic frequencies of the star: -The Lamb frequency S l , with The Lamb frequency is a local characteristic frequency of horizontally-propagating sound waves with a wavenumber k h = √ l(l + 1)/r.
-The buoyancy frequency (or Brunt-Väisälä frequency) N, where g is the local gravitational acceleration, p is pressure and Γ 1 = (∂ ln p/∂ ln ρ) ad , the derivative being for an adiabatic process; note that N 2 is negative in convectively unstable regions. The Brunt-Väisälä frequency is the local frequency of internal gravity waves of short horizontal wavelength.
-The acoustic cut-off frequency ω c , where H = −(d ln ρ/dr) −1 is the density scale height. The acoustic cut-off frequency arises from the inability of modes to propagate when their vertical wavelength is too long compared with the scale of the density variation in the equilibrium structure. This leads to reflection of the waves. We note that since ω 2 c depends on the second derivative of density (see the definition of density scale height and Eq. 25), it shows rapid variations in the region of substantial superadiabaticity just below the surface. This is visible in its curious behaviour in the right-hand panel of Figure 10. On the other hand, in the nearly isothermal atmosphere H ≃ H p is essentially constant, and ω 2 c ≃ c 2 /(4H 2 ). Thus ω c /2π reduces to ν ac (cf. Eq. 3).
According to Equations (21) and (22) the behaviour of the oscillations is determined by the dependence of the characteristic frequencies on position. This is illustrated in Figure 10 for a model approximating the present Sun. To characterize the location of typical frequencies of solar-like oscillations the horizontal line shows the estimated frequency ν max of maximum power, obtained as where ν ac is the isothermal acoustic cut-off frequency (cf. Eq. 3); in the Sun ν max ≃ 3150 µHz. At this frequency ω ≫ N except in the atmosphere, reducing Eq. (22) to . The corresponding modes are acoustic modes, or p modes, where pressure is the restoring force. In the region where K > 0 the solution X oscillates as a function of r (see Eq. 21), whereas X behaves locally exponentially where K < 0. For acoustic oscillations the oscillatory region extends to a distance r t from the centre, approximately given by ω ≃ S l , or Physically this corresponds to total internal reflection of the oscillations, described as a superposition of sound waves. The upper turning point R t is where ω ≃ ω c (R t ), which is satisfied just below the photosphere at high frequency and somewhat deeper at lower frequency (see left panel of Fig. 10). Waves with frequencies exceeding the acoustic cut-off frequency in the atmosphere are free to travel outwards in the atmosphere, resulting in strong damping. A more quantitative analysis of Eq. (21) can be carried out using JWKB theory (see Gough, 2007). In general, this results in an eigenvalue condition on ω given by for integer k, where r 1 and r 2 are adjacent turning points at which K = 0, such that K > 0 between r 1 and r 2 . In the present case, for predominantly acoustic modes in main-sequence stars, this leads to For low-degree modes, r t is close to the centre in main-sequence stars. Equation (29) can be reduced to using an expansion around r = 0 as well as near the surface (Gough, 1986a(Gough, , 1993. For main-sequence stars k can in general be related directly to the radial order n of the mode which was used in Eq. (30). For evolved stars the definition of mode order is more complex; we return to this in Section 4.2.1. Eq. (30) is the basis for the different frequency separations and the phase term ǫ. As discussed in Sections 3.3.3 -3.3.6 these quantities provide important diagnostics of stars based on their acoustic oscillations. 7 From homology scaling c 2 ∝ M/R. It follows from Eq. (31) that ∆ν, and hence from the leading-order first term in Eq. (30) ν nl , scale as From Eq. (23) it follows that the same scaling applies to S l . To investigate the changes in oscillation properties as a star evolves from the main sequence through the subgiant phase to the red-giant branch we consider a 1.3 M ⊙ evolution sequence from Jiang and Christensen-Dalsgaard (2014). The evolution in the HR diagram is illustrated in Fig. 11. The behaviour of the Lamb frequency remains similar as the star evolves, apart from the scaling withρ 1/2 . However, the behaviour of the Brunt-Väisälä frequency is very different for red-giant models compared with main-sequence models. This can be seen by approximating N 2 as where ∇ = d ln T/d ln p, ∇ ad is its adiabatic value and ∇ µ = d ln µ/d ln p, with µ being the mean molecular weight. In a red giant with a compact core g reaches very high values in the deep interior of the star and so therefore does N. This is illustrated  in Fig. 12 (note the logarithmic abscissa). Given the larger radius, both S l and ν max are substantially reduced. However, the most dramatic difference compared with Fig.  10 is the very large value of N in the core. Additional features in N are caused by the composition discontinuity at r ≃ 0.09R left behind after the first dredge-up (see Sect. 2.1.2) and local maximum at r ≃ 0.005R arising from the steep abundance gradient in the hydrogen-burning shell. Note also the very deep convective envelope, where N is imaginary and hence not shown.
For ascending-branch red-giant models such as the one shown in Fig. 12 there are two regions where K > 0, leading to an oscillatory behaviour of the eigenfunctions, i.e., with ω > S l , N or ω < S l , N. Of these, the outer region (ω > S l , N) corresponds essentially to the oscillatory region in main-sequence stars discussed above. Modes trapped in this regions, with the eigenfunction decreasing exponentially below it, satisfy Eq. (29). We note, however, that the analysis leading to the asymptotic approximation for low-degree modes cannot immediately be transferred to more evolved stars, given that the lower turning points lie outside the compact core (see also Section 4.2.3). Even so, both model computations and observations show that the acoustic modes of red giants satisfy a relation very similar to Eq. (30), i.e., the 'universal pattern'. This was discussed in Sect. 3.3.2, where departures from this pattern were also mentioned. Indeed, Dziembowski (2012) found from stellar models that for the most luminous red-giant branch stars the dipolar mode frequencies are shifted substantially relative to the location at the mid-point between the neighbouring radial-mode frequencies which is predicted by the leading-order term in Eq. (29). This was confirmed observationally by Stello et al (2014).
The inner oscillatory region (ω < S l , N) is in the core of the model, where the frequency is below the Brunt-Väisälä frequency. The modes trapped in this region are standing internal gravity waves, g modes, where buoyancy is the restoring force. Their frequencies may be estimated from Eq. (28), with r 1 and r 2 being approximately the points where ω = N. Approximating K by assuming that ω ≪ S l and neglecting ω c we obtain In most of the region ω ≪ N. With a correction for the behaviour near the turning points we can therefore approximate this further. We obtain a relation for the period Π = 2π/ω: where and ǫ g is a phase term accounting for the behaviour near the turning points. 8 Thus for modes trapped in the core of the model we obtain oscillations with uniformly spaced periods. These periods increase with increasing k, and with a spacing Π k+1 l − Π k l ≃ ∆Π l = Π 0 / √ l(l + 1) which depends on the degree of the modes. The diagnostic power of the period spacing is summarized in Section 3.3.7. As discussed in Section 4.2.1 g modes are by convention assigned negative radial orders n g , with frequency tending to zero as n g tends to −∞. Thus at least in simple cases of a pure g-mode spectrum k in Eq. (37) can be identified with |n g |. We return to a more complete discussion of mode order in Section 4.2.1 and beyond.
In the analysis leading to Eqs (29) and (35) we assumed that the modes were completely trapped in the corresponding regions, i.e., acoustic modes with pressure as the restoring force (p modes) in the outer part of the star and gravity (g) modes with buoyancy as the restoring force in the core. However, it is clear from Figure 12 that the evanescent region separating the two trapping regions is quite thin, particularly for l = 1 modes. This leads to substantial coupling between the two regions and hence generally to a mixed character of the modes. These mixed modes are responsible for the diagnostic richness of the solar-like oscillations in evolved stars. We discuss this in the following subsection.

Mixed modes
The first to consider non-radial mixed modes in highly-evolved stars was likely Dziembowski (1971), who analysed the oscillations of Cepheid-type stars. Dziembowski noted that the huge values of the buoyancy frequency in the core of such stars meant that even at high frequencies the modes behaved as standing internal gravity waves in the inner parts of the star. Dziembowski carried out an asymptotic analysis of these properties. This was subsequently followed by a detailed investigation concerning the effects of non-adiabatic properties of the modes by Dziembowski (1977). Scuflaire (1974) analysed the oscillation properties of polytropes of high polytropic index and hence centrally condensed models. Scuflaire similarly noticed the mixed character of the modes and may have been the first to use explicitly the term 'mixed modes'. Osaki (1975) followed the evolution of mixed modes in a massive main-sequence model. He noticed the effect of the increasing frequencies of the gravity waves in the core that led to the characteristic behaviour of the model frequencies with age (see Fig.  13).
The overall evolution of the frequencies of mixed modes is illustrated in Fig. 13, showing results for a 1.3 M ⊙ evolution sequence extending from just after central hydrogen exhaustion to the early part of the red-giant ascent (cf. Fig. 11). To eliminate the dependence of the acoustic-mode frequencies on stellar radius the frequencies have been scaled with (R/R 0 ) 3/2 (cf. Eq. 7), where R 0 is the zero-age main-sequence radius, such that the modes following the acoustic scaling relation, Eq. (33), appear with constant frequency. This is approximately the case for the radial modes (dashed lines in Fig. 13), and for l = 1 (solid lines in Fig. 13) for those modes that are predominantly of acoustic nature. However, there is clearly a second class of modes with scaled frequencies increasing with age. These modes are predominantly of gmode character. Their frequencies increase with age following the strong increase in the buoyancy frequency caused by the contraction of the core (cf. Eq. 34). Where the frequency of such a mode meets a acoustically dominated mode the frequencies do not cross but approach quite closely. This is followed by an exchange of mode character such that the mode previously of g-mode nature becomes predominantly acoustic and vice versa. This behaviour was first noticed for evolving stars by Osaki (1975). It was also investigated by Aizenman et al (1977), who were probably the first to use the term 'avoided crossing' 9 to describe it in astrophysics. Note Fig. 13 Evolution of frequencies for a 1.3 M ⊙ stellar model, from the beginning of the subgiant phase through the early red-giant phase (marked by the red part of the line in Fig. 11; see Sect. 2.1.2). The frequencies have been scaled by (R/R 0 ) 3/2 , where R 0 is the ZAMS radius, to correct for the evolution with mean density of the acoustic-mode frequencies (cf. Eq. 7). The heavy long-dashed line shows an estimate of ν max (cf. Eq. 26), similarly scaled. The dashed lines show radial modes while the solid lines show dipolar modes, exhibiting avoided crossings first mentioned by Aizenman et al (1977). The lower panel shows a small part, marked by a red rectangle, of the upper panel, to illustrate the behaviour in more detail. Adapted from Jiang and Christensen-Dalsgaard (2014). however that similar phenomena have wide applicability to cases of coupled oscillations, including in atomic physics. An early illustrative analysis was provided by von Neuman and Wigner (1929).
Even though the actual frequencies do not cross it is sometimes useful to relate them to fictitious uncoupled gravity and acoustic modes. These are the so-called 'γ' and 'π' modes, respectively, that do cross (see Aizenman et al, 1977). As discussed below the behaviour of the actual modes can then be analysed by introducing coupling between these fictitious modes. This also allows the definition of the numbers N π and N γ that represent the number of π and γ modes in the relevant frequency interval, e.g., corresponding to the range of radial modes (e.g., Benomar et al, 2013). In the earlier phases of evolution N γ ≪ N π , whereas on the red-giant branch (towards the right edge of Fig. 13) N π ≪ N γ . This variation in the overall structure of the oscillation spectrum has a major effect on the observed Fourier spectra (see Fig. 7).
To characterize the relative contributions of the different regions of the star to a mixed mode, a useful quantity is normalized inertia E. This is defined by where the integral is over the volume of the star, and the normalization uses the average squared photospheric displacement. In addition to E we also consider mode mass M mode = ME, where M is the mass of the star. These quantities are defined such that the average kinetic energy of the oscillation is where V rms is the average surface velocity. For p-dominated modes E is largely a function of frequency. This frequency dependence is determined by the depth of the upper turning point R t , and hence the decrease in amplitude of the eigenfunction between the oscillatory region below R t and the surface. The value of E can be much larger for g-dominated modes, with a considerable amplitude in the gravity-wave propagating region in the deep interior, than for p-dominated modes. This is conveniently characterized by the scaled inertia whereĒ 0 (ω) is the radial (purely acoustic) mode inertia at the frequency of the mode considered.
For the analysis of red-giant oscillations it is convenient to consider the fraction of the mode inertia that comes from the inner parts of the star relative to the total mode inertia. Following Goupil et al (2013) where E core is defined as in Eq. (38) but restricting the integral to the region where ω < N, S l . Evidently ζ is small for modes trapped in the envelope, whereas ζ is close to one for modes trapped in the core. We also note that, as a rough approximation, assuming that the envelope contribution to the inertia is similar to the inertia of a radial mode with the same frequency. 10 For the practical evaluation and later analysis we note that the inertia can be expressed in terms of the displacement vector (19) as where R phot is the photospheric radius.

Mode order
The mixed nature of the modes precludes a simple identification of the radial order of a mode based on the number of nodes. Using an earlier analysis by Eckart (1960), Scuflaire (1974) and Osaki (1975) independently proposed a scheme plotting the eigenfunction in a suitable phase diagram, e.g., in terms of (ξ r , ξ h ); the radial order is determined by counting the zero-crossings of ξ r with a positive or negative sign depending on whether the curve crosses the axis in the counter-clockwise or clockwise direction in the phase diagram. These zero-crossings are associated with the regions in the star where the mode has a p-mode or a g-mode character, respectively. If the number of counter-clockwise crossings isn p and the number of clockwise crossings is |n g |, 11 withn g < 0, the mode order is n =n p +n g .
When the perturbation to the gravitational potential is neglected this defines a mode order that is not changed for a given mode as the star evolves, even though its dominant physical character may change. For modes of degree l ≥ 2 this property has also been found to be satisfied for solutions of the full equations of adiabatic oscillation. On the other hand, an application of the Eckart scheme to dipolar modes (l = 1) leads to a poorly defined mode order for evolved models (Lee, 1985;Guenther, 1991) or centrally condensed polytropes (Christensen-Dalsgaard and Mullan, 1994). It was shown by Takata (2005Takata ( , 2006) that a relation satisfied by the eigenfunctions of dipolar modes allows the definition of a scheme for the determination of mode order that is well-defined and invariant under evolution. As in Eq. (44) this is characterized by contributionsn p andn g from the p-mode and g-mode dominated parts of the star. We use this in the later discussion of dipolar modes. Together with the original scheme proposed by Scuflaire (1974) and Osaki (1975) this defines mode orders that are invariant under evolution. Thus the order is unchanged when following a given mode as the star evolves in Fig. 13. Formally modes with positive n may be classified as p modes and modes with negative n as g modes. This classification largely corresponds to the physical nature of the modes for unevolved stars. However, as discussed below the non-radial mixed modes develop a large number of nodes in the g-mode propagation region as the star evolves on the red-giant branch, and hence typically have a large negative value of n g . Thus for such stars all relevant modes have negative n, and other properties of the modes must be used to characterize their physical nature. In any case it should be kept in mind that the radial order defined here is a purely theoretical concept, although very useful in characterizing the oscillation modes of a given stellar model.

Subgiant stars
As a star evolves, the first g-dominated mixed modes start to be observable in subgiants. This provides a useful illustration of the properties of mixed modes. The characteristic frequencies for a subgiant model are shown in Fig. 14. This model has a helium core of radius 0.022R, containing 5 % of the star's mass. The behaviour of the acoustic frequencies differs little from those of a main-sequence model. However, the compact core and resulting high gravitational acceleration give rise to a high peak in the buoyancy frequency in the core, augmented by the sharp composition gradient in the hydrogen-burning shell. Therefore, at the indicated typical oscillation frequency (horizontal dotted line in Fig. 14) there are two trapping regions. The increase in the core buoyancy frequency with age leads to avoided crossings when following the evolution of the modes with age. This is illustrated in the top left panel of Fig. 15, which shows the behaviour of radial and dipolar modes. The resulting changes in the character of the modes in terms of the mode inertia for two of the dipolar modes, compared with a neighbouring radial mode, are shown in the lower left panel of Fig.  15. As long as the dipolar modes are predominantly acoustic their inertia is very similar to that of the radial mode. As a mode undergoes an avoided crossing and takes on a substantial g-mode character its inertia increases. At the next avoided crossing its inertia decreases again and the mode returns to an acoustic character, exchanging character with the next mode. At the point of closest approach the inertias of the two modes are very similar. Fig. 14 shows that for dipolar modes the evanescent region is relatively thin. This leads to strong coupling between the two oscillatory regions, a relatively large minimum separation during the avoided crossings, and a rather modest increase in the mode inertia when the modes are most g-mode like. This should be contrasted with the case of quadrupolar modes, shown in the right panels of Fig. 15. Here the evanescent region is considerably thicker, the coupling consequently weaker and the avoided crossings very sharp. Also, the inertia increases very rapidly as a mode takes on predominantly g-mode character. This shows that in any given model it is unlikely to find a mode that is not either predominantly p-or g-dominated.
The diagnostic potential of mixed modes in subgiant stars was discussed in Section 3.4.1.
The presence of mixed modes affects the asymptotic behaviour of the pure acoustic spectrum, such as that of the Sun (see theéchelle diagram in the bottom panel of Fig. 16). In case the density of g-dominated modes is relatively low (i.e., N γ ≪ N π ; see above), each avoided crossing effectively adds another mode to the frequency spectrum. The resulting changes in the oscillation spectrum were analysed by Deheuvels and Michel (2010) based on a simple physical model of coupled oscillators. They considered the coupling between a single γ mode and several π modes. 12 The effect on the distribution of peaks in theéchelle diagram is illustrated in Fig.  16 for the model marked by a vertical dotted line in Fig. 15. At high frequency the dipolar modes have not yet been affected by g-mode mixing, and the behaviour corresponds to the purely acoustic case. However, as shown in Fig. 15 there is a pair of mixed dipolar modes with nearly the same inertia at a frequency around 930 µHz. These are visible as a pair of modes on either side of the l = 1 ridge in Fig. 16. At lower frequency essentially all the dipolar modes show an effect of the avoided crossings, as also argued by Deheuvels and Michel (2010).
For l = 2 the right hand panels of Fig. 15 also indicate the presence of gdominated modes. However, for these the scaled inertias Q n l are so high that they would not be visible with the proposed scaling of symbol size in Fig. 16. Therefore, they are indicated by plusses. As discussed above the high scaled inertias of mixed l = 2 modes are a consequence of the much weaker coupling between the buoyancy and acoustic cavities for quadrupole modes. Only mixed quadrupole modes with frequencies close to the pure acoustic l = 2 mode can reach observable amplitudes and could thus be detected. However, such modes would be difficult to identify observationally, except if a clear pair of closely spaced modes is detected.

Ascending-branch red giants
As a star evolves up the red-giant branch its internal structure and hence its oscillation spectrum changes dramatically. As an example of an ascending-branch star we consider the most evolved stellar model in Fig. 13, which is indicated with the upper diamond in Fig. 11. The properties of the modes of this model in terms of the mode inertia and scaled mode inertia (cf. Eqs 38 and 40), are shown in Fig. 17. The l = 1 and 2 modes form a very dense spectrum with inertias typically exceeding the inertia of the neighbouring acoustic modes by several orders of magnitude. These high-inertia modes are buoyancy-dominated modes that are predominantly trapped in the g-mode cavity in the deep interior of the star. Additionally, there are acoustic resonances where one or more modes have inertias close to the radial-mode inertia. This is particularly visible in the bottom panel of Fig. 17, where low scaled inertia Q n l indicates the modes that have largest amplitude in the p-mode cavity.  Fig. 15, of age 4.293 Gyr and with T eff = 5887.8 K. Modes with l = 0 are shown as circles, l = 1 as triangles, l = 2 as squares and l = 3 as diamonds. The size of the symbols is proportional to Q −1/2 (cf. Eq. 40), providing a rough estimate of mode amplitude relative to a radial mode of the same frequency (cf. Eq. 117). For four quadrupolar modes very small symbols have been replaced by '+'. The two dipolar modes at the centre of an avoided crossing in Fig. 15 are shown with red triangles. Bottom: for comparison we also show anéchelle diagram for the observed low-degree observations of the Sun (∆ν = 135.4 µHz; Chaplin et al, 2002).
To further illustrate the properties of the modes, Fig. 18 shows the integrands of the inertia for the most p-dominated and the most g-dominated dipolar modes with a frequency near 80 µHz. There is clear similarity in shape between the two curves; 13 the main difference is in the behaviour in the evanescent region near r/R = 0.1, where the eigenfunction increases with depth for the g-dominated mode (black curve) and decreases with depth for the p-dominated mode (red curve). As a result, the region beneath 0.1R contributes 99.5 % to the inertia of the g-dominated mode, while for the most p-dominated mode 63 % of the inertia comes from the g-mode cavity. Hence, even for p-dominated modes the contribution from the g-mode cavity is significant. Additionally, Fig. 18 illustrates the extremely rapid variation of the eigenfunction that results from the high buoyancy frequency in the g-mode cavity. This rapid variation becomes even more extreme in more evolved red giants where the eigenfunctions may have thousands of nodes in the deep interior. This places severe requirements on the  numerical techniques used to compute these modes. Techniques that explicitly take such rapid variations into account have indeed been developed (Gabriel and Noels, 1976;Townsend and Teitler, 2013).

Pressure-dominated modes
The acoustic resonances of non-radial modes together with the radial modes satisfy a frequency pattern similar to the asymptotic behaviour of acoustic modes in main-sequence stars (cf. Eq. 30). This is evident in theéchelle diagram shown in Fig. 19 where in particular the l = 0 and 2 modes have a behaviour very similar to what is seen in the solar case (cf. bottom panel of Fig. 16). In addition to these p-dominated modes a few additional l = 2 modes are visible with the chosen scaling. For l = 1 the pattern is more complicated, i.e., several mixed modes are visible per acoustic-mode order. Nevertheless, a dominant set of modes still follows the asymptotic expression. This regular pattern is a universal feature of observed (and modelled) red-giant oscillations (Mosser et al, 2011b, and Section 3.3.2). However, the reason that the asymptotic relation holds is not entirely clear. The analysis leading to Eq. (30) is fundamentally related to the behaviour of the solution of the oscillation equations near the singularity at r = 0 (see also Gough, 1986aGough, , 1993. However, in the case of red giants the properties of the acoustic resonances are determined by the requirement that the solution decrease with increasing depth in the evanescent region (see also Fig.  18). Here conditions are very different from the conditions in the central regions of a main-sequence star. Despite these different conditions the universal pattern seems to hold. A better physical understanding of this behaviour of the frequency pattern of the acoustically dominated modes may lead to additional diagnostic potential of the observations.

Gravity-dominated modes
The properties of gravity-dominated modes are to a large extent determined by the structure of the core. Hence, these oscillations have a large potential for the study of the internal structure of red giants. It was found by Beck et al (2011) that the observed spectrum of a red giant showed additional peaks which were identified as coming from mixed modes with a stronger g-dominated component (see Fig. 1). The behaviour of these mixed modes is dominated by the g-mode asymptotic relation, Eq. (36). According to this relation the gravity modes are uniformly spaced Fig. 21 Frequency as a function of scaled mode inertia Q (left), period spacing (centre) and period modulo period spacing (∆Π 1 = 75.4 µHz), i.e., periodéchelle diagram (right) for l = 1 modes in a 1.3 M ⊙ stellar model with a radius of 6.2 R ⊙ (see Jiang and Christensen-Dalsgaard, 2014). The red diamonds mark the modes illustrated in Fig. 18 or, in the central panel, the period spacings between those marked modes and the adjacent mode with higher period. in period. This is illustrated by Fig. 20 showing the period spacing ∆Π n l = Π n l −Π n+1 l for l = 1, as a function of frequency. At relatively low frequency there are modes with nearly constant ∆Π n l while the spacing shows a characteristic decrease around the more p-dominated modes. The actual period spacing is in excellent agreement, particularly at low frequency, with the asymptotic period spacing marked by the horizontal dashed line. Thus observations of low-frequency g-dominated modes provide a direct measure of the integral of the buoyancy frequency in the core of the star.
The fact that several dipolar mixed modes are visible in the vicinity of the acoustically dominated modes, as observed by Beck et al (2011), might indeed have been expected from the rough amplitude scaling in anéchelle diagram such as shown in Fig. 19. In fact, as discussed in Section 4.3, even fully g-dominated modes may have sufficient visibility to be observed in very long time series, such as those that were obtained by Kepler's nominal 4-year mission. An example is shown in fig. 1b of Stello et al (2013a) for a red-clump star.
Following Bedding et al (2011) we introduced periodéchelle diagrams in Section 3.3.7, based on the uniform period spacing of g-dominated modes. Figure 21 shows such a diagram for a part of the modes shown in Fig. 20. For comparison we also show the scaled mode inertia and the period spacing. Gravity-dominated modes with high inertia fall approximately on a vertical line in the periodéchelle diagram. This is as expected from the asymptotic behaviour. Around this vertical line departures are visible near the acoustic resonances. It may be noticed that the variations near the resonances are qualitatively similar to the variations induced by g-dominated modes in the frequencyéchelle diagram for a subgiant star (cf. Fig. 16).
Coupling of buoyancy and acoustic cavities Asymptotic analysis of the mixed modes in red giants must take into account the coupling between the buoyancy and acoustic cavities. A full analysis of this problem has been carried out by Shibahashi (1979).
(see also Unno et al, 1989). 14 The analysis is based on simplified equations of the form given in Eq. (21), and has been implemented for the red-giant case by Goupil (private communication) and Mosser et al (2012c). Shibahashi showed that continuous matching of the solutions between the two cavities leads to the resonance condition implicitly determining the eigenfrequencies. In Eq. (45) r a , r b , r c and r d ≃ R are the four turning points where K changes sign (cf. Fig. 12), and φ g and φ p are phases depending on the properties of the turning points. Additionally, is a measure of the strength of the coupling between the two cavities. To analyse Eq. (45) we approximate K by K ≃ l(l + 1)r −2 N 2 /ω 2 in the buoyancy cavity, [r a , r b ] and by K ≃ ω 2 /c 2 in the acoustic cavity, [r c , r d ]. Then Eq. (45) can be written as where (cf. Eq. 31) and (cf. Eq. 37). In Eq. (48) we neglected the small contribution to ∆ν from the core region.
To interpret Eq. (47) we first consider the uncoupled case, i.e. q = 0. This should in principle yield Eq. (30) and Eq. (36) for the uncoupled π and γ modes, respectively. In the uncoupled case Eq. (47) has one set of solutions with tan(ω/ω p − φ p ) = 0, or equivalently ν = ω 2π ≃ ∆ν(n p + ǫ p ) for integer n p , and ǫ p = φ p /π. This indeed superficially recovers the leading-order term in the acoustic-mode asymptotic relation (30) for the uncoupled π modes, although with φ p and ǫ p being independent of l. In contrast, the observed and computed p-dominated frequencies satisfy Eq. (30) which includes an l dependence. However, the analysis leading to Eq. (50) neglected the dependence on l of K in the acoustic cavity and the detailed behaviour near the lower turning point r c . A full understanding 14 Shibahashi used a slightly different form of the asymptotic equation, replacing K in Eq. (22) by This essentially only differs from Eq. (22) in the near-surface layers and hence affects, e.g., the phase ǫ in Eq. (30). The following analysis follows Shibahashi (1979).

Fig. 22
The effective phase Φ (cf. Eq. 60) as a function of frequency for l = 1, in a 1.3 M ⊙ model of radius 6.2 R ⊙ (see Jiang and Christensen-Dalsgaard, 2014). The horizontal dotted lines indicate multiples of π.
of the universal pattern (Eq. 6) requires a more complete analysis extending the one leading to Eq. (45), which has so far not been carried out. For now we simply replace ǫ p by ǫ p l (and similarly for φ p ) with ǫ p l = ǫ p 0 + l/2, such that we recover Eq. (30). The second set of solutions for an uncoupled (q = 0) case satisfy cot(ω g /ω−φ g ) = 0 or equivalently where n g is an integer and ǫ g = φ g /π. In accordance with the discussion in Section 4.2.1 we have chosen n g < 0, using n g rather thann g in this asymptotic expression. We discuss the relation between n g andn g below. In this case we indeed recover Eq. (36), defining the uncoupled γ modes. Now examining the coupled case, Eq. (47), it is convenient to introduce and Here we introduced the uncoupled p-and g-mode frequencies as per Eqs (50) and (51): ν (p) and ν (g) Additionally, in the last equalities in Eqs (52) and (53) we ignored the contributions n p π and |n g |π which have no effect in Eq. (47). It was noted by Deheuvels et al (2015) and Mosser et al (2015) that expressing θ p in terms of ν (p) n p l allows the use of a more complete description of the location of the acoustic resonances than provided by Eq. (54), such as the universal acoustic-mode pattern, Eq. (6). We use this convention in the following, defining θ p by the last equality in Eq. (52).
Mosser et al (2012c) rewrote Eq. (47) by moving the g-mode related part to the right-hand side as tan(ω g /ω − φ g ) and taking the inverse tangent, to obtain This provides an implicit equation for ν, which can be determined iteratively, given the other parameters of the equation. Hence it provides the basis for determining these parameters through fits to observed frequencies. As discussed in Section 5.4 this results in powerful diagnostics for stellar structure and evolution. The detailed relation between the numerical order n and its componentsn p and n g (cf. Section 4.2.1) and the asymptotic indices n p and n g has apparently not been fully explored. Passing through an acoustic resonance ν = ν (p) n p l is associated with the increase by one in the number of nodes in the p-mode region of the star, and hence with an increase by one inn p . In fact it seems that we can identifyn p with n p for the acoustically-dominated modes. Similarly, there is a close relation betweenn g and n g characterizing the asymptotic properties of the g-dominated modes, at least for stars on the ascending red-giant branch. The situation is probably more complicated in more complex stars, possibly with several separate regions of g-mode propagation and/or buoyancy glitches (see also Sections 4.2.4 and 4.2.5).
As noted by Mosser et al (2012c) a useful quantity is the estimate of the number of g-dominated modes associated with a given acoustic resonance, which can be obtained as Adding the mode at the acoustic resonance the total number of modes in a ∆ν interval around the resonance is N + 1. An alternative way to analyse Eq. (47) was proposed by Christensen-Dalsgaard (2012a) 15 and developed further by Jiang and Christensen-Dalsgaard (2014). They rewrote Eq. (47) as sin θ p cos θ g − q cos θ p sin θ g = 0 .
This can be rewritten as where Φ satisfies By choosing the appropriate integer multiples of π we can define Φ as a continuous function of ω. Thus the eigenfrequencies satisfy for integer k, or equivalently The function Φ is illustrated in Fig. 22. To understand its properties it is convenient to rewrite Eq. (60) as Φ = arctan[q cot π(ν − ν (p) n p l )/∆ν]. 16 . The value of q is typically relatively small. Therefore, for most frequencies Φ is close to a multiple of π. This is illustrated by the dotted lines in Fig. 22. However, near the acoustic resonances, where ν ≃ ν (p) n p l , the argument to arctan goes through a singularity. This leads to the rapid variation of Φ.
We note that Eq. (62) provides a straightforward way to determine the frequencies of mixed modes, given ∆ν, ν (p) n p l , q, ∆Π l and ǫ g : Φ can be directly computed from Eq. (60), defining Ψ as a function of frequency from Eq. (62); the frequencies of the mixed modes are then simply obtained as those values where Ψ takes on consecutive integer values. In this way one avoids the iterative solution required to determine the frequencies from Eq. (56).
In Eq. (62) increasing k leads to decreasing frequency, corresponding to decreasing n, where n is the numerical mode order introduced in Section 4.2.1; also, Φ is only defined from Eq. (60) to within an integer multiple of π. It appears that this can be chosen such that k = −n. This reflects the fact that Eq. (62) defines the complete spectrum of mixed modes for the given degree.
The behaviour of Φ near the resonances is directly related to the variation in the period spacing (cf. Fig. 20). Writing Eq. (62) as with ν n l = 1/Π n l , and making a Taylor expansion of Φ(2πν k+1 ) − Φ(2πν k ) we obtain Thus for the g-dominated modes, where Φ is nearly constant, we recover ∆Π l . At the same time the period spacing decreases near the acoustic resonances (note that dΦ/dω < 0). Jiang and Christensen-Dalsgaard (2014) presented an approximation to Eq. (64) and discussed the use of this to estimate q from the variation in the period spacings.
The quality of the asymptotic fit derived above (Eq. 62) is illustrated in Fig. 23 for model M 4 of Jiang and Christensen-Dalsgaard (2014). For the g-dominated modes the relative differences between the computed and asymptotic frequencies are less than 10 −4 . The relative difference increases to 5 × 10 −4 for frequencies undergoing acoustic resonances. However, these differences have a far stronger effect on the period spacings.
The coupling constant q depends on the properties of the evanescent region and hence in principle provides further diagnostics of the stellar interior. This is reflected by Eq. (46) that was already obtained by Shibahashi (1979). The values of q determined from Eq. (46) and from fits to numerical frequencies were compared by Jiang and Christensen-Dalsgaard (2014), who found substantial differences between these values. This appears in part to be due to the neglect of the perturbation to the gravitational potential in the asymptotic analysis. Furthermore, the approximation to Eq. (64) by Jiang and Christensen-Dalsgaard (2014) neglects potential singularities, particularly in the evanescent region. These effects require further analysis.
The phase ǫ g in principle provides further information about stellar structure, particularly the regions at the turning points of the g-mode cavity. Buysschaert et al (2016) carried out fits to data for three red giants obtained by Kepler, including a variable ǫ g as one of the fitted parameters (see Section 3.3.7). They found that ǫ g was only weakly constrained by the observations; however, allowing ǫ g to vary resulted in a clearer definition of the values of the period spacing and its uncertainty obtained in the fit. These issues also deserve further studies.

Asymptotic properties of mode inertia
It is instructive to consider the asymptotic properties of the mode inertia, in particular the relative contribution ζ from the core (cf. Eq. 41). The value of ζ is closely related to the effects of rotation on the oscillation frequencies and the damping of the modes. The following is based on the detailed analysis by Shibahashi (1979, S79), which was further developed by Goupil et al (2013) and Deheuvels et al (2015).

S79 based the asymptotic analysis on the two functions v(r) and w(r) related to ξ r and ξ h by
In the outer parts of the star the modes are predominantly acoustic, with |ξ r | ≫ |ξ h |.
In this case the integral for the inertia (Eq. 43) is dominated by the term in ξ r . From the asymptotic expression for v (eq. 28 of S79, applied in the acoustic outer region) we here obtain where C is a normalization constant and φ ′ p is a phase. As noted above, we use the definition of K given in footnote 14. As in the derivation of Eq. (48), we neglected S 2 l /ω 2 in the amplitude function in Eq. (65) compared with 1, which is a good approximation except close to the turning point r c . Using the same approximation, we obtain from Eq. (66) the contribution of the outer parts of the star to the numerator in Eq. (43) is the acoustic depth, and τ c is the acoustic depth of the turning point r c . For a mode of high acoustic order the integral over cos 2 can be replaced by 1/2 τ c . Using Eq. (48) we obtain In the inner parts of the star the mode is buoyancy dominated, with |ξ h | ≫ |ξ r |. Here, the inertia integral is dominated by the term in ξ h . In this case the relevant asymptotic solution, from the relevant expression for w (eq. 29 of S79, in the gmode cavity) is where L 2 = l(l + 1), A is a normalization constant and φ ′ g is a phase. As in the derivation of Eq. (49) we neglected 1 compared with N 2 /ω 2 . From Eq. (70) we then obtain for the contribution of the g-mode region to the numerator in Eq. (43) where we extended the integrals to r = 0. It should be noticed that the amplitude of the integrand scales like N, when integrating with respect to ln r. This is indeed confirmed by comparing the integrand shown in Fig. 18 with the buoyancy frequency in Fig. 12. Introducing the buoyancy radius we obtain where υ b is the buoyancy radius of the turning point r b , and replacing as before the integral over sin 2 by 1/2 υ b . Using Eq. (49) yields Finally, by combining Eqs (69) and (74), we obtain To complete the analysis we need the ratio C/A between the amplitudes of the eigenfunctions in Eqs (66) and (70). This ratio is determined by the coupling across the evanescent region, [r b , r c ]. From eq. (30) of S79, and using the dispersion relation (Eq. 47), it can be shown that Thus or, using Eqs (52) and (53), Thus the fractional contribution of the g-mode region to the inertia is obtained as  (2014). The heavy dashed line shows the asymptotic period spacing ∆Π 1 = Π 0 / √ 2, and the red dotted line and symbols show ζ∆Π 1 (cf. Eq. 82), where ζ was obtained from the numerically computed eigenfunctions. . They also demonstrated that ζ as provides a remarkably good fit to the numerically computed ζ. By using Eq. (47) this can be further simplified as Mosser et al (2015) pointed out that Eq. (80) shows a very interesting relation between ζ and the actual period spacing between mixed modes. To see this, we return to Eq. (64). Differentiating Eq. (60) we obtain dΦ dω = −qω −1 p 1 (q 2 cot 2 θ p + 1) sin 2 θ p = −qω −1 Substituting this into Eq. (64), using Eqs (48) and (52), yields with ζ as given by Eq. (80). This is the result obtained by Mosser et al (2015). Equation (82) is satisfied to high accuracy by computed values of ζ and period spacing for red-giant models. As an example, Fig. 24 shows part of the period spacing from Fig.  20, but including also ζ∆Π 1 , with ζ computed from the numerical eigenfunctions. We note that the derivation of ζ as , and hence the relation in Eq. (82), assumes the validity of the comparatively simple asymptotics employed by Shibahashi (1979), which requires a slow variation of the equilibrium quantities. This approximation is not valid in the case of sharp features in the buoyancy frequency, such as is the case for Model 1a of Cunha et al (2015) (see also Fig. 30). In these cases the more complex analysis of Cunha et al would be required.
From Eq. (82) Mosser et al (2015) and Vrard et al (2016) noted that the properties of ζ as (or ζ) 17 are largely controlled by the properties of the acoustic resonances, with little dependence on the period spacing. As a result they were able to develop a technique for determining ζ from the observations in the vicinity of the acoustic resonances. Given ζ, they introduced a stretching function P(ν) 18 by and the stretched periodsΠ nl = P(ν nl ). Replacing ζ by ζ as and relating ζ −1 as to dΦ/dω using Eqs (80) and (81) yield or, with suitable choice of integration constant, using the definition of Ψ in Eq. (62). From Eq. (62) it follows that the stretched periods of the modes satisfyΠ n l = ∆Π l (k + ǫ g + 1/2) .
Thus we formally recover the relation (51) for the uncoupled g modes. However, this relations is now valid for all modes. As shown by Mosser et al (2015) this greatly simplifies the periodéchelle diagram when expressed in terms of the stretched period. This is illustrated in Fig. 25; as suggested by Eq. (86) all modes fall essentially on a straight line. Vrard et al (2016) developed the stretching technique into an automated method for analysing the red-giant mixed-mode spectra of large samples of stars; this is particularly powerful in dealing with the complications of rotational splitting. We return to this in Section 4.4.

Effects of glitches
Much of the preceding discussion of the properties of red-giant oscillations was based on the asymptotic properties of the modes. This conveniently allows characterizing the oscillation spectrum by a relatively small number of parameters, such as ∆ν, ∆Π 1 , etc. However, the asymptotic approximation also indicates a limitation of the diagnostic potential of the observations, as far as the detailed internal properties of the stars are concerned. It should be kept in mind that the asymptotic analysis is based on the assumption that the underlying stellar structure varies on a scale that is long compared with the local wavelength of the oscillations. Structure features on a smaller scale introduce perturbations to the frequencies which may provide diagnostics about these structures. This was pointed out in the case of p modes in the Sun by Gough (1990), who noted that at the base of the convection zone in normal solar models the second derivative of the sound speed is discontinuous. In addition, the relatively rapid variation of Γ 1 , and hence the sound speed, in the second helium ionization zone may also cause departures from the asymptotic behaviour. Vorontsov et al (1991) used the signature of helium ionization in observed solar frequencies to infer the solar envelope helium abundance. Such features in the sound speed were denoted acoustic glitches by Gough (2002). Similarly, sharp features in the buoyancy frequency (buoyancy glitches, Cunha et al, 2015) may affect g-dominated modes.
The base of the convective envelope in red giants is typically too deep to induce a significant effect on the p-dominated modes ). An example of the acoustic glitch associated with the helium ionization zone in a red-giant star is illustrated in Fig. 26, showing the derivative of sound speed in terms of acoustic depth τ (cf. Eq. 68). Here the dominant feature is the dip in dc/dτ near τ = 1.5×10 4 s which is associated with the second helium ionization zone. For comparison the lower panel shows the scaled displacement eigenfunction of a radial mode of frequency 98.1 µHz,   Fig. 26 for l = 0 (circles) and the most acoustically dominated modes with l = 1 (triangles) and l = 2 (squares).
with a variation at the helium feature which is substantially slower than the structural variation. Thus the helium dip does act as an acoustic glitch. Gough (1990) showed that the effect of an acoustic glitch is to introduce a variation δν (gl) of the form δν (gl) ∝ sin(4πντ gl + 2φ gl ) (see also Eq. 18), where τ gl is the acoustic depth of the glitch and φ gl is a phase. Thus at fixed τ gl we expect an oscillatory behaviour of δν (gl) as a function of cyclic frequency, with a period of around (2τ gl ) −1 . Gough (1990) furthermore proposed isolating the effect of the glitch from other possible slower variations of the frequencies with radial order by considering second differences at fixed l, This is illustrated in Fig. 27 for the same model as shown in Fig. 26. The effect of the helium glitch is only relevant for the acoustic properties of the modes, and hence we have selected just those modes with l = 1 and 2 with the lowest inertia in each radial-mode frequency interval. The figure shows the expected oscillatory behaviour, although modes with l = 1 that are more strongly affected by the mixed nature show some scatter around the general trend. The period of the oscillation is roughly consistent with the acoustic depth of the glitch. The amplitude of the oscillation depends on the strength of the glitch, which is determined by the magnitude of the variation in Γ 1 and hence provides a diagnostic of the helium abundance in the convective envelope of the star. With increasing frequency the local wavelength of the eigenfunction gets shorter, making the glitch appear less sharp and hence reducing the amplitude of the variation in the second difference (see also Section 3.4.5).
Houdek and Gough (2007) carried out a detailed analysis of the effects of the acoustic glitches in the first and second helium ionization zones and the base of the convective envelope in the Sun. This led to an expression already summarized in Eq. (18). They noted the importance of a proper definition of the acoustic surface, i.e., the zero point of the acoustic depth, in the interpretation of the analysis of acoustic glitches. The potential uncertainty involved in the definition of the acoustic surface, and the effects of the near-surface errors in stellar modelling (cf. Section 3.4.3), can in principle be avoided by converting the acoustic depth to an acoustic distance from the centre (e.g., Ballot et al, 2004); this uses the fact that the asymptotic large frequency separation ∆ν is related to the total acoustic radius of the star (cf. Eq. 31). In practice the acoustic depth of a feature and the large frequency separation may respond differently to the near-surface uncertainties, in particular since ∆ν is typically determined from a suitable fit or average over the observed modes.
Buoyancy glitches, i.e., sharp features in the buoyancy frequency, are of great potential interest owing to their sensitivity to the composition structure. Analysis of buoyancy glitches caused by the composition discontinuities established through gravitational settling has played a major role in the study of white-dwarf oscillations as important diagnostics of the internal structure of the stars (e.g., Winget et al, 1994). In a detailed analysis Miglio et al (2008) demonstrated the diagnostic potential of the buoyancy glitch at the edge of convective cores in more massive main-sequence pulsators, leading to oscillations with mode order in the period spacings of high-order  g modes. Evidence for this behaviour was found by Degroote et al (2010) in a slowly pulsating B star observed by CoRoT, indicating effects of processes smoothing the composition gradient at the edge of the core. Thus it is of great interest to consider the potential for studies of buoyancy glitches in evolved stars. Figure 12 shows the presence of potential glitches in the buoyancy frequency in red-giant models. These were analysed in considerable detail by Cunha et al (2015), including an asymptotic analysis of the effect of a discontinuity in composition and a resulting singularity in the buoyancy frequency (cf. Eq. 34). Here we discuss aspects of these issues using as an example a model corresponding to Model 1a of Cunha et al (2015), but with a density discontinuity at the base of the dredge-up region. The buoyancy frequency for this model is shown in Fig. 28. The model is in a somewhat later evolutionary stage, just before the red-giant 'bump', than the model shown in Fig. 12. Here the convective envelope has retracted in mass since the maximal extent of the first dredge-up (see Section 2.1.2), leaving behind a composition and density discontinuity within the g-mode propagation region, here shown as a sharp spike in N at r/R ≃ 0.02. This is obviously a glitch. 19 The second relatively sharp feature is associated with the hydrogen-burning shell. As illustrated in Fig. 29 this is not a glitch, however, since the eigenfunction varies extremely rapidly in this region. This is an immediate consequence of the asymptotic behaviour in Eq. (70), according to which the local radial wave number scales as N/r and hence becomes very large at the peak in the buoyancy frequency.
The analysis by Cunha et al (2015) shows that, for a pure g mode, the presence of a buoyancy glitch gives rise to an oscillatory behaviour such as was found by, e.g., Miglio et al (2008), with a period depending on the location of the glitch within the g-mode cavity. Taking into account also the mixed nature of the modes, with the coupling to the p-mode behaviour, leads to a complex variation in the frequencies and the period spacings, as illustrated in Fig. 30. 20 The presence of a glitch clearly invalidates the derivation of the asymptotic relation in Eq. (82) between ζ and the period spacing, as also illustrated in Fig. 30. This shows a very interesting diagnostic potential for the characterization of the chemical profile left by the dredge-up in red giants just below the bump; however, since the effects are subtle a more detailed analysis is required to ascertain whether the effect can reliably be isolated in the existing data. Cunha et al (2015) proposed a diagnostic of the irregularity of the period spacings, based on a Fourier transform of period spacing as a function of period, which provides a strong indication in the model results of glitch effects. This also demonstrated that the effects in red-giant models are closely linked to the red-giant bump.
As discussed in the next section buoyancy glitches play a dominant role in the oscillation spectrum of core helium-burning stars. 19 As discussed by Cunha et al (2015) this feature remains sharp compared with the local wavelength even if some numerical (or physical) diffusion is allowed, smoothing the composition gradient. This deserves further study. 20 We note that this is rather more complex than the simple model of the effects of a buoyancy glitch considered by Mosser et al (2015).

Fig. 30
Period spacings for modes with l = 1, in a 1 M ⊙ model at the bump, with buoyancy glitch (adapted from Cunha et al, 2015). The red dotted line shows ζ∆Π 1,as . The lower panel provides a blow-up of the region near the asymptotic period spacing.

Stars on the red clump
A major early realization in the analysis of space-based asteroseismic observations of evolved stars was that there is a sharp distinction between stars on the ascending red-giant branch, with hydrogen fusion around the helium core, and stars in the red clump where in addition there is helium fusion near the centre Mosser et al, 2011a) (see also Section 3.3.7). The g-mode period spacing is substantially smaller in the red giants than in clump stars. This is a straightforward consequence of the difference in the internal structure between these different stars (at fixed surface radius, say), given the asymptotic expression for the period spacing, Fig. 31 Helium abundance as a function of fractional radius in a model in the core helium-burning phase, with initial mass 1.5 M ⊙ , radius 12.2 R ⊙ and effective temperature 4648 K.
Eq. (37), and the expression for the buoyancy frequency, Eq. (34). The core heliumburning stars have a convective core (see Figs 4 and 5), which restricts the g-mode cavity and contributes to reducing the integral over the buoyancy frequency in Eq. (37), hence increasing the asymptotic period spacing (e.g., Christensen-Dalsgaard, 2014). In addition, with the ignition of helium burning the core expands, reducing the local gravitational acceleration and hence the buoyancy frequency (cf. Eq. 34), further increasing the period spacing.
To illustrate these effects we consider models discussed by Christensen-Dalsgaard et al (2014). The evolution sequence, extending from the pre-main sequence through the central helium-burning phase, was computed with the GARSTEC code (Weiss and Schlattl, 2008). This includes full treatment of the helium flash and the subsequent sub-flashes (see Section 2.1.3). The resulting helium-abundance profile Y in a model with initial mass 1.5 M ⊙ near the end of central helium burning is illustrated in Fig. 31. The region of constant Y inside 0.0015 R is convective. As discussed in Section 2.1.4 a discontinuity in Y, and hence in density, may be established at the edge of the convective core. In the rest of the helium core the smaller steps in composition reflect the convective mixing associated with the initial helium flash and the sub-flashes.
The buoyancy frequency of this model is illustrated in Fig. 32 and compared with that of a red-giant model of the same initial mass and approximately the same radius and slightly lower effective temperature. This immediately shows the two effects mentioned above: the presence of a convective core and the lower general level of the buoyancy frequency in the core helium-burning model. As a result the asymptotic dipolar period spacing in the core helium-burning model is ∆Π 1 = 244.2 s while in the red-giant model ∆Π 1 = 56.9 s. This is the effect found by Bedding et al (2011) and Mosser et al (2011a) and extensively used since then to characterize evolved stars (e.g., Stello et al, 2013b;Mosser et al, 2014, see Fig. 47).

Fig. 32
The solid curve shows the buoyancy frequency in the core helium-burning model illustrated in Fig. 31. For comparison, the dashed curve shows the buoyancy frequency in a red-giant model on the same 1.5 M ⊙ initial mass evolution sequence, with a radius of 12.0 R ⊙ and effective temperature 4532 K. The solid and dashed red lines show the Lamb frequency S 1 for the clump and red-giant models, respectively. The dotted line indicates the frequency ν max of maximum oscillation power (cf. Eq. 26). Fig. 32 is the extent of the evanescent region between the buoyancy and acoustic cavities, which is much smaller in the clump than in the red-giant model, leading to a stronger coupling. The effect of this is illustrated in frequencyéchelle diagrams in Fig. 33 where, as in Fig. 16, the symbol size provides an indication of the expected amplitude of the modes. It is evident that the stronger coupling in the clump model leads to a much broader spread of modes, particularly for l = 1. This is probably the reason for the less regular behaviour of the dipolar mode separations found by Montalbán et al (2010) who identified the most acoustic modes as those modes having the smallest inertia in a given radial-mode interval (see also Section 3.3.6). To illustrate this, these modes are shown as red in Fig. 33.

A second clear difference between the red-giant and clump models in
The diagnostics of clump stars has so far predominantly been based on the observed and asymptotically calculated period spacings. However, the actual behaviour of the model frequencies is rather more complicated, owing to the composition discontinuities and resulting spikes in the buoyancy frequency, also visible in Fig. 32. These act as glitches and lead to partial reflection of the g modes. The specific behaviour depends strongly on fine details of the model calculation, such as whether diffusion (numerical or physical) is included and whether the composition profile and buoyancy frequency are properly resolved in the numerical computations. Figure 34 shows the computed dipolar period spacings in the core helium-burning model in Fig. 32, compared with the asymptotic period spacing and the period spacing expected from the asymptotics in Eq. (82), but using the numerically computed ζ. The computed period spacing is clearly related to the expected behaviour in this case, For clarity the frequencies have been shifted by ν 0 = 0.4 µHz in both panels. Circles, triangles and squares show modes of degree 0, 1 and 2, respectively. As in Fig. 16 the symbol size for the non-radial modes have been scaled with Q −1/2 , where Q is the inertia ratio (Eq. 40), providing a rough indication of the expected amplitude of the modes. The dipolar modes with lowest inertia in each radial-mode frequency interval are shown as red in the bottom panel.
reflecting the location of the acoustic resonances shown by ζ, but with substantial variations caused by the glitches in the buoyancy frequency. 21 To illustrate further the sensitivity of the clump-model oscillations to details of the model structure we consider two models with initial mass 1 M ⊙ computed with the LPCODE (Althaus et al, 2005). Their helium-abundance profiles Y are illustrated Fig. 34 Dipolar period spacings in the core helium-burning model illustrated in Fig. 31. The horizontal dashed line shows the asymptotic period spacing while the red dotted line shows ζ∆Π 1 , using ζ computed from the numerical eigenfunctions.

Fig. 35
Helium abundance as a function of fractional radius in two models in the core helium-burning phase, with initial mass 1.0 M ⊙ , radius 12 R ⊙ and effective temperature 4699 K. The inset shows the region near the sharp decrease in Y at r/R ≃ 0.0047 caused by convection associated with the helium flash at the start of helium burning. One model (solid) has a composition discontinuity, whereas in the second model (dashed) this has been slightly smoothed by diffusion. (Models courtesy of Marcelo Miguel Miller Bertolami.)

Fig. 36
Top panel: Buoyancy frequency in the vicinity of the sharp decrease in the helium abundance in the models illustrated in Fig. 35. For the nearly discontinuous helium profile the buoyancy frequency is essentially a delta function (solid line). The buoyancy frequency for the slightly smoothed model is shown by a dashed line. Bottom panel: scaled horizontal displacement (on arbitrary scale) for a mode in the model with the nearly discontinuous helium profile, with l = 1, ν = 19.9 µHz. (Models courtesy of Marcelo Miguel Miller Bertolami.) in Fig. 35. The models only differ in the sharp decrease in Y caused by the off-centre onset of the helium flash and the associated convective mixing. In one model this is essentially discontinuous, whereas in the second model the composition has been slightly smoothed by diffusion. The resulting buoyancy frequencies in the vicinity of this glitch are shown in Fig. 36, which also shows the scaled horizontal displacement for a dipolar mode in the 'sharp' model. For the model with the sharp helium profile this variation in the buoyancy frequency is clearly a glitch. For the slightly smoothed model the scales of the buoyancy frequency and the eigenfunction are comparable, and hence one may expect a smaller effect of the local model structure. 22 This is indeed confirmed by the dipolar period spacings shown in Fig. 37. For the 'sharp' model the period spacings vary wildly, with little indication of the effects of the acoustic resonances, as reflected by the period spacing expected from asymptotics (cf. Eq. 82) shown by the dotted red curve. In the slightly smoothed model (bottom panel) the variations in the period spacing are less dramatic and there are some indications of the decreases in the period spacing expected from the asymptotic behaviour, shown by the red dotted curve. However, even in this case a proper analysis of the variation in the period spacing, and a determination of its asymptotic value, would be difficult.
The behaviour of the internal structure during the phases of the helium flash and sub-flashes was analysed by Bildsten et al (2012). They found strong changes in the asymptotic period spacing and the coupling between the buoyancy and acoustic propagation regions, noting that this may provide opportunities for the detailed diagnostics of this evolution phase. Even though it is brief, the large number of clump stars for which oscillation data are available from Kepler may make it realistic that at least a few stars are in this phase. Cunha et al (2015) pointed out that the expansion of the helium core at the onset of helium burning compresses the hydrogen-burning shell, so that it acts as a buoyancy glitch during this phase (unlike on the red-giant branch; cf. Fig. 29). This would add to the diagnostic potential of the observed oscillations. We note, however, that the early helium-burning model considered by Cunha et al (2015) did not have discontinuities in the composition and hence did not suffer from the strong glitch effects illustrated in Fig. 37. Constantino et al (2015) carried out an extensive analysis of the oscillation properties of stars in the clump phase, varying also convective-core overshoot and other mixing processes at the edge of the convective core. This clearly showed the very rich pulsational behaviour that may be found in these stars. An interesting issue raised was the relation to the g-mode oscillations of the core helium-burning subdwarf B stars which essentially correspond to the naked core of the clump stars (e.g., Østensen et al, 2014). Constantino et al also addressed the fact that computed asymptotic period spacings for clump models tend to be lower than the observed spacings (see also Fig. 47). A possible solution, although perhaps with limited physical justification, would be the effect of what Constantino et al call the 'maximum-overshoot' scheme. Similarly, Bossini et al (2015) considered various mixing schemes at the edge of the helium-burning convective core and their effects on the period spacing. They also pointed out that combining observed period spacings with the location of the so-called AGB bump in stellar clusters, resulting from non-monotonic luminosity evolution at the onset of helium shell burning (see Section 2.1.5), would provide further observational constraints on these mixing processes. Further studies are certainly needed on the apparent discrepancies between the theoretically predicted detailed properties and the observations, and of the consequences for the interpretation of the observations in terms of simple diagnostics such as the asymptotic period spacing.

Energetics of stellar oscillations
It is generally accepted that solar-like oscillations, including those in the Sun, are intrinsically damped (but see Xiong and Deng, 2007) and excited stochastically by the near-surface convection whose near-sonic speed makes the gas motions efficient in generating acoustic noise (Stein, 1968). This is confirmed by analysis of the statistical properties of the variations in solar-oscillation amplitudes (Chaplin et al, 1997), 23 which follow the pattern expected from stochastic excitation (Kumar et al, 1988;Chang and Gough, 1998). Also, there is a striking similarity of the distribution of amplitudes with frequency from the main sequence to red-giant branch, apart from the scaling of the frequency of maximum power with surface gravity (e.g. De Ridder et al, 2009;Stello et al, 2010). Furthermore, there is strong evidence that the variability seen in very evolved stars represents an extension of the red-giant solar-like oscillations (Dziembowski and Soszyński, 2010;Mosser et al, 2013a). This possibly includes semiregular variables, where analysis of up to century-long series of amateur observations has revealed statistical properties matching those of solar-like oscillations . A concise review of energetics of red-giant oscillations was provided by Dupret and Belkacem (2012).
The general theory of stochastic forcing of a damped oscillator was described by Batchelor (1953). It was applied to the analysis of observed solar modes by Christensen-Dalsgaard et al (1989). The resulting average power spectrum of a single mode, with angular frequency ω 0 and damping rate η = −ω i (cf. Eq. 20), has the form where P f is the average power spectrum of the forcing function, which varies relatively slowly with frequency. Consequently the spectrum is a Lorentzian with a full width at half maximum in angular frequency of 2η; in terms of cyclic frequency the full width at half maximum is Measurement of the width of the observed peaks therefore provides a measure of the damping rates of the modes (see also Section 3.4, in particular Eq. 14). In addition to the damping rate it is convenient to characterize the damping by the mode damping time t damp , defined by the time required to reduce the mode amplitude by a factor e and given by

Properties of the damping rate
To study the energetics of the modes the full set of non-adiabatic equations must be solved. These include also perturbations to the energy transport and the energy equation and result in a determination of the complex frequency ω = ω r + iω i as an eigenvalue, and hence the damping rate η = −ω i . However, to analyse η it is convenient to express it in terms of the work integral (e.g., Baker and Kippenhahn, 1962;Cox, 1967, see also Aerts et al 2010). Considering just perturbations to thermodynamic quantities the result is where ε is the rate of energy generation per unit mass, F is the flux of energy, Γ 3 − 1 = (∂ ln T/∂ ln ρ) ad and the star indicates the complex conjugate; also, δ denotes the Lagrangian perturbation, i.e., the perturbation following the motion (see also footnote 24). The integral in the numerator of Eq. (92) reflects the operation of a heat engine, with δρ defining compression, and δ(ρε − div F) defining heating. If the integral is positive, ω i is positive and the mode is excited. However, in the case discussed here of solar-like oscillations the integrated effect is negative and the mode is damped. Equation (92) is generally valid, including for fully non-adiabatic solutions. A major source of uncertainty in the calculation of the work integral is the treatment of the perturbations to the convective contribution to the flux, particularly near the surface. An additional complication in the near-surface region is the effect of turbulent pressure p t which makes a significant contribution to the total pressure in the outermost parts of the convection zone and hence also affects the pulsations. As a consequence, the full expression of the damping rate becomes where η gas is given by Eq. (92) and (e.g., Balmforth, 1992a;Houdek and Dupret, 2015), Im denoting the imaginary part. We return to these convective effects below.
In much of the star where the oscillations are essentially adiabatic, the contribution to the work integral can be estimated from the quasi-adiabatic approximation, calculating all terms from the adiabatic eigenfunctions. In red giants a potentially substantial contribution to the damping comes from the buoyancy-dominated region in the core, where the very high radial order of the g-mode behaviour causes a strong diffusive damping. In the outer parts of the star, on the other hand, the contributions to the damping are dominated by the near-surface layers, where the properties of the oscillations are independent of degree. To analyse this it is therefore convenient to separate the work integral into contributions from the core and the envelope. We introduce such that the full expression for the work integral can be written Here, as in the determination of ζ (Eq. 41), we have separated the star into the envelope and the core, the latter being defined as the region where ω < N, S l . In the envelope contribution we can replace D by the functionD 0 for radial modes, interpolated to the frequency of the mode considered. Also, the denominator is closely related to the mode inertia (Eq. 38); using this we replace the denominator in the first term by Q times the corresponding integral for radial modes, similarly interpolated, where Q is the scaled inertia (Eq. 40). Thus the first term becomes Q −1η 0 , whereη 0 is the interpolated damping rate for radial modes. In the second term we replace the denominator by an integral just over the core, using Eq. (41). Thus we finally obtain the full damping rate as (see also Grosjean et al, 2014). As discussed below,η 0 can be obtained from solving the full non-adiabatic equations for radial oscillations. However, we first consider the estimate of the second term from the asymptotic properties of the eigenfunction. In the core there is no contribution to D from the turbulent pressure. For simplicity we neglect the term in the energy generation in Eq. (95), and we assume that the flux is purely radiative. Also, in the core we can use the quasi-adiabatic approximation, so that the eigenfunctions are real, and no complex conjugate is needed, and we replace (Γ 3 − 1)δρ/ρ by δT/T . Thus we obtain, for the contribution from the core in Eq. (97), The flux is obtained as defining the conductivity K; here a is the radiation density constant,c is the speed of light and κ is the opacity. In δ(div F) the radial component of F dominates, since it involves the second derivative of the rapidly varying eigenfunction. Thus, we neglect the tangential component of the flux and obtain, using also the neglect of the nuclear term, where L = 4πr 2 F r is the luminosity, F r being the radial component of the flux. Here, using Eq. (99), or, expanding the last term 24 In the asymptotic analysis of this equation for the extreme g-mode behaviour in the core we follow Godart et al (2009) (see also Dziembowski, 1977). From the oscillation equations it may be shown that and hence, according to footnote 24, It follows that δp/p is small compared with terms involving derivatives of ξ r and, given the quasi-adiabatic approximation, the same is true of δρ/ρ and δT/T . From the equation of continuity, it then follows that dξ r dr where in addition we neglected 2ξ r /r compared with dξ r /dr. In Eq. (102) 2ξ r /r and δK/K = O(δp/p) can be neglected compared with dξ r /dr. Using also and Eqs (104) and (106) we obtain Using Eqs (100), (104), (107) and the equation of hydrostatic support we can finally write the numerator in Eq. (98) as To approximate the integrand in Eq. (109) we use the g-mode asymptotic eigenfunction in Eq. (70). From Eq. (106) we furthermore have, to leading order, where again we introduced L 2 = l(l + 1). Here we used that, to leading order, ξ r satisfies an equation of the form in Eq. (21); in the g-mode cavity K can be approximated by Differentiating Eq. (70), Eq. (110) also yields ξ r ≃ −Aρ −1/2 ω −1/2 r r −3/2 L 1/2 N −1/2 cos Thus the integrand in Eq. (109) becomes The integral in the denominator of Eq. (98) is evaluated essentially as in Eq. (71), yielding core ρ|δ δ δr| 2 dV ≃ 4πA 2 ω −3 L (114) Fig. 38 Mode damping time t damp = 1/η core corresponding to the asymptotic core damping rate (Eq. 115), evaluated at the frequency ν max of maximum oscillation power estimated as in Eq. (26). The results were obtained for a 1 M ⊙ evolution sequence and are shown against surface radius (lower abscissa) and luminosity (upper abscissa) in solar units, for l = 1 (solid) and l = 2 (dashed).
Thus we finally obtain, replacing cos 2 and sin 2 by their average values 1/2, As a star evolves up the red-giant branch, increasing the luminosity and the mass of the core, the core gravitational acceleration and buoyancy frequency increase. According to Eq. (115) these effects all contribute to increasing the core damping rate. This is illustrated in Fig. 38, showing the corresponding damping time in a 1 M ⊙ evolution sequence. On the low red-giant branch the core damping is small and unlikely to affect the observed properties of the modes, while around the bump, clearly reflected in the figure, the core damping rate has increased to a level where the corresponding damping time is comparable with or smaller than the duration of the nominal Kepler mission, such that significant effects can be expected (see also Fig. 41  below).
Returning to the near-surface contributions to the damping of the modes, a full non-adiabatic treatment of the oscillations is required, taking into account the pulsationinduced perturbations to the convective properties. These issues were discussed in detail by Houdek and Dupret (2015). Time-dependent generalizations of mixing-length theory were developed by Unno (1967) and Gough (1977bGough ( , based on earlier work in 1965, while Xiong (1977) used a Reynolds-stress model to treat convection in pulsating stars.
Unno's theory was further developed by Gabriel (1996) and Grigahcène et al (2005). With a suitable adjustment of parameters it provides a relatively reasonable fit to the detailed observations of solar oscillation linewidths (Dupret et al, 2006).

Fig. 39
Mode damping times t damp = 1/η in the 7 R ⊙ red-giant model from the 1 M ⊙ sequence illustrated in Fig. 38. The dotted curve shows radial-mode damping times computed using a nonlocal time-dependent convection formulation (Houdek et al, 1999), and the solid (l = 1) and dashed (l = 2) curves show the results of combining the radial-mode damping rates with the asymptotic estimate (Eq. 115), as in Eq. (97). The vertical dot-dashed line marks an estimate of the frequency ν max of maximum oscillation power (cf. Eq. 26). (Radial-mode damping rate courtesy Günter Houdek.) It was used in the analysis of red-giant mode energetics by Dupret et al (2009) and Grosjean et al (2014).
Gough's theory was generalized to included nonlocal effects by Balmforth (1992a), based on an analysis by Gough (1977a). Specifically, local mixing-length theory implicitly makes the assumption that the relevant convective scales are much smaller than the scale of variation of stellar structure, whereas in fact the typical mixinglength scale is of order a pressure scale height. The nonlocal analysis involves an average over the extent of the motion of convective eddies and over the ensemble of eddies at any given location; this also gives rise to limited convective overshoot and circumvents mathematical problems that occur in a fully local convection formulation when consistently including turbulent pressure in the equation of hydrostatic support. The theory has so far only been developed for radial oscillations; however, given that the relevant effects are concentrated in the superficial layers of the stars, the results are expected to be representative also for non-radial oscillations of low degree, at least in the frequency range of solar-like oscillations. With suitable, and plausible, choice of parameters characterizing the nonlocality this formulation results in oscillation line widths in reasonable agreement with solar observations (e.g., Balmforth, 1992a;Houdek et al, 2001;Chaplin et al, 2005).
Given the fully non-adiabatic results for radial modes, including the effects of turbulent pressure, we can estimate the damping times of the mixed modes by com-bining the radial-mode damping rate 25 with the asymptotic treatment of the damping in the core (cf. Eq. 97). An example is shown in Fig. 39, for a fairly low-luminosity 1 M ⊙ red giant. Here we obtainedη 0 using the Gough (1977b) and Balmforth (1992a) nonlocal treatment of the perturbations to the convective properties, as described by Houdek et al (1999). We note that the acoustic-mode lifetimes are a factor 2 -3 larger than those found by Grosjean et al (2014) and hence in substantially better agreement with observations Corsaro et al, 2015a, Handberg et al, submitted). Comparison with Fig. 38 indicates that in this case the near-surface damping dominates near ν = ν max , even for the most g-dominated modes.

Mode excitation
Stochastically excited mode amplitudes were first estimated by Goldreich and Keeley (1977). Their results were used by Christensen-Dalsgaard and Frandsen (1983b) in a first estimate of the expected mode amplitudes of stochastically excited oscillations across the HR diagram. This was later summarized by  in a widely used relation according to which the mode amplitudes roughly scale as L/M. In particular, it was already then clear that evolved stars were expected to have higher oscillation amplitudes, as has certainly been confirmed observationally.
The stochastic energy input from convection is driven by turbulent Reynolds stresses and entropy fluctuations, their relative importance depending on the detailed assumptions made in the calculation (e.g., Balmforth, 1992b;Goldreich et al, 1994;Samadi and Goupil, 2001;Samadi et al, 2003;Chaplin et al, 2005). However, regardless of these details the dominant contributions come from the near-surface layers where convection is most vigorous. Here the local properties of the oscillations are essentially independent of the degree of the mode, at least for the low-degree modes that are relevant in distant stars, and the rate of energy input is consequently a function of frequency but not degree. It may be shown (e.g., Chaplin et al, 2005) that the resulting mean square amplitude can be expressed as where F is the rate of energy input. The details of this expression, including the precise form of F (ω), depend on the observed quantity represented by A (see also Section 1.3). However, since the ratio between the physical amplitudes of different observable oscillation properties is generally a function of frequency this does not change the form of the equation. From Eq. (116) the rate of energy input can be determined from η A 2 , with η determined from the observed line width and the mode inertia obtained from computed eigenfunctions. In this manner Stein and Nordlund (2001) showed a remarkable agreement between the predictions from three-dimensional simulations of convection and radial-velocity observations of solar oscillations. Also, Jacoutot et al (2008) used such a comparison to constrain the detailed properties of their simulations of solar convection. A similar analysis based on the results from Kepler, including for red giants, would clearly be very interesting. Samadi et al (2007) compared the rate of energy input obtained from hydrodynamical simulations with the results obtained for mixing-length based models. The analysis was extended to red giants by Samadi et al (2012), comparing the results with CoRoT observations analysed by Baudin et al (2011). They noted the importance of the proper conversion of the predicted velocity amplitudes to the photometric data, such as obtained by CoRoT.
As discussed above, the properties of η can be understood in terms of the work integral, Eq. (92). Here (and in the expression for the contribution η t from the turbulent pressure, Eq. 94) the denominator essentially corresponds to the normalized mode inertia E, apart from the normalization with the surface displacement (cf. Eq. 38). Near the surface the integrands in the numerators of Eqs (92) and (94) depend only on frequency. Thus if the contribution from the core to the work integral can be neglected, Eη is just a function of frequency. Consequently, A 2 ∝ E −1 , and the amplitude of a general non-radial mode, relative to a radial mode of the same frequency, scales as A 2 = Q −1 Ā 2 0 (117) (cf. Eq. 40), where Ā 2 0 is the radial-mode amplitude, interpolated to the frequency of the mode considered. However, in the mixed modes in red giants the contribution η core from the damping in the core can be very significant and must be taken into account. Using Eqs (40) and (97) we obtain here, despite the caveat in footnote 10, we used Eq. (42) to replace Qζ by Q − 1, with sufficient accuracy for the present purpose.
In terms of the observed oscillation power spectrum A 2 corresponds to the area under Lorentzian peak (cf. Eq. 89). Of more relevance to the interpretation of the observations is the peak height H, related to A 2 by (see Baudin et al, 2005;Chaplin et al, 2005, and Section 3.4). Relating as in Eq. (118) the peak height to the peak heightH 0 of radial modes, interpolated to the relevant frequency, we obtain In particular, in the lower part of the red-giant branch the core damping is small compared with the envelope contribution (cf. Fig. 38) and hence can be neglected; consequently the peak heights of all mixed modes are predicted to be comparable to the peak height of the adjacent radial modes, in clear contradiction to the observations.
In an important breakthrough in the interpretation of the observations of red-giant oscillations it was noted by Dupret et al (2009) that the origin of this discrepancy was the assumption that all peaks have a Lorentzian profile, leading to Eq. (119). In fact, this is only strictly true if the oscillations are observed for an infinite period. For observations over a finite period T obs the line profile is a combination of a Lorentzian and a sinc 2 function (cf. footnote 4). The broader peaks result in a reduction in H, at the given A 2 , and a behaviour that is qualitatively consistent with the observations. This analysis was extended, including a more detailed comparison with the observations, by Grosjean et al (2014). Fletcher et al (2006) showed that the transition between unresolved and fully resolved peaks can be approximated by replacing Eq. (119) by As a result, Eq. (120) is replaced by where we introduced the interpolated radial-mode damping timet damp,0 = 1/η 0 and peak heightH 0 . For observations extending over 100 -1000 days, as was the case for CoRoT and the nominal Kepler mission, typicallyt damp,0 ≪ T obs and the correction in the numerator in Eq. (122) can be neglected. However, the same is not true of the corresponding term in the denominator for the more g-dominated mixed modes, with Q much larger than 1 (cf. Fig. 17), whose heights are therefore reduced. A further reduction can clearly come from the core damping, represented by η core . It seems probably that this reduction of peak height is the reason that Frandsen et al (2002) only identified the radial modes in 30-day observations of ξ Hydrae, since even the acoustically dominated dipolar modes have a somewhat higher inertia than the neighbouring radial modes.
To illustrate these effects Fig. 40 shows a simplified modelling of power density spectra, for the 1 M ⊙ , 7 R ⊙ model illustrated in Fig. 39, based on the damping times shown there. Here we have approximatedH 0 by a Gaussian centred on the estimated ν max ≃ 68 µHz and with a maximum of 1. The total power is calculated as the sum over the modes of degree l = 0 − 2, each mode represented by a Lorentzian with the width corrected for the finite observing time, as implicit in Eq. (121), and a height given by Eq. (122). In the top panel, for T obs = 100 d, most mixed modes have a strongly reduced height, owing to the relatively short observing time, although some mixed dipolar modes may in principle be visible. At high frequency the peaks are very substantially broadened by the short mode lifetimes (cf. Fig. 39). With T obs = 1000 d, shown in the lower panel, most dipolar mixed modes are excited to substantial heights, and hence one might expect to detect an almost complete dipolar spectrum, at least at moderate and high frequency. However, at low frequency the dipolar modes are strongly suppressed by the increase in the damping rate (note that, according to Eq. (115), η core ∝ ν −2 ) and in particular the small damping rate compared with 2/T obs . For l = 2 the suppression of the peak height, except for the most acoustic modes, is far higher owing to the larger values of Q, and very few if any mixed modes are predicted to be visible.
The strong increase in the damping rates as the star moves up the red-giant branch (cf. Fig. 38) has a dramatic effect on the predicted power density spectra; this is illustrated in Fig. 41 from Grosjean et al (2014), for a 1.5 M ⊙ evolution sequence here assuming T obs = 360 d. 26 For the 5.2 R ⊙ model all dipolar mixed modes are visible. When R reaches 7.3 R ⊙ the heights of the dipolar mixed modes are substantially reduced, although many modes would still be expected to be visible. Finally, at a radius Fig. 41 Simulated power density, in m 2 s −2 µHz −1 , for models in a 1.5 M ⊙ evolution sequence with radii 5.2 R ⊙ (top), 7.3 R ⊙ (middle) and 11.9 R ⊙ (bottom). Adapted from Grosjean et al (2014). of 11.9 R ⊙ , just below the bump, the mixed modes are no longer visible. Although further investigations along these lines, matching also the observed linewidths, are needed, this clearly indicates a serious limitation on the diagnostic possibilities for very evolved stars.
At even higher luminosity, the damping in the g-mode cavity becomes so strong that gravity waves are damped before being reflected from the centre, essentially eliminating the g-dominated mixed modes. In this case the remaining p-dominated oscillation spectrum can be computed for just the envelope model, applying boundary conditions at the edge of the g-mode cavity which select those waves that propagate towards the centre (Dziembowski, 1977;Osaki, 1977;Van Hoolst et al, 1998;Dziembowski et al, 2001;Dziembowski, 2012). It was noted by Dziembowski (2012) that towards the tip of the red-giant branch the resulting loss of wave energy to the core becomes essentially negligible.
As shown by Eqs (118) and (122) the mode amplitudes and peak heights of the mixed modes depend strongly on the mode inertia, in units of the radial-mode inertia. It was pointed out by Benomar et al (2014) that this provides an opportunity to use the observed amplitudes or peak heights as diagnostics of the stellar interior, supplementing the information obtained from the frequencies. In particular, given the sensitivity of ζ, and equivalently Q, on the coupling strength q (cf. Eq. 80) this may provide information about the evanescent region between the acoustic and g-mode cavities. These possibilities should be further explored, given the extended data available from Kepler.

Rotation
A striking result of the early analysis of red-giant observations from Kepler was the detection by Beck et al (2012) of rotational splitting of the observed frequencies, leading to a first inference of the internal rotation of a red-giant star. Extensive results are now available for both red-giant and clump stars (e.g., Mosser et al, 2012b;Deheuvels et al, 2015;Vrard et al, 2016). Strikingly, the inferred core rotation of red giants is far slower than would be expected from models of angular-momentum evolution (e.g., Eggenberger et al, 2012;Marques et al, 2013;Cantiello et al, 2014, see Section 5.5). Here we discuss the effects of rotation on the stellar oscillation frequencies.
For evolved stars it is generally assumed that rotation is so slow that second-order effects of rotation, including the centrifugal acceleration, can be neglected. Then in particular the hydrostatic structure of the star is not affected. The oscillation frequencies are affected by advection of the pattern of waves propagating in the azimuthal direction, such that the frequencies of prograde waves (travelling in the direction of rotation) increase and frequencies of retrograde waves decrease. In addition, the oscillations are affected locally by the Coriolis force in a frame rotating with the star. The result is the first-order rotational splitting, with the frequencies given by (see also Section 3.4.4).
We cannot assume that stars rotate as a solid body. Hence the rotational splitting is an average over the internal rotation rate, determined by the properties of the mode. For simplicity we only consider the case of the so-called shellular rotation where the angular frequency Ω = Ω(r) only depends on the distance to the centre. Then δω nlm does not depend on m, and we can express it as where the rotational kernel K nl is given by Here the first two terms provide a weighted average of Ω and correspond to the advection, while the last two terms arise from the Coriolis force. It is common also to consider such that for uniform rotation δω nl = β nl Ω. For acoustic modes with |ξ r | ≫ |ξ h | we can neglect the last two terms in the numerator of K nl , so that β nl ≃ 1. For general rotation we obtain in this case δω nl ≃ Ω , i.e., an average angular velocity weighted by the local contribution to the inertia. For high-order g modes, on the other hand, with |ξ h | ≫ |ξ r | we can neglect the terms in ξ r to obtain in particular, for dipolar modes, with l = 1, we obtain β n1 ≃ 1/2. To analyse the rotational splitting for red giants, we follow Goupil et al (2013) and use the asymptotic description of the mode inertia discussed in Section 4.2.3. For simplicity, we suppress the subscript 'nl' in the following. Neglecting the term in ξ r ξ h in Eq. (125) and using Eqs (79) and (128), β may be approximated by β ≃ [1 − (l(l + 1)) −1 ]I g + I p I p + I g = [1 − (l(l + 1)) −1 ]ζ + (1 − ζ) = β core + β env ; (129) here the contributions to β from the core and the envelope are and Fig. 42 Computed rotational splittings as functions of ζ for dipole modes, with l = 1, in a 1.3 M ⊙ model of radius 6.2 R ⊙ from Jiang and Christensen-Dalsgaard (2014). An angular velocity was imposed with Ω = 4.98 × 10 −6 s −1 in the core, corresponding to a rotation period of 15 days, Ω = 3.02 × 10 −7 s −1 in the envelope and a narrow transition at the hydrogen-burning shell, with r = 0.09R. All modes with frequency above 0.4ν ac were included. We note the expected linear dependence (cf. Eq. 135). The dashed line shows a uniformly weighted linear least-squares fit to the results.
We furthermore introduce the kernel-weighted core and envelope angular velocities: and Ω env = env Ω(r)K(r)dr env K(r)dr .
These averages formally depend on the mode; however, the dependence on at least the mode order is weak, given the asymptotic description of the modes and assuming that Ω(r) varies smoothly with r on the scale of the eigenfunctions. With these definitions we obtain where R = Ω env / Ω core . Writing the splitting δν = δω/2π in terms of cyclic frequency and specializing to l = 1 we finally obtain , where δν max = Ω core /4π, assuming that the core is rotating substantially faster that the envelope, so that the maximum splitting occurs for the g-dominated modes with ζ ≃ 1. Thus we find that the splitting varies linearly with ζ. This is illustrated for a stellar model in Fig. 42; from the coefficients of the linear fit the core and envelope rotation can be determined with very good precision: the relative errors in the inferred Ω core and Ω env are 0.6 and 6 %, respectively. In analyses of observed power spectra Mosser et al (2012c,b) assumed that the envelope contribution to the splitting could be neglected, corresponding to taking R = 0, and applied an empirically based approximation to the splitting, which can be expressed as δν = ζ Mosser δν max , where 27 where λ and γ are empirically determined parameters. An expression of this form can in fact be obtained from the full asymptotic expression (80), by expanding the inverse, assuming that ζ is close to 1, and expanding sin 2 on the assumption that ν is close to ν (p) n p l . A similar functional form, based on a Lorentzian departure from the asymptotic period spacing, was used by Stello (2012) in the analysis of computed period spacings.
The relative contribution to the rotational splittings of the different parts of the star can also be illustrated by the partial integrals of the kernels (see Beck et al, 2012, supplementary material). These are shown in Fig. 43 for selected p-and g-dominated modes with l = 1 and 2 in a 1.3 M ⊙ red-giant model, normalized such that the surface values are one. This shows that even for the p-dominated modes a substantial contribution to the splitting comes from the radiative core, particularly for l = 1. This is related to the fact that ζ even for these modes is substantially bigger than 0. Given that generally R ≪ 1 for red giants, this makes it difficult to determine the envelope rotation rate; observations of rotational splittings for modes with l = 2 or, even better, l = 3 would be very helpful.
The analysis of rotationally split modes in red giants is complicated by the possible overlap between neighbouring multiplets in the mixed-mode spectrum, when δν max ∼ ν 2 ∆Π. Even so, Mosser et al (2012b) succeeded in determining the rotational splitting in a large number of stars, including cases of overlap, by identifying the pattern of rotationally split components using the approximation in Eq. (136). As pointed out by Mosser et al (2015) the analysis is greatly simplified by using the stretching with the function P introduced by Eq. (83). We first note that, according to Eq. (135), the frequencies of the rotationally split modes can be approximated by where for simplicity we neglected the envelope contribution and hence assumed that R ≃ 0. From this it follows, using also that P is defined as a function of frequency, that the rotationally split modes for a given m are approximately uniformly spaced in stretched period P(ν nlm ), with a period spacing given by where ζ is a suitable average of ζ. In a periodéchelle diagram based on the stretched periods this corresponds of a sequence of modes on a line inclined relative to the modes with m = 0. Two examples of this are shown in Fig. 44, for l = 1 and two different rotation rates. We note that in the right-hand panel the rotational splitting in terms of period is around 62 s, i.e., comparable with the period spacing of 75 s and hence leading to a complex structure of the power spectrum. Even so, the structure of the rotationally split modes can be unambiguously identified in the stretched period echelle diagram. Mosser et al (2015) argued that ζ can be represented by N(ν)/[N(ν) + 1] where N measures the number of gravity modes in a ∆ν-wide interval around ν (see Eq. 57). Assuming also that the modes considered have frequencies close to the frequency ν max of maximum oscillation power they obtain for the period spacings where As a complication in the treatment of red-giant rotational splitting it was noted by Ouazzani et al (2013) that the core rotation may be so rapid that the perturbation analysis discussed here is inadequate. In such cases Ouazzani et al demonstrated that a two-dimensional solution of the oscillation equations is required, involving coupling between components of different degrees. A complication in these cases is also that properties of the eigenfunctions, in particular the mode inertia, can vary substantially between the different components of the multiplet.

Groundbreaking results
Asteroseismology of red-giant stars has been very successful over the past decade with many publications and ground-breaking results. We anticipate that this will continue in the decade(s) to come (see next section). In this section we discuss results that we think have been seminal for the field.
Early observations and interpretations of oscillations in very large red giants such as Arcturus (Merline, 1995) and α UMa (Buzasi et al, 2000;Dziembowski et al, 2001) have already been conducted during the last century and the beginning of this century. Following these discoveries, the study of solar-like oscillations in red-giant stars all along the red-giant and horizontal branch (Sections 2.1.2 & 2.1.4) has effectively started with the spectroscopic campaign in which ξ Hydrae (Frandsen et al, 2002), ǫ Ophiuchi  and η Serpentis ( Barban et al, 2004) were observed. In these observations oscillation power excesses exhibiting a regular pattern of oscillation modes were observed unambiguously for the first time for early red giants. These discoveries led to many questions concerning e.g. whether non-radial oscillation modes would be observable or damped in the core or what the typical lifetime of the stochastically excited and damped modes would be. Data from the dedicated photometric space missions CoRoT (Convection Rotation and planetary Transits; Baglin et al, 2006) and Kepler (Borucki et al, 2008) have been vital in answering these and subsequent questions.

Non-radial oscillation modes
The theoretical work by Dziembowski et al (2001) for α UMa showed that high up on the red-giant branch non-radial modes are strongly damped in the core. Extrapolating from this result, it was initially thought that non-radial modes would be strongly damped in the cores of stars all along the red-giant branch (Christensen-Dalsgaard, 2004). However, this was contradicted by observations. Firstly, Hekker et al (2006) claimed the detection of non-radial modes in red-giant stars. For this claim they analysed the moments of the cross-correlation functions (CCFs) of the spectra of ξ Hydrae, ǫ Ophiuchi and η Serpentis taken during the previously mentioned spectroscopic campaigns (see Section 3.4 for this diagnostic). The variations in the moments as a function of frequency were inconsistent with radial modes, and hence Hekker et al (2006) concluded that these oscillations are non-radial. This was subsequently followed by a similar claim based on MOST (Microvariability and Oscillations of STars; Matthews et al, 2000) data of ǫ Ophiuchi , where the authors extracted frequencies of non-radial modes from the power spectrum. These claims were finally unambiguously confirmed using CoRoT data (De Ridder et al, 2009). The CoRoT data showed that the oscillations in red giants form a regular pattern similar to that seen for the Sun. This includes oscillation modes with significant amplitudes at the expected locations of dipole (l = 1) and quadrupole (l = 2) modes. Hence, these observations confirmed that non-radial modes reach observable heights at the surface of red-giant stars (see Fig. 45

Mode lifetimes
The mode lifetimes of the individual oscillation frequencies (see Fig. 39) are tied to the excitation and damping of the modes. From the first observations of ξ Hydrae the lifetimes and amplitudes of the oscillation modes were examined by e.g. Houdek and Gough (2002) and Stello et al (2006b). Houdek and Gough (2002) reproduced the amplitudes tolerably well with their computations of a stochastic excitation model in which they use a non-local time-dependent generalisation of the mixing-length formulation of Gough (1977b). Subsequently, Stello et al (2006b) estimated the mode lifetime for ξ Hydrae from the scatter of the measured frequencies about a regular pattern. With this method Stello et al (2006b) found a substantially shorter mode lifetime than Houdek and Gough (2002). Data from the CoRoT space mission were again seminal in showing that lifetimes of oscillation modes can vary from tens of days to order 100 of days (see narrow peaks in the left panel of Fig. 45; De Ridder et al, 2009).
Although the physics responsible for the damping mechanism is not yet fully understood (see also Section 4.3), the mode lifetime (τ ∝ 1/Γ, with Γ indicating the FWHM mode linewidth) follows a trend with temperature (see Fig. 46 and Baudin et al, 2011;Corsaro et al, 2012Corsaro et al, , 2015a, and references therein). Interestingly, there seems to be a steeper temperature gradient for the hotter main-sequence stars than for the cooler red-giant stars, whose origin is yet to be explained.

Glitches
The first evidence for the helium glitch in a red-giant star was found using CoRoT data by Miglio et al (2010), who pointed out that the properties of the glitch could be used as additional seismic diagnostics in the global characterization of the star. Broomhall et al (2014) made an extensive analysis of the effects of the helium glitch in red-giant models. They noted the difficulty in obtaining meaningful inferences of the helium abundance, given the limited number of modes available for the analysis. Encouraging results were obtained by Corsaro et al (2015b) for 18 low-luminosity red giants observed for the full Kepler mission; they showed that the formulation derived by Houdek and Gough (2007) gave an excellent fit to the observed glitch signatures and obtained relatively precise determinations of the amplitude and acoustic depth of the glitch signal. Vrard et al (2015) analysed a large sample of Kepler red giants and found interesting differences in the glitch properties of red-giantbranch and clump stars. This may be related to the difference found by Kallinger et al (2012) between these two evolutionary stages in the phase term ǫ (cf. Section 3.3.4), which appears to be caused by differences in the ionization behaviour of helium .

Mixed modes
A first mention of dense and/or irregular frequency patterns in the solar-like oscillations of red giants was made by Hekker et al (2009). These authors already indicated that this could be explained by the fact that the observed oscillations are influenced by their behaviour in both the p-mode and g-mode cavity. Kepler observations were needed to resolve these modes and identify that non-radial modes in red-giant stars are mixed modes, i.e. they propagate in both the outer acoustic cavity as well as in the inner buoyancy cavity. Results of mixed modes for stars observed by both CoRoT and Kepler have been presented by e.g. Beck et al (2011);Bedding et al (2011);Mosser et al (2011aMosser et al ( , 2014 The detection of mixed modes in timeseries data from the Kepler mission for a few hundred stars showed that the period spacing between the mixed modes is a direct measure of whether a star is in the hydrogen shell burning phase or also burning helium in the core Mosser et al, 2011a). Subsequently, Mosser et al (2014) used the combined period spacing and large frequency separation to also identify stars in short evolutionary phases such as the helium subflash stage, pre-secondary clump stars and stars leaving the red clump moving towards the AGB (see Fig. 47). The differences in the period spacings between stars with an inert helium core and stars with core-helium fusion are in part attributed to the presence of a convective core in stars with helium-core fusion (Christensen-Dalsgaard et al, 2014, see also Section 4.2.5).

Radial differential rotation
Using the fact that different mixed modes probe different radial regions in stars, Beck et al (2012) investigated the mixed modes for rotational splittings. They found that rotational splittings show different behaviour depending on whether the rotationallysplit modes are pressure dominated or gravity dominated. Qualitative comparison with models revealed that subgiants and red-giant stars ascending the red-giant branch exhibit radial differential rotation with the core rotating faster than the surface. Mosser et al (2012b) subsequently presented results for an ensemble of stars. They found that there is a small increase of the core rotation for stars ascending the red-giant branch. However they also found an important slow down for red-clump stars compared to the red-giant-branch stars (see Fig. 48). Deheuvels et al (2012Deheuvels et al ( , 2014; Di Mauro et al (2016) subsequently performed rotational inversions to investigate the radial rotation profiles. They showed that the core rotation rate can be determined with only a weak model dependence. Additionally, an upper limit for the surface rotation could be obtained.
Interestingly, from a theoretical point of view the observed core rotation cannot be explained. Current models include transport of angular momentum due to rotationally induced instabilities and circulations as well as magnetic fields in radial zones (generated by the Tayler-Spruit dynamo; Spruit, 1999). However, these models overpredict the core rotation rate by about one order of magnitude (e.g. Eggenberger et al, 2012;Marques et al, 2013;Cantiello et al, 2014), as illustrated in Fig. 49. This leads to the conclusion that an additional angular momentum transport process must be Fig. 49 Evolution of the average core rotational period as a function of stellar radius for different assumptions of angular momentum transport in a 1.5 M ⊙ model initially rotating at 50 km s −1 . Models without angular momentum transport (green), including transport of angular momentum due to rotational instabilities (purple) and accounting for magnetic torques in radiative regions (red, Tayler-Spruit magnetic fields) are shown. The star symbols indicate the locations of KIC8366239 and KIC5006817 as derived using the maximum observed splitting of their mixed modes (Beck et al, 2012(Beck et al, , 2014. Dashed lines indicate a linear fit to the different curves during the early RGB. The vertical dotted line shows the location of H-core exhaustion. The red dotted line shows the evolution of core rotational period for a model where the resulting Tayler-Spruit diffusion coefficient has been multiplied by a factor of 100. Stars in the red giant sample of Mosser et al (2012b) with R < 7.5 R ⊙ are shown as black dots. The best fit to the core rotation of the sample Mosser et al (2012b) is also shown as a dashed blue line. operating that is currently not included in the models. Belkacem et al (2015b,a) investigated the efficiency of mixed modes in extracting angular momentum from the innermost regions of subgiants and red-giants. They concluded that for evolved red giants, mixed modes are sufficiently efficient to balance and exceed the effect of the core contraction, in particular in the hydrogen-burning shell. However, this is not the case for subgiants and early red giants.

Suppressed dipole modes
For a small subset of stars the dipole (l = 1) modes seem to be suppressed, i.e. they have a low height or visibility (V 2 ) in the power spectrum. This was first flagged by Mosser et al (2012a) and followed by a detailed investigation for a single star by García et al (2014). These authors concluded that the low visibility cannot be ex-

Fig. 50
Left: Average visibility of each star vs. ν max , which correlates closely with surface gravity (shown on the top axis) for stars with masses ranging from 0.9 to 2.1 (top to bottom and colour scale). Stellar mass has a formal 1σ uncertainty of 10%. Stars evolve from right to left in the diagrams, corresponding roughly with the beginning of the red giant phase to the red giant luminosity bump. The solid black line shows the theoretical predicted dipole-mode suppression for 1.1M ⊙ , 1.3M ⊙ , 1.5M ⊙ , 1.7M ⊙ and 1.9M ⊙ (top to bottom panels) and a radial-mode lifetime of 20 days. The fiducial dotted lines separate normal and dipolesuppressed stars. Panel taken from Stello et al (2016b). Right: Schematic representation of the magnetic greenhouse effect. Acoustic waves excited in the envelope couple to gravity waves in the radiative core.
In the presence of a magnetic field in the core, the gravity waves are scattered at regions of high field strength. Because the field cannot be spherically symmetric, the waves are scattered to high angular degree and become trapped within the core, where they eventually dissipate (dashed wave with arrow). Panel taken from Fuller et al (2015).
plained by damping effects nor by a fast-rotating core. Recently, Fuller et al (2015) and Stello et al (2016b) performed a joined observational and theoretical investigation providing the first scenario that may explain these suppressed dipole modes: a magnetic greenhouse effect. Essentially, the mixed-mode oscillations in the core get deflected by the magnetic field in the core and scattered to higher degree or converted into Alfvén waves. This prevents them from coupling with an acoustic mode and hence causes them to be trapped in the core. This effective damping suppresses the amplitudes of the observed dipole modes. Fig. 50 shows a schematic diagram of this scenario as well as a comparison of the predicted visibilities from this scenario with observed visibilities of dipole modes as a function of stellar mass. Subsequently, Stello et al (2016a) also investigated the suppression of quadrupole and octupole modes. They found that mode suppression weakens for higher-degree modes with a reduction in the quadrupole mode visibility of up to 49% and no detectable suppression in octuple modes. This is consistent with predictions based on the theory of the magnetic greenhouse effect (Fuller et al, 2015) applied to these higher-degree modes.

Galactic Archaeology
The large number of dwarfs and intrinsically bright red-giant stars for which accurate stellar parameters are now becoming available from asteroseismology has added a new dimension to the field of Galactic Archaeology, i.e. the study of the formation and evolution of the Milky Way by reconstructing its past from its current constituents. The first study using asteroseismic data to probe the galactic disk was performed on a set of about 800 CoRoT stars. This study showed qualitative agreement between the models and observations, but also flagged differences in ν max distributions (Miglio et al, 2009). This was subsequently followed by a population synthesis study using asteroseismic data of about 500 dwarfs observed with Kepler. Again the models and observations were in qualitative agreement, although there were differences present in the mass distributions . These early investigations have initiated major collaborations between asteroseismic experts and experts in galactic astronomy to combine the efforts of large spectroscopic surveys (see Feltzing, 2015, for an overview of the ongoing spectroscopic surveys and the requirements for them) with the asteroseismic measurements. This has led for instance to a direct measurement of a vertical age gradient in the Milky Way  and to the detection of α-rich young stars, i.e. stars that are young according to their asteroseismic measures and old based on their chemical abundance of α elements Martig et al, 2015). First studies investigating these young α-rich stars hint towards these stars being blue stragglers (Yong et al, 2016;Jofre et al, 2016) 5.8 Extra-solar planets The photometric timeseries obtained by the CoRoT and Kepler space instruments are suitable for both exoplanet studies and asteroseismology. The combination of both has been the basis of many ground-breaking planet discoveries (e.g. Carter et al, 2012;Barclay et al, 2013;Campante et al, 2015). Transiting exoplanets leave larger signatures in the timeseries data when orbiting a smaller star, hence one may expect an observational bias against stars with larger radii. Nevertheless, it is possible to detect transiting planets around early red giants. Huber et al (2013) detected a planetary system consisting of two transiting planets around a low-luminosity red giant. Using the rotationally split (mixed) oscillation modes (Section 3.4.4) they also found that the spin axis of the star is not aligned with the orbital axis of the two planets. This could only be explained by a third wide compagnion, which has indeed been detected. Studies like the one by Huber et al. are essential to unravel planet formation scenarios.

Future
The wealth of data obtained by photometric space missions CoRoT (Baglin et al, 2006), Kepler (Borucki et al, 2008) and K2 (Haas et al, 2014) are currently being explored. Exploiting these as well as data from future complementary (space) telescopes such as TESS (Ricker et al, 2014), PLATO (Rauer et al, 2014) and the SONG network (Grundahl et al, 2014) will be essential to study many of the questions that are still open in stellar structure and evolution of giant stars exhibiting solar-like oscillations. These questions include, but are not limited to convection, rotation, (core) overshoot, additional mixing and stellar ages. We discuss these in some detail here.
Convection and surface effects Deep in stars convection leads to a temperature stratification that is essentially adiabatic. Near the surface, where the convective transport is less efficient, a substantially superadiabatic region is present. In stellar models this is typically described with some form of a mixing-length approach (e.g., Böhm-Vitense, 1958;Canuto and Mazzitelli, 1991). To take the next step forward it is important to improve the models to take convection properly into account. Very promising efforts are in progress to perform 3D hydrodynamical simulations of the atmosphere and near-surface part of the interior in which convection is realistically modelled (see for instance Trampedach et al, 2014a,b;Magic and Weiss, 2016). The results of the simulations can be used in stellar modelling through a calibration of parameters of simpler formulations, such as the mixing length . Alternatively, the outer layers of the model can be replaced by suitably averaged versions of the simulations, interpolated to the parameters of the star.
The inadequate modelling of the structure of the outermost layers is an important contribution to the frequency dependent offset between the model frequencies and observed frequencies, i.e. the surface effect. As discussed in Section 3.4 several approaches exist to mitigate this, although these all require adjustments on a star by star basis. Thus a better theoretical understanding of these effects are highly desirable. Replacing the outer layers of the model by an averaged 3D simulation goes some way towards reducing these effects (Rosenthal et al, 1999;Robinson et al, 2003;Sonoi et al, 2015;Ball et al, 2016). An additional contribution comes from non-adiabatic effects and the influence of turbulent pressure which are typically ignored in the modelling. Interestingly, combining 3D simulations with a non-local treatment of time-dependent convection and non-adiabaticity may very substantially improve the agreement between observed and modelled solar oscillation frequencies (Houdek et al. submitted).
Rotation The current stellar models cannot reproduce the core rotation rates observed using rotationally split mixed modes (see Section 5.5). In these models transport of angular momentum due to rotationally induced instability and circulation, as well as magnetic fields in radiative zones generated by the Tayler-Spruit dynamo are included (Cantiello et al, 2014). The fact that these models can still not match the observations is most likely due to missing physics in the models. Observations of (radially differential) rotation in stars across the HR diagram may provide indications of what additional physical effects play a role. Additionally, the impact on red-giants models of the recent results by van Saders et al (2016) and Metcalfe et al (2016) may need to be investigated.  showed that the effective loss of angular momentum ceases above a critical Rossby number (Ro): Ro = rotation period convective turn-over time ≈ 2.1 Subsequently, Metcalfe et al (2016) aimed to explain the underlying reasons for this by proposing a scenario implying a change in the character of differential rotation that ultimately disrupts the large-scale organisation of the magnetic field in solartype stars. This process begins at Ro ≈ 1, where the rotation period becomes comparable to the convective turn-over time. Metcalfe et al (2016) speculated that this may lead to an accelerated decrease in the surface area of spots while the star goes through a rapid phase of magnetic evolution while crossing the Vaughan-Preston gap (Vaughan and Preston, 1980). Due to the changes in magnetic topology during this fast phase of magnetic evolution stars reach Ro ≈ 2 where magnetic breaking operates with a dramatically reduced efficiency . Although this scenario is still speculative the proposed changes in the rotation of main-sequence stars may impact the rotation profiles of red giants. Finally, in terms of observations of the envelope rotation rate it will be of importance to measure rotational splittings of p-dominated l = 2 and l = 3 modes (see Sections 3.4.4 and 4.4) as these have larger sensitivity to the surface layers (see red triple-dot-dashed line in Fig. 43). Radial-velocity observations are more sensitive to these higher-degree modes and hence SONG could play a major role in this.
Additional mixing The standard theory of stellar evolution fails to explain abundance anomalies observed in stars ascending the red-giant branch. Spectroscopic studies showed that when stars reach the bump (see Section 2.1.2) a drop in the surface carbon isotopic ratio, lithium and carbon abundances is present, while nitrogen increases slightly (e.g. Gilroy and Brown, 1991;Tautvaisiene et al, 2013). This is not seen in standard stellar models and provides evidence that an extra-mixing process should occur when low-mass stars reach the bump (e.g. Lagarde et al, 2015, for a study into the effects of rotation-induced mixing using both asteroseismic and spectroscopic constraints).
Asteroseismic measurements of the strength and locations of boundaries between radiative and convective regions that are observable as glitches (see Sections 3.4.5 and 4.2.4) could be very important to determine the strength, location and efficiency of the additional mixing that needs to be added to the models to match both the observed stellar internal structure and the observed chemical yields.
Core overshoot In models the edge of a convective core can be described by different criteria (Section 2.1.1). Whether semi-convection exists or another form of slow mixing is present can currently not be directly constrained from observations (Constantino et al, 2015(Constantino et al, , 2016. The fact that with the mixed modes it is possible to probe the stellar core (e.g. Bedding et al, 2011;Mosser et al, 2014) and that from model investigations there are prospects to investigate sharp features around the core  will be essential to investigate the edge of a convective core. The glitches (Sections 3.4.5 and 4.2.4) in the mixed modes may be the only source of direct observational evidence that can be obtained for the region where core overshoot may take place. A full exploitation of the archival data from Kepler may provide these observational constraints on the core overshoot processes.
Suppressed dipole modes As described in Section 5.6 for a fraction of the observed red-giant stars the dipole modes are suppressed. Currently there exists one scenario (Fuller et al, 2015) to explain the presence of these suppressed dipole modes. For this scenario, the predictions are consistent with observations (Stello et al, 2016b, see also Fig. 50). However, the presence of large magnetic field in the core can currently not be tested directly, while further theoretical elaboration of the magnetic model is required (e.g. Cantiello et al, 2016). Hence, it is hard to confirm or reject this scenario at the moment. In this respect, it may be interesting if other scenarios that can be confirmed/rejected or correlations with other stellar parameters can be established.
Ages There is no observable that is sensitive to age and age only (Soderblom, 2010), so all measures of stellar ages are either empirical (with a model dependent calibration) or model dependent. Hence, any improvement in the models and/or matching the data with models by mitigating the surface effect will influence the age determinations of stars. This is particularly relevant for evolved stars as uncertainties in earlier evolution phases accumulate. Despite these limitations, it may be possible to find combinations of (asteroseismic) observables that are more sensitive to age and to improve age estimates using asteroseismology.
Red-clump stars Red-clump stars are both from an observational as well as from a modelling point of view not well understood. Red-clump models carry information from the main-sequence stellar structure. Even a small change in the stellar structure at the main sequence can lead to significant differences in terms of effective temperature, luminosity and internal structure after the onset of helium-core burning. Additionally, the subflashes that are predicted to occur at the onset of helium-core burning for low-mass stars with degenerate cores (see Section 2.1.3) leave sharp features in the model internal structure. This results in a very irregular pattern of frequencies and period spacings for stars entering the red clump (Constantino et al, 2015, see also Section 4.2.5). Information on the reality and nature of the subflashes from observations will be essential to better understand the red-clump models. This may be obtained from stars that are descending the red-giant branch after helium ignition. However, this is challenging as this is a very short phase in evolution, reducing the probability that any observed star is in that phase. Additionally, stars already in the red clump may provide the observational evidence for the presence of the flashes.
However, from an observational point of view, the combination of period spacings, rotational splittings and noise that seems to be present in the observed power density spectra make it very difficult to disentangle all individual oscillation modes. Hence, observational evidence of regularity or irregularity is still awaiting for many of these stars. We anticipate that with the long timeseries of Kepler data it will be possible to disentangle these observational features, but a detailed mining of the data with tailormade tools will be required.
The future for asteroseismology of giant stars with solar-like oscillations looks very promising. In addition to the wealth of archival data there are two new space missions planned which will again provide complementary data in terms of nearly allsky coverage with TESS and brighter stars with PLATO. Additionally, the Stellar Oscillations Network Group (SONG) has one node fully operational and one node in development. SONG provides spectroscopic data from the ground, which will be valuable to study higher degree modes, as well as stars with long periods. These observational data together with improved 1-D and 3-D models are bound to improve significantly our knowledge of stellar structure and evolution of giant stars over the next decade(s).