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.

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.

Numerical simulations of the solar chromosphere predict a diverse thermal structure with both hot and cool regions. Observations of plage regions in particular feature broader and brighter chromospheric lines, which suggest that they are formed in hotter and denser conditions than in the quiet-Sun, but also implies a non-thermal component whose source is unclear. We revisit the problem of the stratification of temperature and microturbulence in plage now adding millimeter continuum observations provided by ALMA to inversions of near-ultraviolet IRIS spectra as a powerful new diagnostic to disentangle the two parameters. We fit cool chromospheric holes and track the fast evolution of compact mm brightenings in the plage region. We use the STiC non-LTE inversion code to simultaneously fit real ultraviolet and millimeter spectra in order to infer the thermodynamic parameters of the plasma. We confirm the anticipated constraining potential of ALMA in non-LTE inversions of the solar chromosphere. We find significant differences between the inversion results of IRIS data alone compared to the results of a combination with the mm data: the IRIS+ALMA inversions have increased contrast and temperature range, and tend to favor lower values of microturbulence in the chromosphere of plage. The average brightness temperature of the plage region at 1.25 mm is 8500 K, but the ALMA maps also show much cooler ($\sim3000$ K) and hotter ($\sim11\,000$ K) evolving features partially seen in other diagnostics. To explain the former, the inversions require the existence of localized, low temperature regions in the chromosphere where molecules such as CO could form. The hot features could sustain such high temperatures due to non-equilibrium hydrogen ionization effects in a shocked chromosphere - a scenario that is supported by low-frequency shock wave patterns found in the MgII lines probed by IRIS.

Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest (KOIs) with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013. Aims. We ensure a comprehensive characterisation of the investigated systems by analysing Kepler data combined with new ground-based transit data using a photodynamical model. This method is applied to the Kepler-82 system leading to its first dynamic analysis. Methods. In order to provide a coherent description of all observations simultaneously, we combine the numerical integration of the gravitational dynamics of a system over the time span of observations with a transit light curve model. To explore the model parameter space, this photodynamical model is coupled with a Markov chain Monte Carlo algorithm. Results. The Kepler-82b/c system shows sinusoidal TTVs due to their near 2:1 resonance dynamical interaction. An additional chopping effect in the TTVs of Kepler-82c hints to a further planet near the 3:2 or 3:1 resonance. We photodynamically analysed Kepler long- and short-cadence data and three new transit observations obtained by KOINet between 2014 and 2018. Our result reveals a non-transiting outer planet with a mass of $m_f=20.9\pm1.0\;M_\bigoplus$ near the 3:2 resonance to the outermost known planet, Kepler-82c. Furthermore, we determined the densities of planets b and c to the significantly more precise values $\rho_b=0.98_{-0.14}^{+0.10}\;\text{g cm}^{-3}$ and $\rho_c=0.494_{-0.077}^{+0.066}\;\text{g cm}^{-3}$.

From its discovery, the low density transiting Neptune HAT-P-26b showed a 2.1-sigma detection drift in its spectroscopic data, while photometric data showed a weak curvature in the timing residuals that required further follow-up observations to be confirmed. To investigate this suspected variability, we observed 11 primary transits of HAT-P-26b between March, 2015 and July, 2018. For this, we used the 2.15 meter Jorge Sahade Telescope placed in San Juan, Argentina, and the 1.2 meter STELLA and the 2.5 meter Nordic Optical Telescope, both located in the Canary Islands, Spain. To add upon valuable information on the transmission spectrum of HAT-P-26b, we focused our observations in the R-band only. To contrast the observed timing variability with possible stellar activity, we carried out a photometric follow-up of the host star along three years. We carried out a global fit to the data and determined the individual mid-transit times focusing specifically on the light curves that showed complete transit coverage. Using bibliographic data corresponding to both ground and space-based facilities, plus our new characterized mid-transit times derived from parts-per-thousand precise photometry, we observed indications of transit timing variations in the system, with an amplitude of ~4 minutes and a periodicity of ~270 epochs. The photometric and spectroscopic follow-up observations of this system will be continued in order to rule out any aliasing effects caused by poor sampling and the long-term periodicity.

