The new data delivered by NASA's Juno spacecraft significantly increase our understanding of Jupiter's internal dynamics. The gravity data constrain the depth of the zonal flows observed at cloud level and suggest that they slow down considerably at a depth of about $0.96\,r_J$, where $r_J$ is the mean radius at the one bar level. Juno's magnetometer reveals the planet's internal magnetic field. We combine the new zonal flow and magnetic field models with an updated electrical conductivity profile to assess the zonal wind induced dynamo action, concentrating on the outer part of Jupiter's molecular hydrogen region where the conductivity increases very rapidly with depth. Dynamo action remains quasi-stationary and can thus reasonably be estimated where the magnetic Reynolds number remains smaller than one, which is roughly the region above $0.96\,r_J$. We calculate that the locally induced radial magnetic field reaches rms values of about $10^{-6}\,$T in this region and may just be detectable by the Juno mission. Very localized dynamo action and a distinct pattern that reflects the zonal wind system increases the chance to disentangle this locally induced field from the background field. The estimates of the locally induced currents also allow calculating the zonal flow related Ohmic heating and associated entropy production. The respective quantities remain below new revised predictions for the total dissipative heating and total entropy production in Jupiter for any of the explored model combinations. Thus neither Ohmic heating nor entropy production offer additional constraints on the depth of the zonal winds.

The Juno mission is delivering spectacular data of Jupiter's magnetic field, while the gravity measurements finally allow constraining the depth of the winds observed at cloud level. However, to which degree the zonal winds contribute to the planet's dynamo action remains an open question. Here we explore numerical dynamo simulations that include an Jupiter-like electrical conductivity profile and successfully model the planet's large scale field. We concentrate on analyzing the dynamo action in the Steeply Decaying Conductivity Region (SDCR) where the high conductivity in the metallic Hydrogen region drops to the much lower values caused by ionization effects in the very outer envelope of the planet. Our simulations show that the dynamo action in the SDCR is strongly ruled by diffusive effects and therefore quasi stationary. The locally induced magnetic field is dominated by the horizontal toroidal field, while the locally induced currents flow mainly in the latitudinal direction. The simple dynamics can be exploited to yield estimates of surprisingly high quality for both the induced field and the electric currents in the SDCR. These could be potentially be exploited to predict the dynamo action of the zonal winds in Jupiter's SDCR but also in other planets.

The Juno mission will measure Jupiter's magnetic field with unprecedented precision and provide a wealth of additional data that will allow to constrain the planet's interior structure and dynamics. Here we analyse 66 numerical simulations in order to explore the sensitivity of the dynamo-generated magnetic field to the planets interior properties. The degree l=4 field model VIP4 and up-to-date interior models based on ab initio simulations serve as benchmarks. Our results suggest that VIP4-like magnetic fields can be found for a number of different models. We find that whether we assume an ideal gas or use the more realistic interior model based on ab initio simulations makes no difference. However, two other factors are important. Low Rayleigh number leads to strong axial dipole contribution while the axial dipole dominance is lost when the convective driving is too strong. The required intermediate range that yields Jupiter-like magnetic fields depends on the other system properties. The second factor is the convective magnetic Reynolds number profile Rmc(r), a product of the non-axisymmetric flow velocity and electrical conductivity. We find that the depth where Rmc exceeds about 50 is a good proxy for the top of the dynamo region. When the dynamo region sits too deep, the axial dipole is once more too dominant due to geometric reasons. We extrapolated our results to Jupiter and the result suggests that the Jovian dynamo extends to 95% of the planetary radius. The zonal flows in our simulations are dominated by an equatorial jet largely confined to the molecular layer. Where the jet reaches down to higher electrical conductivities, it is gives rise to a secondary alpha-Omega dynamo that modifies the dipole-dominated field produced deeper in the planet. This secondary dynamo can lead to strong magnetic field patches at low latitudes that seem compatible with the VIP4 model.

