Tidal excitation of auto-resonant oscillations in stars with close-by planets

Context. Close-by planets can excite various kinds of oscillations in their host stars through their time-varying tidal potential. Aims. Magnetostrophic oscillations with a frequency much smaller than the stellar rotation frequency have recently been proposed to account for the spin–orbit commensurability observed in several planet-hosting stars. In principle, they can be resonantly excited in an isolated slender magnetic ﬂux tube by a Fourier component of the time-varying tidal potential with a very low frequency in the reference frame rotating with the host. However, due to the weakness of such high-order tidal components, a mechanism is required to lock the oscillations in phase with the forcing for long time intervals (10 3 –10 7 yr) in order to allow the oscillation amplitude to grow. Methods. We propose that the locking mechanism is an auto-resonance produced by the nonlinear dependence of the oscillation frequency on its amplitude. We suggest that the angular momentum loss rate is remarkably reduced in hosts entering auto-resonance and this contributes to maintaining those systems in that regime for a long time. Results. We apply our model to a sample of ten systems that show spin–orbit commensurability and estimate the maximum drifts of the relevant tidal potential frequencies that allow them to enter the auto-resonant regime. Such drifts are compared with the drifts expected from the tidal evolution of the planetary orbits and the stellar angular momentum loss in the magnetized winds, and we ﬁnd that auto-resonance is a viable mechanism in eight systems, at least in our idealized model. Conclusions.


Introduction
Most of the extrasolar planets known to date orbit close to their host stars because the most effective discovery techniques, the method of the transits and that of the radial velocity, favor the detection of close-by planets. In particular, the probability of transiting a star of radius R s for a planet on a circular orbit of radius a is R s /a, while the amplitude of the radial-velocity wobble of the star is proportional to a −1/2 (e.g., Cochran & Hatzes 1996). As a consequence of the small star-planet separations, these systems are characterized by remarkable tidal interactions (Mathis et al. 2013;Ogilvie 2014).
In a reference frame rotating with the star, the tidal potential due to a close-by planet is time independent only if the orbital plane coincides with the equatorial plane of the star, the orbit is circular, and the stellar rotation period is equal to the orbital period. In most of the known systems, at least one of these conditions is not verified and the time-dependent potential can excite waves in the fluid body of the star, which is called the "dynamical tide." The Kepler space telescope has revealed several close stellar binary systems with eccentric orbits where dynamical tides produce multi-period oscillations of the optical flux. They have been called heartbeat stars (e.g., Fuller 2017). The phenomenology produced by close-by planets is much less extreme because their masses are ∼ 10 −5 −10 −3 times those of their host stars; nevertheless, the possibility of exciting inertial waves and/or gravity waves in the hosts has been recognized and ana-lyzed in several works (e.g., Ogilvie & Lin 2007;Barker 2020;Ahuir et al. 2021). Friedland (2009) and references therein. An application to celestial mechanics is presented by Friedland (2001).
In the present work, the auto-resonant physical system is introduced in Sect. 2 together with a description of its relevant nonlinearities. An application to the star-planet systems in close spin-orbit commensurability considered by Lanza (2022) is presented in Sect. 3, and some consequences and the general relevance of the excited auto-resonant oscillations in the host stars are discussed in Sect. 4.

Tidal potential in star-planet systems
The gravitational tidal potential, Ψ, produced by a planet inside its host star is a solution of the Laplace equation ∇ 2 Ψ = 0. Considering a reference frame with the origin in the barycenter of the star, the polar axis along its spin axis, and rotating with the angular velocity Ω s of the star, Ψ can be expressed as (e.g., Ogilvie 2014, Sect. 2.1) where {z} is the real part of a complex quantity z and the components of the Fourier series of the potential are given by where j = √ −1, G is the gravitation constant, m p the mass of the planet, a the semimajor axis of the relative orbit, A lmn (e, i) a coefficient depending on the eccentricity of the orbit e and its obliquity i (that is, the angle between the stellar spin and the orbital angular momentum; see Kaula 1961;Mathis & Le Poncin-Lafitte 2009), r the distance from the center of the star of radius R s , Y m l (θ, φ) the spherical harmonic of degree l and azimuthal order m (with |m| ≤ l) that is a function of the colatitude θ measured from the polar (spin) axis and the azimuthal coordinate φ, ω mn ≡ (mΩ s − nΩ 0 ), with Ω s ≡ 2π/P rot the rotation frequency of the star having a rotation period P rot and Ω 0 = 2π/P orb the orbital mean motion, where P orb is the orbital period of the planet; n is an integer specifying the time harmonic of the potential oscillation, and t is the time. The series development given in Eq. (2) is valid in the interior of the star, that is, for r ≤ R s .
The frequency of the magnetostrophic waves, ω, that we consider in this work is much smaller than the stellar spin frequency Ω s . They can resonate with a component of the tidal potential, the frequency of which |ω mn | = |mΩ s − nΩ 0 | Ω s , which implies P rot /P orb m/n. In other words, a close commensurability between the stellar rotation period and the orbital period of the planet is required for a tidal excitation of those waves. A sample of stars where such a commensurability between P rot and P orb has been observed is introduced and discussed by Lanza (2022).
Owing to the weakness of the resonant tidal potential component, Ψ lmn , a long time interval is required to amplify the magnetostrophic oscillations up to a significant amplitude implying that the magnetic structure where they are excited must be very long lived (see Sect. 2.2).

