-
The solar dynamo begins near the surface
Authors:
Geoffrey M Vasil,
Daniel Lecoanet,
Kyle Augustson,
Keaton J Burns,
Jeffrey S Oishi,
Benjamin P Brown,
Nicholas Brummell,
Keith Julien
Abstract:
The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called "torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common cause. Contrary to theories suggesting deep origins for these phenomena, helioseismology pin…
▽ More
The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called "torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common cause. Contrary to theories suggesting deep origins for these phenomena, helioseismology pinpoints low-latitude torsional oscillations to the Sun's outer 5-10%, the "Near-Surface Shear Layer". Within this zone, inwardly increasing differential rotation coupled with a poloidal magnetic field strongly implicates the Magneto-Rotational Instability prominent in accretion-disk theory and observed in laboratory experiments. Together, these two facts prompt the general question: Is it possible that the solar dynamo is a near-surface instability? Here, we report strong affirmative evidence in stark contrast to traditional paradigms focusing on the deeper tachocline. Simple analytic estimates show that the near-surface magneto-rotational instability better explains the spatiotemporal scales of the torsional oscillations and inferred subsurface magnetic field amplitudes. State-of-the-art numerical simulations corroborate these estimates and, strikingly, reproduce hemispherical magnetic current helicity laws. The dynamo resulting from a well-understood near-surface phenomenon improves prospects for accurate predictions of full magnetic cycles and space weather, impacting Earth's electromagnetic infrastructure.
△ Less
Submitted 13 April, 2024; v1 submitted 11 April, 2024;
originally announced April 2024.
-
The Puzzling Structure of Solar Convection: Window into the Dynamo
Authors:
Nicholas A. Featherstone,
Kyle C. Augustson,
Jonathan M. Aurnou,
Catherine Blume,
Benjamin P. Brown,
Nicholas Brummell,
Keaton J. Burns,
Michael A. Calkins,
Maria Camisassa,
Mausumi Dikpati,
Yuhong Fan,
J. R. Fuentes,
Gustavo Guerrero,
Bradley W. Hindman,
Keith Julien,
Irina N. Kitiashvili,
Lydia Korre,
Daniel Lecoanet,
Bhishek Manek,
Loren Matilsky,
Mark Miesch,
Nicholas J. Nelson,
Jeffrey S. Oishi,
Whitney T. Powers,
Matthias Rempel
, et al. (3 additional authors not shown)
Abstract:
The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective…
▽ More
The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements.
△ Less
Submitted 15 May, 2023;
originally announced May 2023.
-
How do tidal waves interact with convective vortices in rapidly-rotating planets and stars?
Authors:
Virgile Dandoy,
Junho Park,
Kyle Augustson,
Aurélie Astoul,
Stéphane Mathis
Abstract:
The dissipation of tidal inertial waves in planetary and stellar convective regions is one of the key mechanisms that drive the evolution of star-planet/planet-moon systems. In this context, the interaction between tidal inertial waves and turbulent convective flows must be modelled in a realistic and robust way. In the state-of-the-art simulations, the friction applied by convection on tidal wave…
▽ More
The dissipation of tidal inertial waves in planetary and stellar convective regions is one of the key mechanisms that drive the evolution of star-planet/planet-moon systems. In this context, the interaction between tidal inertial waves and turbulent convective flows must be modelled in a realistic and robust way. In the state-of-the-art simulations, the friction applied by convection on tidal waves is modelled most of the time by an effective eddy-viscosity. This approach may be valid when the characteristic length scales of convective eddies are smaller than those of tidal waves. However, it becomes highly questionable in the case where tidal waves interact with potentially stable large-scale vortices, as those observed at the pole of Jupiter and Saturn. They are potentially triggered by convection in rapidly-rotating bodies in which the Coriolis acceleration forms the flow in columnar vortical structures along the direction of the rotation axis. In this paper, we investigate the complex interactions between a tidal inertial wave and a columnar convective vortex. We use a quasi-geostrophic semi-analytical model of a convective columnar vortex. We perform linear stability analysis to identify the unstable regime and conduct linear numerical simulations for the interactions between the convective vortex and an incoming tidal inertial wave. We verify that in the unstable regime, an incoming tidal inertial wave triggers the most unstable mode of the vortex leading to turbulent dissipation. For stable vortices, the wave-vortex interaction leads to the momentum mixing while it creates a low-velocity region around the vortex core and a new wave-like perturbation in the form of a progressive wave radiating in the far field. The emission of this secondary wave is the strongest when the wavelength of the incoming wave is close to the characteristic size of the vortex.
△ Less
Submitted 8 March, 2023; v1 submitted 10 November, 2022;
originally announced November 2022.
-
Rossby numbers and stiffness values inferred from gravity-mode asteroseismology of rotating F- and B-type dwarfs: consequences for mixing, transport, magnetism, and convective penetration
Authors:
C. Aerts,
K. Augustson,
S. Mathis,
M. G. Pedersen,
J. S. G. Mombarg,
V. Vanlaer,
J. Van Beeck,
T. Van Reeth
Abstract:
Multi-dimensional (magneto-)hydrodynamical simulations of physical processes in stellar interiors depend on a multitude of uncalibrated free parameters, which set the spatial and time scales of their computations. We aim to provide an asteroseismic calibration of the wave and convective Rossby numbers, and of the stiffness at the interface between the convective core and radiative envelope of inte…
▽ More
Multi-dimensional (magneto-)hydrodynamical simulations of physical processes in stellar interiors depend on a multitude of uncalibrated free parameters, which set the spatial and time scales of their computations. We aim to provide an asteroseismic calibration of the wave and convective Rossby numbers, and of the stiffness at the interface between the convective core and radiative envelope of intermediate-mass stars. We deduce these quantities for rotating dwarfs from the observed properties of their identified gravity and gravito-inertial modes. We rely on near-core rotation rates and asteroseismic models of 26 B- and 37 F-type dwarf pulsators derived from 4-year Kepler space photometry, high-resolution spectroscopy and Gaia astrometry in the literature to deduce their convective and wave Rossby numbers. We compute the stiffness at the convection/radiation interface from the inferred maximum buoyancy frequency at the interface and the convective turnover frequency in the core. We use those asteroseismically inferred quantities to make predictions of convective penetration levels, local flux levels of gravito-inertial waves triggered by the convective core, and of the cores' potential rotational and magnetic states. Our sample of 63 gravito-inertial mode pulsators covers near-core rotation rates from almost zero up to the critical rate. The frequencies of their identified modes lead to models with stiffness values between $10^{2.69}$ and $10^{3.60}$ for the B-type pulsators, while those of F-type stars cover the range from $10^{3.47}$ to $10^{4.52}$. The convective Rossby numbers derived from the maximum convective diffusion coefficient in the convective core, based on mixing length theory and a value of the mixing length coefficient relevant for these pulsators, vary between $10^{-2.3}$ and $10^{-0.8}$ for B-type stars and $10^{-3}$ and $10^{-1.5}$ for F-type stars. (abridged)
△ Less
Submitted 12 December, 2021; v1 submitted 12 October, 2021;
originally announced October 2021.
-
Magnetic signatures on mixed-mode frequencies. I. An axisymmetric fossil field inside the core of red giants
Authors:
L. Bugnet,
V. Prat,
S. Mathis,
A. Astoul,
K. Augustson,
R. A. García,
S. Mathur,
L. Amard,
C. Neiner
Abstract:
The discovery of the moderate differential rotation between the core and the envelope of evolved solar-like stars could be the signature of a strong magnetic field trapped inside the radiative interior. The population of intermediate-mass red giants presenting a surprisingly low-amplitude of their mixed modes could also arise from the effect of an internal magnetic field. Indeed, stars more massiv…
▽ More
The discovery of the moderate differential rotation between the core and the envelope of evolved solar-like stars could be the signature of a strong magnetic field trapped inside the radiative interior. The population of intermediate-mass red giants presenting a surprisingly low-amplitude of their mixed modes could also arise from the effect of an internal magnetic field. Indeed, stars more massive than about 1.1Ms are known to develop a convective core during their main sequence, which could relax into a strong fossil magnetic field trapped inside the core of the star for the rest of its evolution. The observations of mixed modes can constitute an excellent probe of the deepest layers of evolved solar-like stars. The magnetic perturbation on mixed modes may thus be visible in asteroseismic data. To unravel which constraints can be obtained from observations, we theoretically investigate the effects of a plausible mixed axisymmetric magnetic field with various amplitudes on the mixed-mode frequencies of evolved solar-like stars. The first-order frequency perturbations are computed for dipolar and quadrupolar mixed modes. These computations are carried out for a range of stellar ages, masses, and metallicities. We show that typical fossil-field strengths of 0.1-1 MG, consistent with the presence of a dynamo in the convective core during the main sequence, provoke significant asymmetries on mixed-mode frequency multiplets during the red-giant branch. We show that these signatures may be detectable in asteroseismic data for field amplitudes small enough for the amplitude of the modes not to be affected by the conversion of gravity into Alfven waves inside the magnetised interior. Finally, we infer an upper limit for the strength of the field, and the associated lower limit for the timescale of its action, to redistribute angular momentum in stellar interiors.
△ Less
Submitted 1 February, 2021;
originally announced February 2021.
-
Probing the internal magnetism of stars using asymptotic magneto-asteroseismology
Authors:
Stéphane Mathis,
Lisa Bugnet,
Vincent Prat,
Kyle Augustson,
Savita Mathur,
Rafael A. Garcia
Abstract:
Our knowledge of the dynamics of stars has undergone a revolution thanks to the simultaneous large amount of high-quality photometric observations collected by space-based asteroseismology and ground-based high-precision spectropolarimetry. They allowed us to probe the internal rotation of stars and their surface magnetism in the whole Hertzsprung-Russell diagram. However, new methods should still…
▽ More
Our knowledge of the dynamics of stars has undergone a revolution thanks to the simultaneous large amount of high-quality photometric observations collected by space-based asteroseismology and ground-based high-precision spectropolarimetry. They allowed us to probe the internal rotation of stars and their surface magnetism in the whole Hertzsprung-Russell diagram. However, new methods should still be developed to probe the deep magnetic fields in those stars. Our goal is to provide seismic diagnoses that allow us to sound the internal magnetism of stars. Here, we focus on asymptotic low-frequency gravity modes and high-frequency acoustic modes. Using a first-order perturbative theory, we derive magnetic splittings of their frequencies as explicit functions of stellar parameters. As in the case of rotation, we show how asymptotic gravity and acoustic modes can allow us to probe the different components of the magnetic field in the cavities where they propagate. This demonstrates again the high potential of using mixed-modes when this is possible.
△ Less
Submitted 20 December, 2020;
originally announced December 2020.
-
The impact of a fossil magnetic field on dipolar mixed-mode frequencies in sub- and red-giant stars
Authors:
L. Bugnet,
V. Prat,
S. Mathis,
R. A. García,
S. Mathur,
K. Augustson,
C. Neiner,
M. J. Thompson
Abstract:
Stars more massive than $\sim 1.3$ M$_\odot$ are known to develop a convective core during the main-sequence: the dynamo process triggered by this convection could be the origin of a strong magnetic field inside the core of the star, trapped when it becomes stably stratified and for the rest of its evolution. The presence of highly magnetized white dwarfs strengthens the hypothesis of buried fossi…
▽ More
Stars more massive than $\sim 1.3$ M$_\odot$ are known to develop a convective core during the main-sequence: the dynamo process triggered by this convection could be the origin of a strong magnetic field inside the core of the star, trapped when it becomes stably stratified and for the rest of its evolution. The presence of highly magnetized white dwarfs strengthens the hypothesis of buried fossil magnetic fields inside the core of evolved low-mass stars. If such a fossil field exists, it should affect the mixed modes of red giants as they are sensitive to processes affecting the deepest layers of these stars. The impact of a magnetic field on dipolar oscillations modes was one of Pr. Michael J. Thompson's research topics during the 90s when preparing the helioseismic SoHO space mission. As the detection of gravity modes in the Sun is still controversial, the investigation of the solar oscillation modes did not provide any hint of the existence of a magnetic field in the solar radiative core. Today we have access to the core of evolved stars thanks to the asteroseismic observation of mixed modes from CoRoT, Kepler, K2 and TESS missions. The idea of applying and generalizing the work done for the Sun came from discussions with Pr. Michael Thompson in early 2018 before we loss him. Following the path we drew together, we theoretically investigate the effect of a stable axisymmetric mixed poloidal and toroidal magnetic field, aligned with the rotation axis of the star, on the mixed modes frequencies of a typical evolved low-mass star. This enables us to estimate the magnetic perturbations to the eigenfrequencies of mixed dipolar modes, depending on the magnetic field strength and the evolutionary state of the star. We conclude that strong magnetic fields of $\sim$ 1MG should perturbe the mixed-mode frequency pattern enough for its effects to be detectable inside current asteroseismic data.
△ Less
Submitted 15 December, 2020;
originally announced December 2020.
-
Lithium depletion and angular momentum transport in solar-type stars
Authors:
T. Dumont,
A. Palacios,
C. Charbonnel,
O. Richard,
L. Amard,
K. Augustson,
S. Mathis
Abstract:
Transport processes occurring in the radiative interior of solar-type stars are evidenced by the surface variation of light elements, in particular Li, and the evolution of their rotation rates. For the Sun, inversions of helioseismic data indicate that the radial profile of angular velocity in its radiative zone is nearly uniform, which implies the existence of angular momentum transport mechanis…
▽ More
Transport processes occurring in the radiative interior of solar-type stars are evidenced by the surface variation of light elements, in particular Li, and the evolution of their rotation rates. For the Sun, inversions of helioseismic data indicate that the radial profile of angular velocity in its radiative zone is nearly uniform, which implies the existence of angular momentum transport mechanisms. While there are many independent transport models for angular momentum and chemical species, there is a lack of self-consistent theories that permit stellar evolution models to simultaneously match the present-day observations of solar lithium abundances and radial rotation profiles. We explore how additional transport processes can improve the agreement between evolutionary models of rotating stars and observations. We constrain the resulting models by simultaneously using the evolution of the surface rotation rate and Li abundance in the solar-type stars of open clusters, and the solar surface and internal rotation profile as inverted from helioseismology. We show the relevance of penetrative convection for the depletion of Li. The rotational dependence of the depth of penetrative convection yields an anti-correlation between the initial rotation rate and Li depletion in our models of solar-type stars that is in agreement with the observed trend. Simultaneously, the addition of an ad hoc vertical viscosity leads to efficient transport of angular momentum between the core and the envelope. We also self-consistently compute for the first time the thickness of the tachocline and find that it is compatible with helioseismic estimations. However, the main sequence depletion of Li in solar-type stars is only reproduced when adding a parametric turbulent mixing below the convective envelope. The need for additional transport processes in stellar evolution models is confirmed.
△ Less
Submitted 7 December, 2020;
originally announced December 2020.
-
A model of rotating convection in stellar and planetary interiors: II -- gravito-inertial wave generation
Authors:
K. C. Augustson,
S. Mathis,
A. Astoul
Abstract:
Gravito-inertial waves are excited at the interface of convective and radiative regions and by the Reynolds stresses in the bulk of the convection zones of rotating stars and planets. Such waves have notable asteroseismic signatures in the frequency spectra of rotating stars, particularly among rapidly rotating early-type stars, which provides a means of probing their internal structure and dynami…
▽ More
Gravito-inertial waves are excited at the interface of convective and radiative regions and by the Reynolds stresses in the bulk of the convection zones of rotating stars and planets. Such waves have notable asteroseismic signatures in the frequency spectra of rotating stars, particularly among rapidly rotating early-type stars, which provides a means of probing their internal structure and dynamics. They can also transport angular momentum, chemical species, and energy from the excitation region to where they dissipate in radiative regions. To estimate the excitation and convective parameter dependence of the amplitude of those waves, a monomodal model for stellar and planetary convection as described in Paper I is employed, which provides the magnitude of the rms convective velocity as a function of rotation rate. With this convection model, two channels for wave driving are considered: excitation at a boundary between convectively stable and unstable regions and excitation due to Reynolds-stresses. Parameter regimes are found where the sub-inertial waves may carry a significant energy flux, depending upon the convective Rossby number, the interface stiffness, and the wave frequency. The super-inertial waves can also be enhanced, but only for convective Rossby numbers near unity. Interfacially excited waves have a peak energy flux near the lower cutoff frequency when the convective Rossby number of the flows that excite them are below a critical Rossby number that depends upon the stiffness of the interface, whereas that flux decreases when the convective Rossby number is larger than this critical Rossby number.
△ Less
Submitted 22 September, 2020;
originally announced September 2020.
-
Transport of angular momentum by stochastically excited waves as an explanation for the outburst of the rapidly rotating Be star HD49330
Authors:
C. Neiner,
U. Lee,
S. Mathis,
H. Saio,
C. C. Lovekin,
K. C. Augustson
Abstract:
HD49330 is a Be star that underwent an outburst during its five-month observation with the CoRoT satellite. An analysis of its light curve revealed several independent p and g pulsation modes, in addition to showing that the amplitude of the modes is directly correlated with the outburst. We modelled the results obtained with CoRoT. We modelled the flattening of the structure of the star due to ra…
▽ More
HD49330 is a Be star that underwent an outburst during its five-month observation with the CoRoT satellite. An analysis of its light curve revealed several independent p and g pulsation modes, in addition to showing that the amplitude of the modes is directly correlated with the outburst. We modelled the results obtained with CoRoT. We modelled the flattening of the structure of the star due to rapid rotation in two ways: Chandrasekhar-Milne's expansion and 2D structure computed with ROTORC. We then modelled kappa-driven pulsations. We also adapted the formalism of the excitation and amplitude of stochastically excited gravito-inertial modes to rapidly rotating stars, and we modelled those pulsations as well. We find that while pulsation p modes are excited by the kappa mechanism, the observed g modes are a result of stochastic excitation. In contrast, g and r waves are stochastically excited in the convective core and transport angular momentum to the surface, increasing its rotation rate. This destabilises the external layers of the star, which then emits transient stochastically excited g waves. These transient waves produce most of the low-frequency signal detected in the CoRoT data and ignite the outburst. During this unstable phase, p modes disappear at the surface because their cavity is broken. Following the outburst and ejection of the surface layer, relaxation occurs, making the transient g waves disappear and p modes reappear. This work includes the first coherent model of stochastically excited gravito-inertial pulsation modes in a rapidly rotating Be star. It provides an explanation for the correlation between the variation in the amplitude of frequencies detected in the CoRoT data and the occurrence of an outburst. This scenario could apply to other pulsating Be stars, providing an explanation to the long-standing questions surrounding Be outbursts and disks.
△ Less
Submitted 17 July, 2020;
originally announced July 2020.
-
Magnetism in Massive Stars
Authors:
Kyle C. Augustson
Abstract:
Massive stars are the drivers of star formation and galactic dynamics due to their relatively short lives and explosive demises, thus impacting all of astrophysics. Since they are so impactful on their environments, through their winds on the main sequence and their ultimate supernovae, it is crucial to understand how they evolve. Recent photometric observations with space-based platforms such as…
▽ More
Massive stars are the drivers of star formation and galactic dynamics due to their relatively short lives and explosive demises, thus impacting all of astrophysics. Since they are so impactful on their environments, through their winds on the main sequence and their ultimate supernovae, it is crucial to understand how they evolve. Recent photometric observations with space-based platforms such as CoRoT, K2, and now TESS have permitted access to their interior dynamics through asteroseismology, while ground-based spectropolarimetric measurements such as those of ESPaDOnS have given us a glimpse at their surface magnetic fields. The dynamics of massive stars involve a vast range of scales. Extant methods can either capture the long-term structural evolution or the short-term dynamics such as convection, magnetic dynamos, and waves due to computational costs. Thus, many mysteries remain regarding the impact of such dynamics on stellar evolution, but they can have strong implications both for how they evolve and what they leave behind when they die. Some of these dynamics including rotation, tides, and magnetic fields have been addressed in recent work, which is reviewed in this paper.
△ Less
Submitted 13 December, 2019;
originally announced December 2019.
-
Does magnetic field impact tidal dynamics inside the convective zone of low-mass stars along their evolution?
Authors:
A. Astoul,
S. Mathis,
C. Baruteau,
F. Gallet,
A. Strugarek,
K. C. Augustson,
A. S. Brun,
E. Bolmont
Abstract:
The dissipation of the kinetic energy of wave-like tidal flows within the convective envelope of low-mass stars is one of the key physical mechanisms that shapes the orbital and rotational dynamics of short-period exoplanetary systems. Although low-mass stars are magnetically active objects, the question of how the star's magnetic field impacts large-scale tidal flows and the excitation, propagati…
▽ More
The dissipation of the kinetic energy of wave-like tidal flows within the convective envelope of low-mass stars is one of the key physical mechanisms that shapes the orbital and rotational dynamics of short-period exoplanetary systems. Although low-mass stars are magnetically active objects, the question of how the star's magnetic field impacts large-scale tidal flows and the excitation, propagation and dissipation of tidal waves still remains open. Our goal is to investigate the impact of stellar magnetism on the forcing of tidal waves, and their propagation and dissipation in the convective envelope of low-mass stars as they evolve. We have estimated the amplitude of the magnetic contribution to the forcing and dissipation of tidally induced magneto-inertial waves throughout the structural and rotational evolution of low-mass stars (from M to F-type). For this purpose, we have used detailed grids of rotating stellar models computed with the stellar evolution code STAREVOL. The amplitude of dynamo-generated magnetic fields is estimated via physical scaling laws at the base and the top of the convective envelope. We find that the large-scale magnetic field of the star has little influence on the excitation of tidal waves in the case of nearly-circular orbits and coplanar hot-Jupiter planetary systems, but that it has a major impact on the way waves are dissipated. Our results therefore indicate that a full magneto-hydrodynamical treatment of the propagation and dissipation of tidal waves is needed to properly assess the impact of star-planet tidal interactions throughout the evolutionary history of low-mass stars hosting short-period massive planets.
△ Less
Submitted 23 September, 2019;
originally announced September 2019.
-
Probing the shape of the mixing profile and of the thermal structure at the convective core boundary through asteroseismology
Authors:
Mathias Michielsen,
May G. Pedersen,
Kyle C. Augustson,
Stéphane Mathis,
Conny Aerts
Abstract:
Aims: We investigate from a theoretical perspective if space asteroseismology can be used to distinguish between different thermal structures and shapes of the near-core mixing profiles for different types of coherent oscillation modes in massive stars with convective cores, and if this capacity depends on the evolutionary stage of the models along the main sequence. Methods: We compute 1D stellar…
▽ More
Aims: We investigate from a theoretical perspective if space asteroseismology can be used to distinguish between different thermal structures and shapes of the near-core mixing profiles for different types of coherent oscillation modes in massive stars with convective cores, and if this capacity depends on the evolutionary stage of the models along the main sequence. Methods: We compute 1D stellar structure and evolution models for four different prescriptions of the mixing and temperature gradient in the near-core region. Their effect on the frequencies of dipole prograde gravity modes in both Slowly Pulsating B and $β$ Cep stars is investigated, as well as for pressure modes in $β$ Cep stars. Results: A comparison between the mode frequencies of the different models at various stages during the main sequence evolution reveals that they are more sensitive to a change in temperature gradient than to the exact shape of the mixing profile in the near-core region. Depending on the duration of the observed light curve, one can distinguish between either just the temperature gradient, or also between the shapes of the mixing coefficient. The relative frequency differences are in general larger for more evolved models, and are largest for the higher-frequency pressure modes in $β$ Cep stars. Conclusions:In order to unravel the core boundary mixing and thermal structure of the near-core region, one must have asteroseismic masses and radii with $\sim 1\%$ relative precision for hundreds of stars.
△ Less
Submitted 18 August, 2019; v1 submitted 12 June, 2019;
originally announced June 2019.
-
Rossby and Magnetic Prandtl Number Scaling of Stellar Dynamos
Authors:
K. ~C. Augustson,
A. ~S. Brun,
J. Toomre
Abstract:
Rotational scaling relationships are examined for the degree of equipartition between magnetic and kinetic energies in stellar convection zones. These scaling relationships are approached from two paradigms, with first a glance at scaling relationship built upon an energy-balance argument and second a look at a force-based scaling. The latter implies a transition between a nearly-constant inertial…
▽ More
Rotational scaling relationships are examined for the degree of equipartition between magnetic and kinetic energies in stellar convection zones. These scaling relationships are approached from two paradigms, with first a glance at scaling relationship built upon an energy-balance argument and second a look at a force-based scaling. The latter implies a transition between a nearly-constant inertial scaling when in the asymptotic limit of minimal diffusion and magnetostrophy, whereas the former implies a weaker scaling with convective Rossby number. Both scaling relationships are then compared to a suite of 3D convective dynamo simulations with a wide variety of domain geometries, stratifications, and range of convective Rossby numbers.
△ Less
Submitted 30 March, 2019;
originally announced April 2019.
-
Rotating Convection and Gravito-Inertial Wave Generation in Stellar Interiors
Authors:
Kyle C. Augustson,
Stéphane Mathis
Abstract:
Gravito-inertial waves can be excited at the interface of convective and radiative regions and by the Reynolds stresses in the bulk of the convection zone. The magnitude of their energy flux will therefore vary with the properties of the convection. To assess how convection changes with rotation, a simplified local monomodal model for rotating convection is presented that provides the magnitude of…
▽ More
Gravito-inertial waves can be excited at the interface of convective and radiative regions and by the Reynolds stresses in the bulk of the convection zone. The magnitude of their energy flux will therefore vary with the properties of the convection. To assess how convection changes with rotation, a simplified local monomodal model for rotating convection is presented that provides the magnitude of the rms velocity, degree of superadiabaticity, and characteristic length scale as a function of the convective Rossby number as well as with thermal and viscous diffusivities. In the context of this convection model, two models for assessing the gravito-inertial wave flux are considered: an interfacial model and a full treatment of the Reynolds stress impact on the waves. It is found that there are regimes where the sub-inertial waves may carry a significant energy flux relative to pure gravity waves that depend upon the convective Rossby number, the ratio of the buoyancy time-scale in the stable region to the convective overturning time, and the wave frequency.
△ Less
Submitted 27 February, 2019;
originally announced February 2019.
-
A model of rotating convection in stellar and planetary interiors: I - convective penetration
Authors:
Kyle C. Augustson,
Stéphane Mathis
Abstract:
A monomodal model for stellar and planetary convection is derived for the magnitude of the rms velocity, degree of superadiabaticity, and characteristic length scale as a function of rotation rate as well as with thermal and viscous diffusivities. The convection model is used as a boundary condition for a linearization of the equations of motion in the transition region between convectively unstab…
▽ More
A monomodal model for stellar and planetary convection is derived for the magnitude of the rms velocity, degree of superadiabaticity, and characteristic length scale as a function of rotation rate as well as with thermal and viscous diffusivities. The convection model is used as a boundary condition for a linearization of the equations of motion in the transition region between convectively unstable and stably-stratified regions, yielding the depth to which convection penetrates into the stable region and establishing a relationship between that depth and the local convective Rossby number, diffusivity, and pressure scale height of those flows. Upward and downward penetrative convection have a similar scaling with rotation rate and diffusivities, but they depend differently upon the pressure scale height due to the differing energetic processes occurring in convective cores of early-type stars versus convective envelopes of late-type stars.
△ Less
Submitted 27 February, 2019;
originally announced February 2019.
-
Impact of general differential rotation on gravity waves in rapidly rotating stars
Authors:
Vincent Prat,
Stéphane Mathis,
Kyle Augustson,
François Lignières,
Jérôme Ballot,
Lucie Alvan,
Allan Sacha Brun
Abstract:
Differential rotation plays a key role in stellar evolution by triggering hydrodynamical instabilities and large-scale motions that induce transport of chemicals and angular momentum and by modifying the propagation and the frequency spectrum of gravito-inertial waves. It is thus crucial to investigate its effect on the propagation of gravity waves to build reliable seismic diagnostic tools, espec…
▽ More
Differential rotation plays a key role in stellar evolution by triggering hydrodynamical instabilities and large-scale motions that induce transport of chemicals and angular momentum and by modifying the propagation and the frequency spectrum of gravito-inertial waves. It is thus crucial to investigate its effect on the propagation of gravity waves to build reliable seismic diagnostic tools, especially for fast rotating stars, where perturbative treatments of rotation fail. Generalising a previous work done in the case of uniform rotation, we derived a local dispersion relation for gravity waves in a differentially rotating star, taking the full effect of rotation (both Coriolis and centrifugal accelerations) into account. Then we modelled the propagation of axisymmetric waves as the propagation of rays. This allowed us to efficiently probe the properties of the waves in various regimes of differential rotation.
△ Less
Submitted 7 December, 2018;
originally announced December 2018.
-
Does magnetic field modify tidal dynamics in the convective envelope of Solar mass stars?
Authors:
Aurélie Astoul,
Stéphane Mathis,
Clément Baruteau,
Florian Gallet,
Kyle C. Augustson,
Emeline Bolmont,
Allan S. Brun,
Antoine Strugarek
Abstract:
The energy dissipation of wave-like tidal flows in the convective envelope of low-mass stars is one of the key physical mechanisms that shape the orbital and rotational dynamics of short-period planetary systems. Tidal flows, and the excitation, propagation, and dissipation of tidally-induced inertial waves can be modified by stellar magnetic fields (e.g., Wei 2016, 2018, Lin and Ogilvie 2018). It…
▽ More
The energy dissipation of wave-like tidal flows in the convective envelope of low-mass stars is one of the key physical mechanisms that shape the orbital and rotational dynamics of short-period planetary systems. Tidal flows, and the excitation, propagation, and dissipation of tidally-induced inertial waves can be modified by stellar magnetic fields (e.g., Wei 2016, 2018, Lin and Ogilvie 2018). It is thus important to assess for which stars, at which location of their internal structure, and at which phase of their evolution, one needs to take into account the effects of magnetic fields on tidal waves. Using scaling laws that provide the amplitude of dynamo-generated magnetic fields along the rotational evolution of these stars (e.g., Christensen et al. 2009, Brun et al. 2015), combined with detailed grids of stellar rotation models (e.g., Amard et al. 2016), we examine the influence of a magnetic field on tidal forcing and dissipation near the tachocline of solar-like stars. We show that full consideration of magnetic fields is required to compute tidal dissipation, but not necessarily for tidal forcing.
△ Less
Submitted 21 November, 2018;
originally announced November 2018.
-
Perturbations of Stellar Oscillations: a General Magnetic Field
Authors:
Kyle Augustson,
Stéphane Mathis
Abstract:
The purpose of this short paper is to give a summary of the derivation of a perturbative method to assess the impact of an arbitrary 3D magnetic field configuration, whose energy is everywhere small relative to the gravitational binding energy, on the form and value of the eigenfrequencies and eigenfunctions of stellar oscillations.
The purpose of this short paper is to give a summary of the derivation of a perturbative method to assess the impact of an arbitrary 3D magnetic field configuration, whose energy is everywhere small relative to the gravitational binding energy, on the form and value of the eigenfrequencies and eigenfunctions of stellar oscillations.
△ Less
Submitted 11 November, 2018;
originally announced November 2018.
-
Penetration of Rotating Convection
Authors:
Kyle Augustson,
Stéphane Mathis
Abstract:
A simplified model for stellar and planetary convection is derived for the magnitude of the rms velocity, degree of superadiabaticity, and characteristic length scale with Rossby number as well as with thermal and viscous diffusivities. Integrating the convection model into a linearization of the dynamics in the transition region between convectively unstable and stably-stratified region yields a…
▽ More
A simplified model for stellar and planetary convection is derived for the magnitude of the rms velocity, degree of superadiabaticity, and characteristic length scale with Rossby number as well as with thermal and viscous diffusivities. Integrating the convection model into a linearization of the dynamics in the transition region between convectively unstable and stably-stratified region yields a Rossby number, diffusivity, and pressure scale height dependent convective penetration depth into the stable region. This may have important consequences for mixing along the evolution of rotating stars.
△ Less
Submitted 11 November, 2018;
originally announced November 2018.
-
Waves in the radiative zones of rotating, magnetized stars
Authors:
Aurélien Valade,
Vincent Prat,
Stéphane Mathis,
Kyle Augustson
Abstract:
Asteroseismology has reached a level of accuracy that may allow us to detect the effect of a deep magnetic field on oscillation modes. We thus aim to develop an asymptotic theory for short-wavelength waves in radiative zones of rotating stars with a general magnetic field topology (toroidal and poloidal). A dispersion relation is derived for a uniformly rotating star with a low-amplitude, axisymme…
▽ More
Asteroseismology has reached a level of accuracy that may allow us to detect the effect of a deep magnetic field on oscillation modes. We thus aim to develop an asymptotic theory for short-wavelength waves in radiative zones of rotating stars with a general magnetic field topology (toroidal and poloidal). A dispersion relation is derived for a uniformly rotating star with a low-amplitude, axisymmetric magnetic field. The parameter space of this model is explored with a Hamiltonian ray-tracing method. The features of the gravito-inertial waves modified by the magnetic field are studied with different magnetic topologies.
△ Less
Submitted 12 October, 2018; v1 submitted 4 October, 2018;
originally announced October 2018.
-
Asymptotic theory of gravity modes in rotating stars II. Impact of general differential rotation
Authors:
Vincent Prat,
Stéphane Mathis,
Kyle Augustson,
François Lignières,
Jérôme Ballot,
Lucie Alvan,
Allan Sacha Brun
Abstract:
Context. Differential rotation has a strong influence on stellar internal dynamics and evolution, notably by triggering hydrodynamical instabilities, by interacting with the magnetic field, and more generally by inducing transport of angular momentum and chemical elements. Moreover, it modifies the way waves propagate in stellar interiors and thus the frequency spectrum of these waves, the regions…
▽ More
Context. Differential rotation has a strong influence on stellar internal dynamics and evolution, notably by triggering hydrodynamical instabilities, by interacting with the magnetic field, and more generally by inducing transport of angular momentum and chemical elements. Moreover, it modifies the way waves propagate in stellar interiors and thus the frequency spectrum of these waves, the regions they probe, and the transport they generate.
Aims. We investigate the impact of a general differential rotation (both in radius and latitude) on the propagation of axisymmetric gravito-inertial waves.
Methods. We use a small-wavelength approximation to obtain a local dispersion relation for these waves. We then describe the propagation of waves thanks to a ray model that follows a Hamiltonian formalism. Finally, we numerically probe the properties of these gravito-inertial rays for different regimes of radial and latitudinal differential rotation.
Results. We derive a local dispersion relation that includes the effect of a general differential rotation. Subsequently, considering a polytropic stellar model, we observe that differential rotation allows for a large variety of resonant cavities that can be probed by gravito-inertial waves. We identify that for some regimes of frequency and differential rotation, the properties of gravito-inertial rays are similar to those found in the uniformly rotating case. Furthermore, we also find new regimes specific to differential rotation, where the dynamics of rays is chaotic.
Conclusions. As a consequence, we expect modes to follow the same trend. Some parts of oscillation spectra corresponding to regimes similar to those of the uniformly rotating case would exhibit regular patterns, while parts corresponding to the new regimes would be mostly constituted of chaotic modes with a spectrum rather characterised by a generic statistical distribution.
△ Less
Submitted 21 March, 2018; v1 submitted 12 March, 2018;
originally announced March 2018.
-
On differential rotation and overshooting in solar-like stars
Authors:
Allan Sacha Brun,
Antoine Strugarek,
Jacobo Varela,
Sean P. Matt,
Kyle C. Augustson,
Constance Emeriau,
Olivier Long DoCao,
Benjamin Brown,
Juri Toomre
Abstract:
We seek to characterize how the change of global rotation rate influences the overall dynamics and large scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like…
▽ More
We seek to characterize how the change of global rotation rate influences the overall dynamics and large scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like stars spin down over the course of their main sequence evolution, such change must have a direct impact on their dynamics and rotation state. We indeed find that three main states of rotation may exist for a given star: anti-solar-like (fast poles, slow equator), solar-like (fast equator, slow poles), or a cylindrical rotation profile. Under increasingly strict rotational constraints, the latter profile can further evolve into a Jupiter-like profile, with alternating prograde and retrograde zonal jets. We have further assessed how far the convection and meridional flows overshoot into the radiative zone and investigated the morphology of the established tachocline. Using simple mixing length arguments, we are able to construct a scaling of the fluid Rossby number $R_{of} = \tildeω/2Ω_* \sim \tilde{v}/2Ω_* R_*$, which we calibrate based on our 3-D ASH simulations. We can use this scaling to map the behavior of differential rotation versus the global parameters of stellar mass and rotation rate. Finally, we isolate a region on this map ($R_{of} \gtrsim 1.5-2$) where we posit that stars with an anti-solar differential rotation may exist in order to encourage observers to hunt for such targets.
△ Less
Submitted 21 February, 2017;
originally announced February 2017.
-
Dynamo Scaling Relationships
Authors:
Kyle Augustson,
Stéphane Mathis,
Allan Sacha Brun,
Juri Toomre
Abstract:
This paper provides a brief look at dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Two simple models are examined, where one that assumes magnetostrophy and another that includes the effects of inertia. These models are then compared to a suite of convective dynamo simulations of the convective core of a main-sequence B-type star and applied to…
▽ More
This paper provides a brief look at dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Two simple models are examined, where one that assumes magnetostrophy and another that includes the effects of inertia. These models are then compared to a suite of convective dynamo simulations of the convective core of a main-sequence B-type star and applied to its later evolutionary stages.
△ Less
Submitted 14 February, 2017;
originally announced February 2017.
-
Dynamos and Differential Rotation: Advances at the Crossroads of Analytics, Numerics, and Observations
Authors:
Kyle Augustson
Abstract:
The recent observational, theoretical, and numerical progress made in understanding stellar magnetism is discussed. Particularly, this review will cover the physical processes thought to be at the origin of these magnetic fields and their variability, namely dynamo action arising from the interaction between convection, rotation, radiation and magnetic fields. Some care will be taken to cover rece…
▽ More
The recent observational, theoretical, and numerical progress made in understanding stellar magnetism is discussed. Particularly, this review will cover the physical processes thought to be at the origin of these magnetic fields and their variability, namely dynamo action arising from the interaction between convection, rotation, radiation and magnetic fields. Some care will be taken to cover recent analytical advances regarding the dynamics and magnetism of radiative interiors, including some thoughts on the role of a tachocline. Moreover, recent and rapidly advancing numerical modeling of convective dynamos will be discussed, looking at rapidly rotating convective systems, grand minima and scaling laws for magnetic field strength. These topics are linked to observations or their observational implications.
△ Less
Submitted 10 January, 2017;
originally announced January 2017.
-
Simple Scaling Relationships for Stellar Dynamos
Authors:
Kyle Augustson,
Stéphane Mathis,
Allan Sacha Brun
Abstract:
This paper provides a brief overview of dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Three basic approaches are adopted to explore these scaling relationships, with a first look at two simple models: one assuming magnetostrophy and another that includes the effects of inertia. Next, a third scaling relationship is derived that utilizes the ass…
▽ More
This paper provides a brief overview of dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Three basic approaches are adopted to explore these scaling relationships, with a first look at two simple models: one assuming magnetostrophy and another that includes the effects of inertia. Next, a third scaling relationship is derived that utilizes the assumptions that the dynamo possesses two integral spatial scales and that it is driven by the balance of buoyancy work and ohmic dissipation as studied in Davidson 2013. The results of which are then compared to a suite of convective dynamo simulations that possess a fully convective domain with a weak density stratification and that captured the behavior of the resulting dynamo for a range of convective Rossby numbers (Augustson et al. 2016).
△ Less
Submitted 10 January, 2017;
originally announced January 2017.
-
On the Stability of a General Magnetic Field Topology in Stellar Radiative Zones
Authors:
Kyle Augustson,
Stephane Mathis,
Antoine Strugarek
Abstract:
This paper provides a brief overview of the formation of stellar fossil magnetic fields and what potential instabilities may occur given certain configurations of the magnetic field. In particular, a purely magnetic instability can occur for poloidal, toroidal, and mixed poloidal-toroidal axisymmetric magnetic field configurations as originally studied in Tayler (1973), Markey & Tayler (1973), and…
▽ More
This paper provides a brief overview of the formation of stellar fossil magnetic fields and what potential instabilities may occur given certain configurations of the magnetic field. In particular, a purely magnetic instability can occur for poloidal, toroidal, and mixed poloidal-toroidal axisymmetric magnetic field configurations as originally studied in Tayler (1973), Markey & Tayler (1973), and Tayler (1980). However, most of the magnetic field configurations observed at the surface of massive stars are non- axisymmetric. Thus, extending earlier studies of the axisymmetric Tayler instability in spherical geometry (Goossens, 1980), we introduce a formulation for the global change in the potential energy contained in a convectively-stable region given an arbitrary Lagrangian perturbation, which permits the inclusion of both axisymmetric and non-axisymmetric magnetic fields. With this tool in hand, a path is shown by which more general stability criterion can be established.
△ Less
Submitted 27 October, 2016;
originally announced October 2016.
-
Magnetic Helicity Reversals in a Cyclic Convective Dynamo
Authors:
Mark S. Miesch,
Mei Zhang,
Kyle C. Augustson
Abstract:
We investigate the role of magnetic helicity in promoting cyclic magnetic activity in a global, 3D, magnetohydrodynamic (MHD) simulation of a convective dynamo. This simulation is characterized by coherent bands of toroidal field that exist within the convection zone, with opposite polarities in the northern and southern hemispheres. Throughout most of the cycle, the magnetic helicity in these ban…
▽ More
We investigate the role of magnetic helicity in promoting cyclic magnetic activity in a global, 3D, magnetohydrodynamic (MHD) simulation of a convective dynamo. This simulation is characterized by coherent bands of toroidal field that exist within the convection zone, with opposite polarities in the northern and southern hemispheres. Throughout most of the cycle, the magnetic helicity in these bands is negative in the northern hemisphere and positive in the southern hemisphere. However, during the declining phase of each cycle, this hemispheric rule reverses. We attribute this to a global restructuring of the magnetic topology that is induced by the interaction of the bands across the equator. This band interaction appears to be ultimately responsible for, or at least associated with, the decay and subsequent reversal of both the toroidal bands and the polar fields. We briefly discuss the implications of these results within the context of solar observations, which also show some potential evidence for toroidal band interactions and helicity reversals.
△ Less
Submitted 31 May, 2016;
originally announced May 2016.
-
Velocity Amplitudes in Global Convection Simulations: The Role of the Prandtl Number and Near-Surface Driving
Authors:
Bridget O'Mara,
Mark S. Miesch,
Nicholas A. Featherstone,
K. C. Augustson
Abstract:
Several lines of evidence suggest that the velocity amplitude in global simulations of solar convection, U, may be systematically over-estimated. Motivated by these recent results, we explore the factors that determine U and we consider how these might scale to solar parameter regimes. To this end, we decrease the thermal diffusivity $κ$ along two paths in parameter space. If the kinematic viscosi…
▽ More
Several lines of evidence suggest that the velocity amplitude in global simulations of solar convection, U, may be systematically over-estimated. Motivated by these recent results, we explore the factors that determine U and we consider how these might scale to solar parameter regimes. To this end, we decrease the thermal diffusivity $κ$ along two paths in parameter space. If the kinematic viscosity $ν$ is decreased proportionally with $κ$ (fixing the Prandtl number $P_r = ν/κ$), we find that U increases but asymptotes toward a constant value, as found by Featherstone & Hindman (2016). However, if $ν$ is held fixed while decreasing $κ$ (increasing $P_r$), we find that U systematically decreases. We attribute this to an enhancement of the thermal content of downflow plumes, which allows them to carry the solar luminosity with slower flow speeds. We contrast this with the case of Rayleigh-Benard convection which is not subject to this luminosity constraint. This dramatic difference in behavior for the two paths in parameter space (fixed $P_r$ or fixed $ν$) persists whether the heat transport by unresolved, near-surface convection is modeled as a thermal conduction or as a fixed flux. The results suggest that if solar convection can operate in a high-$P_r$ regime, then this might effectively limit the velocity amplitude. Small-scale magnetism is a possible source of enhanced viscosity that may serve to achieve this high-$P_r$ regime.
△ Less
Submitted 19 March, 2016;
originally announced March 2016.
-
The Magnetic Furnace: Intense Core Dynamos in B-stars
Authors:
Kyle C. Augustson,
Allan Sacha Brun,
Juri Toomre
Abstract:
The dynamo action achieved in the convective cores of main-sequence massive stars is explored here through 3D global simulations of convective core dynamos operating within a young $\Msun{10}$ B-type star, using the anelastic spherical harmonic (ASH) code. These simulations capture the inner 65\% of this star by radius, encompassing the convective nuclear-burning core (about 23\% by radius) and a…
▽ More
The dynamo action achieved in the convective cores of main-sequence massive stars is explored here through 3D global simulations of convective core dynamos operating within a young $\Msun{10}$ B-type star, using the anelastic spherical harmonic (ASH) code. These simulations capture the inner 65\% of this star by radius, encompassing the convective nuclear-burning core (about 23\% by radius) and a portion of the overlying radiative envelope. Eight rotation rates are considered, ranging from 0.05\% to 16\% of the surface breakup velocity, thereby capturing both convection that barely senses the effects of rotation and other situations in which the Coriolis forces are prominent. The vigorous dynamo action realized within all of these turbulent convective cores builds magnetic fields with peak strengths exceeding a megagauss, with the overall magnetic energy (ME) in the faster rotators reaching super-equipartition levels compared to the convective kinetic energy (KE). The core convection typically involves turbulent columnar velocity structures roughly aligned with the rotation axis, with magnetic fields threading through these rolls and possessing complex linkages throughout the core. The very strong fields are able to coexist with the flows without quenching them through Lorentz forces. The velocity and magnetic fields achieve such a state by being nearly co-aligned, and with peak magnetic islands being somewhat displaced from the fastest flows as the intricate evolution proceeds. As the rotation rate is increased, the primary force balance shifts from nonlinear advection balancing Lorentz forces to a magnetostrophic balance between Coriolis and Lorentz forces.
△ Less
Submitted 1 July, 2016; v1 submitted 11 March, 2016;
originally announced March 2016.
-
Convective Dynamo Simulation with a Grand Minimum
Authors:
Kyle Augustson,
Sacha Brun,
Mark Miesch,
Juri Toomre
Abstract:
The global-scale dynamo action achieved in a simulation of a Sun-like star rotating at thrice the solar rate is assessed. The 3-D MHD Anelastic Spherical Harmonic (ASH) code, augmented with a viscosity minimization scheme, is employed to capture convection and dynamo processes in this G-type star. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the…
▽ More
The global-scale dynamo action achieved in a simulation of a Sun-like star rotating at thrice the solar rate is assessed. The 3-D MHD Anelastic Spherical Harmonic (ASH) code, augmented with a viscosity minimization scheme, is employed to capture convection and dynamo processes in this G-type star. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. It is found that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals occurring roughly every 6.2 years. The magnetic energy also rises and falls with a regular period. The magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the processes leading to polarity reversals are more complex, appearing to arise from the interaction of convection with the mean toroidal fields. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation, and potentially to a nonlinear dynamo wave. This simulation also enters a grand minimum lasting roughly 20~years, after which the dynamo recovers its regular polarity cycles.
△ Less
Submitted 13 March, 2015;
originally announced March 2015.
-
Grand Minima and Equatorward Propagation in a Cycling Stellar Convective Dynamo
Authors:
Kyle Augustson,
Allan Sacha Brun,
Mark Miesch,
Juri Toomre
Abstract:
The 3-D magnetohydrodynamic (MHD) Anelastic Spherical Harmonic (ASH) code, using slope-limited diffusion, is employed to capture convective and dynamo processes achieved in a global-scale stellar convection simulation for a model solar-mass star rotating at three times the solar rate. The dynamo generated magnetic fields possesses many time scales, with a prominent polarity cycle occurring roughly…
▽ More
The 3-D magnetohydrodynamic (MHD) Anelastic Spherical Harmonic (ASH) code, using slope-limited diffusion, is employed to capture convective and dynamo processes achieved in a global-scale stellar convection simulation for a model solar-mass star rotating at three times the solar rate. The dynamo generated magnetic fields possesses many time scales, with a prominent polarity cycle occurring roughly every 6.2 years. The magnetic field forms large-scale toroidal wreaths, whose formation is tied to the low Rossby number of the convection in this simulation. The polarity reversals are linked to the weakened differential rotation and a resistive collapse of the large-scale magnetic field. An equatorial migration of the magnetic field is seen, which is due to the strong modulation of the differential rotation rather than a dynamo wave. A poleward migration of magnetic flux from the equator eventually leads to the reversal of the polarity of the high-latitude magnetic field. This simulation also enters an interval with reduced magnetic energy at low latitudes lasting roughly 16 years (about 2.5 polarity cycles), during which the polarity cycles are disrupted and after which the dynamo recovers its regular polarity cycles. An analysis of this grand minimum reveals that it likely arises through the interplay of symmetric and antisymmetric dynamo families. This intermittent dynamo state potentially results from the simulations relatively low magnetic Prandtl number. A mean-field-based analysis of this dynamo simulation demonstrates that it is of the $α$-$Ω$ type. The time scales that appear to be relevant to the magnetic polarity reversal are also identified.
△ Less
Submitted 10 July, 2015; v1 submitted 23 October, 2014;
originally announced October 2014.
-
Dynamo modeling of the Kepler F star KIC 12009504
Authors:
S. Mathur,
K. C. Augustson,
A. S. Brun,
R. A. Garcia,
T. S. Metcalfe
Abstract:
The Kepler mission has collected light curves for almost 4 years. The excellent quality of these data has allowed us to probe the structure and the dynamics of the stars using asteroseismology. With the length of data available, we can start to look for magnetic activity cycles. The Kepler data obtained for the F star, KIC 12009504, shows a rotation period of 9.5 days and additional variability th…
▽ More
The Kepler mission has collected light curves for almost 4 years. The excellent quality of these data has allowed us to probe the structure and the dynamics of the stars using asteroseismology. With the length of data available, we can start to look for magnetic activity cycles. The Kepler data obtained for the F star, KIC 12009504, shows a rotation period of 9.5 days and additional variability that could be due to the magnetic activity of the star. Here we present recent and preliminary 3D global-scale dynamo simulations of this star with the ASH and STELEM codes, capturing a substantial portion of the convection and the stable radiation zone below it. These simulations reveal a multi-year activity cycle whose length tentatively depends upon the width of the tachocline present in the simulation. Furthermore, the presence of a magnetic field and the dynamo action taking place in the convection zone appears to help confine the tachocline, but longer simulations will be required to confirm this.
△ Less
Submitted 25 August, 2014;
originally announced August 2014.
-
Cycling Dynamo in a Young Sun: Grand Minima and Equatorward Propagation
Authors:
Kyle Augustson,
Allan Sacha Brun,
Mark Steven Miesch,
Juri Toomre
Abstract:
We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zon…
▽ More
We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. We find that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals every 6.2 years. The magnetic energy also rises and falls with a regular period, with two magnetic energy cycles required to complete a polarity cycle. These magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the polarity reversals are present due to the spatial separation of the equatorial and polar dynamos. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation and to a near-surface shear layer. This simulation also enters a grand minimum lasting roughly 20 years, after which the dynamo recovers its regular polarity cycles.
△ Less
Submitted 31 October, 2013;
originally announced October 2013.
-
Effects of Granulation upon Larger-Scale Convection
Authors:
Neal Hurlburt,
Marc DeRosa,
Kyle Augustson,
Juri Toomre
Abstract:
We examine the role of small-scale granulation in helping to drive supergranulation and even larger scales of convection. The granulation is modeled as localized cooling events introduced at the upper boundary of a 3-D simulation of compressible convection in a rotating spherical shell segment. With a sufficient number of stochastic cooling events compared to uniform cooling, we find that supergra…
▽ More
We examine the role of small-scale granulation in helping to drive supergranulation and even larger scales of convection. The granulation is modeled as localized cooling events introduced at the upper boundary of a 3-D simulation of compressible convection in a rotating spherical shell segment. With a sufficient number of stochastic cooling events compared to uniform cooling, we find that supergranular scales are realized, along with a differential rotation that becomes increasingly solar-like.
△ Less
Submitted 23 January, 2012;
originally announced January 2012.
-
Modeling the Near-Surface Shear Layer: Diffusion Schemes Studied With CSS
Authors:
Kyle Augustson,
Mark Rast,
Regner Trampedach,
Juri Toomre
Abstract:
As we approach solar convection simulations that seek to model the interaction of small-scale granulation and supergranulation and even larger scales of convection within the near-surface shear layer (NSSL), the treatment of the boundary conditions and minimization of sub-grid scale diffusive processes become increasingly crucial. We here assess changes in the dynamics and the energy flux balance…
▽ More
As we approach solar convection simulations that seek to model the interaction of small-scale granulation and supergranulation and even larger scales of convection within the near-surface shear layer (NSSL), the treatment of the boundary conditions and minimization of sub-grid scale diffusive processes become increasingly crucial. We here assess changes in the dynamics and the energy flux balance of the flows established in rotating spherical shell segments that capture much of the NSSL with the Curved Spherical Segment (CSS) code using two different diffusion schemes. The CSS code is a new massively parallel modeling tool capable of simulating 3-D compressible MHD convection with a realistic solar stratification in rotating spherical shell segments.
△ Less
Submitted 21 December, 2010;
originally announced December 2010.
-
Convection and dynamo action in B stars
Authors:
Kyle C. Augustson,
Allan S. Brun,
Juri Toomre
Abstract:
Main-sequence massive stars possess convective cores that likely harbor strong dynamo action. To assess the role of core convection in building magnetic fields within these stars, we employ the 3-D anelastic spherical harmonic (ASH) code to model turbulent dynamics within a 10 solar mass main-sequence (MS) B-type star rotating at 4 times the solar rate. We find that strong (900 kG) magnetic fields…
▽ More
Main-sequence massive stars possess convective cores that likely harbor strong dynamo action. To assess the role of core convection in building magnetic fields within these stars, we employ the 3-D anelastic spherical harmonic (ASH) code to model turbulent dynamics within a 10 solar mass main-sequence (MS) B-type star rotating at 4 times the solar rate. We find that strong (900 kG) magnetic fields arise within the turbulence of the core and penetrate into the stably stratified radiative zone. These fields exhibit complex, time-dependent behavior including reversals in magnetic polarity and shifts between which hemisphere dominates the total magnetic energy.
△ Less
Submitted 3 November, 2010;
originally announced November 2010.