We present a non-supersymmetric scenario in which the R-parity symmetry $R_P = (-1)^{3(B-L)+2s}$ arises as a result of spontaneous gauge symmetry breaking, leading to a viable Dirac fermion WIMP dark matter candidate. Direct detection in nuclear recoil experiments probes dark matter masses around $2-5$ TeV for $M_{Z^{\prime}} \sim 3-4$ TeV consistent with searches at the LHC, while lepton flavor violation rates and flavor changing neutral currents in neutral meson systems lie within reach of upcoming experiments.

We discuss here a purely hydrodynamical mechanism to invert the sign of the kinetic helicity, which plays a key role in determining the direction of propagation of cyclical magnetism in most models of dynamo action by rotating convection. Such propagation provides a prominent, and puzzling constraint on dynamo models. In the Sun, active regions emerge first at mid-latitudes, then appear nearer the equator over the course of a cycle, but most previous global-scale dynamo simulations have exhibited poleward propagation (if they were cyclical at all). Here, we highlight some simulations in which the direction of propagation of dynamo waves is altered primarily by an inversion of the kinetic helicity throughout much of the interior, rather than by changes in the differential rotation. This tends to occur in cases with a low Prandtl number and internal heating, in regions where the local density gradient is relatively small. We analyse how this inversion arises, and contrast it to the case of convection that is either highly columnar (i.e., rapidly rotating) or locally very stratified; in both of those situations, the typical profile of kinetic helicity (negative throughout most of the northern hemisphere) instead prevails.

We study rotating thermal convection in spherical shells. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from $10^{-3}$ to $10^{-5}$. The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh-B\'enard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Furthermore, we explore the influence of the magnetic field on the {\em non-zonal} flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at $E=10^{-5}$ magnetic field significantly affects the convection and a maximum gain of about 30\% (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at $E=10^{-5}$ is probably in a `magnetostrophic' regime.

Spacecraft data reveal a very Earth-like Jovian magnetic field. This is surprising since numerical simulations have shown that the vastly different interiors of terrestrial and gas planets can strongly affect the internal dynamo process. Here we present the first numerical dynamo that manages to match the structure and strength of the observed magnetic field by embracing the newest models for Jupiter's interior. Simulated dynamo action primarily occurs in the deep high electrical conductivity region while zonal flows are dynamically constrained to a strong equatorial jet in the outer envelope of low conductivity. Our model reproduces the structure and strength of the observed global magnetic field and predicts that secondary dynamo action associated to the equatorial jet produces banded magnetic features likely observable by the Juno mission. Secular variation in our model scales to about 2000 nT per year and should also be observable during the one year nominal mission duration.

Dynamos operating in the interiors of rapidly rotating planets and low-mass stars might belong to a similar category where rotation plays a vital role. We quantify this similarity using scaling laws. We analyse direct numerical simulations of Boussinesq and anelastic spherical shell dynamos. These dynamos represent simplified models which span from Earth-like planets to rapidly rotating low-mass stars. We find that magnetic field and velocity in these dynamos are related to the available buoyancy power via a simple power law which holds over wide variety of control parameters.

Observations of active M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. We detail the analogy between some anelastic dynamo simulations and spectropolarimetric observations of 23 M stars. In numerical models, the relative contribution of inertia and Coriolis force in the global force balance -estimated by the so-called local Rossby number- is known to have a strong impact on the magnetic field geometry. We discuss the relevance of this parameter in setting the large-scale magnetic field of M dwarfs.

