Asteroseismology of solar-type stars

Until the last few decades, investigations of stellar interiors had been restricted to theoretical studies only constrained by observations of their global properties and external characteristics. However, in the last thirty years the field has been revolutionized by the ability to perform seismic investigations of stellar interiors. This revolution begun with the Sun, where helioseismology has been yielding information competing with what can be inferred about the Earth's interior from geoseismology. The last two decades have witnessed the advent of asteroseismology of solar-like stars, thanks to a dramatic development of new observing facilities providing the first reliable results on the interiors of distant stars. The coming years will see a huge development in this field. In this review we focus on solar-type stars, i.e., cool main-sequence stars where oscillations are stochastically excited by surface convection. After a short introduction and a historical overview of the discipline, we review the observational techniques generally used, and we describe the theory behind stellar oscillations in cool main-sequence stars. We continue with a complete description of the normal mode analyses through which it is possible to extract the physical information about the structure and dynamics of the stars. We then summarize the lessons that we have learned and discuss unsolved issues and questions that are still unanswered.


Introduction
Helio-and asteroseismology allow us to study the internal structure and dynamics of the Sun and other stars by means of their resonant oscillations (e.g. Turck-Chièze et al., 1993;Christensen-Dalsgaard, 2002;Aerts et al., 2010;Basu, 2016, and references therein). These vibrations manifest themselves by motions of the stellar photosphere and by temperature and density changes implying modulations of the positions of the Fraunhofer lines and of the stellar luminosity respectively.
Repeated sequences of stochastic excitation and damping by turbulent motions in the external convective layers lead to a suite of resonant modes in the Sun (Goldreich and Keeley, 1977;Goldreich and Kumar, 1988;Balmforth, 1992;Goldreich et al., 1994;Samadi and Goupil, 2001;Belkacem et al., 2008). The stars where their modes are excited in this way are usually called "solar-like pulsators" or simply "solar-like" stars even though their structure and dynamics could be different compared to the actual Sun, covering main-sequence (MS), sub-giant and red-giant stars (e.g. De Ridder et al., 2009;Bedding et al., 2010a;García and Stello, 2015;Hekker and Christensen-Dalsgaard, 2017). The oscillation periods of solar-like stars range from minutes to years (e.g. Mosser et al., 2010Mosser et al., , 2013Stello et al., 2013Stello et al., , 2014Chaplin et al., 2014a).
In this review, we focus on "Solar-type" stars defined as cool main-sequence dwarfs located below the red edge of the classical instability strip (spectral types from late F, G and K dwarfs, see the zone encircled by the red circle in Figure 1) with a near-surface convective zone that excites the oscillation modes. However, to put these stars in context, we will sometimes discuss sub-giants and early red giants. In such a way, the continuity towards more evolved stars will be ensured.
There are other mechanisms exciting stellar oscillations in more massive and luminous mainsequence stars: a) the heat-engine mechanism (also known as κ or opacity-driven mechanism), related to the changes in the opacity profile due to temperature variations, and responsible for the pulsation in the instability strip and white dwarfs (e.g. Eddington, 1926); b) the mechanism, where the nuclear reaction rate changes as a consequence of the contraction and expansion of the star (e.g. Lenain et al., 2006).; c) tidal effects, where non-radial oscillations can be forced in stars belonging to multiple systems (e.g. Welsh et al., 2011). All of these pulsating stars are usually referred to by the generic term "classical pulsators" (e.g. δ Scuti, γ Doradus, RR Lyrae, Cepheids, etc) and their study is out of the scope of this review (for more information on these variable stars see, for example, Aerts et al., 2010).

Asteroseismology of solar-type stars in a helioseismic context
The best known star representative of solar-like stars is the Sun. Over the last 30 years, helioseismology has proven its ability to study the structure and dynamics of the solar interior in a stratified (layer-by-layer) and latitudinally varying way (see Figure 2). Starting from the photosphere, the internal structure is divided into a convective envelope followed by a radiative zone that includes the inner core where the nuclear reactions to transform Hydrogen in Helium take place. In more massive stars ( 1.2 M ) a convective core develops with a total size that is not well constrained yet (e.g. Saslaw and Schwarzschild, 1965;Zahn, 1991;Deheuvels et al., 2016).
Seismic analysis tools were first applied to our closest star, the Sun, in order to infer its radial and latitudinal internal structure and dynamics. Therefore, the sound-speed profile (e.g. Basu et al., 1997;Turck-Chièze et al., 1997), the density profile (e.g. Basu et al., 2009), the internal rotation in the convective zone (e.g. Thompson et al., 1996) and in the radiative region (e.g. Elsworth et al., 1995;Basu et al., 1997;Chaplin et al., 1999a;Couvidat et al., 2003;García et al., 2004a;Eff-Darwich et al., 2008;García et al., 2008c) or the conditions and properties of the solar core (e.g. Turck-Chièze et al., 2001García et al., 2007García et al., , 2008aBasu et al., 2009;Fig. 1 Kiel diagram, i.e. the logarithm of the surface gravity, log g, as a function of the effective temperature, T eff . The diagram is color coded by the logarithm of the number of stars, N, observed by Kepler . Open red circles locate MS and sub-giant stars where pulsations have been measured by Kepler (Chaplin et al., 2014a). In this review, we focus on cool MS solar-like stars that will be referred as solar-type stars. Appourchaux et al., 2010) have been studied and well determined. Moreover, the characterization of the p-mode properties has led to the determination of other quantities such as the position of the base of the convection zone (e.g. Christensen-Dalsgaard et al., 1985;Ballot et al., 2004b) or the Helium abundance (e.g. Gough, 1983;Vorontsov et al., 1991) with high precision. With all of these observational constraints, the standard solar and stellar evolution models have been significnatly improved, reducing the uncertainties in the calculation of the stellar ages when individual p-mode frequencies are considered (see for more detail the reviews by Lebreton et al. (2014a) and Lebreton et al. (2014b)). However, new asteroseismic observations of many other stars (e.g. Chaplin et al., 2011b;Huber et al., 2011;Lund et al., 2017) covering a larger fraction of the H-R diagram, allow us to test stellar evolution under many different conditions (e.g. Christensen-Dalsgaard and Houdek, 2010) while putting the Sun in its evolutionary context. In asteroseismology, due to the absence of spatial resolution in the observations, only lowdegree modes (those with a small number of nodal lines at the surface of the star) are measured. Therefore compared to the Sun, less detailed information is available on stellar interiors. On the other hand, some pulsating solar-like stars offer the possibility to observe mixed modes, i.e., modes with mixed character resulting from the coupling between p and g modes (Arentoft et al., 2008;Bedding et al., 2010b;Chaplin et al., 2010;Deheuvels et al., 2010a;Beck et al., 2011). The study of these modes allows us to better constrain the structure and dynamics of the deep radiative interiors (e.g. Deheuvels et al., 2010b;Metcalfe et al., 2010;Bedding et al., 2011). Unfortunately, neither mixed-modes nor pure g modes have been identified individually in main-sequence solar-like stars so far because they become evanescent in the convective region and their surface amplitudes are Fig. 2 Artistic view of the structure of a typical 1 M cool main-sequence solar-like pulsating star where two modes are propagating, a low-degree p mode and a g mode. The structure of the star is divided in two main regions: the inner radiative zone (including the core where the nuclear reactions take place) and the external convective zone. In the Sun, the convective envelope extends inward from the surface 30% of its radius.
small compared to the granulation signal. Thus, all of the information that we are obtaining for these stars comes from the characterization of p modes. However, it is important to note that for the special case of the Sun, the global signatures of the dipolar g modes have been measured  with GOLF/SoHO, as well as some individual low-frequency g modes through the study of the perturbations induced by the g modes on the acoustic modes (Fossat et al., 2017). Both results are still controversial as shown for example by Schunker et al. (2018), who demonstrated that the latter detection of individual modes is highly dependent on the selection of the parameters used in the analysis.
Stars are also known to be rotating magnetic objects. Such dynamical processes affect the internal structure and evolution of stars (e.g. Brun et al., 2004;Zahn et al., 2008;Duez et al., 2010;Eggenberger et al., 2010;van Saders and Pinsonneault, 2013;Eggenberger, 2013). Thus, it is necessary to go beyond the classical modeling of stellar interiors and evolution by taking into account transport and mixing mechanisms both on dynamical and secular time-scales (e.g. Zahn, 2004, 2005;Turck-Chièze et al., 2010;Tayar and Pinsonneault, 2013;Eggenberger et al., 2017). A new generation of stellar models is fundamental to understand current and future observations of stellar surfaces and interiors. Indeed, more and more constraints on the stellar rotation profiles are being obtained (e.g. Mathur et al., 2008;Turck-Chièze et al., 2010;Deheuvels et al., 2012Deheuvels et al., , 2014Benomar et al., 2015;Di Mauro et al., 2016;Nielsen et al., 2017), while anomalies at stellar surfaces are found (e.g. Mathur et al., 2007;Zaatri et al., 2007). Furthermore, hints on magnetic fields either at the surface (e.g. Jouve et al., 2010) or in the inner cores of the stars (Stello et al., 2016) have been suggested and could be due to on-going dynamos developing in the convective cores of stars above 1.2 to 1.4 M (see Fig. 3).

