The almost unexplored frequency window from submillimeter to mid-infrared (mid-IR) may bring new clues about the particle acceleration and transport processes and the atmospheric thermal response during solar flares. Because of its technical complexity and the special atmospheric environment needed, observations at these frequencies are very sparse. The High Altitude THz Solar Photometer (HATS) is a full-Sun ground-based telescope designed to observe the continuum from the submillimeter to the mid-IR. It has a 457-mm spherical mirror with the sensor in its primary focus. The sensor is a Golay cell with high sensitivity in a very wide frequency range. The telescope has a polar mount, and a custom-built data acquisition system based on a 32 ksamples per second, 24 bits (72 dB dynamic range), 8 channels analog-to-digital board. Changing only the composition of the low- and band-pass filters in front of the Golay cell, the telescope can be setup to detect very different frequency bands; making the instrument very versatile. In this article we describe the telescope characteristics and its development status. Moreover, we give estimates of the expected fluxes during flares.

To better understand the influence of the activity cycle on the solar atmosphere, we report the time variation of the radius observed at 37 GHz ($\lambda$=8.1 mm) obtained by the Mets\"ahovi Radio Observatory (MRO) through Solar Cycles 22 to 24 (1989-2015). Almost 5800 maps were analyzed, however, due to instrumental setups changes the data set showed four distinct behaviors, which requested a normalisation process to allow the whole interval analysis. When the whole period was considered, the results showed a positive correlation index of 0.17 between the monthly means of the solar radius at 37 GHz and solar flux obtained at 10.7 cm (F10.7). This correlation index increased to 0.44, when only the data obtained during the last period without instrumental changes were considered (1999-2015). The solar radius correlation with the solar cycle agrees with the previous results obtained at mm/cm wavelengths (17 and 48 GHz), nevertheless, this result is the opposite of that reported at submillimetre wavelengths (212 and 405 GHz).

The Neupert effect is the empirical observation that the time evolution of non-thermal emission (e.g. hard X-rays) is frequently proportional to the time derivative of the thermal emission flux (soft X-rays), or, vice versa, that time integrated non-thermal flux is proportional to thermal flux. We analyzed the GOES M2.2 event SOL2011-02-14T17:25, and found that the 212 GHz emission plays quite well the role of the thermal component of the Neupert effect. We show that the maximum of the hard X-ray flux for energies above 50 keV is coincident in time with the time-derivative of the 212 GHz flux, within the uncertainties. The microwave flux density at 15.4 GHz, produced by optically thin gyrosynchrotron mechanism, and hard-X rays above 25 keV mark the typical impulsive phase, and have similar time evolution. On the other hand, the 12 GHz emission is delayed by about 25 seconds with respect of the microwave and hard X-ray peak. We argue that this delay cannot be explained by magnetic trapping of non-thermal electrons. With all the observational evidence, we suggest that the 212 GHz emission is produced by thermal bremsstrahlung, initially in the chromosphere, and shifting to optically thin emission from thehot coronal loops at the end of the gradual phase.

We use the Monte Carlo particle physics code FLUKA (Fluktuierende Kaskade) to calculate $\gamma$-ray spectra expected from solar flare energetic ion distributions. The FLUKA code includes robust physics-based models for electromagnetic, hadronic and nuclear interactions, sufficiently detailed for it to be a useful tool for calculating nuclear de-excitation, positron annihilation and neutron capture line fluxes and shapes, as well as $\approx \; {\rm GeV}$ continuum radiation from pion decay products. We show nuclear de-excitation $\gamma$-ray line model spectra from a range of assumed primary accelerated ion distributions and find them to be in good agreement with those found using the code of Murphy et al. (2009). We also show full $\gamma$-ray model spectra which exhibit all the typical structures of $\gamma$-ray spectra observed in solar flares. From these model spectra we build templates which are incorporated into the software package Objective Spectral Executive (OSPEX) and used to fit the combined Fermi Gamma-ray Burst Monitor (GBM)/Large Area Telescope (LAT) spectrum of the 2010 June 12 solar flare, providing a statistically acceptable result. To the best of our knowledge, the fit carried out with the FLUKA templates for the full $\gamma$-ray spectrum can be regarded as the first attempt to use a single code to implement a self-consistent treatment of the several spectral components in the photon energy range from $\approx 100$s ${\rm keV}$ to $\approx 100$s ${\rm MeV}$.

