Search

In this paper, we introduce a new family of spatially co-located field solvers for particle-in-cell applications which evolve the potential formulation of Maxwell's equations under the Lorenz gauge. Our recent work introduced the concept of time-consistency, which connects charge conservation to the preservation of the gauge at the semi-discrete level. It will be shown that there exists a large family of time discretizations which satisfy this property. Additionally, it will be further shown that for large classes of time marching methods, the satisfaction of the gauge condition automatically implies the satisfaction of Gauss's law for electricity, with the potential formulation ensuring that that Gauss's law for magnetism is satisfied by definition. We focus on popular time marching methods including centered differences, backward differences, and diagonally-implicit Runge-Kutta methods, which are coupled to a spectral discretization in space. We demonstrate the theory by testing the methods on a relativistic Weibel instability and a drifting cloud of electrons.

Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still many open questions that cannot be answered with current facilities. In this report we summarise the science that is guiding the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST). We demonstrate how tranformational advances in topics including star formation in high redshift galaxies, the diffuse circumgalactic medium, Galactic ecology, cometary compositions and solar flares motivate the need for a 50m, single-dish telescope with a 1-2 degree field of view and a new generation of highly multiplexed continuum and spectral cameras. AtLAST will have the resolution to drastically lower the confusion limit compared to current single-dish facilities, whilst also being able to rapidly map large areas of the sky and detect extended, diffuse structures. Its high sensitivity and large field of view will open up the field of submillimeter transient science by increasing the probability of serendipitous detections. Finally, the science cases listed here motivate the need for a highly flexible operations model capable of short observations of individual targets, large surveys, monitoring programmes, target of opportunity observations and coordinated observations with other observatories. AtLAST aims to be a sustainable, upgradeable, multipurpose facility that will deliver orders of magnitude increases in sensitivity and mapping speeds over current and planned submillimeter observatories., Comment: 47 pages, 12 figures. For further details on AtLAST see https://atlast.uio.no

We report for the first time the detection of thermal free-free emission from post-flare loops at 34GHz in images from the Nobeyama Radioheliograph (NoRH). We studied 8 loops, 7 of which were from regions with extremely strong coronal magnetic field reported by Fedenev et al. (2023). Loop emission was observed in a wide range of wavelength bands, up to soft X-rays, confirming their multi-temperature structure and was associated with noise storm emission in metric wavelengths. The comparison of the 17GHz emission with that at 34GHz, after a calibration correction of the latter, showed that the emission was optically thin at both frequencies. We describe the structure and evolution of the loops and we computed their density, obtaining values for the top of the loops between 1 and 6 x 10^10 cm^-3, noticeably varying from one loop to another and in the course of the evolution of the same loop system; these values have only a weak dependence on the assumed temperature, 2 x 10^6 K in our case, as we are in the optically thin regime. Our density values are above those reported from EUV observations, which go up to about 10^10 cm^-3. This difference could be due to the fact that different emitting regions are sampled in the two domains and/or due to the more accurate diagnostics in the radio range, which do not suffer from inherent uncertainties arising from abundances and non-LTE excitation/ionization equilibria. We also estimated the magnetic field in the loop tops to be in the range of 10 to 30G., Comment: Accepted for publication in the Astrophysical Journal

This paper investigates the incidence of coherent emission in solar radio bursts, using a revised catalog of 3800 solar radio bursts observed by the Nobeyama Radio Polarimeters from 1988 to 2023. We focus on the 1.0 and 2.0 GHz data, where radio fluxes of order 10 billion Jansky have been observed. Previous work has suggested that these bursts are due to electron cyclotron maser (ECM) emission. In at least one well studied case, the bright emission at 1 GHz consists of narrowband spikes of millisecond duration. Coherent emission at 1 GHz can be distinguished from traditional incoherent gyrosynchrotron flare emission based on the radio spectrum: gyrosynchrotron emission at 1 GHz usually has a spectrum rising with frequency, so bursts in which 1 GHz is stronger than higher frequency measurements are unlikely to be incoherent gyrosynchrotron. Based on this criterion it is found that, for bursts exceeding 100 sfu, three-quarters of all bursts at 1 GHz and half of all 2 GHz bursts have a dominant coherent emission component, assumed to be ECM. The majority of the very bright bursts at 1 GHz are highly circularly polarized, consistent with a coherent emission mechanism, but not always 100% polarized. The frequency range from 1 to 2 GHz is heavily utilized for terrestrial applications, and these results are relevant for understanding the extreme flux levels that may impact such applications. Further, they provide a reference for comparison with the study of ECM emission from other stars and potentially exoplanets., Comment: 24 pages, 6 figures, in press

