Explore the vast range of titles available.

Semiconvection occurs in regions of stars and planets that are unstable to overturning convection according to the Schwarzschild criterion, but stable according to the Ledoux criterion. Previous simulations in Cartesian boxes have advanced our understanding of the semiconvective instability, layer formation, and transport properties. However, much less is known about semiconvection in spherical geometry and under the influence of rotation or magnetic fields. We present 3D simulations of semiconvection in the full sphere (including r = 0) and accounting for rotation. We find that the formation and evolution of semiconvective layers in nonrotating spheres occurs in a similar way to that in nonrotating Cartesian boxes, in the sense that the critical density ratio at which layers are expected to form is approximately the same in both geometries. Layers rapidly merge once they form, ultimately leading to a fully mixed convective sphere. The transport properties measured through the Nusselt numbers and the buoyancy flux ratio are also similar to results from previous studies in boxes. When rotation is added to the system, layer formation and evolution proceed in a similar fashion to the nonrotating runs. However, rotation hampers the radial transport of heat and composition, and, as a result, the time that it takes for the sphere to become fully mixed gets longer as the flow becomes more rotationally constrained. We also find that semiconvective layers exhibit spherical mixing in nonrotating cases, whereas in rotating cases the mixing becomes more cylindrical. We discuss what is needed for future work to build more realistic models.

We investigate crystallization-driven convection in carbon–oxygen white dwarfs. We present a version of the mixing length theory that self-consistently includes the effects of thermal diffusion and composition gradients, and provides solutions for the convective parameters based on the local heat and composition fluxes. Our formulation smoothly transitions between the regimes of fast adiabatic convection at large Peclet number and slow thermohaline convection at low Peclet number. It also allows for both thermally driven and compositionally driven convection, including correctly accounting for the direction of heat transport for compositionally driven convection in a thermally stable background. We use the MESA stellar evolution code to calculate the composition and heat fluxes during crystallization in different models of cooling white dwarfs, and determine the regime of convection and the convective velocity. We find that convection occurs in the regime of slow thermohaline convection during most of the cooling history of the star. However, at the onset of crystallization, the composition flux is large enough to drive fast overturning convection for a short time (∼10 Myr). We estimate the convective velocities in both of these phases and discuss the implications for explaining observed white dwarf magnetic fields with crystallization-driven dynamos.

The origin of large magnetic fields (≳106 G) in isolated white dwarfs is not clear. One possible explanation is that crystallization of the star’s core drives compositional convection, which when combined with the star’s rotation, can drive a dynamo. However, whether convection is efficient enough to explain the large intensity of the observed magnetic fields is still under debate. Recent work has shown that convection in cooling white dwarfs spans two regimes: efficient convection at the onset of crystallization, and thermohaline convection during most of the star’s cooling history. Here, we calculate the properties of crystallization-driven convection for cooling models of several white dwarfs of different masses. We combine mixing-length theory with scalings from magnetorotational convection to estimate the typical magnitude of the convective velocity and induced magnetic field for both scenarios. In the thermohaline regime, we find velocities ∼10−6–10−5 cm s−1, with fields restricted to ≲ 100 G. However, when convection is efficient, the flow velocity can reach magnitudes of ∼102–103 cm s−1, with fields of ∼106–108 G, independent of the star’s rotation rate. Thus, dynamos driven at the onset of crystallization could explain the large intensity magnetic fields measured for single white dwarfs.

