Tidal friction is thought to be important in determining the long-term spin-orbit evolution of short-period extrasolar planetary systems. Using a simple model of the orbit-averaged effects of tidal friction, we study the evolution of close-in planets on inclined orbits, due to tides. We analyse the effects of the inclusion of stellar magnetic braking by performing a phase-plane analysis of a simplified system of equations, including the braking torque. The inclusion of magnetic braking is found to be important, and its neglect can result in a very different system history. We then present the results of numerical integrations of the tidal evolution equations, where we find that it is essential to consider coupled evolution of the orbital and rotational elements, including dissipation in both the star and planet, to accurately model the evolution. The main result of our integrations is that for typical Hot Jupiters, tidal friction aligns the stellar spin with the orbit on a similar time as it causes the orbit to decay. This means that if a planet is observed to be aligned, then it probably formed coplanar. This reinforces the importance of Rossiter-McLaughlin effect observations in determining the degree of spin-orbit alignment in transiting systems. We apply these results to the XO-3 system, and constrain the tidal quality factors Q' in both the star and planet in this system. Using a model in which inertial waves are excited by tidal forcing in the outer convective envelope and dissipated by turbulent viscosity, we calculate Q' for a range of F-star models, and find it to vary considerably within this class of stars. This means that assuming a single Q' applies to all stars is probably incorrect. We propose an explanation for the survival of WASP-12 b & OGLE-TR-56 b, in terms of weak dissipation in the star.Comment: 19 pages, 8 figures, accepted in MNRA

We study tidal dissipation in stars with masses in the range 0.1 − 1.6M⊙ throughout their evolution, including turbulent effective viscosity acting on equilibrium tides and inertial waves in convection zones, and internal gravity waves in radiation zones. We consider a range of stellar evolutionary models and incorporate the frequency-dependent effective viscosity acting on equilibrium tides based on the latest simulations. We compare the tidal flow and dissipation obtained with the conventional equilibrium tide, which is strictly invalid in convection zones, finding that the latter typically over-predicts the dissipation by a factor of 2-3. Dissipation of inertial waves is computed using a frequency-averaged formalism accounting for realistic stellar structure for the first time, and is the dominant mechanism for binary circularization and synchronization on the main sequence. Dissipation of gravity waves in the radiation zone assumes these waves to be fully damped (e.g. by wave breaking), and is the dominant mechanism for planetary orbital decay. We calculate the critical planetary mass required for wave breaking as a function of stellar mass and age, and show that this mechanism predicts destruction of many hot Jupiters but probably not Earth-mass planets on the main sequence. We apply our results to compute tidal quality factors following stellar evolution, and tidal evolutionary timescales, for the orbital decay of hot Jupiters, and the spin synchronization and circularization of binary stars. We also provide predictions for shifts in transit arrival times due to tidally-driven orbital decay of hot Jupiters that may be detected with NGTS, TESS or PLATO.

Tidally distorted rotating stars and gaseous planets are subject to a well-known linear fluid instability-the elliptical instability. It has been proposed that this instability might drive enough energy dissipation to solve the long-standing problem of the origin of tidal dissipation in stars and planets. But the nonlinear outcome of the elliptical instability has yet to be investigated in the parameter regime of interest, and the resulting turbulent energy dissipation has not yet been quantified. We do so by performing three dimensional hydrodynamical simulations of a small patch of a tidally deformed fluid planet or star subject to the elliptical instability. We show that when the tidal deformation is weak, the nonlinear outcome of the instability leads to the formation of long-lived columnar vortices aligned with the axis of rotation. These vortices shut off the elliptical instability, and the net result is insufficient energy dissipation to account for tidal dissipation. However, further work is required to account for effects neglected here, including magnetic fields, turbulent convection, and realistic boundary conditions.

We study the fate of internal gravity waves approaching the centre of an initially non‐rotating solar‐type star, by performing three‐dimensional numerical simulations using a Boussinesq‐type model. These waves are excited at the top of the radiation zone by the tidal forcing of a short‐period planet on a circular, coplanar orbit. This extends previous work done in two dimensions by Barker & Ogilvie. We first derive a linear wave solution, which is not exact in three dimensions; however, the reflection of ingoing waves from the centre is close to perfect for moderate amplitude waves. Waves with sufficient amplitude to cause isentropic overturning break, and deposit their angular momentum near the centre. This forms a critical layer, at which the angular velocity of the flow matches the orbital angular frequency of the planet. This efficiently absorbs ingoing waves, and spins up the star from the inside out, while the planet spirals into the star. We also perform numerical integrations to determine the linearized adiabatic tidal response throughout the star, in a wide range of solar‐type stellar models with masses in the range 0.5 ≤m★/M⊙≤ 1.1, throughout their main‐sequence lifetimes. The aim is to study the influence of the launching region for these waves at the top of the radiation zone in more detail, and to determine the accuracy of a semi‐analytic approximation for the tidal torque on the star, which was derived under the assumption that all ingoing wave angular momentum is absorbed in a critical layer. The main conclusion of this work is that this non‐linear mechanism of tidal dissipation could provide an explanation for the survival of all short‐period extrasolar planets observed around FGK stars, while it predicts the destruction of more massive planets. This work provides further support for the model outlined in a previous paper by Barker & Ogilvie, and makes predictions that will be tested by ongoing observational studies, such as WASP and Kepler.

