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.

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.

The altitude of the ionospheric lower layer (D-region) is highly influenced by the solar UV flux affecting in turn, the propagation of Very Low Frequency (VLF) signals inside the waveguide formed between this layer and the Earth surface. A rapid change of the solar irradiance, as during a solar eclipse, can help to understand the details of the energy transfer of the solar radiation onto the ionospheric D-layer. Using the "Latin American VLF Network" (LAVNet-Mex) receiver station in Mexico City, Mexico, we detected the phase and amplitude changes of the VLF signals transmitted by the NDK station at 25.2 kHz in North Dakota, USA during the August 21, 2017, solar eclipse. As the Sunlight was eclipsed, the rate of ionization in the ionosphere (D-region) was reduced and the effective reflection height increased, causing a considerable drop of the phase and amplitude of the observed VLF waves. The corresponding waveguide path is 3007.15 km long and crossed almost perpendicularly the total eclipse path. Circumstantially, at the time of the total eclipse, a C3 flare took place allowing us to isolate the flare flux from the background flux of a large portion of the disk. In this work, we report the observations and present a new model of the ionospheric effects of the eclipse and flare. The model is based on a detailed setup of the degree of Moon shadow that affects the entire Great Circle Path (GCP). During the eclipse, the maximum phase variation was -63.36$^{\circ}$ at 18:05 UT which, according to our model, accounts for a maximum increase of the reflection height of 9.3 km.

We report on the 22-23 June 2015 geomagnetic storm. There has been a shortage of intense geomagnetic storms during the current solar cycle 24 in relation to the previous cycle. This situation changed after mid-June 2015 when one of the biggest solar active regions (AR 2371) of current solar cycle 24, close to the central meridian produced several coronal mass ejections (CMEs) associated with M-class flares. The CMEs impact on the Earth's magnetosphere resulted in a moderately-severe G4-class geomagnetic storm on 22-23 June 2015 and a G2 (moderate) geomagnetic storms on 24 June. The G4 solstice storm was the second biggest (so far) geomagnetic storms of cycle 24. We highlight the ground level observations made by New-Tupi, Muonca and the CARPET El Leoncito cosmic ray detectors that are located within the South Atlantic Anomaly (SAA) region. These observations are studied in correlation with data obtained by space-borne detectors and other ground-based experiments. The CME designations are from the Computer Aided CME Tracking (CACTus) automated catalog. As expected, Forbush Decreases (FD) associated with the passing CMEs were recorded by these detectors. We noticed a peculiar feature linked to a severe geomagnetic storm event. The 21 June 2015 CME 0091 was likely associated with the 22 June summer solstice FD event. The angular width of CME 0091 was very narrow and measured 56 degrees seen from Earth. In most cases, only CME halos and partial halos, lead to severe geomagnetic storms. We performed a cross-check analysis of the FD events detected during the rise phase of the current solar cycle 24, the geomagnetic parameters, and the CACTus CME catalog. Our study suggests that narrow angular-width CMEs that erupt in the western region of the ecliptic plane can lead to moderate and severe geomagnetic storms.

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 recent discovery of impulsive solar burst emission in the 30 THz band is raising new interpretation challenges. One event associated with a GOES M2 class flare has been observed simultaneously in microwaves, H-alpha, EUV, and soft X-ray bands. Although these new observations confirm some features found in the two prior known events, they exhibit time profile structure discrepancies between 30 THz, microwaves, and hard X-rays (as inferred from the Neupert effect). These results suggest a more complex relationship between 30 THz emission and radiation produced at other wavelength ranges. The multiple frequency emissions in the impulsive phase are likely to be produced at a common flaring site lower in the chromosphere. The 30 THz burst emission may be either part of a nonthermal radiation mechanism or due to the rapid thermal response to a beam of high-energy particles bombarding the dense solar atmosphere.

