-
A spectroscopic investigation of thermal instability for cylindrical equilibria with background flow
Authors:
Joris Hermans,
Rony Keppens
Abstract:
Flows are omnipresent and govern the dynamics of plasma. Solar tornadoes are a class of apparently rotating prominences, that might be formed by thermal instability. In spectroscopic studies on thermal instability background flow is commonly neglected. We here determine the effect of background flow on thermal instability in cylindrical magnetic field configurations as the influence of various par…
▽ More
Flows are omnipresent and govern the dynamics of plasma. Solar tornadoes are a class of apparently rotating prominences, that might be formed by thermal instability. In spectroscopic studies on thermal instability background flow is commonly neglected. We here determine the effect of background flow on thermal instability in cylindrical magnetic field configurations as the influence of various parameters on the MHD spectrum. We investigate discrete thermal modes. In an analytical study, we extend upon the literature by including a generic background flow in a cylindrical coordinate system. The non-adiabatic MHD equations are linearised, Fourier-analysed, and are examined to understand how a background flow changes the continua. An approximate expression for discrete thermal modes is derived using a WKB analysis. The analytical results are then verified for a benchmark equilibrium using the eigenvalue code Legolas. The eigenfunctions of discrete thermal modes are visualised in 2D and 3D. The thermal continuum is Doppler-shifted due to the background flow, just like the slow and Alfvén continua. Discrete modes are altered because the governing equations contain flow-related terms. An approximate expression to predict the appearance of discrete thermal modes based on the equilibrium parameters is derived. All analytical expressions match the numerical results. The distribution of the density perturbations of the discrete thermal modes is not a singular condensation. 3D visualisation of the total velocity field shows that the helical field is heavily influenced by the radial velocity perturbation. We derived analytic expressions for non-adiabatic MHD modes of a cylindrical equilibrium with background flow and verified them using a coronal equilibrium. However, the equations are valid for and can be applied in other astrophysical environments.
△ Less
Submitted 29 February, 2024;
originally announced March 2024.
-
Effect of optically thin cooling curves on condensation formation: Case study using thermal instability
Authors:
Joris Hermans,
Rony Keppens
Abstract:
Non-gravitationally induced condensations are observed in many astrophysical environments. Such structures are formed due to energy loss by optically thin radiative emission. Instead of solving the full radiative transfer equations, precomputed cooling curves are typically used in numerical simulations. In the literature, there exists a wide variety of cooling curves and they are quite often used…
▽ More
Non-gravitationally induced condensations are observed in many astrophysical environments. Such structures are formed due to energy loss by optically thin radiative emission. Instead of solving the full radiative transfer equations, precomputed cooling curves are typically used in numerical simulations. In the literature, there exists a wide variety of cooling curves and they are quite often used as unquestionable ingredients. We determine the effect of the optically thin cooling curves on the formation and evolution of condensations. We perform a case study using thermal instability as a mechanism to form in situ condensations. We compare 2D numerical simulations with different cooling curves using interacting slow magnetohydrodynamic (MHD) waves as trigger for the thermal instability. Furthermore, we discuss a bootstrap measure to investigate the far non-linear regime of thermal instability. In the appendix, we include the details of all cooling curves implemented in MPI-AMRVAC and briefly discuss a hydrodynamic variant of the slow MHD waves setup for thermal instability. For all tested cooling curves, condensations are formed. However, the growth rates of the thermal instability are different. Also, the morphology of the formed condensation widely varies. We find fragmentation that is affected by the low-temperature treatment of the cooling curves. Condensations formed using cooling curves that vanish for temperatures lower than 20 000 K seem to be more stable against dynamical instabilities. The non-linear regime and fragmentation in the hydrodynamic case differ greatly from the MHD case. We advocate the use of modern cooling curves, based on accurate computations and up-to-date atomic parameters and solar abundances. Our bootstrap procedure can be used in future multi-dimensional simulations, to study fine-structure dynamics in solar prominences.
△ Less
Submitted 2 August, 2021; v1 submitted 15 July, 2021;
originally announced July 2021.
-
KIC5950759: a high-amplitude $δ$ Sct star with amplitude and frequency modulation near the terminal age main sequence
Authors:
D. M. Bowman,
J. Hermans,
J. Daszyńska-Daszkiewicz,
D. L. Holdsworth,
A. Tkachenko,
S. J. Murphy,
B. Smalley,
D. W. Kurtz
Abstract:
Amongst the intermediate mass pulsating stars known as $δ$ Sct stars is a subset of high-amplitude and predominantly radial-mode pulsators known as high-amplitude $δ$ Sct (HADS) stars. From more than 2000 $δ$ Sct stars observed by the Kepler space mission, only two HADS stars were detected. We investigate the more perplexing of these two HADS stars, KIC5950759. We study its variability using groun…
▽ More
Amongst the intermediate mass pulsating stars known as $δ$ Sct stars is a subset of high-amplitude and predominantly radial-mode pulsators known as high-amplitude $δ$ Sct (HADS) stars. From more than 2000 $δ$ Sct stars observed by the Kepler space mission, only two HADS stars were detected. We investigate the more perplexing of these two HADS stars, KIC5950759. We study its variability using ground- and space-based photometry, determine its atmospheric parameters from spectroscopy and perform asteroseismic modelling to constrain its mass and evolutionary stage. From spectroscopy, we find that KIC5950759 is a metal-poor star, which is in agreement with the inferred metallicity needed to reproduce its pulsation mode frequencies from non-adiabatic pulsation models. Furthermore, we combine ground-based WASP and Kepler space photometry, and measure a linear change in period of order $\dot{P}/P \simeq 10^{-6}$ yr$^{-1}$ for both the fundamental and first overtone radial modes across a time base of several years, which is at least two orders of magnitude larger than predicted by evolution models, and is the largest measured period change in a $δ$ Sct star to date. Our analysis indicates that KIC5950759 is a metal-poor HADS star near the short-lived contraction phase and the terminal-age main sequence, with its sub-solar metallicity making it a candidate SX Phe star. KIC5950759 is a unique object amongst the thousands of known $δ$ Sct stars and warrants further study to ascertain why its pulsation modes are evolving remarkably faster than predicted by stellar evolution.
△ Less
Submitted 19 April, 2021;
originally announced April 2021.
-
Towards constraining warm dark matter with stellar streams through neural simulation-based inference
Authors:
Joeri Hermans,
Nilanjan Banik,
Christoph Weniger,
Gianfranco Bertone,
Gilles Louppe
Abstract:
A statistical analysis of the observed perturbations in the density of stellar streams can in principle set stringent contraints on the mass function of dark matter subhaloes, which in turn can be used to constrain the mass of the dark matter particle. However, the likelihood of a stellar density with respect to the stream and subhaloes parameters involves solving an intractable inverse problem wh…
▽ More
A statistical analysis of the observed perturbations in the density of stellar streams can in principle set stringent contraints on the mass function of dark matter subhaloes, which in turn can be used to constrain the mass of the dark matter particle. However, the likelihood of a stellar density with respect to the stream and subhaloes parameters involves solving an intractable inverse problem which rests on the integration of all possible forward realisations implicitly defined by the simulation model. In order to infer the subhalo abundance, previous analyses have relied on Approximate Bayesian Computation (ABC) together with domain-motivated but handcrafted summary statistics. Here, we introduce a likelihood-free Bayesian inference pipeline based on Amortised Approximate Likelihood Ratios (AALR), which automatically learns a mapping between the data and the simulator parameters and obviates the need to handcraft a possibly insufficient summary statistic. We apply the method to the simplified case where stellar streams are only perturbed by dark matter subhaloes, thus neglecting baryonic substructures, and describe several diagnostics that demonstrate the effectiveness of the new method and the statistical quality of the learned estimator.
△ Less
Submitted 30 November, 2020;
originally announced November 2020.
-
Mining for Dark Matter Substructure: Inferring subhalo population properties from strong lenses with machine learning
Authors:
Johann Brehmer,
Siddharth Mishra-Sharma,
Joeri Hermans,
Gilles Louppe,
Kyle Cranmer
Abstract:
The subtle and unique imprint of dark matter substructure on extended arcs in strong lensing systems contains a wealth of information about the properties and distribution of dark matter on small scales and, consequently, about the underlying particle physics. However, teasing out this effect poses a significant challenge since the likelihood function for realistic simulations of population-level…
▽ More
The subtle and unique imprint of dark matter substructure on extended arcs in strong lensing systems contains a wealth of information about the properties and distribution of dark matter on small scales and, consequently, about the underlying particle physics. However, teasing out this effect poses a significant challenge since the likelihood function for realistic simulations of population-level parameters is intractable. We apply recently-developed simulation-based inference techniques to the problem of substructure inference in galaxy-galaxy strong lenses. By leveraging additional information extracted from the simulator, neural networks are efficiently trained to estimate likelihood ratios associated with population-level parameters characterizing substructure. Through proof-of-principle application to simulated data, we show that these methods can provide an efficient and principled way to simultaneously analyze an ensemble of strong lenses, and can be used to mine the large sample of lensing images deliverable by near-future surveys for signatures of dark matter substructure.
△ Less
Submitted 17 October, 2019; v1 submitted 4 September, 2019;
originally announced September 2019.