We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of a quiet-Sun region at a wavelength of 3 mm, obtained during the first solar ALMA cycle on April 27, 2017, and compare them with available chromospheric observations in the UV and visible as well as with photospheric magnetograms. ALMA images clearly reveal the presence of distinct particularly dark/cool areas in the millimeter maps having temperatures as low as 60% of the normal quiet Sun at 3 mm, which are not seen in the other data. We speculate that ALMA is sensing cool chromospheric gas, whose presence had earlier been inferred from infrared CO spectra.

3D picture of the coronal magnetic field remains an outstanding problem in solar physics, particularly, in active regions. Nonlinear force-free field reconstructions that employ routinely available full-disk photospheric vector magnetograms represent state-of-the-art coronal magnetic field modeling. Such reconstructions, however, suffer from an inconsistency between a force-free coronal magnetic field and non-force-free photospheric boundary condition, from which the coronal reconstruction is performed. In this study we focus on integrating the additional chromospheric and / or coronal magnetic field data with the vector photospheric magnetograms with the goal of improving the reliability of the magnetic field reconstructions. We develop a corresponding modification of the available optimization codes described in Fleishman et al. (2017) and test their performance using a full-fledged MHD model obtained from the Bifrost code by performing a `voxel-by-voxel' comparison between the reconstructed and the model magnetic fields. We demonstrate that adding even an incomplete set of chromospheric magnetic field data can measurably improve the reconstruction of the coronal magnetic field, greatly improve reconstructions of the magnetic connectivity and of the coronal electric current.

Quantifying coronal magnetic field remains a central problem in solar physics. Nowadays the coronal magnetic field is often modelled using nonlinear force-free field (NLFFF) reconstructions, whose accuracy has not yet been comprehensively assessed. Here we perform a detailed casting of the NLFFF reconstruction tools, such as pi-disambiguation, photospheric field preprocessing, and volume reconstruction methods using a 3D snapshot of the publicly available full-fledged radiative MHD model. Specifically, from the MHD model we know the magnetic field vector in the entire 3D domain, which enables us to perform "voxel-by-voxel" comparison of the restored and the true magnetic field in the 3D model volume. Our tests show that the available pi-disambiguation methods often fail at the quiet sun areas dominated by small-scale magnetic elements, while they work well at the AR photosphere and (even better) chromosphere. The preprocessing of the photospheric magnetic field, although does produce a more force-free boundary condition, also results in some effective `elevation' of the magnetic field components. This `elevation' height is different for the longitudinal and transverse components, which results in a systematic error in absolute heights in the reconstructed magnetic data cube. The extrapolations performed starting from actual AR photospheric magnetogram are free from this systematic error, while have other metrics comparable with those for extrapolations from the preprocessed magnetograms. This finding favors the use of extrapolations from the original photospheric magnetogram without preprocessing. Our tests further suggest that extrapolations from a force-free chromospheric boundary produce measurably better results, than those from the photospheric boundary.

The spectrum of gyrosynchrotron emission from solar flares generally peaks in the microwave range. Its optically-thin, high-frequency component, above the spectral peak, is often used for diagnostics of the nonthermal electrons and the magnetic field in the radio source. Under favorable conditions, its low-frequency counterpart brings additional, complementary information about these parameters as well as thermal plasma diagnostics, either through gyrosynchrotron self-absorption, free-free absorption by the thermal plasma, or the suppression of emission through the so-called Razin effect. However, their effects on the low-frequency spectrum are often masked by spatial nonuniformity. To disentangle the various contributions to low-frequency gyrosynchrotron emission, a combination of spectral and imaging data is needed. To this end, we have investigated Owens Valley Solar Array (OVSA) multi-frequency images for 26 solar bursts observed jointly with Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) during the first half of 2002. For each, we examined dynamic spectra, time- and frequency-synthesis maps, RHESSI images with overlaid OVSA contours, and a few representative single-frequency snapshot OVSA images. We focus on the frequency dependence of microwave source sizes derived from the OVSA images and their effect on the low-frequency microwave spectral slope. We succeed in categorizing 18 analyzed events into several groups. Four events demonstrate clear evidence of being dominated by gyrosynchrotron self-absorption, with an inferred brightness temperature of $\geq10^8$~K. The low-frequency spectra in the remaining events are affected to varying degree by Razin suppression. We find that many radio sources are rather large at low frequencies, which can have important implications for solar energetic particle production and escape.