Polar brightening of the Sun at radio frequencies has been studied for almost fifty years and yet a disagreement persists between solar atmospheric models and observations. Some observations reported brightening values much smaller than the expected values obtained from the models, with discrepancies being particularly large at millimeter wavelengths. New clues to calibrate the atmospheric models can be obtained with the advent of the Atacama Large Millimeter/submillimeter Array (ALMA) radio interferometer. In this work, we analyzed the lower limit of the polar brightening observed at 100 and 230 GHz by ALMA, during its Science Verification period, 2015 December 16-20. We find that the average polar intensity is higher than the disk intensity at 100 and 230 GHz, with larger brightness intensities at the South pole in eight of the nine maps analyzed. The observational results were compared with calculations of the millimetric limb brightnening emission for two semi-empirical atmospheric models, FAL- C (Fontenla et al. 1993) and SSC (Selhorst et al. 2005a). Both models presented larger limb intensities than the averaged observed values. The intensities obtained with the SSC model were closer to the observations, with polar brightenings of 10.5% and 17.8% at 100 and 230 GHz, respectively. This discrepancy may be due to the presence of chromospheric features (like spicules) at regions close to the limb.

Radio-bright regions near the solar poles are frequently observed in Nobeyama Radioheliograph (NoRH) maps at 17 GHz, and often in association with coronal holes. However, the origin of these polar brightening has not been established yet. We propose that small magnetic loops are the source of these bright patches, and present modeling results that reproduce the main observational characteristics of the polar brightening within coronal holes at 17 GHz. The simulations were carried out by calculating the radio emission of the small loops, with several temperature and density profiles, within a 2D coronal hole atmospheric model. If located at high latitudes, the size of the simulated bright patches are much smaller than the beam size and they present the instrument beam size when observed. The larger bright patches can be generated by a great number of small magnetic loops unresolved by the NoRH beam. Loop models that reproduce bright patches contain denser and hotter plasma near the upper chromosphere and lower corona. On the other hand, loops with increased plasma density and temperature only in the corona do not contribute to the emission at 17 GHz. This could explain the absence of a one-to-one association between the 17 GHz bright patches and those observed in extreme ultraviolet. Moreover, the emission arising from small magnetic loops located close to the limb may merge with the usual limb brightening profile, increasing its brightness temperature and width.

Observations of the Sun with the Atacama Large Millimeter Array have now started, and the thermal infrared will regularly be accessible from the NSF's Daniel K. Inouye Solar Telescope. Motivated by the prospect of these new data, and by recent flare observations in the mid infrared, we set out here to model and understand the source of the infrared continuum in flares, and to explore its diagnostic capability for the physical conditions in the flare atmosphere. We use the 1D radiation hydrodynamics code RADYN to calculate mid-infrared continuum emission from model atmospheres undergoing sudden deposition of energy by non-thermal electrons. We identify and characterise the main continuum thermal emission processes relevant to flare intensity enhancement in the mid- to far-infrared (2-200 $\mu$m) spectral range as free-free emission on neutrals and ions. We find that the infrared intensity evolution tracks the energy input to within a second, albeit with a lingering intensity enhancement, and provides a very direct indication of the evolution of the atmospheric ionization. The prediction of highly impulsive emission means that, on these timescales, the atmospheric hydrodynamics need not be considered in analysing the mid-IR signatures.

