skip to main content


Title: Prospects and challenges of numerical modeling of the Sun at millimeter wavelengths
The Atacama Large Millimeter/submillimeter Array (ALMA) offers new diagnostic possibilities that complement other commonly used diagnostics for the study of the Sun. In particular, ALMA’s ability to serve as an essentially linear thermometer of the chromospheric gas at unprecedented spatial resolution at millimeter wavelengths and future polarization measurements has great diagnostic potential. Solar ALMA observations are therefore expected to contribute significantly to answering long-standing questions about the structure, dynamics, and energy balance of the outer layers of the solar atmosphere. In this regard, current and future ALMA data are also important for constraining and further developing numerical models of the solar atmosphere, which in turn are often vital for the interpretation of observations. The latter is particularly important given the Sun’s highly intermittent and dynamic nature that involves a plethora of processes occurring over extended ranges in spatial and temporal scales. Realistic forward modeling of the Sun therefore requires time-dependent three-dimensional radiation magnetohydrodynamics that account for non-equilibrium effects and, typically as a separate step, detailed radiative transfer calculations, resulting in synthetic observables that can be compared to observations. Such artificial observations sometimes also account for instrumental and seeing effects, which, in addition to aiding the interpretation of observations, provide instructive tools for designing and optimizing ALMA’s solar observing modes. In the other direction, ALMA data in combination with other simultaneous observations enable the reconstruction of the solar atmospheric structure via data inversion techniques. This article highlights central aspects of the impact of ALMA for numerical modeling of the Sun and their potential and challenges, together with selected examples.  more » « less
Award ID(s):
1820613 2121632 2206424
NSF-PAR ID:
10387561
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Frontiers in Astronomy and Space Sciences
Volume:
9
ISSN:
2296-987X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    By direct measurements of the gas temperature, the Atacama Large Millimeter/submillimeter Array (ALMA) has yielded a new diagnostic tool to study the solar chromosphere. Here, we present an overview of the brightness-temperature fluctuations from several high-quality and high-temporal-resolution (i.e. 1 and 2 s cadence) time series of images obtained during the first 2 years of solar observations with ALMA, in Band 3 and Band 6, centred at around 3 mm (100 GHz) and 1.25 mm (239 GHz), respectively. The various datasets represent solar regions with different levels of magnetic flux. We perform fast Fourier and Lomb–Scargle transforms to measure both the spatial structuring of dominant frequencies and the average global frequency distributions of the oscillations (i.e. averaged over the entire field of view). We find that the observed frequencies significantly vary from one dataset to another, which is discussed in terms of the solar regions captured by the observations (i.e. linked to their underlying magnetic topology). While the presence of enhanced power within the frequency range 3–5 mHz is found for the most magnetically quiescent datasets, lower frequencies dominate when there is significant influence from strong underlying magnetic field concentrations (present inside and/or in the immediate vicinity of the observed field of view). We discuss here a number of reasons which could possibly contribute to the power suppression at around 5.5 mHz in the ALMA observations. However, it remains unclear how other chromospheric diagnostics (with an exception of H α line-core intensity) are unaffected by similar effects, i.e. they show very pronounced 3-min oscillations dominating the dynamics of the chromosphere, whereas only a very small fraction of all the pixels in the 10 ALMA datasets analysed here show peak power near 5.5 mHz. This article is part of the Theo Murphy meeting issue ‘High-resolution wave dynamics in the lower solar atmosphere’. 
    more » « less
  2. Abstract

    Atmospheric gravity waves (AGWs) are low-frequency, buoyancy-driven waves that are generated by turbulent convection and propagate obliquely throughout the solar atmosphere. Their proposed energy contribution to the lower solar atmosphere and sensitivity to atmospheric parameters (e.g., magnetic fields and radiative damping) highlight their diagnostic potential. We investigate AGWs near a quiet-Sun disk center region using multiwavelength data from the Interferometric Bidimensional Spectrometer and the Solar Dynamics Observatory. These observations showcase the complex wave behavior present in the entire acoustic-gravity wave spectrum. Using Fourier spectral analysis and local helioseismology techniques on simultaneously observed line core Doppler velocity and intensity fluctuations, we study both the vertical and horizontal properties of AGWs. Propagating AGWs with perpendicular group and phase velocities are detected at the expected temporal and spatial scales throughout the lower solar atmosphere. We also find previously unobserved, varied phase difference distributions among our velocity and intensity diagnostic combinations. Time–distance analysis indicates that AGWs travel with an average group speed of 4.5 km s−1, which is only partially described by a simple simulation, suggesting that high-frequency AGWs dominate the signal. Analysis of the median magnetic field (4.2 G) suggests that propagating AGWs are not significantly affected by quiet-Sun photospheric magnetic fields. Our results illustrate the importance of multiheight observations and the necessity of future work to properly characterize this observed behavior.

     
    more » « less
  3. Abstract Waves and oscillations are important solar phenomena, not only because they can propagate and dissipate energy in the chromosphere, but also because they carry information about the structure of the atmosphere in which they propagate. The nature of the 3 minute oscillations observed in the umbral region of sunspots is considered to be an effect of propagation of magnetohydrodynamic waves upward from below the photosphere. We present a study of sunspot oscillations and wave propagation in NOAA Active Region 12470 using an approximately 1 hr long data set acquired on 2015 December 17 by the Atacama Large Millimeter/submillimeter Array (ALMA), the Goode Solar Telescope (GST) operating at the Big Bear Solar Observatory, the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, and the Interface Region Imaging Spectrograph. The ALMA data are unique in providing a time series of direct temperature measurements in the sunspot chromosphere. The 2 s cadence of ALMA images allows us to well resolve the 3 minute periods typical of sunspot oscillations in the chromosphere. Fourier analysis is applied to ALMA Band 3 (∼100 GHz, ∼3 mm) and GST H α data sets to obtain power spectra as well as oscillation phase information. We analyzed properties of the wave propagation by combining multiple wavelengths that probe physical parameters of solar atmosphere at different heights. We find that the ALMA temperature fluctuations are consistent with that expected for a propagating acoustic wave, with a slight asymmetry indicating nonlinear steepening. 
    more » « less
  4. We report the detection of three large millimeter flaring events from the nearby Sun-like,ϵ Eridani, found in archival Atacama Large Millimeter/submillimeter Array (ALMA) 12 m and Atacama Compact Array observations at 1.33 mm taken from 2015 January 17 to 18 and 2016 October 24 to November 23, respectively. This is the first time that flares have been detected from a Sun-like star at millimeter wavelengths. The largest flare among our data was detected in the ALMA observations on 2015 January 17 from 20:09:10.4–21:02:49.3 UT with a peak flux density of 28 ± 7 mJy and a duration of 9 s. The peak brightness of the largest flare is 3.4 ± 0.9 × 10^14 erg s^−1 Hz^−1, a factor of >50× times brighter than the star’s quiescent luminosity and >10× brighter than solar flares observed at comparable wavelengths. We find changes in the spectral index (F ν ∝ ν α ) at the flare peak, with α = 1.81 ± 1.94 and a lower limit on the fractional linear polarization ∣Q/I∣ = 0.08 ± 0.12. This positive spectral index is more similar to millimeter solar flares, differing from M-dwarf flares also detected at millimeter wavelengths that exhibit steeply negative spectral indices. 
    more » « less
  5. Introduction: With the capture of the first high- resolution, in-situ images of Near-Earth Objects (NEOs) a couple of decades ago [1–4], the ubiquity of regolith and the granular nature of small objects in the Solar System became apparent. Benefiting from an increased access to high computing power, new numerical studies emerged, modeling granular structures forming and evolving as small bodies in the Solar System [5–7]. Now adding laboratory studies on granular material strength for asteroid and other small body applications [8,9], we are steadily progressing in our understanding of how regolith is shaping the interiors and surfaces of these worlds. In addition, our ever-more powerful observation capabilities are uncovering interesting dust-related phenomena in the outer skirts of our Solar System, in the form of activity at large heliocentric distances and rings [10–12]. We find that our recent progress in understanding the behavior of granular material in small body environments also has applications to the more distant worlds of Centaurs and Trans-Neptunian Objects (TNOs). Internal Strength: We currently deduce internal friction of rubble piles from the observation of large numbers of small asteroids and their rotation rates, combined with the associated numerical simulations [13,14]. In the laboratory, we study internal friction of simulant materials using shear strength measurements [8]. Combining observations, modeling, and laboratory work, the picture emerges of rubble pile interiors being composed of coarse grains in the mm to cm range. The irregular shapes of the grains lead to mechanical interlocking, thus generating the internal friction required to match observations of the asteroid population [8,9]. We find that the presence of a fine fraction in the confined interior of a rubble pile actually leads weaker internal strength [9]. Surface Strength: Deducing surface regolith strength for NEOs is usually performed via average slope measurements [15–17] or, most notably, observing the outcome of an impact of known energy [18]. In the laboratory, we measure the angle of repose of simulant material via pouring tests, as well as its bulk cohesion using shear strength measurements [8]. In some cases, this allows us to infer grain size ranges for various regions of the surface and subsurface of pictured NEOs, beyond the resolution of their in-situ images. Surface Activity: The Rosetta mission revealed that a number of activity events on comet 67P/Churyumov–Gerasimenko were linked to active surface geology, most notably avalanches and cliff collapses [19]. In addition, the role of regolith strength in asteroid disruption patterns has been inferred from numerical simulations of rotating rubble piles [20]. By studying strength differences in simulant samples, it becomes apparent that a difference in cohesion between a surface and its subsurface layer can lead to activity events with surface mass shedding, without the presence of volatiles sublimating as a driver [8]. We show that such differences in surface strength can be brought upon by a depletion in fine grains or a change in composition (e.g. depletion in water ice) and could account for regular activity patterns on small bodies, independently of their distance to the Sun. This is of particular interest to the study of Centaur activity and a potential mechanism for feeding ring systems. 
    more » « less