Space Weather is the study of the conditions in the solar wind that can affect life on the surface of the Earth, particularly the increasingly technologically sophisticated devices that are part of modern life. Solar radio observations are relevant to such phenomena because they generally originate as events in the solar atmosphere, including flares, coronal mass ejections and shocks, that produce electromagnetic and particle radiations that impact the Earth. Low frequency solar radio emission arises in the solar atmosphere at the levels where these events occur: we can use frequency as a direct measure of density, and an indirect measure of height, in the atmosphere. The main radio burst types are described and illustrated using data from the Green Bank Solar Radio Burst Spectrometer, and their potential use as diagnostics of Space Weather is discussed., Comment: 18 pages, 12 figures. Added to arxiv to provide appropriate reference

Observations at (sub-)millimeter wavelengths offer a complementary perspective on our Sun and other stars, offering significant insights into both the thermal and magnetic composition of their chromospheres. Despite the fundamental progress in (sub-)millimeter observations of the Sun, some important aspects require diagnostic capabilities that are not offered by existing observatories. In particular, simultaneous observations of the radiation continuum across an extended frequency range would facilitate the mapping of different layers and thus ultimately the 3D structure of the solar atmosphere. Mapping large regions on the Sun or even the whole solar disk at a very high temporal cadence would be crucial for systematically detecting and following the temporal evolution of flares, while synoptic observations, i.e., daily maps, over periods of years would provide an unprecedented view of the solar activity cycle in this wavelength regime. As our Sun is a fundamental reference for studying the atmospheres of active main sequence stars, observing the Sun and other stars with the same instrument would unlock the enormous diagnostic potential for understanding stellar activity and its impact on exoplanets. The Atacama Large Aperture Submillimeter Telescope (AtLAST), a single-dish telescope with 50\,m aperture proposed to be built in the Atacama desert in Chile, would be able to provide these observational capabilities. Equipped with a large number of detector elements for probing the radiation continuum across a wide frequency range, AtLAST would address a wide range of scientific topics including the thermal structure and heating of the solar chromosphere, flares and prominences, and the solar activity cycle. In this white paper, the key science cases and their technical requirements for AtLAST are discussed., Comment: 18 pages, 4 figures, submitted to Open Research Europe as part of a collection on the Atacama Large Aperture Submillimeter Telescope (AtLAST) -- revised version

In a previous paper, we developed a new particle-in-cell method for the Vlasov-Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell's equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in an non-staggered PIC formulation., Comment: 18 pages, 4 figures, 3 tables,1 algorithm. arXiv admin note: text overlap with arXiv:2208.11291