The aim of this study is to demonstrate how the logarithmic millimeter continuum gradient observed using the Atacama Large Millimeter/submillimeter Array (ALMA) may be used to estimate optical thickness in the solar atmosphere. We discuss how using multi-wavelength millimeter measurements can refine plasma analysis through knowledge of the absorption mechanisms. Here we use sub-band observations from the publicly available science verification (SV) data, whilst our methodology will also be applicable to regular ALMA data. The spectral resolving capacity of ALMA SV data is tested using the enhancement coincident with an X-ray Bright Point (XBP) and from a plasmoid ejection event near active region NOAA12470 observed in Band 3 (84-116 GHz) on 17/12/2015. We compute the interferometric brightness temperature light-curve for both features at each of the four constituent sub-bands to find the logarithmic millimetre spectrum. We compared the observed logarithmic spectral gradient with the derived relationship with optical thickness for an isothermal plasma to estimate the structure's optical thicknesses. We conclude, within 90% confidence, that the stationary enhancement has an optical thickness between $0.02 \leq \tau \leq 2.78$, and that the moving enhancement has $0.11 \leq \tau \leq 2.78$, thus both lie near to the transition between optically thin and thick plasma at 100 GHz. From these estimates, isothermal plasmas with typical Band 3 background brightness temperatures would be expected to have electron temperatures of $\sim 7370 - 15300$ K for the stationary enhancement and between $\sim 7440 - 9560$ K for the moving enhancement, thus demonstrating the benefit of sub-band ALMA spectral analysis.

Solar observations with the Atacama Large Millimeter/submillimeter Array (ALMA) provide us with direct measurements of the brightness temperature in the solar chromosphere. We study the temperature distributions obtained with ALMA Band 6 (in four sub-bands at 1.21, 1.22, 1.29, and 1.3 mm) for various areas at, and in the vicinity of, a sunspot, comprising quasi-quiet and active regions with different amounts of underlying magnetic fields. We compare these temperatures with those obtained at near- and far-ultraviolet (UV) wavelengths (and with the line-core intensities of the optically-thin far-UV spectra), co-observed with the Interface Region Imaging Spectrograph (IRIS) explorer. These include the emission peaks and cores of the Mg II k 279.6 nm and Mg II h 280.4 nm lines as well as the line cores of C II 133.4 nm, O I 135.6 nm, and Si IV 139.4 nm, sampling the mid-to-high chromosphere and the low transition region. Splitting the ALMA sub-bands resulted in an slight increase of spatial resolution in individual temperature maps, thus, resolving smaller-scale structures compared to those produced with the standard averaging routines. We find that the radiation temperatures have different, though somewhat overlapping, distributions in different wavelengths and in the various magnetic regions. Comparison of the ALMA temperatures with those of the UV diagnostics should, however, be interpreted with great caution, the former is formed under the local thermodynamic equilibrium (LTE) conditions, the latter under non-LTE. The mean radiation temperature of the ALMA Band 6 is similar to that extracted from the IRIS C II line in all areas with exception of the sunspot and pores where the C II poses higher radiation temperatures. In all magnetic regions, the Mg II lines associate with the lowest mean radiation temperatures in our sample. These will provide constraints for future numerical models.

In this work we use solar observations with the ALMA radio telescope at the wavelength of 1.21 mm. The aim of the analysis is to improve understanding of the solar chromosphere, a dynamic layer in the solar atmosphere between the photosphere and corona. The study has an observational and a modeling part. In the observational part full-disc solar images are analyzed. Based on a modified FAL atmospheric model, radiation models for various observed solar structures are developed. Finally, the observational and modeling results are compared and discussed.

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.

