The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

ALMA observations of the Sun at mm-$\lambda$ offer a unique opportunity to investigate the temperature structure of the solar chromosphere. In this article we expand our previous work on modeling the chromospheric temperature of the quiet Sun, by including measurements of the brightness temperature in the network and cell interiors, from high resolution ALMA images at 3 mm (Band 3) and 1.26 mm (Band 6). We also examine the absolute calibration of ALMA full-disk images. We suggest that the brightness temperature at the center of the solar disk in Band 6 is $\sim440$ K above the value recommended by White et al. (2017) and we give improved results for the electron temperature variation of the average quiet Sun with optical depth, as well as the derived spectrum at the center of the disk. We found that the electron temperature in the network is considerably lower than predicted by model F of Fontenla et al. (1993) and that of the cell interior considerably higher than predicted by model A. Depending upon the network/cell segregation scheme, the electron temperature difference between network and cell at $\tau=1$ (100 GHz) is from $\sim$660 to $\sim$1550 K, compared to $\sim$3280 K predicted by the models; similarly, the $T_e$ ratio is from $\sim$1.10, to 1.24, against $\sim$1.55 of the model prediction. We also found that the network/cell $T_e(\tau)$ curves diverge as $\tau$ decreases, indicating an increase of contrast with height and possibly a steeper temperature rise in the network than in the cell interior.

We present observations of the solar chromosphere obtained simultaneously with the Atacama Large Millimeter/submillimeter Array (\emph{ALMA}) and the Interface Region Imaging Spectrograph (\emph{IRIS}). The observatories targeted a chromospheric plage region of which the spatial distribution (split between strongly and weakly magnetized regions) allowed the study of linear-like structures in isolation, free of contamination from background emission. Using these observations in conjunction with a radiative magnetohydrodynamic 2.5D model covering the upper convection zone all the way to the corona that considers non-equilibrium ionization effects, we report the detection of an on-disk chromospheric spicule with \emph{ALMA} and confirm its multithermal nature. In addition, we discuss the strikingly high degree of similarity between chromospheric plage features observed in \emph{ALMA}/Band6 and \emph{IRIS}/\ion{Si}{4} (also reproduced in our model) suggesting that \emph{ALMA}/Band6 does not observe in the low chromosphere as previously thought but rather observes the upper chromospheric parts of structures such as spicules and other bright structures above plage at geometric heights near transition region temperatures. We also show that \emph{IRIS}/\ion{Mg}{2} is not as well correlated with \emph{ALMA}/Band6 as was previously thought. For these comparisons, we propose and employ a novel empirical method for the determination of plage regions, which seems to better isolate plage from its surrounding regions as compared to other methods commonly used. We caution that isolating plage from its immediate surroundings must be done with care to mitigate stastistical bias in quantitative comparisons between different chromospheric observables. Lastly, we report indications for chromospheric heating due to traveling shocks supported by the \emph{ALMA}/Band6 observations.

Using ALMA observations, we performed the first systematic survey for transient brightenings (i.e. weak, small-scale episodes of energy release) in the quiet solar chromosphere at 3 mm. Our dataset included images of six 87'' x 87'' regions of the quiet Sun obtained with angular resolution of a few arcsec at a cadence of 2 s. The transient brightenings were detected as weak enhancements above the average intensity after we removed the effect of the p-mode oscillations. A similar analysis, over the same regions, was performed for simultaneous 304 and 1600 \AA\ data obtained with the Atmospheric Imaging Assembly. We detected 184 3 mm transient brightening events with brightness temperatures from 70 K to more than 500 K above backgrounds of $\sim 7200-7450$ K. Their mean duration and maximum area were 51.1 s and 12.3 Mm$^2$, respectively, with a weak preference of appearing at network boundaries rather than in cell interiors. Both parameters exhibited power-law behavior with indices of 2.35 and 2.71, respectively. Only a small fraction of ALMA events had either 304 or 1600 \AA\ counterparts but the properties of these events were not significantly different from those of the general population except that they lacked their low-end energy values. The total thermal energies of the ALMA transient brightenings were between $1.5 \times 10^{24}$ and $9.9 \times 10^{25}$ erg and their frequency distribution versus energy was a power law with an index of 1.67. We found that the power per unit area provided by the ALMA events could account for only 1\% of the chromospheric radiative losses (10\% of the coronal ones). Therefore, their energy budget falls short of meeting the requirements for the heating of the upper layers of the solar atmosphere and this conclusion does not change even if we use the least restrictive criteria possible for the detection of transient brightenings.

