Current state-of-the-art models of the solar convection zone consist of solutions to the Navier-Stokes equations in rotating, 3D spherical shells. Such models are highly sensitive to the choice of boundary conditions. Here, we present two suites of simulations differing only in their outer thermal boundary condition, which is either one of fixed-entropy or fixed-entropy-gradient. We find that the resulting differential rotation is markedly different between the two sets. The fixed-entropy-gradient simulations have strong differential rotation contrast and isocontours tilted along radial lines (in good agreement with the Sun's interior rotation revealed by helioseismology), whereas the fixed-entropy simulations have weaker contrast and contours tilted in the opposite sense. We examine in detail the force balances in our models and find that the poleward transport of heat by Busse columns drives a thermal wind responsible for the different rotation profiles. We conclude that the Sun's strong differential rotation along radial lines may result from the solar emissivity being invariant with latitude (which is similar to the fixed-entropy-gradient condition in our models) and the poleward transport of heat by Busse columns. In future work on convection in the solar context, we strongly advise modelers to use a fixed-gradient outer boundary condition.

The calling card of solar magnetism is the sunspot cycle, during which sunspots regularly reverse their polarity sense every 11 years. However, a number of more complicated time-dependent behaviors have also been identified. In particular, there are temporal modulations associated with active longitudes and hemispheric asymmetry, when sunspots appear at certain solar longitudes or else in one hemisphere preferentially. So far, a direct link between between this asymmetric temporal behavior and the underlying solar dynamo has remained elusive. In this work, we present results from global, 3D magnetohydrodynamic (MHD) simulations, which for the first time display both behavior reminiscent of the sunspot cycle (regular polarity reversals and equatorward migration of internal magnetic field) and asymmetric, irregular behavior that in the simulations we interpret as active longitudes and hemispheric asymmetry. The simulations are thus bistable, in that the turbulent convection can stably support two distinct flavors of magnetism at different times, in superposition, or with smooth transitions from one state to the other. We discuss this new family of dynamo models in the context of the extensive observations of the Sun's surface magnetic field with the Solar and Heliospheric Observatory (SOHO) and the Solar Dynamics Observatory (SDO), as well as earlier observations of sunspot number and synoptic maps. We suggest that the solar dynamo itself may be bistable in nature, exhibiting two types of temporal behavior in the magnetic field.

The advent of helioseismology has determined in detail the average rotation rate of the Sun as a function of radius and latitude. These data immediately reveal two striking boundary layers of shear in the solar convection zone (CZ): a tachocline at the base, where the differential rotation of the CZ transitions to solid-body rotation in the radiative zone, and a 35-Mm-thick near-surface shear layer (NSSL) at the top, where the rotation rate slows by about 5% with increasing radius. Though asteroseismology cannot probe the differential rotation of distant stars to the same level of detail that helioseismology can achieve for the Sun, it is possible that many cool stars with outer convective envelopes possess similar differential rotation characteristics, including both a tachocline and a NSSL. Here we present the results of 3D global hydrodynamic simulations of spherical-shell convection for a Sun-like star at different levels of density contrast across the shell. The simulations with high stratification possess characteristics of near-surface shear, especially at low latitudes. We discuss in detail the dynamical balance of torques giving rise to the NSSL in our models and interpret what these balances imply for the real Sun. We further discuss the dynamical causes that may serve to wipe out near-surface shear at high latitudes, and conclude by offering some theories as to how this problem might be tackled in future work.

The dynamical origins of the Sun's tachocline and near-surface shear layer (NSSL) are still not well understood. We have attempted to self-consistently reproduce a NSSL in numerical simulations of a solar-like convection zone by increasing the density contrast across rotating, 3D spherical shells. We explore the hypothesis that high density contrast leads to near-surface shear by creating a rotationally unconstrained layer of fast flows near the outer surface. Although our high-contrast models do have near-surface shear, it is confined primarily to low latitudes (between $\pm15^\circ$). Two distinct types of flow structures maintain the shear dynamically: rotationally $\textit{constrained}$ Busse columns aligned with the rotation axis and fast, rotationally $\textit{unconstrained}$ downflow plumes that deplete angular momentum from the outer fluid layers. The plumes form at all latitudes, and in fact are more efficient at transporting angular momentum inward at high latitudes. The presence of Busse columns at low latitudes thus appears essential to creating near-surface shear in our models. We conclude that a solar-like NSSL is unobtainable from a rotationally unconstrained outer fluid layer alone. In numerical models, the shear is eliminated through the advection of angular momentum by the meridional circulation. Therefore, a detailed understanding how the solar meridional circulation is dynamically achieved will be necessary to elucidate the origin of the Sun's NSSL.

We study the Rossby wave instability (RWI) in a thin accretion disc threaded by an initially toroidal magnetic field using the magnetohydrodynamics (MHD) code PLUTO. We find that for plasma $10^1 < \beta < 10^3$ , the growth rate and late-time density are suppressed, whereas for plasma $\beta > 10^3$, the magnetic field has negligible effect. The initially toroidal field is twisted inside the vortex and expelled at late times. This creates a radially directed field outside the vortex region, which can be understood via a simple kinematic model of a magnetic field in a rotating fluid and may be observable via polarised dust emission.