The Atacama Large Millimeter-Submillimeter Array (ALMA) has opened a new window for studying the Sun via high-resolution high-sensitivity imaging at millimeter wavelengths. In this contribution I review the capabilities of the instrument for solar observing and describe the extensive effort taken to bring the possibility of solar observing with ALMA to the scientific community. The first solar ALMA observations were carried out during 2014 and 2015 in two ALMA bands, Band 3 (3 mm) and Band 6 (1.3 mm), in single-dish and interferometric modes, using single pointing and mosaicing observing techniques, with spatial resolution up to 2arcsec and 1arcsec in the two bands, respectively. I overview several recently published studies which made use of the first solar ALMA observations, describe current status of solar observing with ALMA and briefly discuss the future capabilities of the instrument.

We present the first high-resolution Atacama Large Millimeter/Submillimeter Array (ALMA) observations of a sunspot at wavelengths of 1.3 mm and 3 mm, obtained during the solar ALMA Science Verification campaign in 2015, and compare them with the predictions of semi-empirical sunspot umbral/penumbral atmosphere models. For the first time millimeter observations of sunspots have resolved umbral/penumbral brightness structure at the chromospheric heights, where the emission at these wavelengths is formed. We find that the sunspot umbra exhibits a radically different appearance at 1.3 mm and 3 mm, whereas the penumbral brightness structure is similar at the two wavelengths. The inner part of the umbra is ~600 K brighter than the surrounding quiet Sun (QS) at 3 mm and is ~700 K cooler than the QS at 1.3 mm, being the coolest part of sunspot at this wavelength. On average, the brightness of the penumbra at 3 mm is comparable to the QS brightness, while at 1.3 mm it is ~1000 K brighter than the QS. Penumbral brightness increases towards the outer boundary in both ALMA bands. Among the tested umbral models, that of Severino et al. (1994) provides the best fit to the observational data, including both the ALMA data analyzed in this study and data from earlier works. No penumbral model amongst those considered here gives a satisfactory fit to the currently available measurements. ALMA observations at multiple mm wavelengths can be used for testing existing sunspot models, and serve as an important input to constrain new empirical models.

We present joint observations of the Sun by the Atacama Large Millimeter/submillimeter Array (ALMA) and the Interface Region Imaging Spectrograph (IRIS). The observations were made of a solar active region on 2015 December 18 as part of the ALMA science verification effort. A map of the Sun's continuum emission of size $2.4' \times 2.3'$ was obtained by ALMA at a wavelength of 1.25 mm (239 GHz) using mosaicing techniques. A contemporaneous map of size $1.9'\times 2.9'$ was obtained in the Mg II h doublet line at 2803.5\AA\ by IRIS. Both mm/submm$-\lambda$ continuum emission and ultraviolet (UV) line emission are believed to originate from the solar chromosphere and both have the potential to serve as powerful and complementary diagnostics of physical conditions in this poorly understood layer of the solar atmosphere. While a clear correlation between mm-$\lambda$ brightness temperature $T_B$ and the Mg II h line radiation temperature $T_{rad}$ is observed the slope is $<1$, perhaps as a result of the fact that these diagnostics are sensitive to different parts of the chromosphere and/or the Mg II h line source function includes a scattering component. There is a significant offset between the mean $T_B$(1.25 mm) and mean $T_{rad}$(Mg II), the former being $\approx 35\%$ greater than the latter. Partitioning the maps into "sunspot", "quiet regions", and "plage regions" we find that the slope of the scatter plots between the IRIS Mg II h line $T_{rad}$ and the ALMA brightness temperature $T_B$ is 0.4 (sunspot), 0.56 (quiet regions), and 0.66 (plage regions). We suggest that this change may be caused by the regional dependence of the formation heights of the IRIS and ALMA diagnostics, and/or the increased degree of coupling between the UV source function and the local gas temperature in the hotter, denser gas in plage regions.

