-
An impulsive geomagnetic effect from an early-impulsive flare
Authors:
Hugh S. Hudson,
Edward. W. Cliver,
Lyndsay Fletcher,
Declan A. Diver,
Peter T. Gallagher,
Ying Li,
Christopher M. J. Osborne,
Craig Stark,
Yang Su
Abstract:
The geomagnetic "solar flare effect" (SFE) results from excess ionization in the Earth's ionosphere, famously first detected at the time of the Carrington flare in 1859. This indirect detection of a flare constituted one of the first cases of "multimessenger astronomy," whereby solar ionizing radiation stimulates ionospheric currents. Well-observed SFEs have few-minute time scales and perturbation…
▽ More
The geomagnetic "solar flare effect" (SFE) results from excess ionization in the Earth's ionosphere, famously first detected at the time of the Carrington flare in 1859. This indirect detection of a flare constituted one of the first cases of "multimessenger astronomy," whereby solar ionizing radiation stimulates ionospheric currents. Well-observed SFEs have few-minute time scales and perturbations of >10 nT, with the greatest events reaching above 100 nT. In previously reported cases the SFE time profiles tend to resemble those of solar soft X-ray emission, which ionizes the D-region; there is also a less-well-studied contribution from Lyman-alpha. We report here a specific case, from flare SOL2024-03-10 (M7.4), in which an impulsive SFE deviated from this pattern. This flare contained an "early impulsive" component of exceptionally hard radiation, extending up to gamma-ray energies above 1 MeV, distinctly before the bulk of the flare soft X-ray emission. We can characterize the spectral distribution of this early-impulsive component in detail, thanks to the modern extensive wavelength coverage. A more typical gradual SFE occurred during the flare's main phase. We suggest that events of this type warrant exploration of the solar physics in the "impulse response" limit of very short time scales.
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
Doppler Dimming and Brightening Effects in Solar Prominences
Authors:
Aaron W. Peat,
Christopher M. J. Osborne,
Petr Heinzel
Abstract:
We explored the impact that Doppler dimming and brightening effects from bulk motions of solar prominences have on the formation of Lya, Ha, and MgII h line profiles. We compared two schemes in which these effects manifest; when the prominence is moving radially away from the solar surface (radial case), and when the prominence is moving parallel to the solar surface (horizontal case). To do this,…
▽ More
We explored the impact that Doppler dimming and brightening effects from bulk motions of solar prominences have on the formation of Lya, Ha, and MgII h line profiles. We compared two schemes in which these effects manifest; when the prominence is moving radially away from the solar surface (radial case), and when the prominence is moving parallel to the solar surface (horizontal case). To do this, we analysed 13,332 model profiles generated through the use of the 1D NLTE (i.e. departures from Local Thermodynamic equilibrium) radiative transfer (RT) code Promweaver, built on the Lightweaver NLTE RT framework to mimic the behaviour and output of the 1D NLTE RT code PROM. We found that horizontal velocities are just as, or more important than radial velocities. This demonstrates that horizontal velocities need to be accounted for when attempting to do any sort of forward modelling.
△ Less
Submitted 19 June, 2024; v1 submitted 18 June, 2024;
originally announced June 2024.
-
Assessment of SDSS-Derived Galaxy Morphologies Using HST Imaging
Authors:
Chandler Osborne,
Samir Salim
Abstract:
The Sloan Digital Sky Survey (SDSS) was foundational to the study of galaxy evolution, having revealed the bimodality of galaxies and the relationship between their structure and star-forming activity. However, ground-based optical surveys like SDSS are limited in resolution and depth which may lead to biases or poor quality in the derived morphological properties, potentially impacting our unders…
▽ More
The Sloan Digital Sky Survey (SDSS) was foundational to the study of galaxy evolution, having revealed the bimodality of galaxies and the relationship between their structure and star-forming activity. However, ground-based optical surveys like SDSS are limited in resolution and depth which may lead to biases or poor quality in the derived morphological properties, potentially impacting our understanding of how and why galaxies cease their star formation (quench). We use archival HST imaging of ~2,000 SDSS objects to assess the reliability of SDSS-derived morphologies, taking advantage of both SDSS statistical samples and of HST's superior resolution and sensitivity. Single Sersic fitting and bulge-disk decomposition is performed on HST images for direct comparison with SDSS results. Of the three catalogs of SDSS-derived morphologies considered, none are significantly more accurate than the others. For disk-dominated galaxies (n<2.5), global Sersic indices (n) from Meert et al. 2015 (M15) are preferred. For bulge-dominated galaxies (n>2.5), Simard et al. 2011 (S11) and M15 overestimate n by ~20%, and NYU-derived global n are preferred. Global R_eff from S11 are preferred, but overestimate R_eff for the largest galaxies by 0.1 dex. SDSS-derived single-component parameters are generally significantly more robust than SDSS-derived two-component parameters. The bulge Sersic index (n_bulge) cannot be reliably constrained from SDSS imaging. The bulge-to-total (B/T) ratio can be reliably inferred from SDSS for galaxies with SDSS B/T<0.6 provided that n_bulge=4 is enforced. The difference in global n between HST and SDSS depends strongly on B/T; an empirical correction based only on it accounts for most of the systematics in global n.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
The Bright Rim Prominences according to 2.5D Radiative Transfer
Authors:
Jack M. Jenkins,
Christopher M. J. Osborne,
Ye Qiu,
Rony Keppens,
Chuan Li
Abstract:
Solar prominences observed close to the limb commonly include a bright feature that, from the perspective of the observer, runs along the interface between itself and the underlying chromosphere. Despite several idealised models being proposed to explain the underlying physics, a more general approach remains outstanding. In this manuscript we demonstrate as a proof-of-concept the first steps in a…
▽ More
Solar prominences observed close to the limb commonly include a bright feature that, from the perspective of the observer, runs along the interface between itself and the underlying chromosphere. Despite several idealised models being proposed to explain the underlying physics, a more general approach remains outstanding. In this manuscript we demonstrate as a proof-of-concept the first steps in applying the Lightweaver radiative transfer framework's 2.5D extension to a `toy' model prominence + VAL3C chromosphere, inspired by recent 1.5D experiments that demonstrated a significant radiative chromosphere--prominence interaction. We find the radiative connection to be significant enough to enhance both the electron number density within the chromosphere, as well as its emergent intensity across a range of spectral lines in the vicinity of the filament absorption signature. Inclining the viewing angle from the vertical, we find these enhancements to become increasingly asymmetric and merge with a larger secondary enhancement sourced directly from the prominence underside. In wavelength, the enhancements are then found to be the largest in both magnitude and horizontal extent for the spectral line cores, decreasing into the line wings. Similar behaviour is found within new Chinese H$α$ Solar Explorer (CHASE)/H$α$ Imaging Spectrograph (HIS) observations, opening the door for subsequent statistical confirmations of the theoretical basis we develop here.
△ Less
Submitted 14 March, 2024;
originally announced March 2024.
-
Strategies for obtaining robust SED fitting parameters for galaxies at z~1 and z~2 in the absence of IR data
Authors:
Chandler Osborne,
Samir Salim
Abstract:
Robust estimation of star formation rates (SFRs) at higher redshifts (z>1) using UV-optical-NIR photometry is contingent on the ability of spectral energy distribution (SED) fitting to simultaneously constrain the dust attenuation, stellar metallicity, and star formation history (SFH). IR-derived dust luminosities can help break the degeneracy between these parameters, but IR data is often not ava…
▽ More
Robust estimation of star formation rates (SFRs) at higher redshifts (z>1) using UV-optical-NIR photometry is contingent on the ability of spectral energy distribution (SED) fitting to simultaneously constrain the dust attenuation, stellar metallicity, and star formation history (SFH). IR-derived dust luminosities can help break the degeneracy between these parameters, but IR data is often not available. Here, we explore strategies for SED fitting at z>1 in the absence of IR data using a sample of log M*>10.2 star-forming galaxies from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) for which 24mu data are available. We adopt the total IR luminosity (L_TIR) obtained from 24mu as the 'ground truth' that allows us to assess how well it can be recovered (as L_dust) from UV-optical-NIR SED fitting. We test a variety of dust attenuation models, stellar population synthesis models, metallicity assumptions, and SFHs separately to identify which assumptions maximize the agreement (correlation and linearity) between L_TIR and L_dust. We find that a flexible dust attenuation law performs best. For stellar populations, we find that BC03 models are favored over those of BPASS. Fixing the stellar metallicity at solar value is preferred to other fixed values or leaving it as a free parameter. For SFHs, we find that minimizing the variability in the recent (<100 Myr) SFH improves the agreement with L_TIR. Finally, we provide a catalog of galaxy parameters (including M* and SFR) for CANDELS galaxies with log M*>8 and 0.7<z<1.3 obtained using the models we found to be the most robust.
△ Less
Submitted 12 January, 2024;
originally announced January 2024.
-
$Σ_{\mathrm{SFR}}$-M* Diagram: A Valuable Galaxy Evolution Diagnostic to Complement (s)SFR-M* Diagrams
Authors:
Samir Salim,
Sandro Tacchella,
Chandler Osborne,
S. M. Faber,
Janice C. Lee,
Sara L. Ellison
Abstract:
The specific star formation rate (sSFR) is commonly used to describe the level of galaxy star formation (SF) and to select quenched galaxies. However, being a relative measure of the young-to-old population, an ambiguity in its interpretation may arise because a small sSFR can be either because of a substantial previous mass build up, or because SF is low. We show, using large samples spanning 0 <…
▽ More
The specific star formation rate (sSFR) is commonly used to describe the level of galaxy star formation (SF) and to select quenched galaxies. However, being a relative measure of the young-to-old population, an ambiguity in its interpretation may arise because a small sSFR can be either because of a substantial previous mass build up, or because SF is low. We show, using large samples spanning 0 < z < 2, that the normalization of SFR by the physical extent over which SF is taking place (i.e., SFR surface density, $Σ_{\mathrm{SFR}}$) overcomes this ambiguity. $Σ_{\mathrm{SFR}}$ has a strong physical basis, being tied to the molecular gas density and the effectiveness of stellar feedback, so we propose $Σ_{\mathrm{SFR}}$-M* as an important galaxy evolution diagram to complement (s)SFR-M* diagrams. Using the $Σ_{\mathrm{SFR}}$-M* diagram we confirm the Schiminovich et al. (2007) result that the level of SF along the main sequence today is only weakly mass dependent - high-mass galaxies, despite their redder colors, are as active as blue, low-mass ones. At higher redshift, the slope of the "$Σ_{\mathrm{SFR}}$ main sequence" steepens, signaling the epoch of bulge build-up in massive galaxies. We also find that $Σ_{\mathrm{SFR}}$ based on the optical isophotal radius more cleanly selects both the starbursting and the spheroid-dominated (early-type) galaxies than sSFR. One implication of our analysis is that the assessment of the inside-out vs. outside-in quenching scenarios should consider both sSFR and $Σ_{\mathrm{SFR}}$ radial profiles, because ample SF may be present in bulges with low sSFR (red color).
△ Less
Submitted 17 October, 2023;
originally announced October 2023.
-
Improved GALEX UV Photometry for 700,000 SDSS Galaxies
Authors:
Chandler Osborne,
Samir Salim,
Mederic Boquien,
Mark Dickinson,
Stephane Arnouts
Abstract:
The Galaxy Evolution Explorer (GALEX) satellite performed the first and only large-area UV survey, which in tandem with the Sloan Digital Sky Survey (SDSS) has facilitated modeling of the spectral energy distributions (SEDs) of low-redshift galaxies and the determination of various galaxy properties, in particular the star formation rate. However, the relatively crude angular resolution of GALEX (…
▽ More
The Galaxy Evolution Explorer (GALEX) satellite performed the first and only large-area UV survey, which in tandem with the Sloan Digital Sky Survey (SDSS) has facilitated modeling of the spectral energy distributions (SEDs) of low-redshift galaxies and the determination of various galaxy properties, in particular the star formation rate. However, the relatively crude angular resolution of GALEX (5") made its images susceptible to blending of sources, resulting in potentially biased far-UV (FUV) and near-UV (NUV) pipeline photometry. To remedy this issue and take advantage of model-fit photometry, we use the EMphot software to obtain forced GALEX photometry for ~700,000 SDSS galaxies at z<0.3. Positional priors of target galaxies and potentially contaminating neighbors were taken from SDSS. New photometry is based on the best-fitting of three model profiles: optical-like, exponential and flat. New photometry mitigates blending present in the original pipeline catalogs, which affected 16% of galaxies at a level of >0.2 mag and 2% at a level of >1 mag. Pipeline NUV magnitudes are severely affected (>1 mag) when the neighbor is brighter than the target galaxy and within 10", or when the neighbor is fainter and within ~3" of the target. New photometry fixes edge-of-detector bias, which affected pipeline photometry by up to 0.1 mag in NUV. We present catalogs with new photometry for GALEX observations of different depths, corresponding to the all-sky imaging survey (AIS), medium imaging survey (MIS) and deep imaging survey (DIS). Catalogs feature combined magnitudes for multiple detections of the same galaxy in a survey.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
1.5D NLTE spectral synthesis of a 3D filament/prominence simulation
Authors:
J. M. Jenkins,
C. M. J. Osborne,
R. Keppens
Abstract:
Aims. We here demonstrate how the recently developed Lightweaver framework makes non-LTE (NLTE) spectral synthesis feasible on a new 3D ab-initio magnetohydrodynamic (MHD) filament/prominence simulation, in a post-processing step. Methods. We clarify the need to introduce filament/prominent-specific Lightweaver boundary conditions that accurately model incident chromospheric radiation, and include…
▽ More
Aims. We here demonstrate how the recently developed Lightweaver framework makes non-LTE (NLTE) spectral synthesis feasible on a new 3D ab-initio magnetohydrodynamic (MHD) filament/prominence simulation, in a post-processing step. Methods. We clarify the need to introduce filament/prominent-specific Lightweaver boundary conditions that accurately model incident chromospheric radiation, and include a self-consistent and smoothly varying limb darkening function. Results. Progressing from isothermal/isobaric models to the self-consistently generated stratifications within a fully 3D MHD filament/prominence simulation, we find excellent agreement between our 1.5D non local thermodynamic equilibrium Lightweaver synthesis and a popular Hydrogen Hα proxy. We compute additional lines including Ca~\textsc{ii} 8542 alongside the more optically-thick Ca~\textsc{ii} H&K & Mg~\textsc{ii} h&k lines, for which no comparable proxy exists, and explore their formation properties within filament/prominence atmospheres. Conclusions. The versatility of the Lightweaver framework is demonstrated with this extension to 1.5D filament/prominence models, where each vertical column of the instantaneous 3D MHD state is spectrally analysed separately, without accounting for (important) multi-dimensional radiative effects. The general agreement found in the line core contrast of both observations and the Lightweaver-synthesised simulation further validates the current generation of solar filaments/prominences models constructed numerically with MPI-AMRVAC.
△ Less
Submitted 27 November, 2022;
originally announced November 2022.
-
Flare Kernels May be Smaller than You Think: Modelling the Radiative Response of Chromospheric Plasma Adjacent to a Solar Flare
Authors:
Christopher M. J. Osborne,
Lyndsay Fletcher
Abstract:
Numerical models of solar flares typically focus on the behaviour of directly-heated flare models, adopting magnetic field- aligned, plane-parallel methodologies. With high spatial- and spectral-resolution ground-based optical observations of flares, it is essential also to understand the response of the plasma surrounding these strongly heated volumes. We investigate the effects of the extreme ra…
▽ More
Numerical models of solar flares typically focus on the behaviour of directly-heated flare models, adopting magnetic field- aligned, plane-parallel methodologies. With high spatial- and spectral-resolution ground-based optical observations of flares, it is essential also to understand the response of the plasma surrounding these strongly heated volumes. We investigate the effects of the extreme radiation field produced by a heated column of flare plasma on an adjacent slab of chromospheric plasma, using a two-dimensional radiative transfer model and considering the time-dependent solution to the atomic level populations and electron density throughout this model. The outgoing spectra of H$α$ and Ca II 854.2 nm synthesised from our slab show significant spatial-, time-, and wavelength-dependent variations (both enhancements and reductions) in the line cores, extending on order 1 Mm into the non-flaring slab due to the incident transverse radiation field from the flaring boundary. This may lead to significant overestimates of the sizes of directly-heated flare kernels, if line-core observations are used. However, the radiation field alone is insufficient to drive any significant changes in continuum intensity, due to the typical photospheric depths at which they forms, so continuum sources will not have an apparent increase in size. We show that the line formation regions near the flaring boundary can be driven upwards in altitude by over 1 Mm despite the primary thermodynamic parameters (other than electron density) being held horizontally uniform. This work shows that in simple models these effects are significant and should be considered further in future flare modelling and interpretation.
△ Less
Submitted 7 September, 2022;
originally announced September 2022.
-
On the Importance of Ca II Photoionisation by the Hydrogen Lyman Transitions in Solar Flare Models
Authors:
Christopher M. J. Osborne,
Petr Heinzel,
Jana Kašparová,
Lyndsay Fletcher
Abstract:
The forward fitting of solar flare observations with radiation-hydrodynamic simulations is a common technique for learning about energy deposition and atmospheric evolution during these explosive events. A frequent spectral line choice for this process is Ca II 854.2 nm due to its formation in the chromosphere and substantial variability. It is important to ensure that this line is accurately mode…
▽ More
The forward fitting of solar flare observations with radiation-hydrodynamic simulations is a common technique for learning about energy deposition and atmospheric evolution during these explosive events. A frequent spectral line choice for this process is Ca II 854.2 nm due to its formation in the chromosphere and substantial variability. It is important to ensure that this line is accurately modeled to obtain the correct interpretation of observations. Here we investigate the importance of photoionisation of Ca II to Ca III by the hydrogen Lyman transitions; whilst the Lyman continuum is typically considered in this context in simulations, the associated bound-bound transitions are not. This investigation uses two RADYN flare simulations and reprocesses the radiative transfer using the Lightweaver framework which accounts for the overlapping of all active transitions. The Ca II 854.2 nm line profiles are found to vary significantly due to photoionisation by the Lyman lines, showing notably different shapes and even reversed asymmetries. Finally, we investigate to what extent these effects modify the energy balance of the simulation and the implications on future radiation-hydrodynamic simulations. There is found to be a 10-15% change in detailed optically thick radiative losses from considering these photoionisation effects on the calcium lines in the two simulations presented, demonstrating the importance of considering these effects in a self-consistent way.
△ Less
Submitted 23 July, 2021;
originally announced July 2021.
-
The Lightweaver Framework for NLTE Radiative Transfer in Python
Authors:
Christopher M J Osborne,
Ivan Milić
Abstract:
Tools for computing detailed optically thick spectral line profiles out of local thermodynamic equilibrium have always been focused on speed, due to the large computational effort involved. With the Lightweaver framework, we have produced a more flexible, modular toolkit for building custom tools in a high-level language, Python, without sacrificing speed against the current state of the art. The…
▽ More
Tools for computing detailed optically thick spectral line profiles out of local thermodynamic equilibrium have always been focused on speed, due to the large computational effort involved. With the Lightweaver framework, we have produced a more flexible, modular toolkit for building custom tools in a high-level language, Python, without sacrificing speed against the current state of the art. The goal of providing a more flexible method for constructing these complex simulations is to decrease the barrier to entry and allow more rapid exploration of the field.
In this paper we present an overview of the theory of optically thick NLTE radiative transfer, the numerical methods implemented in Lightweaver including the problems of time-dependent populations and charge-conservation, as well as an overview of the components most users will interact with, to demonstrate their flexibility.
△ Less
Submitted 30 June, 2021;
originally announced July 2021.
-
CANDELS Meets GSWLC: Evolution of the Relationship Between Morphology and Star Formation Since z = 2
Authors:
Chandler Osborne,
Samir Salim,
Ivana Damjanov,
S. M. Faber,
Marc Huertas-Company,
David C. Koo,
Kameswara Bharadwaj Mantha,
Daniel H. McIntosh,
Joel R. Primack,
Sandro Tacchella
Abstract:
Galaxy morphology and its evolution over the cosmic epoch hold important clues for understanding the regulation of star formation (SF). However, studying the relationship between morphology and SF has been hindered by the availability of consistent data at different redshifts. Our sample, combining CANDELS (0.8 < z < 2.5) and the GALEX-SDSS-WISE Legacy Catalog (GSWLC; z ~ 0), has physical paramete…
▽ More
Galaxy morphology and its evolution over the cosmic epoch hold important clues for understanding the regulation of star formation (SF). However, studying the relationship between morphology and SF has been hindered by the availability of consistent data at different redshifts. Our sample, combining CANDELS (0.8 < z < 2.5) and the GALEX-SDSS-WISE Legacy Catalog (GSWLC; z ~ 0), has physical parameters derived using consistent SED fitting with flexible dust attenuation laws. We adopt visual classifications from Kartaltepe et al. 2015 and expand them to z ~ 0 using SDSS images matching the physical resolution of CANDELS rest-frame optical images and deep FUV GALEX images matching the physical resolution of the CANDELS rest-frame FUV images. Our main finding is that disks with SF clumps at z ~ 0 make a similar fraction (~15%) of star-forming galaxies as at z ~ 2. The clumpy disk contribution to the SF budget peaks at z ~ 1, rather than z ~ 2, suggesting that the principal epoch of disk assembly continues to lower redshifts. Star-forming spheroids ("blue nuggets"), though less centrally concentrated than quenched spheroids, contribute significantly (~15%) to the SF budget at z ~ 1-2, suggesting that compaction precedes quenching. Among green valley and quiescent galaxies, the pure spheroid fraction drops since z ~ 1, whereas spheroids with disks (S0-like) become dominant. Mergers at or nearing coalescence are enhanced in SFR relative to the main sequence at all redshifts by a factor of ~2, but contribute $\lesssim$5% to the SF budget, with their contribution remaining small above the main sequence.
△ Less
Submitted 2 September, 2020;
originally announced September 2020.
-
Spectral deconvolution with deep learning: removing the effects of spectral PSF broadening
Authors:
Momchil Molnar,
Kevin Reardon,
Christopher Osborne,
Ivan Milić
Abstract:
We explore novel methods of recovering the original spectral line profiles from data obtained by instruments that sample those profiles with an extended or multipeaked spectral transmission profile. The techniques are tested on data obtained at high spatial resolution from the Fast Imaging Solar Spectrograph (FISS) grating spectrograph at the Big Bear Solar Observatory and from the Interferometric…
▽ More
We explore novel methods of recovering the original spectral line profiles from data obtained by instruments that sample those profiles with an extended or multipeaked spectral transmission profile. The techniques are tested on data obtained at high spatial resolution from the Fast Imaging Solar Spectrograph (FISS) grating spectrograph at the Big Bear Solar Observatory and from the Interferometric Bidimensional Spectrometer (IBIS) instrument at the Dunn Solar Telescope. The method robustly deconvolves wide spectral transmission profiles for fields of view sampling a variety of solar structures (granulation, plage and pores) with a photometrical precision of less than 1%. The results and fidelity of the method are tested on data from IBIS obtained using several different spectral resolution modes.
The method, based on convolutional neural networks (CNN), is extremely fast, performing about $10^5$ deconvolutions per second on a CPU and $10^6$ deconvolutions per second on NVIDIA TITAN RTX GPU for a spectrum with 40 wavelength samples. This approach is applicable for deconvolving large amounts of data from instruments with wide spectral transmission profiles, such as the Visible Tunable Filter (VTF) on the DKI Solar Telescope (DKIST). We also investigate its application to future instruments by recovering spectral line profiles obtained with a theoretical multi-peaked spectral transmission profile.
We further discuss the limitations of this deconvolutional approach through the analysis of the dimensionality of the original and multiplexed data.
△ Less
Submitted 11 May, 2020;
originally announced May 2020.
-
Thyr: A Volumetric Ray-Marching Tool for Simulating Microwave Emission
Authors:
Christopher M. J. Osborne,
Paulo J. A. Simões
Abstract:
Gyrosynchrotron radiation is produced by solar flares, and can be used to infer properties of the accelerated electrons and magnetic field of the flaring region. This microwave emission is highly dependent on many local plasma parameters, and the viewing angle. To correctly interpret observations, detailed simulations of the emission are required. Additionally, gyrosynchrotron emission from the ch…
▽ More
Gyrosynchrotron radiation is produced by solar flares, and can be used to infer properties of the accelerated electrons and magnetic field of the flaring region. This microwave emission is highly dependent on many local plasma parameters, and the viewing angle. To correctly interpret observations, detailed simulations of the emission are required. Additionally, gyrosynchrotron emission from the chromosphere has been largely ignored in modelling efforts, and recent studies have shown the importance of thermal emission at millimetric wavelengths. Thyr is a new tool for modelling microwave emission from three-dimensional flaring loops with spatially varying atmosphere and increased resolution in the lower corona and chromosphere. Thyr is modular and open-source, consisting of separate components to compute the thermal and non-thermal microwave emission coefficients and perform three-dimensional radiative transfer (in local thermodynamic equilibrium). The radiative transfer integral is computed by a novel ray-marching technique to efficiently compute the contribution of many volume elements. This technique can also be employed on a variety of astrophysics problems. Herein we present a review of the theory of gyrosynchrotron radiation, and two simulations of identical flare loops in low- and high-resolution performed with Thyr, with a spectral imaging analysis of differing regions. The high-resolution simulation presents a spectral hardening at higher frequencies. This hardening originates around the top of the chromosphere due to the strong convergence of the magnetic field, and is not present in previous models due to insufficient resolution. This hardening could be observed with a coordinated flare observation from active radio observatories.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
RADYNVERSION: Learning to Invert a Solar Flare Atmosphere with Invertible Neural Networks
Authors:
Christopher M. J. Osborne,
John A. Armstrong,
Lyndsay Fletcher
Abstract:
During a solar flare, it is believed that reconnection takes place in the corona followed by fast energy transport to the chromosphere. The resulting intense heating strongly disturbs the chromospheric structure, and induces complex radiation hydrodynamic effects. Interpreting the physics of the flaring solar atmosphere is one of the most challenging tasks in solar physics. Here we present a novel…
▽ More
During a solar flare, it is believed that reconnection takes place in the corona followed by fast energy transport to the chromosphere. The resulting intense heating strongly disturbs the chromospheric structure, and induces complex radiation hydrodynamic effects. Interpreting the physics of the flaring solar atmosphere is one of the most challenging tasks in solar physics. Here we present a novel deep learning approach, an invertible neural network, to understanding the chromospheric physics of a flaring solar atmosphere via the inversion of observed solar line profiles in Hα and Ca II λ8542. Our network is trained using flare simulations from the 1D radiation hydrodynamics code RADYN as the expected atmosphere and line profile. This model is then applied to single pixels from an observation of an M1.1 solar flare taken with SST/CRISP instrument just after the flare onset. The inverted atmospheres obtained from observations provide physical information on the electron number density, temperature and bulk velocity flow of the plasma throughout the solar atmosphere ranging from 0-10 Mm in height. The density and temperature profiles appear consistent with the expected atmospheric response, and the bulk plasma velocity provides the gradients needed to produce the broad spectral lines whilst also predicting the expected chromospheric evaporation from flare heating. We conclude that we have taught our novel algorithm the physics of a solar flare according to RADYN and that this can be confidently used for the analysis of flare data taken in these two wavelengths. This algorithm can also be adapted for a menagerie of inverse problems providing extremely fast ($\sim$10 μs) inversion samples.
△ Less
Submitted 29 April, 2019; v1 submitted 24 January, 2019;
originally announced January 2019.
