The Cassini and Juno probes have revealed large coherent cyclonic vortices in the polar regions of Saturn and Jupiter, a dramatic contrast from the east-west banded jet structure seen at lower latitudes. Debate has centered on whether the jets are shallow, or extend to greater depths in the planetary envelope. Recent experiments and observations have demonstrated the relevance of deep convection models to a successful explanation of jet structure and cyclonic coherent vortices away from the polar regions have been simulated recently including an additional stratified shallow layer. Here we present new convective models able to produce long-lived polar vortices. Using simulation parameters relevant for giant planet atmospheres we find flow regimes that are in agreement with geostrophic turbulence (GT) theory in rotating convection for the formation of large scale coherent structures via an upscale energy transfer fully three-dimensional. Our simulations generate polar characteristics qualitatively similar to those seen by Juno and Cassini: they match the structure of cyclonic vortices seen on Jupiter; or can account for the existence of a strong polar vortex extending downwards to lower latitudes with a marked spiral morphology and the hexagonal pattern seen on Saturn. Our findings indicate that these vortices can be generated deep in the planetary interior. A transition differentiating these two polar flows regimes is described, interpreted in terms of different force balances and compared with previous shallow atmospheric models which characterised polar vortex dynamics in giant planets. In addition, the heat transport properties are investigated confirming recent scaling laws obtained in the context of reduced models of GT.

Waves in thin fluid layers are important in various stellar and planetary problems. Due to rapid rotation such systems will become oblate, with a latitudinal variation in the gravitational acceleration across the surface of the object. In the case of accreting neutron stars, rapid rotation could lead to a polar radius smaller than the equatorial radius by a factor $\sim 0.8$. We investigate how the oblateness and a changing gravitational acceleration affect different hydrodynamic modes that exist in such fluid layers through analytic approximations and numerical calculations. The wave vectors of $g$-modes and Yanai modes increase for more oblate systems compared to spherical counterparts, although the impact of variations in the changing gravitational acceleration is effectively negligible. We find that for increased oblateness, Kelvin modes show less equatorial confinement and little change in their wave vector. For $r$-modes, we find that for more oblate systems the wave vector decreases. The exact manner of these changes for the $r$-modes depends on the model for the gravitational acceleration across the surface.

Accreting neutron stars (NS) can exhibit high frequency modulations in their lightcurves during thermonuclear X-ray bursts, known as burst oscillations. These frequencies can be offset from the NS spin frequency by several Hz (where known independently) and can drift by 1-3 Hz. One plausible explanation is that a wave is present in the bursting ocean, the rotating frame frequency of which is the offset. The frequency of the wave should decrease (in the rotating frame) as the burst cools hence explaining the drift. A strong candidate is a buoyant $r$-mode. To date, models that calculated the frequency of this mode taking into account the radial structure neglected relativistic effects and predicted rotating frame frequencies of $\sim$ 4 Hz and frequency drifts of > 5 Hz; too large to be consistent with observations. We present a calculation that includes frame-dragging and gravitational redshift that reduces the rotating frame frequency by up to 30 % and frequency drift by up to 20 %. Updating previous models for the ocean cooling in the aftermath of the burst to a model more representative of detailed calculations of thermonuclear X-ray bursts reduces the frequency of the mode still further. This model, combined with relativistic effects, can reduce the rotating frequency of the mode to $\sim$ 2 Hz and frequency drift to $\sim$ 2 Hz, which is closer to the observed values.

Convection in rotating spherical geometries is an important physical process in planetary and stellar systems. Using continuation methods at low Prandtl number, we find both strong equatorially asymmetric and symmetric polar nonlinear rotating waves in a model of thermal convection in thin rotating spherical shells with stress-free boundary conditions. For the symmetric waves convection is confined to high latitude in both hemispheres but is only restricted to one hemisphere close to the pole in the case of asymmetric waves. This is in contrast to what is previously known from studies in the field. These periodic flows, in which the pattern is rotating steadily in the azimuthal direction, develop a strong axisymmetric component very close to onset. Using stability analysis of periodic orbits the regions of stability are determined and the topology of the stable/unstable oscillatory flows bifurcated from the branches of rotating waves is described. By means of direct numerical simulations of these oscillatory chaotic flows, we show that these three-dimensional convective polar flows exhibit characteristics, such as force balance or mean physical properties, which are similar to flows occuring in planetary atmospheres. We show that these results may open a route to understanding unexplained features of gas giant atmospheres, in particular for the case of Jupiter. These include the observed equatorial asymmetry with a pronounced decrease at the equator (the so-called dimple), and the coherent vortices surrounding the poles recently observed by the Juno mission.

