-
Parthenon -- a performance portable block-structured adaptive mesh refinement framework
Authors:
Philipp Grete,
Joshua C. Dolence,
Jonah M. Miller,
Joshua Brown,
Ben Ryan,
Andrew Gaspar,
Forrest Glines,
Sriram Swaminarayan,
Jonas Lippuner,
Clell J. Solomon,
Galen Shipman,
Christoph Junghans,
Daniel Holladay,
James M. Stone,
Luke F. Roberts
Abstract:
On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived…
▽ More
On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multi-dimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of $1.7\times10^{13}$ zone-cycles/s on 9,216 nodes (73,728 logical GPUs) at ~92% weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR.
△ Less
Submitted 21 November, 2022; v1 submitted 24 February, 2022;
originally announced February 2022.
-
Axion Dark Matter eXperiment: Run 1B Analysis Details
Authors:
ADMX Collaboration,
C. Bartram,
T. Braine,
R. Cervantes,
N. Crisosto,
N. Du,
G. Leum,
L. J Rosenberg,
G. Rybka,
J. Yang,
D. Bowring,
A. S. Chou,
R. Khatiwada,
A. Sonnenschein,
W. Wester,
G. Carosi,
N. Woollett,
L. D. Duffy,
M. Goryachev,
B. McAllister,
M. E. Tobar,
C. Boutan,
M. Jones,
B. H. Laroque,
N. S. Oblath
, et al. (23 additional authors not shown)
Abstract:
Searching for axion dark matter, the ADMX collaboration acquired data from January to October 2018, over the mass range 2.81--3.31 $μ$eV, corresponding to the frequency range 680--790 MHz. Using an axion haloscope consisting of a microwave cavity in a strong magnetic field, the ADMX experiment excluded Dine-Fischler-Srednicki-Zhitnisky (DFSZ) axions at 100% dark matter density over this entire fre…
▽ More
Searching for axion dark matter, the ADMX collaboration acquired data from January to October 2018, over the mass range 2.81--3.31 $μ$eV, corresponding to the frequency range 680--790 MHz. Using an axion haloscope consisting of a microwave cavity in a strong magnetic field, the ADMX experiment excluded Dine-Fischler-Srednicki-Zhitnisky (DFSZ) axions at 100% dark matter density over this entire frequency range, except for a few gaps due to mode crossings. This paper explains the full ADMX analysis for Run 1B, motivating analysis choices informed by details specific to this run.
△ Less
Submitted 13 October, 2020;
originally announced October 2020.
-
Axion Dark Matter eXperiment: Detailed Design and Operations
Authors:
R. Khatiwada,
D. Bowring,
A. S. Chou,
A. Sonnenschein,
W. Wester,
D. V. Mitchell,
T. Braine,
C. Bartram,
R. Cervantes,
N. Crisosto,
N. Du,
S. Kimes,
L. J Rosenberg,
G. Rybka,
J. Yang,
D. Will,
G. Carosi,
N. Woollett,
S. Durham,
L. D. Duffy,
R. Bradley,
C. Boutan,
M. Jones,
B. H. LaRoque,
N. S. Oblath
, et al. (26 additional authors not shown)
Abstract:
Axion Dark Matter eXperiment (ADMX) ultra low noise haloscope technology has enabled the successful completion of two science runs (1A and 1B) that looked for dark matter axions in the $2.66$ to $3.1$ $μ$eV mass range with Dine-Fischler-Srednicki-Zhitnisky (DFSZ) sensitivity Ref. [1,2]. Therefore, it is the most sensitive axion search experiment to date in this mass range. We discuss the technolog…
▽ More
Axion Dark Matter eXperiment (ADMX) ultra low noise haloscope technology has enabled the successful completion of two science runs (1A and 1B) that looked for dark matter axions in the $2.66$ to $3.1$ $μ$eV mass range with Dine-Fischler-Srednicki-Zhitnisky (DFSZ) sensitivity Ref. [1,2]. Therefore, it is the most sensitive axion search experiment to date in this mass range. We discuss the technological advances made in the last several years to achieve this sensitivity, which includes the implementation of components, such as state-of-the-art quantum limited amplifiers and a dilution refrigerator. Furthermore, we demonstrate the use of a frequency tunable Microstrip Superconducting Quantum Interference Device (SQUID) Amplifier (MSA), in Run 1A, and a Josephson Parametric Amplifier (JPA), in Run 1B, along with novel analysis tools that characterize the system noise temperature.
△ Less
Submitted 30 September, 2020;
originally announced October 2020.
-
On the occurrence of thermal non-equilibrium in coronal loops
Authors:
C. Froment,
F. Auchère,
Z. Mikić,
G. Aulanier,
K. Bocchialini,
E. Buchlin,
J. Solomon,
E. Soubrié
Abstract:
Long-period EUV pulsations, recently discovered to be common in active regions, are understood to be the coronal manifestation of thermal non-equilibrium (TNE). The active regions previously studied with EIT/SOHO and AIA/SDO indicated that long-period intensity pulsations are localized in only one or two loop bundles. The basic idea of this study is to understand why. For this purpose, we tested t…
▽ More
Long-period EUV pulsations, recently discovered to be common in active regions, are understood to be the coronal manifestation of thermal non-equilibrium (TNE). The active regions previously studied with EIT/SOHO and AIA/SDO indicated that long-period intensity pulsations are localized in only one or two loop bundles. The basic idea of this study is to understand why. For this purpose, we tested the response of different loop systems, using different magnetic configurations, to different stratifications and strengths of the heating. We present an extensive parameter-space study using 1D hydrodynamic simulations (1,020 in total) and conclude that the occurrence of TNE requires specific combinations of parameters. Our study shows that the TNE cycles are confined to specific ranges in parameter space. This naturally explains why only some loops undergo constant periodic pulsations over several days: since the loop geometry and the heating properties generally vary from one loop to another in an active region, only the ones in which these parameters are compatible exhibits TNE cycles. Furthermore, these parameters (heating and geometry) are likely to vary significantly over the duration of a cycle, which potentially limits the possibilities of periodic behavior. This study also confirms that long-period intensity pulsations and coronal rain are two aspects of the same phenomenon: both phenomena can occur for similar heating conditions and can appear simultaneously in the simulations.
△ Less
Submitted 12 February, 2018;
originally announced February 2018.
-
Long-Period Intensity Pulsations in Coronal Loops Explained by Thermal Non-Equilibrium Cycles
Authors:
Clara Froment,
Frédéric Auchère,
Guillaume Aulanier,
Zoran Mikić,
Karine Bocchialini,
Eric Buchlin,
Jacques Solomon
Abstract:
In solar coronal loops, thermal non-equilibrium (TNE) is a phenomenon that can occur when the heating is both highly-stratified and quasi-constant. Unambiguous observational identification of TNE would thus permit to strongly constrain heating scenarios. Up to now, while TNE is the standard interpretation of coronal rain, the long-term periodic evolution predicted by simulations has never been obs…
▽ More
In solar coronal loops, thermal non-equilibrium (TNE) is a phenomenon that can occur when the heating is both highly-stratified and quasi-constant. Unambiguous observational identification of TNE would thus permit to strongly constrain heating scenarios. Up to now, while TNE is the standard interpretation of coronal rain, the long-term periodic evolution predicted by simulations has never been observed yet. However, the detection of long-period intensity pulsations (periods of several hours) has been recently reported with SoHO/EIT, and this phenomenon appears to be very common in loops. Moreover, the three intensity-pulsation events that we recently studied with SDO/AIA show strong evidence for TNE in warm loops. In the present paper, a realistic loop geometry from LFFF extrapolations is used as input to 1D hydrodynamic simulations. Our simulations show that for the present loop geometry, the heating has to be asymmetrical to produce TNE. We analyse in detail one particular simulation that reproduces the average thermal behavior of one of the pulsating loop bundle observed with AIA. We compare the properties of this simulation with the properties deduced from the observations. The magnetic topology of the LFFF extrapolations points to the presence of sites of preferred reconnection at one footpoint, supporting the presence of asymmetric heating. In addition, we can reproduce the temporal large-scale intensity properties of the pulsating loops. This simulation further strengthens the interpretation of the observed pulsations as signatures of TNE. This thus gives important information on the heating localization and time scale for these loops.
△ Less
Submitted 5 January, 2017;
originally announced January 2017.
-
Thermal non-equilibriuem revealed by periodic pulses of random amplitudes in solar coronal loops
Authors:
F. Auchère,
C. Froment,
K. Bocchialini,
E. Buchlin,
J. Solomon
Abstract:
We recently detected variations in extreme ultraviolet intensity in coronal loops repeating with periods of several hours. Models of loops including stratified and quasi-steady heating predict the development of a state of thermal non-equilibrium (TNE): cycles of evaporative upflows at the footpoints followed by falling condensations at the apex. Based on Fourier and wavelet analysis, we demonstra…
▽ More
We recently detected variations in extreme ultraviolet intensity in coronal loops repeating with periods of several hours. Models of loops including stratified and quasi-steady heating predict the development of a state of thermal non-equilibrium (TNE): cycles of evaporative upflows at the footpoints followed by falling condensations at the apex. Based on Fourier and wavelet analysis, we demonstrate that the observed periodic signals are indeed not signatures of vibrational modes. Instead, superimposed on the power law expected from the stochastic background emission, the power spectra of the time series exhibit the discrete harmonics and continua expected from periodic trains of pulses of random amplitudes. These characteristics reinforce our earlier interpretation of these pulsations as being aborted TNE cycles.
△ Less
Submitted 2 February, 2018; v1 submitted 12 August, 2016;
originally announced August 2016.
-
On The Fourier And Wavelet Analysis Of Coronal Time Series
Authors:
F. Auchère,
C. Froment,
K. Bocchialini,
E. Buchlin,
J. Solomon
Abstract:
Using Fourier and wavelet analysis, we critically re-assess the significance of our detection of periodic pulsations in coronal loops. We show that the proper identification of the frequency dependence and statistical properties of the different components of the power spectra provies a strong argument against the common practice of data detrending, which tends to produce spurious detections aroun…
▽ More
Using Fourier and wavelet analysis, we critically re-assess the significance of our detection of periodic pulsations in coronal loops. We show that the proper identification of the frequency dependence and statistical properties of the different components of the power spectra provies a strong argument against the common practice of data detrending, which tends to produce spurious detections around the cut-off frequency of the filter. In addition, the white and red noise models built into the widely used wavelet code of Torrence & Compo cannot, in most cases, adequately represent the power spectra of coronal time series, thus also possibly causing false positives. Both effects suggest that several reports of periodic phenomena should be re-examined. The Torrence & Compo code nonetheless effectively computes rigorous confidence levels if provided with pertinent models of mean power spectra, and we describe the appropriate manner in which to call its core routines. We recall the meaning of the default confidence levels output from the code, and we propose new Monte-Carlo-derived levels that take into account the total number of degrees of freedom in the wavelet spectra. These improvements allow us to confirm that the power peaks that we detected have a very low probability of being caused by noise.
△ Less
Submitted 25 November, 2016; v1 submitted 16 June, 2016;
originally announced June 2016.
-
Energetic characterisation and statistics of solar coronal brightenings
Authors:
Vincent Joulin,
Éric Buchlin,
Jacques Solomon,
Chloé Guennou
Abstract:
To explain the high temperature of the corona, much attention has been paid to the distribution of energy in dissipation events. Indeed, if the event energy distribution is steep enough, the smallest, unobservable events could be the largest contributors to the total energy dissipation in the corona. Previous observations have shown a wide distribution of energies but remain inconclusive about the…
▽ More
To explain the high temperature of the corona, much attention has been paid to the distribution of energy in dissipation events. Indeed, if the event energy distribution is steep enough, the smallest, unobservable events could be the largest contributors to the total energy dissipation in the corona. Previous observations have shown a wide distribution of energies but remain inconclusive about the precise slope. Furthermore, these results rely on a very crude estimate of the energy. On the other hand, more detailed spectroscopic studies of structures such as coronal bright points do not provide enough statistical information to derive their total contribution to heating. We aim at getting a better estimate of the distributions of the energy dissipated in coronal heating events using high-resolution, multi-channel Extreme Ultra-Violet (EUV) data. To estimate the energies corresponding to heating events and deduce their distribution, we detect brightenings in five EUV channels of the Atmospheric Imaging Assembly (AIA) on-board the Solar Dynamics Observatory (SDO). We combine the results of these detections and we use maps of temperature and emission measure derived from the same observations to compute the energies. We obtain distributions of areas, durations, intensities, and energies (thermal, radiative, and conductive) of events. These distributions are power-laws, and we find also power-law correlations between event parameters. The energy distributions indicate that the energy from a population of events like the ones we detect represents a small contribution to the total coronal heating, even when extrapolating to smaller scales. The main explanations for this are how heating events can be extracted from observational data, and the incomplete knowledge of the thermal structure and processes in the coronal plasma attainable from available observations.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
Evidence for evaporation-incomplete condensation cycles in warm solar coronal loops
Authors:
Clara Froment,
Frédéric Auchère,
Karine Bocchialini,
Eric Buchlin,
Chloé Guennou,
Jacques Solomon
Abstract:
Quasi-constant heating at the footpoints of loops leads to evaporation and condensation cycles of the plasma: thermal non-equilibrium (TNE). This phenomenon is believed to play a role in the formation of prominences and coronal rain. However, it is often discarded to be involved in the heating of warm loops as the models do not reproduce observations. Recent simulations have shown that these incon…
▽ More
Quasi-constant heating at the footpoints of loops leads to evaporation and condensation cycles of the plasma: thermal non-equilibrium (TNE). This phenomenon is believed to play a role in the formation of prominences and coronal rain. However, it is often discarded to be involved in the heating of warm loops as the models do not reproduce observations. Recent simulations have shown that these inconsistencies with observations may be due to oversimplifications of the geometries of the models. In addition, our recent observations reveal that long-period intensity pulsations (several hours) are common in solar coronal loops. These periods are consistent with those expected from TNE. The aim of this paper is to derive characteristic physical properties of the plasma for some of these events to test the potential role of TNE in loop heating. We analyzed three events in detail using the six EUV coronal channels of SDO/AIA. We performed both a Differential Emission Measure (DEM) and a time-lag analysis, including a new method to isolate the relevant signal from the foreground and background emission. For the three events, the DEM undergoes long-period pulsations, which is a signature of periodic heating even though the loops are captured in their cooling phase, as is the bulk of the active regions. We link long-period intensity pulsations to new signatures of loop heating with strong evidence for evaporation and condensation cycles. We thus witness simultaneously widespread cooling and TNE. Finally, we discuss the implications of our new observations for both static and impulsive heating models.
△ Less
Submitted 26 August, 2015; v1 submitted 30 April, 2015;
originally announced April 2015.
-
Long-period intensity pulsations in the solar corona during activity cycle 23
Authors:
F. Auchère,
K. Bocchialini,
J. Solomon,
E. Tison
Abstract:
We report on the detection (10 σ) of 917 events of long-period (3 to 16 hours) intensity pulsations in the 19.5 nm passband of the SOHO Extreme ultraviolet Imaging Telescope. The data set spans from January 1997 to July 2010, i.e the entire solar cycle 23 and the beginning of cycle 24. The events can last for up to six days and have relative amplitudes up to 100%. About half of the events (54%) ar…
▽ More
We report on the detection (10 σ) of 917 events of long-period (3 to 16 hours) intensity pulsations in the 19.5 nm passband of the SOHO Extreme ultraviolet Imaging Telescope. The data set spans from January 1997 to July 2010, i.e the entire solar cycle 23 and the beginning of cycle 24. The events can last for up to six days and have relative amplitudes up to 100%. About half of the events (54%) are found to happen in active regions, and 50% of these have been visually associated with coronal loops. The remaining 46% are localized in the quiet Sun. We performed a comprehensive analysis of the possible instrumental artifacts and we conclude that the observed signal is of solar origin. We discuss several scenarios which could explain the main characteristics of the active region events. The long periods and the amplitudes observed rule out any explanation in terms of magnetohydrodynamic waves. Thermal nonequilibrium could produce the right periods, but it fails to explain all the observed properties of coronal loops and the spatial coherence of the events. We propose that moderate temporal variations of the heating term in the energy equation, so as to avoid a thermal nonequilibrium state, could be sufficient to explain those long-period intensity pulsations. The large number of detections suggests that these pulsations are common in active regions. This would imply that the measurement of their properties could provide new constraints on the heating mechanisms of coronal loops.
△ Less
Submitted 17 December, 2013; v1 submitted 13 December, 2013;
originally announced December 2013.
-
Oscillatory motions observed in eruptive filaments
Authors:
K. Bocchialini,
F. Baudin,
S. Koutchmy,
G. Pouget,
J. Solomon
Abstract:
Context: The origin of the variable component of the solar wind is of great intrinsic interest for heliophysics and space-weather, e.g. the initiation of coronal mass ejections, and the problem of mass loss of all stars. It is also related to the physics of coronal neutral sheets and streamers, occurring above lines of magnetic polarity reversal. Filaments and prominences correspond to the cool co…
▽ More
Context: The origin of the variable component of the solar wind is of great intrinsic interest for heliophysics and space-weather, e.g. the initiation of coronal mass ejections, and the problem of mass loss of all stars. It is also related to the physics of coronal neutral sheets and streamers, occurring above lines of magnetic polarity reversal. Filaments and prominences correspond to the cool coronal component of these regions.
Aims: We examine the dynamical behaviour of these structures where reconnection and dissipation of magnetic energy in the turbulent plasma are occurring. The link between the observed oscillatory motions and the eruption occurrence is investigated in detail for two different events.
Method: Two filaments are analysed using two different datasets: time series of spectra using a transition region line (He I at 584.33 A) and a coronal line (Mg X at 609.79 A) measured with CDS on-board SOHO, observed on May 30, 2003, and time series of intensity and velocity images from the NSO/Dunn Solar Telescope in the Halpha line on September 18, 1994 for the other. The oscillatory content is investigated using Fourier transform and wavelet analysis and is compared to different models.
Results: In both filaments, oscillations are clearly observed, in intensity and velocity in the He I and Mg X lines, in velocity in Halpha, with similar periods from a few minutes up to 80 minutes, with a main range from 20 to 30 minutes, cotemporal with eruptions. Both filaments exhibit vertical oscillating motions. For the filament observed in the UV (He I and Mg X lines), we provide evidence of damped velocity oscillations, and for the filament observed in the visible (Halpha line), we provide evidence that parts of the filament are oscillating, while the filament is moving over the solar surface, before its disappearance.
△ Less
Submitted 4 July, 2011;
originally announced July 2011.
-
Solar off-limb line widths: Alfven waves, ion-cyclotron waves, and preferential heating
Authors:
L. Dolla,
J. Solomon
Abstract:
Alfven waves and ion-cyclotron absorption of high-frequency waves are frequently brought into models devoted to coronal heating and fast solar-wind acceleration. Signatures of ion-cyclotron resonance have already been observed in situ in the solar wind (HELIOS spacecrafts) and, recently, in the upper corona (UVCS/SOHO remote-sensing results). We propose a method to constrain both the Alfven wave…
▽ More
Alfven waves and ion-cyclotron absorption of high-frequency waves are frequently brought into models devoted to coronal heating and fast solar-wind acceleration. Signatures of ion-cyclotron resonance have already been observed in situ in the solar wind (HELIOS spacecrafts) and, recently, in the upper corona (UVCS/SOHO remote-sensing results). We propose a method to constrain both the Alfven wave amplitude and the preferential heating induced by ion-cyclotron resonance, above a partially developed polar coronal hole observed with the SUMER/SOHO spectrometer. The instrumental stray light contribution is first substracted from the spectra. By supposing that the non-thermal velocity is related to the Alfven wave amplitude, it is constrained through a density diagnostic and the gradient of the width of the Mg X 625 A line. The temperatures of several coronal ions, as functions of the distance above the limb, are then determined by substracting the non-thermal component to the observed line widths. The effect of stray light explains the apparent decrease with height in the width of several spectral lines, this decrease usually starting about 0.1-0.2 Rs above the limb. This result rules out any direct evidence of damping of the Alfven waves, often suggested by other authors. We also find that the ions with the smallest charge-to-mass ratios are the hottest ones at a fixed altitude and that they are subject to a stronger heating, as compared to the others, between 57" and 102" above the limb. This constitutes a serious clue to ion-cyclotron preferential heating.
△ Less
Submitted 18 April, 2008;
originally announced April 2008.