The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised to follow up transiting planet candidate KOIs with large transit timing variations (TTVs). Its main goal is to complete their TTV curves, as the Kepler telescope no longer observes the original Kepler field. Combining Kepler and new ground-based transit data we improve the modelling of these systems. To this end, we have developed a photodynamical model, and we demonstrate its performance using the Kepler-9 system as an example. Our comprehensive analysis combines the numerical integration of the system's dynamics over the time span of the observations along with the transit light curve model. This model is coupled with a Markov chain Monte Carlo algorithm, allowing the exploration of the model parameter space. Applied to the Kepler-9 long cadence data, short cadence data and 13 new transit observations collected by KOINet between the years 2014 to 2017, our modelling provides well constrained predictions for the next transits and the system's parameters. We have determined the densities of the planets Kepler-9b and 9c to the very precise values of rho_b = 0.439 +/-0.023 g/cm3 and rho_c = 0.322 +/- 0.017 g/cm3. Our analysis reveals that Kepler-9c will stop transiting in about 30 years. This results from strong dynamical interactions between Kepler-9b and 9c, near 2:1 resonance, that leads to a periodic change in inclination. Over the next 30 years the inclination of Kepler-9c (-9b) will decrease (increase) slowly. This should be measurable by a substantial decrease (increase) in the transit duration, in as soon as a few years' time. Observations that contradict this prediction might indicate the presence of additional objects. If this prediction proves true, this behaviour opens up a unique chance to scan the different latitudes of a star.

During its four years of photometric observations, the Kepler space telescope detected thousands of exoplanets and exoplanet candidates. One of Kepler's greatest heritages has been the confirmation and characterization of hundreds of multi-planet systems via Transit Timing Variations (TTVs). However, there are many interesting candidate systems displaying TTVs on such long time scales that the existing Kepler observations are of insufficient length to confirm and characterize them by means of this technique. To continue with Kepler's unique work we have organized the "Kepler Object of Interest Network" (KOINet). The goals of KOINet are, among others, to complete the TTV curves of systems where Kepler did not cover the interaction timescales well. KOINet has been operational since March, 2014. Here we show some promising first results obtained from analyzing seven primary transits of KOI-0410.01, KOI-0525.01, KOI-0760.01, and KOI-0902.01 in addition to Kepler data, acquired during the first and second observing seasons of KOINet. While carefully choosing the targets we set demanding constraints about timing precision (at least 1 minute) and photometric precision (as good as 1 part per thousand) that were achieved by means of our observing strategies and data analysis techniques. For KOI-0410.01, new transit data revealed a turn-over of its TTVs. We carried out an in-depth study of the system, that is identified in the NASA's Data Validation Report as false positive. Among others, we investigated a gravitationally-bound hierarchical triple star system, and a planet-star system. While the simultaneous transit fitting of ground and space-based data allowed for a planet solution, we could not fully reject the three-star scenario. New data, already scheduled in the upcoming 2018 observing season, will set tighter constraints on the nature of the system.

Various solar features can be seen on maps of the Sun in the mm and sub-mm wavelength range. The recently installed Atacama Large Millimeter/submillimeter Array (ALMA) is capable of observing the Sun in that wavelength range with an unprecedented spatial, temporal and spectral resolution. To interpret solar observations with ALMA the first important step is to compare ALMA maps with simultaneous images of the Sun recorded in other spectral ranges. First we identify different structures in the solar atmosphere seen in the optical, IR and EUV parts of the spectrum (quiet Sun (QS), active regions (AR), prominences on the disc, magnetic inversion lines (IL), coronal holes (CH) and coronal bright points (CBPs)) in a full disc solar ALMA image. The second aim is to measure the intensities (brightness temperatures) of those structures and compare them with the corresponding QS level. A full disc solar image at 1.21 mm obtained on December 18, 2015 during a CSV-EOC campaign with ALMA is calibrated and compared with full disc solar images from the same day in H\alpha, in He I 1083 nm core, and with SDO images (AIA at 170 nm, 30.4 nm, 21.1 nm, 19.3 nm, and 17.1 nm and HMI magnetogram). The brightness temperatures of various structures are determined by averaging over corresponding regions of interest in the ALMA image. Positions of the QS, ARs, prominences on the disc, ILs, CHs and CBPs are identified in the ALMA image. At 1.21 mm ARs appear as bright areas (but sunspots are dark), while prominences on the disc and CHs are not discernible from the QS background, although having slightly less intensity than surrounding QS regions. ILs appear as large, elongated dark structures and CBPs correspond to ALMA bright points. These results are in general agreement with sparse earlier measurements at similar wavelengths. The identification of CBPs represents the most important new result.