We present an initial study of one of the first ALMA Band 3 observations of the Sun with the aim to characterise the diagnostic potential of brightness temperatures measured with ALMA on the Sun. The observation covers 48min at a cadence of 2s targeting a Quiet Sun region at disk-centre. Corresponding time series of brightness temperature maps are constructed with the first version of the Solar ALMA Pipeline (SoAP) and compared to simultaneous SDO observations. The angular resolution of the observations is set by the synthesized beam (1.4x2.1as). The ALMA maps exhibit network patches, internetwork regions and also elongated thin features that are connected to large-scale magnetic loops as confirmed by a comparison with SDO maps. The ALMA Band 3 maps correlate best with the SDO/AIA 171, 131 and 304 channels in that they exhibit network features and, although very weak in the ALMA maps, imprints of large-scale loops. A group of compact magnetic loops is very clearly visible in ALMA Band 3. The brightness temperatures in the loop tops reach values of about 8000-9000K and in extreme moments up to 10 000K. ALMA Band 3 interferometric observations from early observing cycles already reveal temperature differences in the solar chromosphere. The weak imprint of magnetic loops and the correlation with the 171, 131, and 304 SDO channels suggests though that the radiation mapped in ALMA Band 3 might have contributions from a larger range of atmospheric heights than previously assumed but the exact formation height of Band 3 needs to be investigated in more detail. The absolute brightness temperature scale as set by Total Power measurements remains less certain and must be improved in the future. Despite these complications and the limited angular resolution, ALMA Band 3 observations have large potential for quantitative studies of the small-scale structure and dynamics of the solar chromosphere.

We studied chromospheric oscillations using Atacama Large millimeter and sub-millimeter Array (ALMA) time-series of interferometric observations of the quiet Sun obtained at 3 mm with a 2-s cadence and a spatial resolution of a few arcsec. The same analysis, over the same fields of view and for the same intervals, was performed for simultaneous Atmospheric Imaging Assembly (AIA) image sequences in 1600 A. Spatially-resolved chromospheric oscillations at 3 mm, with frequencies of $ 4.2 +- 1.7$ mHz are observed in the quiet Sun, in both cell and network. The coherence length-scale of the oscillations is commensurate with the spatial resolution of our ALMA observations. Brightness-temperature fluctuations in individual pixels could reach up to a few hundred K, while the spatially averaged power spectral densities yield rms in the range ~ 55-75 K, i.e., up to ~ 1 % of the averaged brightness temperatures and exhibit a moderate increase towards the limb. For AIA 1600 A, the oscillation frequency is 3.7 +- 1.7 mHz. The relative rms is up to 6 % of the background intensity, with a weak increase towards disk center (cell, average). ALMA 3 mm time-series lag AIA 1600 A by ~ 100 s, which corresponds to a formation-height difference of ~ 1200 km. The ALMA oscillations that we detected exhibit higher amplitudes than those derived from the lower (~ 10 arcsec) resolution observations at 3.5 mm by White et al. (2006). Chromospheric oscillations are, therefore, not fully resolved at the length-scale of the chromospheric network, and possibly not even at the spatial resolution of our ALMA observations. Any study of transient brightenings in the mm-domain should take into account the oscillations.

This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

Remarkable progress has been made in understanding turbulent astrophysical plasmas in past decades including, notably, the solar wind and the interstellar medium. In the case of the solar wind, much of this progress has relied on in situ measurements from space-borne instruments. However, ground-based radio observations also have played a significant role and have the potential to play an even bigger role. In particular, using distant background sources (quasars, pulsars, satellite beacons) to transilluminate the foreground corona and solar wind, a variety of radio propagation phenomena can be used to map plasma properties of the solar corona and heliosphere, as well as the warm interstellar medium. These include angular broadening, interplanetary and interstellar scintillations, and differential Faraday rotation. These observations are highly complementary to in situ observations of the solar wind, and could be a mainstay of investigations into turbulence of the ISM. We point out that the Next Generation Very Large Array (ngVLA) fulfills all the requirements necessary to exploit radio observations of astrophysical turbulence fully.

The Sun is an active star that can have a direct impact on the Earth, its magnetosphere, and the technological infrastructure on which modern society depends. Among the phenomena that drive "space weather" are fast solar wind streams and co-rotating interaction regions, solar flares, coronal mass ejections, the shocks they produce, and the energetic particles they accelerate. Radio emission from these and associated phenomena offer unique diagnostic possibilities that complement those available at other wavelengths. Here, the relevant space weather drivers are briefly described, the potential role of radio observations is outlined, and the requirements of an instrument to provide them are provided: specifically, ultrabroadband imaging spectropolarimetry. The insights provided by radio observations of space weather drivers will not only inform the science of space weather, they will pave the way for new tools for forecasting and "nowcasting" space weather. They will also serve as an important touchstone against which local environment of exoplanets and the impact of "exo-space weather" can be evaluated.