The Atacama Large Millimeter-submillimeter Array (ALMA) radio telescope has commenced science observations of the Sun starting in late 2016. Since the Sun is much larger than the field of view of individual ALMA dishes, the ALMA interferometer is unable to measure the background level of solar emission when observing the solar disk. The absolute temperature scale is a critical measurement for much of ALMA solar science, including the understanding of energy transfer through the solar atmosphere, the properties of prominences, and the study of shock heating in the chromosphere. In order to provide an absolute temperature scale, ALMA solar observing will take advantage of the remarkable fast-scanning capabilities of the ALMA 12m dishes to make single-dish maps of the full Sun. This article reports on the results of an extensive commissioning effort to optimize the mapping procedure, and it describes the nature of the resulting data. Amplitude calibration is discussed in detail: a path that utilizes the two loads in the ALMA calibration system as well as sky measurements is described and applied to commissioning data. Inspection of a large number of single-dish datasets shows significant variation in the resulting temperatures, and based on the temperature distributions we derive quiet-Sun values at disk center of 7300 K at lambda=3 mm and 5900 K at lambda=1.3 mm. These values have statistical uncertainties of order 100 K, but systematic uncertainties in the temperature scale that may be significantly larger. Example images are presented from two periods with very different levels of solar activity. At a resolution of order 25 arcsec, the 1.3 mm wavelength images show temperatures on the disk that vary over about a 2000 K range.

We report the discovery of a brightness enhancement in the center of a large sunspot umbra at a wavelength of 3 mm using the Atacama Large Millimeter/sub-millimeter Array (ALMA). Sunspots are amongst the most prominent features on the solar surface, but many of their aspects are surprisingly poorly understood. We analyzed a {\lambda}=3 mm (100 GHz) mosaic image obtained by ALMA, which includes a large sunspot within the active region AR12470 on December 16, 2015. The 3 mm map has a field-of-view and spatial resolution, which is the highest spatial-resolution map of an entire sunspot in this frequency range. We find a gradient of 3 mm brightness from a high value in the outer penumbra to a low value in the inner penumbra/outer umbra. Within the inner umbra, there is a marked increase in 3mm brightness temperature, which we call an umbral brightness enhancement. This enhanced emission corresponds to a temperature excess of 800 K relative to the surrounding inner penumbral region and coincides with excess brightness in the 1330 and 1400 {\AA} slitjaw images of the Interface Region Imaging Spectrograph (IRIS), adjacent to a partial lightbridge. This {\lambda}=3 mm brightness enhancement may be an intrinsic feature of the sunspot umbra at chromospheric heights, such as a manifestation of umbral flashes, or it could be related to a coronal plume since the brightness enhancement was coincident with the footpoint of a coronal loop observed at 171 {\AA}.

Observations of the Sun at millimeter and submillimeter wavelengths offer a unique probe into the structure, dynamics, and heating of the chromosphere; the structure of sunspots; the formation and eruption of prominences and filaments; and energetic phenomena such as jets and flares. High-resolution observations of the Sun at millimeter and submillimeter wavelengths are challenging due to the intense, extended, low- contrast, and dynamic nature of emission from the quiet Sun, and the extremely intense and variable nature of emissions associated with energetic phenomena. The Atacama Large Millimeter/submillimeter Array (ALMA) was designed with solar observations in mind. The requirements for solar observations are significantly different from observations of sidereal sources and special measures are necessary to successfully carry out this type of observations. We describe the commissioning efforts that enable the use of two frequency bands, the 3 mm band (Band 3) and the 1.25 mm band (Band 6), for continuum interferometric-imaging observations of the Sun with ALMA. Examples of high-resolution synthesized images obtained using the newly commissioned modes during the solar commissioning campaign held in December 2015 are presented. Although only 30 of the eventual 66 ALMA antennas were used for the campaign, the solar images synthesized from the ALMA commissioning data reveal new features of the solar atmosphere that demonstrate the potential power of ALMA solar observations. The ongoing expansion of ALMA and solar-commissioning efforts will continue to enable new and unique solar observing capabilities.