We present an analysis of intermittent processes occurred during the impulsive phase of the flare SOL2012-03-13, using hard X-rays and submillimeter radio data. Intermittency is a key characteristic in turbulent plasmas and have been a analyzed recently for hard X-rays data only. Since in a typical flare the same accelerated electron population is believed to produce both hard X-rays and gyrosynchrotron, we compare both time profiles searching for intermittency signatures. For that we define a cross-wavelet power spectrum, that is used to obtain the Local Intermittency Measure or LIM. When greater than three, the square LIM coefficients indicate a local intermittent process. The LIM$^2$ coefficient distribution in time and scale helps to identify avalanche or cascade energy release processes. We find two different and well separated intermittent behaviors in the submillimeter data: for scales greater than 20 s, a broad distribution during the rising and maximum phases of the emission seems to favor a cascade process; for scales below 1 s, short pulses centered on the peak time, are representative of avalanches. When applying the same analysis to hard X-rays, we find only the scales above 10 s producing a distribution related to a cascade energy fragmentation. Our results suggest that different acceleration mechanisms are responsible for tens of keV and MeV energy ranges of electrons.

The present solar cycle is particular in many aspects: it had a delayed rising phase, it is the weakest of the last 100 years, and it presents two peaks separated by more than one year. To understand the impact of these characteristics on the solar chromosphere and coronal dynamics, images from a wide wavelength range are needed. In this work we use the 17~GHz radio continuum, formed in the upper chromosphere and the EUV lines 304 and 171~{\AA}, that come from the transition region (He II) and the corona (Fe IX, X), respectively. We analyze daily images at 304 and 171~{\AA} obtained by the Atmospheric Imaging Assembly (AIA). The 17~GHz maps were obtained by the Nobeyama Radioheliograph (NoRH). To construct synoptic limb charts, we calculated the mean emission of delimited limb areas with 100" wide and angular separation of $5^\circ$. At the equatorial region, the results show an hemispheric asymmetry of the solar activity. The northern hemisphere dominance is coincident with the first sunspot number peak, whereas the second peak occurs concurrently with the increase in the activity at the south. The polar emission reflects the presence of coronal holes at both EUV wavelengths, moreover, the 17~GHz polar brightenings can be associated with the coronal holes. Until 2013, both EUV coronal holes and radio polar brightenings were more predominant at the south pole. Since then they have not been apparent in the north, but thus appear in the beginning of 2015 in the south as observed in the synoptic charts. This work strengthens the association between coronal holes and the 17~GHz polar brightenings as it is evident in the synoptic limb charts, in agreement with previous case study papers. The enhancement of the radio brightness in coronal holes is explained by the presence of bright patches closely associated with the presence of intense unipolar magnetic fields.

We report the statistics of the number of active regions (NAR) observed at 17 GHz with the Nobeyama Radioheliograph between 1992, near the maximum of cycle 22, and 2013, that also includes the maximum of cycle 24, and we compare with other activity indexes. We find that NAR minima are shorter than those of the sunspot number (SSN) and radio flux at 10.7 cm (F10.7). This shorter NAR minima could reflect the presence of active regions generated by faint magnetic fields or spotless regions, which were a considerable fraction of the counted active regions. The ratio between the solar radio indexes F10.7/NAR shows a similar reduction during the two minima analyzed, which contrasts with the increase of the ratio of both radio indexes in relation to the SSN during the minimum of cycle 23-24. These results indicate that the radio indexes are more sensitive to weaker magnetic fields than those necessary to form sunspots, of the order of 1500 G. The analysis of the monthly averages of the active region brightness temperatures shows that its long term variation mimics the solar cycle, although, due to the gyro-resonance emission, a great number of intense spikes are observed in the maximum temperature study. The decrease, in number, of these spikes is also evident during the current cycle 24, a consequence of the sunspot magnetic field weakening in the last years.