The solar soft X-ray observations from the GOES satellites provide one of the best quantitative records of solar activity, with nearly continuous flare records since 1975. We present a uniform analysis of the entire time series for 1975 to 2022 at NOAA class C1 level or above, to characterize the occurrence distribution function (ODF) of the flares observed in the 1-8 A spectral band. The analysis includes estimations of the peak fluxes of the 12 flares that saturated the 1-8 A time series. Our new estimates include NOAA's recently established correction factor (1.43) to adjust the GOES-1 through GOES-15 data covering 1975-2016. For each of the 12 saturated events we have made new estimates of peak fluxes based on fits to the rise and fall of the flare time profile, and have validated our extrapolation schemes by comparing with artificially truncated but unsaturated X10-class events. SOL2003-11-04 now has a peak flux of 4.32e-3 W/m^2. This corresponds to X43 on the new scale, or X30 on the old scale. We provide a list in the Appendix for peak fluxes of all 38 events above 10^-3 W/m^2, the GOES X10 level, including the 12 saturated events. The full list now gives us a first complete sample from which we obtain an occurrence distribution function (ODF) for peak energy flux , often represented as a power law dF/dE ~ E^-alpha, for which we find alpha = 1.973 +- 0.014 in the range M1 to X3. The power-law description fails at the high end, requiring a downward break in the ODF above the X10 level. We give a tapered powerlaw description of the resulting CCDF (complementary cumulative distribution function) and extrapolate it into the domain of "superflares," i.e. flares with bolometric energies > 10^33 erg. Extrapolation of this fit provides estimates of 100-yr and 1000-yr GOES peak fluxes that agree reasonably well with other such estimates using different data sets and methodology., Comment: Submitted to Solar Physics, September 2023

We identify a set of ~100 "cold" solar flares and perform a statistical analysis of them in the microwave range. Cold flares are characterized by a weak thermal response relative to nonthermal emission. This work is a follow up of a previous statistical study of cold flares, which focused on hard X-ray emission to quantify the flare nonthermal component. Here we focus on the microwave emission. The thermal response is represented by the soft X-ray emission measured by the GOES X-ray sensors. We obtain spectral parameters of the flare gyrosynchrotron emission and investigate patterns of the temporal evolution. The main results of the previous statistical study are confirmed: as compared to a "mean" flare, the cold flares have shorter durations, higher spectral peak frequencies, and harder spectral indices above the spectral peak. Nonetheless, there are some cold flares with moderate and low peak frequencies. In a majority of cold flares, we find evidence suggesting the presence of the Razin effect in the microwave spectra, indicative of rather dense flaring loops. We discuss the results in the context of electron acceleration efficiency.

FOXSI is a direct-imaging, hard X-ray (HXR) telescope optimized for solar flare observations. It detects hot plasma and energetic electrons in and near energy release sites in the solar corona via bremsstrahlung emission, measuring both spatial structure and particle energy distributions. It provides two orders of magnitude faster imaging spectroscopy than previously available, probing physically relevant timescales (<1s) never before accessible to address fundamental questions of energy release and efficient particle acceleration that have importance far beyond their solar application (e.g., planetary magnetospheres, flaring stars, accretion disks). FOXSI measures not only the bright chromospheric X-ray emission where electrons lose most of their energy, but also simultaneous emission from electrons as they are accelerated in the corona and propagate along magnetic field lines. FOXSI detects emission from high in the tenuous corona, where previous instruments have been blinded by nearby bright features and will fully characterizes the accelerated electrons and hottest plasmas as they evolve in energy, space, and time to solve the mystery of how impulsive energy release leads to solar eruptions, the primary drivers of space weather at Earth, and how those eruptions are energized and evolve., Comment: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 14 pages, 4 figures, 1 table

Understanding the nature of energetic particles in the solar atmosphere is one of the most important outstanding problems in heliophysics. Flare-accelerated particles compose a huge fraction of the flare energy budget; they have large influences on how events develop; they are an important source of high-energy particles found in the heliosphere; and they are the single most important corollary to other areas of high-energy astrophysics. Despite the importance of this area of study, this topic has in the past decade received only a small fraction of the resources necessary for a full investigation. For example, NASA has selected no new Explorer-class instrument in the past two decades that is capable of examining this topic. The advances that are currently being made in understanding flare-accelerated electrons are largely undertaken with data from EOVSA (NSF), STIX (ESA), and NuSTAR (NASA Astrophysics). This is despite the inclusion in the previous Heliophysics decadal survey of the FOXSI concept as part of the SEE2020 mission, and also despite NASA's having invested heavily in readying the technology for such an instrument via four flights of the FOXSI sounding rocket experiment. Due to that investment, the instrumentation stands ready to implement a hard X-ray mission to investigate flare-accelerated electrons. This white paper describes the scientific motivation for why this venture should be undertaken soon., Comment: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 15 pages, 5 figures