Solar observations in the infrared domain can bring important clues on the response of the low solar atmosphere to primary energy released during flares. At present the infrared continuum has been detected at 30 THz (10 $\mu$m) in only a few flares. In this work we present a detailed multi-frequency analysis of SOL2012-03-13, including observations at radio millimeter and sub-millimeter wavelengths, in hard X-rays (HXR), gamma-rays (GR), H-alpha, and white-light. HXR/GR spectral analysis shows that the event is a GR line flare and allows estimating the numbers of and energy contents in electrons, protons and alpha particles produced during the flare. The energy spectrum of the electrons producing the HXR/GR continuum is consistent with a broken power-law with an energy break at ~800 keV. It is shown that the high-energy part (above ~800 keV) of this distribution is responsible for the high-frequency radio emission (> 20 GHz) detected during the flare. By comparing the 30 THz emission expected from semi-empirical and time-independent models of the quiet and flare atmospheres, we find that most (~80%) of the observed 30 THz radiation can be attributed to thermal free-free emission of an optically-thin source. Using the F2 flare atmospheric model this thin source is found to be at temperatures T~8000 K and is located well above the minimum temperature region. We argue that the chromospheric heating, which results in 80% of the 30 THz excess radiation, can be due to energy deposition by non-thermal flare accelerated electrons, protons and alpha particles. The remaining 20% of the 30 THz excess emission is found to be radiated from an optically-thick atmospheric layer at T~5000 K, below the temperature minimum region, where direct heating by non-thermal particles is insufficient to account for the observed infrared radiation.

Recent solar flare observations in the sub-THz range have provided evidence of a new spectral component with fluxes increasing for larger frequencies, separated from the well-known microwave emission that maximizes in the GHz range. Suggested interpretations explain the THz spectral component, but do not account for the simultaneous microwave component. We present a mechanism for producing the observed double-spectra. Based on coherent enhancement of synchrotron emission at long wavelengths in laboratory accelerators, we consider how similar processes may occur within a solar flare. The instability known as microbunching arises from perturbations that produce electron beam density modulations, giving rise to broadband coherent synchrotron emission at wavelengths comparable to the characteristic size of the microbunch structure. The spectral intensity of this coherent synchrotron radiation (CSR) can far exceed that of the incoherent synchrotron radiation (ISR), which peaks at higher frequency, thus producing a double-peaked spectrum. Successful CSR simulations are shown to fit actual burst spectral observations, using typical flaring physical parameters and power-law energy distributions for the accelerated electrons. The simulations consider an energy threshold below which microbunching is not possible because of Coulomb repulsion. Only a small fraction of the radiating charges accelerated to energies above the threshold is required to produce the microwave component observed for several events. The ISR-CSR mechanism can occur together with other emission processes producing the microwave component. It may bring an important contribution at microwaves at least for certain events where physical conditions for the occurrence of the ISR-CSR microbunching mechanism are possible.

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 present a detailed theoretical analysis of the local radio emission at the lower part of the solar atmosphere. To accomplish this, we have used a numerical code to simulate the emission and transport of high frequency electromagnetic waves from 2 GHz up to 10 THz. As initial conditions we used three well know chromospheric models. In this way, the generated synthetic spectra allows us to study the local emission and absorption processes with high resolution in both altitude and frequency. Associated with the temperature minimum predicted by these models we found that the local optical depth at millimeter wavelengths remains constant, producing an optically thin layer which is surrounded by two layers of high local emission. We call this structure the Chromospheric Solar Millimeter-wave Cavity (CSMC). The CSMC shows the complexity of the relationship between the theoretical temperature profile and the observed brightness temperature and may help to understand the dispersion of the observed brightness temperature in the millimeter wavelength range.

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.

