Abstract
We present the first survey of quiet Sun features observed in hard X-rays (HXRs), using 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.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The study of the hard X-ray (HXR) emission from the quiet Sun could provide insight into the source of the sustained high temperature of the solar corona, termed the coronal heating problem. It was suggested by Parker (1988) that the source of this heating could be a large number of small-scale energy release events, taking place all through the solar cycle. If such events were weaker versions of flares, they would be expected to produce a HXR signature.
One phenomenon that has been linked to coronal heating is coronal bright points (CBPs), small-scale loop structures located in the lower corona, which are observed in extreme ultraviolet (EUV) and soft X-rays (SXRs) (e.g. Madjarska, 2019). These features are observed throughout the solar cycle, including solar minimum, when the Sun is quiet in the absence of active regions and large flares. Bright points have been studied extensively in EUV and SXRs. In EUV, these features typically have lifetimes of < 20 hours (Alipour and Safari, 2015; Zhang, Kundu, and White, 2001), whereas they have been found to be shorter-lived in SXRs, with lifetimes of ∼ 12 hours (Harvey et al., 1993). These features sometimes have associated transient phenomena, such as “microflares” (Golub et al., 1974; Shimojo and Shibata, 1999), and small-scale eruptions (Mou et al., 2018).
Previous studies have used EUV and SXR observations to investigate the temperatures of bright points. Kariyappa et al. (2011) used SXR data from the Hinode X-ray Telescope (Hinode/XRT) (Kosugi et al., 2007) to study a number of bright points, using filter ratios to determine that their temperatures ranged between 1.1 – 3.4 MK. Another study by Doschek et al. (2010) used EUV data from the Hinode Extreme Ultra-violet Imaging Spectrometer (Hinode/EIS) finding that the bright points reached maximum temperatures between 2 – 3 MK. Alexander, Del Zanna, and Maclean (2011) studied the evolution of a single bright point, using both SXR and EUV data. This study found that throughout the 13 hours of observation, the bright point was almost isothermal, with an average temperature of 1.3 MK. The Hinode/XRT time profile for this bright point showed a steady increase in intensity for the first two hours of observation, followed by several spikes (which the authors speculated were likely due to heating or reconnection events) until it began to decay.
Bright points can result from bipolar flux emergence. When this process occurs, new magnetic flux emerges to create an emerging flux region (EFR), which may subsequently evolve into a small-scale bright point, an example of which was investigated by Kontogiannis et al. (2020). This study tracked the evolution in Hinode/XRT of a bright point associated with an EFR, finding that the bright point exhibited a continuous increase in emission for ∼ 1.5 hours before it began to fade. However, rather than a bright point, an EFR may instead evolve into a large-scale active region if the emergence continues (van Driel-Gesztelyi and Green, 2015).
HXR observations of bright points have been made difficult by the faintness of this emission from these sources, and the lack of a solar dedicated instrument able to observe them individually. Previous HXR studies of the quiet Sun (Hannah et al., 2007, 2010) have been performed using data from the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) (Lin et al., 2002). However, as RHESSI was designed for the observation of bright sources, such as large flare events, only upper limits on the HXR emission from the whole solar disk were obtained. The authors noted that using a sensitive HXR focusing telescope would allow a more detailed study of the energy release mechanisms occurring in quiet Sun features. A more recent paper by Buitrago-Casas et al. (2022) also constrains the quiet Sun HXR emission in the 5 – 10 keV range using data from the Focusing Optics X-ray Solar Imager (FOXSI) sounding rocket. This study found similar upper limits to Hannah et al. (2010), though using only minutes’ worth of data, compared to the days that were required with RHESSI.
The Nuclear Spectroscopic Telescope ARray (NuSTAR) (Harrison et al., 2013) is a HXR focusing telescope which is capable of being pointed at the Sun to provide sensitive observations of faint solar sources (Grefenstette et al., 2016). Since 2014, there have been a number of NuSTAR solar observing campaigns.Footnote 1 Much of the work on the NuSTAR solar observations has focused on active region microflares (Wright et al., 2017; Glesener et al., 2017; Hannah et al., 2019; Cooper et al., 2020; Duncan et al., 2021; Cooper et al., 2021), with Glesener et al. (2020) reporting for the first time on non-thermal emission detected in a microflare observed with NuSTAR. Kuhar et al. (2018) presented work on quiet Sun flares observed by NuSTAR outside of an active region, finding that their temperatures ranged from 3.2 – 4.1 MK.
NuSTAR’s use of focusing optics means that it can directly image very faint HXR sources on the quiet Sun, and perform spectroscopy on regions of interest. During the recent solar minimum between Cycles 24 and 25 (2018 – 2020), when the solar disk was free of active regions, NuSTAR was used to observe the Sun on a number of occasions, providing several bright points and other quiet Sun phenomena to study. The recent solar minimum combined with NuSTAR’s sensitivity has provided a unique opportunity to study the HXR emission from these features and investigate their contribution to the heating of the solar atmosphere by searching for the presence of a high-temperature (> 5 MK) or non-thermal component due to the presence of accelerated electrons.
Here, we present the first survey of small features in the quiet Sun observed in HXRs with NuSTAR. We present analysis of several features from the 28 September 2018 full-disk mosaics, the first of the NuSTAR quiet Sun observations, including the first HXR imaging spectroscopy of such features. We include in our analysis EUV data from the Solar Dynamics Observatory’s Atmospheric Imaging Assembly (SDO/AIA) (Lemen et al., 2012) and SXR data from Hinode/XRT. An overview of this observation is presented in Section 2. The methods used to analyse the quiet Sun features are detailed in Section 3. The detailed analysis of an EFR, X-ray bright points, a jet and bright limb source are presented in Sections 4, 5, 6 and 7, respectively. A comparison of the thermal properties of these features is discussed in Section 8.
2 Overview of Observation
On 28 September 2018, NuSTAR observed the quiet Sun for two orbits (18:25 – 19:25 UT and 20:01 – 21:01 UT), producing a full-disk mosaic for each. As each of NuSTAR’s two focal plane modules (FPMA and FPMB) cover approximately the same \(12'\) × \(12'\) field of view (FOV; Harrison et al., 2013), multiple pointings are required to build up an image of the full-disk. These are formed over an orbit (∼ an hour), comprised of 25 pointings each of duration ∼ 100 s, in a 5 × 5 grid pattern. The first pointing, P1, is in the top left corner, with the subsequent pointings shifting to the right until P5. P6 is shifted down from P5, with the following four pointings being increasingly shifted to the left. This pattern continues until P25, which is in the bottom right corner of the mosaic. The pointings overlap with each other such that a feature on the disk can be captured up to four times within a single orbit.
For the times of these two orbits, there is also full-disk data available from SDO/AIA and Hinode/XRT (though there is a gap in the data from Hinode/XRT between 19:13 – 19:38 UT). The NuSTAR mosaics, and the SDO/AIA 211 Å and Hinode/XRT Be-thin full-disk images from the mid-times of both NuSTAR orbits, are shown in Figure 1. NuSTAR is a photon-counting detector, and the data list the photons’ properties, including time of detection, energy, and position on the detector. These NuSTAR full-disk maps are constructed by correcting each pointing for livetime (the fraction of time that the detectors are recording incoming photons) individually and aligning the NuSTAR image with SDO/AIA, and then summing the corrected counts.
Because of the way in which the mosaic tiles overlap, the disk is sampled more times than the limb (which is captured in only one or two pointings, compared to four pointings on the disk). An additional correction has been applied per detector quadrant to these mosaics to account for this by normalising by the number of times a given region has been sampled over the whole mosaic. Note that these NuSTAR mosaics, as well as all of the NuSTAR images shown here, have been smoothed with a Gaussian filter.
There are several features present in the NuSTAR images shown in Figure 1, with corresponding sources appearing also in SDO/AIA and Hinode/XRT.
-
One feature bright in both orbits is an EFR (just left of disk centre in Figure 1) which went on to become the NOAA active region 12723\(\beta\) a few days after this observation. Its properties and evolution are detailed in Section 4.
-
A number of X-ray bright points are identified, two appearing brightly in the second NuSTAR orbit (first panel of the bottom row in Figure 1 - which we label BP1 and BP2), and a fainter one (we label BP3), which is nearby some quiet Sun loops. Analysis of these features is detailed in Section 5.
-
In just the first NuSTAR orbit (first panel of the top row in Figure 1), there is a faint source that appears in the top right corner of the mosaic, which can be seen to be a jet in SDO/AIA (see Section 6).
-
The brightest source in the NuSTAR images is the large region at the east limb, seen in both NuSTAR orbits. Between the NuSTAR detectors there are chip gaps, and any photons that land on these gaps are not detected. This source is bright and extended enough to cross over the chip gaps between NuSTAR’s detector quadrants, with these detector gaps are visible in Figure 1. The bright loops in this region, clearer in the SDO/AIA and XRT images, are the remains of a decayed active region (DAR). Details and analysis are given in Section 7.
3 Analysis Methods
3.1 Fitting the NuSTAR HXR Spectra
NuSTAR is an imaging spectrometer, allowing the X-ray spectra of these quiet Sun features to be fit to investigate their properties. The spectral fitting of the features presented here was done using XSPEC (Arnaud, 1996). In order to perform spectral fitting for a given source, the spectrum and Spectral Response Matrix (SRM) - via the Response Matrix and Ancillary Response Files (RMF and ARF) - were obtained for a circular region enclosing the source (in this study, all with radii \({>}\,40''\)) using the NuSTAR Data Analysis Software. As the features investigated here produce a low number of counts, Cash statistics (Cash, 1979), the maximum likelihood-based statistic for Poisson data, were used for the fitting. For the APEC thermal model used in XSPEC, coronal abundances were assumed.
For all of the NuSTAR quiet Sun features investigated here, the HXR emission is faint (resulting in high livetimes between 67 – 92%). Although NuSTAR has a high sensitivity, its detector throughput is limited to 400 counts/s/FPM. This limited throughput, combined with HXR spectra that are sharply falling off with increasing photon energy (whether due to thermal or non-thermal continuum sources), results in the low energy counts dominating. This observed over the short duration of the mosaic pointing times (∼ 100 s), means only noisy spectra are observed with few, or no counts, above a few keV. These noisy spectra over a limited energy range are tricky to fit. Our study is helped by the recent update to NuSTAR’s calibration, which makes it possible to fit down to 2.2 keV (Madsen et al., 2021), whereas previously only down to 2.5 keV was recommended for solar observations (Grefenstette et al., 2016).
Our noisy spectra can be improved by simultaneously fitting the FPMA and FPMB spectra, introducing a multiplicative constant to the fits to account for any systematic difference between the responses of the two telescopes. This constant is a fit parameter for FPMB (so relative to FPMA value), which varies depending on this systematic uncertainty and where the source lies on the detector. As most of the features presented in this paper are very faint, using the simultaneously fitted two FPMs’ spectra is not enough. Fortunately, all of the features were captured in more than one pointing due to the overlapping mosaic tiles, so the spectra from multiple FPMs and pointings can be simultaneously fitted to improve the signal-to-noise ratio. However, it is important to find a balance between obtaining a good fit, and being able to investigate how a source evolves over time.
3.2 Reconstructing Differential Emission Measures
To investigate the multi-thermal nature of a given feature, a differential emission measure can be recovered by combining data from NuSTAR, Hinode/XRT, and the six SDO/AIA optically thin coronal-temperature channels. As with the spectral fitting, we can consider multiple pointings jointly by averaging both the NuSTAR data values and responses over the pointings. The Python version of the regularized inversion approach of Hannah and Kontar (2012) was used to reconstruct the DEMs for the features, weighted using the minimum of the EM loci curves (the data divided by the corresponding temperature response function, giving the maximum possible emission at each temperature).
The fluxes for SDO/AIA and Hinode/XRT are obtained from an image from each channel, averaged over the relevant NuSTAR pointing time. A systematic error of 20\(\%\) is assigned to the SDO/AIA, Hinode/XRT, and NuSTAR fluxes in order to account for uncertainties in their temperature responses. When calculated, the photon shot noise for SDO/AIA and Hinode/XRT was found to be negligible for these features as the data values were calculated over relatively large regions. However, in the case of NuSTAR, the calculated shot noise was not negligible, and therefore is added in quadrature with the 20\(\%\) systematic error for all DEMs shown here. For the DEM calculation, the NuSTAR data is split into two energy bands: 2.2 – 2.6 keV and 2.6 – 3.6 keV. The SDO/AIA and Hinode/XRT temperature responses were calculated using the standard SolarSoft routines from the instrument teams: aia_get_response.pro for SDO/AIA and make_xrt_temp_resp.pro for Hinode/XRT. The NuSTAR response was calculated in PythonFootnote 2 using the spectral responses used for the spectral fitting.
DEM analysis combining data from NuSTAR, SDO/AIA and Hinode/XRT has been done before by Wright et al. (2017). In this study, it was found that multiplying the Hinode/XRT responses by a factor of 2 before calculating the DEMs produced a solution with smaller residuals. This choice was made following the suggestion of previous authors who also found discrepancies when using Hinode/XRT data (Schmelz et al., 2015; Testa et al., 2011; Cheung et al., 2015). Introducing this factor was similarly found to improve the DEM results here, and so was used for all of the DEMs presented in this paper.
3.3 NuSTAR Non-thermal Upper Limits
The NuSTAR spectra of these quiet Sun features are well-fitted by an isothermal model, but as they are noisy or have no counts at higher energies (possibly due to NuSTAR’s limited detector throughput), a weak non-thermal component could be present but undetected. We can determine an upper limit on the non-thermal emission that could be present and consistent with a null detection, following the approach of Wright et al. (2017). This is done by adding a non-thermal component to the thermal model obtained from the NuSTAR spectral fitting. This thick target non-thermal model depends on three parameters: the power-law index, \(\delta\), the low-energy cutoff, \(E_{c}\), and the total electron flux, \(N_{N}\). For a chosen \(\delta\), \(E_{c}\) and \(N_{N}\) value, the resulting non-thermal model and the fitted thermal model are folded through the NuSTAR response, producing model count spectra. From these, synthetic spectra are generated through a Monte Carlo process, randomly sampling the model count spectra for the total number of counts (calculated using the livetime and duration of the observation). For a range of different \(\delta\) and \(E_{c}\) combinations, \(N_{N}\) is reduced, until these synthetic spectra lie within each others’ Poisson errors between 2 and 4 keV, and there are < 4 counts above 4 keV – consistent with a null detection to 2\(\sigma\) (Gehrels, 1986). Because the spectra considered here are noisy, this test can be repeated multiple times (1000 times in the cases discussed later), to obtain more accurate results. As well as testing models with different \(\delta\) values, the case of a mono-energetic beam of electrons, with an energy of \(E_{c}\), can also be tested. These simulations were done in Python using the thermal and thick-target models from the new solar X-ray fitting package.Footnote 3
The upper limit on \(N_{N}\) can then be used to determine an upper limit on the power in the non-thermal distribution, via:
This can then be multiplied by the duration of the observation and compared to the thermal energy, calculated as:
where \(k_{B}\) is the Boltzmann constant, V is the volume of the emitting plasma (found from the SDO/AIA or Hinode/XRT image, using source area \(A^{3/2}\)), and T and EM are the temperature and emission measure of the plasma (from the NuSTAR spectral fit), respectively (Hannah et al., 2008). Note that this equation does not take into account a loop filling factor, which means that this is an upper limit on the thermal energy. This calculated thermal energy then provides a heating requirement for the source, which if less than the upper-limits of the non-thermal power, could be the produced by the accelerated electrons, like in larger flares.
4 Emerging Flux Region
The EFR, which later went on to become an active region, was observed by NuSTAR in both orbits. For this EFR, flux first begins to emerge just after 00:00 UT on 28 September. There is some cancellation between opposite polarities, and the positive and negative polarities then spread apart. This is the time during which this feature is observed by NuSTAR. The next day, beginning at ∼ 14:00 UT on 29 September, there is more intense flux emergence in this region, producing an active region.
In each NuSTAR orbit, the EFR was captured in four pointings, split into two pairs of consecutive pointings: 12 and 13, and 18 and 19. Capturing this feature four times in each orbit, for a total number of eight pointings over the whole observation, gives the opportunity to study its temporal evolution despite the short duration of the mosaic pointings.
Hinode/XRT Be-thin and SDO/AIA 211 Å images of the EFR are shown in Figure 2. These images are from the mid-times of all of the mosaic pointings where NuSTAR captured the EFR. The NuSTAR 2.2 – 4.0 keV contours, aligned with SDO/AIA, are also plotted in this figure to show the HXR evolution of the feature. These contours indicate that the feature is generally brighter in the second orbit than the first. In the Hinode/XRT and SDO/AIA images, the EFR appears to be comprised of three separate regions. It is clear that the majority of the NuSTAR emission originates from the bright central region (\(-250'', -250''\)), particularly in the second orbit when the feature has brightened. However, there does appear to be a contribution from the upper region (\(-300'', -200''\)) to the NuSTAR emission in P12 of the first orbit, as shown in the top left panel for each instrument of Figure 2. The lower region (\(-350'', -350''\)) does not appear to contribute significantly to the NuSTAR emission. This is expected in P12 and 13 of each orbit, as this region would be outside NuSTAR’s FOV. In P18 and 19, this feature would lie close to a detector gap, which could explain the lack of NuSTAR emission from this region.
The Hinode/XRT Be-thin and SDO/AIA 211 Å lightcurves for the EFR are plotted in Figure 3, for each of the three regions separately (the boxes shown in the top left panel in Figure 2) and combined. These lightcurves confirm that the central region is the brightest of the three in both SDO/AIA and Hinode/XRT, and the greatest contributor to the emission from the EFR as a whole, as shown in the top row of Figure 3. In the first orbit, both lightcurves for this region decrease between the times of the NuSTAR pointings (18:53 UT and 19:07 UT, shown by the green shaded regions in Figure 3). Interestingly, between the two NuSTAR pointing times in the second orbit (20:29 UT and 20:44 UT) the Hinode/XRT lightcurve decreases where the SDO/AIA 211 Å lightcurve increases. However, the change in the SDO/AIA 211 Å lightcurve is relatively much smaller than the change in Hinode/XRT.
The Hinode/XRT lightcurve for the upper region shows a sharp peak in brightness during P12 of the first NuSTAR orbit, before falling off. This is consistent with the NuSTAR contours in the top left panel for each instrument in Figure 2, where the upper region is significantly contributing to the NuSTAR emission in P12 of orbit 1. The peak in the Hinode/XRT lightcurve, followed later by the peak in the SDO/AIA 211 Å lightcurve, for this region suggests a heating of material followed by cooling.
The lower region is outside NuSTAR’s FOV in P12 and 13, and would be positioned near or on a detector gap in P18 and 19. However, the lightcurves confirm that this feature is relatively faint in Hinode/XRT in comparison to the central region, and therefore would be unlikely to contribute significantly to the NuSTAR emission from the EFR.
4.1 NuSTAR Spectral Analysis
Using the approach detailed in Section 3.1, we fit an isothermal model to the NuSTAR spectra of the EFR. Multiple pointings were simultaneously fit over to reduce the noise, but it was sufficient just to use consecutive pointings (P12,13 and P18,19) so that we could still investigate the time evolution. The spectra for orbit 1 P18,19 and orbit 2 P12,13 are shown in Figure 4, and a summary of the fitting results for the EFR (in comparison with all other features investigated) is given in Table 1. The spectral fits from all four times give a reasonably constant temperature of ∼ 2.5 MK, and emission measures ranging between 1.9 and 6.1 × \(10^{44}\) cm−3. The fits do suggest a slight increase in temperature of the EFR over the two orbits, from 2.54 MK to 2.63 MK. However, taking into account the uncertainties on these temperatures, this increase is not statistically significant. The emission measure from the spectral fits decreases from 3.42 to 1.93 × \(10^{44}\) cm−3 between P12,13 and P18,19 in the first orbit. It then increases up to 6.10 × \(10^{44}\) cm−3 in P12,13 in the second orbit before falling to 2.40 × \(10^{44}\) cm−3 for P18,19.
This matches the behaviour in the Hinode/XRT lightcurve for the central region in Figure 3. In both orbits, the NuSTAR fit temperature remains approximately constant while the emission measure decreases. This is in agreement with the decreases in the Hinode/XRT lightcurve between the NuSTAR pointing times in both orbits, and the highest NuSTAR emission measure corresponds to the highest peak in Hinode/XRT at 20:30 UT. This similar behaviour is expected as Hinode/XRT and NuSTAR should be observing emission at approximately the same temperatures.
4.2 Differential Emission Measures
In the EFR lightcurves in Figure 3, at around 20:30 UT (NuSTAR orbit 2, P12 and 13), there is a peak in the Hinode/XRT Be-thin lightcurve that coincides with a minimum in the SDO/AIA 211 Å lightcurve. Later, at around 20:45 UT (NuSTAR orbit 2, P18 and 19), there is increased SDO/AIA 211 Å emission but decreased Hinode/XRT Be-thin emission. As Hinode/XRT is sensitive to higher temperature emission than SDO/AIA 211 Å, this suggests that there is more higher temperature emission present at 20:30 UT than at 20:45 UT. In order to confirm this, we perform differential emission measure (DEM) analysis using the method outlined in Section 3.2.
A comparison between the DEMs calculated for P12 and 13, and P18 and 19 is shown for each orbit in the top row of Figure 5. In orbit 1, between these two times, the emission in both the Hinode/XRT Be-thin and SDO/AIA 211 Å lightcurves in Figure 3 decreases, though the change in brightness is smaller than the 20\(\%\) systematic error used in the DEM calculation in both cases. This is reflected in the two DEMs for these times, as the DEM for P12 and 13 is higher than the one for P18 and 19 for \(\log_{10} T\gtrsim\) 6.2, though they are not significantly different when taking into account the error bars.
In the case of orbit 2, the Hinode/XRT lightcurve shows a peak at around the time of NuSTAR P12 and 13 before decreasing for NuSTAR P18 and 19. Between these two times, the decrease in Hinode/XRT emission and corresponding increase in SDO/AIA 211 Å emission suggests that there is hotter material present at the earlier time. The change in Hinode/XRT is higher than the 20\(\%\) systematic errors, whereas the change in 211 Å is again very small. From the calculated DEMs, there is more emission above \(\log_{10} T\gtrsim 6.3\) at the earlier time, in agreement with the fall-off in the Hinode/XRT lightcurve. The DEM for P18 and 19 is higher than that for P12 and 13 for \(6.1\lesssim \log_{10} T\lesssim6.3\), though there is no difference between the DEMs at these temperatures outwith the error regions. Material at these temperatures could be responsible for the peak in the SDO/AIA 211 Å at the later time.
Adding X-ray data to a DEM calculation using SDO/AIA is important to constrain the higher temperature emission, as shown in the bottom row of Figure 5. In this figure, the intersection of the NuSTAR EM loci curves, obtained by dividing the data by the NuSTAR temperature response in each energy bin, is consistent with the T and EM values from the NuSTAR spectral fit, as expected. All three DEMs are similar for \(\log_{10} T\lesssim6.4\), but demonstrate that the addition of Hinode/XRT helps to constrain the DEM at temperatures higher than this, and adding NuSTAR strengthens this constraint. From including the X-ray data, rather than relying on SDO/AIA alone, it is clear that there is virtually no emission above 4 MK here. Previous studies looking at non-flaring active regions using different DEM approaches with SDO/AIA data also produced erroneous higher temperature DEM components, removed when HXR data were included (Schmelz et al., 2009; Reale, McTiernan, and Testa, 2009). A previous study of a quiet Sun EFR by Kontogiannis et al. (2020) also included DEM analysis, using data from Hinode/EIS, finding a similar magnitude peak in the DEM at \(\log_{10} T\sim 6.1\) to those shown in Figure 5. As Kontogiannis et al. (2020) only used Hinode/EIS, their DEMs were not well constrained at higher temperatures.
4.3 NuSTAR Non-thermal Upper Limits
Using the approach detailed in Section 3.3, the upper limits on any non-thermal emission present in the EFR were calculated and compared to the required heating power obtained from the thermal energy of the plasma. This calculation was done for the peak time in the Hinode/XRT lightcurve, corresponding to NuSTAR orbit 2 P12 and 13. The central region of the EFR is ∼ \(35''\) square (see Figure 2) and, using \(V\sim A^{3/2}\), has a volume of 1.65 × \(10^{28}\) cm3. Using equation 2 and the NuSTAR spectral fit values (see Figure 4), we find a thermal energy of 3.37 × \(10^{27}\) erg, and hence, by dividing by the NuSTAR observation time (246 s), a heating power of 1.37 × \(10^{25}\) erg s−1.
We find that all of the upper limits on the non-thermal power were smaller than the heating requirement. The area used here may have been an over-estimate, and making this smaller would lower the heating requirement, but only by a small factor – this would still be an order of magnitude larger than the upper limits on the non-thermal power. Therefore, it can be concluded that, in the case of the EFR, if any non-thermal component is present, it is not responsible for the observed heating.
5 Bright Points
Three bright points were identified in the NuSTAR observations and confirmed with Hinode/XRT and SDO/AIA, as was shown in Figure 1. Zoomed-in maps of these bright points, with images in Hinode/XRT and SDO/AIA 211 Å and over-plotted NuSTAR contours, are shown in Figure 6. The top row shows bright points labelled BP1 and BP2 from NuSTAR orbit 2. In both 211 Å and Hinode/XRT, BP2 is a more compact feature than BP1, but it is brighter in NuSTAR at this time. BP1 is observed with NuSTAR in P14, 15, 16, and 17 of both orbits. BP2 is close to BP1, and lies in a region that is captured by NuSTAR also in P14, 15, 16, and 17. However, this feature is extremely faint in the first orbit, making it unusable for spectroscopy. Though it is present in all pointings in orbit 2, BP2 is located over the edge of the detector in P14 and 15, and is therefore only well observed in P16 and 17. In the bottom row of Figure 6, a fainter bright point, labelled BP3, is shown near to some larger loops in the quiet Sun (labelled QS loops). Each of these features were captured in P13, 14, 17, and 18 in both of the NuSTAR orbits.
The SDO/AIA 211 Å and Hinode/XRT Be-thin lightcurves for the BP1 and BP2 are shown in Figure 7, calculated for the boxes shown in Figure 6. In the case of BP1, the lightcurves indicate an increase in brightness in both channels, peaking just before the NuSTAR pointings in the first orbit, and then they continue to increase in brightness until the pointings in the second orbit. The two lightcurves for BP2 also show increasing brightness throughout both orbits of NuSTAR observation. This behaviour explains why BP2 is not observed by NuSTAR in the first orbit; it is not yet producing sufficiently bright emission to be detected by NuSTAR. SDO/AIA 211 and Hinode/XRT lightcurves for the QS loops and for BP3 (not shown here) show only small changes in brightness for these features between the two orbits.
5.1 NuSTAR Spectral Analysis
Figure 8 shows example NuSTAR spectra for BP1, BP2 and BP3 fitted with a thermal model. For each bright point (except BP2), and the QS loops, spectra were obtained and fitted for both orbits, and all the fitted thermal parameters are given in Table 1.
For BP1, the spectral fits give a temperature of ∼ 2.5 MK in both orbits, but with emission measure increasing from 1.12 × \(10^{44}\) cm−3 in orbit 1 (as shown in left panel of Figure 8) to 4.07 × \(10^{44}\) cm−3 in orbit 2, matching the behaviour seen in the Hinode/XRT and SDO/AIA lightcurves (Figure 7). Although, NuSTAR gives an EM increase by a factor of ∼ 4, while the feature brightens by a factor of ∼ 1.2 in SDO/AIA 211 Å and ∼ 2 in Hinode/XRT Be-thin. This could be due to both of these channels, 211 Å in particular, being sensitive to cooler material than NuSTAR, meaning that this brightening might be occurring slightly higher temperatures, and is less significant to the overall emission.
As BP2 was only well observed during NuSTAR orbit 2 P16 and 17, only this spectrum could be fitted, finding 3.22 MK and 5.98 × \(10^{43}\) cm−3, respectively (shown middle panel of Figure 8). BP3 is fainter but was well observed over both NuSTAR orbits, the fitted spectrum in orbit 1 giving 3.22 MK and 1.33 × \(10^{43}\) cm−3 (shown right panel of Figure 8), and 2.56 MK and 5.10 × \(10^{43}\) cm−3 during orbit 2. The QS loops were also well observed over both NuSTAR orbits, the fitted spectra giving similar or slightly cooler temperatures (2.07 – 2.51 MK) than the bright points (see Table 1). All the BP (and QS loop) spectra are well fitted with the isothermal model and, similarly to the previous examples, do not show any evidence of either a higher temperature or non-thermal component.
5.2 Differential Emission Measures
Following the approach of Section 3.2, DEMs were reconstructed for the bright points and are shown in Figure 9. DEMs are shown for both orbits for BP1, but only in orbit 2 for BP2 and several of the Hinode/XRT Be-thin pixels which covered BP2 were saturated at the times of these pointings, so only SDO/AIA and NuSTAR data were included in the DEM calculation for this feature (though including Hinode/XRT was found to not change the shown solution significantly). All three of the resulting DEMs for BP1 and BP2 are shown in the top row of Figure 9. They confirm the behaviour found in the NuSTAR spectral fits, with the increase in emission for BP1 between the two orbits, and confirm the presence of slightly hotter material in BP2, with a higher DEM for \(\log_{10} T\gtrsim 6.2\). The bottom row of Figure 9 plots the DEM for all 3 bright points: BP1 and BP2 during orbit 2 and BP3 during orbit 1. BP3 is fainter than BP1 and BP2, but its DEM has a similarly shaped fall for \(\log_{10} T\gtrsim 6.2\) to BP2, for both of which the NuSTAR spectral fit found a similar temperature (which was higher than that found for BP1). All the bright point DEMs shown have peaks at \(\log_{10} T\sim 5.7\) and 6.15, a result which has also been found in previous DEM analyses of coronal bright points using EUV spectroscopy (Brosius et al., 2008; Doschek et al., 2010). These previous works and the DEMs presented in this paper all show no significant emission is present above 4 MK. In our case, these DEMs benefit from using the X-ray data from NuSTAR and Hinode/XRT to constrain the higher temperature emission.
5.3 NuSTAR Non-thermal Upper Limits
For strongest bright points, BP1 and BP2, the non-thermal upper limits were calculated when these were brightest (from orbit 2), again following the approach from Section 3.3. For BP1, an area of \(35''\) square gives a volume of 1.65 × \(10^{28}\) cm3, a thermal energy of 2.71 × \(10^{27}\) erg, and a heating power of 5.02 × \(10^{24}\) erg s−1 when the energy is divided by the NuSTAR observation time (540 s). The largest non-thermal upper limits for BP1 were about an order of magnitude smaller than this heating requirement. If there was some filling factor < 1 and the area used here was an overestimate, then the heating requirement could be reduced, but at best BP1 would still be at the very limit of being a possible non-thermally heated source. BP2, with an area of \(15''\) square, gives a volume of 1.30 × \(10^{27}\) cm3, thermal energy of 3.71 × \(10^{26}\) erg, and heating power of 1.52 × \(10^{24}\) erg s−1 when this energy is divided by the NuSTAR observation time (245 s). The upper limits on the non-thermal heating power are only slightly lower than this value for a very steep, almost mono-energetic, spectrum with a low energy cutoff of ∼ 3 keV. Again, the heating requirement could be shifted down by using a filling factor < 1 and reducing the area. Therefore, it is possible that this BP2 could have been heated non-thermally, but this result is marginal.
6 Jet
A transient feature is observed at the top right of the orbit 1 NuSTAR mosaic (first panel of Figure 1), caught in pointing P4 through 7, but gone by orbit 2. SDO/AIA images of this confirm that it is a compact jet, which begins to brighten at around 18:27 UT, and has disappeared by 18:40 UT. SDO/AIA 211 Å and Hinode/XRT Be-thin images from the mid-times of each pointing are plotted in Figure 10, with aligned NuSTAR contours. From the SDO/AIA images, it can be seen that the configuration of this jet is atypical, with the jet material being ejected perpendicularly rather than radially outwards, implying that the overlying magnetic field is pushing it sideways.
The jet’s lightcurves are shown in Figure 11, with the NuSTAR count rates for each pointing plotted for comparison. When NuSTAR sees the brightest emission from the jet (P4 in Figure 11) there is also a peak in both SDO/AIA 211 Å and Hinode/XRT Be-thin. Another peak in SDO/AIA 211 Å coincides with P5, but this feature has decreased in brightness in both NuSTAR and Hinode/XRT. Again showing agreement with Hinode/XRT, though not SDO/AIA 211 Å, the NuSTAR signal is at its lowest in P6 (making the feature almost indistinguishable from the background), before increasing in brightness in P7. This behaviour is also apparent in the Hinode/XRT images, where the feature appears brighter during the times of NuSTAR P4, 5, and 7 compared to P6. The jet is positioned far enough away from any chip gaps that this change in brightness in NuSTAR is genuine, as opposed to an effect of it moving in and out of detector gaps.
6.1 NuSTAR Spectral Analysis
The NuSTAR spectra of P4, 5, and 7 were fitted simultaneously and are shown in Figure 12. P6 was not used due to the faintness of the NuSTAR emission. Though this source is clearly evolving in SDO/AIA 211 Å, the NuSTAR spectra for each of the pointings were fit individually and it was found that there was no significant change in the temperature or emission measure throughout. The fit for the NuSTAR jet spectrum gives 2.60 MK and an emission measure of 8.86 × \(10^{43}\) cm−3. The isothermal model fits the spectrum well, with no indication of a higher temperature component or any non-thermal emission. This temperature lies in the range of the sensitivity of Hinode/XRT, but slightly above the peak in SDO/AIA 211 Å, which might explain the different behaviour of SDO/AIA 211 Å compared to Hinode/XRT and NuSTAR in Figure 11.
6.2 NuSTAR Non-thermal Upper Limits
From an SDO/AIA 211 Å image of the jet, we get an area of \(3''\) square and a volume of 1.04 × \(10^{25}\) cm3. Therefore, taking the temperature and emission measure values from the NuSTAR spectral fits, the thermal energy of the jet is 3.26 × \(10^{25}\) erg. Upper limits on the possible non-thermal emission that could be present and remain undetected by NuSTAR were calculated for this event using the method discussed in Section 3.3, with the results plotted in Figure 13. Note that for the case of the mono-energetic beam and for \(\delta = 9\) and \(E_{c}\) < 4 keV the results were more well-defined that for the other cases, and therefore no spread is indicated in Figure 13. The heating power required for the jet over the time range considered here (obtained by dividing the thermal energy by the NuSTAR observation time of 528 s) is 6.18 × \(10^{22}\) erg s−1, which is marked in Figure 13 in comparison to the non-thermal upper limits. The non-thermal power would have to be greater than or equal to the heating requirement in order for the feature to have been heated through by accelerated electrons. Therefore, it can be concluded that only if the non-thermal emission was very steep, almost mono-energetic, between 3 and 4 keV could it power the required heating as determined from the NuSTAR thermal emission.
7 Bright Limb Emission
The brightest and also largest source that appears in the NuSTAR mosaics in Figure 1 is the emission from the east limb. From the SDO/AIA and Hinode/XRT images, this feature is two different sources: a bright loop, and surrounding diffuse emission. Looking back a solar rotation before this observation, this area was the site of an active region that, though too faint to be given a NOAA identification number, was detected with the Spatial Possibilities Clustering Algorithm, as SPoCA 22053 (Delouille et al., 2012). Therefore, this bright emission is likely due to the presence of a decayed active region, and so we label it DAR loop.
This source is captured fully by P11 and partially by P10 and 20 in both NuSTAR orbits. This limb emission is relatively bright and extended, meaning that the noise is not as big of an issue as it is with the other features presented in this paper. However, there are other factors, which complicate the fitting of its NuSTAR spectra. Firstly, in each of the pointings this source is extended enough to be positioned over multiple NuSTAR detector quadrants, which have varying responses. Also, as mentioned previously, the NuSTAR emission is a combination of the bright loop and the surrounding diffuse emission. Therefore, in order to investigate the properties of the emission that originates only from the bright loop in SDO/AIA and Hinode/XRT, the brightest section of the NuSTAR emission was chosen for the fitting. However, this bright loop is not ideally positioned in most of the pointings that capture this area, lying over the detector gap or just off the edge of the detector in several of these pointings. The bright region was best observed in Orbit 1 P10 and 11 with FPMA, and therefore only these pointings were combined for the spectral fitting. The NuSTAR spectrum was again fit with an isothermal model, giving a temperature of 2.53 MK and an emission measure of 9.62 × \(10^{44}\) cm−3. Because it is only possible to do a NuSTAR spectral fit for one point in time for this feature, its HXR evolution cannot be compared to its EUV and SXR evolution in SDO/AIA and Hinode/XRT.
8 Comparison of Thermal Properties
A variety of quiet Sun features have been analysed in this survey, with their X-ray emission detected for the first time by a focusing optics imaging spectrometer, allowing their X-ray spectra to be fitted. The results of fitting an isothermal model to all these features’ NuSTAR spectra are given in Table 1. The temperatures found from the NuSTAR spectral fits for all of the features lie in a narrow temperature range between 2.0 – 3.2 MK. This is due to these events being at the limit of NuSTAR’s temperature sensitivity and there is so little hotter material that the spectra are dominated by these cooler sources, appearing effectively isothermal. This is consistent with previous studies that used EUV and SXR data, which found that quiet Sun features like bright point are generally < 3 MK (Doschek et al., 2010; Alexander, Del Zanna, and Maclean, 2011; Kariyappa et al., 2011). BP2 and BP3 are slightly hotter than have been previously studied, with the NuSTAR spectral fits giving 3.2 MK. Although NuSTAR is more sensitive to higher temperature emission, there is so little of it in these features that the spectra are dominated by these cooler sources, appearing effectively isothermal, capturing the 2 – 3 MK peak of the DEM.
The DEMs for several of the features are plotted together in Figure 14 and all show very similar two peak structures, with peaks around \(\log_{10} T\sim 5.7\) and 6.1 – 6.2. This lower peak has been seen before in DEM analysis using EUV spectroscopy (Doschek et al., 2010; Brosius et al., 2008) and in our observations is dominated by the emission seen in SDO/AIA 131 Å. Increasing the lower temperature limit our DEMs are calculated over (towards \(\log_{10} T\sim 5.9\)) has minimal change on the higher temperature DEM component \(\log_{10} T\gtrsim 6.1\) and only produces a large discrepancy to the observed SDO/AIA 131 Å. Removing this channel from our DEM calculation again has minimal change to the higher temperature peak and tail of the DEMs. The higher temperature peak is slightly lower for the QS loops, jet and BP3, about \(\log_{10} T\sim 6.1\) compared to the EFR and DAR loops. Both the EFR and DAR loop have very similar DEMs, which may be just coincidence, but is curious given that one is the very start of an active region and the other the decayed remains. All DEMs fall off rapidly above this peak, highlighting that very little material has been heated to higher temperatures. The jet DEM is slightly flatter than the others, indicating possibly more hot material than in the other features, but this was not confirmed by the NuSTAR spectral fit. However, this was a small, faint and short duration event so would have been hampered by NuSTAR’s limited detector throughput, which would not have been helped by a longer dwell observation.
Difficulties also arise in this analysis when working with SDO/AIA data, because none of the channels have a peak in sensitivity in the 2 – 3 MK range. As a result, the SDO/AIA 211 Å lightcurves for these sometimes do not show behaviour consistent with Hinode/XRT or NuSTAR. Previous analysis of microflares observed with NuSTAR (for example, Cooper et al., 2020, 2021) has used the SDO/AIA Fe xviii proxy channel (Del Zanna, 2013). Unfortunately, the temperatures of these quiet Sun features are too low for this to be useful. However, Hinode/XRT has sensitivity in a similar temperature range to NuSTAR, and makes a useful comparison.
9 Discussion and Conclusions
In this paper, we have presented the first survey of quiet Sun features in HXRs observed during solar minimum. NuSTAR’s full-disk solar mosaic mode allowed for a range of different types of features to be observed. In these two mosaics, NuSTAR observed steady features, such as bright points and an EFR, but also captured a transient jet. This is the first observation of these types of features using a HXR focusing telescope. The mosaics also reveal large-scale sources (the diffuse sources in Figure 1), diagnostically important for investigating the heating of the diffuse corona.
As summarised in Section 8 and Table 1, we find the features’ temperatures lie in the range 2.0 – 3.2 MK, capturing the sharp fall off in their DEMs. We find no evidence of a higher temperature or non-thermal component present in their X-ray spectra. We have used EUV and SXR data from SDO/AIA and Hinode/XRT in addition to NuSTAR to investigate the temperature evolution of the quiet Sun features, including successfully reconstructing DEMs, which combine data from all three of these instruments. The DEM solutions for these quiet Sun features show no evidence of emission above 4 MK, a result achieved by using X-ray data in the DEM calculation to constrain the solution at high temperatures.
As all of the NuSTAR spectra were adequately fitted with an isothermal model, only non-thermal upper limits were found for some of the features. In most cases, it was found that the possible non-thermal component was not sufficient to produce the required heating. The feature that was the best candidate for non-thermal heating was the jet. However, even this would require a very steep (effectively mono-energetic) non-thermal distribution with a low energy cutoff between 3 – 4 keV.
From the spectral, DEM, and non-thermal upper limit analysis performed here, it can be concluded that no higher temperature or non-thermal sources were found in this quiet Sun data. However, if any such components were present, they would be very faint and therefore NuSTAR does not have the sensitivity required to detect them in the short 100 s mosaic pointings combined with its limited throughput. Higher temperature or non-thermal components would only be detectable if they were relatively strong, or in longer duration observations of non-transient features.
The work presented in this paper used the first NuSTAR quiet Sun campaign from the recent solar minimum. Additional data sets were taken throughout the solar minimum (2018 – 2020) in both the full-disk mosaic mode as well as longer dwells, in which pointing was not changed. In these dwells, any bright points would be observed for several hours over multiple orbits. These longer observing campaigns could increase the chances of detecting more energetic HXR emission from the quiet Sun and of capturing more atypical, harder sources. Having observations of quiet Sun features over a longer period of time will also mean that a more rigorous investigation of their temporal evolution in HXRs will be possible. The NuSTAR quiet Sun dwell data will be used to further the work presented here and will be the subject of future papers. However, shorter time-scale variability in the HXR emission from quiet Sun features such as these may remain difficult to detect until there is a dedicated solar X-ray instrument with higher sensitivity and throughput.
Data Availability
All the data used in this paper are publicly available. In particular, NuSTAR via the NuSTAR Master Catalog with the OBSIDs 90410101001 - 90410125001 and 90410201001 - 90410225001.
Notes
Summary of NuSTAR solar observations can be found at https://ianan.github.io/nsigh_all/.
References
Alexander, C.E., Del Zanna, G., Maclean, R.C.: 2011, Hinode observations and 3D magnetic structure of an X-ray bright point. Astron. Astrophys. 526, A134. DOI. ADS.
Alipour, N., Safari, H.: 2015, Statistical properties of solar coronal bright points. Astrophys. J. 807(2), 175. DOI. ADS.
Arnaud, K.A.: 1996, XSPEC: the first ten years. In: Jacoby, G.H., Barnes, J. (eds.) Astronomical Data Analysis Software and Systems V, Astronomical Society of the Pacific Conference Series 101, 17. ADS.
Brosius, J.W., Rabin, D.M., Thomas, R.J., Landi, E.: 2008, Analysis of a solar coronal bright point extreme ultraviolet spectrum from the EUNIS sounding rocket instrument. Astrophys. J. 677(1), 781. DOI. ADS.
Buitrago-Casas, J.C., Glesener, L., Christe, S., Krucker, S., Vievering, J., Athiray, P., Musset, S., Davis, L., Courtade, S., Dalton, G., et al.: 2022, On the faintest solar coronal hard x-rays observed with foxsi. Astron. Astrophys. 665, A103. DOI.
Cash, W.: 1979, Parameter estimation in astronomy through application of the likelihood ratio. Astrophys. J. 228, 939. DOI. ADS.
Cheung, M.C.M., Boerner, P., Schrijver, C.J., Testa, P., Chen, F., Peter, H., Malanushenko, A.: 2015, Thermal diagnostics with the atmospheric imaging assembly on board the solar dynamics observatory: a validated method for differential emission measure inversions. Astrophys. J. 807(2), 143. DOI. ADS.
Cooper, K., Hannah, I.G., Grefenstette, B.W., Glesener, L., Krucker, S., Hudson, H.S., White, S.M., Smith, D.M.: 2020, NuSTAR observation of a minuscule microflare in a solar active region. Astrophys. J. Lett. 893(2), L40. DOI. ADS.
Cooper, K., Hannah, I.G., Grefenstette, B.W., Glesener, L., Krucker, S., Hudson, H.S., White, S.M., Smith, D.M., Duncan, J.: 2021, NuSTAR observations of a repeatedly microflaring active region. Mon. Not. Roy. Astron. Soc. 507(3), 3936. DOI. ADS.
Del Zanna, G.: 2013, The multi-thermal emission in solar active regions. Astron. Astrophys. 558, A73. DOI. ADS.
Delouille, V., Mampaey, B., Verbeeck, C., de Visscher, R.: 2012, The SPOCA-suite: a software for extraction and tracking of active regions and coronal holes on EUV images. arXiv. ADS.
Doschek, G.A., Landi, E., Warren, H.P., Harra, L.K.: 2010, Bright points and jets in polar coronal holes observed by the extreme-ultraviolet imaging spectrometer on Hinode. Astrophys. J. 710(2), 1806. DOI. ADS.
Duncan, J., Glesener, L., Grefenstette, B.W., Vievering, J., Hannah, I.G., Smith, D.M., Krucker, S., White, S.M., Hudson, H.: 2021, NuSTAR observation of energy release in 11 solar microflares. Astrophys. J. 908(1), 29. DOI. ADS.
Gehrels, N.: 1986, Confidence limits for small numbers of events in astrophysical data. Astrophys. J. 303, 336. DOI. ADS.
Glesener, L., Krucker, S., Hannah, I.G., Hudson, H., Grefenstette, B.W., White, S.M., Smith, D.M., Marsh, A.J.: 2017, NuSTAR hard X-ray observation of a sub-a class solar flare. Astrophys. J. 845(2), 122. DOI. ADS.
Glesener, L., Krucker, S., Duncan, J., Hannah, I.G., Grefenstette, B.W., Chen, B., Smith, D.M., White, S.M., Hudson, H.: 2020, Accelerated electrons observed down to <7 keV in a NuSTAR solar microflare. Astrophys. J. Lett. 891(2), L34. DOI. ADS.
Golub, L., Krieger, A.S., Silk, J.K., Timothy, A.F., Vaiana, G.S.: 1974, Solar X-ray bright points. Astrophys. J. Lett. 189, L93. DOI. ADS.
Grefenstette, B.W., Glesener, L., Krucker, S., Hudson, H., Hannah, I.G., Smith, D.M., Vogel, J.K., White, S.M., Madsen, K.K., Marsh, A.J., Caspi, A., Chen, B., Shih, A., Kuhar, M., Boggs, S.E., Christensen, F.E., Craig, W.W., Forster, K., Hailey, C.J., Harrison, F.A., Miyasaka, H., Stern, D., Zhang, W.W.: 2016, The first focused hard X-ray images of the sun with NuSTAR. Astrophys. J. 826(1), 20. DOI. ADS.
Hannah, I.G., Kontar, E.P.: 2012, Differential emission measures from the regularized inversion of Hinode and SDO data. Astron. Astrophys. 539, A146. DOI. ADS.
Hannah, I.G., Hurford, G.J., Hudson, H.S., Lin, R.P., van Bibber, K.: 2007, First limits on the 3-200 keV X-ray spectrum of the quiet sun using RHESSI. Astrophys. J. Lett. 659(1), L77. DOI. ADS.
Hannah, I.G., Christe, S., Krucker, S., Hurford, G.J., Hudson, H.S., Lin, R.P.: 2008, RHESSI microflare statistics. II. X-ray imaging, spectroscopy, and energy distributions. Astrophys. J. 677(1), 704. DOI. ADS.
Hannah, I.G., Hudson, H.S., Hurford, G.J., Lin, R.P.: 2010, Constraining the hard X-ray properties of the quiet sun with new RHESSI observations. Astrophys. J. 724(1), 487. DOI. ADS.
Hannah, I.G., Kleint, L., Krucker, S., Grefenstette, B.W., Glesener, L., Hudson, H.S., White, S.M., Smith, D.M.: 2019, Joint X-ray, EUV, and UV observations of a small microflare. Astrophys. J. 881(2), 109. DOI. ADS.
Harrison, F.A., Craig, W.W., Christensen, F.E., Hailey, C.J., Zhang, W.W., Boggs, S.E., Stern, D., Cook, W.R., Forster, K., Giommi, P., Grefenstette, B.W., Kim, Y., Kitaguchi, T., Koglin, J.E., Madsen, K.K., Mao, P.H., Miyasaka, H., Mori, K., Perri, M., Pivovaroff, M.J., Puccetti, S., Rana, V.R., Westergaard, N.J., Willis, J., Zoglauer, A., An, H., Bachetti, M., Barrière, N.M., Bellm, E.C., Bhalerao, V., Brejnholt, N.F., Fuerst, F., Liebe, C.C., Markwardt, C.B., Nynka, M., Vogel, J.K., Walton, D.J., Wik, D.R., Alexander, D.M., Cominsky, L.R., Hornschemeier, A.E., Hornstrup, A., Kaspi, V.M., Madejski, G.M., Matt, G., Molendi, S., Smith, D.M., Tomsick, J.A., Ajello, M., Ballantyne, D.R., Baloković, M., Barret, D., Bauer, F.E., Blandford, R.D., Brandt, W.N., Brenneman, L.W., Chiang, J., Chakrabarty, D., Chenevez, J., Comastri, A., Dufour, F., Elvis, M., Fabian, A.C., Farrah, D., Fryer, C.L., Gotthelf, E.V., Grindlay, J.E., Helfand, D.J., Krivonos, R., Meier, D.L., Miller, J.M., Natalucci, L., Ogle, P., Ofek, E.O., Ptak, A., Reynolds, S.P., Rigby, J.R., Tagliaferri, G., Thorsett, S.E., Treister, E., Urry, C.M.: 2013, The nuclear spectroscopic telescope array (NuSTAR) high-energy X-ray mission. Astrophys. J. 770(2), 103. DOI. ADS.
Harvey, K.L., Strong, K.T., Nitta, N., Tsuneta, S.: 1993, Lifetimes and distribution of coronal bright points observed with Yohkoh. Adv. Space Res. 13(9), 27. DOI. ADS.
Kariyappa, R., Deluca, E.E., Saar, S.H., Golub, L., Damé, L., Pevtsov, A.A., Varghese, B.A.: 2011, Temperature variability in X-ray bright points observed with Hinode/XRT. Astron. Astrophys. 526, A78. DOI. ADS.
Kontogiannis, I., Tsiropoula, G., Tziotziou, K., Gontikakis, C., Kuckein, C., Verma, M., Denker, C.: 2020, Emergence of small-scale magnetic flux in the quiet Sun. Astron. Astrophys. 633, A67. DOI. ADS.
Kosugi, T., Matsuzaki, K., Sakao, T., Shimizu, T., Sone, Y., Tachikawa, S., Hashimoto, T., Minesugi, K., Ohnishi, A., Yamada, T., Tsuneta, S., Hara, H., Ichimoto, K., Suematsu, Y., Shimojo, M., Watanabe, T., Shimada, S., Davis, J.M., Hill, L.D., Owens, J.K., Title, A.M., Culhane, J.L., Harra, L.K., Doschek, G.A., Golub, L.: 2007, The Hinode (Solar-B) mission: an overview. Solar Phys. 243(1), 3. DOI. ADS.
Kuhar, M., Krucker, S., Glesener, L., Hannah, I.G., Grefenstette, B.W., Smith, D.M., Hudson, H.S., White, S.M.: 2018, NuSTAR detection of X-ray heating events in the quiet sun. Astrophys. J. Lett. 856(2), L32. DOI. ADS.
Lemen, J.R., Title, A.M., Akin, D.J., Boerner, P.F., Chou, C., Drake, J.F., Duncan, D.W., Edwards, C.G., Friedlaender, F.M., Heyman, G.F., Hurlburt, N.E., Katz, N.L., Kushner, G.D., Levay, M., Lindgren, R.W., Mathur, D.P., McFeaters, E.L., Mitchell, S., Rehse, R.A., Schrijver, C.J., Springer, L.A., Stern, R.A., Tarbell, T.D., Wuelser, J.-P., Wolfson, C.J., Yanari, C., Bookbinder, J.A., Cheimets, P.N., Caldwell, D., Deluca, E.E., Gates, R., Golub, L., Park, S., Podgorski, W.A., Bush, R.I., Scherrer, P.H., Gummin, M.A., Smith, P., Auker, G., Jerram, P., Pool, P., Soufli, R., Windt, D.L., Beardsley, S., Clapp, M., Lang, J., Waltham, N.: 2012, The Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO). Solar Phys. 275(1–2), 17. DOI. ADS.
Lin, R.P., Dennis, B.R., Hurford, G.J., Smith, D.M., Zehnder, A., Harvey, P.R., Curtis, D.W., Pankow, D., Turin, P., Bester, M., Csillaghy, A., Lewis, M., Madden, N., van Beek, H.F., Appleby, M., Raudorf, T., McTiernan, J., Ramaty, R., Schmahl, E., Schwartz, R., Krucker, S., Abiad, R., Quinn, T., Berg, P., Hashii, M., Sterling, R., Jackson, R., Pratt, R., Campbell, R.D., Malone, D., Landis, D., Barrington-Leigh, C.P., Slassi-Sennou, S., Cork, C., Clark, D., Amato, D., Orwig, L., Boyle, R., Banks, I.S., Shirey, K., Tolbert, A.K., Zarro, D., Snow, F., Thomsen, K., Henneck, R., McHedlishvili, A., Ming, P., Fivian, M., Jordan, J., Wanner, R., Crubb, J., Preble, J., Matranga, M., Benz, A., Hudson, H., Canfield, R.C., Holman, G.D., Crannell, C., Kosugi, T., Emslie, A.G., Vilmer, N., Brown, J.C., Johns-Krull, C., Aschwanden, M., Metcalf, T., Conway, A.: 2002, The Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). Solar Phys. 210(1), 3. DOI. ADS.
Madjarska, M.S.: 2019, Coronal bright points. Living Rev. Solar Phys. 16(1), 2. DOI. ADS.
Madsen, K.K., Forster, K., Grefenstette, B.W., Harrison, F.A., Miyasaka, H.: 2021, Effective area calibration of the Nuclear Spectroscopic Telescope ARray (NuSTAR). arXiv. ADS.
Mou, C., Madjarska, M.S., Galsgaard, K., Xia, L.: 2018, Eruptions from quiet Sun coronal bright points. I. Observations. Astron. Astrophys. 619, A55. DOI. ADS.
Parker, E.N.: 1988, Nanoflares and the solar X-ray corona. Astrophys. J. 330, 474. DOI. ADS.
Price-Whelan, A.M., Sipőcz, B.M., Günther, H.M., Lim, P.L., Crawford, S.M., Conseil, S., Shupe, D.L., Craig, M.W., Dencheva, N., Ginsburg, A., VanderPlas, J.T., Bradley, L.D., Pérez-Suárez, D., de Val-Borro, M., Aldcroft, T.L., Cruz, K.L., Robitaille, T.P., Tollerud, E.J., Ardelean, C., Babej, T., Bach, Y.P., Bachetti, M., Bakanov, A.V., Bamford, S.P., Barentsen, G., Barmby, P., Baumbach, A., Berry, K.L., Biscani, F., Boquien, M., Bostroem, K.A., Bouma, L.G., Brammer, G.B., Bray, E.M., Breytenbach, H., Buddelmeijer, H., Burke, D.J., Calderone, G., Cano Rodríguez, J.L., Cara, M., Cardoso, J.V.M., Cheedella, S., Copin, Y., Corrales, L., Crichton, D., D’Avella, D., Deil, C., Depagne, É., Dietrich, J.P., Donath, A., Droettboom, M., Earl, N., Erben, T., Fabbro, S., Ferreira, L.A., Finethy, T., Fox, R.T., Garrison, L.H., Gibbons, S.L.J., Goldstein, D.A., Gommers, R., Greco, J.P., Greenfield, P., Groener, A.M., Grollier, F., Hagen, A., Hirst, P., Homeier, D., Horton, A.J., Hosseinzadeh, G., Hu, L., Hunkeler, J.S., Ivezić, Ž., Jain, A., Jenness, T., Kanarek, G., Kendrew, S., Kern, N.S., Kerzendorf, W.E., Khvalko, A., King, J., Kirkby, D., Kulkarni, A.M., Kumar, A., Lee, A., Lenz, D., Littlefair, S.P., Ma, Z., Macleod, D.M., Mastropietro, M., McCully, C., Montagnac, S., Morris, B.M., Mueller, M., Mumford, S.J., Muna, D., Murphy, N.A., Nelson, S., Nguyen, G.H., Ninan, J.P., Nöthe, M., Ogaz, S., Oh, S., Parejko, J.K., Parley, N., Pascual, S., Patil, R., Patil, A.A., Plunkett, A.L., Prochaska, J.X., Rastogi, T., Reddy Janga, V., Sabater, J., Sakurikar, P., Seifert, M., Sherbert, L.E., Sherwood-Taylor, H., Shih, A.Y., Sick, J., Silbiger, M.T., Singanamalla, S., Singer, L.P., Sladen, P.H., Sooley, K.A., Sornarajah, S., Streicher, O., Teuben, P., Thomas, S.W., Tremblay, G.R., Turner, J.E.H., Terrón, V., van Kerkwijk, M.H., de la Vega, A., Watkins, L.L., Weaver, B.A., Whitmore, J.B., Woillez, J., Zabalza, V. (Astropy Contributors, Astropy Collaboration): 2018, The Astropy project: building an open-science project and status of the v2.0 core package. Astron. J. 156, 123. DOI. ADS.
Reale, F., McTiernan, J.M., Testa, P.: 2009, Comparison of Hinode/XRT and RHESSI detection of hot plasma in the non-flaring solar corona. Astrophys. J. Lett. 704(1), L58. DOI. ADS.
Robitaille, T.P., Tollerud, E.J., Greenfield, P., Droettboom, M., Bray, E., Aldcroft, T., Davis, M., Ginsburg, A., Price-Whelan, A.M., Kerzendorf, W.E., Conley, A., Crighton, N., Barbary, K., Muna, D., Ferguson, H., Grollier, F., Parikh, M.M., Nair, P.H., Unther, H.M., Deil, C., Woillez, J., Conseil, S., Kramer, R., Turner, J.E.H., Singer, L., Fox, R., Weaver, B.A., Zabalza, V., Edwards, Z.I., Azalee Bostroem, K., Burke, D.J., Casey, A.R., Crawford, S.M., Dencheva, N., Ely, J., Jenness, T., Labrie, K., Lim, P.L., Pierfederici, F., Pontzen, A., Ptak, A., Refsdal, B., Servillat, M., Streicher, O. (Astropy Collaboration): 2013, Astropy: a community Python package for astronomy. Astron. Astrophys. 558, A33. DOI. ADS
Schmelz, J.T., Kashyap, V.L., Saar, S.H., Dennis, B.R., Grigis, P.C., Lin, L., De Luca, E.E., Holman, G.D., Golub, L., Weber, M.A.: 2009, Some like it hot: coronal heating observations from Hinode X-ray telescope and RHESSI. Astrophys. J. 704(1), 863. DOI. ADS.
Schmelz, J.T., Asgari-Targhi, M., Christian, G.M., Dhaliwal, R.S., Pathak, S.: 2015, Hot plasma from solar active region cores: a test of AC and DC coronal heating models? Astrophys. J. 806(2), 232. DOI. ADS.
Shimojo, M., Shibata, K.: 1999, Occurrence rate of microflares in an X-ray-bright point within an active region. Astrophys. J. 516(2), 934. DOI. ADS.
SunPy Community, T., Mumford, S.J., Christe, S., Pérez-Suárez, D., Ireland, J., Shih, A.Y., Inglis, A.R., Liedtke, S., Hewett, R.J., Mayer, F., Hughitt, K., Freij, N., Meszaros, T., Bennett, S.M., Malocha, M., Evans, J., Agrawal, A., Leonard, A.J., Robitaille, T.P., Mampaey, B., Iván Campos-Rozo, J., Kirk, M.S.: 2015, SunPy—Python for solar physics. Comput. Sci. Discov. 8(1), 014009. DOI. ADS.
Testa, P., Reale, F., Landi, E., DeLuca, E.E., Kashyap, V.: 2011, Temperature distribution of a non-flaring active region from simultaneous Hinode XRT and EIS observations. Astrophys. J. 728(1), 30. DOI. ADS.
van Driel-Gesztelyi, L., Green, L.M.: 2015, Evolution of active regions. Living Rev. Solar Phys. 12(1), 1. DOI. ADS.
Wright, P.J., Hannah, I.G., Grefenstette, B.W., Glesener, L., Krucker, S., Hudson, H.S., Smith, D.M., Marsh, A.J., White, S.M., Kuhar, M.: 2017, Microflare heating of a solar active region observed with NuSTAR, Hinode/XRT, and SDO/AIA. Astrophys. J. 844(2), 132. DOI. ADS.
Zhang, J., Kundu, M.R., White, S.M.: 2001, Spatial distribution and temporal evolution of coronal bright points. Solar Phys. 198(2), 347. DOI. ADS.
Acknowledgments
This paper made use of data from the NuSTAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, funded by the National Aeronautics and Space Administration. We thank the NuSTAR Operations, Software and Calibration teams for support with the execution and analysis of these observations. This research made use of the NuSTAR Data Analysis Software (NUSTARDAS), jointly developed by the ASI Science Data Center (ASDC, Italy) and the California Institute of Technology (USA). Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as domestic partner and NASA and UKSA as international partners. It is operated by these agencies in co-operation with ESA and NSC (Norway). AIA on the Solar Dynamics Observatory is part of NASA’s Living with a Star program. This research has made use of SunPy, an open-source and free community-developed solar data analysis package written in Python (SunPy Community et al., 2015). This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al., 2018, 2013). The authors would also like to thank Kathy Reeves for co-ordinating the Hinode/XRT observations.
Funding
SP acknowledges support from the UK’s Science and Technology Facilities Council (STFC) doctoral training grant (ST/T506102/1). IGH acknowledges support from a Royal Society University Fellowship (URF/R/180010) and STFC grant (ST/T000422/1).
Author information
Authors and Affiliations
Contributions
SP did the data processing and applied the analysis techniques described and wrote the manuscript text. IGH developed the demreg analysis techniques and helped with the spectral analysis. All authors revised the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Paterson, S., Hannah, I.G., Grefenstette, B.W. et al. The First Survey of Quiet Sun Features Observed in Hard X-Rays with NuSTAR. Sol Phys 298, 47 (2023). https://doi.org/10.1007/s11207-023-02135-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11207-023-02135-4