-
Hot Spot Offset Variability from Magnetohydrodynamical Thermoresistive Instability in Hot Jupiters
Authors:
Raphaël Hardy,
Paul Charbonneau,
Andrew Cumming
Abstract:
Hot Jupiter atmospheres are possibly subject to a thermoresistive instability. Such an instability may develop as the ohmic heating increases the electrical conductivity in a positive feedback loop, which ultimately leads to a runaway of the atmospheric temperature. We extend our previous axisymmetric one-dimensional radial model, by representing the temperature and magnetic diffusivity as a first…
▽ More
Hot Jupiter atmospheres are possibly subject to a thermoresistive instability. Such an instability may develop as the ohmic heating increases the electrical conductivity in a positive feedback loop, which ultimately leads to a runaway of the atmospheric temperature. We extend our previous axisymmetric one-dimensional radial model, by representing the temperature and magnetic diffusivity as a first order Fourier expansion in longitude. This allows us to predict the hot spot offset during the unfolding of the thermoresistive instability and following Alfvénic oscillations. We show a representative simulation undergoing the thermoresistive instability, in which the peak flux offset varies between approximately $\pm 60^{\circ}$ on timescales of a few days with potentially observable brightness variations. Therefore, this thermoresistive instability could be an observable feature of hot Jupiters, given the right timing of observation and transit and the right planetary parameters.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
Energy Definition and Minimization in Avalanche Models for Solar Flares
Authors:
Henri Lamarre,
Paul Charbonneau,
Antoine Strugarek
Abstract:
Self-organised critical avalanche models are a class of cellular automata that, despite their simplicity, can be applied to the modeling of solar (and stellar) flares and generate robust power-law distributions in event size measures. However, bridging the conceptual gap to both magnetohydrodynamics and real flare observations continues to prove challenging. In this paper, we focus on a specific,…
▽ More
Self-organised critical avalanche models are a class of cellular automata that, despite their simplicity, can be applied to the modeling of solar (and stellar) flares and generate robust power-law distributions in event size measures. However, bridging the conceptual gap to both magnetohydrodynamics and real flare observations continues to prove challenging. In this paper, we focus on a specific, key aspect of this endeavour, namely the definition of magnetic energy and its consequences for the model's internal dynamics and energy release statistics. We show that the dual requirement of releasing energy and restoring local stability demands that the instability criterion and boundary conditions be set in a manner internally consistent with a given energy definition, otherwise unphysical behavior ensues, e.g., negative energy release. Working with three energy definitions previously used in the literature, we construct such internally consistent avalanche models and compare/contrast their energy release statistics. Using the same set of models, we also explore a recent proposal by Farhang et al. (2018, 2019), namely that avalanches/flares should maximize the amount of energy released by the lattice when instabilities are triggered. This tends to produce avalanches of shorter duration but higher peak energy release, but adding up to similar total energy release. For the three energy definition we tested, such avalanche models exhibit almost identical distributions of event size measures. Our results indicate that the key to reproduce solar-like power-law slopes in these size measures is lattice configurations in which most nodes remain relatively far from the instability threshold.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
Magnetohydrodynamical torsional oscillations from thermo-resistive instability in hot jupiters
Authors:
Raphaël Hardy,
Paul Charbonneau,
Andrew Cumming
Abstract:
Hot jupiter atmospheres may be subject to a thermo-resistive instability where an increase in the electrical conductivity due to ohmic heating results in runaway of the atmospheric temperature. We introduce a simplified one-dimensional model of the equatorial sub-stellar region of a hot jupiter which includes the temperature-dependence and time-dependence of the electrical conductivity, as well as…
▽ More
Hot jupiter atmospheres may be subject to a thermo-resistive instability where an increase in the electrical conductivity due to ohmic heating results in runaway of the atmospheric temperature. We introduce a simplified one-dimensional model of the equatorial sub-stellar region of a hot jupiter which includes the temperature-dependence and time-dependence of the electrical conductivity, as well as the dynamical back-reaction of the magnetic field on the flow. This model extends our previous one-zone model to include the radial structure of the atmosphere. Spatial gradients of electrical conductivity strongly modify the radial profile of Alfvén oscillations, leading to steepening and downwards transport of magnetic field, enhancing dissipation at depth. We find unstable solutions that lead to self-sustained oscillations for equilibrium temperatures in the range $T_\mathrm{eq}\approx 1000$--$1200$~K, and magnetic field in the range $\approx 10$--$100$~G. For a given set of parameters, self-sustained oscillations occur in a narrow range of equilibrium temperatures which allow the magnetic Reynolds number to alternate between large and small values during an oscillation cycle. Outside of this temperature window, the system reaches a steady state in which the effect of the magnetic field can be approximated as a magnetic drag term. Our results show that thermo-resistive instability is a possible source of variability in magnetized hot jupiters at colder temperatures, and emphasize the importance of including the temperature-dependence of electrical conductivity in models of atmospheric dynamics.
△ Less
Submitted 1 August, 2023;
originally announced August 2023.
-
Evolution of solar and stellar dynamo theory
Authors:
Paul Charbonneau,
Dmitry Sokoloff
Abstract:
In this paper, written as a general historical and technical introduction to the various review papers collected in the special issue ``Solar and Stellar Dynamo: A New Era'', we review the evolution and current state of dynamo theory and modelling, with emphasis on the solar dynamo. Starting with a historical survey, we then focus on a set of ``tension points'' that are still left unresolved despi…
▽ More
In this paper, written as a general historical and technical introduction to the various review papers collected in the special issue ``Solar and Stellar Dynamo: A New Era'', we review the evolution and current state of dynamo theory and modelling, with emphasis on the solar dynamo. Starting with a historical survey, we then focus on a set of ``tension points'' that are still left unresolved despite the remarkable progress of the past century. In our discussion of these tension points we touch upon the physical well-posedness of mean-field electrodynamics; constraints imposed by magnetic helicity conservation; the troublesome role of differential rotation; meridional flows and flux transpost dynamos; competing inductive mechanisms and Babcock-Leighton dynamos; the ambiguous precursor properties of the solar dipole; cycle amplitude regulation and fluctuation through nonlinear backreaction and stochastic forcing, including Grand Minima; and the promises and puzzles offered by global magnethydrodynamical numerical simulations of convection and dynamo action. We close by considering the potential bridges to be constructed between solar dynamo theory and modelling, and observations of magnetic activity in late-type stars.
△ Less
Submitted 25 May, 2023;
originally announced May 2023.
-
Solaris: A Focused Solar Polar Discovery-class Mission to achieve the Highest Priority Heliophysics Science Now
Authors:
Donald M. Hassler,
Sarah E Gibson,
Jeffrey S Newmark,
Nicholas A. Featherstone,
Lisa Upton,
Nicholeen M Viall,
J Todd Hoeksema,
Frederic Auchere,
Aaron Birch,
Doug Braun,
Paul Charbonneau,
Robin Colannino,
Craig DeForest,
Mausumi Dikpati,
Cooper Downs,
Nicole Duncan,
Heather Alison Elliott,
Yuhong Fan,
Silvano Fineschi,
Laurent Gizon,
Sanjay Gosain,
Louise Harra,
Brad Hindman,
David Berghmans,
Susan T Lepri
, et al. (11 additional authors not shown)
Abstract:
Solaris is a transformative Solar Polar Discovery-class mission concept to address crucial outstanding questions that can only be answered from a polar vantage. Solaris will image the Sun's poles from ~75 degree latitude, providing new insight into the workings of the solar dynamo and the solar cycle, which are at the foundation of our understanding of space weather and space climate. Solaris will…
▽ More
Solaris is a transformative Solar Polar Discovery-class mission concept to address crucial outstanding questions that can only be answered from a polar vantage. Solaris will image the Sun's poles from ~75 degree latitude, providing new insight into the workings of the solar dynamo and the solar cycle, which are at the foundation of our understanding of space weather and space climate. Solaris will also provide enabling observations for improved space weather research, modeling and prediction, revealing a unique, new view of the corona, coronal dynamics and CME eruptions from above.
△ Less
Submitted 18 January, 2023;
originally announced January 2023.
-
Variability from thermo-resistive instability in the atmospheres of hot jupiters
Authors:
Raphaël Hardy,
Andrew Cumming,
Paul Charbonneau
Abstract:
The atmosphere of a hot jupiter may be subject to a thermo-resistive instability, in which the increasing electrical conductivity with temperature leads to runaway Ohmic heating. We introduce a simplified model of the local dynamics in the equatorial region of a hot jupiter that incorporates the back reaction on the atmospheric flow as the increasing electrical conductivity leads to flux freezing,…
▽ More
The atmosphere of a hot jupiter may be subject to a thermo-resistive instability, in which the increasing electrical conductivity with temperature leads to runaway Ohmic heating. We introduce a simplified model of the local dynamics in the equatorial region of a hot jupiter that incorporates the back reaction on the atmospheric flow as the increasing electrical conductivity leads to flux freezing, which in turn quenches the flow and therefore the Ohmic heating. We demonstrate a new time-dependent solution that emerges for a temperature-dependent electrical conductivity (whereas a temperature-independent conductivity always evolves to a steady-state). The periodic cycle consists of bursts of Alfven oscillations separated by quiescent intervals, with the magnetic Reynolds number alternating between values smaller than and larger than unity, maintaining the oscillation. We investigate the regions of pressure and temperature in which the instability operates. For the typical equatorial accelerations seen in atmospheric models, we find instability at pressures $\sim 0.1$--$1\ {\rm bar}$ and temperatures $\approx 1300$--$1800\ {\rm K}$ for magnetic fields $\sim 10\ {\rm G}$. Unlike previous studies based on a constant wind velocity, we find that the instability is stronger for weaker magnetic fields. Our results add support to the idea that variability should be a feature of magnetized hot jupiter atmospheres, particularly at intermediate temperatures. The temperature-dependence of the electrical conductivity is an important ingredient that should be included in MHD models of hot jupiter atmospheric dynamics.
△ Less
Submitted 19 October, 2022; v1 submitted 5 August, 2022;
originally announced August 2022.
-
Forecasting Solar Flares by Data Assimilation in Sandpile Models
Authors:
Christian Thibeault,
Antoine Strugarek,
Paul Charbonneau,
Benoit Tremblay
Abstract:
The prediction of solar flares is still a significant challenge in space weather research, with no techniques currently capable of producing reliable forecasts performing significantly above climatology. In this paper, we present a flare forecasting technique using data assimilation coupled with computationally inexpensive cellular automata called sandpile models. Our data assimilation algorithm u…
▽ More
The prediction of solar flares is still a significant challenge in space weather research, with no techniques currently capable of producing reliable forecasts performing significantly above climatology. In this paper, we present a flare forecasting technique using data assimilation coupled with computationally inexpensive cellular automata called sandpile models. Our data assimilation algorithm uses the simulated annealing method to find an optimal initial condition that reproduces well an energy-release time series. We present and empirically analyze the predictive capabilities of three sandpile models, namely the Lu and Hamilton model (LH) and two deterministically-driven models (D). Despite their stochastic elements, we show that deterministically-driven models display temporal correlations between simulated events, a needed condition for data assimilation. We present our new data assimilation algorithm and demonstrate its success in assimilating synthetic observations produced by the avalanche models themselves. We then apply our method to GOES X-Ray time series for 11 active regions having generated multiple X-class flares in the course of their lifetime. We demonstrate that for such large flares, our data assimilation scheme substantially increases the success of ``All-Clear'' forecasts, as compared to model climatology.
△ Less
Submitted 27 June, 2022;
originally announced June 2022.
-
Powering Stellar Magnetism: Energy Transfers in Cyclic Dynamos of Sun-like Stars
Authors:
Allan Sacha Brun,
Antoine Strugarek,
Quentin Noraz,
Barbara Perri,
Jacobo Varela,
Kyle Augsutson,
Paul Charbonneau,
Juri Toomre
Abstract:
We use the ASH code to model the convective dynamo of solar-type stars. Based on a series of 15 3-D MHD simulations spanning 4 bins in rotation and mass, we show what mechanisms are at work in these stellar dynamos with and without magnetic cycles and how global stellar parameters affect the outcome. We also derive scaling laws for the differential rotation and magnetic field based on these simula…
▽ More
We use the ASH code to model the convective dynamo of solar-type stars. Based on a series of 15 3-D MHD simulations spanning 4 bins in rotation and mass, we show what mechanisms are at work in these stellar dynamos with and without magnetic cycles and how global stellar parameters affect the outcome. We also derive scaling laws for the differential rotation and magnetic field based on these simulations. We find a weaker trend between differential rotation and stellar rotation rate, ($ΔΩ\sim (|Ω|/Ω_{\odot})^{0.46}$) in the MHD solutions than in their HD counterpart $(|Ω|/Ω_{\odot})^{0.66})$, yielding a better agreement with the observational trends based on power laws. We find that for a fluid Rossby number between $0.15 \lesssim Ro_f \lesssim 0.65$ the solutions possess long magnetic cycle, if $Ro_f \lesssim 0.42$ a short cycle and if $Ro_f \gtrsim 1$ (anti-solar-like differential rotation) a statistically steady state. We show that short-cycle dynamos follow the classical Parker-Yoshimura rule whereas the long-cycle period ones do not. We further demonstrate that the Rossby number dependency of the large-scale surface magnetic field in the simulation ($B_{L,surf} \sim Ro_{f}^{-1.26}$) agrees better with observations ($B_{V} \sim Ro_{s}^{-1.4 \pm 0.1}$) and differs from dynamo scaling based on the global magnetic energy ($B_{bulk} \sim Ro_{f}^{-0.5}$). We also show that up to few percents of the stellar luminosity can be channelled into the star's magnetism, hence providing a large energy reservoir for possible surface eruptive events.
△ Less
Submitted 31 January, 2022;
originally announced January 2022.
-
Towards an algebraic method of solar cycle prediction II. Reducing the need for detailed input data with ARDoR
Authors:
Melinda Nagy,
Kristóf Petrovay,
Alexandre Lemerle,
Paul Charbonneau
Abstract:
An algebraic method for the reconstruction and potentially prediction of the solar dipole moment value at sunspot minimum (known to be a good predictor of the amplitude of the next solar cycle) was suggested in the first paper in this series. The method sums up the ultimate dipole moment contributions of individual active regions in a solar cycle: for this, detailed and reliable input data would i…
▽ More
An algebraic method for the reconstruction and potentially prediction of the solar dipole moment value at sunspot minimum (known to be a good predictor of the amplitude of the next solar cycle) was suggested in the first paper in this series. The method sums up the ultimate dipole moment contributions of individual active regions in a solar cycle: for this, detailed and reliable input data would in principle be needed for thousands of active regions in a solar cycle. To reduce the need for detailed input data, here we propose a new active region descriptor called ARDoR (Active Region Degree of Rogueness). In a detailed statistical analysis of a large number of activity cycles simulated with the 2$\times$2D dynamo model we demonstrate that ranking active regions by decreasing ARDoR, for a good reproduction of the solar dipole moment at the end of the cycle it is sufficient to consider the top $N$ regions on this list explicitly, where $N$ is a relatively low number, while for the other regions the ARDoR value may be set to zero. E.g., with $N=5$ the fraction of cycles where the dipole moment is reproduced with an error exceeding $\pm$30% is only 12%, significantly reduced with respect to the case $N=0$, i.e. ARDoR set to zero for all active regions, where this fraction is 26%. This indicates that stochastic effects on the intercycle variations of solar activity are dominated by the effect of a low number of large ``rogue'' active regions, rather than the combined effect of numerous small ARs. The method has a potential for future use in solar cycle prediction.
△ Less
Submitted 4 September, 2020;
originally announced September 2020.
-
Principles Of Heliophysics: a textbook on the universal processes behind planetary habitability
Authors:
Karel Schrijver,
Fran Bagenal,
Tim Bastian,
Juerg Beer,
Mario Bisi,
Tom Bogdan,
Steve Bougher,
David Boteler,
Dave Brain,
Guy Brasseur,
Don Brownlee,
Paul Charbonneau,
Ofer Cohen,
Uli Christensen,
Tom Crowley,
Debrah Fischer,
Terry Forbes,
Tim Fuller-Rowell,
Marina Galand,
Joe Giacalone,
George Gloeckler,
Jack Gosling,
Janet Green,
Nick Gross,
Steve Guetersloh
, et al. (37 additional authors not shown)
Abstract:
Heliophysics is the system science of the physical connections between the Sun and the solar system. As the physics of the local cosmos, it embraces space weather and planetary habitability. The wider view of comparative heliophysics forms a template for conditions in exoplanetary systems and provides a view over time of the aging Sun and its magnetic activity, of the heliosphere in different sett…
▽ More
Heliophysics is the system science of the physical connections between the Sun and the solar system. As the physics of the local cosmos, it embraces space weather and planetary habitability. The wider view of comparative heliophysics forms a template for conditions in exoplanetary systems and provides a view over time of the aging Sun and its magnetic activity, of the heliosphere in different settings of the interstellar medium and subject to stellar impacts, of the space physics over evolving planetary dynamos, and of the long-term influence on planetary atmospheres by stellar radiation and wind.
Based on a series of NASA-funded summer schools for early-career researchers, this textbook is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences. The book emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it could thrive. The text includes 200 "Activities" in the form of exercises, explorations, literature readings, "what if" challenges, and group discussion topics; many of the Activities provide additional information complementing the main text. Solutions and discussions are included in an Appendix for a selection of the exercises.
△ Less
Submitted 27 March, 2024; v1 submitted 30 October, 2019;
originally announced October 2019.
-
Impact of rogue active regions on hemispheric asymmetry
Authors:
Melinda Nagy,
Alexandre Lemerle,
Paul Charbonneau
Abstract:
The solar dipole moment at activity minimum is a good predictor of the strength of the subsequent solar cycle. Through a systematic analysis using a state-of-the-art 2$\times$2 D solar dynamo model, we found that bipolar magnetic regions (BMR) with atypical characteristics can modify the strength of the next cycle via their impact on the buildup of the dipole moment as a sunspot cycle unfolds. In…
▽ More
The solar dipole moment at activity minimum is a good predictor of the strength of the subsequent solar cycle. Through a systematic analysis using a state-of-the-art 2$\times$2 D solar dynamo model, we found that bipolar magnetic regions (BMR) with atypical characteristics can modify the strength of the next cycle via their impact on the buildup of the dipole moment as a sunspot cycle unfolds. In addition to summarizing these results, we present further effects of such "rogue" BMRs. These have the ability to generate hemispheric asymmetry in the subsequent sunspot cycle, since they modify the polar cap flux asymmetry of the ongoing cycle. We found strong correlation between the polar cap flux asymmetry of cycle i and the total pseudo sunspot number asymmetry of cycle i + 1 . Good correlation also appears in the case of the time lag of the hemispheres of cycle i + 1.
△ Less
Submitted 17 September, 2019;
originally announced September 2019.
-
On the sensitivity of magnetic cycles in global simulations of solar-like stars
Authors:
A. Strugarek,
P. Beaudoin,
P. Charbonneau,
A. S. Brun
Abstract:
The periods of magnetic activity cycles in the Sun and solar-type stars do not exhibit a simple or even single trend with respect to rotation rate or luminosity. Dynamo models can be used to interpret this diversity, and can ultimately help us understand why some solar-like stars do not exhibit a magnetic cycle, whereas some do, and for the latter what physical mechanisms set their magnetic cycle…
▽ More
The periods of magnetic activity cycles in the Sun and solar-type stars do not exhibit a simple or even single trend with respect to rotation rate or luminosity. Dynamo models can be used to interpret this diversity, and can ultimately help us understand why some solar-like stars do not exhibit a magnetic cycle, whereas some do, and for the latter what physical mechanisms set their magnetic cycle period. Three-dimensional non-linear magnetohydrodynamical simulations present the advantage of having only a small number of tunable parameters, and produce in a dynamically self-consistent manner the flows and the dynamo magnetic fields pervading stellar interiors. We conducted a series of such simulations within the EULAG-MHD framework, varying the rotation rate and luminosity of the modeled solar-like convective envelopes. We find decadal magnetic cycles when the Rossby number near the base of the convection zone is moderate (typically between 0.25 and 1). Secondary, shorter cycles located at the top of the convective envelope close to the equator are also observed in our numerical experiments, when the local Rossby number is lower than 1. The deep-seated dynamo sustained in these numerical experiments is fundamentally non-linear, in that it is the feedback of the large-scale magnetic field on the large-scale differential rotation that sets the magnetic cycle period. The cycle period is found to decrease with the Rossby number, which offers an alternative theoretical explanation to the variety of activity cycles observed in solar-like stars.
△ Less
Submitted 25 June, 2018;
originally announced June 2018.
-
The Effect of "Rogue" Active Regions on the Solar Cycle
Authors:
Melinda Nagy,
Alexandre Lemerle,
François Labonville,
Kristóf Petrovay,
Paul Charbonneau
Abstract:
The origin of cycle-to-cycle variations in solar activity is currently the focus of much interest. It has recently been pointed out that large individual active regions with atypical properties can have a significant impact on the long term behaviour of solar activity. We investigate this possibility in more detail using a recently developed 2$\times$2D dynamo model of the solar magnetic cycle. We…
▽ More
The origin of cycle-to-cycle variations in solar activity is currently the focus of much interest. It has recently been pointed out that large individual active regions with atypical properties can have a significant impact on the long term behaviour of solar activity. We investigate this possibility in more detail using a recently developed 2$\times$2D dynamo model of the solar magnetic cycle. We find that even a single "rogue" bipolar magnetic region (BMR) in the simulations can have a major effect on the further development of solar activity cycles, boosting or suppressing the amplitude of subsequent cycles. In extreme cases an individual BMR can completely halt the dynamo, triggering a grand minimum. Rogue BMRs also have the potential to induce significant hemispheric asymmetries in the solar cycle. To study the effect of rogue BMRs in a more systematic manner, a series of dynamo simulations were conducted, in which a large test BMR was manually introduced in the model at various phases of cycles of different amplitudes. BMRs emerging in the rising phase of a cycle can modify the amplitude of the ongoing cycle while BMRs emerging in later phases will only impact subsequent cycles. In this model, the strongest impact on the subsequent cycle occurs when the rogue BMR emerges around cycle maximum at low latitudes but the BMR does not need to be strictly cross-equatorial. Active regions emerging as far as $20^{\circ}$ from the equator can still have a significant impact. We demonstrate that the combined effect of the magnetic flux, tilt angle and polarity separation of the BMR on the dynamo is via their contribution to the dipole moment, $δD_{\mathrm{BMR}}$. Our results indicate that prediction of the amplitude, starting epoch and duration of a cycle requires an accurate accounting of a broad range of active regions emerging in the previous cycle.
△ Less
Submitted 6 December, 2017;
originally announced December 2017.
-
Reconciling solar and stellar magnetic cycles with nonlinear dynamo simulations
Authors:
A. Strugarek,
P. Beaudoin,
P. Charbonneau,
A. S. Brun,
J. -D. do Nascimento Jr
Abstract:
The magnetic fields of solar-type stars are observed to cycle over decadal periods -11 years in the case of the Sun. The fields originate in the turbulent convective layers of stars and have a complex dependency upon stellar rotation rate. We have performed a set of turbulent global simulations that exhibit magnetic cycles varying systematically with stellar rotation and luminosity. We find that t…
▽ More
The magnetic fields of solar-type stars are observed to cycle over decadal periods -11 years in the case of the Sun. The fields originate in the turbulent convective layers of stars and have a complex dependency upon stellar rotation rate. We have performed a set of turbulent global simulations that exhibit magnetic cycles varying systematically with stellar rotation and luminosity. We find that the magnetic cycle period is inversely proportional to the Rossby number, which quantifies the influence of rotation on turbulent convection. The trend relies on a fundamentally non-linear dynamo process and is compatible with the Sun's cycle and those of other solar-type stars.
△ Less
Submitted 13 July, 2017;
originally announced July 2017.
-
Confinement of the solar tachocline by a cyclic dynamo magnetic field
Authors:
R. Barnabé,
A. Strugarek,
P. Charbonneau,
A. S. Brun,
J. -P. Zahn
Abstract:
The surprising thinness of the solar tachocline is still not understood with certainty today. Among the numerous possible scenarios suggested to explain its radial confinement, one hypothesis is based on Maxwell stresses that are exerted by the cyclic dynamo magnetic field of the Sun penetrating over a skin depth below the turbulent convection zone. Our goal is to assess under which conditions (tu…
▽ More
The surprising thinness of the solar tachocline is still not understood with certainty today. Among the numerous possible scenarios suggested to explain its radial confinement, one hypothesis is based on Maxwell stresses that are exerted by the cyclic dynamo magnetic field of the Sun penetrating over a skin depth below the turbulent convection zone. Our goal is to assess under which conditions (turbulence level in the tachocline, strength of the dynamo-generated field, spreading mechanism) this scenario can be realized in the solar tachocline. We develop a simplified 1D model of the upper tachocline under the influence of an oscillating magnetic field imposed from above. The turbulent transport is parametrized with enhanced turbulent diffusion (or anti-diffusion) coefficients. Two main processes that thicken the tachocline are considered; either turbulent viscous spreading or radiative spreading. An extensive parameter study is carried out to establish the physical parameter regimes under which magnetic confinement of the tachocline that is due to a surface dynamo field can be realized. We have explored a large range of magnetic field amplitudes, viscosities, ohmic diffusivities and thermal diffusivities. We find that, for large but still realistic magnetic field strengths, the differential rotation can be suppressed in the upper radiative zone (and hence the tachocline confined) if weak turbulence is present (with an enhanced ohmic diffusivity of $η> 10^{7-8} \, cm^2/s$), even in the presence of radiative spreading. Our results show that a dynamo magnetic field can, in the presence of weak turbulence, prevent the inward burrowing of a tachocline subject to viscous diffusion or radiative spreading.
△ Less
Submitted 7 March, 2017;
originally announced March 2017.
-
Meridional Circulation Dynamics in a Cyclic Convective Dynamo
Authors:
D. Passos,
M. Miesch,
G. Guerrero,
P. Charbonneau
Abstract:
Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to the slow down…
▽ More
Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to the slow down of the axisymmetric poleward horizontal component. The behavior of this large scale flow deep inside the convection zone remains largely unknown. Present helioseismic techniques are not sensitive enough to capture the dynamics of this weak large scale flow. Moreover, the large time of integration needed to map the meridional circulation inside the convection zone, also masks some of the possible dynamics on shorter timescales. In this work we examine the dynamics of the meridional circulation that emerges from a 3D MHD global simulation of the solar convection zone. Our aim is to assess and quantify the behavior of meridional circulation deep inside the convection zone, where the cyclic large-scale magnetic field can reach considerable strength. Our analyses indicate that the meridional circulation morphology and amplitude are both highly influenced by the magnetic field, via the impact of magnetic torques on the global angular momentum distribution. A dynamic feature induced by these magnetic torques is the development of a prominent upward flow at mid latitudes in the lower convection zone that occurs near the equatorward edge of the toroidal bands and that peaks during cycle maximum. Globally, the dynamo-generated large-scale magnetic field drives variations in the meridional flow, in stark contrast to the conventional kinematic flux transport view of the magnetic field being advected passively by the flow.
△ Less
Submitted 18 July, 2017; v1 submitted 8 February, 2017;
originally announced February 2017.
-
A coupled $2\times2$D Babcock-Leighton solar dynamo model. II. Reference dynamo solutions
Authors:
Alexandre Lemerle,
Paul Charbonneau
Abstract:
In this paper we complete the presentation of a new hybrid $2\times2$D flux transport dynamo (FTD) model of the solar cycle based on the Babcock-Leighton mechanism of poloidal magnetic field regeneration via the surface decay of bipolar magnetic regions (BMRs). This hybrid model is constructed by allowing the surface flux transport (SFT) simulation described in Lemerle et al. 2015 to provide the p…
▽ More
In this paper we complete the presentation of a new hybrid $2\times2$D flux transport dynamo (FTD) model of the solar cycle based on the Babcock-Leighton mechanism of poloidal magnetic field regeneration via the surface decay of bipolar magnetic regions (BMRs). This hybrid model is constructed by allowing the surface flux transport (SFT) simulation described in Lemerle et al. 2015 to provide the poloidal source term to an axisymmetric FTD simulation defined in a meridional plane, which in turn generates the BMRs required by the SFT. A key aspect of this coupling is the definition of an emergence function describing the probability of BMR emergence as a function of the spatial distribution of the internal axisymmetric magnetic field. We use a genetic algorithm to calibrate this function, together with other model parameters, against observed cycle 21 emergence data. We present a reference dynamo solution reproducing many solar cycle characteristics, including good hemispheric coupling, phase relationship between the surface dipole and the BMR-generating internal field, and correlation between dipole strength at cycle maximum and peak amplitude of the next cycle. The saturation of the cycle amplitude takes place through the quenching of the BMR tilt as a function of the internal field. The observed statistical scatter about the mean BMR tilt, built into the model, acts as a source of stochasticity which dominates amplitude fluctuations. The model thus can produce Dalton-like epochs of strongly suppressed cycle amplitude lasting a few cycles and can even shut off entirely following an unfavorable sequence of emergence events.
△ Less
Submitted 12 January, 2017; v1 submitted 23 June, 2016;
originally announced June 2016.
-
Modelling turbulent stellar convection zones: sub-grid scales effects
Authors:
A. Strugarek,
P. Beaudoin,
A. S. Brun,
P. Charbonneau,
S. Mathis,
P. K. Smolarkiewicz
Abstract:
The impressive development of global numerical simulations of turbulent stellar interiors unveiled a variety of possible differential rotation (solar or anti-solar), meridional circulation (single or multi-cellular), and dynamo states (stable large scale toroidal field or periodically reversing magnetic fields). Various numerical schemes, based on the so-called anelastic set of equations, were use…
▽ More
The impressive development of global numerical simulations of turbulent stellar interiors unveiled a variety of possible differential rotation (solar or anti-solar), meridional circulation (single or multi-cellular), and dynamo states (stable large scale toroidal field or periodically reversing magnetic fields). Various numerical schemes, based on the so-called anelastic set of equations, were used to obtain these results. It appears today mandatory to assess their robustness with respect to the details of the numerics, and in particular to the treatment of turbulent sub-grid scales. We report on an ongoing comparison between two global models, the ASH and EULAG codes. In EULAG the sub-grid scales are treated implicitly by the numerical scheme, while in ASH their effect is generally modelled by using enhanced dissipation coefficients. We characterize the sub-grid scales effect in a turbulent convection simulation with EULAG. We assess their effect at each resolved scale with a detailed energy budget. We derive equivalent eddy-diffusion coefficients and use the derived diffusivities in twin ASH numerical simulations. We find a good agreement between the large-scale flows developing in the two codes in the hydrodynamic regime, which encourages further investigation in the magnetohydrodynamic regime for various dynamo solutions.
△ Less
Submitted 27 May, 2016;
originally announced May 2016.
-
Characterisation of the turbulent electromotive force and its magnetically-mediated quenching in a global EULAG-MHD simulation of solar convection
Authors:
C. Simard,
P. Charbonneau,
C. Dube
Abstract:
We perform a mean-field analysis of the EULAG-MHD millenium simulation of global magnetohydrodynamical convection presented in Passos et al. 2014. The turbulent electromotive force operating in the simulation is assumed to be linearly related to the cyclic axisymmetric mean magnetic field and its first spatial derivatives. At every grid point in the simulation's meridional plane, this assumed rela…
▽ More
We perform a mean-field analysis of the EULAG-MHD millenium simulation of global magnetohydrodynamical convection presented in Passos et al. 2014. The turbulent electromotive force operating in the simulation is assumed to be linearly related to the cyclic axisymmetric mean magnetic field and its first spatial derivatives. At every grid point in the simulation's meridional plane, this assumed relationship involves 27 independent tensorial coefficients. Expanding on Racine et al. 2011, we extract these coefficients from the simulation data through a least-squares minimization procedure based on singular value decomposition. The reconstructed alpha-tensor shows good agreement with that obtained by Racine et al. 2011, who did not include derivatives of the mean-field in their fit, as well as with the alpha-tensor extracted by Augustson et al. 2015 from a distinct ASH MHD simulation. The isotropic part of the turbulent magnetic diffusivity tensor beta is positive definite and reaches values of 5.0x10^7 m2s-1 in the middle of the convecting fluid layers. The spatial variations of both alpha_phiphi and beta_phiphi component are well reproduced by expressions obtained under the SOCA, with a good matching of amplitude requiring a turbulent correlation time about five times smaller than the estimated turnover time of the small-scale turbulent flow. We find the magnetic quenching of the alpha-effect to be driven primarily by a reduction of the small-scale flow's kinetic helicity, with variations of the current helicity playing a lesser role in most locations in the simulation domain. Our measurements of turbulent diffusivity quenching are restricted to the beta_phiphi component, but indicate a weaker quenching, by a factor of 1.36, than of the alpha effect, which in our simulation drops by a factor of three between the minimum and maximum phases of the magnetic cycle.
△ Less
Submitted 6 April, 2016;
originally announced April 2016.
-
A coupled 2$\times$2D Babcock-Leighton solar dynamo model. I. Surface magnetic flux evolution
Authors:
Alexandre Lemerle,
Paul Charbonneau,
Arnaud Carignan-Dugas
Abstract:
The need for reliable predictions of the solar activity cycle motivates the development of dynamo models incorporating a representation of surface processes sufficiently detailed to allow assimilation of magnetographic data. In this series of papers we present one such dynamo model, and document its behavior and properties. This first paper focuses on one of the model's key components, namely surf…
▽ More
The need for reliable predictions of the solar activity cycle motivates the development of dynamo models incorporating a representation of surface processes sufficiently detailed to allow assimilation of magnetographic data. In this series of papers we present one such dynamo model, and document its behavior and properties. This first paper focuses on one of the model's key components, namely surface magnetic flux evolution. Using a genetic algorithm, we obtain best-fit parameters of the transport model by least-squares minimization of the differences between the associated synthetic synoptic magnetogram and real magnetographic data for activity cycle 21. Our fitting procedure also returns Monte Carlo-like error estimates. We show that the range of acceptable surface meridional flow profiles is in good agreement with Doppler measurements, even though the latter are not used in the fitting process. Using a synthetic database of bipolar magnetic region (BMR) emergences reproducing the statistical properties of observed emergences, we also ascertain the sensitivity of global cycle properties, such as the strength of the dipole moment and timing of polarity reversal, to distinct realizations of BMR emergence, and on this basis argue that this stochasticity represents a primary source of uncertainty for predicting solar cycle characteristics.
△ Less
Submitted 26 November, 2015;
originally announced November 2015.
-
IAU commission 10 "Solar Activity": Legacy report and triennial report for 2012-2015
Authors:
Carolus J. Schrijver,
Lyndsay Fletcher,
Lidia van Driel-Gesztelyi,
Ayumi Asai,
Paul S. Cally,
Paul Charbonneau,
Sarah E. Gibson,
Daniel Gomez,
Siraj S. Hasan,
Astrid M. Veronig,
Yihua Yan
Abstract:
After more than half a century of community support related to the science of "solar activity'', IAU's Commission 10 was formally discontinued in 2015, to be succeeded by C.E2 with the same area of responsibility. On this occasion, we look back at the growth of the scientific disciplines involved around the world over almost a full century. Solar activity and fields of research looking into the re…
▽ More
After more than half a century of community support related to the science of "solar activity'', IAU's Commission 10 was formally discontinued in 2015, to be succeeded by C.E2 with the same area of responsibility. On this occasion, we look back at the growth of the scientific disciplines involved around the world over almost a full century. Solar activity and fields of research looking into the related physics of the heliosphere continue to be vibrant and growing, with currently over 2,000 refereed publications appearing per year from over 4,000 unique authors, publishing in dozens of distinct journals and meeting in dozens of workshops and conferences each year. The size of the rapidly growing community and of the observational and computational data volumes, along with the multitude of connections into other branches of astrophysics, pose significant challenges; aspects of these challenges are beginning to be addressed through, among others, the development of new systems of literature reviews, machine-searchable archives for data and publications, and virtual observatories. As customary in these reports, we highlight some of the research topics that have seen particular interest over the most recent triennium, specifically active-region magnetic fields, coronal thermal structure, coronal seismology, flares and eruptions, and the variability of solar activity on long time scales. We close with a collection of developments, discoveries, and surprises that illustrate the range and dynamics of the discipline.
△ Less
Submitted 12 October, 2015;
originally announced October 2015.
-
Evidence of active MHD instability in EULAG-MHD simulations of solar convection
Authors:
N. Lawson,
A. Strugarek,
P. Charbonneau
Abstract:
We investigate the possible development of magnetohydrodynamical instabilities in the EULAG-MHD "millenium simulation" of Passos & Charbonneau (2014). This simulation sustains a large-scale magnetic cycle characterized by solar-like polarity reversals taking place on a regular multidecadal cadence, and in which zonally-oriented bands of strong magnetic field accumulate below the convective layers,…
▽ More
We investigate the possible development of magnetohydrodynamical instabilities in the EULAG-MHD "millenium simulation" of Passos & Charbonneau (2014). This simulation sustains a large-scale magnetic cycle characterized by solar-like polarity reversals taking place on a regular multidecadal cadence, and in which zonally-oriented bands of strong magnetic field accumulate below the convective layers, in response to turbulent pumping from above in successive magnetic half-cycles. Key aspects of this simulation include low numerical dissipation and a strongly subadiabatic fluid layer underlying the convectively unstable layers corresponding to the modeled solar convection zone. These properties are conducive to the growth and development of two-dimensional instabilities otherwise suppressed by stronger dissipation. We find evidence for the action of a non-axisymmetric magnetoshear instability operating in the upper portions of the stably stratified fluid layers. We also investigate the possibility that the Tayler instability may be contributing to the destabilization of the large-scale axisymmetric magnetic component at high latitudes. On the basis of our analyses, we propose a global dynamo scenario whereby the magnetic cycle is driven primarily by turbulent dynamo action in the convecting layers, but MHD instabilities accelerate the dissipation of the magnetic field pumped down into the overshoot and stable layers, thus perhaps significantly influencing the magnetic cycle period. Support for this scenario is found in the distinct global dynamo behaviors observed in an otherwise identical EULAG-MHD simulations, using a different degree of subadiabaticity in the stable fluid layers underlying the convection zone.
△ Less
Submitted 24 September, 2015;
originally announced September 2015.
-
Meridional circulation dynamics from 3D MHD global simulations of solar convection
Authors:
Dario Passos,
Paul Charbonneau,
Mark Miesch
Abstract:
The form of the solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. Yet a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work we use results from 3D global simulations of solar convection to infer the dynamical behavior of the established meridiona…
▽ More
The form of the solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. Yet a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work we use results from 3D global simulations of solar convection to infer the dynamical behavior of the established meridional circulation. We make a direct comparison between the meridional circulation that arises in these simulations and the latest observations. Based on our results we argue that there should be an equatorward flow at the base of the convection zone at mid latitudes, below the current maximum depth helioseismic measures can probe (0.75 R). We also provide physical arguments to justify this behaviour. The simulations indicate that the meridional circulation undergoes substantial changes in morphology as the magnetic cycle unfolds. We close by discussing the importance of these dynamical changes for current methods of observation that involve long averaging periods of helioseismic data. Also noteworthy is the fact that these topological changes indicate a rich interaction between magnetic fields and plasma flows, which challenges the ubiquitous kinematic approach used in the vast majority of mean field dynamo simulations.
△ Less
Submitted 4 February, 2015;
originally announced February 2015.
-
Hemispheric Coupling: Comparing Dynamo Simulations and Observations
Authors:
Aimee A. Norton,
Paul Charbonneau,
Dario Passos
Abstract:
Numerical simulations that reproduce solar-like magnetic cycles can be used to generate long-term statistics. The variations in N-S hemispheric cycle synchronicity and amplitude produced in simulations has not been widely compared to observations. The observed limits on asymmetry show that hemispheric sunspot area production is no more than 20% asymmetric for cycles 12-23 and phase lags do not exc…
▽ More
Numerical simulations that reproduce solar-like magnetic cycles can be used to generate long-term statistics. The variations in N-S hemispheric cycle synchronicity and amplitude produced in simulations has not been widely compared to observations. The observed limits on asymmetry show that hemispheric sunspot area production is no more than 20% asymmetric for cycles 12-23 and phase lags do not exceed 20% (2 yrs) of the total cycle period. Independent studies have found a long-term trend in phase values as one hemisphere leads the other for ~four cycles. Such persistence in phase is not indicative of a stochastic phenomenon. We compare the findings to results from a numerical simulation of solar convection recently produced with the EULAG-MHD model. This simulation spans 1600 yrs and generated 40 regular, sunspot-like cycles. While the simulated cycle length is too long and the toroidal bands remain at too high of latitudes, some solar-like aspects of hemispheric asymmetry are reproduced. The model reproduces the synchrony of polarity inversions and onset of cycle as the simulated phase lags do not exceed 20% of the cycle period. Simulated amplitude variations between the N and S hemispheres are larger than observed in the Sun. The simulations show one hemisphere persistently leads the other for several successive cycles, placing an upper bound on the efficiency of transequatorial magnetic coupling mechanisms. These include magnetic diffusion, cross-equatorial mixing within elongated convective rolls and transequatorial meridional flow cells. One or more of these processes may lead to magnetic flux cancellation whereby the oppositely directed fields come in close proximity and cancel each other across the magnetic equator late in the solar cycle. We discuss the discrepancies between model and observations and the constraints they pose on possible mechanisms of hemispheric coupling.
△ Less
Submitted 25 November, 2014;
originally announced November 2014.
-
Predictive Capabilities of Avalanche Models for Solar Flares
Authors:
Antoine Strugarek,
Paul Charbonneau
Abstract:
We assess the predictive capabilities of various classes of avalanche models for solar flares. We demonstrate that avalanche models cannot generally be used to predict specific events due to their high sensitivity to their embedded stochastic process. We show that deterministically driven models can nevertheless alleviate this caveat and be efficiently used for large events predictions. Our result…
▽ More
We assess the predictive capabilities of various classes of avalanche models for solar flares. We demonstrate that avalanche models cannot generally be used to predict specific events due to their high sensitivity to their embedded stochastic process. We show that deterministically driven models can nevertheless alleviate this caveat and be efficiently used for large events predictions. Our results promote a new approach for large (typically X-class) solar flares predictions based on simple and computationally inexpensive avalanche models.
△ Less
Submitted 25 June, 2014;
originally announced June 2014.
-
Deterministically Driven Avalanche Models of Solar Flares
Authors:
Antoine Strugarek,
Paul Charbonneau,
Richard Joseph,
Dorian Pirot
Abstract:
We develop and discuss the properties of a new class of lattice-based avalanche models of solar flares. These models are readily amenable to a relatively unambiguous physical interpretation in terms of slow twisting of a coronal loop. They share similarities with other avalanche models, such as the classical stick--slip self-organized critical model of earthquakes, in that they are driven globally…
▽ More
We develop and discuss the properties of a new class of lattice-based avalanche models of solar flares. These models are readily amenable to a relatively unambiguous physical interpretation in terms of slow twisting of a coronal loop. They share similarities with other avalanche models, such as the classical stick--slip self-organized critical model of earthquakes, in that they are driven globally by a fully deterministic energy loading process. The model design leads to a systematic deficit of small scale avalanches. In some portions of model space, mid-size and large avalanching behavior is scale-free, being characterized by event size distributions that have the form of power-laws with index values, which, in some parameter regimes, compare favorably to those inferred from solar EUV and X-ray flare data. For models using conservative or near-conservative redistribution rules, a population of large, quasiperiodic avalanches can also appear. Although without direct counterparts in the observational global statistics of flare energy release, this latter behavior may be relevant to recurrent flaring in individual coronal loops. This class of models could provide a basis for the prediction of large solar flares.
△ Less
Submitted 19 February, 2014;
originally announced February 2014.
-
Torsional Oscillations in a Global Solar Dynamo
Authors:
P. Beaudoin,
P. Charbonneau,
E. Racine,
P. K. Smolarkiewicz
Abstract:
We characterize and analyze rotational torsional oscillations developing in a large-eddy magnetohydrodynamical simulation of solar convection (Ghizaru, Charbonneau, and Smolarkiewicz, Astrophys. J. Lett., 715, L133 (2010); Racine et al., Astrophys. J., 735, 46 (2011)) producing an axisymmetric large-scale magnetic field undergoing periodic polarity reversals. Motivated by the many solar-like featu…
▽ More
We characterize and analyze rotational torsional oscillations developing in a large-eddy magnetohydrodynamical simulation of solar convection (Ghizaru, Charbonneau, and Smolarkiewicz, Astrophys. J. Lett., 715, L133 (2010); Racine et al., Astrophys. J., 735, 46 (2011)) producing an axisymmetric large-scale magnetic field undergoing periodic polarity reversals. Motivated by the many solar-like features exhibited by these oscillations, we carry out an analysis of the large-scale zonal dynamics. We demonstrate that simulated torsional oscillations are not driven primarily by the periodically-varying large-scale magnetic torque, as one might have expected, but rather via the magnetic modulation of angular-momentum transport by the large-scale meridional flow. This result is confirmed by a straightforward energy analysis. We also detect a fairly sharp transition in rotational dynamics taking place as one moves from the base of the convecting layers to the base of the thin tachocline-like shear layer formed in the stably stratified fluid layers immediately below. We conclude by discussing the implications of our analyses with regards to the mechanism of amplitude saturation in the global dynamo operating in the simulation, and speculate on the possible precursor value of torsional oscillations for the forecast of solar cycle characteristics.
△ Less
Submitted 3 October, 2012;
originally announced October 2012.
-
Predicting Solar Flares by Data Assimilation in Avalanche Models. I. Model Design and Validation
Authors:
Eric Bélanger,
Alain Vincent,
Paul Charbonneau
Abstract:
Data assimilation techniques, developed in the last two decades mainly for weather prediction, produce better forecasts by taking advantage of both theoretical/numerical models and real-time observations. In this paper, we explore the possibility of applying the data-assimilation techniques known as 4D-VAR to the prediction of solar flares. We do so in the context of a continuous version of the…
▽ More
Data assimilation techniques, developed in the last two decades mainly for weather prediction, produce better forecasts by taking advantage of both theoretical/numerical models and real-time observations. In this paper, we explore the possibility of applying the data-assimilation techniques known as 4D-VAR to the prediction of solar flares. We do so in the context of a continuous version of the classical cellular-automaton-based self-organized critical avalanche models of solar flares introduced by Lu and Hamilton (Astrophys. J., 380, L89, 1991). Such models, although a priori far removed from the physics of magnetic reconnection and magneto-hydrodynamical evolution of coronal structures, nonetheless reproduce quite well the observed statistical distribution of flare characteristics. We report here on a large set of data assimilation runs on synthetic energy release time series. Our results indicate that, despite the unpredictable (and unobservable) stochastic nature of the driving/triggering mechanism within the avalanche model, 4D-VAR succeeds in producing optimal initial conditions that reproduce adequately the time series of energy released by avalanches/flares. This is an essential first step towards forecasting real flares.
△ Less
Submitted 14 August, 2007;
originally announced August 2007.
-
Stellar Structure Modeling using a Parallel Genetic Algorithm for Objective Global Optimization
Authors:
T. S. Metcalfe,
P. Charbonneau
Abstract:
Genetic algorithms are a class of heuristic search techniques that apply basic evolutionary operators in a computational setting. We have designed a fully parallel and distributed hardware/software implementation of the generalized optimization subroutine PIKAIA, which utilizes a genetic algorithm to provide an objective determination of the globally optimal parameters for a given model against…
▽ More
Genetic algorithms are a class of heuristic search techniques that apply basic evolutionary operators in a computational setting. We have designed a fully parallel and distributed hardware/software implementation of the generalized optimization subroutine PIKAIA, which utilizes a genetic algorithm to provide an objective determination of the globally optimal parameters for a given model against an observational data set. We have used this modeling tool in the context of white dwarf asteroseismology, i.e., the art and science of extracting physical and structural information about these stars from observations of their oscillation frequencies. The efficient, parallel exploration of parameter-space made possible by genetic-algorithm-based numerical optimization led us to a number of interesting physical results: (1) resolution of a hitherto puzzling discrepancy between stellar evolution models and prior asteroseismic inferences of the surface helium layer mass for a DBV white dwarf; (2) precise determination of the central oxygen mass fraction in a white dwarf star; and (3) a preliminary estimate of the astrophysically important but experimentally uncertain rate for the ^12C(alpha,gamma)^16O nuclear reaction. These successes suggest that a broad class of computationally-intensive modeling applications could also benefit from this approach.
△ Less
Submitted 22 November, 2002; v1 submitted 16 August, 2002;
originally announced August 2002.
-
Preliminary Constraints on 12C(alpha,gamma)16O from White Dwarf Seismology
Authors:
T. S. Metcalfe,
D. E. Winget,
P. Charbonneau
Abstract:
For many years, astronomers have promised that the study of pulsating white dwarfs would ultimately lead to useful information about the physics of matter under extreme conditions of temperature and pressure. In this paper we finally make good on that promise. Using observational data from the Whole Earth Telescope and a new analysis method employing a genetic algorithm, we empirically determine…
▽ More
For many years, astronomers have promised that the study of pulsating white dwarfs would ultimately lead to useful information about the physics of matter under extreme conditions of temperature and pressure. In this paper we finally make good on that promise. Using observational data from the Whole Earth Telescope and a new analysis method employing a genetic algorithm, we empirically determine that the central oxygen abundance in the helium-atmosphere variable white dwarf GD 358 is 84+/-3 percent. We use this value to place preliminary constraints on the 12C(alpha,gamma)16O nuclear reaction cross-section. More precise constraints will be possible with additional detailed simulations. We also show that the pulsation modes of our best-fit model probe down to the inner few percent of the stellar mass. We demonstrate the feasibility of reconstructing the internal chemical profiles of white dwarfs from asteroseismological data, and find an oxygen profile for GD 358 that is qualitatively similar to recent theoretical calculations.
△ Less
Submitted 11 April, 2001;
originally announced April 2001.