Although observational data unequivocally point out to the presence of chromospheres in red giant stars, no attempts have been made so far to model them using 3D hydrodynamical model atmospheres. We therefore compute an exploratory 3D hydrodynamical model atmosphere for a cool red giant in order to study the dynamical and thermodynamic properties of its chromosphere, as well as the influence of the chromosphere on its observable properties. 3D radiation hydrodynamics simulations are carried out with the CO5BOLD model atmosphere code for a star with the atmospheric parameters (Teff=4010 K, log g=1.5, [M/H]=0.0), which are similar to those of the K-type giant star Aldebaran (alpha Tau). ... we compute the emergent continuum intensity maps at different wavelengths, spectral line profiles of Ca II K, the Ca II infrared triplet line at 854.2nm, and H alpha, as well as the spectral energy distribution (SED) of the emergent radiative flux. The initial model quickly develops a dynamical chromosphere characterised by propagating and interacting shock waves. The peak temperatures in the chromospheric shock fronts reach values on the order of up to 5000 K although the shock fronts remain quite narrow. Like for the Sun, the gas temperature distribution in the upper layers is composed of a cool component due to adiabatic cooling in the expanding post-shock regions and a hot component due to shock waves. For this red giant model, the hot component is a rather flat high-temperature tail, which nevertheless affects the resulting average temperatures significantly. The simulations show that the atmospheres of red giant stars are dynamic and intermittent. Consequently, many observable properties cannot be reproduced with one-dimensional static models but demand for advanced 3D HD modelling. Furthermore, including a chromosphere in the models might produce significant contributions to the emergent UV flux.

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.

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.

Solar "magnetic tornadoes" are produced by rotating magnetic field structures that extend from the upper convection zone and the photosphere to the corona of the Sun. Recent studies show that such rotating features are an integral part of atmospheric dynamics and occur on a large range of spatial scales. A systematic statistical study of magnetic tornadoes is a necessary next step towards understanding their formation and their role for the mass and energy transport in the solar atmosphere. For this purpose, we have developed a new automatic detection method for chromospheric swirls, i.e. the observable signature of solar tornadoes or, more generally, chromospheric vortex flows and rotating motions. Unlike the previous studies that relied on visual inspections, our new method combines a line integral convolution (LIC) imaging technique and a scalar quantity which represents a vortex flow on a two-dimensional plane. We have tested two detection algorithms, based on the enhanced vorticity and vorticity strength quantities, by applying them to 3D numerical simulations of the solar atmosphere with CO5BOLD. We conclude that the vorticity strength method is superior compared to the enhanced vorticity method in all aspects. Applying the method to a numerical simulation of the solar atmosphere revealed very abundant small-scale, short-lived chromospheric vortex flows that had not been found by visual inspection before.

Flare-driven coronal rain can manifest from rapidly cooled plasma condensations near coronal loop-tops in thermally unstable post-flare arcades. We detect 5 phases that characterise the post-flare decay: heating, evaporation, conductive cooling dominance for ~120 s, radiative / enthalpy cooling dominance for ~4700 s and finally catastrophic cooling occurring within 35-124 s leading to rain strands with s periodicity of 55-70 s. We find an excellent agreement between the observations and model predictions of the dominant cooling timescales and the onset of catastrophic cooling. At the rain formation site we detect co-moving, multi-thermal rain clumps that undergo catastrophic cooling from ~1 MK to ~22000 K. During catastrophic cooling the plasma cools at a maximum rate of 22700 K s-1 in multiple loop-top sources. We calculated the density of the EUV plasma from the DEM of the multi-thermal source employing regularised inversion. Assuming a pressure balance, we estimate the density of the chromospheric component of rain to be 9.21x10^11 +-1.76x10^11 cm-3 which is comparable with quiescent coronal rain densities. With up to 8 parallel strands in the EUV loop cross section, we calculate the mass loss rate from the post-flare arcade to be as much as 1.98x10^12 +/-4.95x10^11 g s-1. Finally, we reveal a close proximity between the model predictions of 10^5.8 K and the observed properties between 10^5.9 K and 10^6.2 K, that defines the temperature onset of catastrophic cooling. The close correspondence between the observations and numerical models suggests that indeed acoustic waves (with a sound travel time of 68 s) could play an important role in redistributing energy and sustaining the enthalpy-based radiative cooling.