Solar flares and the often associated solar eruptive events serve as an outstanding laboratory to study the magnetic reconnection and the associated energy release and conversion processes under plasma conditions difficult to reproduce in the laboratory, and with considerable spatiotemporal details not possible elsewhere in the universe. In the past decade, thanks to advances in multi-wavelength imaging spectroscopy, as well as developments in theories and numerical modeling, significant progress has been made in improving our understanding of solar flare/eruption energy release. In particular, broadband imaging spectroscopy at microwave wavelengths offered by the Expanded Owens Valley Solar Array (EOVSA) has enabled the revolutionary capability of measuring the time-evolving coronal magnetic fields at or near the flare reconnection region. However, owing to EOVSA's limited dynamic range, imaging fidelity, and angular resolution, such measurements can only be done in a region around the brightest source(s) where the signal-to-noise is sufficiently large. In this white paper, after a brief introduction to the outstanding questions and challenges pertinent to magnetic energy release in solar flares and eruptions, we will demonstrate how a next-generation radio facility with many (~100-200) antenna elements can bring the next revolution by enabling high dynamic range, high fidelity broadband imaging spectropolarimetry along with a sub-second time resolution and arcsecond-level angular resolution. We recommend to prioritize the implementation of such a ground-based instrument within this decade. We also call for facilitating multi-wavelength, multi-messenger observations and advanced numerical modeling in order to achieve a comprehensive understanding of the "system science" of solar flares and eruptions., Comment: Science white paper submitted to the 2024 Solar and Space Physics Decadal Survey. All submitted white papers (including this one) are available at https://www.nationalacademies.org/our-work/decadal-survey-for-solar-and-space-physics-heliophysics-2024-2033

Coronal mass ejections (CMEs) are the most important drivers of space weather. Central to most CMEs is thought to be the eruption of a bundle of highly twisted magnetic field lines known as magnetic flux ropes. A comprehensive understanding of CMEs and their impacts hence requires detailed observations of physical parameters that lead to the formation, destabilization, and eventual eruption of the magnetic flux ropes. Recent advances in remote-sensing observations of coronal cavities, filament channels, sigmoids, EUV "hot channels," white light CMEs, and in situ observations of magnetic clouds points to the possibility of significant progress in understanding CMEs. In this white paper, we provide a brief overview of the potential of radio diagnostics for CMEs and CME progenitors, with a particular focus on the unique means for constraining their magnetic field and energetic electron population. Using synthetic observations based on realistic 3D MHD models, we also demonstrate the transformative potential of advancing such diagnostics by using broadband radio imaging spectropolarimetry with a high image dynamic range and high image fidelity. To achieve this goal, a solar-dedicated radio facility with such capabilities is recommended for implementation in the coming decade., Comment: Science white paper submitted to the 2024 Solar and Space Physics Decadal Survey. All submitted white papers (including this one) are available at https://www.nationalacademies.org/our-work/decadal-survey-for-solar-and-space-physics-heliophysics-2024-2033