Ground- and space-based observations of solar flares from radio wavelengths to gamma-rays have produced considerable insights but raised several unsolved controversies. The last unexplored wavelength frontier for solar flares is in the range of submillimeter and infrared wavelengths. Here we report the detection of an intense impulsive burst at 30 THz using a new imaging system. The 30 THz emission exhibited remarkable time coincidence with peaks observed at microwave, mm/submm, visible, EUV and hard X-ray wavelengths. The emission location coincides with a very weak white-light feature, and is consistent with heating below the temperature minimum in the atmosphere. However, there are problems in attributing the heating to accelerated electrons. The peak 30 THz flux is several times larger than the usual microwave peak near 9 GHz, attributed to non-thermal electrons in the corona. The 30 THz emission could be consistent with an optically thick spectrum increasing from low to high frequencies. It might be part of the same spectral component found at sub-THz frequencies whose nature remains mysterious. Further observations at these wavelengths will provide a new window for flare studies.

We study a solar flare that occurred on September 10, 2002, in active region NOAA 10105 starting around 14:52 UT and lasting approximately 5 minutes in the radio range. The event was classified as M2.9 in X-rays and 1N in H\alpha. Solar Submillimeter Telescope observations, in addition to microwave data give us a good spectral coverage between 1.415 and 212 GHz. We combine these data with ultraviolet images, hard and soft X-rays observations and full-disk magnetograms. Images obtained from Ramaty High Energy Solar Spectroscopic Imaging data are used to identify the locations of X-ray sources at different energies and to determine the X-ray spectrum, while ultra violet images allow us to characterize the coronal flaring region. The magnetic field evolution of the active region is analyzed using Michelson Doppler Imager magnetograms. The burst is detected at all available radio-frequencies. X-ray images (between 12 keV and 300 keV) reveal two compact sources and 212 GHz data, used to estimate the radio source position, show a single compact source displaced by 25" from one of the hard X-ray footpoints. We model the radio spectra using two homogeneous sources, and combine this analysis with that of hard X-rays to understand the dynamics of the particles. Relativistic particles, observed at radio wavelengths above 50 GHz, have an electron index evolving with the typical soft-hard-soft behaviour.

Radio and optical observations of the evolution of flare-associated phenomena have shown an initial and rapid burst at 0.4 THz only followed subsequently by a localized chromospheric heating producing an H{\alpha} brightening with later heating of the whole active region. A major instability occurred several minutes later producing one impulsive burst at microwaves only, associated with an M2.0 GOES X-ray flare that exhibited the main H{\alpha} brightening at the same site as the first flash. The possible association between long-enduring time profiles at soft X-rays, microwaves, H{\alpha} and sub-THz wavelengths is discussed. In the decay phase the H{\alpha} movie shows a disrupting magnetic arch structure ejecting dark, presumably chromospheric, material upwards. The time sequence of events suggests genuine interdependent and possibly non-thermal instabilities triggering phenomena, with concurrent active region plasma heating and material ejection.

The GOES X1.5 class flare that occurred on August 30,2002 at 1327:30 UT is one of the few events detected so far at submillimeter wavelengths. We present a detailed analysis of this flare combining radio observations from 1.5 to 212 GHz (an upper limit of the flux is also provided at 405 GHz) and X-ray. Although the observations of radio emission up to 212 GHz indicates that relativistic electrons with energies of a few MeV were accelerated, no significant hard X-ray emission was detected by RHESSI above ~ 250 keV. Images at 12--20 and 50--100 keV reveal a very compact, but resolved, source of about ~ 10" x 10". EUV TRACE images show a multi-kernel structure suggesting a complex (multipolar) magnetic topology. During the peak time the radio spectrum shows an extended flatness from ~ 7 to 35 GHz. Modeling the optically thin part of the radio spectrum as gyrosynchrotron emission we obtained the electron spectrum (spectral index delta, instantaneous number of emitting electrons). It is shown that in order to keep the expected X-ray emission from the same emitting electrons below the RHESSI background at 250 keV, a magnetic field above 500 G is necessary. On the other hand, the electron spectrum deduced from radio observations >= 50 GHz is harder than that deduced from ~ 70 - 250 keV X-ray data, meaning that there must exist a breaking energy around a few hundred keV. During the decay of the impulsive phase, a hardening of the X-ray spectrum is observed which is interpreted as a hardening of the electron distribution spectrum produced by the diffusion due to Coulomb collisions of the trapped electrons in a medium with an electron density of n_e ~ 3E10 - 5E10 cm-3.