Using radiation magnetohydrodynamic simulations of the solar atmospheric layers from the upper convection zone to the lower corona, we investigate the self-consistent excitation of slow magneto-acoustic body waves (slow modes) in a magnetic flux concentration. We find that the convective downdrafts in the close surroundings of a two-dimensional flux slab "pump" the plasma inside it in the downward direction. This action produces a downflow inside the flux slab, which encompasses ever higher layers, causing an upwardly propagating rarefaction wave. The slow mode, excited by the adiabatic compression of the downflow near the optical surface, travels along the magnetic field in the upward direction at the tube speed. It develops into a shock wave at chromospheric heights, where it dissipates, lifts the transition region, and produces an offspring in the form of a compressive wave that propagates further into the corona. In the wake of downflows and propagating shock waves, the atmosphere inside the flux slab in the chromosphere and higher tends to oscillate with a period of $\nu\approx 4$~mHz. We conclude that this process of "magnetic pumping" is a most plausible mechanism for the direct generation of longitudinal chromospheric and coronal compressive waves within magnetic flux concentrations, and it may provide an important heat source in the chromosphere. It may also be responsible for certain types of dynamic fibrils.

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.

Here, we present a set of time-dependent 3D RMHD simulations of a M-dwarf star representative of AD Leo, which extend from the upper convection zone into the chromosphere. The 3D model atmospheres are characterized by a very dynamic and intermittent structure on small spatial and temporal scales and a wealth of physical processes, which by nature cannot be described by means of 1D static model atmospheres. Artificial observations of these models imply that a combination of complementary diagnostics such as Ca II lines and the continuum intensity from UV to millimeter wavelengths, probe various properties of the dynamics, thermal and magnetic structure of the photosphere and the chromosphere and thus provide measures of stellar activity, which can be compared to observations. The complicated magnetic field structure and its imprint in synthetic diagnostics may have important implications for the understanding and characterization of stellar activity and with it possibly for the evaluation of planetary habitability around active M-dwarf stars.

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.

We present the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarf stars. We demonstrate that convective energy transfer is seriously impeded by magnetic fields when the plasma-beta parameter, the thermal to magnetic pressure ratio, becomes smaller than unity. The critical field strength that inhibits convection in the photosphere of white dwarfs is in the range B = 1-50 kG, which is much smaller than the typical 1-1000 MG field strengths observed in magnetic white dwarfs, implying that these objects have radiative atmospheres. We have then employed evolutionary models to study the cooling process of high-field magnetic white dwarfs, where convection is entirely suppressed during the full evolution (B > 10 MG). We find that the inhibition of convection has no effect on cooling rates until the effective temperature (Teff) reaches a value of around 5500 K. In this regime, the standard convective sequences start to deviate from the ones without convection owing to the convective coupling between the outer layers and the degenerate reservoir of thermal energy. Since no magnetic white dwarfs are currently known at the low temperatures where this coupling significantly changes the evolution, effects of magnetism on cooling rates are not expected to be observed. This result contrasts with a recent suggestion that magnetic white dwarfs with Teff < 10,000 K cool significantly slower than non-magnetic degenerates.

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.

The Atacama Large Millimeter/submillimeter Array offers an unprecedented view of our Sun at sub-/millimeter wavelengths. The high spatial, temporal, and spectral resolution facilitates the measurement of gas temperatures and magnetic fields in the solar chromosphere with high precision. The anticipated results will revolutionize our understanding of the solar atmosphere and may in particular result in major steps towards solving the coronal heating problem. Based on state-of-the-art 3D radiation magnetohydrodynamic simulations, we calculate the emergent continuum intensity (and thus brightness temperature maps) in the wavelength range accessed by ALMA and simulate instrumental effects for different array configurations. First results show that the local gas temperature can be closely mapped with ALMA and that much of the complex small-scale chromospheric pattern can be resolved.