The "middle corona," defined by West et al. (2022) as the region between ~1.5-6 solar radii, is a critical transition region that connects the highly structured lower corona to the outer corona where the magnetic field becomes predominantly radial. At radio wavelengths, remote-sensing of the middle corona falls in the meter-decameter wavelength range where a critical transition of radio emission mechanisms occurs. In addition, plasma properties of the middle corona can be probed by trans-coronal radio propagation methods including radio scintillation and Faraday rotation techniques. Together they offer a wealth of diagnostic tools for the middle corona, complementing current and planned missions at other wavelengths. These diagnostics include unique means for detecting and measuring the magnetic field and energetic electrons associated with coronal mass ejections, mapping coronal shocks and electron beam trajectories, as well as constraining the plasma density, magnetic field, and turbulence of the "young" solar wind. Following a brief overview of pertinent radio diagnostic methods, this white paper will discuss the current state of radio studies on the middle corona, challenges to obtaining a more comprehensive picture, and recommend an outlook in the next decade. Our specific recommendations for advancing the middle coronal sciences from the radio perspective are: (1) Prioritizing solar-dedicated radio facilities in the ~0.1-1 GHz range with broadband, high-dynamic-range imaging spectropolarimetry capabilities. (2) Developing facilities and techniques to perform multi-perspective, multiple lines-of-sight trans-coronal radio Faraday Rotation measurements., Comment: Science white paper submitted to the 2024 Solar and Space Physics Decadal Survey. All submitted white papers (including this one) are available at https://www.nationalacademies.org/our-work/decadal-survey-for-solar-and-space-physics-heliophysics-2024-2033. arXiv admin note: text overlap with arXiv:2208.04485