We present an extensive study of Boussinesq thermal convection including a temperature-dependent internal heating source, based on numerical three-dimensional simulations. The temperature dependence mimics triple-$\alpha$ nuclear reactions and the fluid geometry is a rotating spherical shell. These are key ingredients for the study of convective accreting neutron star oceans. A dimensionless parameter ${{\rm Ra}}_n$, measuring the relevance of nuclear heating, is defined. We explore how flow characteristics change with increasing ${{\rm Ra}}_n$ and give an astrophysical motivation. The onset of convection is investigated with respect to this parameter and periodic, quasiperiodic, chaotic flows with coherent structures, and fully turbulent flows are exhibited as ${{\rm Ra}}_n$ is varied. Several regime transitions are identified and compared with previous results on differentially heated convection. Finally, we explore (tentatively) the potential applicability of our results to the evolution of thermonuclear bursts in accreting neutron star oceans.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

Accreting neutron stars (NS) can exhibit high-frequency modulations, known as burst oscillations, in their lightcurves during thermonuclear X-ray bursts. Their frequencies can be offset from the spin frequency of the NS (known independently) by several Hz, and can drift by 1-3 Hz. One plausible explanation for this phenomenon is that a wave is present in the bursting ocean that decreases in frequency (in the rotating frame) as the burst cools. The strongest candidate is the buoyant $r$-mode; however, models for the burning ocean background used in previous studies over-predict frequency drifts by several Hz. Using new background models (which include shallow heating, and burning in the tail of the burst) the evolution of the buoyant $r$-mode is calculated. The resulting frequency drifts are smaller, in line with observations. This illustrates the importance of accounting for the detailed nuclear physics in these bursts.

This study considers the onset of stress-free Boussinesq thermal convection in rotating spherical shells with aspect ratio $\eta=r_i/r_o=0.9$ ($r_i$ and $r_o$ being the inner and outer radius), Prandtl numbers ${\rm Pr} \in[10^{-4},10^{-1}]$, and Taylor numbers ${\rm Ta}\in[10^{4},10^{12}]$. We are particularly interested in the form of the convective cell pattern that develops, and in its time scales, since this may have observational consequences. For a fixed ${\rm Ta}<10^{9}$ and by decreasing ${\rm Pr}$ from 0.1 to $10^{-4}$ a transition between spiralling columnar (SC) and equatorially-attached (EA) modes, and a transition between EA and equatorially antisymmetric or symmetric polar (AP/SP) weakly multicellular modes are found. The latter modes are preferred at very low ${\rm Pr}$. Surprisingly, for ${\rm Ta}>3\times 10^{9}$ the unicellular polar modes become also preferred at moderate ${\rm Pr}\sim10^{-2}$ because two new transition curves between EA and AP/SP and between AP/SP and SC modes are born at a triple-point bifurcation. The dependence on ${\rm Pr}$ and ${\rm Ta}$ of the transitions is studied to estimate the type of modes, and their critical parameters, preferred at different stellar regimes.

Accreting neutron stars (NS) can exhibit high frequency modulations in their lightcurves during thermonuclear X-ray bursts, known as burst oscillations. The frequencies can be offset from the spin frequency of the NS by several Hz, and can drift by 1-3 Hz. One possible explanation is a mode in the bursting ocean, the frequency of which would decrease (in the rotating frame) as the burst cools, hence explaining the drifts. Most burst oscillations have been observed during H/He triggered bursts, however there has been one observation of oscillations during a superburst; hours' long Type I X-ray bursts caused by unstable carbon burning deeper in the ocean. This paper calculates the frequency evolution of an oceanic r-mode during a superburst. The rotating frame frequency varies during the burst from 4-14 Hz, and is sensitive to the background parameters, in particular the temperature of the ocean and ignition depth. This calculation is compared to the superburst oscillations observed on 4U-1636-536. The predicted mode frequencies ($\sim$ 10 Hz) would require a spin frequency of $\sim$ 592 Hz to match observations; 6 Hz higher than the spin inferred from an oceanic r-mode model for the H/He triggered burst oscillations. This model also over-predicts the frequency drift during the superburst by 90 %.