Asteroseismic observations of solar-like stars
The requirements needed to perform asteroseismic studies of distant stars are shared with helioseismology and any other seismic studies. Stable and uninterrupted observations are ideal because Fig. 3 Artist's view of the structure of an F-type main-sequence solar-like star with a convective core (Mass 1.2 M ) and two possible dynamos, one at the external convective zone and another one in the convective core.
most of the analyses are performed in the frequency domain, requiring long observations to increase the frequency resolution.
When preparing the observations it is also mandatory to choose a sampling rate rapid enough that the Nyquist frequency is well above the acoustic cut-off frequency of the oscillation modes. Conversely, when the stellar oscillations are just above the Nyquist frequency, aliased peaks are reflected from the Nyquist frequency leaking into lower frequencies. In this case, it is still possible to do asteroseismology for "super-Nyquist" oscillations as first shown by Murphy et al. (2013) and then applied to solar-like stars by Chaplin et al. (2014b). Low-mass cool main-sequence and sub-giant stars have a frequency of maximum power above ∼ 500-8000 µHz, i.e., in the range of ∼2 to ∼30 minutes. Therefore, a sampling rate faster than ∼1 minute is recommended to avoid dealing with super-Nyquist asteroseismology.
Continuity is needed to reduce the effect of gaps in the data. In particular, regular gaps -seen as a Dirac Comb function -should be avoided. Regular gaps are typical in single ground-based observations due to the day/night cycle or from a satellite due to regular operations such as angular momentum dumps of the reaction wheels used to stabilize the spacecraft. When regular gaps are present in the time series, the power of every stellar oscillation peak leaks into surrounding sidelobes due to the convolution by the Fourier transform of the signal with the window function. Examples of the impact of the NASA Kepler window function on stellar oscillations can be found in García et al. (2014b). If the gaps are sparse, the level of noise in the spectrum increases as a function of the duty cycle (see examples in Pires et al., 2015) and the signal-to-noise ratio decreases.
Finally, stable instruments are necessary to minimize any possible instrumental modulations that could generate peaks in the same frequency domain as the expected stellar oscillations. In the case of multi-site observations, it is recommended to have instruments as similar as possible. However, global asteroseismic observing campaigns have shown that it is possible to use very different instruments, normalize the data, and recover the stellar pulsations (e.g. Bedding et al., 2010b).
Oscillations on the Sun and stars can be measured using two fundamental techniques: by measuring the Doppler velocity of the surface of the stars or by measuring the luminosity changes induced by the changes in temperature due to the pulsations. By using the instruments on board the ESA/NASA Solar and Heliospheric Observatory (SoHO, Domingo et al., 1995), it is possible to compare the solar power spectrum measured in Sun-as-a-star Doppler velocity by the Global Oscillations at Low Frequencies (GOLF, Gabriel et al., 1995Gabriel et al., , 1997 instrument and using integrated photometry by the Sun photometers of the Variability of solar IRradiance and Gravity Oscillations instrument (VIRGO/SPMs Fröhlich et al., 1995). The power spectral density obtained from the two instruments is shown in Fig. 4. To facilitate the comparison at all frequency ranges, the power spectral density has been normalized in such a way that the maximum of the p-mode bump has the same amplitude in both observables. The convective background is higher in photometry than in Doppler velocity (more than an order of magnitude at 0.5 mHz). Therefore, the signal-to-noise ratio of the acoustic modes is smaller in photometry (a factor ∼30) while it reaches a factor of ∼300 in Doppler velocity. Using Doppler velocity, it is then possible to characterize a larger number of modes at low frequency. Fig. 4 Comparison between the power spectrum density extracted from 21 years of Doppler velocity (by GOLF, in black) and from photometric measurements (by VIRGO/SPM green channel, in green).
Observational asteroseismology of main-sequence cool dwarfs developed during the 1990s and the first years of the twenty-first century. The first solar-like star for which pulsations were observed was α Cen A. It was first observed in photometry by Buzasi (2000, 2001) using the Wide-Field Infrared Explorer satellite (WIRE, Buzasi et al., 2000;Fletcher et al., 2006). Later, it was re-observed using Doppler velocity techniques from the ground (Bouchy and Carrier, 2001;Martić et al., 2001;Bouchy and Carrier, 2002).
Two other more evolved stars (sub-giants) were studied at that time too: α CMi (Procyon) and η Boo. Procyon was observed by several teams from the ground using Doppler-velocity measurements (Brown et al., 1991;Mosser et al., 1998;Martić et al., 1999;Bouchy et al., 2004). Procyon was also studied in photometry from space (Matthews et al., 2004) by the Canadian satellite Microvariability and Oscillations of STars (MOST, Matthews, 1998;Guenther et al., 2007), providing some controversial results (see for example the discussions in Bedding et al., 2005;Régulo and Roca Cortés, 2005;. Today, however, mode detection for Procyon is well established thanks to a multi-site ground-based campaign (Arentoft et al., 2008;Bedding et al., 2010b). Several individual modes were identified and the internal structure of the star was extracted using these asteroseismic constraints. η Boo was first asteroseismically observed in the equivalent width of the Balmer lines by  (re-observed in radial velocity by Kjeldsen et al. (2003)) and its oscillations were independently confirmed by the MOST space photometric observations (Guenther et al., 2005).
Although the highest signal-to-noise asteroseismic results are obtained by observing stars in Doppler velocity, most of the pulsating main-sequence, solar-like stars have been observed using photometric techniques. Indeed, photometric instruments have better performance outside the Earth's atmosphere. They require fewer photons per star allowing a high sampling rate while keeping the telescope size small. Moreover, many objects can be studied at a time. To give an example, the NASA Kepler mission (Borucki et al., 2010) allowed observing 512 stars at any time with a short cadence of 60 s.
Modern space-based asteroseismology of solar-type stars and sub-giants started with the Convection Rotation and Planetary Transits space mission (CoRoT Baglin et al., 2006), which observed around a dozen such targets. Originally, the objective was to study hot F stars because they were expected to have higher amplitude modes when compared to G and K dwarfs. Unfortunately, the widths of the modes in these stars were also very large. Although this was expected, it was found that the modes overlapped each other and it was extremely difficult to properly identify the modes and to extract precise p-mode parameters (see the discussions in Appourchaux et al., 2008;Benomar et al., 2009b). Hence, the observing strategy evolved and cooler stars were prioritized. The same strategy was then followed later with Kepler.
So far, cool main-sequence dwarfs and early sub-giants have been observed during five space missions: WIRE, MOST, CoRoT, Kepler, including its second's life as K2 (Howell et al., 2014), and the Transiting Exoplanet Survey Satellite (TESS, Ricker et al., 2014). TESS is primarily a mission for sub-giant stars as demonstrated by the theoretical studies already done (Campante, 2017;Schofield et al., 2019) and corroborated by the first marginal detection of pulsations in the solar-type star π Mensae (Gandolfi et al., 2018) and the clear detection of pulsations in the late sub-giant TOI-197 (Huber et al., 2019). The small fraction of solar-type stars observed by TESS will be extremely useful as these targets will be very bright and ground-based complementary studies will contribute to better characterizing them. An additional space-based observatory, the ESA M3 Planetary Transits and Oscillations of stars (PLATO) mission (Rauer et al., 2014) is expected to be launched around 2026. From the ground, the SONG network (Grundahl et al., 2011) is already running with two sites, Spain and China, with a site in Australia expected to be operative in 2020. SONG will also be able to study solar-type stars although it will be best suited for sub-giants and red giants (e.g. Grundahl et al., 2017;Arentoft et al., 2019).
3.1 Structure of the power spectrum density of a solar-type star Asteroseismic analyses are mostly performed in the Fourier domain by computing the amplitude or power spectrum density. An example of the power spectrum density of HD 52265, a G0V star observed over 117 days by CoRoT, is shown in Fig. 5. Depending upon the frequency range to be explored, the spectrum is dominated by features related to different physical phenomena.
Starting from the low-frequencies (<10 µHz), the spectrum is dominated by a series of high peaks and their harmonics. These peaks correspond to the surface differential rotation of the star because of the modulation induced in the mean stellar luminosity by dark spots crossing the visible face of the stellar disk (e.g. Berdyugina, 2005). The average surface rotation of this star is 12.3 ± 0.15 days . At higher frequencies, between 50 and 1000 µHz, the spectrum  . Physical phenomena associated with each region of the PSD are indicated: photon noise, oscillations, convection (granulation), activity related slope, and rotation through the spot modulation of the emitted stellar flux. The continuous red line represents the fitted background components. The blue continuous line is the gaussian fit over the p-mode hump.
is dominated by a continuum , which is the result of the turbulent movements at the surface of the star due to convection, such as granulation or supergranulation. At even higher frequencies, the p-mode envelope is visible. For this star, the acoustic modes are centered around 2000 µHz, i.e., around 8 minutes. For reference, the oscillations of the Sun are centered around 3000 µHz (5 minutes). Finally, close to the Nyquist frequency, the spectrum is flat and it is dominated by the photon noise of the instrument. This noise level depends on the properties of the instrument and it would eventually be above the p-mode hump in stars for which the modes have low amplitudes or when the star is distant and faint.
Examining the frequency range near the p-mode envelope, one may see that it is composed of a sequence of peaks following a repetitive pattern as shown in Fig. 6 for the Kepler target 16 Cyg A Davies et al., 2015). Two main regularities can then be extracted: the large frequency separation, ∆ν, and the small frequency separation, δν, or simply the large and small separations (see Sec 4.4 for more details). A third global seismic parameter, ν max , can also be defined as the frequency at maximum power of the p-mode envelope (Brown et al., 1991;.

Theory of oscillations
After describing the different observations, we present in this section the theory developed to interpret the stellar oscillation spectra. We focus here on theoretical concepts that are useful for this review. More detailed descriptions may be found in, e.g., Cox (1980); Unno et al. (1989);Christensen-Dalsgaard (2002); Aerts et al. (2010). Solar-like oscillations are resonant modes, resonances occurring at specific frequencies. These oscillations are small enough to be considered as linear perturbations around the equilibrium state of the star. Thus, studying the oscillations boils down to an eigenvalue problem; by solving it, we get a discrete set of eigenfrequencies, each one being associated with an eigenmode describing the distribution of the perturbation inside the star. Stellar oscillation modes are standing waves in a meridional plane and may be propagating in the azimuthal direction, as in a waveguide. Assuming that solar-like stars are purely spherically symmetrical objects, i.e. all of the quantities describing the equilibrium depend on the radius r only, the problem is separable and the horizontal part of the modes are described by spherical harmonics (denoted Y m ); only the radial part depends on the stellar structure. A mode is then fully characterized by three integers: the radial order, n, indicates the number of nodes along the radius. By convention, we denote the p modes by positive numbers and g modes by negative ones (see Sect. 4.1). -The angular degree, , is a non-negative integer denoting the number of nodal lines at the surface of the sphere. Thus, modes with =0 are radial modes while those with ≥ 1 are the non-radial modes. More specifically =1 are called dipole modes, those with =2 are the quadrupole modes, the =3 are the octupole modes, etc. -The azimuthal order, m, gives the number of nodal lines passing through the poles. It can take values from − to + including zero. Positive and negative values corresponding to retrograde and prograde waves respectively. When m = 0, modes are axisymmetric. These modes are usually called zonal modes; modes with |m| = are called sectoral modes.
In Fig. 7, a few modes are represented: low-degree modes as well as a higher degree mode.
ℓ 0 1 2 m 0 1 2 Fig. 7 Left: Example of modes of degree =0, 1, 2 and azimuthal order m=0, 1, 2. Modes are represented by spherical harmonics. The blue regions are those moving toward the observer, while the red regions represent those that are moving away; half a period latter this is reversed. Right: Mode =20, m=16 and n=14.
We usually denote the frequency ν n, ,m = ω n, ,m /2π, expressed in Hz, where ω n, ,m is the angular frequency in rad s −1 . Due to the spherical symmetry, modes are degenerate relative to m, making the frequencies independent on m. Thus, we can write ν n, ,m = ν n, . By considering such a symmetry, we implicitly neglect the impact of stellar rotation. We will introduce rotation in Sect. 4.2.

Acoustic, gravity and mixed modes
The equations for oscillations in non-rotating spherical stars lead to a collection of 1D eigenvalue problems that can be independently solved for each value of . Assuming adiabatic perturbations, these equations are fourth order and may be solved numerically given a radial model of the stellar structure. Various codes have been developed, such as ADIPLS, LOSC, PULSE, POSC, GraCo, NOSC, OSCROX, FILOU, LNAWNR, or GYRE (e.g. Christensen-Dalsgaard, 2008b;Scuflaire et al., 2008;Brassard and Charpinet, 2008;Monteiro, 2008;Moya and Garrido, 2008;Provost, 2008;Roxburgh, 2008;Suárez and Goupil, 2008;Suran, 2008;Townsend and Teitler, 2013). However, to understand the nature of the oscillations, we make some approximations. By assuming that the modes vary much more rapidly than the equilibrium structure, we may crudely reduce the problem to a classical turning-point wave equation (1) In this equation ξ r is the amplitude of the radial displacement, ω the angular frequency of the wave, c, N and S = ( + 1)c/r are the sound speed, the Brunt-Väisälä frequency and the Lamb frequency. Within this approximation, waves can only propagate in the stellar interior where K(r) > 0, that occurs where either (i) ω > N and ω > S or (ii) ω < N and ω < S . The first case corresponds to acoustic or pressure (p) modes and the second one to gravity (g) modes. P modes are then high-frequency modes whereas g modes are low-frequency (long-period) modes. In a main-sequence solar-like star these two frequency domains are well separated. Figure 8 shows the profiles of N and S in a solar-type star. We deduced from this plot that g modes are confined in the internal radiative region of stars whereas p modes are mainly confined in the outer part of the star. Modes are also characterized by their inertia I defined from the eigenmodes as where ξ n, ,m denotes the total mode displacement normalized by its value at the stellar surface. Thus, more energy is required to excite higher inertia modes. Since the g modes are largely restricted to the deeper layers of the stars, their inertia is far larger than that of p modes. The cutoff-frequency profile ωc near the solar photosphere, computed within the isothermal approximation, is plotted with a dashed line. The horizontal dashed-dot line depicts the cutoff frequency of modes (∼5600 µHz). At higher frequency, waves are not trapped, whereas at lower frequency, waves are trapped below an external turning point (indicated with vertical dotted lines for different frequencies shown with horizontal solid lines).

p modes
As indicated by their name, p modes are pressure, or acoustic, waves for which the restoring force arises from the pressure gradient. They are the most important modes in solar-like star seismology since they are by far the most observed ones, with periods of several minutes. Normally stable regarding non-adiabatic processes such as κ mechanism, these modes are stochastically excited by the turbulent convective envelope. Observed modes correspond to high-order low-degree modes. The asymptotic (in the sense that n ) theory of p modes has initially been developed by Vandakurov (1967) and continued in the 80's by Tassoul (1980) and Deubner and Gough (1984) for example. Asymptotic theory described also the structure of p-mode spectra. The theory predicts regular patterns: modes are organized in a comb structure as observed. These different regularities are detailed in Section 4.4. P modes are confined in the outer region of the stars. Using Equation 1, their resonant cavity is limited at the top by the stellar photosphere and at the bottom by an internal turning point r t defined such as S (r t ) = ω n, , i.e. r t = c(r t )L/ω n, , where considering L = + 1/2 is a better approximation (e.g. Vorontsov, 1991;Lopes and Turck-Chieze, 1994). As a consequence, we see that for the same radial order, lower-degree modes propagate deeper into the interior of the star, while for the same degree, higher frequency modes, which correspond to higher order modes, also penetrate deeper than lower-order, lower-frequency modes (see Fig. 9).

Fig. 9
Fractional radii of the inner turning point for low-degree acoustic modes obtained from GOLF (García et al., 2008c). The error bars of the modes have been multiplied by 5000 to be visible.
The outer turning point, close to the photosphere, also depends on the frequency. It corresponds to the radius r e where an acoustic wave propagating outwards is reflected due to the brutal drop in density and pressure. A wave propagates as long as ω > ω c where ω c is the cut-off frequency. We then define r e such that ω n, = ω c (r e ). In an isothermal atmosphere, where H is the classical density scale height, that is equal to the pressure scale height H p = (−d ln p/dr) −1 in an isothermal atmosphere. Deubner and Gough (1984) proposed a more detailed expression that is well approximated by this one. Figure 8 shows the profile of ω c for the Sun. In the upper layers of solar-like stars, the profile ω c (r) increases with r and reaches a maximum. This has two consequences: (i) above a cutoff frequency (∼ 5600 µHz for the Sun) acoustic waves are no more reflected, then mode trapping is no longer possible. Nevertheless, we may observe a modelike pattern in the oscillation spectrum above the cutoff frequency (Jefferies et al., 1988;Libbrecht, 1988;Duvall et al., 1991). These so-called pseudo-modes or high interference peaks (HIPS) were first observed in Sun-as-a-star solar observations using GOLF by García et al. (1998) and then confirmed by BiSON . A few years later, in 2007, there were also found in solar-like stars from ground-based observations (Karoff, 2007) and from space using Kepler (Jiménez et al., 2015). From disk-integrated observations, these HIPS may be explained as an interference of high-frequency waves partially reflected at the unobserved far side of the observed star (García et al., 1998). (ii) Higher frequency modes reach higher layers in the atmosphere and this does not depend on ; as a consequence any perturbation in the outer layers of the star affects similarly low-degree modes with very close frequencies. These layers near the photosphere are known to be poorly modelled with 1D stellar evolution codes, since they are highly turbulent, nonadiabatic and with a low plasma beta. All of these phenomena generate the so-called near-surface effects visible as a departure between computed and observed frequencies. Such effects have to be taken into account when analysing Solar-type stars. We can also build diagnosis tools that cancel out this effect (see Sect 4.4).

g modes
At the low-frequency part of the spectrum of solar-like stars, we find g modes, for which the restoring force is the buoyancy due to density fluctuations and gravity. These low-frequency modes are confined in the radiative interior of solar-like stars since gravity waves can only propagate in a region where N 2 > 0 (see Eq. 1), i.e. in non-convective zones by definition. They are thus evanescent in solar-like envelopes and reach the surface with very low amplitudes relative to the p modes. For the Sun, expected periods for g modes are of hours and longer than 35 minutes for the shortest period. The g-mode spectrum also has a regular pattern: g modes of the same degree are evenly spaced in period, the period increasing when the absolute value of the radial order increases. Moreover, radial gravity modes do not exist. Up to now, the only claims of g-mode detection in solar-like stars concerns the Sun (e.g. Turck-Chièze et al., 2004;García et al., 2007;Fossat et al., 2017).

Mixed p/g modes
When solar-like stars reach the end of the main sequence, due to the build-up of strong density gradients in the core, the Brunt-Väisälä frequency increases there. As a consequence there exists a frequency range where g modes in the core and p modes in the envelope may coexist. If the evanescent region between the p-and g-mode cavities is small enough, a coupling between them occurs. In such a case we get mixed modes with both p and g characteristics. Their properties were initially discussed theoretically by Scuflaire (1974) and the first observations of mixed modes in solar-like stars were reported from ground-based observations of η Bootis, by  and confirmed later by Kjeldsen et al. (2003), and Carrier et al. (2005). They were in very good agreement with theoretical predictions by, e.g., Christensen-Dalsgaard et al. (1995). More recently, many observations of mixed modes have been made from ground-based observations but also from space thanks to CoRoT (e.g. Deheuvels et al., 2010a) and Kepler observations Campante et al., 2011;Mathur et al., 2011a). Beck et al. (2011) first reported the existence of mixed modes in red-giant stars. Later, Bedding et al. (2011) and Mosser et al. (2011a) showed the power of these modes to measure the evolutionary status of red giants, with a clear difference between stars ascending the red-giant branch (RGB) and those in the so-called "clump". Mixed modes are very useful to put constraints on the internal structure and dynamics of stars since they are very sensitive to the core due to their g-mode behaviour, whereas their p-mode properties make their surface amplitudes high enough to be detected.

Effects of rotation
In previous sections, we have assumed a perfect spherical symmetry implying that modes are degenerate in m. Everything that breaks this symmetry lifts this degeneracy. Rotation is the most common phenomenon breaking the symmetry. When the rotation rate is slow enough, as it is in most solar-type stars, it may be treated as a perturbation of the non-rotating case. Modes are thus split into 2 + 1 m-components forming a multiplet. Each component has a frequency ω n, ,m = ω n, + δω n, ,m where the second term is a small perturbation usually called a rotational splitting (or just a splitting). This quantity will directly depend on the rotation rate of the star. To first order, in the perturbative expansion, the splittings may be expressed as: where K n, ,m is the rotational kernel depending on the eigenmode and Ω is the rotation profile in the star (Hansen et al., 1977). Thus, the kernel indicates how the splitting is sensitive to the rotation inside the star. By assuming the rotation is solid in the mode cavity, this equation is simplified (Ledoux, 1951): where C n, is named the Ledoux (or Coriolis) coefficient. For p modes C n, ≈ 0 for high-order modes (for the Sun, with n 10 we get C n, 10 −2 ), whereas for g modes C n, ≈ 1/[ ( + 1)]. As a consequence, p modes of solar-like stars are generally modelled as symmetrical multiplets of frequencies where ν s = Ω/(2π) is simply called the splitting. If the rotation is nearly solid and observed modes are pure p modes, ν s is the same for all of them. If the rotation changes with the stellar radius, the splittings depend on the cavities probed by the modes and they will be different. The presence of mixed modes in sub-giant and red giant stars gives opportunities to measure their core rotation. Equation 7 induces symmetrical multiplets. However asymmetrical splittings may be observed if differential rotation in latitude is strong enough (e.g. Gizon and Solanki, 2004), or if the star is oblate due to fast rotation, or in mixed modes due to near-degeneracy effects (Deheuvels et al., 2017).
It is important to note that magnetic fields can also produce asymmetries in the multiplets in both amplitudes and frequencies (e.g. Gough and Thompson, 1990;Goode and Thompson, 1992;Shibahashi and Takata, 1993;Kiefer and Roth, 2018;Augustson and Mathis, 2018).

Mode visibility
We focus in this review on low-degree modes, because they are the only ones that we observe in stars without any spatial resolution due to cancellation effects. Indeed, since all of the information (intensity or velocity) is integrated over the visible stellar disc, the contribution of small-scale modes, i.e. high-degree modes, vanishes. We can only measure high-degree modes in the Sun because its surface is resolved. If we denote the fluctuation induced by a mode f n, ,m (θ, φ) = AY m (θ, φ), then we may show (Dziembowski, 1977;Toutain and Gouttebroze, 1993;Gizon and Solanki, 2003) that the observed amplitude can be expressed as: a n, ,m = r ,m (i)V A .
V is called the mode visibility, r ,m (i) is the relative amplitude of mode inside a multiplet and depends only on the inclination angle i between the rotation axis and the line of sight. Mode visibilities may be written as: where P is the -th order Legendre polynomial, W (µ) a weighting function depending only on the distance to the limb µ. It indicates the contribution of each surface element over the stellar disc. For intensity measurements W (µ) is mainly the limb darkening function L(µ). Limb darkening depends on atmosphere properties (effective temperature, surface gravity, metallicity,...) and the observed wavelength. For velocity observations, since the motion induced by low-degree p modes are mainly radial at the surface, we may approximate the weighting function as W (µ) ≈ µL(µ).
Mode visibilities are plotted in Fig. 10 for low degree. We see that they dramatically drop for > 2.
In practice, we observe mainly modes with = 0, 1 and 2 and some = 3. V mainly depends on the atmospheric properties and the mode physics in the upper layers. We notice that = 3 modes are only visible thanks to limb-darkening effects. Specific derivations of V for CoRoT and Kepler observations have been carried out by Michel et al. (2009) and Ballot et al. (2011a). When V depends on the stellar atmosphere and the instrument, the factor r ,m (i) is purely geometric and reads where P m is the associated Legendre functions. Figure 11 shows the variation of these coefficients with i for = 1 and 2. The squared factor r 2 ,m (i) represents the relative power of modes in a multiplet. We notice that m r 2 ,m (i) = 1. To be able to separate r ,m and V factors in Eq. 8, W must depend on µ only. It is, for example, well known that the amplitude ratio in = 2 and 3 multiplets in solar data observed, for instance, by GOLF does not follow Eq. 10 since the spatial response of the instrument does not depend on µ only. Recent detailed measurements have been presented and discussed in Salabert et al. (2011a).

Frequency separations
The regular structure of p-mode spectra predicted by asymptotic theory is well observed as we can see in Fig. 12 for the solar case. A sequence of modes of degrees 2, 0 and 3, 1 is repeated. We then may define two important seismic variables: the large-and the small-frequency separations or simply large and small separations.
The large separation of low-degree p modes is given by (see also Fig. 12): This regularity is well captured by their first-order asymptotic expansion (Tassoul, 1980): where the asymptotic large spacing ∆ν depends inversely on the sound-travel time between the centre and the surface of the star, which is also called acoustic radius: where R is the stellar radius, and c is the sound speed. The small separation of low-degree p modes is given by (see also Fig. 12): δν , +2 (n) = ν n, − ν n−1, +2 .
(14) Fig. 12 Region of the power spectrum of 16 Cyg A (shown in Fig. 6 between 2130 and 2400 µHz. Some modes =0, 1, 2 and 3 are labelled with radial orders between 19 and 21. The horizontal lines ending with arrows indicate the large spacing, ∆ν, and the small spacings, δ 0,2 and δ 1,3 . The small separation is the difference of two modes with nearly identical eigenfunctions at the surface (e.g. those with almost the same frequency, and thus similar outer turning points) and being only different in the deeper layers, with different inner turning points. Using second-order asymptotic theory (Tassoul, 1990;Vorontsov, 1991) it can be shown that: This asymptotic expression shows that the small separation is dominated by the sound-speed gradient near the core (the integral is weighted by 1/r) and, therefore, it is sensitive to the chemical composition in the central regions. In solar-like stars, generally only δν 0,2 is directly observable. As the frequencies of both modes are very close, they have similar near-surface effects. Hence, the small separation is less affected by such effects than the large separation. However, some residuals can still remain. Therefore, it has been demonstrated that the ratio of the small separation to the large separation, defined as r 0,2 ≡ r 0,2 (n) = δν 0,2 (n)/∆ν (n), can exclude such near-surface effects to a great extent (for more details see Roxburgh and Vorontsov, 2003). By extension, we also define small separations between radial and dipole modes, δ 0,1 and δ 1,0 , as the amount by which the modes =1 (resp. = 0) are offset from the midpoint of the modes =0 (resp. = 1): These quantities are also very sensitive to stellar cores and may be used, for example, to probe for the presence of a convective core (see Sect. 6.7).

Spectral analysis
In this section, we introduce various practical tools that are used to analyze seismic observations of solar-like stars. We do not deal with the computation of a power spectrum density from a velocity curve or a light curve.

Theéchelle diagram
A common way to represent the oscillation spectrum is theéchelle diagram. It consists of plotting the mode frequencies as a function of the frequencies modulo the large frequency separation. Similarly, anéchelle diagram may be built from an observed spectrum by cutting the spectrum into segments of multiples of the large-frequency spacing, stacking them one on top of the next, and making then a 2D map. Doing so, modes with the same degree are aligned along almost vertical ridges. It was first used in helioseismology by Grec et al. (1980) to identify the modes in the solar oscillation spectrum observed from the South Pole. It is now commonly used in asteroseismology to correctly identify the degree and the order of the modes. Any departure from the first-order asymptotic relation will produce curvature and/or wiggles in the vertical ridges. For example, the presence of mixed modes in evolved solar-like stars is clearly visible and easy to identify: some dipole modes are displaced, or bumped, from their original position due to the presence of a mixed mode in the vicinity. Figure 13 shows theéchelle diagram of three solar-like stars observed by Kepler.
The ridges corresponding to the even ( = 0 and 2) modes are clearly shown on the left-hand side of the diagrams, while the ridge of the =1 is visible onto the right. ( =3 mode amplitudes are too small to be unambiguously observed). From left to right, stars are increasingly evolved. Indeed theéchelle digram of KIC 11026764 (right-hand panel) shows a bumped =1 mode (a mode that has been displaced from its original position by the presence of a mixed mode or by the mode immediately below) at ∼ 900 µHz. This is a clear signature of a more evolved star, probably a sub-giant star.

Modelled spectrum
The observed power spectrum density is the superimposition of several components (see Fig. 5) composing a background on top of which the oscillation modes are visible. We denote B(ν) the background and P (ν) the p-mode spectrum. The full model is thus This section provides models for B and P that are commonly used in the literature.

Background model
At high frequency, the spectrum is dominated by the photon noise W ; this is a white noise, independent of the frequency. At low frequency, the spectrum is dominated by long-term trends, including both instrumental (drifts, thermal variations, etc.) and stellar effects (especially activity and rotation). Between these two extreme ranges, the background originates from the surface convection, which is dominated by its granular scales. Low-frequency trends and convective components may be described by empirical laws -initially proposed by ) -corresponding to exponentially decaying temporal variations. The background is then modelled as where the different components are Here, σ i is the rms amplitude of the component, τ i its characteristic time scale, ζ i = 2α i sin(π/α i ) is a normalisation constant. The exponent α i measures the amount of memory in the physical process responsible for the component. A larger exponent means less memory in the process. Exponential decay gives α = 2. In this case, ζ i = 4. Sometimes, following Mathur et al. (2011b), authors define a new timescale, τ eff , as the e-folding time of the autocorrelation function of the signal; for α = 2 both timescales are identical (τ = τ eff ). In practical cases, one granulation component may be enough to model the convective contribution, but a second one, possibly due to faculae, is sometimes required (Karoff et al., 2013a).
This model B(ν) is the background limit spectrum that would be obtained after an infinite observing time, which would average all statistical fluctuations. The observed background is then this limiting spectrum multiplied by a random noise following a two-degree-of-freedom (2-dof) χ 2 statistic. Random processes in time series tend to produce normal (Gaussian) noises in the Fourier domain, both for the real and imaginary parts of the Fourier transform, due to the central limit theorem. Thus, a power spectrum being the sum of squared real and imaginary parts follows a 2-dof χ 2 statistic. It is true for the background but also for stochastically excited modes (see next section). This statistical distribution has been well verified in observations.

Mode model
The last components of the spectrum are the modes themselves. As previously mentioned, solar-like oscillations are not unstable modes but are excited stochastically and damped by turbulence in the outer layers of the convection zone (e.g. Goldreich and Keeley, 1977;Goldreich and Kumar, 1988). Following these authors, each mode can be simply modelled as a randomly excited and damped harmonic oscillator following the equation where ξ(t) is the displacement of the oscillator, η its damping rate, ω o = 2πν o the frequency of the undamped oscillator and f (t) the random forcing function. Assuming that the mean value of the square of the Fourier transform of f (t), |F (ν)| 2 , is a slowly-varying function of ν, and that η ν o , then the power spectrum density of a mode is modelled as a Lorentzian profile multiplied by a stochastic noise following a 2-dof χ 2 statistic. More specifically, the limit Lorentzian profile reads where H is the mode height and Γ is the width at half-height. They are expressed as and For a single mode, the integrated power, P , is given by: which corresponds to the mean square of the mode amplitude (see also the discussion in Jiménez-Reyes et al., 2003). The total energy, E, is taken to be the sum of both the kinetic and the potential energy: where I is the corresponding mode inertia defined in Sect. 4.1. The rate at which energy is dissipated in the modes (e.g. Houdek et al., 1999) can be derived by using the harmonic damped oscillator analogy (e.g. Chaplin et al., 2000): From the study of these equations, the linewidth provides a direct measurement of the damping rate. The mode power (or mode energy) provides a measure of the balance between the excitation and the damping of the modes. Finally, the energy supply rate provides information about the excitation or the forcing of the oscillator. Such a Lorentzian model is sufficient to reproduce the observed modes, even if asymmetries were reported for the Sun since Duvall et al. (1993). These asymmetries are small enough for lowdegree modes -typically of the order of a few percent (e.g. Toutain et al., 1998;Chaplin et al., 1999b;Thiery et al., 2000) -to consider such Lorentzian description as sufficiently accurate for analysing observations shorter than a year. For longer time series, including asymmetries may be relevant (Benomar et al., 2018).
This Lorentzian model is also valid as long as the modes are resolved, i.e. the spectra resolution is finer than the mode width. It means that the observing time is longer than the mode lifetime. This point is verified for observed p modes in solar-like stars, but it may be invalidated at very low frequency.
The oscillation spectrum is the sum of a full sequence of modes: where M is the profile of a multiplet where i is the inclination angle and ν s the rotation splitting. The coefficients a ,m are given by Eq. 10 (see. Sect. 4). To obtain this expression, we assume that the intrinsic mode heights and widths inside a multiplet are the same. Since the components have similar frequencies (ν s ν o ), the energy injected by stochastic excitation is the same (see Eq. 23), and we also assume that damping processes only depend smoothly on frequency.

Fitting spectra
This model S(ν; p) described in the previous section is parametrized and we aim in determining the most likely values for these parameters p (e.g. mode frequencies, widths, heights, splittings...), given an observed spectrum. The observed spectrum Y = {Y i } is sampled at n frequencies ν i and we assume that the realization noise for each point is independent. This last point is absolutely verified when the spectrum is estimated through a Fast Fourier Transform of an evenly spaced time series without gaps. Since the observed power density is distributed around the limit spectrum with 2-dof χ 2 statistics (e.g. Woodard, 1984;Appourchaux et al., 1998), the likelihood is Given a model S, one seeks to find, for a spectrum Y observed over a chosen frequency interval, the parameters p that maximize L(Y ; p). The way to choose the frequency interval to consider and the details of the models depend on the fitting strategy and will be discussed in 5.5. Maximum likelihood estimators (MLE) are frequently used to analyse seismic data of solar-like stars. This is an inheritance of helioseismology where this approach was intensively used . In practice, the negative logarithm of the likelihood function − ln L is minimized with a standard algorithm such as a modified Newton's method.

Errors and correlations
To estimate the covariance matrix of parameters p, a typical method is to approximate it by the inverse of the Hessian matrix. The uncertainties on fitted parameters are therefore taken as the square roots of the diagonal elements of the inverted matrix. These estimates are based on the Kramer-Rao theorem. By using them, we have to keep in mind a few crucial points: (i) these error estimates are only lower limits of the statistical errors and (ii) it is only asymptotically valid: the statistical distribution of parameters are assumed to be normal, that is not necessarily the case (e.g. Ballot, 2010). Thus, the variables are to be carefully chosen, for example h = ln H and γ = ln Γ are more suitable variables than H and Γ to estimate the errors through the Hessian (see discussion in Toutain and Appourchaux, 1994). The prefered way to estimate the errors remains in performing Monte Carlo simulations. Nevertheless, the Kramer-Rao theorem may be used to understand the evolution of errors with various parameters. For example, following Duvall (1990), Libbrecht (1992) showed that the lower limit for the error of the frequency of a radial mode is where T is the observation time, β = B/H is the local background to mode height ratio, and More general expressions have been proposed for multiplets in the solar case (i ≈ 90 • ) by Toutain and Appourchaux (1994) and for any inclination angle by Ballot et al. (2008). Such relations show that for resolved modes, the precision of frequencies varies with T −1/2 : to reduce the error by a factor of 2, we need observations 4 time longer. It is also important to notice that some parameter estimates are strongly correlated. It is specially the case between the mode height H and width Γ , which are strongly anticorrelated, making the quantity HΓ , hence the mode energy P , better determined than H and Γ individually (e.g. Toutain and Appourchaux, 1994). Similar correlations exist between the inclination angle i and the splitting ν s when ν s is not significantly larger than the mode width. The correlation is well visible in a likelihood mapped in the (ν s , i) plane. The likelihood maximum show a banana shape structure (Fig. 14). In this case, the projected splitting ν s sin i is better determined than the splitting itself, even when the inclination is poorly measured (see Ballot et al., 2006Ballot et al., , 2008.

Bayesian inference
During the last decade, "Frequentist" approaches based on MLE have been completed with Bayesian approaches to analyse oscillation spectra. First introduced for the seismology of solar-like stars to correctly interpret the first CoRoT data (Benomar et al., 2009a;Gaulme et al., 2009;Gruberbauer et al., 2009), it has been extensively developed and used ever since (e.g. Handberg and Campante, 2011;Gruberbauer et al., 2013;Corsaro et al., 2013;Davies et al., 2016;Lund et al., 2017). The likelihood is nothing but the probability to observe the data set Y assuming a parameter set p, under given a priori information I (including for example a spectral model S), we denote it p(Y|p, I)(≡ L(Y; p)). We maximize this probability with MLE. However, we would prefer to assign a probability to a parameter set p for the observation Y, i.e. p(p|Y, I). The is called the posterior probability. It is linked to the likelihood through Bayes' theorem: p(p|I) is the prior probability, i.e. the probability of getting these parameters before looking in the data; it includes our current knowledge (physical properties, information coming from other data or other parts of the spectrum) and current ignorance. The probability p(Y|I) is called the global likelihood: p(Y|I) = p(p|I)p(Y|p, I) dp.
Using a marginalization procedure, we can derived the marginal posterior probability distribution for a subset of parameters of interest p I by integrating out the remaining parameters p N (p = {p I , p N }), called nuisance parameters: Using Bayesian methods, we then obtain the statistics for all of the parameters of our model: not only an estimated mean and variance, but the full distribution. For spectra with low signal to noise ratio, Gaulme et al. (2009) showed that MLE may be biased and Bayesian approaches are more robust.

Priors: knowledge and ignorance
Bayesian methods sample the function p(p|Y, I) to give a global picture of the problem. Moreover, Bayesian approaches allow us to include relevant priors in fitting procedures. Even more, priors are needed: posterior probabilities only make sense when prior probabilities are set. Sometimes, when our knowledge is limited we look for priors that take that ignorance into account. As an example, let us consider the prior for the inclination angle i. We may certainly consider our prior independent of the other stellar quantities. When we do not have any complementary observations (or we do not want to use them), we would naturally assumed that all orientations are evenly probable. This does not mean that prior probability is uniform between 0 • and 90 • . If we assume isotropy, the prior probability distribution for i is p(i) di = sin i di. A uniform distribution for i would favour a rotation axis oriented toward us. Ignorance priors for frequency (or splitting) are uniform probability distributions whereas it is uniform in logarithm for heights and widths (see, e.g., Benomar et al., 2009a;Handberg and Campante, 2011, for more detailed discussions).

Markov Chain Monte Carlo
The most common method to sample the posterior probability is to use a Markov Chain Monte Carlo (MCMC) approach. Even if they are rather slow, their implementation is easy and flexible. The Metropolis-Hastings algorithm (Metropolis et al., 1953;Hastings, 1970) is the most commonly used method. Its application to solar-like star oscillation spectra are detailed in Benomar et al. (2009a) and Handberg and Campante (2011). The distribution that we want to sample, called the target distribution, is p(p|Y, I). To do so, we construct a pseudo-random walk in the parameter space such that the density of drawn points in a given region of the parameter space is proportional to p(p|Y, I). A Markov Chain achieves such a random walk. The chain is built such that a new point p t+1 is added to the chain depending on the previous point p t according to a time-independent transition kernel p(p t+1 |p t ). In the Metropolis-Hastings algorithm, a new point p new is randomly drawn from a proposed probability distribution centred on p t , p p (p new |p t ) (In practice we consider multivariate normal distributions). We then accept p new as a new point (p t+1 = p new ) with an acceptance probability If we reject the new point, then p t+1 = p t . The expression for α(p new , p t ) is simpler than the ones found in literature, because we have taken into account the symmetry of normal distributions that we used as proposed distributions (p p (p 1 |p 2 ) = p p (p 2 |p 1 )).
From a mathematical point of view, the Markov chain asymptotically represent the target function (i.e. p(p|Y, I)) independently from the choice of the proposal probability distribution p p . However, since we will get a chain with a limited size, the p p law must be chosen to ensure an efficient sampling of p(p|Y, I) in reasonable computing time. As we consider p p as a multivariate normal law, we must find an suited covariance matrix for this law. Different automated algorithms are used to constructed such covariance matrices, for example in Benomar et al. (2009a) and Handberg and Campante (2011).
Finally, improved versions of MCMC using parallel tempering (Earl and Deem, 2005) are often implemented in asteroseismology. It is very similar, except that we do not use only one chain, but several ones corresponding to different "temperature" (by analogy with statistical physics), characterised by a parameter β. Each chain samples the following probability: for β = 1 (the "cool" distribution) we recover the usual posterior probability. By decreasing β (i.e., by using "hotter and hotter"), the function is flatter and flatter, and is reduced to the prior distribution for β = 0. Doing so, we understand that it is easier to explore regions of a hotter distribution that would never been visited for a cooler one. The different chains are independent, but randomly we swap the position of the two chains. It is very useful, for example, when the posterior distribution possesses several well separated local maximums. With a classical MCMC, the chain may get stuck in one of them and it may take a very long time to sample the different maxima. Parallel tempering is an efficient way to avoid this (e.g. Handberg and Campante, 2011). MCMC is not the only sampling method developed for seismology of solar-like stars. For example, Corsaro and De Ridder (2014) proposed a spectrum analysis method based on a nested sampling Monte Carlo algorithm.

Fitting strategy: local vs. global fits
Several strategies can be used to fit an oscillation spectrum. Historically, to analyse helioseismic data of the Sun seen as a star, a local strategy was developed: modes are fitted within windows narrow enough to isolate a single mode or a pair of modes = 0, 2 or = 1, 3. For intensity measurements, = 3 modes are generally small enough to be ignored, and = 1 modes are often fitted alone. A second approach is to perform a global fit of all modes above a given amplitude threshold around the maximum of the p-mode hump, simultaneously. This approach, pioneered by Roca Cortés et al. (1999) on solar data, became more popular in asteroseismology when first CoRoT data had to be analysed .

Fig. 15
Example of fits (green line) =2, 0 and =3, 1 modes, left-and right-hand panels respectively at low (top) and high (bottom) frequencies of a typical GOLF spectrum. In the first case, the lifetimes of the modes are longer and therefore the linewidths are smaller than for the modes at high frequency.
Local fitting can be performed when mode pairs are well separated from each other, i.e. when the large separation is large enough to ensure that the wings of the Lorentzian profiles of other modes do not contaminate the fitting window. In a local approach, we assume the background is flat and is only parametrised with a constant. We use a single linewidth for the two modes ( = 0 and 2), and we also may use a single free height, H 0 , which is the height of the = 0 mode, and impose H 2 = V 2 2 /V 2 0 H 0 for the = 2 mode by fixing the visibilities from theory. This last constraint is useful when = 2 and 0 modes are blended. Even if this local fitting scheme is easier to implement, with a reduced number of free parameters for each fitting procedure, it has been proven that a global approach is generally better suited ). An example of such global fit is presented in Fig. 16 corresponding to the analysis of 11 orders of the CoRoT target HD 169392 . A main advantage of this approach is to impose that some parameters be the same for all modes, especially the inclination angle i, and for some cases the splitting ν s . Due to the correlation previously mentioned between i and ν s , i may be poorly determined for each mode independently, we thus increase its precision by using a common value for all modes. Since all parameters are correlated, a better determination of i improve the determination of ν s , hence linewidths, hence heights... However, by doing this, we assume that all of the stellar layers rotate with the same axis. This reasonable assumption may be questioned for the deep core (e.g. Bai and Sturrock, 1993, for the Sun). Free parameters to parametrize widths and heights can also be reduced assuming that they vary only slowly with frequency; one free parameter per large separation may be enough. Mode visibilities can also be fixed or fitted as global parameters. Concerning the background, it cannot be taken as constant in the fitting window, a more complete profile (Sect. 5.2.1) must be used. Its parameters may be fixed by a previous overall fitting or left (partially) free (see, e.g., Appourchaux et al., 2008;Ballot et al., 2011b;Mathur et al., 2011a).
A global fit of the background is generally previously performed by ignoring the frequency range of p modes, or modelling the p-mode hump as a Gaussian profile or two Lorentzian profiles (as has been shown for the solar case by Lefebvre et al., 2008). An example of such a fit for the two solar analogs 16 Cyg A and B observed by Kepler is shown in Fig. 17  .

Hypothesis tests and model comparisons
The fitting techniques described in the previous section provide the best parameters of a model, assuming that the hypotheses are correct. Hypotheses include especially the choice of the model, and some a priori information such as the mode identification. It is however frequently necessary to test several competiting hypotheses to make a final decision. The most outstanding one is to test the mode identification. The first solar-like star observed by CoRoT was the F-type star HD49933. Early-type solar-like stars have broad modes due to strong mode damping, thus mode widths are not smaller than the small separation, making the even modes = 0, 2 overlap. As a consequence, identifying the = 1 ridge from the = 0, 2 ridge is not obvious in anéchelle diagram. Both Index Fig. 17 Background fitting of the two solar analogs 16 Cyg A (left) and B (right) observed by Kepler. The PSD has been smoothed by a 20 µHz boxcar (grey), with best-fitting background components attributed to granulation (dashed lines), stellar activity and/or larger scales of granulation (dot-dashed lines) and shot noise (dotted lines), with the sum of the background components plotted as solid black lines .
hypotheses have to be tested with a frequentist  or Bayesian approach (Benomar et al., 2009a). Nevertheless to complement statistical approaches, White et al. (2012) propose a method based on physical properties of the stars to disentangle both scenarios using the variation of (see Eq. 12) with effective temperature.
It is also very useful to use model comparisons to validate the significance of a splitting. Thus, we can verify if a mode is significantly better fitted with a multiplet than with a single Lorentzian profile (e.g. Deheuvels et al., 2015).
When the two hypotheses H 0 and H 1 have the same number of free parameters (typically when we want to compare the two different mode identifications), the ratio of their likelihoods L 0 and L 1 gives a direct comparison of the two hypotheses. The p-value of favouring H 0 over H 1 is then p = (1 + L 1 /L 0 ) −1 . However, as already mentioned, since MLE may be biased for low SNR, such a comparison may be skewed.
Assuming now that the two hypotheses H 0 and H 1 have respectively n 0 and n 1 free parameters (n 1 > n 0 ), we want to determine whether H 1 may be favoured. We thus need to assess the likelihood improvement of using more free parameters. Wilks (1938) showed that the quantity Λ = 2(ln L 1 − ln L 0 ) follows a χ 2 distribution with n 1 − n 0 degrees of freedom. The p-value is thus the probability A low p-value rejects the simplest model in favour of the most complex. Using a Bayesian approach to compare two different competing models M i and M j (including different set of priors), we need to compute the evidence of the model We then introduce the odds ratio By assuming that the two models are equally probable (p(M i |I) = p(M j |I) = 1/2), odds ratios reduces to Bayes factor The relative probability to favour model M j over model M i is p = (1+O ij ) −1 . It may be generalized to more competing models. The main difficulty is to compute the evidences p(Y|M i , I). Parallel tempering MCMC allows us to compute them using the different chains with different temperatures. It reads (e.g. Gregory, 2005) p where the term · · · is the average of the logarithm of the likelihood for all points of a given tempered chain, characterised by β. This shows another interest to perform parallel tempering MCMC.

Global seismic parameters
We mainly detailed in the previous sections how individual mode properties can be extracted. Of course, global parameters, especially the mean large separation ∆ν, the frequency at maximum amplitude ν max and the maximum amplitude of radial mode A max can be derived from a detailed fit. However, there exist quick methods to recover these main features. Such methods are useful when we must deal with a large number of stars. This is the reason why they have been massively used to analyse several thousand red giants observed by Kepler, but they have also been applied to main-sequence stars. Various pipelines have been developed during the last decade by several teams around the world (e.g. Roxburgh, 2009b;Hekker et al., 2010;Huber et al., 2010;Kallinger et al., 2010;Mathur et al., 2010b;Mosser et al., 2011b). The most common technique to measure ν max is to fit a bell-shape curve to the p-mode hump and define its maximum as ν max . A max can be derived from the total mode power around the maximum. By measuring the total power of the modes over a large separation interval around ν max , we measure the total power of one mode of each degree. Using mode visibilities and assuming that all modes in the interval have the same intrinsic amplitude, we recover the power of a radial mode. The square root of this quantity provides the rms amplitude (see Kjeldsen et al., 2008a). Generally, we convert the observed amplitude into a bolometric amplitude using bolometric corrections derived for each instrument (see, e.g., for CoRoT and Kepler Michel et al., 2009;Ballot et al., 2011a).
To determine ∆ν, two techniques may be used to recover the regular pattern of p modes. This first one is to find a maximum of the autocorrelation of the spectrum in the p-mode region. The autocorrelation lag providing the largest peak is the large separation, another peak (generally slightly smaller) occurs at ∆ν/2 when the = 1 modes coincide with the = 0, 2 modes. The second technique, which has appeared to be more robust in practice, is to consider the Fourier transform of the power spectrum in the p-mode region. Doing so, the largest peak of the Fourier transform is τ = 2/∆ν, which corresponds to the main regularity visible in spectra produced by the alternation of even and odd modes. Since products in Fourier domain are convolutions in the time domain, we can easily show that this technique is equivalent to looking at the autocorrelation of the time series. We can refer to Roxburgh and Vorontsov (2006) for detailed discussions of its use and Mosser and Appourchaux (2009) for discussions on measurement errors with this technique.
6 Inferences on stellar structure

Scaling relation for masses and radii
If detailed modelling of a star to reproduce the observed frequencies is the most accurate way to determine its mass and radius, a simpler approach using scaling relations from solar values have been massively used these last few years. These scaling relations link the mass M and radius R of a star to the large separation ∆ν, the frequency at maximum amplitude ν max , and the effective temperature T eff .
A scaling relation for ∆ν can easily be derived assuming an homology relation between stellar structures . Assuming two homologous stars (with radii R and R and masses M and M ) means that m(r)/M = m (r )/M for all r/R = r /R , where m(r) is the mass of the star encompassed within the radius r. In homologous stars, Kippenhahn and Weigert (1990, §20.1) show that the density and pressure profiles from one star is deduced from the other through the following scaling relations: ρ (r ) = (M /M )(R /R) −3 ρ(r) and p (r ) = (M /M ) 2 (R /R) −4 p(r). Hence we deduce for the sound speed profiles that c (r ) = (M /M ) 1/2 (R /R) −1/2 c(r). Using Eq. 13, we then show that The large separation then scales with the square root of the mean stellar density. Of course, real stars are not homologous, however, it is well verified with models (e.g. Ulrich, 1986). White et al. (2011b) show that the agreement is better than 5%. Fixing the Sun as a reference, any star can be scaled from it as : where ∆ν =135.1 ± 0.1 µHz, as derived by Huber et al. (2011) using 111 subseries of 30-day each collected by the Variability of solar IRradiance and Gravity Oscillations (VIRGO) instrument (Fröhlich et al., 1995) aboard the Solar and Heliospheric Observatory spacecraft (SoHO, Domingo et al., 1995) spanning from 1996 to 2005 and analyzed them in the same way as asteroseismic data. Nevertheless, the choice of the reference values is not that straight forward and may be the subject of debate. One thing is clear, each pipeline producing ∆ν and ν max to be used in these scaling relations has their own solar reference values computing strictly with the same procedure as the stellar values. The second quantity scaling with stellar global parameter is ν max . Brown et al. (1991) suggested that this quantity scales as the acoustic cutoff frequency ν c (see Sect. 4.1.1). This assumption has been justified by Belkacem et al. (2011). Using Eq. 4 and noting that in the stellar atmosphere, approximated by an isothermal atmosphere of temperature T eff , H p ∝ T eff /g and c 2 ∝ T eff , we deduce that ν c ∝ g/ √ T eff . Hence, ν max can be scaled from the solar value as follows : where T eff, = 5770 K, and ν max, = 3090 ± 30 µHz . Scaling relations from solar values have been extensively tested, in particular, by comparing for example the expected and observed values of ν max for some well-studied stars with accurate parallaxes . Figure 18 shows the comparison of the observed and predicted values of ν max for 14 stars. The general agreement is very good excepting the two lowmass stars τ Cet and 70 Oph A. From a theoretical point of view, some predictions of lighter 0.7 and 0.9 M models show a double bump in the p-mode hump, which complicates the correct extraction of this parameter (e.g. Chaplin et al., 2008b). However, one of the maxima usually lies close to the ν = ν max line. Moreover, ground-based observations of Procyon (Arentoft et al., 2008), as well as other CoRoT (e.g. Mathur et al., 2010b) and Kepler stars have shown that more massive stars could also present a double bump in the p-mode envelope. The question of the definition for ν max arises in this case. Combining Eqs. 43 and 44, we deduced that Because of the relatively weak change in T eff among solar-like stars and the weak dependency on M , we recover that solar-like oscillations in main-sequence stars to a good approximation follow:  Figure 19 (top) shows the relation ∆ν versus ν max computed Index using 1700 stars from the main sequence (black diamonds) completed with red clump (red triangles appearing in a diagonal branch between 20 to 50 µHz in the bottom panel) observed by the Kepler mission. Although the relation appears to be constant for all of these stars, several authors Huber et al., 2010) have suggested that the slope is different for red-giant and mainsequence stars. To enhance such a difference, we have subtracted the luminosity dependence by raising ν max to the power of b = 0.75. A fit to the residuals below and above ν max = 300 µHz -which roughly marks the transition from low-luminosity red giants to sub-giants -enhances a steeper slope as ν max increases. It is important to note that, for ν max close to the solar value for example, the use of a power-law relation calibrated to red-giant stars would lead to an underestimation in ∆ν by ≈ 10%. Recent comparisons of different global seismic pipelines (mostly for red giants) can be found in Serenelli et al. (2017) and Pinsonneault et al. (2018). In a similar way, it is possible to derive a scaling relation from the solar values for the maximum oscillation full -not RMS -intensity amplitude of radial modes, A max , of the radial modes (following, e.g., Samadi et al., 2007): where the exponent s usually has values between 0.7 and 1.0. The β-coefficient has been introduced by Chaplin et al. (2011a) because without any further correction the above relation is known to overestimate the amplitudes for the hottest solar-like stars (Houdek, 2006). A max, is 4.4 ± 0.3 ppm for the VIRGO green channel (λ = 500nm) as calculated by Huber et al. (2011). The inferred bolometric amplitude of A max, ,bol = 3.5 ± 0.2 ppm is in good agreement with the value of 3.6 ppm (Michel et al., 2009) commonly adopted in asteroseismology.
6.2 Model-independent determination of masses and radii Combining Eq. 43 and 44, we derive masses and radii of stars as a function of the seismic variables ν max and ∆ ν , assuming that T eff is known from photometry or spectroscopy: These scalings provide a quick way to measure mass and radius without any models. The errors on these measurements will depend first on the precision of ∆ν and ν max determination. For Kepler data, global methods such as the A2Z or the Octave pipelines (Mathur et al., 2010b;Hekker et al., 2010) provide precisions around 10% and 20% for masses and radii according to their authors. However, using individual frequencies to compute averaged ∆ν and ν max gives uncertainties of about 3 % in radius and 9 % in mass (Mathur et al., 2012). Nevertheless, these relations are only approximations and to test their accuracy we can compare the values inferred from the scaling relations with the best modeling performed on 22 Kepler targets for which we have been able to determine the individual frequencies of a few dozens modes (Mathur et al., 2012). The modeling was done using the Asteroseismic Modeling Portal (AMP), which uses a parallel genetic algorithm (Metcalfe and Charbonneau, 2003) to optimize the match between the model output and the observational constraints. In Fig. 20 we have compared the estimation of both parameters. Stellar radii from both methods are in a very good agreement, while we have a higher scatter in the comparison of the masses, as was expected because the exponents involved on the seismic parameters were 2 to 3 times bigger. According to this work, the scaling relations tend to overestimate the radius by +0.3σ and the mass by +0.4σ relative to the values from AMP (σ 2 being the quadratic sum of the uncertainties from the two methods). This corresponds to radii overestimated by about 1% and masses by about 4% typically. These results suggest that observations of the global oscillation properties combined with an effective temperature through the scaling relations from the solar values can provide reliable estimates of the radius and mass but with a lower precision than AMP when individual frequencies are available. Indeed, the precision reached in that latter case by AMP, in the mass and radius, is much better: 0.8 % in radius and 1.1% in the mass.
It has also been possible to test the global seismic scaling relations using independent measurements of radii for main-sequence and giant stars using interferometry. Huber et al. (2012) performed such an analysis comparing asteroseismic radii obtained with Kepler with interferometric radii obtained with the CHARA array and found an agreement of around 4% (see also White et al., 2013). Unfortunately, because main-sequence stars have small angular diameters, it is extremely difficult to properly extract their radius, these stars being more prone to systematic errors in the adopted calibrator diameters than sub-giants and red giants. This is the reason why it has been recommended to restrict any comparison of asteroseismic diameters with interferometry to stars with angular diameters larger than > 0.3 mas excluding many main-sequence solar-like stars . Another independent validation of the asteroseismic radius can be done using astrometric results (e.g. Silva Aguirre et al., 2012) from Hipparcos (van Leeuwen, 2007) or Gaia (Perryman et al., 2001). Indeed, Huber et al. (2017) compared asteroseismic radii of 2200 oscillating stars observed by Kepler (including 440 main-sequence stars and sub-giants from Chaplin et al., 2011b) with Gaia DR1 results included in the Tycho-Gaia Astrometric Solution (TGAS, Michalik et al., 2015). The overall agreement found was excellent, which helped to empirically demonstrate that asteroseismic radii computed using global seismic scaling relations were accurate to ≈ 10% for stars ranging from ≈0.8 to 10 R without any visible offset between the two radii determinations for main-sequence stars (1 to 1.5 R ). Moreover, no significant trends were found with metallicity [Fe/H] = -0.8 to +0.4 dex.