A new solar burst emission spectral component has been found showing sub-THz fluxes increasing with frequency, spectrally separated from the well known microwave component. Rapid pulsations are found present in all events observed at the two frequencies of the solar submillimeter-wave telescope (SST): 212 and 405 GHz. They were studied in greater detail for three solar bursts exhibiting the new THz spectral component. The pulse amplitudes are of about 5-8% of the mean flux throughout the bursts durations, being comparable for both frequencies. Pulsations range from one pulse every few seconds to 8-10 per second. The pulse repetition rates (R) are linearly proportional to the mean burst fluxes (S), following the simple relationship S = k R, suggesting that the pulsations might be the response to discrete flare particle accelerator injections quantized in energy. Although this result is consistent with qualitative trends previously found in the GHz range, the pulse amplitude relative to the mean fluxes at the sub-THz frequencies appear to be nearly ten times smaller than expected from the extrapolation of the trends found in the GHz range. However there are difficulties to reconcile the nearly simultaneous GHz and THz burst emission spectrally separated components, exhibiting rapid pulsations with considerably larger relative intensities in the GHz range.

The presence of a solar burst spectral component with flux density increasing with frequency in the sub-terahertz range, spectrally separated from the well-known microwave spectral component, bring new possibilities to explore the flaring physical processes, both observational and theoretical. The solar event of 6 December 2006, starting at about 18:30 UT, exhibited a particularly well-defined double spectral structure, with the sub-THz spectral component detected at 212 and 405 GHz by SST and microwaves (1-18 GHz) observed by the Owens Valley Solar Array (OVSA). Emissions obtained by instruments in satellites are discussed with emphasis to ultra-violet (UV) obtained by the Transition Region And Coronal Explorer (TRACE), soft X-rays from the Geostationary Operational Environmental Satellites (GOES) and X- and gamma-rays from the Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The sub-THz impulsive component had its closer temporal counterpart only in the higher energy X- and gamma-rays ranges. The spatial positions of the centers of emission at 212 GHz for the first flux enhancement were clearly displaced by more than one arc-minute from positions at the following phases. The observed sub-THz fluxes and burst source plasma parameters were found difficult to be reconciled to a purely thermal emission component. We discuss possible mechanisms to explain the double spectral components at microwaves and in the THz ranges.

Aims. To determine the solar transition region and coronal radius at EUV wavelengths and its time evolution during Solar Cycle XXIII. Methods. We use daily 30.4 and 17.1 nm images obtained by the Extreme Ultraviolet Imager (EIT) aboard the SoHO satellite and derive the solar radius by fitting a circle to the limb brightness ring. Results. The weighted mean of the temporal series gives (967''.56 +/- 0''.04) and (969''.54 +/- 0''.02) at 30.4 and 17.1 nm respectively. No significant correlation was found with the solar cycle at any of the two wavelengths. Conclusions. Since the temperature formation of the 30.4 nm line is between (60 - 80) 10^3 K (Transition Region), the obtained result is bigger than that derived from present atmospheric models. On the contrary this height is compatible with radio models.

We analyse the 30 October, 2004, X1.2/SF solar event that occurred in AR 10691 (N13 W18) at around 11:44 UT. Observations at 212 and 405 GHz of the Solar Submillimeter Telescope (SST), with high time resolution (5 ms), show an intense impulsive burst followed by a long-lasting thermal phase. EUV images from the Extreme Ultraviolet Imaging Telescope (SOHO/EIT) are used to identify the possible emitting sources. Data from the Radio Solar Telescope Network (RSTN) complement our spectral observations below 15 GHz. During the impulsive phase the turnover frequency is above 15.4 GHz. The long-lasting phase is analysed in terms of thermal emission and compared with GOES observations. From the ratio between the two GOES soft X-ray bands, we derive the temperature and emission measure, which is used to estimate the free-free submillimeter flux density. Good temporal agreement is found between the estimated and observed profiles, however the former is larger than the latter.