This review discusses the complementary relationship between radio and hard X-ray observations of the Sun using primarily results from the era of the Reuven Ramaty High Energy Solar Spectroscopic Imager satellite. A primary focus of joint radio and hard X-ray studies of solar flares uses observations of nonthermal gyrosynchrotron emission at radio wavelengths and bremsstrahlung hard X-rays to study the properties of electrons accelerated in the main flare site, since it is well established that these two emissions show very similar temporal behavior. A quantitative prescription is given for comparing the electron energy distributions derived separately from the two wavelength ranges: this is an important application with the potential for measuring the magnetic field strength in the flaring region, and reveals significant differences between the electrons in different energy ranges. Examples of the use of simultaneous data from the two wavelength ranges to derive physical conditions are then discussed, including the case of microflares, and the comparison of images at radio and hard X-ray wavelengths is presented. There have been puzzling results obtained from observations of solar flares at millimeter and submillimeter wavelengths, and the comparison of these results with corresponding hard X-ray data is presented. Finally, the review discusses the association of hard X-ray releases with radio emission at decimeter and meter wavelengths, which is dominated by plasma emission (at lower frequencies) and electron cyclotron maser emission (at higher frequencies), both coherent emission mechanisms that require small numbers of energetic electrons. These comparisons show broad general associations but detailed correspondence remains more elusive.

Solar flares observed in the 200-400 GHz radio domain may exhibit a slowly varying and time-extended component which follows a short (few minutes) impulsive phase and which lasts for a few tens of minutes to more than one hour. The few examples discussed in the literature indicate that such long-lasting submillimeter emission is most likely thermal bremsstrahlung. We present a detailed analysis of the time-extended phase of the 2003 October 27 (M6.7) flare, combining 1-345 GHz total-flux radio measurements with X-ray, EUV, and H{\alpha} observations. We find that the time-extended radio emission is, as expected, radiated by thermal bremsstrahlung. Up to 230 GHz, it is entirely produced in the corona by hot and cool materials at 7-16 MK and 1-3 MK, respectively. At 345 GHz, there is an additional contribution from chromospheric material at a few 10^4 K. These results, which may also apply to other millimeter-submillimeter radio events, are not consistent with the expectations from standard semi-empirical models of the chromosphere and transition region during flares, which predict observable radio emission from the chromosphere at all frequencies where the corona is transparent.

We use a numerical code called PAKALMPI to compute synthetic spectra of the solar emission in quiet conditions at millimeter, sub-millimeter and infrared wavelengths. PAKALMPI solves the radiative transfer equation, with Non Local Thermodynamic Equilibrium (NLTE), in a three dimensional geometry using a multiprocessor environment. The code is able to use three opacity functions: classical bremsstrahlung, H-and inverse bremsstrahlung. In this work we have computed and compared two synthetic spectra, one in the common way: using bremsstrahlung opacity function and considering a fully ionized atmosphere; and a new one considering bremsstrahlung, inverse bremsstrahlung and H- opacity functions in NLTE. We analyzed in detail the local behavior of the low atmospheric emission at 17, 212, and 405 GHz (frequencies used by the Nobeyama Radio Heliograph and the Solar Submillimeter Telescope). We found that the H- is the major emission mechanism at low altitudes (below 500 km) and that at higher altitudes the classical bremsstrahlung becomes the major mechanism of emission. However the brightness temperature remains unalterable. Finally, we found that the inverse bremsstrahlung process is not important for the radio emission at these heights.

We report on the first detection of ionospheric disturbances caused by short repeated gamma-ray bursts from the magnetar SGR J1550-5418. Very low frequency (VLF) radio wave data obtained in South America clearly show sudden amplitude and phase changes at the corresponding times of eight SGR bursts. Maximum amplitude and phase changes of the VLF signals appear to be correlated with the gamma-ray fluence. On the other hand, VLF recovery timescales do not show any significant correlation with the fluence, possibly suggesting that the bursts' spectra are not similar to each other. In summary, the Earth's ionosphere can be used as a very large gamma-ray detector and the VLF observations provide us with a new method to monitor high energy astrophysical phenomena without interruption such as Earth Occultation.

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.