Model dependent determination of masses and radii
The ultimate precision of asteroseismology can only be reached when individual mode frequencies are combined with spectroscopy and astrometric observations to infer the best stellar model of each star. Tens of individual mode frequencies of the best 66 main-sequence and sub-giant solarlike stars seismically characterized by Lund et al. (2017) were distributed among seven different modeling teams to determine radii, masses, and ages for all the stars in this sample (Silva Aguirre et al., 2017). As it is shown in Fig. 21, masses and radii agree within the error bars for most of the cases with average uncertainties better than ∼ 2% and ∼ 4% in radius and mass respectively. The differences found between the different results could be explained in terms of the different physics used in the codes and the way in which the uncertainties were computed during the minimization process (that explains the large uncertainties found in radius computed by Cesam2k Stellar Model Optimization -C2kSMO, Lebreton and Goupil (2014)). For further details see Silva Aguirre et al.

Ensemble observational asteroseismology
We can take advantage of the ensemble analysis of several hundreds MS solar-like pulsating stars to reveal general trends on the properties of stellar pulsations. In Fig. 22 the ensemble of solar-like stars with measured pulsations is shown in a seismic HR diagram, i.e., plotting the large frequency separation as a function of the effective temperature. In the top panel are represented groundbased Doppler-velocity observations (black triangles) (Bouchy and Carrier, 2002;Kjeldsen et al., 2003;Martic et al., 2004;Bouchy et al., 2005;Kjeldsen et al., 2005;Bedding et al., 2006;Carrier and Eggenberger, 2006;Kjeldsen et al., 2008a;Mosser et al., 2008;Bonanno et al., 2008;Teixeira et al., 2009;Bedding et al., 2010b;Bazot et al., 2011, and  Using different solar-like stars located along the evolutionary tracks traced in Fig. 22, it is today possible to build observational evolutionary sequences of dozens of stars of similar masses (and whenever possible, similar metallicities), using scaling relation discussed in previous sections. Therefore, ensemble asteroseismology studies open new ways of performing differential analysis on field stars along evolutionary sequences and thus constraining how the internal properties of stars change with mass, metallicity, and evolution (e.g. Silva Aguirre et al., 2011b; Ozel et al., 2013). An example of an observational evolutionary sequence of solar analogues is shown in Fig. 23.
Mode amplitudes, heights, and linewidths are more difficult to measure.  showed that systematic effects between 8 different groups of "fitters" were mainly due to the way that the convective background was treated, as well as on the fitted values of the rotational splittings and inclination angles. Following a correction scheme based on the one-fit approach of Toutain et al. (2005), they demonstrated that the systematic effects could be reduced to less than ± 15 % for the linewidths and heights, and to less than ± 5 % for the amplitudes. It is worth mentioning that different convective background models introduce frequency-dependent systematic errors that could bias any comparison with theoretical predictions and between different groups of fitters. Once frequencies, amplitudes, and linewidths are measured with precision and accuracy, it is possible to look for global trends with, for example, age, effective temperature, or mass.
Linewidths also depend strongly on effective temperature. They increase with temperature as is shown in Fig. 25. Several scaling relations of the linewidths at ν max with the effective temperature have been proposed in the literature Corsaro et al., 2012;Appourchaux et al., 2012Appourchaux et al., , 2016. These relations are qualitatively in good agreement with theoretical predictions, although it would be necessary to introduce other dependencies such as with log g (Belkacem et al., 2012). The main differences arises at low temperatures where it is necessary to introduce sub-giants and giants to improve the fit. Moreover, to go further in the comparison with the theoretical predictions, it would be necessary to improve the precision in the calculation of the effective temperatures whose large errors limit the final precision of the fitted scaling relation . Fig. 22 Seismic Hertzsprung-Russell diagram showing ∆ν as a function of T eff of main-sequence, sub-giants, and some early red giant stars. In the top panel a selection of stars is represented while in the bottom panel stars observed in short-cadence by Kepler and K2 have been added. Triangles represent stars analyzed individually or in groups of a few stars. Circles represent stars studied as an ensemble. Colors represent observations made from the ground (black) or with spacial missions: green for CoRoT, red or magenta for Kepler and blue for K2. WIRE observations are included as ground-based because these targets were also observed from ground. The Sun is represented with its usual label ( ). Evolutionary tracks, computed with the ASTEC code (Christensen-Dalsgaard, 2008a), are shown for masses ranging between 0.9 and 1.6 M at solar composition (Z = 0.0246). The values used to plot these figures are listed in Table 1.

The C-D diagram
The Christensen-Dalsgaard diagram (or simply C-D) is the representation of the small separation δν 0,2 as a function of the large separation ∆ν (Christensen-Dalsgaard, 1988). By placing stars on this diagram (see Fig. 26, it is possible to directly determine their masses and ages. Indeed the tracks for different masses and ages are well separated. However, it is important to note that this diagram depends on metallicity. Nevertheless, if the metallicity of a star is known, such a diagram is an interesting tool to derive its age. A C-D diagram has been constructed for 76 Kepler stars and several CoRoT targets by White et al. (2011a). Results are shown in Fig. 26. We can see how evolution tracks are well split during the main sequence but converge at the terminal age main sequence. It also gives an estimate of the mass within 4 -7% when the metallicity  is determined within 0.1 dex. The age-determination accuracy depends strongly on the smallseparation determination.

Modelling a star and surface effects
To go beyond global parameters and global diagnostics, we can model a star in detail to match the observed frequencies. Several approaches are then possible. We can used a precomputed grid of models to find the model with the closest frequencies (e.g. Stello et al., 2009b;Quirion et al., 2010;Gai et al., 2011), we can also perform on-the-fly modelling to recover the best model. For this kind of direct approach, we need to minimize a cost function, generally a χ 2 , between observed and modelled frequencies. This cost function may integrate external constraints such as spectroscopic or photometric determinations of temperature, metallicity, and eventually the mass if the star belongs to a multiple system. Nevertheless, the use of frequencies without specific care provides biased results due to surface effects. It is known for the Sun (e.g. Christensen-Dalsgaard and Berthomieu, 1991) that significant discrepancies between observed frequencies and those computed from 1D model occur above ≈ 2000 µHz, as shown in Fig. 27. This is due to the poor modelling of upper layers. Low frequency modes are almost unaffected because their outer turning point are deeper in the star, as discussed in Sect. 4.1.1 (Fig. 8). Upper layers are affected by dynamical processes missing in 1D modelling: 3D model atmospheres show that turbulent pressure cannot be neglected (at the photosphere it reaches about 15% of the total pressure) and the structure is also affected by convective back-warming (e.g. Trampedach et al., 2013). Moreover the β plasma becomes small, letting the magnetic field play a role (it is the reason why its variations during stellar cycles are visible on frequencies, see Sect. 8). There are not only structural effects, but also modal effects: non adiabatic effects must be included as well as fluctuations of the turbulent pressure. To improve the modelling of outer layers, 1D models patched with 3D simulations have been proposed for decades (e.g. Stein and Nordlund, 1991;Rosenthal et al., 1999;Yang and Li, 2007;Piau et al., 2014;Bhattacharya et al., 2015;Sonoi et al., 2015Sonoi et al., , 2017Ball et al., 2016;Magic and Weiss, 2016;Houdek et al., 2017;Trampedach et al., 2017). To compute the oscillations primarily two different assumptions were proposed to model the Lagrangian fluctuations δp t of turbulent pressure. The first one, called the Gas Γ model (GGM), consists in assuming that they vary as the total pressure, i.e. δp t /p t ≈ δp tot /p tot ≈ δp gas /p gas = Γ 1 δρ/ρ. In the second model, called the Reduced Γ model (RGM), the turbulent pressure vanishes δp t /p t ≈ 0, hence δp tot /p tot ≈ (Γ 1 p gas /p tot )δρ/ρ. Rosenthal et al. (1999) shows that GGM gives better results than RGM. Using patched models with GGM improves the computed high frequencies, but a residual of several µHz remains. Recently Sonoi et al. (2017) and Houdek et al. (2017) used time-dependent convection models to compute the oscillations of patched models. The first team used prescriptions from Grigahcène et al. (2005) and Dupret et al. (2005), whereas the second one used a model developed by Gough (1977). Doing so, the discrepancies do not almost depend on frequency and remain limited to a few µHz (Fig. 27).
Whereas it was easy to get rid of surface effects for the Sun because its mass and radius are extremely well known, dealing with surface effects is more difficult for other stars. A solution is to correct the model frequencies ν model with an ad-hoc term δν to get corrected frequencies ν corr = ν model + δν. The first correction that was broadly used was proposed by Kjeldsen et al. (2008b) it takes the form where ν 0 is a reference frequency, b is calibrated from the Sun, and a has to be fitted. This relation has been extensively used during the last 10 years (e.g. Tang and Gai, 2011;Mathur et al., 2012). Of course, it assumes that the mode inertia varies slowly with frequency, that is the case for main-sequence stars. For sub-giant stars, corrections must be divided by the mode inertia I. Other prescriptions have been proposed later. Following a suggestion by , Ball and Gizon (2014) proposed a cubic correction δν = a 3 (ν/ν 0 ) 3 /I where ν 0 = ν c and a composite correction with an additional inverse term is These new prescriptions appear to give better fits than the old prescription for main-sequence stars (Ball and Gizon, 2014) and sub-giants . A recent study performed on a large sample of 66 Kepler stars shows a significant improvement in the fit by using the composite form including the inverse term (Compton et al., 2018). It is also worth mentioning an alternative prescription derived from patched models by Sonoi et al. (2015): Finally, another way to get rid of surface effects is to construct seismic variables where surface effects cancel out. Roxburgh and Vorontsov (2003) proposed to use the ratio of the small separation (δν 02 or δν 01 ) to the large separation to remove a large part of the surface effects (see also Otí Floranes et al., 2005;Roxburgh, 2005). Such a solution has been used, for example, by Silva Aguirre et al. (2013) on two Kepler targets. In such cases, direct modelling (including surface effect corrections) led to several models with frequencies compatible with the observations but significantly different internal structure. Using ratios pins down the number of plausible models, reducing the errors on radius, mass, and age down to 2%, 4%, and 10%. A large sample of 66 solar-like stars observed by Kepler has been homogeneously analysed by Lund et al. (2017) and modelled by Silva  in the LEGACY project. The average uncertainties on the model parameters for this sample are again 2%, 4%, and 10% for radius, mass, and age.

Constraints on the internal structure
Beyond global parameters, asteroseismology allow us to derive constraints on the internal structure that can be compared with stellar models. Acoustic glitches provide interesting diagnosis. An acoustic glitch is a perturbation in the frequencies due to the presence of discontinuities (or rapid variations) in sound speed (or its derivatives) in the stellar structure. Such discontinuities generate partial reflections of propagating waves, resulting in a change in the resonant frequencies of a cavity (see  for a comprehensive toy model). For the Sun and stars, the impact of glitches have been theoretically studied since late 1980s (see, e.g., Vorontsov, 1988;. It slightly modifies mode frequencies with an oscillatory behaviour such as δν n, ∝ sin(4πτ g ν n, + φ g ), where τ g is the acoustic depth of the glitch, i.e. the travel time of a sound wave to reach the glitch from the surface: r g being the radial position of the glitch and φ g a phase shift. For low degree modes, this phase φ g is the same for all modes. The amplitude and the envelope of the perturbation change with the nature and the abruptness of the transition (the shaper, the larger). Sharp variations of the sound speed (or derivative) occur in stars where (i) the adiabatic exponent Γ 1 quickly decreases, which occurs in the ionisation of abundant elements (H, He I, He II) (ii) when the thermal gradient changes, especially when energy transport processes change from radiative to convective (or the opposite). Measuring glitches allow us to probe (i) the helium ionisation zone, possibly giving constraints on the helium abundance in the envelope, (ii) the position and the structure of the base of the convective envelope (BCE). Glitches generated by convective cores are too deep -in terms of acoustic depth -to be seen as oscillations in frequency. We will discuss convective cores later on.
To search for these oscillations we can look directly at the frequency (e.g. Monteiro et al., 1994Monteiro et al., , 2000. It is also possible to look for them in the surface phase shift or its derivatives (e.g. Lopes and Turck-Chieze, 1994;Roxburgh and Vorontsov, 1994;Lopes et al., 1997;Roxburgh and Vorontsov, 2001). They are also visible in the large separation, or even better in second differences, or even higher differences (e.g. Basu et al., 1994;Basu, 1997;Mazumdar and Antia, 2001;Ballot et al., 2004a). The second differences should almost vanish without glitches. Variations in δ 2 ν n, can be modelled as the sum of a small trend and two glitches: This parametrization of the glitch envelopes is adapted from Houdek and Gough (2007). Glitches are also visible in small separations and ratios, but oscillation periods of signatures are not linked to the acoustic depth but to the acoustic radius (time for a sound wave to propagate to the glitch from the stellar centre). As a consequence, the helium ionisation signature almost disappears in these variables (for more details, see Roxburgh and Vorontsov, 2003;Roxburgh, 2005Roxburgh, , 2009a. Measurement of glitches requires high precision on individual mode frequencies. In CoRoT observations of HD 52265, a BCE glitch was visible  and has been analysed by Lebreton and Goupil (2012). Their analysis suggests that the penetrative distance of the convective envelope in the radiative zone reaches 0.95 H p (6% of the total radius), which is significantly larger than what is found for the Sun. Measuring glitches in F-type stars is challenging because uncertainties of frequencies are high due to mode broadening (Eq. 31), reinforced by the overlap of = 0 and 2 modes. However, Brito and Lopes (2017) succeeded in measuring the position of helium ionisation zone in HD 49933 using CoRoT data. Glitches have also been intensively studied in Kepler solar-type stars. A first work on 19 stars with one-year observations has been carried out by Mazumdar et al. (2014), then extended to the LEGACY sample (66 stars observed for up to 4 years) by Verma et al. (2017). Figure 28 shows an example of second differences fitted with oscillatory components due to BCE and He II ionisation zone. Helium glitches can be robustly measured in many stars, whereas BCE is more difficult to measure for super-solar mass stars. These measurements can now be used to constrain properties of the helium ionisation zone and helium abundance in the envelope.
The question of the extent of convective cores in solar-type stars is still debated. It depends on the mixing processes in the core and also on the details of the nuclear reaction rates. A recent study of α Cen A combining observations in astrometry, spectroscopy, interferometry and ground-based seismology shows that the nature of its core is still uncertain (Bazot et al., 2016). Signatures of convective cores can be found in seismic variables. Small separations δν 01 and δν 10 and the ratio of small to large separation r 01 and r 10 are very sensitive to the core structure (e.g. Provost et al., 2005;Popielski and Dziembowski, 2005;Deheuvels et al., 2010b;Silva Aguirre et al., 2011a;Cunha and Brandão, 2011). Convective cores create glitches with very long periods in δν 01 , δν 10 , mainly visible as a slope. As opposed to the BCE position directly deduced from observations, interpreting the slope as a core extent is not fully model independent. Indeed, several effects contribute to the slope: the amount of hydrogen in the core (no matter whether a convective core is present or not), the size of the convective core and the sharpness of the discontinuities. Using this idea Deheuvels et al. (2016) have been able to measure the extents of convective cores for eight main-sequence stars. From them, they were able to calibrate the quantity of extra mixing needed in the core to reproduce the seismic observations.
In the case of helioseismology, fine constraints on the internal structure were achieved by inversion techniques. This was possible thanks to the observations of many intermediate-and highdegree modes. Structure inversions for other stars, based only on low-degree modes without strong external constraints on mass and radius are challenging. In this context, Reese et al. (2012) proposed an inverse method to accurately (0.5% accuracy) estimate stellar mean densities. Inversion efforts have been pursued for several Kepler targets, especially 16 Cyg A and B by Buldgen et al. (2016a,b). Doing so, they reduced the uncertainties to 2% on mass, 1% on radius, and 3% on the age for 16 Cyg A. New inversion methods to constrain convective regions have been proposed by Buldgen et al. (2018) and should be tested soon on Kepler observations. Finally, a new method called "Inversion for agreement" has been proposed by Bellinger et al. (2017). This method takes into account imprecise estimates of stellar mass and radius as well as the relatively small amount of modes available. Because the result is independent of models, it can be used to test their inferences. Thus, the results obtained on 16 Cyg A and B showed that the core sound speed in both stars exceeds that of the models.

Stellar rotation
One of the main results obtained with helioseismology is the accurate determination of the solar differential rotation in the convective zone (e.g. Thompson et al., 1996Thompson et al., , 2003Howe, 2009), as well as the nearly constant rotation rate in the radiative interior down to ≈ 0.2 R , where the measured low-degree p modes do not provide enough sensitivity below this radius (e.g. Chaplin et al., 1999a;García et al., 2004aGarcía et al., , 2008c. Gravity modes are needed to properly constrain the rotation within the core (e.g. Mathur et al., 2008). However, although the detection of individual g modes in the Sun is still controversial (e.g. Appourchaux et al., 2010;Schunker et al., 2018) the latest results obtained considering these modes or their period spacing suggest a faster rotation rate in the core (García et al., , 2008a(García et al., , 2011Fossat et al., 2017).
When extending the analysis of internal rotation to other stars, the precision and the extension of the region to be explored depends on our ability to probe their inner regions. In MS stars, only pure acoustic modes have been unambiguously characterized so far and, therefore, only the outer convective zone and the outer part of the inner radiative zone can be probed by asteroseismology. Only when stars enter the sub-giant region, can mixed modes be measured (e.g. Deheuvels et al., 2010a;Benomar et al., 2012) and rotation of the inner zones and the core can be obtained (e.g. Deheuvels et al., 2012Deheuvels et al., , 2014. The first unambiguous determination of the average internal rotation and the rotation inclination axis of a MS solar-like star using asteroseismic techniques was performed by Gizon et al. (2013) on the CoRoT target HD 52265. To prove the validity of their measurements, Gizon et al. (2013) compared the asteroseismic rotation with the spectroscopic value, the surface rotation period, and starspot modeling obtained from the analysis of this CoRoT light curve. The result is shown in Fig. 29.
This example illustrates all of the information provided by the study of continuous highprecision photometry. On the one hand, a direct determination of the rotation period can be obtained by the analysis of the light curve, either in the time domain, or by studying the lowfrequency part of the temporal power spectrum. In addition, starspot modeling can also provide additional information such as the rotation inclination angle (e.g. Mosser et al., 2009a;Lanza et al., 2014). On the other hand, seismology yields a weighted average of the internal rotation heavily biased towards the surface during the MS, as well as the rotation inclination angle. In the next two sections, we will describe in more detail these two types of analyses.

Photospheric rotation from the study of the photometric light curves
The crossing of cool spots over the visible disk of a star produces a modulation in the photometric signal which is proportional to the rotation period of the star, P rot , at the latitudes where spots and active regions exist. Several works have been dedicated to the study of the extraction of these rotation periods for CoRoT, Kepler and the K2 missions. Thus, in the last decade, it has been possible to retrieve the rotation periods of thousands of stars (e.g. McQuillan et al., 2013;Nielsen et al., 2013;Walkowicz and Basri, 2013;McQuillan et al., 2014;Leão et al., 2015).
There are several different -but complementary -ways to extract the information on the rotation period from light curves. In this section, we do not pretend to provide an exhaustive review of all of the techniques and results on this topic; we focus mostly on what has been done related to asteroseismic studies of MS solar-like stars.
The careful study of the low-frequency part of the power spectrum can be used to extract the rotation period by selecting the highest peak in this frequency region (e.g. Barban et al., 2009;Campante et al., 2011;Nielsen et al., 2013). However, sometimes the highest peak can be the Index Fig. 29 Angular velocity of HD 52265, Ω, normalized by the angular velocity of the Sun (Ω Sun /2π=0.424 µHz) as a function of the sine of the inclination of the stellar rotation axis to the line of sight, sin i. The black diamond with error bars is the best-fit seismic result. The red shaded region represent the 1-σ seismic constraint in both parameters. The two horizontal lines represent the angular velocity obtained form the analysis of the low-frequency range of the power spectrum. The green ellipse is the 1-σ result obtained from the starspot modeling of the light curve. Finally, the dashed (observations) and the solid (1-σ errors) blue curves are the constraints obtained from ground-based spectroscopy. Figure adapted from Gizon et al. (2013).
second or even the third harmonic of the rotation period instead of the first corresponding to the true rotation period. An example is given in Fig. 30, where the light curve and the longer periods of the period-power spectrum of KIC 4918333 is shown. On the left-hand panel, the first 310 days of the light curve are analyzed and the tallest peak corresponding to a rotation period of 9.8 ± 0.8 days is shown. A second peak is also visible at a period of 19.5 ± 1.2 days but with a smaller amplitude. When this analysis is extended to 1250 days, the period of 19.5 ± 1.2 days is the most prominent one without ambiguity. Hence, when analyzing the rotation period from the study of the power spectrum, it is important to check any possible signal at twice the period corresponding to the highest peak.
This situation occurs when two active regions develop in stars with approximately 180 • separation in longitude. This is clearly visible in the light curve shown in the left panel of Fig. 30. Two modulations of different amplitudes (one large and one small) are separated by approximately 10 days. As a consequence, the largest peak in the spectrum is not the ≈ 20-day periodicity but the ≈ 10 days.
To avoid this possible ambiguity, McQuillan et al. (2013) used the autocorrelation function, ACF, as a robust method to infer the rotation period directly from the Kepler light curves in the time domain. They showed that in stars such as KIC 4918333, the highest peak in the ACF was the first harmonic and not the second and thus the retrieved rotation periods were robust. In Fig. 31 the ACFs of the analysis of the first 310 and 1250 days of KIC 4918333 are shown. With this methodology -contrary to the direct analysis of the power spectrum -the highest peak is the one at ≈ 20 days for both length of the time series. Because the photometric light curves are not free from instrumental effects, time-period diagrams have been used since the CoRoT mission to look for the rotation periods in stars. In the mean time, we also checked any possible "glitches" in the signal that could produce a spike in the low-frequency part of the spectrum that could perturb the determination of the rotation (e.g. García et al., 2009;Mathur et al., 2010b). This way, segments where instabilities appear in the light curve could be easily inspected and removed from the analysis. To extract the rotation period from the time-period diagrams, a projection on the period axis is done, and a fitting of the highest peak is then performed as if it was a normal power spectrum. The application of this time-period methodology to KIC 4918333 is shown in Fig. 32.
Since the different methods used to calibrate the data can filter the stellar signal in different ways, while leaving other instrumental signals in the final light curves, it is recommended to use several types of data calibrations when studying the rotation of a star (e.g. García et al., 2013aGarcía et al., , 2014aCeillier et al., 2016;Buzasi et al., 2016). It has also been shown that each method to extract the rotation (e.g. ACF, time-period diagrams, the direct analysis of the low-frequency part of the spectrum) works better in some circumstances or for some type of stars. Hence, the most reliable procedures to retrieve the rotation periods are those that combine different calibration procedures and analysis techniques (Aigrain et al., 2015).
The average rotation is not the only quantity that can be inferred from the analysis of the light curves. Differential rotation in latitude ∆Ω can also be obtained in some cases (e.g. Fröhlich et al., 2012;Lanza et al., 2014;Reinhold and Gizon, 2015), by directly studying the number and structure of the peaks in the low-frequency part of the spectrum or by performing starspot modeling (e.g. Croll, 2006;Fröhlich, 2007;Nielsen et al., 2013;Walkowicz et al., 2014;Lanza et al., 2016). Thanks to these analyses, general trends have been found. For example, the dependence of ∆Ω with the rotation period is weak and it slightly increases with T eff Index Fig. 32 time-period analysis of the light curves showed in Fig. 30 computed using a Morlet wavelet between 0.5 and 100 days on a logarithmic scale following Torrence and Compo (1998). The green-crossed area is the cone of influence corresponding to the region where it is not possible to obtain reliable results. As in Fig. 30, the global wavelet power spectrum is shown in the right of each of the time-period panels. The red-shaded regions depict the tallest peak in this diagram. The green line represents the best Lorentzian fit to this peak. The dotted line represent the 95% confidence level.
in the range 3500 to 6000 K , while the relative differential rotation, ∆Ω/Ω, increases with the rotation period (Reinhold and Gizon, 2015). In addition, Reinhold and Arlt (2015) were able to discriminate solar and antisolar differential rotation (i.e., to identify the sign of the differential rotation at the stellar surface) using peak-height ratios of the first harmonics of the differential rotation peaks. However, a theoretical and numerical study by Santos et al. (2017) showed that the peak-height ratios are essentially a function of the fraction of time the spots are visible. This time is related to how strongly the spot modulation follows a sinusoidal form. Hence, depending on the rotation inclination angle with respect to the line of sight and on the location of the spots, the inferred sign of the differential rotation can be wrong.
Signatures of the change in latitude of the spots have also been measured for example in KIC 3733735 (Mathur et al., 2014a). During the periods with low activity, the mean rotation period is ∼ 3 days, while during the periods of maximum activity the rotation period decreases to 2.54 days (see Fig. 37). Similarly as for the Sun, this behavior can be explained by the existence of surface differential rotation and a change in the latitude where the spots emerge as the magnetic cycle progresses.

Internal rotation through asteroseismic measurements
As we have already said in this review, only acoustic modes have been characterized in mainsequence solar-like dwarfs so far. Therefore, the information that one can obtain from seismology comes from the careful analysis of acoustic modes. Unfortunately, although low-degree p modes penetrate deep in the stellar interior, they spend only a small fraction of their time in the deep interior and the amount of information that they can provide from these regions is small. This is illustrated in Fig. 33 where we show the rotational kernels, K n, , of the dipolar and quadripolar modes of radial orders n = 10 and n = 25 for three planet-hosting stars: Kepler-25 (KIC 4349452, M=1.26 ± 0.03 M , Benomar et al., 2014), HAT-P7 (KIC 10666592, M∼1.59 ± 0.03 M , Benomar et al., 2014) and the Sun. The two radial orders n = 10 and 25 cover the typical observational range for theses stars. In spite of the different masses of the three stars and the different position of the base of the convective zone, the kernels of these modes are nearly identical (see also Lund et al., 2014). However, the position of the base of the convective zone depends on the mass of the star. Therefore, more massive stars have shallower convective zones and thus, they have a larger contribution from the radiative zone to the average internal rotation rate. It is, however, possible to extract some general properties. These kernels have larger amplitudes near the surface. They are also denser in these outer regions implying that the waves spend more time there than in the inner layers. As a consequence, the sensitivity to the rotation is larger towards the surface of the star. Moreover, the kernels in the radiative zone are almost linear functions of the radius above ∼ 0.15 R up to the base of the convective zone. Therefore, each observed mode probes the radiative zone uniformly.
The first consequence is that the rotation rate extracted from seismology and the rotation period extracted from the surface (either from the study of the photometric variability or from spectroscopy once the inclination angle is known) should be similar as already explained when described the first asteroseismic detection of the rotation in HD 52265 (see Fig 29). Interestingly, any significant difference between the surface and the asteroseismic rotation rate should indicate the existence of differential rotation, either from the surface or from the external convective zone inside the stars. Gizon et al. (2013) found that both determinations of the rotation agreed within one sigma in all cases, suggesting that the differential rotation should be weak in these stars.
A comparison of the asteroseismic rotation rate (from Nielsen et al., 2014) and their surface rotation period was performed by Nielsen et al. (2015) for five Kepler targets with masses in the range 1.02 to 1.2 M . Because of the large errors in any inference of the rotation period, the asteroseismic rotation rate can be considered a good proxy of the surface rotation only under the assumption that there is no large differential rotation neither at the surface nor in the convective envelope. Hence, asteroseismic rotation rates can be used for example in gyrochronology studies when no other precise rotation rates are available (e.g. Davies et al., 2015;Nielsen et al., 2015).
Based on similar studies of the rotational kernels but using inversion methods, Schunker et al. (2016a) suggested that it could be interesting to use many stars in order to reduce the observational errors and be able to, for example, constrain the sign of the radial differential rotation. Moreover, it has been demonstrated using an "ensemble-fit" of 15 stars across the main sequence, that it would be possible to distinguish between solid rotation and radial differential rotation of around 200 nHz using observable splittings of angular degrees 1 and 2 (Schunker et al., 2016b).
Another simplified, complementary approach to rotational inversion was proposed by Benomar et al. (2015) taking advantage of the relative simplicity of the stellar interiors of MS solar-like dwarfs. In principle, it is possible to consider that the region probed by the measurable acoustic modes is composed of two different zones: the outer convective and the inner radiative zone. Therefore, the rotational splitting δν n,l can be expressed as a weighted average of the rotation rate, f rad , f conv , in each zone as follows: where I rad and I conv are the integrals of the rotational kernels in the radiative and convective zones defined as: with I rad + I conv = 1. The position of the base of the convective zone is denoted by r bcz and R is the radius of the star. In this simplified model, assuming that f rad is a good approximation of the rotation rate in the radiative zone, it is necessary to suppose that: a) the two zones rotate uniformly with different rates (about the same axis), b) f conv is approximately equal to the surface rotation (as in the solar case), c) rotational splittings remain nearly constant over the observed range of modes. Hence, the average rotational splitting < δν n, >, can be used as a representative value of the seismically measured internal rotation rate, f seis . This last assumption is validated because the integrals of  the observed mode kernels shown in Fig. 33 are nearly identical. Therefore, it is possible to express Eq. 58 as follows: where <> denotes that the corresponding parameter is averaged over the observed range of modes. Knowing the kernel integrals from a model and the average internal rotation rate from the seismic splittings, it is necessary to have a good determination of the surface rotation rate in order to infer the rotation in the internal radiative region with Eq. 61. Two possibilities are available: on the one hand, the surface rotation obtained from photometric variability can be different form the average surface rotation because it is weighted towards the active latitudes (30 • in the solar case). On the other hand, the spectroscopic value is averaging the full visible disk and it can be considered a better approximation. However, in this latter case, what is measured is the projected velocity into the line of sight velocity: Thus, to extract f surf it is necessary to know the rotation inclination angle, i, from seismology and the radius of the star from the same models used to compute the rotational kernels.
Using 22 stars observed by CoRoT and Kepler -of masses in a range ∼1.07 to ∼ 1.56 M for which f surf is available from spectroscopy and photometric variability - Benomar et al. (2015) deduced that nearly uniform internal rotation (between the radiative zone and the surface) is common in other solar-like stars from masses up to ∼ 1.5 M . Thus, the Sun is not an isolated case in stellar evolution and efficient angular momentum transport mechanisms are required in the Fig. 35 Ratio of the average radial differential rotation in the radiative zone obtained from Eq. 61 and the surface rotation (obtained from spectroscopy) as a function of the stellar inclination angle (obtained from seismology). Colors represent the stellar mass. Grey lines indicate a differential rotation of a factor of 2. The dashed area indicates a negative rotation rate in the interior. Figure adapted from Benomar et al. (2015). main sequence to achieve such rotation profile (e.g. Mathis, 2013;Marques et al., 2013). The ratio < f rad > /f surf − 1 is shown in Fig. 35.
It is important to note that only one star, KIC 9139163, significantly departed from the one-toone relation shown in Fig. 35, offering a scenario where the interior could be spinning much faster than the surface. Interestingly, this star is among the youngest and more massive of the sample suggesting a scenario where the angular momentum transport processes that are responsible for the quasi-uniform internal rotation might not have had enough time to complete their work.

Stellar magnetic activity and magnetic cycles
In distant stars, as was the case for rotation, asteroseismic observations provide two different, but complementary, ways to study magnetic activity in general and magnetic activity cycles in particular. On the one hand, magnetic variability can be measured on different time scales by directly analyzing the average luminosity flux modulation in the light curve or its fluctuation as a function of time (e.g. Basri et al., 2010Basri et al., , 2011García et al., 2014a). It is out of the scope of this review to describe in details the methodologies and results associated to these studies. However, in the next paragraph we will provide a brief overview of them. On the other hand, magnetic variability can be studied through the characterization of the temporal evolution of the oscillation modes, i.e., measuring the variations of the frequencies, amplitudes, and line widths of the acoustic modes.
As discussed in the section about internal rotation, to reach the core, mixed or gravity modes are required. Unfortunately, those modes have not been unambiguously detected in main-sequence solar-type dwarfs. Therefore, nothing can be said about the magnetic field in the core of these stars. However, studying red giant stars, Stello et al. (2016) evoked the possibility of having magnetic fields of dynamo origin in the convective core of stars with masses greater than ∼ 1.2 M . The origin of the depressed dipolar modes seen in Kepler observations (e.g. Mosser et al., 2012;García et al., 2014c) has been suggested to be magnetic, although the exact nature of phenomena is still debated (e.g. Fuller et al., 2015;Mosser et al., 2017;Papaloizou, 2017, 2018).

From the direct analysis of the light curves
High-quality, photometric time series of several months to years duration are now available from CoRoT, Kepler, K2 and TESS missions. They allow a detailed study of the photospheric stellar variability at different time scales and to compare them with the Sun (e.g. Gilliland et al., 2011;Basri et al., 2013, and references therein). However, not all this variability can be directly associated with magnetism. Different physical phenomena (pulsations, convection or rotation) can indeed coexist on the same time scales. To overcome this problem, it has been proposed to study the time scales associated with the rotation period of the star as an indicator of a magnetic origin (e.g. Mathur et al., 2014b,a;García et al., 2014a). For magnetically active stars developing starspots at their surface, the global brightness is modulated as a function of the position and the size of the spots over the visible stellar disk with a period that is related to the rotation period of the star at the active latitude where the spot appears (e.g. Lanza, 2010;Fröhlich et al., 2012;Lanza et al., 2014;McQuillan et al., 2014). As the magnetic cycle evolves, the number and size of the spots change. An example of the solar photospheric and Doppler-velocity time series is shown in Fig. 36. Hence, a photospheric magnetic activity proxy, S ph , can then be constructed by measuring the dispersion of the light curve on time scales of five times the rotation period . The same procedure can be applied to Doppler velocity observations and another proxy, called S vel , has been defined. The GOLF red-wing period denoted by the two vertical dashed lines where the instrument is sensitive to higher layers in the atmosphere is represented in red (see for details García et al. (2004b), García et al. (2005), and Jiménez-Reyes et al. (2007)). Figure adapted from Salabert et al. (2017).
The validity of this proxy as a function of magnetic activity has been verified using the Sun (Salabert et al., 2017) and compared to the chromospheric activity of other Kepler targets including 18 solar analogs (Salabert et al., 2016a;García, 2017). This proxy (or other similar metrics) is currently being used to search for magnetic activity cycles or magnetically-related trends that could be part of magnetic cycles longer than the duration of the observations for stars of many different spectral types (e.g. Oláh et al., 2009;Mathur et al., 2013bMathur et al., , 2014bVida et al., 2014;Ferreira Lopes et al., 2015). Finally, by combining magnetic proxies and seismology it has been possible to revisit the age-rotation-activity relations (Karoff et al., 2013b;García et al., 2014a) and to establish that the photospheric magnetic activity of the Sun is comparable to other mainsequence solar-like pulsating stars , in particular, compared to other solar analogs of similar age (do Nascimento et al., 2014;Salabert et al., 2016a;Beck et al., 2017). In this sense, we can conclude that the Sun is a normal star in terms of its surface magnetism when compared to its siblings. An example of the study of the surface dynamics (rotation and magnetism) using Kepler photometry is given in Fig. 37. In the top panel the relative measured flux is given. The envelope of the surface brightness shows an increase between the day ∼ 550 and ∼ 1200 of the mission that can be interpreted as an increase of the magnetic activity of the star. In the middle panel of Fig. 37 a time-period diagram is computed. It shows two main bands of power (depicted by the yellow horizontal dashed lines) at 3 and 2.54 days corresponding to the main average rotation rate of the star during these two seasons. Projecting the time-period diagram onto the time domain around the periods given before (from to 2 to 6 days), a magnetic proxy can be built (see bottom panel of Fig. 37), confirming the existence of an on-going stellar activity cycle, with a season of maximum activity where the starspots are located at longitudes of faster rotation (2.54 days), while the band at a slower rotation rate is dominating the location of the spots during the minimum of magnetic activity. This behavior is similar to that observed in the Sun which produces the so-called butterfly diagram (e.g. Hathaway, 2015).

From asteroseismology
To study magnetic activity cycles with asteroseismic techniques, it is important to use the Sun as a reference because we can characterize in detail its surface magnetism and look for correlations with different observed features that are today impossible to obtain for distant stars. In this way, we can then apply this knowledge to other stars. At the very beginning of helioseismology, a correlation between the acoustic-mode frequencies and solar magnetic activity was found (van der Raay, 1984;Woodard and Noyes, 1985;Fossat et al., 1987;Pallé et al., 1989;Elsworth et al., 1990), i.e., the frequencies shifted towards higher values as the 11-year magnetic cycle progressed. Later, it was discovered that the frequency shifts increased with frequency for intermediate-and low-degree modes (Libbrecht and Woodard, 1990;Anguera Gubau et al., 1992), i.e. modes at higher frequency had a larger frequency shift than modes at low frequency. This frequency dependencehigh-frequency modes have outer turning points compared to low-frequency modes -is the main change of the mode parameters when the effect of the mode inertia is removed. Considering also that there was no significant dependence of the frequency shifts with the degree of the modes led to the conclusion that the perturbations related with the magnetic activity cycle were confined to a thin layer very close to the photosphere (e.g. Libbrecht and Woodard, 1990;Goldreich et al., 1991;Nishizawa and Shibahashi, 1995;Basu et al., 2012). Moreover, the Sun has not only a 11-year periodicity. Shorter -quasi-biennial-modulations have been measured in the Sun (e.g. Benevolenskaya, 1995) in several magnetic activity proxies and also confirmed by seismology (Fletcher et al., 2010;Broomhall et al., 2011b;Simoniello et al., 2012). The existence of several time scales in the modulation of the magnetic proxies is not peculiar to the Sun and many other stars show several long and short periods that could be interpreted as magnetic cycles with an "active" and "inactive" phase (e.g. Baliunas and Soon, 1995;Brandenburg et al., 1998;Saar and Brandenburg, 1999;Böhm-Vitense, 2007;Brandenburg et al., 2017). Last but not least, important differences were found between the surface magnetic proxies and the frequency shifts during the last extended minima between Solar Cycles 23 and 24 (Broomhall et al., 2009;Salabert et al., 2009). While no activity was measured in standard magnetic proxies, the frequency shifts showed a quasi normal behavior leading to the conclusion that the magnetic perturbations in the subsurface layers were still strong (Salabert et al., 2015).
All of these early findings on the modulation of the frequency shifts with the magnetic cycle were soon followed by the study of the evolution of all of the other oscillation-mode parameters. Hence, the measurement of a reduction of the mode heights, an increase of the mode linewidths and a reduction of the velocity power with the magnetic activity cycle were found. Finally, the energy that is supplied to the modes remained constant along the magnetic cycle within the uncertainties (e.g. Pallé et al., 1990;Chaplin et al., 2000;Komm et al., 2002;Jiménez-Reyes et al., 2003;García et al., 2013b;Broomhall et al., 2015). It is therefore important to take these changes into account when performing the fitting of the acoustic modes in active stars. All of these variations modify the simple Lorentzian profiles of the modes introducing a bias of the retrieved parameters when simple Lorentzian profiles are used to fit spectra obtained from time series that are long compared to the timescale of the variations introduced by the activity (Chaplin et al., 2008a;Chaplin and Basu, 2008). In particular, frequency shifts modify the central frequencies of the Lorentzians. Thus, when the resultant profile is fitted by a single Lorentzian function, the inferred seismic parameters are biased such as the large frequency separations as demonstrated by Broomhall et al. (2011a). An example of this distorted Lorentzian profile is given in Fig. 38 for a mode with a linewidth of Γ =1.5 µHz and four different maximum-to-minimum frequency shifts: 0, 0.6, 1.5, and 2.25 µHz. Fig. 38 Lorentzian profile of a mode with a linewidth Γ =1.5 µHz (grey), and three other pseudo-Lorentzian profiles resulting from the analysis of a long time series with three different maximum-to-minimum frequency shifts of 0.6, 1.5, and 2.25 µHz, black, red and blue lines respectively. Figure adapted from (Chaplin and Basu, 2014).
One exception to this is the characterization of low-degree low-order p modes for which the variations induced by the magnetic activity cycle are too small, around ∆ ν /ν ∼ 10 −5 nHz for the solar case (e.g. García et al., 2001).
To study stellar magnetic activity cycles it is necessary to have long time series covering all -or a significant part -of the magnetic cycle. The longest time series available today for asteroseismology were obtained by the Kepler mission covering slightly more than four continuous years. Due to the rather well-known relation between the rotation period of stars and the length of the magnetic cycle (e.g. Noyes et al., 1984b,a;Brandenburg et al., 1998;Böhm-Vitense, 2007;Brandenburg et al., 2017), the study of full magnetic cycles is limited to fast rotating stars (P rot ≤ 7 days) for the active, "A", branch and for stars with intermediate rotation rates (P rot ≤ 18 days) for the less-active, "I", sequence (see Fig. 39).
The relation between magnetic activity cycle periods and stellar evolution is also well known (e.g. Brandenburg et al., 1998). Recently, Brandenburg et al. (2017) re-visited this relation in light of the most recent and longest spectroscopic observations from Egeland (2017) who combined datasets from several instruments. Unlike previous interpretations where young stars would evolve along the active A branch, they now believe that all stars younger than 2.3 Gyr are capable of exhibiting longer and shorter cycle periods. If their calculations are correct, for the G pulsating dwarfs HD 76151 and KIC 10644253 shown in Fig. 39, longer periods in the 12 to 16 year range are expected, and may have already been found in HD 76151 (Egeland, 2017) and in ι Horologii, Fig. 39 Magnetic activity cycle period versus the rotation period of stars showing the existence of two branches an "active" (red line denoted by A and red symbols) and an "inactive" sequence (blue line denoted by B and blue symbols). Several stars have been added to the original diagram from Böhm-Vitense (2007) such as ι Horlogii (Metcalfe et al., 2013), HD 76151 (Egeland, 2017), and the seismic targets 18 Sco (Hall et al., 2007), HD 49933 (García et al., 2010) and the Kepler targets KIC 8006161 and KIC 10644253 (Salabert et al., 2016b;Kiefer et al., 2017;Karoff et al., 2018). The green shaded region represents the region where full stellar cycles can be studied with the Kepler data. The violet shaded region represents the region where longer activity cycles could be uncovered by Kepler observations if they are in a regime where the magnetic activity changes significantly. The letter H indicates stars in the Hyades, crosses indicate stars on the "active" A sequence, and asterisks indicate stars on the "inactive" B branch. Squares around the crosses show stars with a color B-V < 0.62. Triangles indicate secondary periods for some stars on the active sequence.
HD 17051, where Flores et al. (2017) measured a long-term activity cycle of about five years fitting the "active" branch in the Böhm-Vitense diagram. For the solar analogue 18 Sco, it would be interesting to discover either a second shorter modulation as in the Sun or a longer cycle period.
An interesting result was obtained by Mathur et al. (2014a) from the analysis of the temporal evolution of S ph . They selected 22 solar-like stars with detected solar-like oscillations and with rotation periods shorter than 12 days, i.e. the limit to observe at least half of a cycle from Fig. 39. Only two of the stars show a cyclic-like variation, while two others showed a decreasing or increasing trend in the S ph temporal variation. Five stars show modulations (or beating) due to long-lived spots at two different active longitudes with different rotation rates. The rest of the stars show no cycle-like behavior although they showed surface magnetic activity. Therefore, although it is premature to infer firm conclusions from a so small sample of stars, it seems that fast rotating stars can exhibit magnetic activity but without any magnetic cycle (or at least much longer than it could be expected from Fig. 39. The only correlation found was between the S ph and the rotation period for stars showing a beating between long lived spots at different rotation rates, in the direction of higher S ph for longer P rot . To perform asteroseismic studies of magnetic activity, it is then necessary to first measure the oscillation properties and then study their evolution with time. As it has been said, solar magnetism reduces the amplitude of the solar modes. Therefore, a very active star would probably have oscillation modes with very small or even not detectable amplitudes. This effect was first observed by CoRoT in two main-sequence targets: HD 175726 (Mosser et al., 2009b) and in HD 49933 (García et al., 2010). Later, using a Kepler sample of main-sequence stars, Chaplin et al. (2011a) showed that there was a clear correlation between stellar magnetic activity and the amplitude of the stellar oscillations. The most active stars in the sample did not show measurable oscillation modes.
To look for the evolution with time of the seismic properties, the light curve is divided into small segments for which the characterization of the oscillation modes is carried out. The length of the sub-series is found as the best trade-off between frequency resolution and the number of sub-series to be analyzed. The longer the series the better the precision on the extracted parameters (e.g. Régulo et al., 2016). However, it can be very challenging to obtain individual p-mode frequency shifts that can be as small as half a µHz for short time series. This is why a global method was developed in the early days of helioseismology by Pallé et al. (1989) to obtain averaged frequency shifts by computing the cross correlation of the p-mode hump computed from the PSD of each sub-series in comparison to a reference one. This reference can be either the PSD of one of the sub-series or the average spectrum of all of them (e.g. Pallé et al., 1989;García et al., 2010;Régulo et al., 2016;Salabert et al., 2016b;Kiefer et al., 2017;Santos et al., 2018a).
The first attempts to detect differences in the seismic parameters associated with magnetic activity were reported by Fletcher et al. (2006). They analyzed observations made with the startracker on the WIRE satellite of α Cen A and compared the obtained frequencies with previously obtained frequencies measured by Bouchy and Carrier (2002) and Bedding et al. (2004). Fletcher et al. (2006) conjectured that the average difference of about 0.6 ± 0.3 µHz (∼ 2σ) in the oscillation frequencies could be due to an ongoing activity cycle in α Cen as the WIRE observations were taken 19 months before the others.
Successful asteroseismic measurements of magnetic activity started in 2010 when García et al. (2010) unveiled the presence of a magnetic activity cycle in a star other than the Sun: the CoRoT target HD 49933, which was observed during the first two CoRoT runs with a six month interval Benomar et al., 2009b). They found a modulation of more than 120 days in the three indicators considered: frequency shifts -measured globally through cross-correlation techniques as well as by extracting the individual p-mode frequencies -, amplitude modulation of the p-mode bump and in the dispersion of the light curve. Moreover, as also observed in the Sun, the modulation detected in the amplitude and the frequency shifts of the p modes was anticorrelated (see Fig. 40).
The analogies between the magnetic cycles of both stars continue. As in our Sun, the frequency shifts measured in HD 49933 present a frequency dependence with a clear increase with frequency. However, the maximum frequency shift is about 2 µHz around 2100 µHz, 4 times bigger than in the solar case . Similar variations are obtained between the = 0 and = 1 modes computed independently (see Fig. 41). At higher frequencies, the frequency shifts show indications of a downturn followed by an upturn for both low-degree modes = 1 and 2. Therefore, the frequency variation of the p-mode frequency shifts of this star has a comparable shape to tthat observed in the Sun, which is understood to arise from changes just beneath the photosphere (e.g. Basu, 2016).
Complementary observations of HD 49933 were taken in the Calcium H and K lines since 13 April, 2010, showing that this is an active star with a Mount Wilson S-index of 0.3 confirming previous conclusions (Mosser et al., 2005). The temporal evolution of the Mount Wilson S-index showed that a magnetic activity cycle is ongoing in this star. However, longer observations are needed to clearly establish the length of the cycle and whether or not the magnetic cycle is regular. The same methodology was applied to three other CoRoT targets complemented by groundbased analysis of the chromospheric magnetic activity done with the NARVAL spectropolarimeter located at the Bernard Lyot 2 m telescope at the Pic du Midi Observatory (Aurière, 2003) in order to study any possible hint of magnetic activity cycles . Interestingly one star, HD 181420 ), seems to be in a stationary regime without any visible change of the activity during the observations. This is an unexpected result because this star rotates rapidly (2.6 days) and we were expecting to see some indications of a magnetic activity cycle. For the other two stars, HD 49385 (Deheuvels et al., 2010a) shows a small increase of activity at a 1-σ level but not confirmed by our spectroscopic measurements, while HD 55265 (Soriano et al., 2007;Ballot et al., 2011b) presents a small variation of the seismic parameters, also at a 1-σ level, and in the spectroscopic observations performed by NARVAL. That could indicate that this star was observed during the rising phase of a long magnetic activity cycle.

Influence of metallicity on magnetic activity
The longer observation period of the Kepler main mission (up to four continuous years) allows one to better track temporal changes in the seismic parameters for main-sequence solar-like stars. This is illustrated in Fig 43 where the temporal evolution of the seismic parameters and the photospheric and chromospheric activity proxies are shown for KIC 8006161 (HD 173701) and for the Sun. All of the figures are in the same scale for both stars. KIC 8006161 is a very interesting target because it is a solar analogue (Karoff et al., 2018) with M = 1.00 ± 0.03 M , R = 0.93 ± 0.009 R , T eff =5488 ±77 K, P rot = 21 ± 2 days, age = 4.57 ± 0.36 Gyr, and a metallicity that is twice the solar value (0.3 ± 0.1 dex).
KIC 8006161 has a comparable magnetic activity cycle period (∼ 7.4 yr) deduced from more than thirty years of chromospheric observations (top plot in Fig. 43), being in the rising phase  of its cycle during the four years of Kepler observations. Interestingly, the amplitude of all the temporal variations recorded are much larger than the corresponding solar values. Karoff et al. (2018) conjectured that these differences could be a consequence of the higher metallicity of the star. An increase of the stellar metallicity produces larger opacities and thus a larger internal temperature gradient. Therefore, the Schwarzschild criterion for convection (Schwarzschild, 1906) is satisfied deeper inside the star leading to a deeper convective zone (van Saders and Pinsonneault, 2012). Theoretical studies and numerical simulations have shown that larger convective zones induce larger differential rotation (Brun et al., 2017) and thus a stronger magnetic dynamo (Bessolaz and Brun, 2011). Although it is not possible with the present observations to infer how strong the dynamo is, the surface differential rotation is larger than that in the Sun, reinforcing the conclusions. Unfortunately, firm conclusions cannot be derived with only one observations and the study of a larger sample of stars with different metallicities is necessary to completely understand the influence of metallicity on the length and strength of the magnetic activity cycles.

Relation between the frequency shift strength with effective temperature and age
We have just presented an explanation for a possible relation between metallicity and the magnetic cycle strength. Using the ensemble analysis of frequency shifts from asteroseismology, it is possible to depict other correlations between different parameters. The relation between the effective temperature and the frequency shifts is particularly interesting as two different scaling relations have been proposed by Chaplin et al. (2007a) and Metcalfe et al. (2007). On the one hand, Chaplin et al. (2007a) proposed that the frequency shifts linearly scale with the amplitude of the activity cycles ∆R HK (defined by Saar and Brandenburg, 2002). This is validated in the Sun because the frequency shifts change linearly with the temporal evolution in the Mg II lines (Chaplin et al., 2007b). With this scaling, the frequency shifts generally increase towards lower temperature and decrease with age as seen in Fig. 44  . Fig. 44 Frequency shifts, δνcyc, as a function of effective temperature computed for different Padova isochrones (different line styles, Girardi et al., 2002Girardi et al., , 2004. Black lines follow the scaling relation of Chaplin et al. (2007a) while blue lines follow Metcalfe et al. (2007). On the other hand, Metcalfe et al. (2007), proposed that the frequency shifts are also proportional to ∆R HK but scaled with a factor that depends on the depth of the perturbation normalized by the mode inertia. In this case, the frequency shifts generally increase towards higher temperature and decrease with age (see Fig. 44). While both scalings yield a decrease of the frequency shifts with age, they have different predictions with effective temperature. The result of the hot F star HD 49933 and the Sun validates the second methodology. Unfortunately, the ensemble analysis of 24 Kepler targets by Kiefer et al. (2017) (see Fig. 45) was unable to distinguish between both scaling relations due to the large uncertainties in the frequency shifts. However the authors made a comment about a weak correlation with T eff supporting the scaling relation by Metcalfe et al. (2007) but with one clear exception, KIC 8006161 the high metallic solar analogue already discussed in detail in previous paragraphs. Another small trend was found between the frequency shifts and age (excepting again KIC 8006161) in the direction predicted by both scaling relations, i.e., reduction of the frequency shifts (activity) with age as found by Skumanich (1972). Fig. 45 Average frequency shifts as a function of effective temperature for 24 Kepler targets. Black diamonds, red triangles, orange circles and purple asterisks represent stars younger than 4 Gyr, between 4 to 5 Gyr, between 5 to 6 Gyr, and between 6 to 7 Gyr respectively. Figure

Relation between frequency shifts and amplitude shifts for a large sample of stars
The analysis of the temporal variations of the seismic parameters and S ph of a larger sample of 87 solar-type Kepler stars (Santos et al., 2018b), showed similar results to those of Kiefer et al. (2017). They found that about 60% of the stars in the sample show "(quasi-)periodic variations" in the frequency shifts. Moreover, 20% of the stars show frequency and amplitude shifts correlated instead of anticorrelated. Although these results seem to be puzzling, they could be explained in a simple way. First, Salabert et al. (2018) showed that small variations in frequency shiftsthat could be interpreted as due to magnetic origin -can be explained by different realizations of stochastic noise. Second, the presence of hysteresis effects between different magnetic proxies implies that at several stages of the cycle, two indexes could be in phase while they would be most of the time in anti-phase. Longer time series covering several stellar magnetic cycles would be required to properly understand all of these observations.

On the variation of the frequency shifts with frequency
The linear dependence of the frequency shifts with frequency found in the Sun and in HD 49933 was also found in the young solar analogue KIC 10644253 (Salabert et al., 2016b). Therefore, in these three stars the perturbation inducing the variation of the mode parameters needs to be located outside of the resonant mode cavity of the modes, i.e., in a thin layer very close to the photosphere. Otherwise, an oscillatory signal would be expected as was first discussed by Goldreich et al. (1991). They explained the apparent oscillatory signal superimposed on the linear dependence of the frequency shifts with frequency in the Sun as the consequence of a perturbation located near the He I ionization zone. This depth was found by analyzing the periodicity of this oscillatory signal. Later, Gough (1994) proposed that the perturbation was much deeper and that it was due to changes in the acoustic glitch of the He II ionization layer. Indeed several authors have found variations in the amplitude of the depression in the adiabatic index, Γ 1 , at this ionization layer that could be the consequence of the changing activity on the equation of state of the gas in that layer (Basu and Mandel, 2004;Verner et al., 2006). Salabert et al. (2018) found that in four stars of the Kepler sample (KIC 5184732, KIC 8006161, KIC 8379927, and KIC 11081729), the frequency shifts normalized by the mode inertia show an oscillatory behavior instead of a linear one (see Fig. 46). Before deriving firm conclusions about the positions and mechanisms responsible for the frequency shifts in stars, it is necessary to better determine these frequency shifts as a function of frequency (with smaller uncertainties) as well as in a larger number of stars. However, it seems that the picture is much more complicated than was outlined from the analysis of the Sun.

Conclusions and Perspectives
In this work we have reviewed the theory behind the asteroseismic techniques applied to study main-sequence solar-like dwarfs, as well as the latest results obtained from ground-based and space-borne missions such as CoRoT, Kepler, K2 and TESS. This is a young research topic that has reached its constant pace during the last decade.
It has been shown that asteroseismology of MS solar-like stars is providing strong constraints about the structure and dynamics from the surface to the core of stars. It is possible today to infer masses, radii and ages with a precision and accuracy never reached before. This is impacting many fields in stellar astronomy. For example, precise and accurate ages of field stars at all stages of evolution in the main sequence have shown that stars seem to stop braking when they reach a given Rossby number and from there, they seem to continue with a quite constant surface rotation . This could be the consequence of a change in the properties of external magnetism, which could also impact stellar magnetic cycles (Metcalfe and van Saders, 2017) and the possibility of giving a correct age to middle aged stars through gyrochronology. An extended sample of stars and complementary ground-based observations of the magnetic activity will be necessary to progress in this area.
Asteroseismology is also helping to improve exoplanet research (for a full review on the synergies between asteroseismology and exoplanetary science see Huber, 2018). For example, precise stellar ages are a key parameter to date the full planetary systems and thus, better understand the theory of formation and evolution of planet-hosting stars and the extrasolar planet systems as a whole. The asteroseismic improvement in the determination of stellar radius is directly impacting the precision of planet radius extracted with the transit method. Thanks to this new increased precision, it is now possible to properly characterize the so-called "radius valley" or "photoevaporation desert" at around 2 R ⊕ (Lundkvist et al., 2016;Van Eylen et al., 2018a). Finally we mention that the combination of transit photometry with asteroseismology allows a systematic measurement of orbital eccentricities of transiting planets (e.g. Van Eylen and Albrecht, 2015), which was only possible before in relatively large gas-giant planets, or for multiplanet systems where the effects of eccentricities and masses could be successfully distinguished, New internal rotation profiles have encouraged stellar astrophysicists to study angular momentum transport and efficient mixing processes and develop new mechanisms explaining the quasi-uniform rotation found in the outer part of the radiative zone, the external convective zone and the surface.
In conclusion, asteroseismology of solar-type stars is in very good shape. At the time of writing these conclusions, the community is actively analyzing K2 data, where dozens of new pulsators are foreseen, as well as the first two sectors of TESS data that are already available. New missions will contribute to enhance our known sample of pulsating star. The future is already here because of the work engaged to prepare the ESA's M3 PLATO mission, which will be able to characterize tens of thousand of these MS cool dwarfs after 2026 for which many of them will be planet hosts. Asteroseismology of MS solar-like dwarfs is just at the dawn of its potential.

Acknowledgements
The authors wish to thank the entire SoHO, CoRoT and Kepler teams, without whom many of the results presented in this review would not be possible. The authors received funding from the European Community seventh programme ([FP7/2007([FP7/ -2013) under grant agreement no. 312844 (SPACEINN) and under grant agreement no. 269194 (IRSES/ASK). The authors also acknowledges funding from the CNES. RAG acknowledges the ANR (Agence Nationale de la Recherche, France) program IDEE (n. ANR-12-BS05-0008) "Interaction DesÉtoiles et des Exoplanètes". The authors also want to thank Dr. J. Leibacher who did a thorough lecture of the manuscript helping us to improve the manuscript.

Appendix A
In table 1, the values of the large frequency spacing and the effective temperature used to plot Fig. 22 are listed, as well as the references where these values are from. A complete version of this table is available on line.