We use state-of-the-art, three-dimensional non-local thermodynamic equilibrium (non-LTE) radiative magnetohydrodynamic simulations of the quiet solar atmosphere to carry out detailed tests of chromospheric magnetic field diagnostics from free-free radiation at millimeter and submillimeter wavelengths (mm/submm). The vertical component of the magnetic field was deduced from the mm/submm brightness spectra and the degree of circular polarization synthesized at millimeter frequencies. We used the frequency bands observed by the Atacama Large Millimeter/Submillimeter Array (ALMA) as a convenient reference. The magnetic field maps obtained describe the longitudinal magnetic field at the effective formation heights of the relevant wavelengths in the solar chromosphere. The comparison of the deduced and model chromospheric magnetic fields at the spatial resolution of both the model and current observations demonstrates a good correlation, but has a tendency to underestimate the model field. The systematic discrepancy of about 10 percent is probably due to averaging of the restored field over the heights contributing to the radiation, weighted by the strength of the contribution. On the whole, the method of probing the longitudinal component of the magnetic field with free-free emission at mm/submm wavelengths is found to be applicable to measurements of the weak quiet-Sun magnetic fields. However, successful exploitation of this technique requires very accurate measurements of the polarization properties (primary beam and receiver polarization response) of the antennas, which will be the principal factor that determines the level to which chromospheric magnetic fields can be measured. Consequently, high-resolution and high-precision observations of circularly polarized radiation at millimeter wavelengths can be a powerful tool for producing chromospheric longitudinal magnetograms.

This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interesting observations could be made already with the still limited solar observing modes foreseen for Cycle 4 and that ALMA has the potential to make important contributions to answering long-standing scientific questions in solar physics. With the proposal deadline for ALMA Cycle 4 in April 2016 and the Commissioning and Science Verification campaign in December 2015 in sight, several of the SSALMON Expert Teams composed strategic documents in which they outlined potential solar observations that could be feasible given the anticipated technical capabilities in Cycle 4. These documents have been combined and supplemented with an analysis, resulting in recommendations for solar observing with ALMA in Cycle 4. In addition, the detailed science cases also demonstrate the scientific priorities of the solar physics community and which capabilities are wanted for the next observing cycles. The work on this White Paper effort was coordinated in close cooperation with the two international solar ALMA development studies led by T. Bastian (NRAO, USA) and R. Brajsa, (ESO). This document will be further updated until the beginning of Cycle 4 in October 2016. In particular, we plan to adjust the technical capabilities of the solar observing modes once finally decided and to further demonstrate the feasibility and scientific potential of the included science cases by means of numerical simulations of the solar atmosphere and corresponding simulated ALMA observations.

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

ALMA will open a new chapter in the study of the Sun by providing a leap in spatial resolution and sensitivity compared to currently available mm wave- length observations. In preparation of ALMA, we have carried out a large number of observational tests and state-of-the-art radiation MHD simulations. Here we review the best available observations of the Sun at millimeter wavelengths.Using state of the art radiation MHD simulations of the solar atmosphere we demonstrate the huge potential of ALMA observations for uncovering the nature of the solar chromosphere. We show that ALMA will not only provide a reliable probe of the thermal structure and dynamics of the chromosphere, it will also open up a powerful new diagnostic of magnetic field at chromospheric heights, a fundamentally important, but so far poorly known parameter.

We developed an efficient algorithm integrated in our 3D modeling tool, GX Simulator (Nita et al. 2015), allowing quick computation of the synthetic intensity and polarization maps of solar active regions (AR) in the ALMA spectral range. The algorithm analyzes the photospheric input (white light and magnetogram) to classify a given photospheric pixel to belong to a given photospheric structure. Then, a 1D chromospheric model (Fontenla et al. 2009) is added on top of each pixel, which forms a chromospheric model of the AR. Next step is computation of the mm and sub-mm emission produced from this chromosphere model. A huge advantage of this approach is that emission from any given AR can be synthesized very fast, on the order of a few minutes after the AR selection. Using the GX Simulator tool it is also possible to produce synthetic maps of the microwave (gyroresonance) and EUV emission from the same AR model and compare them with the ALMA synthetic maps and with the corresponding observed microwave and/or EUV data.

The Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) was initiated in 2014 in connection with two ALMA development studies. The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful new tool, which can also observe the Sun at high spatial, temporal, and spectral resolution. The international SSALMONetwork aims at coordinating the further development of solar observing modes for ALMA and at promoting scientific opportunities for solar physics with particular focus on numerical simulations, which can provide important constraints for the observing modes and can aid the interpretation of future observations. The radiation detected by ALMA originates mostly in the solar chromosphere - a complex and dynamic layer between the photosphere and corona, which plays an important role in the transport of energy and matter and the heating of the outer layers of the solar atmosphere. Potential targets include active regions, prominences, quiet Sun regions, flares. Here, we give a brief overview over the network and potential science cases for future solar observations with ALMA.