The sudden release of magnetic energy on the Sun drives powerful solar flares and coronal mass ejections. The key issue is the difficulty in predicting the occurrence time and location of strong solar eruptions, i.e., those leading to the high impact space weather disturbances at the near-Earth environment. Solar radio imaging helps identify the magnetic field characteristics of active regions susceptible to intense flares and energetic coronal mass ejections. Mapping of the Sun at X-band (8.1 -- 9.3 GHz) with the 12-m radio telescope at the Arecibo Observatory allows monitoring of the evolution of the brightness temperature of active regions in association with the development of magnetic complexity, which can lead to strong eruptions. For a better forecasting strategy in the future, such ground-based radio observations of high-spatial and temporal resolution, along with a full polarization capability, would have tremendous potential not only to understand the magnetic activity of solar eruptions, but also for revealing the particle acceleration mechanism and additional exciting science., Comment: 11 pages and 3 figures (White Paper Submitted to `Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033')

Dictionary learning, the problem of recovering a sparsely used matrix $\mathbf{D} \in \mathbb{R}^{M \times K}$ and $N$ $s$-sparse vectors $\mathbf{x}_i \in \mathbb{R}^{K}$ from samples of the form $\mathbf{y}_i = \mathbf{D}\mathbf{x}_i$, is of increasing importance to applications in signal processing and data science. When the dictionary is known, recovery of $\mathbf{x}_i$ is possible even for sparsity linear in dimension $M$, yet to date, the only algorithms which provably succeed in the linear sparsity regime are Riemannian trust-region methods, which are limited to orthogonal dictionaries, and methods based on the sum-of-squares hierarchy, which requires super-polynomial time in order to obtain an error which decays in $M$. In this work, we introduce SPORADIC (SPectral ORAcle DICtionary Learning), an efficient spectral method on family of reweighted covariance matrices. We prove that in high enough dimensions, SPORADIC can recover overcomplete ($K > M$) dictionaries satisfying the well-known restricted isometry property (RIP) even when sparsity is linear in dimension up to logarithmic factors. Moreover, these accuracy guarantees have an ``oracle property" that the support and signs of the unknown sparse vectors $\mathbf{x}_i$ can be recovered exactly with high probability, allowing for arbitrarily close estimation of $\mathbf{D}$ with enough samples in polynomial time. To the author's knowledge, SPORADIC is the first polynomial-time algorithm which provably enjoys such convergence guarantees for overcomplete RIP matrices in the near-linear sparsity regime.

The Frequency Agile Solar Radiotelescope (FASR) has been strongly endorsed as a top community priority by both Astronomy & Astrophysics Decadal Surveys and Solar & Space Physics Decadal Surveys in the past two decades. Although it was developed to a high state of readiness in previous years (it went through a CATE analysis and was declared ``doable now"), the NSF has not had the funding mechanisms in place to fund this mid-scale program. Now it does, and the community must seize this opportunity to modernize the FASR design and build the instrument in this decade. The concept and its science potential have been abundantly proven by the pathfinding Expanded Owens Valley Solar Array (EOVSA), which has demonstrated a small subset of FASR's key capabilities such as dynamically measuring the evolving magnetic field in eruptive flares, the temporal and spatial evolution of the electron energy distribution in flares, and the extensive coupling among dynamic components (flare, flux rope, current sheet). The FASR concept, which is orders of magnitude more powerful than EOVSA, is low-risk and extremely high reward, exploiting a fundamentally new research domain in solar and space weather physics. Utilizing dynamic broadband imaging spectropolarimetry at radio wavelengths, with its unique sensitivity to coronal magnetic fields and to both thermal plasma and nonthermal electrons from large flares to extremely weak transients, the ground-based FASR will make synoptic measurements of the coronal magnetic field and map emissions from the chromosphere to the middle corona in 3D. With its high spatial, spectral, and temporal resolution, as well as its superior imaging fidelity and dynamic range, FASR will be a highly complementary and synergistic component of solar and heliospheric capabilities needed for the next generation of solar science., Comment: White Paper submitted to the Solar and Space Physics 2024 Decadal Survey

We present the first survey of quiet Sun features observed in hard X-rays (HXRs), using the the Nuclear Spectroscopic Telescope ARray (NuSTAR), a HXR focusing optics telescope. The recent solar minimum combined with NuSTAR's high sensitivity has presented a unique opportunity to perform the first HXR imaging spectroscopy on a range of features in the quiet Sun. By studying the HXR emission of these features we can detect or constrain the presence of high temperature (>5 MK) or non-thermal sources, to help understand how they relate to larger more energetic solar phenomena, and determine their contribution to heating the solar atmosphere. We report on several features observed in the 28 September 2018 NuSTAR full-disk quiet Sun mosaics, the first of the NuSTAR quiet Sun observing campaigns, which mostly include steady features of X-ray bright points and an emerging flux region which later evolved into an active region, as well as a short-lived jet. We find that the features' HXR spectra are well fitted with isothermal models with temperatures ranging between 2.0-3.2 MK. Combining the NuSTAR data with softer X-ray emission from Hinode/XRT and EUV from SDO/AIA we recover the differential emission measures, confirming little significant emission above 4 MK. The NuSTAR HXR spectra allow us to constrain the possible non-thermal emission that would still be consistent with a null HXR detection. We found that for only one of the features (the jet) was there a potential non-thermal upper limit capable of powering the heating observed. However, even here the non-thermal electron distribution had to be very steep (effectively mono-energetic) with a low energy cut-off between 3-4 keV. The higher temperature or non-thermal sources in the typical quiet Sun features found in this September 2018 data are therefore found to be very weak, if present at all., Comment: Accepted for publication in Solar Physics

The Atacama Large Millimeter-submillimeter Array (ALMA), sited on the high desert plains of Chajnantor in Chile, has opened a new window onto solar physics in 2016 by providing continuum observations at millimeter and sub-millimeter wavelengths with an angular resolution comparable to that available at optical (O), ultraviolet (UV), extreme ultraviolet (EUV), and X-ray wavelengths, and with superior time resolution. In the intervening years, progress has been made testing and commissioning new observing modes and capabilities, in developing data calibration strategies, and in data imaging and restoration techniques. Here we review ALMA current solar observing capabilities, the process by which a user may propose to use the instrument, and summarize the observing process and work flow. We then discuss some of the challenges users may encounter in imaging and analyzing their data. We conclude with a discussion of additional solar observing capabilities and modes under consideration that are intended to further exploit the unique spectral coverage provided by ALMA.

This paper formulates a new particle-in-cell method for the Vlasov-Maxwell system. Under the Lorenz gauge condition, Maxwell's equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized momentum. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are shown to converge at the same rate as the fields in the both time and space. The field solver we consider in this work is first-order accurate in time and fifth-order accurate in space and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems, and the efficacy of the proposed formulation is demonstrated through a comparison with standard methods presented in the literature. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is close to the grid spacing. The use of high-order spatial approximations in the new method means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to standard explicit methods, which permits further computational savings., Comment: 47 pages, 24 figures, 6 tables, 1 algorithm

The physical conditions within our heliosphere are driven by the Sun's motion through an evolving interstellar environment that remains largely unexplored. The next generation of outer heliosphere and interstellar explorers will answer fundamental questions about the heliosphere's relationship with the very local interstellar medium (VLISM) by diving deeper into the Sun's interstellar surroundings. The impact of these future missions will be vastly enhanced by concurrent, interdisciplinary studies that examine the direct connections between conditions within the heliosphere, the heliosphere's immediate interstellar environment, and the larger-scale Galactic ISM. Comparisons of the heliosphere and VLISM to their analogs across the Galaxy will constrain the global processes shaping both stellar astrospheres and their sustained impact on the ISM., Comment: White paper submitted to the Heliophysics 2024 Decadal Survey

In recent decades our understanding of solar active regions (ARs) has improved substantially due to observations made with better angular resolution and wider spectral coverage. While prior AR observations have shown that these structures were always brighter than the quiet Sun at centimeter wavelengths, recent observations at millimeter and submillimeter wavelengths have shown ARs with well defined dark umbrae. Given this new information, it is now necessary to update our understanding and models of the solar atmosphere in active regions. In this work, we present a data-constrained model of the AR solar atmosphere, in which we use brightness temperature measurements of NOAA 12470 at three radio frequencies: 17 (NoRH), 100 and 230 GHz (ALMA). Based on our model, which assumes that the radio emission originates from thermal free-free and gyroresonance processes, we calculate radio brightness temperature maps that can be compared with the observations. The magnetic field at distinct atmospheric heights was determined in our modelling process by force-free field extrapolation using photospheric magnetograms taken by HMI/SDO. In order to determine the best plasma temperature and density height profiles necessary to match the observations, the model uses a genetic algorithm that modifies a standard quiet Sun atmospheric model. Our results show that the height of the transition region (TR) of the modelled atmosphere varies with the type of region being modelled: for umbrae the TR is located at 1080 +/- 20 km above the solar surface; for penumbrae, the TR is located at 1800 +/- 50 km; and for bright regions outside sunspots, the TR is located at 2000 +/- 100 km. With these results, we find good agreement with the observed AR brightness temperature maps. Our modelled AR can be used to estimate the emission at frequencies without observational coverage., Comment: 20 pages, 8 figures, 4 tables, accepted for publication in Frontiers in Astronomy and Space Sciences

Accurate and rapid detection of COVID-19 pneumonia is crucial for optimal patient treatment. Chest X-Ray (CXR) is the first line imaging test for COVID-19 pneumonia diagnosis as it is fast, cheap and easily accessible. Inspired by the success of deep learning (DL) in computer vision, many DL-models have been proposed to detect COVID-19 pneumonia using CXR images. Unfortunately, these deep classifiers lack the transparency in interpreting findings, which may limit their applications in clinical practice. The existing commonly used visual explanation methods are either too noisy or imprecise, with low resolution, and hence are unsuitable for diagnostic purposes. In this work, we propose a novel explainable deep learning framework (CXRNet) for accurate COVID-19 pneumonia detection with an enhanced pixel-level visual explanation from CXR images. The proposed framework is based on a new Encoder-Decoder-Encoder multitask architecture, allowing for both disease classification and visual explanation. The method has been evaluated on real world CXR datasets from both public and private data sources, including: healthy, bacterial pneumonia, viral pneumonia and COVID-19 pneumonia cases The experimental results demonstrate that the proposed method can achieve a satisfactory level of accuracy and provide fine-resolution classification activation maps for visual explanation in lung disease detection. The Average Accuracy, the Precision, Recall and F1-score of COVID-19 pneumonia reached 0.879, 0.985, 0.992 and 0.989, respectively. We have also found that using lung segmented (CXR) images can help improve the performance of the model. The proposed method can provide more detailed high resolution visual explanation for the classification decision, compared to current state-of-the-art visual explanation methods and has a great potential to be used in clinical practice for COVID-19 pneumonia diagnosis.