We investigate heat transport associated with compositionally driven convection driven by crystallization at the ocean–crust interface in accreting neutron stars, or growth of the solid core in cooling white dwarfs. We study the effect of thermal diffusion and rapid rotation on the convective heat transport, using both mixing length theory and numerical simulations of Boussinesq convection. We determine the heat flux, composition gradient, and Péclet number, Pe (the ratio of thermal diffusion time to convective turnover time) as a function of the composition flux. We find two regimes of convection with a rapid transition between them as the composition flux increases. At small Pe, the ratio between the heat flux and composition flux is independent of Pe, because the loss of heat from convecting fluid elements due to thermal diffusion is offset by the smaller composition gradient needed to overcome the reduced thermal buoyancy. At large Pe, the temperature gradient approaches the adiabatic gradient, saturating the heat flux. We discuss the implications for cooling of neutron stars and white dwarfs. Convection in neutron stars spans both regimes. We find rapid mixing of neutron star oceans, with a convective turnover time of the order of weeks to minutes depending on rotation. Except during the early stages of core crystallization, white dwarf convection is in the thermal-diffusion-dominated fingering regime. We find convective velocities much smaller than recent estimates for crystallization-driven dynamos. The small fraction of energy carried as kinetic energy calls into question the effectiveness of crystallization-driven dynamos as an explanation for observed magnetic fields in white dwarfs.

Neutron stars show a steady decrease in their rotational frequency, occasionally interrupted by sudden spin-up events called glitches. The dynamics of a neutron star after a glitch involve the transfer of angular momentum from the crust (where the glitch is presumed to originate) to the liquid core, causing the core to spin up. The crust–core coupling, which determines how quickly this spin-up proceeds, can be achieved through various physical processes, including Ekman pumping, superfluid vortex-mediated mutual friction, and magnetic fields. Although the complex nature of these mechanisms has made it difficult to study their combined effects, analytical estimations for individual processes reveal that spin-up timescales vary according to the relative strength of Coriolis, viscous, and mutual friction forces, as well as the magnetic field. However, experimental and numerical validations of those analytical predictions are limited. In this paper, we focus on viscous effects and mutual friction. We conduct nonlinear hydrodynamical simulations of the spin-up problem in a two-component fluid by solving the incompressible Hall–Vinen–Bekarevich–Khalatnikov equations in the full sphere (i.e., including r = 0) for the first time. We find that the viscous (normal) component accelerates due to Ekman pumping, although the mutual friction coupling to the superfluid component alters the spin-up dynamics compared to the single-fluid scenario. Close to the sphere’s surface, the response of the superfluid is accurately described by the mutual friction timescale irrespective of its coupling strength with the normal component. However, as we move deeper into the sphere, the superfluid accelerates on different timescales due to the slow viscous spin-up of the internal normal fluid layers. We discuss potential implications for neutron stars, and requirements for future work to build more realistic models.

Thermal convection in rotating stars and planets drives anisotropic turbulence and differential rotation, both capable of feeding energy into global oscillations. Using 3D simulations of rotating convection in spherical shells, we show that inertial modes—oscillations restored by the Coriolis force—emerge naturally in rotationally constrained turbulence, without imposing any external forcing other than thermal/buoyancy driving. By varying the rotation rate at fixed Rayleigh number, we find that coherent modes appear only when the convective Rossby number, the ratio of the rotation period to the convective turnover time, falls below about one-half, where rotation dominates the dynamics. These modes are mostly retrograde in the rotating frame, equatorially symmetric, and confined to mid and high latitudes, with discrete frequencies well below twice the background rotation rate. At lower viscosities, or a smaller Prandtl number, mode excitation becomes more efficient and a broader spectrum of inertial modes emerges. While the precise excitation mechanism remains uncertain, our results suggest that the modes are driven by instabilities due to differential rotation rather than stochastic forcing by convection. We conclude that similar inertial modes are likely to exist in the interiors of giant planets and stars, though their low frequencies will make them difficult to detect.

Radio images of protoplanetary disks demonstrate that dust grains tend to organize themselves into rings. These rings may be a consequence of dust trapping within gas pressure maxima, wherein the local high dust-to-gas ratio is expected to trigger the formation of planetesimals and eventually planets. We revisit the behavior of dust near gas pressure perturbations enforced by a planet in two-dimensional, shearing-box simulations. While dust grains collect into generally long-lived rings, particles with a small Stokes parameter τ s < 0.1 tend to advect out of the ring within a few drift timescales. Scaled to the properties of ALMA disks, we find that rings composed of larger particles (τ s ≥ 0.1) can nucleate a dust clump massive enough to trigger pebble accretion, which proceeds to ingest the entire dust ring well within ∼1 Myr. To ensure the survival of the dust rings, we favor a nonplanetary origin and typical grain size τ s ≲ 0.05–0.1. Planet-driven rings may still be possible but if so we would expect the orbital distance of the dust rings to be larger for older systems.