Resonant slender magnetic flux tubes
The interior structure of late-type stars consists of an inner radiative zone surrounded by a convective envelope. Magnetic fields are subject to a buoyancy force in the super-adiabatically stratified layers of the convection zones, while they can reside for long times in the sub-adiabatic radiative zones because magnetic buoyancy becomes relevant only when the magnetized plasma comes in thermal equilibrium with the surrounding medium. The boundary layer between the radiative zone and the overlying convective envelope is characterized by some penetration of the convective motions into the stably stratified layers because the downward directed convective columns do not stop where the Schwarzschild criterion is no longer verified, but continue their motion below the convection zone due to their inertia. Such a boundary layer is characterized by a slightly sub-adiabatic stratification and is called the overshoot layer by analogy with the case of the convective cores of upper main-sequence stars because of such undershooting convective motions. A description of the physical properties of the layer and of the convective downdrafts was originally given by Zahn (1991) and Rieutord & Zahn (1995).
The overshoot layer allows the storage of strong azimuthal magnetic fields thanks to its sub-adiabatic stratification parameterized by δ ≡ ∇ − ∇ ad < 0, where ∇ is the gradient in the layer and ∇ ad is the adiabatic gradient (see Kippenhahn et al. 2012). The strong gradient of the rotation angular velocity observed in the Sun by means of helioseismic techniques (Schou et al. 1998), makes the overshoot layer apt to amplify toroidal fields so that it has been considered a possible seat for the solar and stellar hydromagnetic dynamos (e.g., Parker 1993). Unfortunately, its physical properties are not precisely known because of the lack of a detailed description of the undershooting convection, the strong shear, and the possible molecular weight gradient there present (e.g., Spiegel & Zahn 1992).
A detailed treatment of the azimuthal magnetic fields stored inside the overshoot layer is hampered by this lack of information, so that it is useful to model the azimuthal magnetic field as consisting of several slender azimuthal flux tubes separated by an unmagnetized plasma. For simplicity, we consider one of such axisymmetric flux tubes located in the equatorial plane of the star, so that its magnetic field is B = Be φ , where e φ is the unit vector in the azimuthal direction and B is the intensity of the magnetic field (see Fig. 1).
We adopt the slender flux tube approximation that assumes that the field is uniform across the section of the flux tube and that its section radius is much smaller than the curvature radius, assumed equal to the radius r 0 at the base of the convection zone, and the local pressure scale height, H. Such assumptions make the equations governing the motion of the slender flux tube particularly simple as shown by, for example, Ferriz-Mas & Schüssler (1993) and Ferriz-Mas & Schuessler (1994).
In our model, resonant magnetostrophic oscillations are excited in a slender toroidal flux tube and are slowly amplified along a timescale ranging from ∼ 10 3 to ∼ 10 7 years. This implies that the flux tube must remain stably stored inside the overshoot layer for those long time intervals; therefore, it is distinct from the flux tubes that emerge along the stellar activity cycle giving rise to active regions in the photosphere (Caligari et al. 1995). The emergence of such slender flux tubes is produced by an undulatory instability that requires their magnetic field B to exceed a threshold that depends mainly on the sub-adiabatic gradient δ, the radius r 0 , the acceleration of gravity g, and the stellar angular velocity Ω s (e.g., Spruit & van Ballegooijen 1982;van Ballegooijen 1983). It is about 10 5 G in the Sun and can reach several times 10 6 G in rapidly rotating late-type stars with a deep convection zone (cf. Granzer et al. 2000 Fig. 1. Sketch of a slender magnetic flux tube in equilibrium in the equatorial plane of a late-type main-sequence star (solid green line). The xy plane is the equatorial plane of the star as seen from the north pole. The barycenter of the star is O, its radius is R s , and r 0 is the radius of the unperturbed flux tube. The convective zone of the star is rendered in yellow, and the radiative zone, including the overshoot layer where the flux tube is located, is rendered in white. The flux tube perturbed by an oscillation mode with azimuthal wavenumber m = 8 is plotted as a solid black line.
The magnetic field strength of the resonant flux tubes considered in our model is below the threshold for the onset of the undulatory instability being between 10 G and a few times 10 4 G in the systems investigated by Lanza (2022). Therefore, such resonant flux tubes are stable to the undulatory instability and can be stored for long time intervals inside the sub-adiabatic overshoot layer. We assume that the field inside the flux tube is strong enough to overcome the effects of the convective downdrafts that tend to break it apart. If the velocity of the downdrafts, d , is on the order of 10 m s −1 in the layer where the tube is stored (cf. Zahn 1991), this requires that the tube field B > ∼ d √ µρ, where µ is the magnetic permeability of the plasma and ρ the density in the layer. In the case of the Sun, this implies a field on the order of 100 − 500 G, which is much lower than the undulatory instability threshold. The overshoot layer can be the seat of a mild turbulence characterized by motions mainly in the horizontal direction because the stable stratification strongly hampers motions in the vertical (radial) direction. This can give rise to an anisotropic turbulent viscosity that can play a fundamental role in avoiding the spreading of the radial-velocity gradient into the radiative zone (Spiegel & Zahn 1992;Brun & Zahn 2006). Nevertheless, by assuming that the magnetic field inside a resonant magnetic flux tube is larger than the equipartition field with the downdrafts and the associated turbulent motions, we can assume molecular values for the magnetic and thermal diffusivities across the flux tube and consider the molecular kinematic viscosity for the damping of the wave motions excited inside the flux tube itself (Lanza 2022).
In view of such considerations, we assume the molecular values adopted by Brun & Zahn (2006) for the top of the solar radiative zone, that is, a thermal diffusivity κ T = 10 3 m 2 s −1 , a magnetic diffusivity η M ∼ 0.1 m 2 s −1 , and a kinematic viscosity ν ∼ 3 × 10 −3 m 2 s −1 . With a reference diameter of the resonant flux tube of d f = 5 × 10 6 m, such values imply that the magnetic diffusion timescale of the flux tube is d 2 f /η M ∼ 2.5 × 10 14 s, that is, on the order of ∼ 10 7 years. In other words, the flux tube can have a very long lifetime, a necessary condition for the slow auto-resonant wave excitation to operate. On the other hand, the thermal diffusion timescale is shorter by a factor of 10 4 , which implies that the resonant flux tube must be assumed in thermal equilibrium with its surroundings. The corresponding magnetic buoyancy can be counterbalanced by the aerodynamic drag force produced by the convective downdrafts, thus allowing the flux tube to be stably stored inside the lower part of the overshoot layer (see Sect. 3.2 of Lanza 2022, for more details).
We stress that the considered resonant flux tubes do not take part in the stellar activity cycles having periods of a few decades. They should be regarded as the low-field tail of the magnetic flux tube distribution produced by the hydromagnetic stellar dynamo. They can be pumped into the lower part of the overshoot layer by the convective downdrafts and do not emerge to the surface, but are subject to a very slow decay under the action of their internal molecular magnetic diffusivity. The processes producing slender flux tubes in the overshoot layers of the Sun and late-type stars are beyond the scope of this work and have been discussed, for example, by Schüssler (2005) and Fan (2021).

Tidally forced oscillations in stars
Our idealized model for the magnetic field configuration in the overshoot layer allows us to consider a simple approach to study the tidal excitation of resonant oscillations in our model star. A more realistic model should consider the configuration of the internal magnetic field and the effects it has on the oscillation modes that propagate across the whole star adopting an approach similar to that of, for example, Schenk et al. (2001), Lai & Wu (2006), Fuller (2017), or Lai (2021. This approach shows that the tidal excitation of a stellar oscillation mode requires not only a close match between the mode frequency and the frequency of a Fourier component of the tidal potential, but also a significant spatial overlap between the mode oscillation and the gradient of the exciting potential component as quantified by the so-called tidal overlap coefficient where ρ(r) is the density of the plasma at position r, ξ α (r) the mode spatial eigenfunction with α an index identifying the mode and the asterisk indicating complex conjugation, while r l Y m l gives the spatial dependence of the potential of the exciting tidal component specified by the degree l and the azimuthal number m; the integral is extended over the volume V of the star (Lai & Wu 2006).
In our simplified approach, the mode eigenfunctions coincide with the eigenfunctions of the slender flux tube oscillations. They are different from zero only inside the flux tube centered at r = r 0 in the equatorial plane of the star and have a spatial dependence proportional to exp( jmφ) that overlaps with the azimuthal dependence of the tidal potential component with the same azimuthal number m (cf. sect. 2.4, Eq. 11). In other words, our highly idealized model allows an easy excitation of a slender flux tube mode provided that its frequency comes close to one of the frequencies of the tidal potential components and both the oscillation and the exciting tidal component have the same azimuthal wavenumber m to warrant their spatial overlap. This is possible because we consider a single slender flux tube in isolation inside the star. Moreover, we assume the magnetic field and the density to be constant across the flux tube section to avoid any phase mixing that could damp the excited oscillations. Small gradients of the field or of the density can be tolerated if the flux tube magnetic field has a small twist as usually assumed in dynamical models to keep flux tubes coherent during their motion (cf. Fan 2021).
In a real star the situation can be rather different if several slender flux tubes can interact with each other. In this case, the tidal potential can excite oscillations only in the resonant flux tube, while no oscillations are excited in the neighbor flux tubes owing to the sharpness of the resonance (see Sect. 2.4). Once the amplitude of the oscillations in the resonant tube becomes significant, its oscillations are strongly damped by magnetic reconnection with the unperturbed fields of the neighbor flux tubes whose oscillations are not excited. Therefore, our model can work only if the resonant flux tube is sufficiently isolated from neighbor flux tubes as to avoid any dynamical and magnetic interactions with them. Fortunately, such a condition is verified if the separation between neighbor flux tubes is at least on the order of 10 6 m, as we shall see in Sect. 2.8, where we find that the strongest coupling is due to the turbulent viscosity in the overshoot layer. On the other hand, magnetic diffusivity can be neglected in our approximation because the relatively strong fields of the flux tubes damp any turbulence inside them, making the magnetic diffusivity equal to its molecular value, that is, much smaller than the turbulent viscosity.

Forced oscillations in a magnetic flux tube
The equations describing the small free oscillations of a toroidal magnetic flux tube around its equilibrium configuration were derived by Ferriz-Mas & Schüssler (1993) and Ferriz-Mas & Schuessler (1994) who applied them to study the stability of the flux tube itself. A useful property consists in the decoupling of the oscillations in the equatorial plane from those in the latitudinal direction, which allows us to simplify our analysis by considering only oscillation modes in the equatorial plane. We generalize Ferriz-Mas & Schüssler's nondimensional equations of motion to the case of a flux tube subject to an external tidal force −ρ 0i ∇Ψ per unit volume: where ρ 0i is the internal density in the unperturbed flux tube, a dot indicates the time derivative, ξ = (ξ φ , ξ r ) is the small Lagrangian displacement in the equatorial plane with respect to the position of equilibrium; φ 0 the azimuthal coordinate measured along the unperturbed flux tube; τ a characteristic timescale defined as with β ≡ 2µp i0 /B 2 , p i0 being the unperturbed internal pressure in the flux tube, H the pressure scale height, g 0 the acceleration of gravity at the tube radius r 0 ; c T is the tube wave speed with c s = γp i0 /ρ i0 being the internal sound speed and γ = 5/3 the adiabatic exponent; f ≡ H/r 0 ; x ≡ r 0 Ω 2 s /g 0 ; x e ≡ r 0 Ω 2 e0 /g 0 ; and T is given in the limit β 1 by where σ ≡ d log g(r)/d log r r 0 expresses the dependence of the acceleration of gravity on the radial distance and q ≡ r 0 Ω e (r 0 )/[2Ω e (r 0 )] measures the shear due to the radial differential rotation at r = r 0 with the prime indicating the derivative with respect to the radial coordinate. The forcing produced by the mode of the tidal potential with frequency ω mn can be written as where the complex quantities (p φ ,p r ) can be deduced from the expression of the potential in Eq.
(2), while the tidal frequency ω mn = mΩ s − nΩ 0 is real. We consider only the stationary solution of the system of equations, assuming that the effects of the initial conditions have had time to decay under the action of the viscous dissipation. Alternatively, the same solution is obtained if the oscillations are excited from an initial state of vanishing amplitude.
To solve Eqs. (4) and (5), we consider an individual mode of oscillation of the form whereξ is a constant vector, t the time, and the periodicity of the solution in the azimuthal direction imposes that m be an integer (see Fig. 1). Substituting Eqs. (11), (9), and (10) into Eqs. (4) and (5), it is possible to reduce the latter equations to an algebraic system with ω = ω mn and the same m of the component of the tidal potential (cf. Ferriz-Mas & Schuessler 1994, Eqs. (7a) and (7b), for the corresponding freely oscillating system): whereω mn = τ ω mn ,Ω s = τ Ω s , The matrix of the algebraic system (12) is Hermitian. Therefore, we can diagonalize it and decouple the oscillations into two normal modes orthogonal to each other, each satisfying the equation where i = 1, 2 and P i are the components of the vector of the forcing terms (9) and (10) after applying the coordinate transformation that leads to the normal modes (see Eqs. 17 and 18).
In the present application, we focus on forced oscillations having a frequency much smaller than the stellar rotation frequency in the so-called magnetostrophic regime that is characterized by a flux tube magnetic field, B, such that the Alfven frequency is also much smaller than the rotation frequency. Following Ferriz-Mas & Schuessler (1994), such a condition can be written as A (m/r 0 ) 2Ω s , where A = B/ √ µρ 0i is the Alfven velocity along the magnetic flux tube, m the azimuthal wavenumber of the oscillation mode and r 0 the radius of the flux tube. As shown in Sect. 4 of Lanza (2022), such an approximation is well verified in the systems we investigate.
For a magnetic field intensity on the order of a few hundred gauss, as considered in our model, β ∼ 10 9 so that |T | ∼ β|δ| ∼ 10 3 when we assume |δ| ∼ 10 −6 , considered to be a typical value for the overshoot layer of the Sun or solar-like stars (cf. Sect. 4 of Lanza 2022, or Table 3). This implies that |A| |T | in the limit ω mn Ω s , f ∼ 0.1, and m ≤ 10. The eigenvalues Λ i of the matrix in Eq. (12), each associated with the normal mode Θ i , verify |Λ 1 | |Λ 2 | ∼ |T |, so that the amplitudes of the normal oscillation modesΘ i becomeŝ where the final equalities imply a suitable normalization of the eigenvectors and we have made use of A/D ∼ ω mn /Ω s 1 in the magnetostrophic regime. We explore the validity of these approximations from a quantitative viewpoint in Sect. 3 for the eight systems investigated in the present study (cf. Table 4 for the values of the coefficients appearing in Eqs. 17 and 18).
The physical motivation for the almost orthogonal oscillations of ξ φ and ξ r as expressed by Eqs. (17) and (18) is the stable stratification of the overshoot layer that hampers any radial motion, making the Coriolis force term 2Ω sξr in Eq. (4) negligible. The stabilizing effect of such a stratification is quantified by the large |T |. Therefore, the oscillations occur mainly in the azimuthal direction, as already shown by the solutions of the equations of motion for the oscillations found by Lanza (2022) that implyξ φ ∼ |T |ξ r .
In view of the above considerations, we make the approximation Θ 2 ∼ ξ φ , focusing on the oscillations of the toroidal flux tube in the azimuthal direction and writing its equation of motion in the form where 2ζ = τ −1 d is the inverse of the viscous damping timescale given by τ d = d 2 f /ν with d f the diameter of the flux tube and ν the kinematic viscosity of the plasma inside it as introduced in Sect. 2.2, while ω is the frequency of the magnetostrophic mode. This frequency is given by the dispersion relation obtained by equating the determinant of the matrix in Eq. (12) to zero because the damping term is so small that it has no practical effect on the eigenfrequencies of the free oscillations. The presence of the small damping term avoids a singularity when the frequency of the forcing tidal potential becomes equal to the eigenfrequency of the magnetostrophic oscillations (its expression is given in Eq. 29 of Lanza 2022, and is not repeated here).
In principle, the thermal and Ohmic damping of the oscillations are also relevant because the molecular values of the thermal and magnetic diffusivities verify the inequality κ T η M ν in a stellar radiative zone. Nevertheless, the thermal and Ohmic diffusion timescales are on the order of ∼ 10 4 − 10 7 years (see Sect. 2.2), that is, they are comparable or longer than the typical growth timescales of the resonant oscillations that range between 10 3 and 10 5 years (see Sect. 3). Therefore, we can neglect thermal and Ohmic damping of the oscillations in Eq. (19) and consider only the smaller viscous damping because it gives the narrowest resonance width, making their excitation the most difficult to achieve when the tidal frequency varies along with the evolution of the stellar spin and the orbital semimajor axis of the system. In the linear regime, we can neglect the dependence of the frequency of the magnetostrophic waves on their amplitude. However, such a dependence plays a fundamental role in the auto-resonance process; therefore, we go beyond the linear regime in the next subsection.

Nonlinear dependence of the oscillation frequency
The frequency ω of the magnetostrophic oscillations in Eq. (19) depends on the magnetic field B inside the oscillating flux tube. (Lanza 2022) showed that it is a monotonously increasing function of B and found that it depends on the amplitude of the oscillations through the perturbation of the mean magnetic field along the flux tube induced by the oscillations themselves. Specifically, he considered only the perturbation of the mean magnetic field induced by the radial displacement, ξ r (see Sect. 3.9 and, in particular, Eqs. (45) and (46) of Lanza 2022) and found that it is capable of keeping the flux tube in resonance with the tidal potential despite small variations in its frequency, ω mn , owing to tidal changes in the orbit of the planet and slow variations in the stellar rotation rate.
In addition to the nonlinear dependence of the oscillation frequency found by Lanza (2022), two other dependences can be pointed out. The first depends both on the radial displacement ξ r and the radial differential rotation and is introduced in Appendix A. There we show that it is negligible in comparison to the second nonlinearity, which depends on the azimuthal displacement ξ φ ξ r , and plays the most relevant role in the present model. Therefore, we introduce it in the following.
The component of the momentum equation of the plasma along the instantaneous axisl of the flux tube, that is, along its magnetic field lines, can be written as where is the velocity of the plasma inside the flux tube with l = · e l its component in the axial direction of unit vector e l , t the time, ρ the density, p the pressure, and Φ G the gravitational potential. The component of the Lorentz force along the flux tube axis is zero; therefore, it does not appear in the right-hand side of Eq. (20). Given that the deformation of the flux tube is small (ξ r ξ φ ), the deviation of the tube axis e l from the unperturbed azimuthal direction e φ is a second-order effect when we regard the oscillation velocity field as a firstorder quantity in the hypothesis of small oscillations. Therefore, we approximate l φ , where φ = · e φ is the azimuthal velocity component, that is, the component along the axis of the unperturbed flux tube e φ . In the same way, the derivative of any physical quantity along the direction e l is approximated by its derivative in the azimuthal direction. In other words, we assume ∂/∂l ∼ r −1 0 ∂/∂φ 0 , where φ 0 is the azimuthal angle along the unperturbed flux tube. In particular, the stellar gravitational po-tential Φ G can be assumed to be axisymmetric, thus allowing us to drop the term ∂Φ G /∂l ∼ r −1 0 ∂Φ G /∂φ 0 in Eq. (20). Equation (20) can be put in a simpler form considering the identity ∇( 2 /2) = ( · ∇) + × (∇ × ) and that the perturbations produced by the magnetostrophic oscillations are adiabatic because the thermal adjustment timescale of the flux tube is much longer than the oscillation period (cf. Sect. 2.2). This condition allows us to neglect the density perturbation δρ produced by the wave in comparison with the perturbation of the pressure δp because it is smaller by a factor δρ/δp ∼ 2 /c 2 s 1, where is the wave velocity and c s the sound speed in the layer.
By applying the above results, we can recast Eq. (20) in the form where the azimuthal component of the velocity is where 0 is the uniform and constant azimuthal velocity that enters into the equation of the mechanical equilibrium of the unperturbed flux tube because the angular velocity of rotation of the plasma inside the flux tube is in general slightly different from that of the surrounding plasma (see Moreno-Insertis et al. 1992;Lanza 2022, Eq. 5). Substituting Eq. (22) into Eq. (21) and taking into account that ∂/∂φ 0 → jm, we can integrate the resulting equation with respect to φ 0 and obtain where f (t) is a function of the time that plays the role of the integration constant that is independent of φ 0 and we have approximated φ because the radial component of the velocity is much smaller than the azimuthal component. Considering the nodes of the wave, that is, the points where its velocity is zero at any time t, we have where p 0 is the plasma pressure in the nodes that we assume equal to the uniform pressure inside the flux tube in equilibrium and in the absence of the wave. In this way, we have that f (t) = F 0 is a constant independent of the time t and the azimuthal coordinate φ 0 .
Making an average of Eq. (23) over one oscillation period of the wave, the first-order terms disappear and the average perturbation of the pressure inside the flux tube is found to be where δp ≡ p − p 0 is the pressure perturbation produced by the wave, while the angle brackets indicate the time average over one period of the sinusoidal oscillation. The lateral pressure balance across the flux tube is given by where p e is the constant external pressure that is independent of φ 0 . Substituting Eq.(25) into Eq. (26), we find the average variation of the magnetic field along the slender flux tube when a magnetostrophic wave of frequency ω is propagating along it: In the regime of magnetic field strength characteristic of our resonant flux tubes, the frequency of the magnetostrophic waves is linearly dependent on the average magnetic field, that is, Therefore, as a consequence of Eq. (27), the frequency of the wave varies nonlinearly with its amplitude as where and ω 0 is the frequency of the wave for a vanishing amplitude. Neglecting the effect of the stratification on the frequency of the magnetostrophic mode, ω 0 ≈ A m/r 0 , which gives the approximate relationship γ 0 ≈ (1/4)(m/r 0 ) 2 , independent of the magnetic field strength. The first-order perturbation analysis applied to derive Eq. (28) is valid only when γ 0ξ 2 φ 1, which is equivalent to (ξ φ /r 0 ) 2 1 because of the above approximation.

Auto-resonance
Equation (19) of the magnetostrophic oscillations, taking into account the nonlinear dependence of their frequency on their amplitude found in Sect. 2.5 (cf. Eq. 28), can be written as where γ D = −2γ 0 and we have considered that γ 0ξ 2 φ 1. Equation (30) is the equation of a forced Duffing oscillator with a small frictional damping (e.g., Wawrzynski 2021), that is, one of the nonlinear systems that shows auto-resonance (Friedland 2009).
The dependence of the oscillation frequency on the amplitude shown by Eq. (30) allows it to adjust itself in order to remain in resonance with the external forcing, even when this is varying its frequency, provided that such a variation is sufficiently slow (see below). In such a way, the amplitude of the oscillations can continuously grow by extracting energy from the external forcing, even when its strength is very small because the resonance condition can be maintained for a time interval much longer than the oscillation period.
A quantitative model of auto-resonance can be formulated following different approaches developed to study nonlinear oscillators (e.g., Chirikov 1979;Friedland 2001). We follow the multiple timescale method as detailed in Appendix B. Defining a nondimensional time t ≡ ω 0 t, the forced oscillator equation can be recast as where λ = ζω −1 0 and the forcing term in the right hand side of Eq. 30, depending on the tidal potential component Ψ lmn , has been written as Article number, page 6 of 21 where η > 0 is the small amplitude of the oscillations of the potential the frequency of which is assumed to be slowly varying as where κ ω 0 is the inverse of the timescale of variation of ω mn = mΩ s −nΩ 0 due to the slow evolution of the orbit of the planet and the braking of the stellar rotation under the action of the stellar wind and the tides. The multiple timescale method seeks for a solution of the forced oscillator equation of the form where A is the slowly varying amplitude and φ A is the initial phase of the oscillation. The lag angle between the forcing and the oscillation can be defined as ψ ≡ θ − ϕ + π, where θ ≡ t + φ A is the phase of the oscillation. The oscillator is in resonance when ψ = π/2 implying that the forcing makes work to increase the amplitude of the oscillation during all the oscillation cycle. When the system is not too far away from resonance, ψ is a slowly varying function of the time such as the amplitude A.
The multiple timescale method allows us to derive the equations describing the slow variations in A and ψ as (cf. Appendix B) This system of ordinary differential equations can be numerically integrated to find the amplitude and the lag angle as functions of the time. Auto-resonance corresponds to a regime where the angle ψ is fixed or oscillating in an interval centered on ψ = π/2 that corresponds to the resonance condition. In such a regime, the variations of A and ψ are characterized by oscillations that keep the system always in phase with the external forcing, thus allowing the mean amplitude of the oscillations to grow. For given values of the parameters γ D and η, the autoresonant regime is possible only if |κ| is smaller than a threshold value that depends on those parameters (Friedland 2009). To find the maximum value of |κ| compatible with autoresonance, we followed the approach in Sect. 3 of Friedland (2001). We assume that the system crosses the linear resonance keeping a continuous phase locking with an almost constant angle ψ 1, so that cos ψ ∼ 1 can be assumed to be constant in Eq. (36). We differentiate this equation with respect to t and substitute Eq. (35) neglecting the very small dissipation term because λ 1. In such a way, we obtain Equation (37) describes the motion of a quasi-particle moving in a tilted cosine potential: In other words, the time variation in the angle ψ is described by the unidimensional motion of a point mass ruled by the Lagrangian L = (1/2)ψ 2 − V qp (ψ). The system is in auto-resonance if and only if the angle ψ can oscillate around a fixed point ψ 0 that corresponds to a minimum of the potential V qp . This implies that auto-resonance is possible only if the potential has a relative minimum where dV qp /dψ = 0, in other words, otherwise, the potential V qp has no relative minima and increases monotonically, leading to a monotone variation of ψ in time. The function S (A) has a minimum at the point where dS /dA = 0, that is, at A = A min = [(4/3)(η/γ D )] 1/3 . The corresponding minimum value of S is S min ≡ S (A min ) and it can be substituted into Eq. (41) to obtain the condition to be verified by κ for the system to enter into an auto-resonant regime, that is, where the threshold value κ th is given by The value of κ th given by Eq. (43) is an order-of-magnitude estimate because of the simplifying assumption adopted to derive Eq. 37, but it is adequate for our purposes in view of the comparable or larger uncertainties in the current tidal and magnetic braking effects in planet-hosting stars (cf. Sect. 2.7). External perturbations acting on the flux tube may produce a phase shift of its oscillations and destroy the resonance locking. An order of magnitude estimate based on the energy integral of the equation of motion (37) requires that the maximum value of the nondimensional time derivative of the phase lag verifieṡ ψ max < ∼ (2k th ) 1/2 to keep the oscillations of ψ confined around a minimum of V qp . This translates into a condition on the characteristic timescale of any external perturbation t pert that must be t pert > ∼ (2k th ) −1/2 ω −1 0 , that is, the perturbation must be sufficiently slow to allow the system to maintain the phase locking. Considering typical values to be k th ∼ 10 −7 and ω 0 ∼ 5 × 10 −8 s −1 (see Sect. 3), we find t pert > ∼ 1400 years, which is a very long timescale in comparison to those of the known perturbations that can act on the flux tube. Therefore, the effects of such perturbations cannot be neglected and should be included as an additional term on the right-hand side of Eq. (36). Nevertheless, the oscillation period of the lag angle ψ in the absence of external perturbations ranges between 10 2 and 10 5 years as we shall see in Sect. 3. The known perturbations are periodic with a much shorter period, for example, the deformation due to the l = 2 component of the tidal potential has a frequency 2(Ω 0 − Ω s ), or are random with typical timescales of ∼ 10 6 − 10 7 s as in the case of the convective downdrafts. Therefore, we can average the additional term accounting for those external perturbations in Eq. (36) over a timescale remarkably longer than the perturbation timescales, yet much shorter than the period of the unperturbed oscillations of ψ, making its contribution vanish. In other words, the large difference between the timescale of the ψ oscillations and those of the known perturbations acting on the flux tube makes it possible to neglect the effects of the latter on the auto-resonant excitation of the former.
The asymptotic increase of the mean amplitude in autoresonance can be derived by considering that the term proportional to A −1 in the right-hand side of Eq. (36) becomes negligible in the asymptotic regime, while dψ/dt ∼ 0 on the left-hand side because the angle ψ makes small oscillations around a minimum of the pseudo-potential. Therefore, we obtain the asymptotic scaling A 2 (t ) ∼ (8/3)κt /γ D . Considering two time instants t 0 and t , this gives where A 0 ≡ A(t 0 ) is the value of A at the initial time t 0 when the auto-resonant asymptotic regime can be assumed to begin. Formally, the asymptotic solution given by Eq. (44) is valid until the system reaches its maximum amplitudeξ φ max when the external forcing is balanced by dissipation, that is,ξ φ max = A max ∼ η sin ψ/(2λ) ≤ η/(2λ) (cf. Eq. 35). It corresponds to a great amplification of the forcing amplitude because λ −1 = ζ −1 ω 0 is large, thus allowing the small tidal potential oscillations to produce a sizeable effect as already discussed by Lanza (2022).
In the high-Reynolds-number environment of the overshoot layer, the oscillations are likely to produce turbulence that limits their amplitude in the asymptotic regime. Inside the flux tube, the magnetic field can effectively dump any turbulent motion; therefore, turbulence will not directly limit the ξ φ displacement, but the radial displacement ξ r is affected by the turbulence in the surrounding medium produced by the oscillations themselves. Such a turbulence can effectively damp ξ r , if its amplitude grows too large. The maximum amplitude achievable in that case has been estimated in Sect. 3.7, Eq. (39), of Lanza (2022) and can be expressed in our formulation asξ r max ∼ ω −2 mn (∂Ψ lmn /∂r). The maximum azimuthal amplitudeξ φ max is linked to the maximum radial amplitude byξ φ max ∼ |T |ξ r max because of the equations of motion governing the flux tube oscillations (cf. Eqs. 4 and 5 and Sect. 2.4), so that turbulence in the surrounding medium will limit also the azimuthal displacement.
In some cases, the value ofξ φ max obtained from the above conditions is so large that the small-displacement approximation |γ D | ξ 2 φ 1 of our first-order theory is no longer valid. When this happens, we conservatively assume a maximumξ φ max = 0.1|γ D | −1/2 in such a way that the small-displacement approximation remains valid with |γ D | ξ 2 φ ≤ 0.01. In summary, the value ofξ φ max is computed as the minimum between i) the equilibrium amplitude as given by Eq. (35), that is, η/(2λ); ii) the asymptotic amplitude as given by Eq. (44) with κ = κ th and t = 10 7 years corresponding to the diffusion timescale of the resonant flux tube; iii) the maximum amplitude in the presence of self-generated turbulence; and iv) the requirement of small oscillations, that is,ξ φ max ≤ 0.1|γ D | −1/2 . The value ofξ φ max determined in this way has to be regarded as an order-of-magnitude estimate in view of the uncertainties on κ th , λ, and the diffusion timescale of the resonant flux tubes, or the possibility of adopting a less strict smallness condition on |γ D |ξ 2 φ .