New and advanced space-based observing facilities continue to lower the resolution limit and detect solar coronal loops in greater detail. We continue to discover even finer sub-structures within coronal loop cross sections, in order to understand the nature of the solar corona. Here, we push this lower limit further to search for the finest coronal loop sub-structures, through taking advantage of the resolving power of the Swedish 1- m Solar Telescope (SST) / CRisp Imaging Spectro-Polarimeter (CRISP), together with co-observations from the Solar Dynamics Observatory (SDO) / Atmospheric Image Assembly (AIA). High resolution imaging of the chromospheric H-alpha 656.28 nm spectral line core and wings can, under certain circumstances, allow one to deduce the topology of the local magnetic environment of the solar atmosphere where its observed. Here, we study post-flare coronal loops, which become filled with evaporated chromosphere that rapidly condenses into chromospheric clumps of plasma (detectable in H-alpha) known as a coronal rain, to investigate their fine-scale structure. We identify, through analysis of three datasets, large-scale catastrophic cooling in coronal loop-tops and the existence of multi-thermal, multi-stranded sub-structures. Many cool strands even extend fully-intact from loop-top to foot-point. We discover that coronal loop fine-scale strands can appear bunched with as many as 8 parallel strands, within an AIA coronal loop cross-section. The strand number density vs cross-sectional width distribution, as detected by CRISP within AIA-defined coronal loops, most-likely peaks at well below 100 km and currently 69% of the sub-structure strands are statistically unresolved in AIA coronal loops.

There have been ubiquitous observations of wave-like motions in the solar atmosphere for decades. Recent improvements to space- and ground-based observatories have allowed the focus to shift to smaller magnetic structures on the solar surface. In this paper, high-resolution ground-based data taken using the Swedish 1 m Solar Telescope is combined with co-spatial and co-temporal data from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) satellite to analyze running penumbral waves (RPWs). RPWs have always been thought to be radial wave propagation that occurs within sunspots. Recent research has suggested that they are in fact upwardly propagating field-aligned waves (UPWs). Here, RPWs within a solar pore are observed for the first time and are interpreted as UPWs due to the lack of a penumbra that is required to support RPWs. These UPWs are also observed co-spatially and co-temporally within several SDO/AIA elemental lines that sample the transition region and low corona. The observed UPWs are traveling at a horizontal velocity of around 17 +- 0.5 km s-1 and a minimum vertical velocity of 42 +- 21 km s-1. The estimated energy of the waves is around 150 W m-2, which is on the lower bound required to heat the quiet-Sun corona. This is a new, yet unconsidered source of wave energy within the solar chromosphere and low corona.

High-resolution observations with the Swedish 1-m Solar Telescope (SST) and the Solar Dynamics Observatory (SDO) reveal rotating magnetic field structures that extend from the solar surface into the chromosphere and the corona. These so-called magnetic tornadoes are primarily detected as rings or spirals of rotating plasma in the Ca II 854.2 nm line core (also known as chromospheric swirls). Detailed numerical simulations show that the observed chromospheric plasma motion is caused by the rotation of magnetic field structures, which again are driven by photospheric vortex flows at their footpoints. Under the right conditions, two vortex flow systems are stacked on top of each other. We refer to the lower vortex, which extends from the low photosphere into the convection zone, as intergranular vortex flow (IVF). Once a magnetic field structure is co-located with an IVF, the rotation is mediated into the upper atmospheric layers and an atmospheric vortex flow (AVF, or magnetic tornado) is generated. In contrast to the recent work by Shelyag et al., we demonstrate that particle trajectories in a simulated magnetic tornado indeed follow spirals and argue that the properties of the trajectories decisively depend on the location in the atmosphere and the strength of the magnetic field.

Observations in the 171 AA channel of the Atmospheric Imaging Assembly of the space-borne Solar Dynamics Observatory show tornadoes-like features in the atmosphere of the Sun. These giant tornadoes appear as dark, elongated and apparently rotating structures in front of a brighter background. This phenomenon is thought to be produced by rotating magnetic field structures that extend throughout the atmosphere. We characterize giant tornadoes through a statistical analysis of properties like spatial distribution, lifetimes, and sizes. A total number of 201 giant tornadoes are detected in a period of 25 days, suggesting that on average about 30 events are present across the whole Sun at a time close to solar maximum. Most tornadoes appear in groups and seem to form the legs of prominences, thus serving as plasma sources/sinks. Additional Halpha observations with the Swedish 1-m Solar Telescope imply that giant tornadoes rotate as a structure although clearly exhibiting a thread-like structure. We observe tornado groups that grow prior to the eruption of the connected prominence. The rotation of the tornadoes may progressively twist the magnetic structure of the prominence until it becomes unstable and erupts. Finally, we investigate the potential relation of giant tornadoes to other phenomena, which may also be produced by rotating magnetic field structures. A comparison to cyclones, magnetic tornadoes and spicules implies that such events are more abundant and short-lived the smaller they are. This comparison might help to construct a power law for the effective atmospheric heating contribution as function of spatial scale.