We study thermal convection in a rotating fluid in order to better understand the properties of convection zones in rotating stars and planets. We first derive mixing-length theory for rapidlyrotating convection, arriving at the results of Stevenson (1979) via simple physical arguments. The theory predicts the properties of convection as a function of the imposed heat flux and rotation rate, independent of microscopic diffusivities. In particular, it predicts the mean temperature gradient; the rms velocity and temperature fluctuations; and the size of the eddies that dominate heat transport. We test all of these predictions with high resolution three-dimensional hydrodynamical simulations of Boussinesq convection in a Cartesian box. The results agree remarkably well with the theory across more than two orders of magnitude in rotation rate. For example, the temperature gradient is predicted to scale as the rotation rate to the 4/5th power at fixed flux, and the simulations yield 0.75 ± 0.06. We conclude that the mixing length theory is a solid foundation for understanding the properties of convection zones in rotating stars and planets.

We perform one of the first studies into the nonlinear evolution of tidally excited inertial waves in a uniformly rotating fluid body, exploring a simplified model of the fluid envelope of a planet (or the convective envelope of a solar-type star) subject to the gravitational tidal perturbations of an orbiting companion. Our model contains a perfectly rigid spherical core, which is surrounded by an envelope of incompressible uniform density fluid. The corresponding linear problem was studied in previous papers which this work extends into the nonlinear regime, at moderate Ekman numbers (the ratio of viscous to Coriolis accelerations). By performing high-resolution numerical simulations, using a combination of pseudo-spectral and spectral element methods, we investigate the effects of nonlinearities, which lead to time-dependence of the flow and the corresponding dissipation rate. Angular momentum is deposited non-uniformly, leading to the generation of significant differential rotation in the initially uniformly rotating fluid, i.e. the body does not evolve towards synchronism as a simple solid body rotator. This differential rotation modifies the properties of tidally excited inertial waves, changes the dissipative properties of the flow, and eventually becomes unstable to a secondary shear instability provided that the Ekman number is sufficiently small. Our main result is that the inclusion of nonlinearities eventually modifies the flow and the resulting dissipation from what linear calculations would predict, which has important implications for tidal dissipation in fluid bodies. We finally discuss some limitations of our simplified model, and propose avenues for future research to better understand the tidal evolution of rotating planets and stars.

We study the fate of internal gravity waves approaching the centre of an initially non-rotating solar-type star, primarily using two-dimensional numerical simulations based on a cylindrical model. A train of internal gravity waves is excited by tidal forcing at the interface between the convection and radiation zones of such a star. We derive a Boussinesq-type model of the central region of a star and find a non-linear wave solution that is steady in the frame rotating with the angular pattern speed of the tidal forcing. We then use spectral methods to integrate the equations numerically, with the aim of studying at what amplitude the wave is subject to instabilities. These instabilities are found to lead to wave breaking whenever the amplitude exceeds a critical value. Below this critical value, the wave reflects perfectly from the centre of the star. Wave breaking leads to mean flow acceleration, which corresponds to a spin-up of the central region of the star, and the formation of a critical layer, which acts as an absorbing barrier for subsequent ingoing waves. As these waves continue to be absorbed near the critical layer, the star is spun up from the inside out.Our results point to an important amplitude dependence of the (modified) tidal quality factor Q , since non-linear effects are responsible for dissipation at the centre of the star. If the amplitude of the tidal forcing exceeds the critical amplitude for wave breaking to occur, then this mechanism produces efficient dissipation over a continuous range of tidal frequencies. This requires (3.3, for a planet of mass m p in an orbit of period P around the current Sun, neglecting stellar rotation. However, this criterion depends strongly on the strength of the stable stratification at the centre of the star, and so it depends on stellar mass and main-sequence age. If breaking occurs, we find Q ≈ 1.5 × 10 5 (3 , for the current Sun. This varies by no more than a factor of 5 throughout the range of solar-type stars with masses between 0.5 and 1.1 M , for fixed orbital parameters. This estimate of Q is therefore quite robust and can be reasonably considered to apply to all solar-type main-sequence stars, if this mechanism operates. The strong frequency dependence of the resulting dissipation means that this effect could be very important in determining the fate of close-in giant planets around G and K stars. We predict fewer giant planets with orbital periods of less than about 2 d around such stars if they cause breaking at the centre, due to the efficiency of this process.Even if the waves are of too low amplitude to initiate breaking, radiative damping could, in principle, lead to a gradual spin-up of the stellar centre and to the formation of a critical layer. This process could provide efficient tidal dissipation in solar-type stars perturbed by less massive companions, but it may be prevented by effects that resist the development of differential rotation.

We explore the linear stability of astrophysical discs exhibiting vertical shear, which arises when there is a radial variation in the temperature or entropy. Such discs are subject to a "vertical-shear instability", which recent nonlinear simulations have shown to drive hydrodynamic activity in the MRI-stable regions of protoplanetary discs. We first revisit locally isothermal discs using the quasi-global reduced model derived by Nelson et al. (2013). This analysis is then extended to global axisymmetric perturbations in a cylindrical domain. We also derive and study a reduced model describing discs with power law radial entropy profiles ("locally polytropic discs"), which are somewhat more realistic in that they possess physical (as opposed to numerical) surfaces. In all cases the fastest growing modes have very short wavelengths and are localised at the disc surfaces (if present), where the vertical shear is maximal. An additional class of modestly growing vertically global body modes is excited, corresponding to destabilised classical inertial waves ("r-modes"). We discuss the properties of both types of modes, and stress that those that grow fastest occur on the shortest available length scales (determined either by the numerical grid or the physical viscous length). This ill-posedness makes simulations of the instability difficult to interpret. We end with some brief speculation on the nonlinear saturation and resulting angular momentum transport.