Frequency variation in the exciting tidal potential ω mn
The auto-resonance regime can be accessed only if the drift of the excitation frequency |κ| is sufficiently small as shown by Eq. (42). The sharpness of the resonance peak due to the smallness of the dissipation parameter λ (see Lanza 2022), implies that ω 0 is virtually equal to ω mn at resonance. In other words, κ, as defined in Eq. (33), is given by The variation in the frequency of the forcing tidal potential ω mn , responsible for the excitation of the oscillations, can be written as where the variation in the angular velocity of the stellar rotation is produced by the action of the stellar magnetized wind and the tides raised by the planet, while the variation in the orbital frequency (mean motion) is due to the tides. Specifically, adopting the tidal model of Mardling & Lin (2002), we can write where I s is the moment of inertia of the star, dL W /dt the angular momentum loss rate due to the stellar wind, Q s the stellar modified tidal quality factor, G the gravitation constant, m p the mass of the planet, a the orbit semimajor axis, and R s the star radius. The moment of inertia is scaled from the solar value because the planet hosts we consider have an internal stratification similar to that of the Sun, that is, I s = 0.07M s R 2 s , where M s is the mass of the star. When Ω s is a function of the position inside the star, the stellar angular velocity entering into the specification of the tidal torque is the angular velocity of the convection zone because the tidal dissipation in the radiative interior is remarkably smaller due to the smaller radius of the radiative zone (cf. Ogilvie & Lin 2007;Barker 2020).
The angular momentum loss rate due to the stellar wind that appears in Eq. (47) can be parameterized as where K W = 2.7 × 10 40 kg m 2 s −1 is a constant calibrated with the present Sun, Ω sat is the saturation angular velocity that is a function of the mass of the star (Bouvier et al. 1997). We assume Ω sat = 5, 9, 14, 16 Ω for a star of mass M s = 0.5, 0.8, 1.0, 1.2 M , respectively, and linearly interpolate for masses in between. Stars with close-by giant planets may experience a reduced angular momentum loss rate because a hot Jupiter could reduce the Alfven radius of the stellar wind and the mass loss rate by modifying the large-scale configuration of the coronal field. The decrease in the angular momentum loss rate may reach two orders of magnitude with respect to Eq. (48), depending on the parameters of the system (Cohen et al. 2010;Lanza 2010Lanza , 2014. The modified tidal quality factor of the star, Q s , depends on the processes responsible for the dissipation of the tidal kinetic energy inside the star. When the rotation period is shorter than twice the orbital period, tidal inertial waves can be excited in a main-sequence star, leading to a remarkable increase of the tidal dissipation that corresponds to a decrease of Q s with respect to the regime where P rot > 2P orb in which tidal dissipation is produced only by the equilibrium tide (Ogilvie & Lin 2007). Following Collier Cameron & Jardine (2018), we assume log 10 Q s = 7.3 when inertial waves can be excited and log 10 Q s = 8.26 when only the equilibrium tide is present.
The variation in the orbital mean motion, Ω 0 , due to the tides, is given by where the symbols have already been introduced above. Neglecting the tides and the possible reduction in the Alfven radius due to a close-by massive planet, the rate of angular momentum loss in the stellar wind would in several cases produce a value of |κ| larger than the maximum allowed by Eq. (42) when the frequency ω mn crosses the resonance. The discrepancy is by at least one or two orders of magnitude in several systems (cf. Sect. 3). Therefore, we must invoke the effects of tides raised by the planet or effects not included in the above simplified model to account for the initial locking of those star-planet systems into auto-resonance.
A necessary condition for tides to counteract the magnetic braking is that Ω 0 > Ω s because only in such a case they transfer angular momentum from the orbit to the stellar spin counteracting the angular momentum loss produced by the stellar wind. Conversely, when Ω s > Ω 0 tides contribute to brake the star by transferring angular momentum from the stellar spin to the orbit (cf. Eq. 47). Moreover, tidal effects can be significant only when the planet has a mass at least comparable with that of Jupiter and a relative orbital separation a/R s < ∼ 10. When a/R s > ∼ 10, the tidal torque on the stellar spin is too small, even for a very massive hot Jupiter (m p ∼ 10 M J ) and assuming a reduced Alfven radius, because the torque scales as (a/R s ) −5 (cf. Eq. 49) and the minimum value we can adopt for Q s ∼ 10 6 for the late-type main-sequence hosts under consideration (e.g., Ogilvie & Lin 2007;Barker 2020).
In the regime where tides due to the planet are not efficient, the variation in ω mn is dominated by the evolution of the stellar spin. Therefore, the possibility of making |κ| smaller than the threshold value for triggering auto-resonance could be related to a temporary stalling of the evolution of the stellar spin. Since the flux tube oscillations predicted by our model require a timescale shorter than 10 6 − 10 7 years for their amplification, a stalled rotation phase with a duration much shorter than the typical timescale of stellar spin evolution (10 8 −10 9 years) is sufficient to allow the system to enter into auto-resonance. Once established, the auto-resonant regime can maintain itself for long time intervals, in some cases comparable with the main-sequence lifetime of the system itself, as discussed by Lanza (2022) and as we shall show in Sect. 2.10.
Observations of stellar rotation in open clusters in the age range ≈ 0.7 − 2.7 Gyr suggest that a stalling of the rotational evolution can indeed occur for main-sequence stars with M s < ∼ 0.8 M with its beginning and duration principally depending on the stellar mass (e.g., Agüeros et al. 2018;Curtis et al. 2019;Douglas et al. 2019). Specifically, while the rotational evolution of stars with a mass of ∼ 0.8 M can be stalled for ≈ 0.3 Gyr, their lower mass siblings of ≈ 0.55 M appear to be stalled for ∼ 1.3 Gyr (Curtis et al. 2020). Stars with a mass higher than 0.8 M could also experience a phase of stalled rotation, but it may have gone unnoticed because of a duration on the order of 30 − 100 Myr, that is, shorter than the shortest separation in age between the clusters analyzed to detect the phenomenon. Spada & Lanzafame (2020) interpreted the stalled rotation phase by invoking an exchange of angular momentum between the radiative interior and the convective envelope in late-type main-sequence stars. Their two-zone model predicts a phase when the resurfacing of angular momentum stored in the radiative interior can counteract the braking of the convection zone by the stellar wind, temporarily stalling the rotational evolution of a star. During such a phase, the value of |κ| may become small enough to lock the system into auto-resonance. This will lead to a significant amplification of the flux tube oscillations in the systems for which one of the tidal potential frequencies crosses the oscillation frequency of one of the slender flux tubes in the overshoot layer. Given that the magnetic field intensity of the potentially resonant flux tubes has a rather wide distribution, the possibility for one of the frequencies of the tidal potential components to come into resonance with one of the flux tubes is not a remote one. The observed duration of the phase of stalled rotation, ranging between ∼ 0.3 and 1.3 Gyr depending on the star mass, is certainly much longer than the amplification timescale of the auto-resonance process, which makes it even more likely to be triggered (cf. Sect. 3).

Wind torque in auto-resonant systems
Once the auto-resonant oscillations are excited and grow to a significant amplitude, the azimuthal velocity field of the wave inside the flux tube propagates in its immediate neighborhood. The propagation of the oscillating azimuthal flow is ruled by the Navier-Stokes equation because the medium around the flux tube is unmagnetized. The radial component of the oscillation velocity field is on the order of |T | −1 ∼ 10 −3 of the azimuthal component and is in general smaller than the velocity field of the convective downdrafts that is on the order of 10 − 100 m s −1 in the mid of the overshoot layer (Zahn 1991;Rieutord & Zahn 1995). The downdrafts are braked and their flow is turned horizontal close to the base of the overshoot layer. Their velocity field is a random field with a coherence timescale on the order of 10 6 − 10 7 s, corresponding to the convective turnover time close to the base of the overlying convection zone, and a spatial coherence scale on the order of 10 7 m or smaller (cf. Eq. 13 and Sect. 3.2 in Pinçon et al. 2021). On the other hand, the azimuthal component of the oscillation velocity field has typical coherence timescales and length scales one order of magnitude larger than those of the downdraft field and is periodic (cf. Tables 3 and 5 in Sect. 3). Therefore, it will dominate the development of the undulatory instability in the neighbor flux tubes.
Specifically, the coherently oscillating velocity field of an auto-resonant flux tube may trigger the development of the undulatory instability in neighbor flux tubes with a field strength ≥ 10 5 G that are regarded to be responsible for the formation of active regions at the surface of the star (Caligari et al. 1995). The azimuthal wavenumber of the most unstable mode is m = 1 or m = 2 (Ferriz-Mas & Schuessler 1994; Granzer et al. 2000); therefore, the growth of the instability can be initiated close to one of the nodes of the azimuthal velocity field of a neighbor resonant flux tube because its azimuthal velocity field propagates to the flux tube whose magnetic field is close to the threshold for the development of the undulatory instability. The amplitude of the perturbation does not need to be large because the unstable state of the strong-field flux tube will exponentially amplify it with a characteristic growth timescale on the order of 10 7 s (Caligari et al. 1995;Schüssler 2005;Fan 2021). This allows the resonant flux tube to be still considered as virtually isolated from the strong-field flux tube and its oscillations are not damped by magnetic reconnection as discussed in Sect. 2.4 also thanks to their very small radial displacement ξ r ∼ |T | −1 ξ φ ∼ 10 −3 ξ φ .
The propagation of the azimuthal velocity field into the fluid around an oscillating flux tube is due to the small viscosity present in the overshoot layer. Assuming that the propagation length scale is small in comparison with the radius r 0 and the diameter of the cross section of the oscillating flux tube, d ft , the Navier-Stokes equation can be approximated as where ν is the kinematical viscosity in the overshoot layer. We seek for a solution of the kind φ (r, φ, t) ∝ ξ φ (φ, t)F(r), assuming a radial dependence F(r) ∝ exp[−(1 + j)λ F r], where λ F is the inverse length scale ruling the radial propagation of the flow in the neighborhood of the oscillating flux tube (cf. Batchelor 1968, Sect. 4.3). By substituting such an ansatz into Eq. (50), we find the inverse propagation length scale: Considering a frequency of the tidal potential in the rotating reference frame of the star ω mn ∼ 10 −4 Ω s ∼ 10 −10 s −1 and a molecular viscosity ν ∼ 3 × 10 −3 m 2 s −1 , the propagation length scale λ −1 F is smaller than 10 4 m. Nevertheless, considering that the viscosity in the radiative zone should be larger by a factor of at least 10 4 than the molecular value to account for the rigid rotation of the solar radiative zone (e.g., Ruediger & Kitchatinov 1996;Spada et al. 2010), such a length scale may increase up to 10 6 m, even without invoking the turbulence produced by the oscillations of the resonant flux tube once these reach a sizeble amplitude. If the magnetic field in the overshoot layer has a spatially intermittent structure with the resonant flux tube close to the strong-field flux tubes with a mean separation comparable with or smaller than λ −1 F , the ensemble of active regions formed in the photosphere will reproduce the same pattern of the auto-resonant oscillation mode. Specifically, different strongfield flux tubes may become unstable with their m = 1 or m = 2 growing oscillations in spatial coincidence with the nodes of the oscillations of the auto-resonant flux tube. The pattern of photospheric active regions will statistically display the same azimuthal wavenumber of the oscillating flux tube. The longitude of the starspots observed in the hot Jupiter host Kepler-17 may be accounted for by such a mechanism (cf. Sect. 2 of Lanza 2022, for a description of the observations).
The potential magnetic field in the stellar corona corresponding to a photospheric pattern dominated by a mode with azimuthal order m will have a degree l ≥ m. Since m can be as large as 8 (cf. Lanza 2022), the corresponding increase in the degree l of the coronal field will strongly reduce the Alfven radius of the stellar wind and consequently the angular momentum loss rate (Réville et al. 2015). Specifically, considering the fitting formula provided by Réville et al. (2015) and based on numerical simulations, the torque for given stellar radius, angular velocity, and mass loss rate scales with an exponent approximately of 2/(2l + 2.7). For solar-like parameters, their models find a decrease by a factor of ∼ 4 in passing from the dipolar (l = 1) to the quadrupolar (l = 2) configuration. Going from the dipolar to the configuration with l = 8, their analytic formula gives a reduction by a factor ∼ 32. Considering that the actual coronal field configuration consists of the contributions of components with different values of l, we prudently adopt a reduction by a factor of 10 of the stellar wind torque when auto-resonant oscillations are excited. The corresponding slower evolution of the stellar spin allows a system to enter and remain locked into an auto-resonant regime for long time intervals as already envisaged in Sect. 3.9.2 of Lanza (2022) and as we shall see in Sect. 2.10.
Finally, we note that the nonlinear development of the instability of the strong-field flux tube produces a perturbation on the resonant tube itself having a duration significantly shorter than the period of its oscillations because the unstable flux tube emerge to the photosphere in approximately one month (Caligari et al. 1995). In other words, the perturbation effects are averaged to zero on the timescales characteristic of the amplification of the oscillations and thus have virtually no effect on the resonant flux tube and the auto-resonant process (cf. Sects. 2.6 and 3).

Tidal dissipation and torque in auto-resonant systems
The equations governing the oscillations of the auto-resonant flux tubes are Eqs. (4) and (5). Following the method introduced in Sect. 3.4 of Lanza (2022), it is possible to write their first integral, that is, the equation of conservation of the mechanical energy, in the form where E, the total mechanical energy density of the oscillating flux tube, consists of the kinetic energy, the magnetic tension energy, and the potential energy coming from the work done against the stable stratification of the overshoot layer by the radial displacement of the oscillations (see Lanza 2022, for details). The power dissipated per unit volume is given by the dissipation function because ξ r ξ φ in resonance. The total dissipated power is obtained by integrating p diss over the volume of the oscillating flux tube. Since its density can be regarded as uniform and its volume is V ft = (1/2)π 2 r 0 d 2 ft , where d ft is the diameter of the cross section of the flux tube, we find the total dissipated power, averaged over one period of the oscillations, where we substitutedξ φ = ω mn ξ φ and considered that the average of ξ 2 φ over one oscillation period isξ 2 φ /2. Equation (54) gives the dissipated power in the reference frame rotating with the star due to the oscillations excited by the tidal potential (cf. Sect. 2.2 of Ogilvie 2014). The corresponding tidal torque, T , acting on the stellar spin is given by (cf. Sect. 2.2 of Ogilvie 2014) The maximum tidal torque corresponds to the maximum amplitude of the oscillations determined according to the method illustrated in the final paragraph of Sect. 2.6, that is, T max = T (ξ φ max ). The corresponding timescale for the variation in the stellar angular velocity is

Evolution of the tidal potential frequency ω mn
In principle, the evolution of the tidal potential frequency ω mn as given by Eq. (46) can be computed if the effects of the stellar wind and the tides on the stellar spin and the orbital mean motion can be evaluated. These require the knowledge of the angular momentum loss rate and the modified stellar tidal quality factor Q s together with a numerical integration in time of the differential equations (47) and (49) for each particular system. Such a detailed investigation is beyond the scope of the present work that focuses on the excitation of the oscillations in the resonant magnetic flux tubes and will be postponed to future analyses. Nevertheless, some consideration of the evolution of the tidal frequency ω mn is in order here, especially in view of a preliminary estimation of the duration of the auto-resonant phase in the observed systems. For this reason, a theory for the evolution of ω mn is developed in this section.
To apply the results of the previous sections to the evolution of the tidal potential frequency ω mn , we write the equation giving the evolution of the orbital angular momentum as and recast Eqs. (47) as where which is valid for m p M s with r 2 g = 0.07 being the square of the nondimensional gyration radius of the host star. From Kepler third law and Eq. (57), we have where m red = M s m p /(M s + m p ) ∼ m p is the reduced mass of the system. In such a way, Eq. (46) can be recast as where In the systems were tides are negligible, Ω 0 is constant and we have and the possibility of having |κ| smaller than the threshold for the excitation of auto-resonant oscillations is linked to a phase of very slow variation in the stellar spin, in other words, a phase of stalled rotational evolution. Once auto-resonant oscillations are triggered and reach a sizeable amplitude, they can slow down the rotational evolution by reducing the coronal Alfven radius, as proposed in Sect. 2.8. This contributes to increase the duration of the phase of stalled rotation because the reduced loss rate of angular momentum from the convection zone makes the reservoir of angular momentum in the radiative zone sufficient to balance the losses for a longer time. By decreasing dL w /dt by a factor of, say, 10, the duration of the phase of stalled rotation may be increased by a corresponding factor reaching 1 − 10 Gyr according to the mass of the host star, that is, comparable with the mainsequence lifetime of the system. Therefore, we expect a long duration of the commensurability phase in such a regime, up to several gigayears, especially for lower main-sequence hosts for which the stalled rotation phase is longer.
On the other hand, in systems where the tidal torque is nonnegligible, the angular momentum loss in the wind may be balanced by the torque itself, making dω mn /dt ∼ 0 and |κ| smaller than the threshold for triggering auto-resonance. Considering that dL w /dt < 0 and K/I s > 0 in Eq. (61), when Ω 0 > Ω s this implies H(m, n) > 0, or H(m, n) < 0 when Ω 0 < Ω s . We show examples of star-planet systems that verify both kinds of prescriptions in Sect. 3.2.
By writing ω mn = n(Ω s − Ω 0 ) + pΩ, where p ≡ m − n is a positive or negative integer, we can recast Eq. (61) as where with |Θ| being the inverse timescale for the evolution of ω mn and L s ≡ I s Ω s the angular momentum of the stellar spin. The general solution of the nonhomogeneous differential equation (64) can be found by the method of the variation of the constant and is where C 0 is a constant that is determined by the initial condition. The frequency of the auto-resonant magnetostrophic waves considered in our model is much smaller than the stellar spin frequency (|ω mn | Ω s ); therefore, the initial condition can be well approximated as ω mn (0) 0 that allows us to find C 0 . The indefinite integral appearing in Eq. (66) can be evaluated only when the evolution of the stellar spin, of the tidal torque, and the angular momentum loss rate are known. Nevertheless, if the evolution of Ω s is very slow because of the balance between the tidal torque and the angular momentum loss rate of the wind (with Ω 0 > Ω s ) and |dL w /dt| is made small because of the reduced Alfven radius in the auto-resonant regime (cf. Sect. 2.8) or the presence of a massive close-by planet (cf. Sect. 2.7), the integrand in Eq. (66) can be regarded as approximately constant and the solution becomes where τ spin ≡ −L s /(dL w /dt) > 0 is the timescale for the evolution of the stellar spin under the action of the reduced angular momentum loss in the stellar wind without considering the action of tides. Restricting our analysis to systems with Ω 0 > Ω s and assuming that the tidal torque is balancing the angular momentum loss in the wind, we have H(m, n) > 0 that implies Θ > 0 according to Eq. (65). When the system is close to commensurability, |ω mn | Ω s and the second term dominates in the right-hand side of Eq. (64). In such a case, the condition for the system to enter the auto-resonant regime as given by Eq. (42) becomes which is satisfied for p > 0 provided that In some systems, pΘ ≤ κ th ω 2 mn /Ω s . In those cases, τ spin max is not defined and we have the only condition τ spin ≥ τ spin min .
The asymptotic value of ω mn , that is, its value when t Θ −1 , is given by [p − m/(Θτ spin )]Ω s and is very small if τ spin is close to m/(pΘ), which falls in the middle of the interval of the allowed values as given by Eq. (69) and, in any case, is larger than τ spin min . In such a case, ω mn remains close to its initial zero value all along the system evolution. In other words, if the host star has τ spin ∼ m/(pΘ), the system remains close to commensurability for a very long time interval comparable with its main-sequence lifetime. This is more likely to occur if pΘ/ω mn κ th (ω mn /Ω s ). When condition (69) is violated in a system with Ω 0 > Ω s , for example, if τ spin > τ spin max , the wind braking will become too weak and the system may exit the auto-resonant regime because of the tidal effects. As a consequence, the angular momentum loss rate by the wind will increase because a lower-l coronal field will become dominant (cf. Sect. 2.8). This will restore condition (69) on τ spin bringing the system back into auto-resonance. Similarly, if for some reason τ spin < τ spin min , the system will exit the auto-resonant regime, the wind braking will strongly increase and the stellar rotation will be slowed down increasing the difference Ω 0 − Ω s , that is, the tidal torque that tends to bring back the system into auto-resonance.
We assume that the value of τ spin may be adjusted in order to verify condition (69) because it depends on the amplitude of the auto-resonant oscillations that are responsible for the perturbation of the unstable flux tubes in the overshoot layer (cf. Sect. 2.8). When the oscillation amplitude increases, the perturbation will be greater, making the coronal field more strongly dominated by the high-l mode that makes the wind braking less efficient. In such a way, the value of τ spin can be increased or decreased according to the amplitude of the excited oscillations.
In conclusion, we foresee a self-regulating mechanism that can keep systems with Ω 0 > Ω s and m > n in auto-resonance with comparable tidal and wind torques for long time intervals, comparable with the main-sequence lifetimes of their host stars.

Applications
The theory introduced in Sect. 2 will now be applied to ten starplanet systems that show a commensurability between the star spin period and the orbital period of their closest planet already investigated by Lanza (2022). The purpose of that work was to account for the observed P rot /P orb = m/n commensurability as the result of the tidal excitation of magnetostrophic oscillations inside a toroidal flux tube stored in the overshoot layer of the host star. The results of the present investigation allow us to study in a quantitative way the excitation of those oscillations and the maintainance of each system in auto-resonance for a sufficiently long time interval as explained in Sect. 2.
The relevant parameters of the considered systems are listed in Table 1 that is extracted from Table 1 of Lanza (2022) with an update of the mass of AU Mic b. In Table 1, we list from the left to the right, the name of the system, the mass, M s , and the radius, R s , of the host star, its effective temperature, T eff , the orbit semimajor axis, a, and its ratio to the host radius, a/R s , the eccentricity of the orbit, e, the mass of the planet, m p , the skyprojected spin orbit angle, λ so , the star rotation period, P rot , the orbital period, P orb , the integers m and n the ratio of which best approximate the P rot /P orb ratio, and the relevant references.
Following the same approach as in Lanza (2022), for the sake of simplicity we neglect the (small) eccentricity of the planetary orbits in the computation of the components Ψ lmn of the tidal potential and keep only the effects of the obliquity i that we assume to be 20 • for all the systems, except for Kepler-13A (i = 25 • ), Kepler-63 (i = −120 • ), HAT-P-11 (i = 70 • ), and WASP-107 (i = 40 • ). Their specific obliquity values are derived from the measurements of their sky-projected spin-orbit obliquity λ so in Table 1 together with an estimate of the inclination of their spin axis to the line of sight based on the stellar rotation period and the inclination of the orbital plane given by the transit best fit. The value of the degree l of the tidal potential component is constrained to be such that l ≥ |m| and l − n be even; otherwise, the corresponding Ψ lmn component has a zero A lmn (0, i) coefficient and does not contribute to the excitation of the oscillations (cf. Ogilvie 2014). In Table 2, we list, from the left to the right, the name of the system, the fractionary radius at the base of the convection zone, r 0 /R s , the density in that layer, ρ, the minimum l compatible with these constraints together with the corresponding values of m, n, and the absolute value of the tidal coefficient, A lmn , computed as in Sect. 3.1 of Lanza (2022).
To apply our theory, we assume that the magnetic field of the resonant flux tube in the overshoot layer is in equipartition with a plasma velocity field of amplitude d = 10 m s −1 . The corresponding field intensity is B eq = d √ µρ. This ensure that the flux tube can neither be disrupted nor significantly distorted in the lower part of the overshoot layer where we assume that the velocity of the downdrafts has been braked well below 10 m s −1 . The equilibrium azimuthal velocity 0 inside the resonant flux tube is assumed to be equal to the Alfven velocity (cf. Sect. 2.5), that is, the maximum value compatible with stability against the Kelvin-Helmholtz instability (Ferriz-Mas & Schüssler 1993). However, the effect of the internal velocity 0 on the flux tube parameters is very small because it is much lower than the rotation velocity, that is, 0 = A Ω s r 0 . The parameters governing the flux tube oscillations are reported in Table 3 where we list, from the left to the right, the name of the system, the equipartition magnetic field, B eq , of the oscillating flux tube, the parameter β ≡ 2µp/B 2 eq , the nondimensional stratification parameter, T , as given by Eq. (8) of Ferriz-Mas & Schuessler (1994) (its negative values indicate a stable stratification of the overshoot layer, that is, δ < 0), the characteristic timescale τ = √ 2H/ A with H being the pressure scale height at r = r 0 (see Sect. 3 of Ferriz-Mas & Schuessler 1994), the frequency, ω 0 , of the magnetostrophic oscillations, its ratio to the spin frequency, ω 0 /Ω s , the coefficient γ D as appearing in our Duffing oscillator equation (30), the threshold frequency drift, κ th , as given by Eq. (43), the frequency drift expected on the basis of the present parameters of the system, κ exp , as computed from Eqs. (45), (46), (47), (48), and (49), and the parameter giving the asymptotic amplitude A versus the time at the threshold drift as A ∼ {[8κ th ω 0 /(3|γ D |)]t} 1/2 . To compute the parameters listed in Table 3, we made use of the parameters listed in Table 3 of Lanza (2022) and assumed a superadiabaticity δ = −2.0 × 10 −6 , an adiabatic exponent γ = 5/3, and a differential rotation parameter q = 0.06 (cf. Ferriz-Mas & Schuessler 1994).
As a preliminary step to the application of our model, we check that the normal modesΘ 1 andΘ 2 verify the approximations in Eqs. (17) and (18) by listing in Table 4 the values of the ratios (Λ i − C)/D, with i = 1, 2. We see thatΘ 2 ξ φ is an acceptable approximation for all our investigated systems.
The parameters describing the evolution of the frequency of the tidal potential ω mn in resonance with the flux tube magne-    System CoRoT-2 0.726 3.29e+02 9 8 3 5.0489e-03 CoRoT-4 0.805 2.26e+01 3 1 1 4.0068e-01 CoRoT-6 0.744 1.86e+02 7 5 7 1.0028e-01 Kepler-13A 0.120 3.31e+04 5 3 5 1.3753e-01 Kepler-17 0.795 4.66e+01 9 8 1 3.2176e-04 Kepler-63 0.726 3.29e+02 5 3 5 1.0420e-02 HAT-P-11 0.675 7.78e+02 7 6 1 8.3935e-02 τ Boo 0.868 8.65e-01 3 1 1 4.0068e-01 WASP-107 0.673 1.79e+03 3 3 1 2.0495e-01 AU Mic 0.303 6.36e+03 7 4 7 3.5363e-02 tostrophic frequency ω 0 are listed in Table 5 where we report, from the left to the right, the name of the system, H(m, n), as defined by Eq. (62), the modified tidal quality factor, Q s eq , corresponding to the equilibrium between the tidal torque and the wind braking torque in the nonsynchronous systems where Ω 0 > Ω s and H(m, n) > 0, the maximum amplitude of the azimuthal oscillations,ξ φ max , estimated with the method in the final paragraph of Sect. 2.6, the tidal timescale, τ tide , associated with the tidal torque due to the auto-resonant oscillations as given by Eq. (56) and 2ζ = 1.5 × 10 −15 s −1 , the braking timescale of the stellar wind, τ red AMLR , in the case of a magnetic braking reduced by a factor of 10 with respect to Eq. (48), the timescale Θ −1 as derived from Eq. (65), the minimum and maximum spin timescales, τ spin min , and τ spin max , as derived from Eqs. (70) and 71, respectively, and the spin timescale τ spin ne = m/(pΘ) giving dω mn /dt ∼ 0 for ω mn Ω s . The missing Q s eq values in Table 5 correspond to synchronous systems, systems with Ω 0 < Ω s , or without significant tidal interactions, that is, a/R s > 10 in our sample. The missing values of τ spin max or τ spin ne correspond to systems where such quantities are not defined because the system is synchronous, or pΘ < κ th ω 2 mn /Ω s , or pΘ < 0. The tidal timescales due to the auto-resonant oscillations, τ tide , are always much longer than the main-sequence lifetimes of the systems because of the oscillations of the phase lag angle ψ that reduce the average tidal dissipation inside the stars (see Figs. 2-7, lower panels). On the other hand, previous estimates of the same tidal timescale for CoRoT-4, τ Boo, and WASP-107 in Table 4 of Lanza (2022) were shorter than the main-sequence lifetimes of these stars because in that study a constant maximal value sin ψ = 1 was assumed. The confinement of the tidal oscillations to the resonant slender flux tube makes the tidal overlap coefficient as defined in Eq. (3) necessarily small, that is, we confirm the general weakness of the tidal star-planet angular momentum exchanges in our systems also from the point of view of the more general model cited in Sect. 2.3. Now we discuss the different kinds of systems by classifying them according to our theory.
The future evolution of the CoRoT-4 and τ Boo systems is ruled by the loss of angular momentum in their stellar winds. Considering a reduced angular momentum loss rate in the autoresonant regime, both systems can remain into the present synchronous state because τ red AMLR ∼ 2 × 10 10 years and ∼ 4.6 × 10 9 years for CoRoT-4 and τ Boo, respectively, when dL w /dt is reduced by a factor of 10 with respect to the values given by Eq. (48). Both timescales are longer than or comparable with the main-sequence lifetimes of those systems that have F-type hosts (M s > ∼ (1.2 − 1.4) M ); therefore, they are expected to maintain their synchronous state for very long time. A reduced angular momentum loss is crucial in maintaining the synchronization of τ Boo because it would be difficult to explain with a normal dL w /dt value given that the timescale for the evolution of the stellar spin would be of ∼ 5 × 10 8 years in that case.

Nonsynchronous systems with significant tidal
interactions: Kepler-13A,  In these systems, the value of a/R s < ∼ 6.5 and m p > ∼ 2.5 M J indicating a remarkable tidal interaction. CoRoT-2 and Kepler-17 are examples of systems where Ω 0 > Ω s , Θ > 0, and p ≡ m − n > 0, so that the tidal torque acting on the star can counterbalance the wind braking, while Kepler-13A is an example of a system with Ω 0 < Ω s , Θ < 0, and p < 0 where both tides and stellar wind contribute to brake stellar rotation. In Kepler-13A, the predicted value of the frequency drift κ exp computed on the basis of our model is ∼ 4 times smaller than the threshold value κ th for entering auto-resonance (cf. Table 3). Given the present uncertainties on the tidal and wind braking processes, we regard this system as capable of entering autoresonance. It can last for all the main-sequence lifetime of the system given that τ red AMLR ∼ 10 9 years. In this particular system, the host star has a mass of ∼ 1.7 M , making the angular momentum loss intrinsically weak as observed in similar early F-type stars. Assuming a tidal quality factor Q s ∼ 2 × 10 7 , in agreement with the estimate by Lanza et al. (2011b) for the Ftype star CoRoT-11, we obtain a tidal braking timescale of the stellar spin of |Θ| −1 ∼ 0.7 Gyr in Table 5. Nevertheless, CoRoT-11 has a mass of 1.2 M ; thus, Q s is probably larger for the  more massive host Kepler-13A, making such a timescale correspondingly longer. In conclusion, a slow evolution away from the present commensurability state is likely to occur in this system, independently of any decrease in the angular momentum loss rate due to the auto-resonant phenomenon. As an illustration of the auto-resonant regime in Kepler-13A, we plot in Fig. 3 the amplitude and the phase lag computed for κ = 0.03κ exp and 2ζ = 1.5 × 10 −15 s −1 to warrant a faster excitation of the oscillations.
In CoRoT-2, the estimated value of the tidal frequency drift κ exp is about 40 times larger than the threshold value κ th , while the discrepancy is only within a factor of 2 in the case of Kepler-17 (cf . Table 3). Therefore, the triggering of the auto-resonance requires a temporarily remarkable reduction of the spin evolution in CoRoT-2 and/or a balance between the tidal torque and the wind angular momentum loss in this host in order to reduce the drift of the frequency ω mn below the threshold. A balance between the tidal torque and the intense wind produced by the high level of magnetic activity of CoRoT-2 requires Q s ∼ 10 7 , assuming a reduction of the braking efficiency by a factor of 10 due to its close-by massive planet (cf. Table 5). Such a value of Q s is rather difficult to justify with the equilibrium tide only, but it cannot be excluded given the present uncertainties in the tidal theory in fast-rotating stars (e.g., Barker 2020; Wei 2022). On the other hand, Kepler-17 could have entered the auto-resonant regime with a value of Q s as large as 6 × 10 8 when we assume a reduction by a factor of 10 in the angular momentum loss rate with respect to Eq. (48) (cf. Table 5). Such a value of Q s is perfectly in line with the expectations based on the theory by Ogilvie & Lin (2007) and the results of Collier Cameron & Jardine (2018) for equilibrium tide.
The evolution of the amplitude A and the phase lag ψ of the auto-resonant oscillations are illustrated in Figs. 4 and 5 for CoRoT-2 and Kepler-17, respectively. For illustrative purposes, we adopted κ values corresponding to a reduction by a factor of 200 and 0.3 with respect to the values of κ exp listed in Table 3, respectively, 2ζ = 1.5 × 10 −15 s −1 , and do not include an upper bound on the amplitude.
Our model predicts the possibility of a long duration of the auto-resonant regime and of the P rot -P orb commensurability in both systems, if the tidal torque is capable of balancing the wind torque, due to the reduced angular momentum loss rate in the wind. While a detailed model of the rotational evolution of CoRoT-2 is beyond the scope of the present work, an application of our simplified model in Sect. 2.10 shows that a value of 3.6 × 10 9 < ∼ τ spin < ∼ 5.3 × 10 9 years may keep the system in auto-Article number, page 15 of 21 A&A proofs: manuscript no. ms43700 resonance all along its main-sequence lifetime. The corresponding reduction in the angular momentum loss rate with respect to the value computed with Eq. 48 is by a factor of ∼ 14 − 20. This result gives support to Poppenhaeger & Wolk (2014), who suggest that CoRoT-2 has an age greater than 4 − 5 Gyr, at variance with its fast rotation that indicates an age of < 0.5 Gyr by comparison with similar single stars in open clusters.
In the case of Kepler-17, τ spin > ∼ 4.5 × 10 9 years is predicted by our model. It is longer by a factor of only ∼ 1.5 than predicted by Eq. (48), implying a modest reduction in the angular momentum loss rate, maybe because Kepler-17 is rotating about three times slower than CoRoT-2. If τ spin is close to τ spin ne ∼ 1.7 × 10 10 years, the initial commensurability is maintained all along the main-sequence evolution given the long timescale Θ −1 ∼ 15 Gyr.
An important consequence of these results is that gyrochronology may not be applicable to such systems owing to their slowed spin-orbit evolution in the auto-resonance regime.

Systems with negligible tidal interactions
Systems with a negligible tidal interaction are characterized by a/R s > 10 and include CoRoT-6, Kepler-63, HAT-P-11, WASP-107, and AU Mic. Kepler-63 and WASP-107 have an expected frequency drift κ exp comparable with and much smaller than their threshold values κ th , respectively; thus, those systems can be expected to have entered the auto-resonance phase. For illustration purpose, we plot the amplitude and phase lag in auto-resonance for Kepler-63 in Fig. 6 for κ = 0.005κ exp and 2ζ = 1.5×10 −15 s −1 .
For the other systems, the estimated values of κ are larger than the thresold values by factors of ∼ 400 for CoRoT-6, ∼ 5000 for HAT-P-11, and 8 × 10 8 for AU Mic, the latter two representing cases that cannot be accounted for by the model proposed by Lanza (2022). A phase of stalled rotational evolution could have allowed CoRoT-6 to enter auto-resonance because a constant stellar rotation lasting for < ∼ 10 5 years is enough to trigger auto-resonance in such a system, as illustrated in Fig. 7   where we plot amplitude and phase lag for κ = 5 × 10 −4 κ exp and 2ζ = 1.5 × 10 −15 s −1 .
In the case of HAT-P-11 and AU Mic, the maximum amplitude of the azimuthal resonant oscillations is onlyξ φ max ∼ 150 − 300 m, that is, at least four orders of magnitude smaller than for the other systems listed in Table 5. This is a consequence of the weakness of the tidal potential in these two systems owing to a/R s > 15 and m p < 0.1 M J . Such small oscillation amplitudes are unlikely to significantly perturb the neighbor magnetic flux tubes responsible for stellar photospheric activity, thus indicating that also the present model may not explain the observed commensurability in HAT-P-11 and AU Mic. Therefore, other models should be considered to account for a potential influence of the close-by planets on the activity of these hosts (e.g., Lanza 2008;Abreu et al. 2012;Stefani et al. 2018) 1 .
The possibility that a low-frequency resonant oscillation mode of a star evolves with the orbital and stellar spin frequencies in order to maintain the commensurability condition has been proposed in stellar binaries and has been called resonance locking (see, e.g., Fuller 2017; Zanazzi & Wu 2021). Its applicability to the star-planet case remains to be investigated. In the specific case of AU Mic and HAT-P-11, this mechanism is highly unlikely because of the extremely long tidal evolution timescales in those systems, which is much longer than the evolution timescales of their stellar structures: the latter are required to be larger than the former in order to make the mechanism work (cf. Eq. 34 in Zanazzi & Wu 2021).
The duration of the auto-resonant regime in the systems considered in this subsection is predicted to be comparable with their main-sequence lifetime because of the weakness of the tidal interactions making the timescales |Θ| −1 much longer than their host star nuclear timescales. With the possible exception of WASP-107, the minimum spin evolution timescales τ spin min are also very long because of the reduction in the angular momentum loss rate in auto-resonance. Therefore, no significant evolution is predicted in their P rot /P orb ratio during their main-sequence lifetime according to our model, except for the pre-main-sequence star AU Mic that is still contracting to reach the zero-age main sequence.

Discussion and conclusions
We have introduced a model that can interpret the close spinorbit commensurability observed in several planetary systems in the framework of the approach proposed by Lanza (2022) that considers the tidal excitation of resonant oscillations in a magnetic slender flux tube stored in the overshoot layer of their host stars. Specifically, our model accounts for the excitation and the amplification of such oscillations despite the weakness of the tidal potential by considering that the system can enter an autoresonant regime. In such a regime, the flux tube oscillations can be kept in phase with the oscillations of a weak tidal potential component for a time interval that is sufficiently long as to amplify them. A consequence of the oscillations of such resonant flux tubes is the triggering of the undulatory instability in other flux tubes stored in the overshoot layer that have a stronger magnetic field; this makes them emerge into the photosphere, where they form a pattern characterized by a high-l coronal field topology. This leads to a remarkable reduction in the angular momentum loss rate, which may keep these systems in the auto-resonant regime for timescales comparable with their main-sequence lifetime.
We analyzed the same ten systems considered by Lanza (2022), finding that auto-resonance is a viable amplification mechanism for eight of them. The two exceptions are HAT-P-11 and AU Mic, for which the mechanism produces a maximum oscillation amplitude at least four orders of magnitude smaller than in the other systems.
It is important to stress that the introduced model, being based on an individual resonant magnetic flux tube in isolation, is highly idealized (cf. Sect. 2.3). The configuration of the magnetic field in the overshoot layers of late-type stars is largely unknown; thus, it is not possible to assess the relevance of this simplifying assumption in real stars. A future direction of study can be the consideration of low-frequency oscillations in a magnetized layer, as proposed by, for example, Zaqarashvili (2018), to see whether they can be tidally excited in the case of a close spin-orbit commensurability.
Most of the previous studies on the effect of magnetic fields on dynamical tides focused on the dissipation of the tidal energy that was found to be generally increased because of the added effect of Ohmic dissipation (cf. Wei 2016Wei , 2018Lin & Ogilvie 2018). On the other hand, Astoul et al. (2019) considered the role of magnetic fields in the excitation of dynamical tides and found that it can become relevant in the case of a commensurability between the stellar spin and the orbital mean motion of a planet on an oblique or eccentric orbit. Such a result, although obtained from an order-of-magnitude comparison between the Lorentz force and the hydrodynamical effects in the excitation of the dynamical tides, provides support for the proposal by Lanza (2022) and the present idealized excitation model. In other words, the possible tidal excitation of magnetostrophic oscillations in stars that host planets in spin-orbit commensurability is worthy of further investigation.
Our treatment allows us to compute the tidal dissipation inside stars in the auto-resonant regime. We find values remarkably smaller than those predicted by Lanza (2022) with the simplifying assumption of a perfect and continuous phase locking between the oscillations and the forcing tidal potential. By computing the actual phase lag between forcing and auto-resonant oscillations, we find that the average torque associated with the dissipation of the flux tube oscillations is in general negligible in comparison with those of the other dynamical tides proposed to operate in late-type planetary hosts (Ogilvie & Lin 2007;Mathis 2015;Barker 2020).
The fraction of close-by planetary systems that display a spin-orbit commensurability can be roughly estimated by considering, for example, the compilation by Albrecht et al. (2021), where the rotation periods for a total of 51 host stars are reported. Those systems were selected because they have a transiting close-by planet with a measurement of its projected orbit obliquity, a property that should not, to our knowledge, be correlated with the presence of a spin-orbit commensurability. Among those stars, we find seven of our systems, that is, ∼ 15% of their sample. Assuming that all those cases are not due to a chance coincidence, this would correspond to a duration of the commensurability regime of at least ∼ 1.5 Gyr, assuming an average main-sequence lifetime of ∼ 10 Gyr and that all the systems enter such a regime sooner or later. On the other hand, if only some fraction of the systems experience a phase of commensurability during their evolution, the duration of this phase can be longer, as suggested in Sect. 2.10. As a consequence, we suggest that gyrochronology cannot be trusted in the case of stars with closeby massive planets, notably when they show a commensurability between their rotation and orbital periods.
Similarly, we expand x 1 into a Fourier series as where A 0 is the initial value of A at T 1 = 0. Substituting into Eq. (B.24), and integrating with respect to T 1 , we find where φ A0 is the initial value of φ A . Finally, we consider the case of the damped Duffing oscillator with an external cosinusoidal forcing. Its equation of motion is (cf. Sect. 2.6) x + 2λẋ + x − γ D x 3 = η cos ϕ, (B.27) where ϕ = t − 1 2 κ t 2 (B.28) is the phase of the forcing, the frequency of which is varying on the slow timescale due to the term containing κ (see below). Equation ( The zeroth order solution can be assumed again x 0 = A(T 1 ) cos T 0 + φ A (T 1 ) = A cos θ. We introduce the angle ψ defined as Therefore, substituting this into Eq. (B.10) and considering that The left-hand side of Eq. (B.33) is ∂ϕ/∂t = 1 − κt . Comparing with the right-hand side and separating the zero-order term from the first-order term, we find where we defined κ 0 ≡ κ/ . Making use of these relationships, developing in a Fourier series h(A cos θ, −A sin θ), and grouping the terms in sin θ and cos θ in Eq. (B.32), we find Considering that ∂/∂T 1 = −1 ∂/∂T 0 , for T 1 = T 0 , and taking into account the definitions of the coefficients κ 0 , λ 0 , γ 0 , and η 0 given above, we finally find (B.40) which are Eqs. (35) and (36), respectively, if we identify the normal time T 0 with the time t as defined in Sect. 2.6.