Chromospheric swirls are the observational signatures of rotating magnetic field structures in the solar atmosphere, also known as magnetic tornadoes. Swirls appear as dark rotating features in the core of the spectral line of singly ionized calcium at a wavelength of 854.2 nm. This signature can be very subtle and difficult to detect given the dynamic changes in the solar chromosphere. Important steps towards a systematic and objective detection method are the compilation and characterization of a statistically significant sample of observed and simulated chromospheric swirls. Here, we provide a more exact definition of the chromospheric swirl phenomenon and also present a first morphological classification of swirls with three types: (I) Ring, (II) Split, (III) Spiral. We also discuss the nature of the magnetic field structures connected to tornadoes and the influence of limited spatial resolution on the appearance of their photospheric footpoints.

We present first results from three-dimensional radiation magnetohydrodynamic simulations of M-type dwarf stars with CO5BOLD. The local models include the top of the convection zone, the photosphere, and the chromosphere. The results are illustrated for models with an effective temperature of 3240 K and a gravitational acceleration of log g = 4.5, which represent analogues of AD Leo. The models have different initial magnetic field strengths and field topologies. This first generation of models demonstrates that the atmospheres of M-dwarfs are highly dynamic and intermittent. Magnetic fields and propagating shock waves produce a complicated fine-structure, which is clearly visible in synthetic intensity maps in the core of the Ca II K spectral line and also at millimeter wavelengths. The dynamic small-scale pattern cannot be described by means of one-dimensional models, which has important implications for the construction of semi-empirical model atmospheres and thus for the interpretation of observations in general. Detailed three-dimensional numerical simulations are valuable in this respect. Furthermore, such models facilitate the analysis of small-scale processes, which cannot be observed on stars but nevertheless might be essential for understanding M-dwarf atmospheres and their activity. An example are so-called "magnetic tornadoes", which have recently been found on the Sun and are presented here in M-dwarf models for the first time.

The interactive IDL-based CO5BOLD Analysis Tool (CAT) was developed to facilitate an easy and quick analysis of numerical simulation data produced with the 2D/3D radiation magnetohydrodynamics code CO5BOLD. The basic mode of operation is the display and analysis of cross-sections through a model either as 2D slices or 1D graphs. A wide range of physical quantities can be selected. Further features include the export of models into VAPOR format or the output of images and animations. A short overview including scientific analysis examples is given.

Three-dimensional numerical simulations with CO5BOLD, a new radiation hydrodynamics code, result in a dynamic, thermally bifurcated model of the non-magnetic chromosphere of the quiet Sun. The 3-D model includes the middle and low chromosphere, the photosphere, and the top of the convection zone, where acoustic waves are excited by convective motions. While the waves propagate upwards, they steepen into shocks, dissipate, and deposit their mechanical energy as heat in the chromosphere. Our numerical simulations show for the first time a complex 3-D structure of the chromospheric layers, formed by the interaction of shock waves. Horizontal temperature cross-sections of the model chromosphere exhibit a network of hot filaments and enclosed cool regions. The horizontal pattern evolves on short time-scales of the order of typically 20 - 25 seconds, and has spatial scales comparable to those of the underlying granulation. The resulting thermal bifurcation, i.e., the co-existence of cold and hot regions, provides temperatures high enough to produce the observed chromospheric UV emission and -- at the same time -- temperatures cold enough to allow the formation of molecules (e.g., carbon monoxide). Our 3-D model corroborates the finding by Carlsson & Stein (1994) that the chromospheric temperature rise of semi-empirical models does not necessarily imply an increase in the average gas temperature but can be explained by the presence of substantial spatial and temporal temperature inhomogeneities.