Recent measurements of Jupiter's gravitational field (by Juno) and seismology of Saturn's rings (by Cassini) strongly suggest that both planets have a stably stratified core that still possesses a primordial gradient in the concentration of heavy elements. The existence of such a “diffusely” stratified core has been a surprise as it was long expected that the Jovian planets should be fully convective and hence fully mixed. A vigorous zone of convection, driven by surface cooling, forms at the surface and deepens through entrainment of fluid from underneath. In fact, it was believed that this convection zone should grow so rapidly that the entire planet would be consumed in less than a million years. Here we suggest that two processes, acting in concert, present a solution to this puzzle. All of the giant planets are rapidly rotating and have a cooling rate that declines with time. Both of these effects reduce the rate of fluid entrainment into the convection zone. Through the use of an analytic prescription of entrainment in giant planets, we demonstrate that these two effects, rotation and dwindling surface cooling, result in a convection zone that initially grows but eventually stalls. The depth to which the convective interface asymptotes depends on the rotation rate and on the stratification of the stable interior. Conversely, in a nonrotating planet, or in a planet that maintains a higher level of cooling than current models suggest, the convection zone deepens forever, eventually spanning the entire planet.

Recent observational constraints on the internal structure of Jupiter and Saturn suggest that these planets have “fuzzy” cores, i.e., gradients of the concentration of heavy elements that might span a large fraction of the planet’s radius. These cores could be composed of a semiconvective staircase, i.e., multiple convective layers separated by diffusive interfaces arising from double-diffusive instabilities. However, to date, no study has demonstrated how such staircases can avoid layer mergers and persist over evolutionary timescales. In fact, previous work has found that these mergers occur rapidly, leading to only a single convective layer. Using 3D simulations, we demonstrate that rotation prolongs the lifetime of a convective staircase by increasing the timescale for both layer merger and erosion of the interface between the final two layers. We present an analytic model for the erosion phase, predicting that rotation increases the erosion time by a factor of approximately Ro−1/2, where Ro is the Rossby number of the convective flows (the ratio of the rotation period to the convective turnover time). For Jovian conditions at early times after formation (when convection is vigorous enough to mix a large fraction of the planet), we find the erosion time to be roughly 109 yr in the nonrotating case and 1011 yr in the rotating case. If these timescales are confirmed with a larger suite of numerical simulations, the existence of convective staircases within the deep interiors of giant planets is a strong possibility, and rotation could be an important factor in the preservation of their fuzzy cores.

Recent measurements of Jupiter’s gravitational moments by the Juno spacecraft and seismology of Saturn’s rings suggest that the primordial composition gradients in the deep interior of these planets have persisted since their formation. One possible explanation is the presence of a double-diffusive staircase below the planet’s outer convection zone, which inhibits mixing across the deeper layers. However, hydrodynamic simulations have shown that these staircases are not long-lasting and can be disrupted by overshooting convection. In this Letter, we suggests that planetary rotation could be another factor for the longevity of primordial composition gradients. Using rotational mixing-length theory and 3D hydrodynamic simulations, we demonstrate that rotation significantly reduces both the convective velocity and the mixing of primordial composition gradients. In particular, for Jovian conditions at t ∼ 108 yr after formation, rotation reduces the convective velocity by a factor of 6, and in turn, the kinetic energy flux available for mixing gets reduced by a factor of 63 ∼ 200. This leads to an entrainment timescale that is more than 2 orders of magnitude longer than without rotation. We encourage future hydrodynamic models of Jupiter and other gas giants to include rapid rotation because the decrease in the mixing efficiency could explain why Jupiter and Saturn are not fully mixed.