The observed surface dynamics of Jupiter and Saturn is dominated by a banded system of zonal winds. Their depth remains unclear but they are thought to be confined to the very outer envelopes where hydrogen remains molecular and the electrical conductivity is small. The dynamo maintaining the dipole-dominated magnetic fields of both gas giants likely operates in the deeper interior where hydrogen assumes a metallic state. Here, we present numerical simulations that attempt to model both the zonal winds and the interior dynamo action in an integrated approach. Using the anelastic version of the MHD code MagIC, we explore the effects of density stratification and radial electrical conductivity variation. The electrical conductivity is mostly assumed to remain constant in the thicker inner metallic region and it decays exponentially towards the outer boundary throughout the molecular envelope. Our results show that the combination of stronger density stratification and weaker conducting outer layer is essential for reconciling dipole dominated dynamo action and a fierce equatorial zonal jet. Previous simulations with homogeneous electrical conductivity show that both are merely exclusive, with solutions either having strong zonal winds and multipolar magnetic fields or weak zonal winds and dipole-dominated magnetic fields. All jets tend to be geostrophic and therefore reach right through the convective shell in our simulations. The particular setup explored here allows a strong equatorial jet to remain confined to the weaker conducting outer region where it does not interfere with the deeper seated dynamo action. The flanking mid to high latitude jets, on the other hand, have to remain faint to yield a strongly dipolar magnetic field. The fiercer jets on Jupiter and Saturn only seem compatible with the observed dipolar fields when they remain confined to a weaker conducting outer layer.

Numerical dynamo models always employ parameter values that differ by orders of magnitude from the values expected in natural objects. However, such models have been successful in qualitatively reproducing properties of planetary and stellar dynamos. This qualitative agreement fuels the idea that both numerical models and astrophysical objects may operate in the same asymptotic regime of dynamics. This can be tested by exploring the scaling behavior of the models. For convection-driven incompressible spherical shell dynamos with constant material properties, scaling laws had been established previously that relate flow velocity and magnetic field strength to the available power. Here we analyze 273 direct numerical simulations using the anelastic approximation, involving also cases with radius-dependent magnetic, thermal and viscous diffusivities. These better represent conditions in gas giant planets and low-mass stars compared to Boussinesq models. Our study provides strong support for the hypothesis that both mean velocity and mean magnetic field strength scale as a function of power generated by buoyancy forces in the same way for a wide range of conditions.

Context: observations of rapidly rotating M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. In dynamo models, the relative importance of inertia in the force balance -- quantified by the local Rossby number -- is known to have a strong impact on the magnetic field geometry. Aims: we aim to assess the relevance of the local Rossby number in controlling the large-scale magnetic field geometry of M dwarfs. Methods: we explore the similarities between anelastic dynamo models in spherical shells and observations of active M-dwarfs, focusing on field geometries derived from spectropolarimetric studies. To do so, we construct observation-based quantities aimed to reflect the diagnostic parameters employed in numerical models. Results: the transition between dipole-dominated and multipolar large-scale fields in early to mid M dwarfs is tentatively attributed to a Rossby number threshold. We interpret late M dwarfs magnetism to result from a dynamo bistability occurring at low Rossby number. By analogy with numerical models, we expect different amplitudes of differential rotation on the two dynamo branches.

Context: dynamo action in giant planets and rapidly rotating stars leads to a broad variety of magnetic field geometries including small scale multipolar and large scale dipole-dominated topologies. Previous dynamo models suggest that solutions become multipolar once inertia becomes influential. Being tailored for terrestrial planets, most of these models neglected the background density stratification. Aims: we investigate the influence of the density stratification on convection-driven dynamo models. Methods: three-dimensional nonlinear simulations of rapidly rotating spherical shells are employed using the anelastic approximation to incorporate density stratification. A systematic parametric study for various density stratifications and Rayleigh numbers allows to explore the dependence of the magnetic field topology on these parameters. Results: anelastic dynamo models tend to produce a broad range of magnetic field geometries that fall on two distinct branches with either strong dipole-dominated or weak multipolar fields. As long as inertia is weak, both branches can coexist but the dipolar branch vanishes once inertia becomes influential. The dipolar branch also vanishes for stronger density stratifications. The reason is the concentration of the convective columns in a narrow region close to the outer boundary equator, a configuration that favors non-axisymmetric solutions. In multipolar solutions, zonal flows can become significant and participate in the toroidal field generation. Parker dynamo waves may then play an important role close to onset of dynamo action leading to a cyclic magnetic field behavior. Conclusion: Our simulations also suggest that the fact that late M dwarfs have dipolar or multipolar magnetic fields can be explained in two ways. They may differ either by the relative influence of inertia or fall into the regime where both types of solutions coexist.