Aims. We use advanced 3D NLTE radiative magnetohydrodynamic simulations of the solar atmosphere to carry out detailed tests of chromospheric diagnostics at millimeter and submillimeter wavelengths. Methods. We focused on the diagnostics of the thermal structure of the chromosphere in the wavelength bands from 0.4 mm up to 9.6 mm that can be accessed with the Atacama Large Millimeter/Submillimeter Array (ALMA) and investigated how these diagnostics are affected by the instrumental resolution. Results. We find that the formation height range of the millimeter radiation depends on the location in the simulation domain and is related to the underlying magnetic structure. Nonetheless, the brightness temperature is a reasonable measure of the gas temperature at the effective formation height at a given location on the solar surface. There is considerable scatter in this relationship, but this is significantly reduced when very weak magnetic fields are avoided. Our results indicate that although instrumental smearing reduces the correlation between brightness and temperature, millimeter brightness can still be used to reliably diagnose electron temperature up to a resolution of 1". If the resolution is more degraded, then the value of the diagnostic diminishes rapidly. Conclusions. We conclude that millimeter brightness can image the chromospheric thermal structure at the height at which the radiation is formed. Thus multiwavelength observations with ALMA with a narrow step in wavelength should provide sufficient information for a tomographic imaging of the chromosphere.

Aims: The aim of this paper is to demonstrate that millimeter wave data can be used to distinguish between various atmospheric models of sunspots, whose temperature structure in the upper photosphere and chromosphere has been the source of some controversy. Methods: We use observations of the temperature contrast (relative to the quiet Sun) above a sunspot umbra at 3.5 mm obtained with the Berkeley-Illinois-Maryland Array (BIMA), complemented by submm observations from Lindsey & Kopp (1995) and 2 cm observations with the Very Large Array. These are compared with the umbral contrast calculated from various atmospheric models of sunspots. Results: Current mm and submm observational data suggest that the brightness observed at these wavelengths is low compared to the most widely used sunspot models. These data impose strong constraints on the temperature and density stratifications of the sunspot umbral atmosphere, in particular on the location and depth of the temperature minimum and the location of the transition region. Conclusions: A successful model that is in agreement with millimeter umbral brightness should have an extended and deep temperature minimum (below 3000 K). Better spatial resolution as well as better wavelength coverage are needed for a more complete determination of the chromospheric temperature stratification above sunspot umbrae.

Aims. We analyze observational data from 4 instruments to study the correlations between chromospheric emission, spanning the heights from the temperature minimum region to the middle chromosphere, and photospheric magnetic field. Methods: The data consist of radio images at 3.5 mm from the Berkeley-Illinois-Maryland Array (BIMA), UV images at 1600 A from TRACE, Ca II K-line filtergrams from BBSO, and MDI/SOHO longitudinal photospheric magnetograms. For the first time interferometric millimeter data with the highest currently available resolution are included in such an analysis. We determine various parameters of the intensity maps and correlate the intensities with each other and with the magnetic field. Results: The chromospheric diagnostics studied here show a pronounced similarity in their brightness structures and map out the underlying photospheric magnetic field relatively well. We find a power law to be a good representation of the relationship between photospheric magnetic field and emission from chromospheric diagnostics at all wavelengths. The dependence of chromospheric brightness on magnetic field is found to be different for network and internetwork regions.

The very nature of the solar chromosphere, its structuring and dynamics, remains far from being properly understood, in spite of intensive research. Here we point out the potential of chromospheric observations at millimeter wavelengths to resolve this long-standing problem. Computations carried out with a sophisticated dynamic model of the solar chromosphere due to Carlsson and Stein demonstrate that millimeter emission is extremely sensitive to dynamic processes in the chromosphere and the appropriate wavelengths to look for dynamic signatures are in the range 0.8-5.0 mm. The model also suggests that high resolution observations at mm wavelengths, as will be provided by ALMA, will have the unique property of reacting to both the hot and the cool gas, and thus will have the potential of distinguishing between rival models of the solar atmosphere. Thus, initial results obtained from the observations of the quiet Sun at 3.5 mm with the BIMA array (resolution of 12 arcsec) reveal significant oscillations with amplitudes of 50-150 K and frequencies of 1.5-8 mHz with a tendency toward short-period oscillations in internetwork and longer periods in network regions. However higher spatial resolution, such as that provided by ALMA, is required for a clean separation between the features within the solar atmosphere and for an adequate comparison with the output of the comprehensive dynamic simulations.