Identification of the mechanisms responsible for heating the solar chromosphere and corona remains an outstanding problem, one of great relevance to late-type stars as well. There has been tremendous progress in the past decade, largely driven by new instruments, new observations, and sophisticated modeling efforts. Despite this progress, gaps remain. We briefly discuss the need for radio coverage of the 3D solar atmosphere and discuss the requirements.

Magnetic reconnection is a fundamental physical process in many laboratory, space, and astrophysical plasma contexts. Solar flares serve as an outstanding laboratory to study the magnetic reconnection and the associated energy release and conversion processes under plasma conditions difficult to reproduce in the laboratory, and with considerable spatiotemporal details not possible elsewhere in astrophysics. Here we emphasize the unique power of remote-sensing observations of solar flares at radio wavelengths. In particular, we discuss the transformative technique of broadband radio dynamic imaging spectroscopy in making significant contributions to addressing several outstanding challenges in magnetic reconnection, including the capability of pinpointing magnetic reconnection sites, measuring the time-evolving reconnecting magnetic fields, and deriving the spatially and temporally resolved distribution function of flare-accelerated electrons.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

We present Very Large Array observations at 33.0 GHz that detect emission coincident with $\epsilon$ Eridani to within $0\rlap.{"}07$ (0.2 AU at the distance of this star), with a positional accuracy of $0\rlap.{"}05$. This result strongly supports the suggestion of previous authors that the quiescent centimeter emission comes from the star and not from a proposed giant exoplanet with a semi-major axis of $\sim1\rlap.{"}0$ (3.4 AU). The centimeter emission is remarkably flat and is consistent with optically thin free-free emission. In particular, it can be modeled as a stellar wind with a mass loss rate of the order of $6.6 \times 10^{-11}~ M_\odot ~yr^{-1}$, which is 3,300 times the solar value, exceeding other estimates of this star's wind. However, interpretation of the emission in terms of other thermal mechanisms like coronal free-free and gyroresonance emission cannot be discarded.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

We present an overview of high resolution quiet Sun observations, from disk center to the limb, obtained with the Atacama Large mm and sub-mm Array (ALMA) at 3 mm. Seven quiet Sun regions were observed with resolution of up to 2.5" by 4.5". We produced both average and snapshot images by self-calibrating the ALMA visibilities and combining the interferometric images with full disk solar images. The images show well the chromospheric network, which, based on the unique segregation method we used, is brighter than the average over the fields of view of the observed regions by $\sim 305$ K while the intranetwork is less bright by $\sim 280$ K, with a slight decrease of the network/intranetwork contrast toward the limb. At 3 mm the network is very similar to the 1600 \AA\ images, with somewhat larger size. We detected for the first time spicular structures, rising up to 15" above the limb with a width down to the image resolution and brightness temperature of $\sim$ 1800 K above the local background. No trace of spicules, either in emission or absorption, was found on the disk. Our results highlight ALMA's potential for the study of the quiet chromosphere.

Solar flares are due to the catastrophic release of magnetic energy in the Sun's corona, resulting in plasma heating, mass motions, particle acceleration, and radiation emitted from radio to $\gamma$-ray wavelengths. They are associated with global coronal eruptions of plasma into the interplanetary medium---coronal mass ejections---that can result in a variety of "space weather" phenomena. Flares release energy over a vast range of energies, from $\sim\!10^{23}$ ergs (nanoflares) to more than $10^{32}$ ergs. Solar flares are a phenomenon of general astrophysical interest, allowing detailed study of magnetic energy release, eruptive processes, shock formation and propagation, particle acceleration and transport, and radiative processes. Observations at radio wavelengths offer unique diagnostics of the physics of flares. To fully exploit these diagnostics requires the means of performing time-resolved imaging spectropolarimetry. Recent observations with the Jansky Very Large Array (JVLA) and the Expanded Owens Valley Solar Array (EOVSA), supported by extensive development in forward modeling, have demonstrated the power of the approach. The ngVLA has the potential to bring our understanding of flare processes to a new level through its combination of high spatial resolution, broad frequency range, and imaging dynamic range---especially when used in concert with multi-wavelength observations and data at hard X-ray energies.

The ngVLA has the potential to play a significant role in characterizing properties of the outer corona and the heating and acceleration of the solar wind into the inner heliosphere. In particular, using distant background sources to transilluminate the foreground corona and solar wind, a variety of radio propagation phenomena can be used to map plasma properties as a function of solar elongation and position angle throughout the solar cycle. These include angular broadening, interplanetary scintillations, and differential Faraday rotation, which can be used to map the solar wind velocity, determine properties of solar wind turbulence, and constrain the solar wind magnetic field. These observations will provide a global characterization of the solar wind that will be highly complementary to in situ observations made by various spacecraft. In addition, such observations can be used to probe disturbances in the solar wind -- coronal mass ejections, for example -- that may impact the near-Earth environment.

Energy stored in the magnetic field in the solar atmosphere above active regions is a key driver of all solar activity (e.g., solar flares and coronal mass ejections), some of which can affect life on Earth. Radio observations provide a unique diagnostic of the coronal magnetic fields that make them a critical tool for the study of these phenomena, using the technique of broadband radio imaging spectropolarimetry. Observations with the ngVLA will provide unique observations of coronal magnetic fields and their evolution, key inputs and constraints for MHD numerical models of the solar atmosphere and eruptive processes, and a key link between lower layers of the solar atmosphere and the heliosphere. In doing so they will also provide practical "research to operations" guidance for space weather forecasting.

Magnetic reconnection, the central engine that powers explosive phenomena throughout the Universe, is also perceived as one of the principal mechanisms for accelerating particles to high energies. Although various signatures of magnetic reconnection have been frequently reported, observational evidence that links particle acceleration directly to the reconnection site has been rare, especially for space plasma environments currently inaccessible to $\textit{in situ}$ measurements. Here we utilize broadband radio dynamic imaging spectroscopy available from the Karl G. Jansky Very Large Array to observe decimetric type III radio bursts in a solar jet with high angular ($\sim$20$''$), spectral ($\sim$1 %), and temporal resolution (50 milliseconds). These observations allow us to derive detailed trajectories of semi-relativistic (tens of keV) electron beams in the low solar corona with unprecedentedly high angular precision ($<0''.65$). We found that each group of electron beams, which corresponds to a cluster of type III bursts with 1-2-second duration, diverges from an extremely compact region ($\sim$600 km$^2$) in the low solar corona. The beam-diverging sites are located behind the erupting jet spire and above the closed arcades, coinciding with the presumed location of magnetic reconnection in the jet eruption picture supported by extreme ultraviolet/X-ray data and magnetic modeling. We interpret each beam-diverging site as a reconnection null point where multitudes of magnetic flux tubes join and reconnect. Our data suggest that the null points likely consist of a high level of density inhomogeneities possibly down to 10-km scales. These results, at least in the present case, strongly favor a reconnection-driven electron acceleration scenario.

As part of a wider search for radio emission from nearby systems known or suspected to contain extrasolar planets $\epsilon$ Eridani was observed by the Jansky Very Large Array (VLA) in the 2-4 GHz and 4-8 GHz frequency bands. In addition, as part of a separate survey of thermal emission from solar-like stars, $\epsilon$ Eri was observed in the 8-12 GHz and the 12-18 GHz bands of the VLA. Quasi-steady continuum radio emission from $\epsilon$ Eri was detected in the three high-frequency bands at levels ranging from 67-83 $\mu$Jy. No significant variability is seen in the quasi-steady emission. The emission in the 2-4 GHz emission, however, is shown to be the result of a circularly polarized (up to 50\%) radio pulse or flare of a few minutes duration that occurred at the beginning of the observation. We consider the astrometric position of the radio source in each frequency band relative to the expected position of the K2V star and the purported planet. The quasi-steady radio emission at frequencies $\ge \!8$ GHz is consistent with a stellar origin. The quality of the 4-8 GHz astrometry provides no meaningful constraint on the origin of the emission. The location of the 2-4 GHz radio pulse is $>2.5\sigma$ from the star yet, based on the ephemeris of Benedict et al. (2006), it is not consistent with the expected location of the planet either. If the radio pulse has a planetary origin, then either the planetary ephemeris is incorrect or the emission originates from another planet.

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.

We measured the center-to-limb variation of the brightness temperature, $T_b$, from ALMA full-disk images at two frequencies and inverted the solution of the transfer equation to obtain the electron temperature, $T_e$ as a function of optical depth, $\tau$. The ALMA images are very similar to AIA images at 1600\AA. The brightness temperature at the center of the disk is 6180 and 7250 K at 239 and 100 GHz respectively, with dispersions of 100 and 170 K. Plage regions stand out clearly in the 239/100 GHz intensity ratio, while faculae and filament lanes do not. The solar disk radius, reduced to 1 AU, is $961.1\pm2.5$ arcsec and $964.1\pm4.5$ arcsec at 239 and 100 GHz respectively. A slight but statistically significant limb brightening is observed at both frequencies. The inversion of the center-to-limb curves shows that $T_e$ varies linearly with the logarithm of optical depth for $0.34<\tau_{100\,GHz}<12$, with a slope $d\ln T_e/d\tau=-608$ K. Our results are 5% lower than predicted by the average quiet sun model C of Fontenla et al. (1993), but do not confirm previous reports that the mm-$\lambda$ solar spectrum is better fitted with models of the cell interior.

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.

Eruptive phenomena such as plasmoid ejections or jets are an important feature of solar activity with the potential for improving our understanding of the dynamics of the solar atmosphere. Such ejections are often thought to be signatures of the outflows expected in regions of fast magnetic reconnection. The 304 A EUV line of Helium, formed at around 10^5 K, is found to be a reliable tracer of such phenomena, but the determination of physical parameters from such observations is not straightforward. We have observed a plasmoid ejection from an X-ray bright point simultaneously at millimeter wavelengths with ALMA, at EUV wavelengths with AIA, in soft X-rays with Hinode/XRT. This paper reports the physical parameters of the plasmoid obtained by combining the radio, EUV and X-ray data. As a result, we conclude that the plasmoid can consist either of (approximately) isothermal 10^5 K plasma that is optically thin at 100 GHz, or else a 10^4 K core with a hot envelope. The analysis demonstrates the value of the additional temperature and density constraints that ALMA provides, and future science observations with ALMA will be able to match the spatial resolution of space-borne and other high-resolution telescopes.

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.

The absolute brightness temperature of the Sun at millimeter wavelengths is an important diagnostic of the solar chromosphere. Because the Sun is so bright, measurement of this property usually involves the operation of telescopes under extreme conditions and requires a rigorous performance assessment of the telescope. In this study, we establish solar observation and calibration techniques at 2.6-mm wavelength for the Nobeyama 45-m telescope and derive the absolute solar brightness temperature accurately. We tune the superconductor-insulator-superconductor (SIS) receiver by inducing different bias voltages onto the SIS mixer to prevent saturation. Then, we examine the linearity of the receiver system by comparing outputs derived from different tuning conditions. Further, we measure the lunar filled beam efficiency of the telescope using the New Moon, and then derive the absolute brightness temperature of the Sun. The derived solar brightness temperature is 7700+-310 K at 115 GHz. The telescope beam pattern is modeled as a summation of three Gaussian functions and derived using the solar limb. The real shape of the Sun is determined via deconvolution of the beam pattern from the observed map. Such well-calibrated single-dish observations are important for high-resolution chromospheric studies because they provide the absolute temperature scale missing from interferometer observations.

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.

Solar flares - the most powerful explosions in the solar system - are also efficient particle accelerators, capable of energizing a large number of charged particles to relativistic speeds. A termination shock is often invoked in the standard model of solar flares as a possible driver for particle acceleration, yet its existence and role have remained controversial. We present observations of a solar flare termination shock and trace its morphology and dynamics using high-cadence radio imaging spectroscopy. We show that a disruption of the shock coincides with an abrupt reduction of the energetic electron population. The observed properties of the shock are well-reproduced by simulations. These results strongly suggest that a termination shock is responsible, at least in part, for accelerating energetic electrons in solar flares.

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.

The Atacama Large Millimeter/submillimeter Array (ALMA) will be a valuable tool for observing the chromosphere of our Sun at (sub-)millimeter wavelengths at high spatial, temporal and spectral resolution and as such has great potential to address long-standing scientific questions in solar physics. In order to make the best use of this scientific opportunity, the Solar Simulations for the Atacama Large Millimeter Observatory Network has been initiated. A key goal of this international collaboration is to support the preparation and interpretation of future observations of the Sun with ALMA.

ALMA provides the necessary spatial, temporal and spectral resolution to explore central questions in contemporary solar physics with potentially far-reaching implications for stellar atmospheres and plasma physics. It can uniquely constraint the thermal and magnetic field structure in the solar chromosphere with measurements that are highly complementary to simultaneous observations with other ground-based and space-borne instruments. Here, we highlight selected science cases.

We have implemented control and data-taking software that makes it possible to scan the beams of individual ALMA antennas to perform quite complex patterns while recording the signals at high rates. We conducted test observations of the Sun in September and December, 2014. The data returned have excellent quality; in particular they allow us to characterize the noise and signal fluctuations present in this kind of observation. The fast-scan experiments included both Lissajous patterns covering rectangular areas, and double-circle patterns of the whole disk of the Sun and smaller repeated maps of specific disk-shaped targets. With the latter we find that we can achieve roughly Nyquist sampling of the Band~6 (230~GHz) beam in 60~s over a region 300$"$ in diameter. These maps show a peak-to-peak brightness-temperature range of up to 1000~K, while the time-series variability at any given point appears to be of order 0.5 percent RMS over times of a few minutes. We thus expect to be able to separate the noise contributions due to transparency fluctuations from variations in the Sun itself. Such timeseries have many advantages, in spite of the non-interferometric observations. In particular such data should make it possible to observe microflares in active regions and nanoflares in any part of the solar disk and low corona.

Magnetic flux ropes (MFRs) are believed to be at the heart of solar coronal mass ejections (CMEs). A well-known example is the prominence cavity in the low corona that sometimes makes up a three-part white-light (WL) CME upon its eruption. Such a system, which is usually observed in quiet-Sun regions, has long been suggested to be the manifestation of an MFR with relatively cool filament material collecting near its bottom. However, observational evidence of eruptive, filament-hosting MFR systems has been elusive for those originating in active regions. By utilizing multi-passband extreme-ultra-violet (EUV) observations from SDO/AIA, we present direct evidence of an eruptive MFR in the low corona that exhibits a hot envelope and a cooler core; the latter is likely the upper part of a filament that undergoes a partial eruption, which is later observed in the upper corona as the coiled kernel of a fast, WL CME. This MFR-like structure exists more than one hour prior to its eruption, and displays successive stages of dynamical evolution, in which both ideal and non-ideal physical processes may be involved. The timing of the MFR kinematics is found to be well correlated with the energy release of the associated long-duration C1.9 flare. We suggest that the long-duration flare is the result of prolonged energy release associated with the vertical current sheet induced by the erupting MFR.

We report observations of type III radio bursts at decimeter wavelengths (type IIIdm bursts) -- signatures of suprathermal electron beams propagating in the low corona -- using the new technique of radio dynamic imaging spectroscopy provided by the recently upgraded Karl G. Jansky Very Large Array (VLA). For the first time, type IIIdm bursts were imaged with high time and frequency resolution over a broad frequency band, allowing electron beam trajectories in the corona to be deduced. Together with simultaneous hard X-ray (HXR) and extreme ultraviolet (EUV) observations, we show these beams emanate from an energy release site located in the low corona at a height below ~15 Mm, and propagate along a bundle of discrete magnetic loops upward into the corona. Our observations enable direct measurements of the plasma density along the magnetic loops, and allow us to constrain the diameter of these loops to be less than 100 km. These over-dense and ultra-thin loops reveal the fundamentally fibrous structure of the Sun's corona. The impulsive nature of the electron beams, their accessibility to different magnetic field lines, and the detailed structure of the magnetic release site revealed by the radio observations indicate that the localized energy release is highly fragmentary in time and space, supporting a bursty reconnection model that involves secondary magnetic structures for magnetic energy release and particle acceleration.

Coronal hard X-ray (HXR) and continuum gamma-ray sources associated with the impulsive phase of solar flares have been the subject of renewed interest in recent years. They have been interpreted in terms of thin-target, nonthermal bremsstrahlung emission. This interpretation has led to rather extreme physical requirements in some cases. For example, in one case, essentially all of the electrons in the source must be accelerated to nonthermal energies to account for the coronal HXR source. In other cases, the extremely hard photon spectra of the coronal continuum gamma-ray emission suggest that the low energy cutoff of the electron energy distribution lies in the MeV energy range. Here we consider the role of inverse Compton scattering (ICS) as an alternate emission mechanism in both the ultra- and mildly relativistic regimes. It is known that relativistic electrons are produced during powerful flares; these are capable of up-scattering soft photospheric photons to HXR and gamma-ray energies. Previously overlooked is the fact that mildly relativistic electrons, generally produced in much greater numbers in flares of all sizes, can up-scatter EUV/SXR photons to HXR energies. We also explore ICS on anisotropic electron distributions and show that the resulting emission can be significantly enhanced over an isotropic electron distribution for favorable viewing geometries. We briefly review results from bremsstrahlung emission and reconsider circumstances under which nonthermal bremsstrahlung or ICS would be favored. Finally, we consider a selection of coronal HXR and gamma-ray events and find that in some cases the ICS is a viable alternative emission mechanism.

We present the first interferometric observation of a zebra-pattern radio burst with simultaneous high spectral (~ 1 MHz) and high time (20 ms) resolution. The Frequency-Agile Solar Radiotelescope (FASR) Subsystem Testbed (FST) and the Owens Valley Solar Array (OVSA) were used in parallel to observe the X1.5 flare on 14 December 2006. By using OVSA to calibrate the FST the source position of the zebra pattern can be located on the solar disk. With the help of multi-wavelength observations and a nonlinear force-free field (NLFFF) extrapolation, the zebra source is explored in relation to the magnetic field configuration. New constraints are placed on the source size and position as a function of frequency and time. We conclude that the zebra burst is consistent with a double-plasma resonance (DPR) model in which the radio emission occurs in resonance layers where the upper hybrid frequency is harmonically related to the electron cyclotron frequency in a coronal magnetic loop.

The coming decade will see the routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness. To capitalize on the new (or soon to be available) facilities such as SDO, ATST and FASR, and the challenges they present in the visualization and synthesis of multi-wavelength datasets, we propose that realistic, sophisticated, 3D active region and flare modeling is timely and critical, and will be a forefront of coronal studies over the coming decade. To make such modeling a reality, a broad, concerted effort is needed to capture the wealth of information resulting from the data, develop a synergistic modeling effort, and generate the necessary visualization, interpretation and model-data comparison tools to accurately extract the key physics.

Time-to-science is an important figure of merit for digital instrumentation serving the astronomical community. A digital signal processing (DSP) community is forming that uses shared hardware development, signal processing libraries, and instrument architectures to reduce development time of digital instrumentation and to improve time-to-science for a wide variety of projects. We suggest prioritizing technological development supporting the needs of this nascent DSP community. After outlining several instrument classes that are relying on digital instrumentation development to achieve new science objectives, we identify key areas where technologies pertaining to interoperability and processing flexibility will reduce the time, risk, and cost of developing the digital instrumentation for radio astronomy. These areas represent focus points where support of general-purpose, open-source development for a DSP community should be prioritized in the next decade. Contributors to such technological development may be centers of support for this DSP community, science groups that contribute general-purpose DSP solutions as part of their own instrumentation needs, or engineering groups engaging in research that may be applied to next-generation DSP instrumentation.

The magnetospheric emissions from extrasolar planets represent a science frontier for the next decade. All of the solar system giant planets and the Earth produce radio emissions as a result of interactions between their magnetic fields and the solar wind. In the case of the Earth, its magnetic field may contribute to its habitability by protecting its atmosphere from solar wind erosion and by preventing energetic particles from reaching its surface. Indirect evidence for at least some extrasolar giant planets also having magnetic fields includes the modulation of emission lines of their host stars phased with the planetary orbits, likely due to interactions between the stellar and planetary magnetic fields. If magnetic fields are a generic property of giant planets, then extrasolar giant planets should emit at radio wavelengths allowing for their direct detection. Existing observations place limits comparable to the flux densities expected from the strongest emissions. Additional sensitivity at low radio frequencies coupled with algorithmic improvements likely will enable a new means of detection and characterization of extrasolar planets within the next decade.

We describe microwave and hard X-ray observations of strong quasiperiodic pulsations from the GOES X1.3 solar flare on 15 June 2003. The radio observations were made jointly by the Owens Valley Solar Array (OVSA), the Nobeyama Polarimeter (NoRP), and the Nobeyama Radioheliograph (NoRH). Hard X-ray observations were made by the Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Using Fourier analysis, we study the frequency- and energy-dependent oscillation periods, differential phase, and modulation amplitudes of the radio and X-ray pulsations. Focusing on the more complete radio observations, we also examine the modulation of the degree of circular polarization and of the radio spectral index. The observed properties of the oscillations are compared with those derived from two simple models for the radio emission. In particular, we explicitly fit the observed modulation amplitude data to the two competing models. The first model considers the effects of MHD oscillations on the radio emission. The second model considers the quasi-periodic injection of fast electrons. We demonstrate that quasiperiodic acceleration and injection of fast electrons is the more likely cause of the quasiperiodic oscillations observed in the radio and hard X-ray emission, which has important implications for particle acceleration and transport in the flaring sources.

We report observations of a radio burst that occurred on the flare star AD Leonis over a frequency range of 1120-1620 MHz (lambda ~18-27 cm). These observations, made by the 305m telescope of the Arecibo Observatory, are unique in providing the highest time resolution (1 ms) and broadest spectral coverage (Delta nu/nu =0.36) of a stellar radio burst yet obtained. The burst was observed on 2005 April 9. It produced a peak flux density of ~500 mJy and it was essentially 100% right circularly polarized. The dynamic spectrum shows a rich variety of structure: patchy emission, diffuse bands, and narrowband, fast-drift striae. Focusing our attention on the fast-drift striae, we consider the possible role of dispersion and find that it requires rather special conditions in the source to be a significant factor. We suggest that the emission may be due to the cyclotron maser instability, a mechanism known to occur in planetary magnetospheres. We briefly explore implications of this possibility.

An interplanetary (IP) type-II-like radio burst is analyzed. It occurred on 2003 June 17-18 in association with a fast halo coronal mass ejection (CME), an M6.8 soft-X-ray (SXR) flare, and produced a solar proton event. Unlike coronal type II bursts and the ma jority of IP type II radio emissions, the IP type-II-like event associated with the fast halo CME on June 17-18 varies smoothly in time and frequency and has a frequency bandwidth that is several times larger than is typical for coronal and IP type II emissions. Moreover, the frequency change with time is inconsistent with that expected from plasma radiation associated with a CME-driven shock. I suggest that this IP type-II-like event, referred to here as an IP type II-S event, is not due to plasma radiation but, rather, incoherent synchrotron radiation from near-relativistic electrons entrained in the CME magnetic field, or in the sheath region between the shock and the CME driver. This event may be an example of a new and distinct class of interplanetary radio phenomenon.

We present radio and X-ray observations of an impulsive solar flare that was moderately intense in microwaves, yet showed very meager EUV and X-ray emission. The flare occurred on 2001 Oct 24 and was well-observed at radio wavelengths by the Nobeyama Radioheliograph (NoRH), the Nobeyama Radio Polarimeters (NoRP), and by the Owens Valley Solar Array (OVSA). It was also observed in EUV and X-ray wavelength bands by the TRACE, GOES, and Yohkoh satellites. We find that the impulsive onset of the radio emission is progressively delayed with increasing frequency relative to the onset of hard X-ray emission. In contrast, the time of flux density maximum is progressively delayed with decreasing frequency. The decay phase is independent of radio frequency. The simple source morphology and the excellent spectral coverage at radio wavelengths allowed us to employ a nonlinear chi-squared minimization scheme to fit the time series of radio spectra to a source model that accounts for the observed radio emission in terms of gyrosynchrotron radiation from MeV-energy electrons in a relatively dense thermal plasma. We discuss plasma heating and electron acceleration in view of the parametric trends implied by the model fitting. We suggest that stochastic acceleration likely plays a role in accelerating the radio-emitting electrons.

We report high-time-resolution, broadband spectroscopic observations of two radio bursts on the classical flare star AD Leonis. The observations were acquired by the 305 m telescope at Arecibo Observatory on 2003 June 13-14. Using the Wideband Arecibo Pulsar Processor, these observations sampled a total bandwidth of 400 MHz, distributed over a 500 MHz frequency range, 1120--1620 MHz, with a frequency resolution of 0.78 MHz and a time resolution of 10 ms. A radio burst observed on June 13 is characterized by the presence of multitudes of short duration ($\Delta t \sim$30 ms), high brightness temperature ($T_{b}>10^{14}$K), highly circularly polarized, fast-drift radio sub-bursts, with median bandwidths $\Delta \nu/\nu \sim$5%. The inverse drift rates are small, and have a symmetric distribution (both positive and negative frequency drifts) with a Gaussian FWHM inverse drift rate of 4.5$\times10^{-4}$ s/MHz. The fast-drift sub-bursts occur at a mean rate of 13 s$^{-1}$ and show no evidence for periodic recurrence. The fast-drift radio events on AD Leo are highly reminiscent of solar decimetric spike bursts. We suggest the emission is due to fundamental plasma radiation. A second highly circularly polarized radio burst, recorded June 14, has markedly different properties: a smoothly varying intensity profile characterized by a slow drift in frequency with time (-52 MHz s$^{-1}$). Under the assumption that the source is due to a disturbance propagating through the low corona, a source size of 0.1--1 R$_{\star}$ is inferred, implying a brightness temperature range 6$\times10^{11}$--6$\times10^{13}$K: another example of a coherent radio burst.

We explore the radio emission from the M9 dwarf, TVLM513-46546, at multiple radio frequencies, determining the flux spectrum of persistent radio emission, as well as constraining the levels of circular polarization. Detections at both 3.6 and 6 cm provide spectral index measurement $\alpha$ (where S$_{\nu} \propto \nu^{\alpha}$) of $-0.4\pm0.1$. A detection at 20 cm suggests that the spectral peak is between 1.4 and 5 GHz. The most stringent upper limits on circular polarization are at 3.6 and 6 cm, with $V/I <$15%. These characteristics agree well with those of typical parameters for early to mid M dwarfs, confirming that magnetic activity is present at levels comparable with those extrapolated from earlier M dwarfs. We apply analytic models to investigate the coronal properties under simple assumptions of dipole magnetic field geometry and radially varying nonthermal electron density distributions. Requiring the spectrum to be optically thin at frequencies higher than 5 GHz and reproducing the observed 3.6 cm fluxes constrains the magnetic field at the base to be less than about 500 G. There is no statistically significant periodicity in the 3.6 cm light curve, but it is consistent with low-level variability.

At sufficiently low frequencies, no ground-based radio array will be able to produce high resolution images while looking through the ionosphere. A space-based array will be needed to explore the objects and processes which dominate the sky at the lowest radio frequencies. An imaging radio interferometer based on a large number of small, inexpensive satellites would be able to track solar radio bursts associated with coronal mass ejections out to the distance of Earth, determine the frequency and duration of early epochs of nonthermal activity in galaxies, and provide unique information about the interstellar medium. This would be a "space-space" VLBI mission, as only baselines between satellites would be used. Angular resolution would be limited only by interstellar and interplanetary scattering.