-
Supernova Pointing Capabilities of DUNE
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
M. Adamowski,
D. Adams,
M. Adinolfi,
C. Adriano,
A. Aduszkiewicz,
J. Aguilar,
B. Aimard,
F. Akbar,
K. Allison,
S. Alonso Monsalve,
M. Alrashed,
A. Alton,
R. Alvarez,
T. Alves,
H. Amar,
P. Amedo,
J. Anderson,
D. A. Andrade
, et al. (1340 additional authors not shown)
Abstract:
The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electr…
▽ More
The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on $^{40}$Ar and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called ``brems flipping'', as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE's burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage.
△ Less
Submitted 14 July, 2024;
originally announced July 2024.
-
Search for Dark Sector by Repurposing the UVX Brazilian Synchrotron
Authors:
L. Duarte,
L. Lin,
M. Lindner,
V. Kozhuharov,
S. V. Kuleshov,
A. S. de Jesus,
F. S. Queiroz,
Y. Villamizar,
H. Westfahl Jr
Abstract:
We propose the first Search for Dark Sector at the Brazilian Synchrotron Light Laboratory, site of Sirius, a fourth-generation storage ring. We show that UVX, Sirius predecessor, can be a promising dark sector detector, SeDS, with unprecedented sensitivity. The search is based on a 1-3 GeV positron beam impinging on a thick target leading the $e^+ e^- \rightarrow γA'$ reaction, followed by a missi…
▽ More
We propose the first Search for Dark Sector at the Brazilian Synchrotron Light Laboratory, site of Sirius, a fourth-generation storage ring. We show that UVX, Sirius predecessor, can be a promising dark sector detector, SeDS, with unprecedented sensitivity. The search is based on a 1-3 GeV positron beam impinging on a thick target leading the $e^+ e^- \rightarrow γA'$ reaction, followed by a missing mass spectrum event reconstruction. We show that SeDS has the potential to probe dark photons with masses up to 55 MeV and kinetic coupling down to $ε\sim 10^{-14}$ within months of data. Therefore, such experiment would constitute the best dark photon probe worldwide in the 10-55 MeV mass range, being able to probe an unexplored region of parameter space.
△ Less
Submitted 10 June, 2022;
originally announced June 2022.
-
Searching for solar KDAR with DUNE
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
D. Adams,
M. Adinolfi,
A. Aduszkiewicz,
J. Aguilar,
Z. Ahmad,
J. Ahmed,
B. Ali-Mohammadzadeh,
T. Alion,
K. Allison,
S. Alonso Monsalve,
M. Alrashed,
C. Alt,
A. Alton,
P. Amedo,
J. Anderson,
C. Andreopoulos,
M. Andreotti,
M. P. Andrews
, et al. (1157 additional authors not shown)
Abstract:
The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search.…
▽ More
The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search. In this work, we evaluate the proposed KDAR neutrino search strategies by realistically modeling both neutrino-nucleus interactions and the response of DUNE. We find that, although reconstruction of the neutrino energy and direction is difficult with current techniques in the relevant energy range, the superb energy resolution, angular resolution, and particle identification offered by DUNE can still permit great signal/background discrimination. Moreover, there are non-standard scenarios in which searches at DUNE for KDAR in the Sun can probe dark matter interactions.
△ Less
Submitted 26 October, 2021; v1 submitted 19 July, 2021;
originally announced July 2021.
-
Hunting down the cause of solar magnetism
Authors:
M. Viviani,
A. Prabhu,
J. Warnecke,
L. Duarte,
J. Pekkilä,
M. Rheinhardt,
M. J. Käpylä
Abstract:
To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions…
▽ More
To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions and state-of-the-art observations. We present here the first test-field measurements from our higher-resolution runs with improved heat conduction description. They indicate significant changes in the profiles of the most crucial inductive effect related to solar and stellar dynamo mechanisms. Higher resolution runs, currently undertaken, will bring us into an even more turbulent regime, in which we will be able to study, for the first time, the interaction of small- and large-scale dynamos in a quantitative way.
△ Less
Submitted 4 February, 2021;
originally announced February 2021.
-
Dynamo Action of Jupiter's Zonal Winds
Authors:
Johannes Wicht,
Thomas Gastine,
Lucia D. V. Duarte,
Wieland Dietrich
Abstract:
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. W…
▽ More
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.
△ Less
Submitted 21 June, 2019;
originally announced June 2019.
-
Dynamo Action in the Steeply Decaying Conductivity Region of Jupiter-like Dynamo Models
Authors:
Johannes Wicht,
Thomas Gastine,
Lucia D. V. Duarte
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 August, 2018;
originally announced August 2018.
-
R-parity as a residual gauge symmetry : probing a theory of cosmological dark matter
Authors:
Alexandre Alves,
Giorgio Arcadi,
P. V. Dong,
Laura Duarte,
Farinaldo S. Queiroz,
José W. F. Valle
Abstract:
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 v…
▽ More
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.
△ Less
Submitted 13 December, 2016;
originally announced December 2016.
-
Physical conditions for Jupiter-like dynamo models
Authors:
Lucia D. V. Duarte,
Johannes Wicht,
Thomas Gastine
Abstract:
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 in…
▽ More
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.
△ Less
Submitted 12 September, 2017; v1 submitted 8 December, 2016;
originally announced December 2016.
-
Helicity inversion in spherical convection as a means for equatorward dynamo wave propagation
Authors:
Lúcia D. V. Duarte,
Johannes Wicht,
Matthew K. Browning,
Thomas Gastine
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 December, 2015; v1 submitted 18 November, 2015;
originally announced November 2015.
-
Effect of shear and magnetic field on the heat-transfer efficiency of convection in rotating spherical shells
Authors:
Rakesh K. Yadav,
Thomas Gastine,
Ulrich R. Christensen,
Lucia Duarte,
Ansgar Reiners
Abstract:
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 an…
▽ More
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énard 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.
△ Less
Submitted 17 November, 2015; v1 submitted 13 July, 2015;
originally announced July 2015.
-
Explaining Jupiter's magnetic field and equatorial jet dynamics
Authors:
T. Gastine,
J. Wicht,
L. Duarte,
M. Heimpel,
A. Becker
Abstract:
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 Jupi…
▽ More
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.
△ Less
Submitted 25 July, 2014; v1 submitted 22 July, 2014;
originally announced July 2014.
-
Bridging planets and stars using scaling laws in anelastic spherical shell dynamos
Authors:
Rakesh K. Yadav,
Thomas Gastine,
Ulrich R. Christensen,
Lucia D. V. Duarte
Abstract:
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 s…
▽ More
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.
△ Less
Submitted 1 October, 2013;
originally announced October 2013.
-
What controls the large-scale magnetic fields of M dwarfs?
Authors:
T. Gastine,
J. Morin,
L. Duarte,
A. Reiners,
U. Christensen,
J. Wicht
Abstract:
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 Ross…
▽ More
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.
△ Less
Submitted 24 September, 2013;
originally announced September 2013.
-
Consistent scaling laws in anelastic spherical shell dynamos
Authors:
Rakesh K. Yadav,
Thomas Gastine,
Ulrich R. Christensen,
Lúcia D. V. Duarte
Abstract:
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.…
▽ More
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.
△ Less
Submitted 25 June, 2013; v1 submitted 23 April, 2013;
originally announced April 2013.
-
What controls the magnetic geometry of M dwarfs?
Authors:
T. Gastine,
J. Morin,
L. Duarte,
A. Reiners,
U. R. Christensen,
J. Wicht
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 December, 2012;
originally announced December 2012.
-
Anelastic dynamo models with variable electrical conductivity: an application to gas giants
Authors:
Lúcia D. V. Duarte,
Thomas Gastine,
Johannes Wicht
Abstract:
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…
▽ More
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.
△ Less
Submitted 26 June, 2013; v1 submitted 11 October, 2012;
originally announced October 2012.
-
Dipolar versus multipolar dynamos: the influence of the background density stratification
Authors:
T. Gastine,
L. Duarte,
J. Wicht
Abstract:
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 stratific…
▽ More
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.
△ Less
Submitted 30 August, 2012;
originally announced August 2012.
