Lingezicht Astrophysics Reports 2 – ADS – August 26, 2021 – newest Solar Activities and their influences in the Heliosphere and Planetary Atmospheres Workshop NIT Calicut, India, March 8-14 2021

Compendium solar spectrum formation Robert J. Rutten1,2,3

1 Lingezicht Astrophysics, Deil, The Netherlands 2 Institute of Theoretical Astrophysics, University in Oslo, Oslo, Norway 3 Rosseland Centre for Solar Physics, University in Oslo, Oslo, Norway

Myself when young did eagerly frequent doctor and saint, and heard great argument about it and about; but evermore came out by the same door as in I went. Omar Khayyam (1048–1131)

Abstract. The solar spectrum conveys most of our diagnostics to find out how our star works. They must be understood for utilization but solar spectrum formation is complex because the interaction of matter and radiation within the solar atmosphere suffers non-local control in space, wavelength, and time. These complexities are summarized and illustrated with classic literature. They combine in chromospheric spectrum formation.

Contents 51 coronal heating by acoustic or gravity waves 51 k − ω diagram 1 Preamble: on-line format and on-line material2 53 solar p-mode history 2 Introduction3 3 RT basics5 4 spectral bandwidth and vision Intensity ...... 5 6 equivalent temperatures Matter–radiation interactions ...... 5 7 intensity calibration Bound-bound, bound-free, free-free interactions . . . .6 54 Schuster-Schwarzschild approximation Pair combinations ...... 6 37 Milne-Eddington approximation Equilibria ...... 6 30 contribution, response, sensitivity Bound-bound redistribution ...... 7 16 extinction diagrams Bound-free redistribution ...... 7 19 deepest wavelengths Free-free redistribution ...... 8 15 white-light corona Thomson and Rayleigh scattering ...... 8 17 emission measure, DEM Extinction and emissivity ...... 8 10 dielectronic√ ionization and recombination Source function ...... 8 20 √ε law Line source function ...... 9 21 ε validity Transfer equation along the beam ...... 10 23 epsilon is small: ALI, acceleration, net radiative bracket Transfer equation in optical depth ...... 10 24 no Thomas photoelectric control 4 LTE: Holweger models 10 22 Menzel departure coefficients 5 NLS: Auer & Mihalas modeling 11 25 photon pumping, lasering, Wien transcription 6 NLS+NLW: Avrett stars 11 26 photon loss, no nanoflame heating, photon suction 7 NLS+NLW+NLT: Oslo simulations 13 38 spectral interlocking 1D HD RADYN simulation ...... 13 39 line haze modeling 2D MHD Stagger simulation ...... 13 13 detour and cross redistribution 3D MHD Bifrost simulations ...... 14 8 NLS – NLW – NLT chromosphere spectrum 14 27 coronal clouds Quiet chromosphere ...... 14 28 sophisticated clouds Active chromosphere ...... 15 44 RE upper photosphere Chromospheric and coronal heating ...... 16 33 classic abundance determination 9 Conclusion 16 42 LTE – NLTE schism arXiv:2103.02369v5 [astro-ph.SR] 25 Aug 2021 Acknowledgments 16 36 NLTE masking Version history 16 41 granulation contrast Notes 17 45 ALC7 line formation plots References 37 55 Ha–Ha scattering 49 Carlsson-Stein shocks Detailing notes, roughly grouped 50 Carlsson four-panel breakdown 58 internetwork shocks in Ca II H and Hα 1 personal background 52 bright-point shifts 1600–1700 Å 2 books 3 ADS page opening, silly publishers 60 non-E quiet chromosphere 61 internetwork fields and basal flux 5 solar telescopes 62 spicules-II 9 solar neutrinos 67 RBE–RRE–fibril ionization-recombination Compendium solar spectrum formation 2 1. Preamble: on-line format and on-line material

63 SDO quiet-chromosphere images The many links make this compendium Wikipedia-style mul- 60 transition region versus heated chromosphere tilayer: you may delve deeper by opening linked graphs, equa- 14 dark EUV features tions, references that suit your interest. In addition this is a 57 dynamic fibrils two-tier text with a relatively straightforward presentation in the 66 SolO campfires main text (contents above) but detailing many specific topics in 65 CBP foot visibility yellow-page endnotes4 starting at page 17 (list above). These add 69 EBs and FAFs further explanation including linked literature, comments, opin- 70 active-region moss 5 72 around and upward heating ions, and even history. 73 chromospheric FIP fractionation The links and two-tier format invite study at different levels: 71 IBIS–SDO active-region blinker 75 Lyα dream – read only the main text of Sect.3 and go through the corre- 76 Lyα features sponding SSF equation displays for a theory overview (as I 78 filament blobs did in my zoom lecturing). Readers only interested in opti- cally thin spectrum formation may stop already at the sep- 54 flash-spectrum color aration point on page8 ; for thicker readers the graphs of 40 Fe I 6302 Å polarimetry doublet 40 Mg I 4571 Å line formation in the didactic ALC7 star (Sect.6 ) starting 12 Ba II 4554 Å at [SSX (pdf 79)] are suited exercise material, with exam in 12 Mn I lines as activity monitor [SSF (pdf 113)]; 29 Na I D, Mg I b, Ca II 8542 Å heights – also read the main text of Sects.4–7 with SSX example in- 11 Ca II H & K core reversals spections for an overview of classic developments; 11 Wilson-Bappu effect 34 Mg I 12 micron emission lines – also study detailing endnotes and their sources on topics of 38 optical Fe I lines with S ≈ B interest to you; 38 weak Ce II and Fe II emission lines 54 He I D3 off-limb – read everything and open each linked page and study all 74 optical He I lines references for a full-fledged full-semester comprehensive 56 Rydberg H I lines course on the many intricacies of solar spectrum formation 77 Rydberg H I candidate for ALMA and some underlying solar physics, with this text serving as 48 quiet-Sun Hα scenes 55 quiet-Sun Hα – Ca II 8542 Å comparison convenient pointer and page and reference opener; 8 optical spectrum atlases and line lists – when no longer a student but a researcher (doctor?) and/or 35 plane-parallel atmosphere files teacher (saint?) this text may serve as page- and reference- 33 Kieler program for 1D LTE spectrum synthesis opening resource in taking your pick in Khayyam-wise re- 43 Pandora program for 1D NLTE spectrum synthesis freshing arguments you heard before. 31 MULTI program for 1D NLTE spectrum synthesis 32 MULTI3D program for 3D NLTE spectrum synthesis My on-line teaching material used here amounts to nine files 47 RH program for 1D-2D-3D NLTE–PRD spectrum synthesis listed and linked at the top of Table1 [1 personal background]. 59 Bifrost program for 3D radiative MHD simulation The first three are lecture notes. ISSF was a bachelor-level 18 CHIANTI package for CE analysis 46 my IDL programs introduction to Sac Peak summer students in 1993. IART was 64 SDO co-alignment programs my Utrecht course for second-year bachelors students, written 68 ALMA–GONG–SDO alignment in Dutch in the 1980s, translated by R.C. Peterson in 1992 and revived by L.H.M. Rouppe van der Voort in 2015. It follows Rybicki and Lightman (1986) with most emphasis on explain- 1. Preamble: on-line format and on-line material ing their first chapter summarizing basic RT up to NLS treat- This is a compact graduate-level course on solar spectrum for- ment. RTSA was initially written in 1995 for my Utrecht masters mation.1 It is display-oriented adding many figures to basic course, intended as an easier-to-read rehash of√ Mihalas (1970) equations but for brevity and to avoid repetition this text con- [2 textbooks] with emphasis on solar NLS ε scattering and tains almost no figures and equations itself: instead it links to adding NLW in particular for solar bound-free continua. The webposted teaching material and ADS-accessible literature.2 I 2003 update (still incomplete) is linked at ADS and therefore therefore recommend on-screen reading with web access.3 most frequently used – but at Oslo and elsewhere IART is pre- ferred as broader course material. None of these courses includes 1 Except high-energy spectra, radio spectra, polarization spectra, parti- NLT nor multi-level detours while these are now appreciated as cle spectra. For high-energy radiation start with Rybicki and Lightman key ingredients of chromospheric spectrum formation. I included (1986). For radio diagnostics start with Gary and Hurford (2004); they were treated by P.K.Manoharan in this school. For polarization (treated them in more recent summary tutorials (Rutten 2017a, Rutten by Rohan Eugene Louis and Debi Prasad Choudhary) del Toro Ini- 2019) and do so here. esta (2007) is a good introduction, Landi Degl’Innocenti and Landolfi 4 Linked with maroon bold-italics topic definers. Each endnote ends (2004) the comprehensive bible. For high-energy particles see Reames with a return link [Main call] to avoid getting marooned there and a (2021). Neutrinos are largely done (endnote9 ). 2 link [Back] that goes there also or back to another link that you used – Weblinks are blue, internal crosslinks are maroon. Acrobat may need working in Firefox, Ubuntu acroread and xpdf, macOS Preview but not web access in Preferences > Trust Manager. The suggestion is to open in Ubuntu evince, chrome, macOS Safari (try the back button instead). blue links (best in a browser parallel to your pdf reader) till you are blue Footnotes are just that (my trademark I’m told). Their number jumps in the face – meaning lack of sub-skin scattering (’tis all scattering here) back to their call. by blood drained far from your brain. 5 3 Naturally a personal selection favoring topics that I have been in- An off-line or slow-connection alternative is to download my files volved in and references and figures I knew best and came first to mind first, open them in parallel, and use the pdf page number specified in in this quick writeup for the Calicut school. every link. This may especially suit SSX which is the largest – but is no good for the many other page and abstract openers. Compendium solar spectrum formation 3 2. Introduction

Table 1. Acronyms and abbreviations, roughly grouped. Here I refer and link to these nine files for equations and ISSF = introduction to solar spectrum formation graphs that otherwise would swamp this text. In oral 3D bodily- IART = introduction to astrophysical radiative transfer on-the-spot teaching I would project and discuss these; here you RTSA = radiative transfer in stellar atmospheres should open and scrutinize them yourself. SSI = introduction to solar spectrum formation Navigation: links to these nine files (square-bracketed blue) SSF = solar spectrum formation theory should open the particular page in your browser. Similarly, blue SSX = solar spectrum formation examples links to ADS-accessible literature should open the cited page7 [3 ADS page opening, silly publishers]. The citations themselves SSA = stellar spectra A (Cannon – Payne – Minnaert) 8 SSB = stellar spectra B (Avrett – Chandrasekhar – Unsöld) should open the ADS abstract page. SSC = stellar spectra C ([Mihalas – Judge – Feautrier]) The many endnotes detail specific matters that in actual teaching I would treat in scheduled example lectures or as asides. RT = radiative transfer ∼ fermion-boson changes of Iν TE = thermodynamic equilibrium ∼ detailed balance in all LTE = local thermodynamic equilibrium ≡ SB populations 2. Introduction NLTE = non-local thermodynamic equilibrium ∼ SE non-E = non-equilibrium ∼ no SE Understanding solar spectrum formation is a necessity to under- R RE = radiative equilibrium ∼ αν [Sν − Jν] dν = 0 stand our star. The state of its matter is relatively simple: just SB = Saha-Boltzmann population equilibrium (→ Sν = Bν) gas, no fluids nor solids, even close to ideal. However, complex- SE = statistical equilibrium ∼ all populations stable ity is brought at wide-ranging spatiotemporal scales by dynam- CE = coronal equilibrium ∼ SE + thin collisional creation ics and electrodynamics with ionization and charging varying DEM = differential emission measure = CE density part j from none to full. The resulting structures and features may vary R ν EM = emission measure ≡ DEM(T) dT from utterly transparent to extremely opaque between different FIP = first ionization potential = neutral-atom ionization energy spectral samples. Using the emergent spectrum as diagnostic to NLS = non-local in space ∼ J important solar structure and solar happenings requires full elucidation of ν radiative interactions within the solar atmosphere. Observation, NLW = non-local in wavelength ∼ other Jν important NLT = non-local in time ∼ non-SE with memory theory and simulation constitute the three-fold way to enlighten- ment. Shortcuts as 1D, 2D, LTE, CE, SE may be misleading into CS = coherent scattering ∼ monofrequent, monochromatic CRD = complete redistribution ∼ resample extinction profile dead ends. PRD = partial redistribution ∼ mix CS and CRD Brief history: the solar spectrum caused the dawn of astro- physics in the nineteenth century in the hands of Fraunhofer, ∼ HD = hydrodynamics no magnetism Herschel, Kirchhoff and Bunsen [SSF (pdf 3)], and others. Fol- MHD = magnetohydrodynamics ∼ single-fluid HD + B lowed by stellar spectrum classification [SSF (pdf 6)] at Har- UV = ultraviolet e.g., AIA 1600 & 1700 Å vard by many ladies [SSF (pdf 5)] but foremost Cannon [SSF EUV = extreme ultraviolet ∼ shorter AIA wavelengths (pdf 7)] – I invite you to now do practical SSA 1.1 (pdf 5)9 DOT = Dutch Open Telescope @ La Palma, 0.45 m before reading further. In the twentieth century stellar spec- SST = Swedish 1-m Solar Telescope @ La Palma, 1 m trum interpretation evolved from basic concepts to understand- DKIST = Daniel K. Inouye Solar Telescope @ Maui, 4 m ing with Schwarzschild, Russell, Milne, Eddington, Minnaert, SDO = Solar Dynamics Observatory @ space Chandrasekhar and others. Major breakthroughs were Payne’s AIA = Atmospheric Imaging Assembly @ SDO demonstration with the Saha law that Cannon’s classification HMI = Helioseismic & Magnetic Imager @ SDO represents temperature ordering [SSF (pdf 57)], the Pannekoek– IRIS = Interface Region Imaging Spectrograph @ space − ALMA = Atacama Large Millimeter Array @ Chajnantor Wildt–Chandrasekhar identification of H as major continuous opacity provider (Fig. 3 (pdf 10) of Chandrasekhar and Breen MC = small kilogauss magnetic concentration (“fluxtube”) 1946), and Grotrian’s establishing the outrageous high temper- RBE = rapid blue-shifted excursion in blue wing Hα ature of the corona (endnote 15). Solar spectrum interpretation RRE = rapid red-shifted excursion in red wing Hα CBP = coronal bright point ∼ small hot quiet-Sun loop group went from LTE (Unsöld, endnote 42) to NLTE with Menzel, EB = Ellerman bomb (or burst: footnote 92) Thomas, Athay, Hummer, Jefferies, Mihalas, Avrett and oth- FAF = flaring active-region fibril e.g., in AIA 1600 Å ers. Then started (M)HD simulation modeling by Nordlund, Schüssler, Stein, Carlsson, Steiner and others with spectral syn- thesis programs by Carlsson, Uitenbroek, Heinzel, Leenaarts, The next three files in Table1 are sets of projection displays Pereira and others, presently emphasizing the new frontier of that I use in teaching since the passing of the viewgraph era.6 non-E interpretation. ADS serves nearly all their works. Figure2 shows an example. SSI is bachelors level, SSF and SSX 7 With most pdf viewers, but some may instead shunt you to the first are masters–graduate level. These also contain newer aspects. I page. This may also happen for pay-walled publications even with a showed a selection in my workshop lectures. license (endnote3 ). For such I add the pdf page number for manually The final three linked files in Table1 are practicals. They date going to the cited page. back to the 1990s and therefore use IDL but are easily coded 8 Following these recipes. They make a classical AAS-style reference in Python. I recommend SSA and SSB, designed to accompany list superfluous except for the few non-ADS publications. I do add one IART, for fresh RT students to gain hands-on insights. The third but with titles, bibcodes for copy-paste and ADS links as in my solar is a practical for RTSA developed by former graduate students physicist publication lists, which with my more useful solar abstract that still awaits conversion into proper SSC format. collection are a byproduct of my wholesome bibtexing with my entire bibfile collection. No index here because your viewer offers search. 9 6 This equation compendium still supplies 90 viewgraph pages copy- This scissors-only (pdf 7) exercise from Zwaan is ing all 545 labeled equations of RTSA. The first few are empty. illuminating. Young kids often do better than solar colleagues by sort- ing spectra without prejudice. Mrs. Cannon’s “early – late” suggests counting lines as wrinkles or by element number. [Spoiler: answer.] Compendium solar spectrum formation 4 2. Introduction

Fig. 1. Snapshot course summary, October 2018, Shandong University, Weihai. First column: particle–photon processes. Second column: corresponding extinction profiles. Third column: process pairs. Last column: key equations. The first word is already wrong but nevertheless I publicly signed this painting.

BASIC QUANTITIES

Monochromatic emissivity dE j dV dt dν dΩ dI (s) = j (s) ds ν ≡ ν ν ν 3 1 1 1 2 1 1 1 units jν: erg cm− s− Hz− ster− Iν: erg cm− s− Hz− ster−

Monochromatic extinction coefficient dI σ n I ds dI α I ds dI κ ρ I ds ν ≡ − ν ν ν ≡ − ν ν ν ≡ − ν ν units: cm2 per particle (physics) cm2 per cm3 = per cm (RTSA) cm2 per gram (astronomy)

Monochromatic source function j jc + jl Sc + η Sl S j /α = j /κ ρ Stot = ν Stot = ν ν = ν ν ν η αl /αc ν ≡ ν ν ν ν ν α ν αc + αl 1 + η ν ≡ ν ν P ν ν ν ν thick: (αν, Sν) more independentP than (αν, jν) stimulated emission negatively into αν, κν

Transport equation with τν as optical thickness along the beam dI dI D dI ν = j α I ν = S I dτ α ds τ (D) = α ds ν = S I ds ν− ν ν α ds ν− ν ν ≡ ν ν ν dτ ν− ν ν Z0 ν

Plane-parallel transport equation with τν as radial optical depth and µ as viewing angle z0 dIν dτν αν dz τν(z0) = αν dz µ cos θ µ = Iν Sν ≡ − − ≡ dτν − Z∞ start all talks index

Fig. 2. Screenshot of [SSF (pdf 23)]. SSFo ffers 119 similarly-styled theory displays (with a negative evaluation in [SSF (pdf 120)]). SSXo ffers 185 example displays. SSIo ffers only 16 introductory displays – you might start with these. The lightblue button start at bottom left should open the contents overview, index at bottom right a thumbnail index. The display title may return to the previously shown one. For me most figures open full-page when I click on them and many blue linkers open figures, images and movies elsewhere in my laptop – but not for you. Compendium solar spectrum formation 5 3. RT basics

This text gives an overview dividing radiative transfer (RT) 3. RT basics complexities between nonlocal in space (NLS), nonlocal in The figures on the preceding page show didactic overviews of wavelength (NLW), nonlocal in time (NLT). Solar physicists 13 generally appreciate NLS better than NLW and NLT least. this section. Table1 specifies acronyms I bombard you with .A NLS means that radiation received from the location we are summary of main RT quantities and equations is in the [RTSA studying is influenced by local radiation there that came from RT Rap (pdf 275)] and [SSF (pdf 107)]. elsewhere. It is daily familiar to us since our daytime outdoors is NLS illuminated. Whether sunshine or overcast, all photons we Intensity. The basic RT quantity is specific intensity Iν describ- see were made in the Sun and made it into our eyes by scattering, ing energy in a beam of photons at a given location at a given mostly multiple. As NLS as it can get! Very much out of LTE time at a given frequency in a given direction, measured per since the radiation temperature is about 6000 K, higher than our units of area, time, frequency bandwidth and beam spreading local temperature[ 4 spectral bandwidth and vision]. [ISSF Sect. 3.1 (pdf 9)] [IART Sect. 2.1 (pdf 22)] [RTSA Eq. 2.1 NLW means that radiation received from the location we (pdf 29)]. The last is hardest, [ISSF Fig. 2 (pdf 8)] may help. It are studying is influenced by radiation at other wavelengths at ensures that intensity does not vary with travel distance [SSF that location. We do not suffer NLW because sunshine scatter- (pdf 25)], an ingenious macroscopic way to express that photons ing around us is monochromatic (monofrequent): every detected are bosons that do not decay and may travel through the whole photon still has the energy (frequency, color) with which it left universe without change and convey the whole universe to us for the Sun. All solar Fraunhofer lines also reach our retinas unmo- inspection. lested (but disk-averaged and unnoticed10). We see the sky blue Common sense suggests that the Sun is brighter than a dis- and the setting Sun red not because of underway color change tant star, but that is irradiance which is flux [SSF (pdf 24)] but because Rayleigh scattering off molecules has higher prob- [RTSA Eq. 2.4 (pdf 30)] that is not measured per unit of beam ability at shorter wavelength. The green flash11 combines that spreading (solid angle [RTSA Fig. 2.1 (pdf 30)]). If the distant with larger refraction in the blue. star is solar-type it shows the same intensity, also when your NLT means that radiation received from the location we are time-travel starship parks you just above its surface.14 If your gi- studying is influenced by what happened there or thereabouts be- ant groundbased telescope resolves it the stellar intensity in the forehand. We suffer outdoors NLT because our sunshine is eight focus is the same as the solar intensity in the focus of photon- minutes retarded – but the same for all photons so we don’t care. starved DKIST[ 5 solar telescopes]. Solar NLS defines outward decline of line source functions12 At far-infrared, mm and radio wavelengths the Planck func- that is well exemplified by the Na I D lines and most exemplified tion simplifies to the Rayleigh-Jeans limit [SSI (pdf 3)] so that by Lyα, but it also brightens ultraviolet continua in and from the the intensity of an optically thick object with an LTE source solar atmosphere [SSF (pdf 109)]. More below. function equals its temperature – if calibrated into absolute units. Solar NLW is exemplified by most atomic lines having opac- In the optical and shorter-wavelength domains the Wien non- ity deficiencies up to an order of magnitude from sensing the linearity implies that intensity contrast comparisons between dif- scattering ultraviolet continua in ionization. This is less com- ferent wavelengths in terms of temperature need conversion into monly appreciated, for example not in profile-fitting “inversion” brightness temperature[ 6 equivalent temperatures] which re- codes nor in spectral irradiance modeling. More below. √ quires absolute units[ 7 intensity calibration] as in the optical Solar NLT is exemplified by solar Hα which is a NLS ε spectrum atlas of Neckel (1999)[ 8 optical spectrum atlases and scattering line but it is also NLW because Lyα defines its extinc- line lists] and the ultraviolet spectrum atlas of Curdt et al. (2001). tion (opacity) and a loop ionizing per scattering Balmer contin- uum from n=2 with cascade recombination back to n=2 includ- 15 ing Hα photon losses contributes to its source function. It is also Matter–radiation interactions. Photons live forever after their creation unless they meet fermions and interact. Also their cre- NLT because this loop is controlled by the n2 population which lags badly in gas that was first heated and then cools, as hap- ation is fermion business. For atoms, ions and molecules the pens continuously in the dynamic structures constituting the so- photon-involving processes are bound-bound transitions produc- lar chromosphere and even in long-lived filaments/prominences ing spectral lines and bound-free and free-free transitions pro- (endnote 78). Solar NLT can reach many orders of magnitude but ducing continua. In addition there are Thomson scattering of is ignored in almost all modeling by invoking statistical equilib- existing photons by free electrons and Rayleigh scattering of rium (SE). More below. existing photons by bound electrons. These processes can pro- I summarize the theory in Sect.3 , followed by classic de- vide local emissivity or extinction changing the intensity in velopments in order of complexity in Sects.4 –7. The extremes a given beam, in particular the intensity along the “line-of- of solar spectrum formation are easiest: LTE in the deep photo- sight” towards our telescope. Other photon-producing or photon- sphere, CE in the corona. In between is the hardest nut to crack: affecting processes may be ignored for the solar atmosphere un- the NLS – NLW – NLT spectrum of the chromosphere discussed less you are a radio astronomer (cyclotron, synchrotron, plasma in Sect.8 . radiation) or high-energy astronomer (pair annihilation). Photons being bosons means that they like to sit in the same 10 Even when you look at your beloveds although they get so richly place and therefore that any interaction producing or scattering adorned: 27500 colorful lines (not wrinkles) visible on their face! 11 13 My most beautiful√ was from the Neemach Mata temple in Udaipur. It would be good to hover your cursor over an acronym to pop up its 12 I call this ε scattering after the blue-bannered equation in Fig.1 meaning. I have tried latex command pdftooltip in pdfcomment.sty for which results for isothermal constant-ε atmospheres from the 2-level such popups but they worked too different between pdf viewers. simplifications of Eq.2 and Eq.3 below. More in endnote 20. 14 If the universe were infinite, eternal and homogeneously filled with stars our night sky would be solar bright everywhere, also in neutrinos. The fraction of astronomers studying this text would also be higher. 15 Heeck (2013) concludes (pdf 4) that they live at least 3 years – in their rest frame, meaning 1018 years for us. Compendium solar spectrum formation 6 3. RT basics two−level pairs detour pairs

2
2 2 u






extinction



1
1 1 l a b c h

2




















2
2



















2



















u u









emissivity





























































1 1



















1



















1 l l d e f g i j

Fig. 3. Atomic transitions governing line formation arranged in photon-involving pairs. The intensity in the beam of interest travels from left to right. All photons drawn here are at the line frequency = jump energy. The pairing categorizes what may happen to an atom after excitation −8 −2 (at left in each pair). It won’t stay excited long, for a 2-level atom only between 1/Aul ≈ 10 s (resonance) and 1/Aul ≈ 10 s (forbidden) until b, d or f occurs spontaneously, yet shorter if the atom while excited meets a deexciting collider a or an inviting akin photon to join on its way (c, e, g). Multi-level detour paths (righthand cartoons) combine other-wavelength transitions involving other levels and may include analogous bound–free transitions. There may be radiative or collisional deexcitation to another lower level or radiative or collisional excitation to a higher level including ionization and recombination. Other transitions eventually return the atom to the start-off level l. There may be endless up-down scattering in between so the detour possibilities are infinite. The upper row shows pair combinations contributing to line extinction: collisional photon destruction (a), scattering out of the beam (b and c), photon conversion out of the beam (h, into other-wavelength photons and/or kinetic energy). The lower row shows pairs contributing to line emissivity: collisional photon creation (d and e), scattering into the beam (f and g), detour photon production into the beam (i and j). Pairs c and g have equal probability by requiring one photon in the beam and one with arbitrary direction. a photon has a stimulated addition in which the new or scat- Equilibria. [SSF (pdf 69)] summarizes bound-bound equilibria tered photon coherently accompanies similar passing radiation with the 2-level simplification at left in Fig.3 . It assumes that that enhances the interaction probability. This contribution scales any excitation is followed by de-excitation in the same transi- with exp (−hν/kT), negligible in the Wien limit [SSI (pdf 3)]. tion, ignoring multi-level NLW and so restricting RT to a single For scattering it has the same probability between extinction and wavelength (CS) or the narrow spectral band covering a single emissivity (pairs c and g in Fig.3 ). transition (CRD). TE = “thermodynamic equilibrium” describes a homoge- neous isothermal fully-enclosed gas where nothing ever happens Bound-bound, bound-free, free-free interactions. The five and where every type of excitation occurs just as frequently as bound-bound interactions are spontaneous photo-deexcitation, its de-exciting counterpart (“detailed balance”): as many col- photo-excitation, induced photo-deexcitation, collisional exci- lisions down as up, as many radiative transitions down as up, tation and collisional deexcitation measured with the Einstein with all rates defined by the temperature. Maybe this boring par- coefficients Aul, Blu, Bul, Clu, Cul. The radiative first three are adise occurs approximately at the center of a black hole, but not transition properties that do not depend on circumstances (by within the Sun since non-enclosed implies net outward energy keeping the local photon supply outside the B definitions). Ein- leak from the hot core to the cold universe and the accompany- stein therefore derived their ratios assuming strict detailed bal- ing outward temperature decline through the entire Sun making ance (valid only in ideal TE) and declared these valid anywhere the fusion-produced energy surplus eventually escape[ 9 solar A f so that ul (or the classical oscillator strength ) is the single neutrinos]. radiative probability parameter [SSF (pdf 70)], [ISSF Sect. 5.1 SE = “statistical equilibrium” foregoes detailed balancing (pdf 23)], [IART Chapt. 5 (pdf 60)], [RTSA Sect. 2.3.1 (pdf 39)]. but requires overall balancing, meaning that all level populations The same five hold for bound-free ionization and recombi- do not vary with time however they are made. For each level the nation with appropriate nomenclature, and the same five hold rate equation [SSF (pdf 75)] [RTSA Eq. 2.100 (pdf 52)] govern- for free-free interactions (but collisional three-body up or down ing changes to its population ends up zero. Excesses between transitions have no radiative interest). going up or down some particular way are compensated by oth- ers going down or up. Simplest example: CE = only pair d, all Pair combinations. Combining the above interactions in suc- ups collisional, all downs spontaneous radiative for lack of im- cessive pairs is so important that I do not refer to Fig. 1 (pdf 3) pinging particles and photons. The SE assumption underlies all of Rutten (2019) but copy it here in Fig.3 . It came from [SSF following equilibria except non-E. (pdf 105)] and is an extension of [RTSA Fig. 3.3 (pdf 85)] with LTE = “local thermodynamic equilibrium” assumes TE at multi-level detours. Ca II and Hα detour examples are shown in the local temperature while permitting temperature and density Fig. 2 (pdf 4) of Rutten (2019). With these added this is a com- gradients and small leaks, clearly a cheap cheat. Its definition plete inventory of how line photons can be locally created and is that Saha-Boltzmann (SB) partitioning applies to all popula- modified in bound-bound manner.16 Below I refer to these pairs tions, as for element Schadeenium17 in my recommended SSA 2 by their alphabetic labels. Similar pair diagrams can be drawn Cecilia Payne exercise (pdf 9). It is often thought that LTE is for bound-free transitions (some in Fig.1 ). 17 Which you can’t pronounce. My Utrecht colleague Aert Schadee in- 16 In 3D teaching I draw these cartoons on the whiteboard as in Fig.1 vented this didactic element. Utrecht Nobel-winner Tini Veltman called and point to them all the time. In the blackboard era I drew them with his revolutionary algebra solver “schoonschip” (shipshape) to annoy his chalk and threw pieces at inattentive students until hitting Minnaert. CERN colleagues on pronunciation. Compendium solar spectrum formation 7 3. RT basics defined as equality of the source function Sν to the Planck func- of Burgess (1964) shows the high-temperature importance for tion Bν but this is a corollary [RTSA Sect. 2.3.2 (pdf 43)] [RTSA He I ionization. The radiative rate drops as 1/T but the dielec- Eq. 2.7.3 (pdf 45)]. LTE holds closely within the Sun and to some tronic rate first peaks two-orders-of-magnitude higher before its extent for many lines and continua in the low photosphere. The 1/T decay sets in. Classical Fig. 3 (pdf 20) of Jordan (1969)18 condition is that collisions both up and down heavily outweigh shows its importance for coronal iron partitioning (the lower photons up and down in Fig.3 and its bound-free counterparts. panel shows partitioning without it from House 1964). The “Jor- Photon destruction a and photon creation d + e must dominate dan versus Payne” CE–SB comparisons on [SSX (pdf 174)] in- over spontaneous and stimulated scattering, also for all steps in clude it also. A nice summary following Burgess (1964) is given the photon conversion sequences at right. Thus, the photons must in Sect. 6.2.5 (pdf 144) of Jefferies (1968)[ 11 Ca II H & K core behave as “honorary gas particles” (Castor), likewise boxed in reversals] to local circumstances; what we observe must be only a small Non-E = “non equilibrium” means non-SE = NLT, admit- leak. In SB partitioning the level population ratios sense only ting temporal variations in populations, i.e., having the SE pop- the temperature [RTSA Eq. 2.86 (pdf 49)] because the dominat- ulation equations [RTSA Eq. 2.100 (pdf 52)] not equaling zero ing up and down collisions both require one collider so that their but a function of time. For hydrogen it is often called “non- density drops out. The ionization stage ratios sense temperature equilibrium ionization” but the actual agent is NLT slow colli- more complicatedly by folding the Boltzmann factor with the sional settling of Lyα in cooling gas whereas H ionization occurs Maxwell distribution for the caught electron into the Saha dis- in an SE Balmer loop (more below). tribution [RTSA Eq. 2.88 (pdf 50)]; they scale inversely with the electron density from the need to catch one. NLTE = “non-local thermodynamic equilibrium” is a mis- Bound-bound redistribution. Eddington asked on page 2 (pdf 2) nomer. Usually it means assuming SE with accounting for NLS, of Eddington (1929) the “crucial question” to what extent a less often NLW, much less often NLT. resonance-scattered photon remembers the frequency it had in b, c, f, g CE = “coronal equilibrium” assumes SE but low electron the preceding photo-excitation (pairs ). The one extreme density and absence of local irradiation so that only d and its is that it does precisely in “coherent” scattering (CS). In Ed- dington’s days this was thought to be the rule, but in his land- bound-free counterpart remain. Every photon is locally created 19 by a collision obeying the local kinetic temperature. They occur mark thesis Houtgast (1942) showed that most solar lines in infrequently but in their ground state ions can wait endlessly for the visible obey the other extreme, suffering complete redistribu- suited colliders. After eventual excitation no collider or stimulat- tion (CRD) meaning complete loss of memory and representing a new sample of the line extinction profile rather than δ-function ing photon arrives within the 1/Aul decay probability for spon- taneous deexcitation. The resulting photon escapes outward to frequency conservation. CRD is indeed valid for most lines ex- space and may enter our telescope (fat chance), or it drowns in- cept the strongest (large extinction) whose photons escape high ward in the Sun, or it gets bound-free scattered in other gas (more in the atmosphere where collisions governing collisional redis- below). In this case the level ratios sense both temperature and tribution (“damping”) are rare. These lines must be described collider density since there is only the required up collision with- by partial redistribution (PRD) which combines Doppler re- out frequent down collision required for LTE. The stage ratios distribution in the core (CS in the frame of the atom but an sense only temperature because each ionization is collisional and redistributing average over the observed individually-Doppler- each recombination is radiative, both requiring a single electron shifted photon ensemble) with coherency in the inner wings and catch. CE and SB are thin versus thick partitioning extremes. collisional redistribution in deeply-formed hence collision-rich [SSX (pdf 174)] compares them for photospheric electron den- outer wings [SSF (pdf 89)]. The principal example is Lyα. Other sity; the SB peaks shift left by about -0.05 in log(T) per tenfold well-known PRD lines are Mg II h & k and Ca II H & K ([SSF (pdf 92)], [SSF (pdf 93)]). My personal example is Ba II 4554 Å Ne reduction. A complication and a boon is that dielectronic excitation [12 Ba II 4554 Å and Mn I lines]. Other strong lines suffer less must be included beyond the single-electron jump pairings in or no PRD[ 13 detour and cross redistribution]. Fig.3 . The reason is that at coronal temperatures the mean I must concede that PRD theory remains “not yet” in [RTSA Maxwellian energy is far above the desired standstill value for Sect. 3.4.3 (pdf 92)]. Classic: Hummer (1962) [SSF (pdf 90)]. the free electron to be caught at ionization threshold in one- Clearest explanation so far: Chapter 5 (pdf 110) of Jefferies electron recombination (peak probability in hydrogenic Kramers (1968) summarized in [SSF (pdf 91)]. The introduction to bf −3 Sukhorukov and Leenaarts (2017) is a good literature overview. αν ∼ ν edge decay). The shift of the Maxwell peak to larger energy for higher temperature gives 1/T dependence to radia- Recent developments are listed in [SSF (pdf 94)]. Next to fre- tive recombination rates. With dielectronic excitation this too quency redistribution there is also angle redistribution to worry large thermal energy is cut down by using a large part for si- about for moving gas due to its Doppler anisotropy – what gas in multaneous bound-bound excitation of another electron in both the solar atmosphere doesn’t move? See Leenaarts et al. (2012b). dielectronic ionization up and dielectronic recombination down. They are multi-level additions behaving in CE as pair d: up only Bound-free redistribution. Bound-free transitions are similar to collisions, down only spontaneous photons[ 10 dielectronic ion- bound-bound transitions and may be described by the same rate ization and recombination]. The temperature-only sensitivity of CE stage partitioning is therefore maintained; the rate ratio 18 Carole Jordan is the CE counterpart of SB Cecilia Payne, moving to between direct and dielectronic depends only on temperature Oxford from Culham rather then to Cambridge USA from Cambridge (Maxwell peak shift). The complication is that many possible UK. My SSA exercises should become Annie Cannon – Cecilia Payne bound-bound transitions to very high levels must be included and – Carole Jordan, likely using CHIANTI, and I should also add a Walter evaluated with CHIANTI (endnote 18). The boon is that the di- Grotrian [SSX (pdf 6)] exercise to make you appreciate that during to- electronic variants produce photons in many lines including the tality you bask in ordinary photospheric sunlight – but NLS weak, NLW line-free, and polarized: your cleanest illumination ever (endnote 15). diagnostics in the AIA EUV passbands and help keep the corona 19 Scanned by A.V. Sukhorukov and put on ADS by me. cool through substantial radiative losses. Original Fig. 1 (pdf 4) Compendium solar spectrum formation 8 3. RT basics equations [RTSA Sect. 3.2.3 (pdf 68)]. Resonance scattering (ra- sity [SSF (pdf 23)], [ISSF Sect. 3.2 (pdf 9)], [IART Chapter 3 diative up and down in the same transition, pairs b, c, f, g) is (pdf 34)], [RTSA Sect. 2.1.2 (pdf 32)]. Extinction is measured as also the same except that the bound-free spectral extent is much cross-section (cm−1 = cm2/cm3) so it simply adds up for differ- wider and that recombination catches an electron as a fresh sam- ent processes acting at one wavelength: (multiple overlapping) ple of the kinetic energy distribution without remembering the lines and continua. A photospheric overview is in [SSF (pdf 68)] energy of the electron kicked out in preceding ionization, hence [16 extinction diagrams]. The ultraviolet bound-free edges of pairwise only memory of ionization to the threshold but not of Si I, Mg I, Fe I and Al I contribute the electrons needed for H− the excess above it. Therefore there is always CRD over the ion- by pairing relatively large abundance with low ionization energy −4 ization edge. [SSF (pdf 56)], with their joint abundance giving Ne ≈ 10 NH A particular case is bound-free scattering of EUV line pho- and keeping the gas near-neutral.22 tons imaged in narrow bands as by EIT, TRACE, AIA. They Line broadening and the Voigt function are summarized in may be extincted in bound-free photo-ionization of H I, He I or [SSF (pdf 64)] and detailed in [RTSA Sect.3.3˙ (pdf 72)]. He II, most likely followed by spontaneous photo-recombination The strongest disk-center optical lines are Ca II H & K and because the colliding-particle and suited-photon densities are the Balmer lines (Fig. 10 (pdf 18) of Rutten 2019; extinction too low for other bound-free equivalents than scattering pair b. comparisons in [SSF (pdf 60)] and [SSF (pdf 62)]). R Through CRD such a bound-free scattered photon not only gets a Optical thickness τν ≡ αν ds is summed extinction along random redirection but also a new wavelength, most likely near the beam, with τν = 1 the 1/e decay separator between thin and the ionization threshold (probability for H I given by Kramers thick (transparent or intransparent). This is a property of the gas ν−3 decay [RTSA Eq. 2.74 (pdf 45)]), hence far from the line. parcel (structure, cloud, layer, shell, object), not of a spectral Bound-free scattering so redirects radiation out of the line of feature. Effective thickness likewise measures transparency for sight and also out of the passband; it may locally darken the monofrequently scattering radiation [IART Sect. 7.3.3 (pdf 103)] narrow-band image (often misnamed “absorption”)[ 14 dark [RTSA Sect. 2.6.3 (pdf 58)]. EUV features]. A particular point is that photons contributed by stimulated emission are not counted positively as contribution to emissiv- ity but negatively as contribution to extinction23. It facilitates Free-free redistribution. Every collider is a fresh sample of the the description because extinction (a + b + c + h) and stimu- kinetic energy distribution giving redistribution over the proba- lated emission (e + g + j) both need beam photons and therefore bility ν−3 decay [RTSA Eq. 2.76 (pdf 46)], a continuum without similar quantification (Einstein B values [SSF (pdf 70)]). More- threshold. In addition there is no memory of a preceding transi- over, with this subtraction the cancelation of stimulated scatter- tion so that pair combinations as in Fig.3 make no sense. If the ing pairs c and g is effectuated. If this isn’t done the extinction kinetic energy distribution is the Maxwell distribution defined would be overestimated by including c and optical depth scales by temperature then that fully defines both emissivity (Brehm- in the Rayleigh-Jeans domain would suggest larger opaqueness strahlung) and extinction so that their ratio (source function) and higher escape than actually the case.24 is the Planck function. Thus, free-free transitions always have For coronal spectrum formation one usually assumes CE S = B because they are intrinsically collisional. ν ν throughout which implies that all extinction is negligible and A particular case is ALMA for which at longer wavelengths emissivities can be evaluated without RT. It then suffices to com- the main agent is hydrogen free-free (electrons meeting protons) pute and sum all emissivities along the line of sight per wave- [SSX (pdf 32)]. ALMA is therefore often appreciated as a ther- length[ 17 emission measure, DEM]. This is done with the CHI- mometer, but it is as often not appreciated that the extinction is ANTI package[ 18 CHIANTI]. likely NLT controlled and may be orders of magnitude above the local LTE value, just as for Hα. ALMA is a thermometer but without easy knowledge where you stick it (more below). Thin – thick separation point. If your problem is so optically thin in all its aspects that assuming CE suffices you may stop reading here. Endnotes 10, 14, 15, 17, 18, 26, 27 may cater to Thomson and Rayleigh scattering. In terms of Fig.3 these you (but beware of my scorn at the end of 34). scatterings are described by pairs b, c and f, g, none of the oth- ers. In the Thomson limit (short of Compton) they still suffer Doppler redistribution. Pairs c and g cancel in occurrence prob- Source function. I wonder who gave this key non-thin quantity ability, averaged for Doppler redistribution. its well-suited name in describing weighted local addition of new A particular case is Thomson scattering by free electrons photons to a given beam by dividing emissivity and extinction making the white-light corona so beautifully visible during [SSF (pdf 23)]. It seems like dividing cows and horses since their eclipses and with LASCO[ 15 white-light corona]. dimensions are [erg cm−3 s−1 Hz−1 sr−1] versus just [cm−1], but the resulting dimension is that of intensity [SSF (pdf 23)] which ’tis all about. Extinction and emissivity. The extinction coefficient20 α and ν The source function is a key quantity for optically thick col- the emissivity21 j are the macroscopic ways of quantifying the ν leagues (mostly studying the photosphere and chromosphere) above processes locally along a given beam with given inten- because it separates local physics and surround physics, or 20 Naming: often called “absorption coefficient” but I follow Zwaan in and I realized that that is a better name because it is not a coefficient, using “extinction” because the photon is taken out of the particular beam not a fraction of something but measuring new radiation. (direction and frequency) but may still exist in another direction (scat- 22 Gaseous matter is mostly plasma throughout the universe but in cool- tering pairs b and c) or at another frequency (conversion pair h). Ab- star photospheres almost all electrons and in planetary atmospheres all sorption is then only pair a, sometimes called “true absorption”. 21 electrons are tied to nuclei. We live in an extraordinary neutral environ- Naming: Zwaan, Jefferies, Mihalas and others called this “emission ment energized by an extraordinary near-neutral environment. coefficient”. I followed them in IART and also while writing RTSA in 23 I wonder whether this was first done by Einstein himself. 1995, but at an Oslo school that summer P.G.Judge used “emissivity” 24 More in this R.N. Thomas memorial. Compendium solar spectrum formation 9 3. RT basics atomic perspective and environmental quality. Take a strong (pdf 128)][ 20 square root√ epsilon law]. Isothermal and 2-level spectral line, meaning large line extinction coefficient. Its local seems√ far-fetched but the ε law is often applicable[ 21 validity addition of many bound-bound transitions to continuous transi- of ε law] and yet more often instructive. An equation summary tions at its wavelength adds a high narrow spike in extinction ver- is shown in [SSF (pdf 79)] and as [RTSA concluding “RT rap” sus wavelength. This is local physics: the required particles must (pdf 275)]. be present in sufficient number set by elemental abundance, par- Two-level scattering source function evaluation became a titioning over ionization stages and atomic levels within stages key part of numerical spectrum synthesis with many programs that is primarily set by temperature and density (closer to SB following the key concepts of opposite-direction averaging of than to CE), and transition probability. In the emissivity there is Feautrier (1964) [RTSA Sect. 5.2 (pdf 137)] and of operator a similar spike given by lookalike particles in the upper level. splitting of Cannon (1973), leading to approximate lambda it- In the ratio, however, this spike may vanish completely – as for eration (ALI) and similar methods [RTSA Sect. 5.3.2 (pdf 145)]. LTE with the featureless Planck function as ratio. Then the ra- Some of these programs are detailed below in Sects.4 –6. 25 tio is completely set by the temperature, an environmental rather The extension with multi-level detours is still lacking in my than an atomic physics quantity. When LTE is not valid the local courses, so I add that here. The best descriptions so far are in environment also senses nonlocalness, surely in space (NLS), Sect. 8.1 (pdf 199) of Jefferies (1968)26 and in Sect. II (pdf 1) of likely also in wavelength (NLW), possibly also in time (NLT). Canfield (1971a), both using ratios rather than fractions as prob- Thus, optically thick colleagues evaluate local atomic properties ability parameters. Split the line extinction coefficient into the to define the extinction which defines where their signal comes destruction (a for absorption), scattering (s), and detour (d) con- from, then evaluate the source function to see what environmen- tributions of Fig.3 : tal effects may make it depart from the local temperature. This αl ≡ αa + αs + αd ε ≡ αa/αl η ≡ αd/αl (1) split is often described as “accounting for NLTE”. ν ν ν ν ν ν ν ν ν ν However, such split interpretation gets mixed when envi- where ε is the collisional destruction probability of an extincted ronmental effects also affect the extinction. This happens when photon and η is its detour conversion probability. With these, the that is controlled by other transitions sensitive to source func- general line source function becomes tion nonlocalness. The extinction of optical atomic lines being Sl = (1 − ε − η ) J + ε B (T) + η Sd (2) defined by scattering bound-free ultraviolet continua is one ex- ν ν ν ν ν ν ν ν ample, the dominance of Lyα in controlling Hα extinction and for CS and ALMA H I free-free extinction is another. In such cases, opac- S l − ε − η J ε B T η S d ν0 = (1 ν0 ν0 ) ν0 + ν0 ν0 ( ) + ν0 ν0 (3) ity NLTE gets more important than source function NLTE. Col- RR leagues that see NLTE only as a source function issue may miss for CRD with Jν0 ≡ (1/4π) Iνϕ(ν−ν0) dν dΩ the “mean mean” the crux of their problem. intensity averaged over all directions and the line profile, with ν0 Also, the source function is a composite because it is the line-center frequency and also used as line identifier. a process property. Multiple processes operating at a given Equations2 and3 describe physical contributions to local wavelength each have their own source function and combine NLTE divergence between S l and B. In each the first term repre- weighted [SSF (pdf 23)]. In the solar photosphere optical lines sents the reservoir of photons that contribute new photons to the often combine a scattering line source function with a near-LTE beam by direct scattering (f and g), the second describes colli- continuous source function because the latter is dominated by H− sional beam-photon creation (d and e), the third the contribution bound-free interactions. These are not necessarily LTE as is the of new beam photons via detours (i and j) which themselves may case for free-free interactions. Ultraviolet bound-free continua contain collisional, scattering and also such bound-free steps. do get far from LTE but in the optical and infrared the ambi- A convenient formal way to quantify S l , B inequality is ent thermal kinetic energy is similar to or exceeds the transition to use NLTE population departure coefficients b ≡ n/nLTE for energy and at sufficiently deep formation collisions reign, mak- the lower (bl) and upper (bu) level[ 22 Menzel departure coeffi- ing LTE a good approximation[ 19 deepest wavelengths]. The cients]. In the Wien approximation (negligible stimulated emis- Na I D lines reach deepest in the√ optical spectrum (endnote 29) sion) these yield the simple relations [SSF (pdf 73)]: because they are photospheric ε scatterers (endnote 29). In hot- l l l α ≈ bl α S ≈ (bu/bl) Bν (T). (4) ter stars where hydrogen is ionized the reverse occurs: weak LTE ν0 0 LTE lines appear in emission on the dark Thomson-scattering The formation graphs for many lines from ALC7 starting at continuum (Fig. 1 (pdf 215) of Rybicki and Hummer 1992). The [SSX (pdf 79)] employ these (endnote 45). solar Ca II H & K wings also contain weak emission lines but Equations2 and3 show that LTE S = B equality holds when these are from NLTE interlocked scattering (endnote 38) with ε = 1, η = 0 and/or J = S d = B, both closely correct below the c Fe II 3969.4 Å close to H the weirdest: it responds deeper standard h = 0 surface at τ5000 = 1. Above it ε becomes small than its quasi-continuum near-LTE background (Fig. 8 (pdf 9) of from lower electron density[ 23 epsilon is small] while η is usu- Cram et al. 1980). ally smaller[ 24 no Thomas photoelectric control]. Most chro- mospheric lines and bound–free continua are heavily scattering with S ≈ J (but free–free continua always have S = B because Line source function. Most classic treatments of NLTE scat- each interaction is collisional). For such scatterers the two-level tering in the line source function assume a gas of 2-level atoms simplification η =0 is a good approximation but weaker spectral to keep the description monochromatic (CS) or limited to a sin- features may suffer detour complications[ 25 photon pumping, gle line (CRD). I did the same in RTSA, with a long deriva- lasering, Wien transcription][ 26 photon loss, no nanoflame tion inventing “sharp-line atoms” enabling monofrequent use of heating, photon suction]. the Einstein coefficients in [RTSA√ Sect. 3.4 (pdf 84)], an ele- gant derivation by Zwaan of the ε law in [RTSA Sect. 4.3.1 25 An online tool for hands-on experimentation with various methods is (pdf 112)], and the classic numerical CRD results for this law in at http://rttools.irap.omp.eu and described by Lambert et al. (2016). It an isothermal atmosphere of Avrett (1965) in [RTSA Fig. 4.12 produces plots as [RTSA Fig. 5.2 (pdf 147)] on your screen. 26 Erratum: remove minus in equation below Eq. 8.8. Compendium solar spectrum formation 10 4. LTE: Holweger models

Transfer equation along the beam. See [ISSF Sect. 3.3 Be aware that use of the Eddington-Barbier approximation is of- (pdf 10)] for a gentle introduction. Here τν is optical thickness ten a fallacy[ 30 contribution, response, sensitivity]. measured in the beam direction. [ISSF Fig. 5 (pdf 14)] treats a Be also aware that solar-atmosphere structures may be opti- spectral line from a homogeneous isothermal cloud with LTE cally thick without being effectively thick so that radiation scat- processes, say a galactic H I cloud showing the 21-cm line to a ters through. Filaments or fibrils may appear intransparent in radio astronomer. Of my cartoon diagrams this one is most of- Hα while illumination from behind is still important. This is ten shown, nowadays as Python animation. Appreciate how the even the case for Hα in the “chromosphere” plateau of 1D stan- line grows or sinks with increasing cloud opacity until it hits the dard models, as in [SSX (pdf 91)] which not only shows S ≈ Iν = Sν = Bν optically thick limit (more formal graph in [IART J > B across the model temperature minimum from backscat- Fig. 3.3 (pdf 42)] and [ISSF Fig. 4 (pdf 13)]). Thick isothermal tering (endnote 24) but also has photospheric photons scatter- LTE clouds show no lines (as “black holes have no hair”). ing through the model chromosphere. This caveat was beauti- Such simple clouds[ 27 coronal clouds] became the basis fully demonstrated in Fig. 7 (pdf 7) of Leenaarts et al. (2012a) for “cloud modeling” started by Beckers (1964)27 for Hα mot- by comparing 1D Hα synthesis for a Bifrost (endnote 59) simu- tles on the disk and also much used for filaments. Assuming lation with MULTI[ 31 MULTI program] and 3D Hα synthesis the cloud homogeneous with Gaussian (“microturbulent”) line with MULTI3D[ 32 MULTI3D program]. In the columnar 1D broadening yields the classic cloud parameters [SSX (pdf 2)]: line synthesis the deep-photosphere granulation is better visible line-center optical thickness, line-of-sight velocity, source func- than lower-contrast chromospheric fibrils but in 3D synthesis the tion, Doppler width, incident intensity profile from below/behind granular contrast is smoothed away by lateral scattering. (review Tziotziou 2007)[ 28 sophisticated clouds]. The most important equation of classical 1D 2-level scatter- ing is the Schwarzschild equation defining the Λ operator [RTSA Eq. 4.14 (pdf 98)]; a synopsis is given after Eq. 10 (pdf 9) of Rut- Transfer equation in optical depth. For very optically thick ten (2019). Its cutoff at the surface produces outward J < S di- objects (as the Sun) there is no interest in what impinges the ν ν vergence for small inward increase of Sν(τν) but outward Jν > Sν other side since nothing comes through. We then flip τ to mea- divergence for steep increase [RTSA Fig. 4.2 (pdf 100)], [RTSA sure optical depth instead, either along the line of sight (τνµ) or Fig. 4.4 (pdf 103)], [RTSA Fig. 4.9 (pdf 121)]). Be aware that more often against height radially into the Sun (τν). The standard 28 steep horizontal gradients are similarly important in 3D radiative recipe is then to apply the Eddington-Barbier approximation transfer. Steep horizontal gradients occur for example already I+ ≈ S S ν (τν = 0, µ = 1) ν(τν = 1) with τν and ν the totals [ISSF in the photosphere in and above granulation and in and around Eq. 35 (pdf 13)], [IART 321 (pdf 45)], [RTSA Eq. 2.44 (pdf 38)], 29 MCs constituting network and plage. Obviously they occur in [RTSA Eq. 4.32 (pdf 106)], [SSF (pdf 29)], [SSF (pdf 51)]. and around the multitude of small dynamic structures constitut- The recipe is easiest for LTE in which emergent line pro- ing the higher atmosphere. files map temperature variation with height through Planck- [SSF (pdf 108)] gives a summary of the key line formation function sensitivity conversion and with extinction-coefficient equations. profile weighting [SSI (pdf 16)]. Absorption lines map outward declining source function, emission lines increasing source func- tion. Avrett sketched this sign change for solar lines in [SSX 4. LTE: Holweger models (pdf 36)]. The flip from absorption to emission near 1600 Å in [SSX (pdf 34)] was already described by de Jager (1963). The The classic Holweger models (Holweger 1967 and HOLMUL of Holweger and Müller 1974) were a pinnacle of solar LTE inter- similar flip at 160 µm is closer to LTE. 31 Absorption lines go together with limb darkening [SSF pretation. They replaced curve of growth abundance determi- (pdf 51)] [RTSA Fig. 2.4 (pdf 39)], emission lines with limb nation by individual line-by-line fitting which became an abun- brightening. Just beyond the limb the first flash-reverse into the dance industry[ 33 classic abundance determination]. It em- latter [RTSA Fig. 7.2 (pdf 168)].30 ployed critical ad-hoc adjustment parameters: fake line broad- The “four-panel” diagrams in [SSF (pdf 31)] ff are a cartoon ening for duplicating observed line broadening with “microtur- course on Eddington-Barbier spectrum formation with increas- bulence” [RTSA Eq. 3.80 (pdf 82)] (adding height-dependence ing complexity. You should blink the absorption – emission pairs as in Fig. 11 (pdf 19) of Vernazza et al. 1981) and “macro- by full-page flipping. Then do the exam on [SSF (pdf 40)], first turbulence” [RTSA Eq.3.81 (pdf 83)] (with “radial-tangential” anisotropy added by Gray 1977), and fake wing extension with a for Na I D2 (answer h ≈ 500 km) and then for the blend. If you decide h ≈ 150 km for the blend by reading off where its center collisional damping “enhancement factor”. The Holweger mod- intensity equals the source function just as for Na I D you are a els ignored the existence of the solar chromosphere to avoid 2 34 Mg I 12-micron emis- million times wrong[ 29 Na I D, Mg I b, Ca II 8542 Å heights]. self-reversals in strong optical lines[ sion lines], but nevertheless HOLMUL became highly popular32 27 J.M. Beckers (1934 – 2021) defended this thesis at Utrecht Univer- with abundance determiners because it worked so well, giving sity with M.G.J.Minnaert˙ as formal adviser but did the work in Sydney them smaller spread for multiple lines than other standard mod- with R. Giovanelli as effective adviser. At his request I put it on ADS els of the time[ 35 plane-parallel atmosphere files]. No wonder (Beckers 1964). It was printed by the US airforce because he had al- because the model was made as best fit to photospheric lines, ready moved to Sunspot where I got to know him – then a solar physics in particular of Fe I, so that fitting similar lines with the model paradise, now a ghost town (Rutten 2018). 28 was a self-fulfilling prophecy, in particular for the then much Unjustly named (Paletou 2018) – as often with Eddington. debated Fe abundance itself (e.g., Kostik et al. 1996). Earlier 29 ≈ When you see Iν Sν(τν = 2/3) you read text from an unresolved Rutten and Kostik (1982) showed that such model fitting is also nighttime colleague confusing flux and intensity. They also add a telltale self-correcting[ 36 NLTE masking]. factor π in their Bν definition [SSF (pdf 29)] [RTSA Eq. 2.7 (pdf 31)]. 30 Ca II H & K and the Balmer lines are also the strongest in the optical 31 Treated in practical SSA 3 (pdf 28). Fig. 4 (pdf 6) of Rutten and flash spectrum, but there out of the blue joined by He I D3 which is Zwaan (1983) may have been the last solar Fe I one. endnote 74 virtually absent on the disk ( ). 32 ADS shows 900+ citations but declining rate this century. Compendium solar spectrum formation 11 5. NLS: Auer & Mihalas modeling

So-called “inversion” codes (I find “best-fit” better suited) 5. NLS: Auer & Mihalas modeling are automated versions of Holweger’s trial-and-error model derivation procedure[ 37 Milne-Eddington approximation]. NLS is the realm of 2-level modeling treated extensively in Copying Holweger’s fads also copies the fallacies including [RTSA Chapt. 4 (pdf 95)], numerical developments around the NLTE masking, an important failure because NLTE scattering in Λ operator summarized in [RTSA Chapt. 5 (pdf 133)], and nu- the ultraviolet bound-free edges controls the opacity of minority- merical PRD formalisms (not yet in RTSA). [SSF (pdf 111)] species atomic lines throughout the spectrum. In the optical these shows a cartoon summary also shown as Fig. 4 (pdf 10) of Rutten often have source functions close to LTE through interlocking (2019) and discussed in Sect. 4 (pdf 11) there. Non-local scat- [38 interlocking Fe I, Ce II, Fe II lines] but their interlocked tering makes J depart from B in setting S according to Eq.2 opacities have NLTE deficits up to an order of magnitude from and Eq.3 , for bound-free continua and for lines of increasing ultraviolet overionization set by steep temperature gradients in strength (or different parts of strong PRD lines). The ultraviolet ≈ deeper layers. This fallacy has been long ignored but was at last continua gain S J excess over the photospheric B decline√ per addressed by Smitha et al. (2020). Worse: in the actual 3D so- Λ [SSF (pdf 86)]. Strong lines obtain low S per Λ from ε scat- lar photosphere matter and radiation gradients also occur and act tering by seeing actual temperature gradients as near-isothermal horizontally, requiring 3D line synthesis (see Smitha et al. 2021). from τ scale compression (endnote 21). Proper quantification of such ultraviolet extinction control is The classic warm-star H I analyses of Auer & Mihalas severely hampered by the dense ultraviolet line haze with its own (1969a, 1969b) advertised here were an early pinnacle of 2- NLTE properties[ 39 line haze modeling]. Stronger lines need level modeling[ 42 LTE – NLTE schism]. Not purely 2-level NLTE modeling not only for their extinction but also for their since adding the H I continuum level and so including the Ly- source function set by their own scattering or by non-2-level de- man and Balmer continua to Lyα in the first paper and to Hα in tour transitions. When “inversion” codes are proudly advertised the second. The key results were copied in Mihalas (1970) and as NLTE usually only the modeled-line resonance scattering is in [RTSA Fig. 7.10 (pdf 187)]. The great didactic value of their meant, not the extinction-affecting scattering in the bound-free results is that the effects of these transitions are well separated in continua nor resonance scattering in the line haze affecting those height. “Big whopper” [RTSA exercise12 (pdf 254)] poses five pages of hard questions on their graphs that I challenge you to [40 Fe I 6302 Å doublet, Mg I 4571 Å]. answer.34 The hardest is to understand: “Amazingly, the pho- Holweger’s assumptions of 1D modeling, LTE and no chro- ton losses in the subordinate Balmer-α line, located in the low- mosphere were also revived in SATIRE irradiance modeling for energy red part of the spectrum, cause heating of the whole outer network and plage (Unruh et al. 1999) which therefore suffers atmosphere of this hot star.” [RTSA (pdf 189)]. Because their similarly from mistreating the opacity-affecting NLTE ultravi- graphs did not include J we did a multi-level re-do including olet continua and their NLTE line haze. More explanation in Avrett-style S, B, J graphs (pdf 5) and radiative heating/cooling Sect. 4 (pdf 11) of Rutten (2019). graphs (pdf 6) in Wiersma et al. (2003), giving as answer: “The The Asplund revolution (Asplund et al. 2009)33 ended the photon losses in Baα suck population from the proton reservoir HOLMUL popularity, primarily by including simulated granu- through a collisionally-dominated Rydberg recombination flow lation (Nordlund et al. 2009 with the foundation laid in Nord- and so boost the outward temperature rise.” (pdf 7). All yours to lund 1984a)[ 41 granulation contrast]. Its main effect is that understand. its temperature inhomogeneity upsets the notion that spatially- averaged intensity may be modeled as spatially-averaged tem- perature because at optical and shorter wavelengths this notion is 6. NLS+NLW: Avrett stars undone by the Wien nonlinearity of the Planck function (Uiten- broek and Criscuoli 2011). The simulated granulation undid the The various quiet-Sun models masterminded by E.H. Avrett with need for turbulence fudging and better atomic physics (Barklem the[ 43 Pandora program] (VALIIM of Vernazza et al. 1976, et al. 2000) undid the need for damping fudging. VALIIIC of Vernazza et al. 1981, MACKKL of Maltby et al. Another Asplund ingredient was NLTE spectrum synthesis. (1986), FALC of Fontenla et al. 1993 and ALC7 of Avrett and However, yet without detailed spectral synthesis of ultraviolet Loeser (2008) with an update in Avrett et al. 2015) are all classic + overionization weakening atomic minority-species lines as re- pinnacles of plane-parallel (1D) NLTE modeling (2100 cita- 35 viewed for Fe I by Rutten (1988) and quantified from granu- tions for VALIIIC at increasing rate). Their strength is com- lar/intergranular 1D models by Shchukina and Trujillo Bueno prehensive inclusion of ultraviolet transitions. Their weakness is (2001). Such synthesis should be 3D to include also the steep ignoring Wien nonlinearity in assigning observed intensities to horizontal temperature gradients (Fig. 4 (pdf 6)) by which hot- averaged temperatures of inhomogeneities, the same mistake as ter granules irradiate cooler intergranular lanes in the ultravio- in HOLMUL abundance determination rectified in the hydrody- let where the Balmer continuum escapes. This irradiation also namic Asplund revolution. causes line weakening and therefore abundance underestimation NLW was still lacking in VALIIM which only accounted when neglected; correction may give a shift back towards classic for scattering in the Si I continuum so that its Fig. 23 (pdf 38) “high” values desired by helioseismologists. multi-page analogon to multi-page VALIII Fig. 36 shows erro- neous Sν ≈ Bν in the ultraviolet above the Si I threshold at 1682 Å 33 Likely at over 5000 citations the most cited publication of solar where other electron-donor bound-free contributions dominate physics. If you desire citations you must cater to non-solar colleagues; [SSX (pdf 68)]. This was remedied in VALIIIC but overestimat- my most-cited (and most boring) article (Carlsson et al. 1994) was non- ing their effect because not enough line-haze blending lines were solar: stellar lithium abundance determiners outnumber solar spectrum analysts. 34 [Spoiler: my answers to most RTSA problems including this one.] 35 Practical SSB 1 (pdf 5) dissects the FALC stratifications with com- parison to our own atmosphere which is transparent to sunlight by miss- ing H− treated in the continuum part of SSI and in practical SSB 2 (pdf 13). Practical SSB 3 (pdf 22) makes you synthesize the Na I D lines from FALC assuming LTE and hence misrepresenting their cores. Compendium solar spectrum formation 12 6. NLS+NLW: Avrett stars included, resulting in a too steep upper-photosphere temperature a double champion of solar physics: the most-shown figure but decline. Including more and more lines from the Kurucz tabu- also the most misinterpreted one. It misled many solar physi- lations (e.g., Kurucz 2009) then brought the upper-photosphere cists into thinking it describes the Sun instead of an unrealistic back up in MACKKL and since then it has remained the same. plane-parallel-layer star. The solar atmosphere is not layered but It is remarkably similar to HOLMUL and even to Kurucz or Up- dynamically structured. psala LTE–RE models [RTSA Fig. 7.3 (pdf 169)][ 44 RE upper Take Hα as example. In his VALIIIC diagram [SSX (pdf 60)] photosphere]. Avrett, admirably careful, drew its core coming from the VALI- I love the Avrett models, in particular VALIIIC because of IIC chromosphere, its wing from the low VALIIIC photosphere the many diagrams in Vernazza et al. (1981) and especially its with an extended gap in between that is not present in his forma- many-page Fig. 36 continuum formation diagrams (some copied tion spans of Ca II K and Mg II k. He did so because at low tem- in [RTSA (pdf 204)]– [RTSA (pdf 206)] and in [SSX (pdf 61)]– perature the Boltzmann excitation and therefore the extinction of [SSX (pdf 67)]). Since VALIIIC obeys all the standard RT equa- Hα are negligible[ 48 quiet-Sun Hα scenes]. Correct – but the tions, meaning all those in RTSA, these graphs describing spec- actual Hα core in actual quiet Sun mostly samples opaque fibril tral continua arising from the VALIIIC density and tempera- canopies reaching much higher. ture stratifications provide valuable insights on spectrum for- Worse, the VALIIIC chromosphere and similar Avrett-star mation and its intricacies. They are of enormous didactic value chromospheres do not represent average canopy-fibril tempera- because everything in Avrett’s output graphs is fully understand- tures but ultraviolet radiation temperatures reached in acoustic able from the input physics ≈ RTSA physics. The same holds for shocks in internetwork regions underneath Hα fibril canopies. the newer Avrett models and his latest ALC7. Therefore I made The shocks result from upward propagating waves excited by many ALC7 demonstration figures: [SSX (pdf 73)] showing the 3-min components of the photospheric p-mode interference ALC7 stratifications, [SSX (pdf 76)] ultraviolet continua, [SSX pattern (endnote 53). Their radiation temperatures are Wien- (pdf 77)] hydrogen line formation, [SSX (pdf 77)] overview of weighted to the shock maxima and are not linear tempera- strong-line formation, then [SSX (pdf 78)] ff explaining line- ture averages. This was shown in Fig. 4 (pdf 8) of Carlsson formation plots[ 45 ALC7 line formation plots] followed by and Stein (1994) where the upgoing waves are still wave-like such plots for thirteen ALC7 lines in [SSX (pdf 82)] ff. They in the Al I edge and fully shocking in the Lyman continuum. were all made[ 46 my IDL programs] with the 1D version of Their model-building experiment in Fig. 5 (pdf 11) (summary in github-publicRH[ 47 RH program]. They represent a detailed [SSX (pdf 157)]) demonstrated that the nonlinear Wien weight- view of how all these lines form in the ALC7 atmosphere. In ing skews the apparent average over internetwork shocks into an addition [SSX (pdf 98)] ff show ALC7 formation comparisons Avrett-type apparent chromosphere[ 49 Carlsson-Stein shocks] between comparable pairs of lines and [SSX (pdf 102)] ff com- [50 Carlsson four-panel breakdown]. pare the formation of twelve lines between three Avrett stars.36 In AIA 1700 and 1600 Å movies these acoustic internetwork Any student of thick solar spectrum formation should be able to waves appear as fast-moving fast-evolving wisps of emissivity understand and appreciate all these curves in all their gory detail. amid stable network evident as bright grains representing kilo- I invite you to study and explain: gauss “fluxtube” MCs. The 1600 Å images sample them slightly – Ne stratification; higher because the summed photo-ionization extinction of elec- – Eddington-Barbier validity in all lines; tron donors Mg I, Fe I, Si I, and Al I increases for shorter wave- – ε decline in scattering lines; length (compare emergence heights and main contributors on the successive pages of VALIII Fig. 36 from page 670 (pdf 37)). – PRD-split source functions of Ba II 4554, Ca II K, Mg II k, This increase produces the Fourier phase differences in the up- Lyα; per panel of Fig. 18 (pdf 18) of Krijger et al. (2001) reproduced – thermalization depth for all lines, in particular Hα; per RADYN simulation by Fossum and Carlsson (2005b)[ 51 – NLTE extinction dip of Mg I 4571, Fe I 6301.5, Mg I b2 but Acoustic and gravity wave heating]. not for Na I D1, He I 10830; The MCs outlining network appear with brighter contrast – LTE extinction of Mg II k, Hα,Hβ, He I 584; at 1600 Å then at 1700 Å. Plane-parallel colleagues attribute – slight rise in NLTE over-extinction for Lyα; this enhancement also to increased donor-ionization opacity and – line formation differences between different Avrett stars. standard “network” and “plage” models therefore have less steep upper-photosphere temperature decay than their “quiet” com- However – the stern caveat is that I love the Avrett models only panions for internetwork. However, the actual fluxtube bright- for representing educational solar-analog “Avrett stars” [SSX ness is photospheric hole-in-the-surface radiation, with deeper (pdf 58)]. These computationally existing stars extend infinitely holes at shorter wavelength. This is also known for decades; I plane-parallel (1D) without any inhomogeneities or any mag- put a review as caption in Fig. 73 (pdf 91) of Rutten (2020) and netism or any sort of waves or any type of reconnection – no refer to that[ 52 MC shifts 1600–1700 Å]. Actual MCs are close granules, spicules, floccules, let be filaments, flares, CMEs, and to RE throughout the photosphere without evidence of kinetic whatever else that makes our kind Sun non-plane-parallel inter- heating as expected for structures that inhibit convection (Fig. 7 esting. The Avrett stars are a marvelous boon to a RT teacher (pdf 7) of Sheminova et al. 2005). like me but they also had and have detrimental effects on far Thus, the standard-model “chromospheres” non-linearly rep- too many colleagues misled into accepting their stratifications as resent under-the-chromosphere internetwork shocks (a domain I a valid average over small fluctuations around realistic means. called “clapotisphere” in Fig. 12 (pdf 7) of Rutten 1995)[ 53 so- The first figure of Vernazza et al. (1981) in [SSX (pdf 60)] is lar p-mode history], plus misinterpreted network surface holes. 36 ALC7 – FALC in [SSX (pdf 102)] ff, FALC – FALP in [SSX Nothing to do with the on-disk chromosphere in Hα that car- (pdf 115)] ff, ALC7 – FALP in [SSX (pdf 128)] ff, each with twelve line ries the name originally given by Lockyer (1868)[ 54 Lockyer buttons at bottom. The display titles may return to the previously shown chromosphere, He I D3, flash color] for the colorful off-limb display. The start button at bottom-left should return to the SSX con- appearance of the Balmer lines and He I D3. tents overview. Compendium solar spectrum formation 13 7. NLS+NLW+NLT: Oslo simulations

Yet worse, actual Hα canopies in quiet areas are mostly made the impinging Balmer continuum originates in the deep photo- by spicules-II (endnote 62) that reach heights around 7000 km sphere, with radiation temperature about 5300 K [SSX (pdf 76)]. (Pereira et al. 2014) far beyond any Avrett atmosphere. On the Where this value exceeds the local temperature photon pumping disk they are observed as outer Hα-wing RBEs and RREs (end- enhances ionization. There are also contributing photon losses note 62) – so much for ascribing outer-wing intensities to deep- in heavily scattering Hα [SSX (pdf 91)] which operates within photosphere sampling and so much for assigning brightness to the low-energy 3.6 eV hydrogen top (footnote 60) representing heating since they heat to EUV visibility (Henriques et al. 2016) a minority-species photon-suction loop as in the neutral alkali’s but are dark. Their subsequent cooling recombining return flows (endnote 26).38 produce most of the dark Hα line-core fibrils (“coarse mottles”) The key hydrogen balance for NLT is the Lyα collision bal- around network (Rutten et al. 2019). Both types of ubiquitous ance. The collisional up and down rates differ much in temper- dynamic Hα structure are transparent in the ultraviolet continua ature sensitivity due to the Boltzmann ratio in their Einstein re- used for constructing standard models[ 55 Ha–Ha scattering]. lation [SSX (pdf 166)] [RTSA p. 23 (pdf 43)] [RTSA Sect. 3.2.5 Yet worse, there even exist publications (not to be named) (pdf 70)]. The exponential Clu sensitivity determines the settling that used the VALIIIC diagram in [SSX (pdf 60)] to define time scale at which Lyα reaches Boltzmann equilibrium. A few “height differences” as the 400 km there between Ca II K3 and cases are plotted in Fig. 6 (pdf 7) of Carlsson and Stein (2002), the core of Hα to identify and measure wave propagation and reproduced in [SSX (pdf 166)]. In hot shocks the frequent Lyα explain coronal heating. Maybe there exist plane-parallel stars up- and down collisions make the balance reach Boltzmann equi- in a parallel plain universe where this is a viable tactic (but do librium in seconds, but in subsequent post-shock cooling gas the those have coronae?). settling time scale increases to multiple or many minutes. During this time the large n2 overpopulation also governs large non-E overionization since the whole top of the hydrogen atom includ- 7. NLS+NLW+NLT: Oslo simulations ing the proton population is defined by this lower boundary con- NLT is the current frontier in solar spectrum interpretation, dition. The Balmer-top SE loop just adds yet more overioniza- not just in underlying MHD phenomena as the Carlsson-Stein tion with respect to SB if the temperature drops below 5300 K. shocks (endnote 50) but also by not assuming SE in spectrum synthesis. Anything you see may not be so much how it appears 37 2D MHD Stagger simulation. The next step in non-E simu- now but rather what it was before. lation was Leenaarts et al. (2007), going from 1D to 2D and The key publications I summarize here are the 1D(t) HD RA- from HD to MHD. [SSX (pdf 167)] is a summary. The three pan- DYN analysis of acoustic shocks by Carlsson and Stein (2002) els at its top are cutouts showing three of the nine quantities in and the 2D(t) MHD Stagger simulation of network and internet- Fig. 1 (pdf 4) at a sample time step. They demonstrate that the work shocks by Leenaarts et al. (2007). The first showcases the simulation (evolved after starting from an LTE initializing sim- physics, the second demonstrates the effects. Neither includes ulation) portrays a quiet solar scene with two opposite-polarity spectrum synthesis but the hydrogen rate analyses of the first MCs 8 Mm apart that resemble network fluxtubes, with inter- and the hydrogen population results of the second make these network in between showing a canopy-like dome with cool to classics of solar spectrum formation. very cool gas and less steep upward density decay underneath. For me the 9-panel thumbnail in [SSX (pdf 167)] is a button that 1D HD RADYN simulation. Carlsson and Stein (2002) studied opens a movie running the whole simulation. Please download hydrogen partitioning in acoustic shocks as those of Carlsson it here and play it. It vividly shows how shocks (thin blue fila- and Stein (1997) in an exemplary analysis. Key results are sum- ments in temperature) travel up along the MCs and more erratic marized in [SSX (pdf 165)]. The main agent is Lyα. They found between them, the latter kicking up the apparent canopy. The that detailed radiative balance is generally a close approxima- first represent field-guided dynamic fibrils[ 57 dynamic fibrils] tion for Lyα at chromospheric heights (the discussion on VALIII as evident from their characteristic time-slice parabolas in the p. 662 (pdf 29) concerns larger heights): about as many photons lefthand column of Fig. 3 (pdf 6). The internetwork shocks are go up as down, many more than collisions up or down. Most pho- acoustic Carlsson-Stein shocks (endnote 49). tons are created per pair d but then resonance scatter many times The last panel of the three in [SSX (pdf 167)] and of the nine without coming far because their mean free path is small (1 cm in Fig. 1 (pdf 4) and its movie version show the NLTE departure in the center panel of [SSX (pdf 32)]) with scattering thermal- coefficient b2 of H I level n = 2, upper level of Lyα and lower ization length about 10 km due to the immense n ≈ N density level of Hα. Its variations are enormous, reaching values up to 1 H 12 39 while the probabilities of destruction a and Balmer detour h are b2 = 10 in cooling post-shock blobs . This is not NLTE at small [SSX (pdf 90)]. This approximation was imposed as input the 0.1-1 dex level traditionally regarded as significant in solar by Leenaarts et al. (2007) to avoid quantifying large numbers and stellar abundance studies – this is 12 dex of retarded over- that effectively cancel. population! The next to last panel of Fig. 1 (pdf 4) and its movie Fig. 3 (pdf 4) of Carlsson and Stein (2002), reproduced in the version shows the departure for the H I continuum level (n=6 in α [SSX (pdf 165)] summary, illustrates this Ly photon balanc- 38 In the VALIIIC“chromosphere” Balmer continuum pumping and Hα ing as negligible net 1-2 rate. The fat transitions there (large net photon losses represent large but balancing terms in the radiative en- photon rate) are photo-ionization in the Balmer continuum and ergy budget [SSX (pdf 69)], so that losses in Ca II H & K, the Ca II in- recombination with ∆n = 1 steps[ 56 Rydberg H I lines] end- frared lines and Mg II h & k dominate [SSX (pdf 70)]. The total curve in ing with fat Hα. For given n2 population (highly non-E in post- [SSX (pdf 70)] defines the unspecified local heating needed to maintain shock cooling) this loop obeys SE but it is far out of LTE because the VALIIIC atmosphere (endnote 26), much more for its chromosphere than for its 10× thinner transition region. 37 Outdoors example: airliner contrails on our sky signify that aircraft 39 This panel also shows high-up green arches. I think that these are ar- engines passed before. If you don’t know that aircraft exist you may tifacts stemming from the input tractability assumption that Lyα obeys have a problem to explain white clouds appearing in long isolated criss- detailed radiative balance everywhere (it does lower down). crossing lines. Compendium solar spectrum formation 14 8. NLS – NLW – NLT chromosphere spectrum the atom model) and illustrates that Balmer continuum pumping (endnote 39). Various Bifrost runs included NLT for hydrogen adds up to another 3 dex for the coolest clouds (instantaneously). (as in the classics above), but not in subsequent spectral synthe- Fig. 2 (pdf 5) of Leenaarts et al. (2007) is also best studied in sis except the NLT analysis of EUV helium lines by Golding its movie version. Its panels show behavior of the y-axis quan- et al. (2017). tity along the dashed vertical lines in the fourth panel of Fig. 1 (pdf 4), at left sampling the lefthand MC and at center sampling internetwork. In the lower three rows of the Fig. 2 movie the 8. NLS – NLW – NLT chromosphere spectrum thick curves are the simulation results, the thin curves the val- SE is usually a sound assumption for thick spectral modeling of ues obtained by assuming LTE instead so that the curve sepa- photospheric structures (SE NLTE) and for thin spectral mod- rations correspond to the departure coefficients in Fig. 1 and its eling of coronal structures (SE CE). Chromospheric structures movie version. Playing this movie shows shocks running up with in between are the hardest to model because NLT reigns next to time in the temperature curves in the lefthand top panel. The ac- NLS and NLW in making them cycle through thick and thin[ 60 tual shocks in the Fig. 1 movie tend to go slanted in the inter- transition region versus heated chromosphere]. For a review network, so that the righthand Fig. 2 movie column samplings see Carlsson et al. (2019). are not along-the-shock path histories[ 58 internetwork shocks in Ca II H and Hα] but sample different shocks at different times while the lefthand movie column samples dynamic fibril Quiet chromosphere. Away from active regions the chromo- shocks better along their near-vertical paths. The movie strik- sphere consists of clusters of MCs making up plage and network, ingly demonstrates that the ionization fraction, proton density the latter partially outlining clapotispheric internetwork “cell in- and n2 population density in the lower rows reach LTE values in teriors” with occasional MCs on their way to cell boundaries[ 61 shocks but thereafter do not drop steeply along with the tempera- internetwork fields and basal flux]. ture. Instead they hang roughly at the value they got in the shock. In the optical the quiet chromosphere is seen only in Balmer There should be eventual settling to the LTE value, but long be- lines, Ca II lines and He I lines. The Lockyer definition implies fore that the next shock arrives: shocks repeat faster than the it is made up by the fibrils seen around network in the Hα core settling time scale. The persistently-high n2 population results and with lesser extent in Ca II 8542 Å. The scene in Ca II H & K in the dark appearance of dynamic fibrils throughout the down- is more complex due to larger Doppler sensitivity enhanced by stroke phase of their timeslice parabolas40 whereas they would PRD and not skipping the clapotisphere in inner wings. Hα fib- fade fast with the lowering temperature if SE LTE were valid. rils around network are largely made by spicules-II[ 62 spicules- The second row shows that hydrogen generally does not ion- II]. ize fully in the shocks but only a few percent. However, this al- In AIA UV images the quiet chromosphere is transparent. ready means two orders of magnitude increase of the electron They sample clapotisphere underneath and yet deeper in MC density which in neutral-hydrogen gas sits at the 10−4 abundance holes in the photospheric surface. These are cooler than their sur- fraction of the electron donor metals (VALIII Fig.47˙ (pdf 67) roundings but appear brighter (Fig. 75 (pdf 91) of Rutten 2020). [SSX (pdf 72)]). This increase in collision frequency appears In AIA 304 Å images quiet areas show two major chromo- sufficient to obtain LTE Boltzmann equilibrium in Lyα within sphere constituents: ubiquitous extended heating patches and the shocks, giving LTE Balmer extinction (bottom panels of the sparser brighter feet of coronal bright points (CBP). For a quick Fig. 2 movie). Another condition for this equilibriating is that a overview inspect the 10 “triple” image sets in Figs. 17 (pdf 35)– feature should be thick enough to contain Lyα, optically thicker 46 (pdf 64) of Rutten (2020)[ 63 SDO quiet-chromosphere im- than the thermalization length from its surface. At partial ioniza- ages]. The first of each triple is a large-field 193 Å image show- tion hydrogen retains sufficient n1 opacity to keep these shocks ing coronal activity, coronal holes, and CBPs. The second is effectively thick in Lyα. a clipped “fire detector” multiplication of precisely co-aligned 304 Å and 131 Å images in which grey patches show heated 3D MHD Bifrost simulations. Since these classic simulations quiet chromosphere and cyan-colored pixels mark CBP feet. The the Oslo simulators went to 3D(t) MHD Bifrost[ 59 Bifrost pro- third is the corresponding HMI magnetogram showing what un- gram] of Gudiksen et al. (2011) [SSX (pdf 159)]. The Bifrost derlies and probably causes the CBPs. The 10 triples sample low stars are far more solar-like than the Avrett stars by being 3D(t) and high activity and low and high latitude and permit inspecting with dynamism and magnetism. Unlike the Avrett stars they con- multiple coronal holes. tain realistic granulation, overshooting granulation and gravity The heated chromosphere patches in the detector images ap- waves, acoustic waves (with box modes as p-mode surrogate) pear everywhere around network but fewer in coronal holes. and shocks, MCs arranged in network patterns, dynamic fib- Their surface patterns are cospatial with the fibrilar Hα canopies rils, bipolar-network-connecting chromospheric fibrils, and re- around network (Appendix B (pdf 11) of Rutten 2020) and serve connective Ellerman bombs (EB). However, they still lack in in my SDO alignment pipeline[ 64 SDO alignment programs] spicules-II and attendant fibrils around network as major quiet- to cross-align AIA 304 Å images directly to HMI magnetograms. chromosphere ingredients, and their cube sizes are yet too small I attribute this component to the heating phase of spicules-II for supergranulation, active regions, and larger-scale coronal fine ejected from network. These also occur in coronal holes; I sus- structure. pect that the lack of corresponding 304 Å heating patches there Bifrost spectral synthesis has mostly concentrated on the is from lack of EUV irradiation from above that elsewhere en- IRIS diagnostics. Obviously it beats Avrett-star spectral synthe- hances He I ionization. AIA 193 Å images do not show these sis in NLS, although most analyses rely for tractability on colum- quiet-chromosphere heating patches but outside coronal holes nar modeling with the RH 1.5D program of Pereira and Uiten- they do show wide-spread diffuse brightness providing such ir- broek (2015). The Avrett stars handle NLW better with respect radiation as a diffuse blanket. to ultraviolet ionization and the problematic ultraviolet line haze CBPs are sparser but not rare in quiet areas. They are best seen in 193 Å images as small ensembles of EUV-bright loops 40 Seen well in the SST/SOUP Hα movie of endnote 57. Compendium solar spectrum formation 15 8. NLS – NLW – NLT chromosphere spectrum located at and connecting bipolar network (unarXived LRSP section to two exemplary active region scenes in different diag- by Madjarska 2019.) They appear as mini active regions, are nostics. likely due to ongoing flux emergence and flux assembly into net- The first scene is the DST/IBIS – SDO multi-wavelength work, and are probably heated by reconnection near their loop image mosaic of Cauzzi and Reardon (2012) and Reardon & tops with subsequent heating of their chromospheric feet[ 65 Cauzzi (2012a, 2012b) showing AR 1092 near disk center on CBP foot visibility](Madjarska et al. 2021), likely by particle 2010-08-03 in the early days of SDO41. I received these concur- beams (Frogner et al. 2020, cf. Rutten 2020). These features can rent co-aligned images from K.P.Reardon and combined them live fairly long, the smallest[ 66 SolO campfires] the shortest, into this pdf blinker.42 The IBIS selection samples Ca II 8542 Å, but AIA 193 Å movies show frequent change suggesting con- Hα (also Hα core width and core-minimum Dopplershift as de- tinuous renewal. The representative AIA 193 Å images in the fined in Cauzzi et al. 2009), and the (small) equivalent width of triple blinkers suggest that CBPs are the main agent in quiet- He I D3. Sun coronal heating and cause the diffuse AIA 193 Å brightness The SDO selection adds comparable HMI and AIA diagnos- spread widely around them, together and with contributions from tics. This IBIS–SDO blinker represents an active-region coun- prior CBPs. Such diffuse brightness is absent in coronal holes terpart to the final showex inspection of the quiet-Sun SST and where heating does not spread to surrounding surface areas along SDO sequences in my alignment practical. These two solar dis- closed loops but instead follows open-field up and out, likely as plays are utterly different. The SST set is a time sequence of coronal plumes (Wang et al. 1997). a small very quiet coronal-hole field with only some network. To reproduce such dynamic nature numerically, even for just The IBIS set is only a single snapshot, not permitting time-delay “quiet” chromosphere, poses enormous challenges because the cause-and-effect sequence inspection, but shows a larger active- MHD simulations must produce spicules-II as copiously as ac- region field containing dense unipolar plage, mossy plage[ 70 tual solar network does and must then be combined with 3D(t) moss], a large regular spot, and even a small flare originating NLS–NLW–NLT (“non-E”) spectral synthesis including PRD at in the spot which behaves as a CBP with loop brightening and least for Lyα to recover the observed Hα chromosphere, and in feet heating. Blinking these images represents a vivid course addition must produce small-scale bipolar activity causing and in active-region phenomenology and also exhibits aspects dis- governing CBPs as observed. cussed above[ 71 IBIS–SDO active-region blinker]. ALMA has the potentiality to sample and measure the chro- One striking difference between quiet and active scenes is mosphere sharper and faster than DKIST. Plane-parallel mod- that quiet closed-field areas permit Hα and He II 304 Å co- eling and various simulations have suggested that ALMA sam- alignment by reversing Hα intensity (endnote 60, endnote 68, ples the clapotisphere dominated by acoustic waves on their way Appendix B (pdf 11) of Rutten 2020) whereas active areas per- to become shocks but these efforts erred in assuming SE. Be- mit co-alignment using normal Hα intensity (blinker p6 – p10). cause Lyα-controlled NLT enhancement of the H I n2 population I attribute this difference to different ratios between roundabout also gives enhanced hydrogen ionization I have instead predicted (“around”) chromosphere heating that dominates in quiet closed- that NLT-opaque Hα canopies are similarly NLT-opaque in the field areas and upward chromosphere heating that dominates in H I free-free continua sampled by ALMA (Rutten 2017b with the hearts of active network and plage from more numerous and a tutorial in Rutten 2017a)[ 67 RBE–RRE–fibril ionization- denser MC clustering[ 72 around and upward heating]. I also recombination]. This means that, while ALMA is a thermome- wonder whether or how upward chromosphere heating in the ter thanks to free-free source function LTE, due to NLT extinc- hearts of active network and plage contributes to the slow-wind tion NLTE it more likely samples Hα canopies than internet- FIP composition bias[ 73 FIP effect]. work shocks underneath – a blessing because we understand in- The enigmatic appearance of the chromosphere in optical ternetwork shocks for decades but do not understand the chro- He I lines may also have mostly to do with upward heating[ 74 mosphere. The good correlation between ALMA intensity and optical He I lines]. Hα width (both gauges of temperature, see Cauzzi et al. 2009) The second display I link to is an appetizer for things hope- in Fig. 4 (pdf 6) of Molnar et al. (2019) confirms this view[ 68 fully to come: comparable active regions in Hα and Lyα in Fig. 2 ALMA–GONG–SDO alignment]. (pdf 6) of Rutten (2017b)[ 75 Lyα dream]. These scenes both Recently Eklund et al. (2021) modeled ALMA response us- sample hydrogen atoms but differ strikingly in appearance. The ing the public Bifrost simulation of Carlsson et al. (2016) which lines have very discordant formation in the ALC7 atmosphere included non-E hydrogen balancing. It emulates active network, [SSX (pdf 28)] but nevertheless should sample common ioniza- but unlike what the actual Sun would show for such area it does tion features, I think with bright Lyα grains showing dynamic not contain spicules-II emanating outward from its network nor heating and dark Hα fibrils showing NLT cooling[ 76 Lyα fea- Hα return fibrils constituting retardedly opaque canopies around tures]. ALMA samples hydrogen ions and may be able to mea- its network. This suggests that some process lacks in the Bifrost sure both types of feature at fast enough cadence to establish and MHD; the 2D simulation of ALMA diagnostics by Martínez- correlate cause–effect delays. It may also sample hydrogen lines

Sykora et al. (2020) points to small-scale charge separation pro- 41 cesses. Hence the modeling of Eklund et al. (2021) applies only The deterioration of AIA sensitivities, worst for chromosphere and to very quiet non-shielded areas as the ones offered in endnote 48 prominence telltale 304 Å, during the immensely productive first SDO decade makes me long for AIA-v2 – with yet more pixels and yet faster? and described in endnote 58. AIA 304 Å images can assist in lo- 42 In a sequence suited to pairwise page blinking. Most pdf viewers fur- cating such canopy-free locations. They are rare. nish blinking by flipping pages in-place with the cursor arrow or page up-down keys. Full-page, fit-to-page, presentation-mode viewing may help. Multiple-page jumping may be possible with a back key. Myself Active chromosphere. For chromospheric activity phenomena I use my showex blinker (endnote 64, from the command line as in I refer to Jongchul Chae’s contribution in this school because I my SDO manual) because it auto-blinks between selectable pairs at se- have not worked on these except EBs and FAFs[ 69 Ellerman lectable speed and offers zoom-in magnification with scatter-diagram bombs and flaring active-region fibrils]. Hence no discussion correlation and more. of spot chromospheres, flares, filaments, etc. here. I restrict this Compendium solar spectrum formation 16 9. Conclusion permitting super-sensitive chromospheric magnetometry[ 77 Ry- with the proviso that “simulations are great tools but still toys” dberg H I candidate for ALMA]. (Martínez-Sykora). Full 3D(t) synthesis of all spectral features Finally, there are dark filaments in these active-region scenes for every voxel at all times of a RADYN-style 3D(t) MHD-etc suggesting “chromospheric” neutral hydrogen gas in the corona. simulation is a bridge too far, so that the simulation synthesis Filaments generally appear long-lived but may nevertheless be frontier lies in clever shortcuts (reviews by Pereira 2019 and made dynamically[ 78 filament blobs]. At high resolution and Leenaarts 2020; non-E simplification on p. 20 (pdf 20) of Rut- fast cadence ALMA may also diagnose their formation includ- ten 2019). ing frequent refurbishment and accompanying retarded-opacity Shortcuts as adopting 1.5D columnar synthesis ignoring non- visibility. radial matter and radiation gradients, Holweger-style “inversion” with NLTE masking, neglect of NLTE ultraviolet scattering and neglect of NLT memory by assuming SE may be as misleading Chromospheric and coronal heating. A text as this should ad- as believing plane-parallel models – but just as those they may dress “one of the longest standing unsolved mysteries in all of as- 43 be educational steps along the way. trophysics” , certainly a text on the chromosphere since gener- The tantalizing trophy is spectral understanding of the chro- ally the blame for coronal heating is sought in the chromosphere mosphere which I now define as the hectic domain separating – while heating that to observed temperatures is less glamorous dense near-LTE near-neutral photospheric gas and tenuous near- but requires more energy for its larger mass. CE ionized coronal gas, where hydrogen intermittently ionizes The literature on outer-atmosphere heating is bewilderingly and recombines cycling between neutral and plasma state with varied. The traditional emphasis for the coronal part is on loop full NLS – NLW – NLT and HD ↔ MHD ↔ plasma complexity: modeling following Rosner et al. (1978) debating AC (Alfvénic the chromosphere ain’t stacked in layers but is dynamically wave) versus DC (reconnection) mechanisms, both dictated by structured and unstuck in time (endnote 60). dynamical fluxtube topography in the photosphere. Chromo- Key property: frequent dynamic refurbishment. spheric heating is attributed to a yet wider variety of agents. Prime diagnostics: H I Balmer lines, Lyman lines, mm con- However, regarding the chromosphere as a layer (not dynam- tinua, helium lines – all enriched by retarded opacity enhance- ically structured and pervaded) and coronal loops as homoge- ment in cooling gas. neous structures (not dynamically multistrand and multithermal) Lowest-hanging fruit: spicule-II heating of quiet chromo- both seem as misleading as regarding the photosphere a plane- sphere and CBP heating of quiet corona. In activity: moss, the parallel layer without convection, waves, magnetism. FIP effect a ripening plum, the Wilson-Bappu relation overdue. The theme here is spectral evidence. My impression from High-hanging fruit: 150+ years since Lockyer (1868) [SSX blinking the SDO triples (endnote 63) and the IBIS – SDO im- (pdf 12)] coined “chromosphere” in the dawn of astrophysics ages (endnote 71) is that a more important split than “corona ([SSF (pdf 2)] ff) we Khayyam-wise have no explanation for his versus chromosphere” is “quiet versus active”. bright yellow “helios” line. Cracking this nut likely helps to un- The SDO triples suggest that quiet chromosphere heating, derstand prominence visibilities. mapped best in the AIA fire detector product, is AC-type with Every aspect treated here must be considered: the Sun is a RBE and RRE spicules-II fired from network and producing most beautiful complex being. We are fortunate that our kind star fibrilar Hα and He II 304 Å canopies. This is “around” heating offers such great inspiring complexity to keep us on our toes. (endnote 72) and does not affect the corona. Quiet corona heat- ing arises chiefly from CBPs, especially the larger ones, with DC loop-top reconnection resulting in diffuser coronal connectivity Acknowledgements. I thank Ram Ajor Maurya and M.K. Ravi Varma for inviting me to teach this subject at their Calicut school. I enjoyed my teaching to nearby CBP areas. Their heated feet are also well-mapped sessions but suspect that the (unseen) students came out by the same door as in the AIA fire detector product. These two quiet-Sun heating in they went [SSF (pdf 120)]. Earlier invitations to teach, in Bandung, Mitaka, agents appear remarkably unrelated, even though both stem from Palo Alto, Seoul, La Laguna, Freiburg, Newcastle and Weihai (usually a full MC dynamics in network and plage outside active regions. week for this material) led to displays SSF and SSX; I hope more will come in a post-covid future. I also thank the organizers for requesting this writeup and The IBIS – SDO images suggest that active regions add endorsing hyperlinks. F. Paletou, J. Leenaarts, T.M.D. Pereira and P.R.Young unidentified “upward” heating (endnote 72) feeding long-loop contributed improvements. As always I relied much on ADS. ADS also taught connectivity to distant other active regions, likely with FIP abun- me pdf page opening and keeps these classic linkers active since its 2019 dance offsets (endnote 73). revamp. Stackexchangers solved latex problems. M.J. Rutten started hosting robrutten.nl for link speed and persistence.

9. Conclusion Version history March 3, 2021 – initial pre-school arXiv post. Understanding the solar spectrum has come a long way but isn’t April 29, 2021 – second arXiv post honoring C. de Jager’s 100th birthday: more content, more links, faster serving. there yet. The basic physics was well understood by the 1970s – May 7, 2021 – notes → endnotes. no dark matter nor dark energy in the solar atmosphere.44 Nu- June 1, 2021 – endnotes on yellow pages, third arXiv post. merical modeling started about then with the advent of elec- July 20, 2021 – more informative reference format. July 29, 2021 – reached deep-thought page limit 42. tronic computing which was welcome because analytic theory August 3, 2021 – fourth arXiv post. had reached its limit in Kourganoff (1952). Numerical study August 26, 2021 – improvements. provided great insights even though initially mostly within the misleading straitjacket of plane-parallel atmospheres. Nowadays the major challenge is to obtain realistic spectrum synthesis from realistic solar-atmosphere simulations to provide guidance to interpret high-quality observations conveying solar truth – 43 Joan Schmelz, ISSI talk July 2021 for a covid-zoom audience. 44 For a conflicting speculation search ADS with abs:"axion quark" and "solar orbiter" (not to be cited). Compendium solar spectrum formation 17 Notes

Notes C. (Kees) Zwaan50 took over the Utrecht radiative transfer courses in the 1970s, modernized RTSA and set up IART. I Welcome to the yellow pages! These endnotes are take-or-leave copied much of his (neatly hand-written Dutch) material when additions adding further explanation including linked literature, I took over IART in 1985, RTSA in 1994. comments, opinions, and even history (also personal). I then also started teaching abroad, appreciating getting to Topic selection: personal, favoring topics that I have been know an appreciable fraction of later colleagues worldwide.51 involved in and references and figures I knew best and came first I quit teaching at Utrecht at my mandatory retirement in to mind in this quick writeup for the Calicut school. With prefer- 45 2007. In 2011 Utrecht University abruptly quit its astronomy ence for my own displays when these open faster than ADS. (even its A&A and ApJ subscriptions) and effectively killed Navigation: each note ends with a link [Main call] returning Dutch solar physics. The DOT is mothballed since. I still teach to its calling location in the main text and a link [Back] that goes abroad when invited but without Utrecht affiliation. there when you came from there or back to another link that you During my time the advent of electronic computing, elec- used – working in Firefox, Ubuntu acroread and xpdf, macOS tronic writing, the internet and ADS drastically changed the way Preview but not in Ubuntu evince, chrome, macOS Safari (try we work. Because these advances may appear granted to you I the back button instead). Footnote numbers jump back to their sprinkle some reminiscences through these notes. [Back][ Main call. call] 1. Personal background 2. Textbooks I developed this material in decades of teaching at Utrecht Uni- versity and abroad. In the 1960s I learned solar spectrum for- An inventory of relevant textbooks is in the [RTSA bibliography mation from C. (Kees) de Jager46. Our bible was Unsöld’s book (pdf 17)]. I found Mihalas (1970) easier to read then Mihalas (endnote 42) but NLTE theory had started and was loudly advo- (1978). Newer books are Oslo-triggered Castor (2007) and com- cated to us by A.B. Underhill, then also professor at Utrecht. The prehensive Hubený and Mihalas (2014), both more advanced Back Main call guru was R.N. Thomas but I found his writings incomprehen- than this course. [ ][ ] sible and Kees much clearer47 (see also footnote 93). Jefferies’ 3. ADS page opening, silly publishers (1968) book48 was a welcome improvement on Thomas-speak49. I had started on solar radio diagnostics (Fokker and Rutten I describe linking citations to the ADS abstract page and open- 1967) but switched to optical spectrum formation thanks to De ing particular pages in ADS-accessible pdf files on this webpage. Jager’s lecturing and then concentrated on the chromosphere The latex macros and usage examples are given in my report thanks to J. Houtgast who took me as junior help to both 1966 recipe for astronomy students. ADS pdf opening does not work eclipses to record the flash spectrum. Half a century and three for publications at what I call “silly publishers” who do not fur- Saros 142 periods later the chromospheric spectrum still fasci- nish direct pdf access on the ADS abstract page; I avoid citing nates me, as evident in Sect.8 which is more a progress report these altogether. More detail in my recipes for publications. than a course. ApJ and A&A furnish ADS pdf openers already after the first year; during the first year ADS links to their publisher html 45 Example: [SSX (pdf 60)] versus Fig. 1 (pdf 3) for the most-shown access. For these Lean Library will silently use your remote- figure of solar physics. Both open relatively slow (but consider foot- institute license to pass directly to the pdf but then opens the note 52). 46 first page, making proper-page arXiv pdf opening a better alter- Cornelis (Kees) de Jager (1921 – 2021) achieved large and fast ex- native. Thus, the immense top-quality astrophysics literature in pansion of Utrecht astronomy after M.G.J. Minnaert (including my job), ApJ and A&A from their start to the present is directly accessi- started Flemish astrophysics at Brussels, started Dutch space research, was instrumental in starting ESRO and its transition to ESA, was the ble: any particular page with a figure or equation or paragraph worth showing is only a single click/tap away in the very sen- founding editor of Space Science Reviews and co-founding editor of So- 52 lar Physics, was IAU General Secretary, twice president of COSPAR, tence where it is cited. This text has many such openers. president of ICSU, and much more – including running marathons and 50 Cornelis (Kees) Zwaan (1928–1999) was my thesis adviser and a giving public lectures up to record age. I honored his 100th birthday close friend. I not only took over his teaching but also his DOT stew- with my second arXiv post and an updated Kees-day SDO album. He ardship. The revolutionary “open principle” of the DOT (no vacuum passed away four weeks later. I put his thesis on solar hydrogen lines but rely on wind flushing to avoid convective turbulence from heat- (de Jager 1952) on ADS with help from Oslo and am co-authoring an ing in the convergent beam to focus) was his idea and proposal to “In Memoriam” for Solar Physics. R.H. Hammerschlag. DKIST is sort of open. 47 Just as R.M. Bonnet who wrote “As seen through the writings of 51 I much missed face-to-face contact with the students in this zoom Thomas, radiative transfer for me would have remained an unintelli- school. Just talking to your own screen is no good. I also missed who- gible theoretical exercise forever” and found De Jager much clearer, what-where student presentations and joint dinner and pub sessions. acknowledging “the clarity of Kees’ views” (page 2 (pdf 2) of Bonnet 52 1996). Unfortunately browsers reload the full pdf for every page, taking a while. But during these long seconds you might comfortingly contem- 48 Jefferies (1968) was scanned by A.V.Sukhorukov and put on ADS plate the much slower pre-internet/ADS chore: go to your institute li- by me with Jefferies’ consent, enabling page openers here. brary or some dark storage for older publications, locate the journal 49 Telltale specimen: NASA “orange book” Stellar Atmospheric Struc- shelf, find a ladder, get the volume, find the article, find the page, make tural Patterns (Thomas 1983). Jordan (1996) called it “hard-to-read notes or get to the nearest Xerox, write the full reference on an index vintage Thomas”; L.E. Cram suggested to read it ten times for proper card (for alphabetic reshuffle by the manuscript typist), return the vol- (Thomas speak) appreciation. ADS lists 26 citations versus 250 for ume to its place, return the ladder, get back to your office. Jefferies (1968), 408 for Mihalas (1970), 2097 for Mihalas (1978). Current on-line procedure: click/tap the citation in the text to jump to the However, Thomas (1921 – 1996) was truly a game-changer in stellar- reference (or go there), open the ADS abstract page or DOI-linked pub- atmosphere RT. JILA (which he co-founded) celebrated his centennial lisher page (or find it), locate a pdf button and load the pdf, turn pages while I wrote this which made me reread my 1997 Thomas memorial to the desired one. Still too many actions and prone to silly-publisher and realize that in publishing this link-strewn compendium I finally ful- mazes. You also may have trouble returning where you were reading. fill my “On publishing” (pdf 9) promise there. Single click/tap page opening in parallel to your reading is much better. Compendium solar spectrum formation 18 Notes

In contrast, Springer promises “documents53 at your finger- tips” but requires many clicks/taps that may impose hefty pay- 6. Equivalent temperatures ment and do not bring you to the right page. Solar Physics has These serve for deWiening Planck function temperature sensi- no ADS pdf openers since Springer bought it in 1996 whereas tivity in intensity-related and population-related quantities. They ADS still page-opens scans from the preceding Reidel-Kluwer are formal by being defined as the temperature value that entered era. Living Reviews in Solar Physics is open access but has no into the pertinent TE expression equals the quantity: brightness 54 pdf openers on ADS since Springer bought it. [Back][ Main temperature Tb for Iν into the Planck function, radiation tem- call] perature Trad for Jν into the Planck function, excitation temper- ature Texc for actual level population ratio into the Boltzmann l 4. Spectral bandwidth and vision ratio making Sν = Bν(Texc), ionization temperature Tion for actual Only our eyesight is adapted to solar radiation temperature but level/continuum population ratio into the SB expression making c not simply to the Planck function for the solar effective tempera- Sν = Bν(Tion) [SSF (pdf 72)] [RTSA Eqs.2.128 – 2/132 (pdf 57)] ture. Stating that needs specification whether our retinas measure after Wijbenga and Zwaan (1972). Thus, the difference Texc − Te l per wavelength, frequency, wavenumber or count photons – see is a deWiened alternative to the Sν − Bν measure of NLTE affect- [SSF (pdf 50)] and remember dν=−(c/λ2) dλ. ing a source function. Scattering lines with small ε likely have The third panel shows why J.W. Brault said in 1978 that the Texc ≈Trad. − intensity at Ca II H & K was too low for his Kitt Peak Fourier The required [B] 1 operator (Eq. 1 (pdf 1) of Samain 1979) Transform Spectrometer: a photon-counting wavenumber de- is straightforward but needs attention to units. My IDL vice (but we got the only center-limb H & K scans in the FTS brighttemp.pro handles it for Bν and Bλ and different band- archive). The last panel shows why high-energy astrophysicists width units. plot λBλ or νBν. These conversions should always be applied when compar- History; our eyesight was optimized in our aquatic past for ing intensity-related quantities at different wavelengths. Exam- the penetration of sunlight in sea water, as for giant squids ples: Fig. 4 (pdf 5) of Rutten et al. (2011) comparing brightness whose eyes have similar but better design (Feynman I Fig. 36- scenes in Na I D1, Mg I b2 and Ca II 8542 Å, their atlas profile 10) and might suit nighttime astronomers better. As early mam- comparisons in Figs. 10 – 12 (pdf 10) there, and the intercompa- mals ashore we became dichromats being nocturnal to hide rable S, B, J panels for many lines from the ALC7 star in the from predating dinosaurs. After these were blasted into obliv- displays starting at [SSX (pdf 82)].[Back][ Main call] ion early hominids restored red-detecting cones (perhaps to spot low-hanging fruit) giving us color-triangle vision (Feynman I 7. Intensity calibration Sect. 35) to appreciate Kodachrome Basin. [Back][ Main call] For the optical spectrum Neckel (2003, 2005) provides limb- darkening expressions based on Brault’s FTS observations used 5. Solar telescopes by Neckel and Labs (1984) that Neckel (1999) converted into a These are intensity telescopes and therefore photon-starved, calibrated atlas (endnote8 ). worse when larger.55 Reason: for pixels that resolve the diffrac- In the ultraviolet the SUMER atlas of Curdt et al. (2001) tion limit the flux of photons per pixel does not vary with aper- provides calibrated disk-center spectra for the 668 – 1611 Å ture size while smaller pixels sample smaller solar surface area range, with center-to-limb observations available in the SUMER (that’s the game) for which the solar change time (say crossing archive. Earlier 1400 – 2100 Å absolute-intensity and center-to- at solar sound speed) is faster. limb estimates are in Samain (1979). Calibrated 1100 – 1900 Å DKIST therefore needs four times faster data-taking ca- flare spectra are shown and provided by Simões et al. (2019) in dences than the SST unless it is used as (expensive) light bucket. their inventory of spectral contributions to flare intensities in the DKIST thus has four times less time (only 1 s at Hα) for multi- AIA UV passbands. [Back][ Main call] frame collection for the required speckle or MOMFBD restora- tion. The latter technique (van Noort et al. 2005) needs fewer 8. Optical spectrum atlases and line lists frames which is an important boon. At the DOT we used two- Solar spectrum atlases started with the beautiful engravings of channel speckle reconstruction (Keller and von der Lühe 1992) Fraunhofer and Kirchhoff, Rowland’s photographs, and then the to reduce speckle burst sizes in scanning through Hα by adding graphical Utrecht Atlas of Minnaert et al. (1940)56, all four sam- wide-band imaging as is also done for SST MOMFBD. pled in [RTSA Fig. 3.1 (pdf 26)]. All lines in the Utrecht Atlas It is easier to forego diffraction-limit resolution so that larger were laboriously measured57 to produce the solar line table of aperture light-buckets more photons, the ELT approach. On Moore et al. (1966). I have it as ascii file courtesy J.W. Harvey the other hand, fortunately our resolvable star kindly offers (in the magnetic-tape era) and consult this regularly although I wavefront-encoding granulation everywhere on her surface so also have the book. that we do not need a laser star for every isoplanatic patch as ELT should have but won’t. [Back][ Main call] 56 History: made by scanning photographic spectrum plates taken at Mt. Wilson with an ingenious microdensitometer using analog nonlin- 53 Springerese for publications, not passports. “Reviewer” means ref- ear “calibration curve” conversion of emulsion opacities into solar in- eree, “ticket” means into a black hole unless “escalated”. tensities, illustrated in the prefaces (English and Esperanto). 54 Example: opening some of the 62 informative figures in LRSP Nord- 57 History: manually by counting mm2 area for each line dip to get its lund et al. (2009)( endnote 41) would suit well in this display-oriented equivalent width, for which the atlas was printed on mm-grid paper. course. The ADS abstract page shows a pdf-promising green-dotted pdf While I was a student De Jager introduced a planimeter for faster mea- button but it serves the hassle of navigating the Springer maze. Among surement. Later I co-developed a sequence of electronic tracing readers solar physics publications the LRSP reviews are top candidate to page- until replacing the chart recorders of Houtgast’s microdensitometer and link but notwithstanding their open access this is obstructed by Springer. the Utrecht solar spectrograph with a digital paper-tape puncher (Rut- Hence no LRSP page opening in my courses and reviews unless arX- ten and van Amerongen 1975). Both instruments are now part of the ived. This one wasn’t. Other such suited but unarXived LRSP’s there- De Jager-started Museum Sterrenwacht Sonnenborgh but regrettably fore remain uncited here. the impressively noisy tape-spewing (120 bytes/s!) punch got lost. 55 Nighttime colleagues may not believe this. Compendium solar spectrum formation 19 Notes

While writing this note I realized that I might supply this ergy excess of the unstable state above that, but if the fresh ion file and similar leftovers from the past and started a solar file is excited it decays, in CE always per line photon. archive supplying them. They include the legendary Revised The other path out of the unstable level is decay of one (or Multiplet Table of Moore (1959) (RMT) used as input to MMH both) of the excited electrons, in CE also always per line photon, = Moore et al. (1966), the more precise wavelength list of Pierce into a configuration with the second electron in a regular bound and Breckenridge (1974), the yet more precise wavelength list level below the nominal ionization limit (“radiative stabiliza- of Allende Prieto and García López (1998), the measurements tion”). for 750 clean optical lines in the Jungfraujoch disk-center atlas Autoionization often has higher probability than spontaneous (Delbouille et al. 1973) starting at Table IV (pdf 11) of Rutten stabilization. and van der Zalm (1984b), the corresponding measurements for Dielectronic ionization: an energetic collider, usually hitting 602 clean optical lines in the Sacramento Peak irradiance atlas of the lower ion in its ground state, kicks up two electrons with Beckers et al. (1976) starting at Table I (pdf 4) of Rutten and van the first one taking a good part of its kinetic energy. With the der Zalm (1984a), and the flash-spectrum line intensities starting remainder the second may ionize bound-free but is more likely at Table 3 (pdf 51) of Dunn et al. (1968) that I got laboriously put into a high level. If this is still below the ground-state limit punched on Hollerith cards in the 1980s. both electrons decay radiatively, the second likely in a cascade, The next printed atlases were the Jungfraujoch photoelec- together a two-electron multi-level version of pair d converting tric scans of Delbouille et al. (1973) and the Sacramento Peak kinetic energy into line radiation without ionization. However, flux scans of Beckers et al. (1976). The first was unwieldy from if the second electron reaches an unstable level above the nom- too large page size, the second had practical letter size and also inal limit then it may be freed by autoionization before its ra- MMH line identifications in the margin. I toted my dog-eared diative decay, a sequence called “excitation-autoionization” by copy to all my spectroscopic observing ever since, annotating Dere (2007). The other still bound excited electron decays to the every line of special interest. Both atlases also became available ground state of the new ion, also producing one or more line on magnetic tape (storing unprecedented millions of bytes) but photons. In this double-excitation manner the energetic collider the digital market was taken over by J.W. Brault’s Fourier Trans- yields higher ionization probability than for one-electron colli- form Spectrometer at the Kitt Peak McMath telescope delivering sional ionization by losing energy that escapes as line photons. precise atlases still available on the NSO website if you master Dielectronic recombination: in the reverse sequence the ener- the non-trivial art of navigating that. getic colliding free electron kicks a bound electron up to a high The FTS atlas that I use habitually, not available at NSO, level and is itself captured into an unstable state which stabilizes is the “Neckel” disk-center intensity atlas named after Neckel already before re-auto-ionization through spontaneous radiative (1999) announcing its availability. Its advantage is that Neckel decay of one of the two excited electrons, often the lower one, to converted Brault’s FTS intensities into absolute units follow- a stable level of the lower stage. If that is still excited it decays ing Neckel and Labs (1984), enabling conversion into brightness radiatively to the ground state; the other electron also, usually temperature for comparing differing spectral regions (endnote6 ). per multi-level cascade. Depending on the branching ratios some I also supply this atlas in my solar file archive.[Back][ Main decays may temporarily end up in a metastable level comparable call] to Ca II 8542 Å in the Ca II detour example in Fig. 2 (pdf 4) of Rutten (2019) but then not go down from there collisionally as 9. Solar neutrinos there but CE-wise radiatively in a forbidden transition such as Fortunately the fusion-released gamma rays escape by slow out- the optical coronal emission lines [SSX (pdf 8)]. ward diffusion converting them into more agreeable sunlight Radiative result: no continuum but line photons ga- [RTSA Shu quote (pdf 110)] plus enthralling dynamo action lore. [Back][ Main call] causing magnetic flux escape and harmonious humming en- abling internal sounding (Athay et al. 2008). In contrast, the neu- 11. Ca II H & K core reversals trinos that are also released in hydrogen fusion are virtually un- On page 128 (pdf 146) Jefferies cited the suggestion of Goldberg stoppable and escape directly from the solar core58. They posed (1964) that the double emission peaks of the Ca II H & K core a big problem last century when too few were detected. It was reversals may be due to dielectronic recombination and also the not a fault of the measurement nor of stellar structure theory demonstration by Noyes (1965) that this does not work for CRD but the blame was on particle physics: part of the solar neutri- in H & K. Noyes used Avrett’s√ brandnew CRD scattering code nos “oscillate” into another flavor on their way here. Photons, producing the canonical ε graphs of [SSF (pdf 80)] in the be- once they escape, fulfill their task of information carrying bet- lief that CRD is a better approximation than CS which he had ter.59 [Back][ Main call] initially assumed in Goldberg and Noyes (1964), but actually H & K are PRD lines. Goldberg’s mechanism seems to be colli- 10. Dielectronic ionization and recombination sional dielectronic ionization of Ca I leaving Ca II ions excited Unstable levels: both processes excite one electron in a normal in the upper levels of H & K – but I note that Ca I is a minority bound-bound transition and put a second electron into a level species and think the idea was a red herring. above the nominal ionization limit of the lower-stage ground It does illustrate how during decades the Ca II H & K rever- state. This is possible because the orbitals are boosted by the sals were the most enigmatic features in the optical spectrum first excitation; they may also be non-hydrogenic affected. Such by being the only Fraunhofer-line departures from regular bell a state is unstable because autoionization may occur, freeing the shape. This extended literature was outstandingly reviewed by electron from there into escape. This usually leaves the fresh ion Linsky and Avrett (1970) (but not mentioning the Goldberg– in its ground state and then the freed electron carries off the en- Noyes dielectronic angle). Athay and Skumanich (1967) had suggested that the H and K dips map a temperature minimum 58 11 2 2 For them we do not need RT theory. Nor worry about the 10 pass- as in NLTE cartoon [SSF (pdf 37)], but then Shine et al. (1975) ing every second through our thumbnails. 59 showed that PRD upsets this notion as in cartoon [SSF (pdf 89)] Solar neutrinos may convey information on scotogenic dark matter with their Fig. 9 (pdf 7) copied in [SSF (pdf 92)]. beyond my and your solar physics interest (de Boer et al. 2021). Compendium solar spectrum formation 20 Notes

However, such plane-parallel modeling could not explain the CRD one and sits in the dip outside the Mg II peak even in that the violet-side peaks in H & K are higher than the red-side strong plage. The apparent Mn I over Fe I brightening of plage peaks (the red Ca II H2R peak is barely visible in the atlas pro- in Fig. 1 (pdf 2) in scans from Malanushenko et al. (2004) re- file in Fig. 11 (pdf 10) of Rutten et al. 2011). In addition, high- sults from mean-field brightness normalization. The surround- quality narrow-passband imaging as in the “Selected spectrohe- ing granulation is actually darker in the Mn I line than in the liograms” of Title(1966) (which I have on my desk but cannot turbulence-sensitive Fe I line while these lines show the same ac- find in ADS) and the remarkableK 2V spectroheliogram taken by tual plage brightness because the fluxtube “holes” schematized B. Gillespie at Kitt Peak that adorned the cover of Lites (1985) in Fig. 4 (pdf 5) reach as deep because both Fe I and Mn I ionize and was partly shown in Fig. 2 (pdf 4) of Rutten (2007) and away. [Back][ Main call] elsewhere showed that the asymmetries come from small bright grains that occur intermittently in internetwork cell interiors and 13. Detour and cross redistribution are not magnetic (although some authors persistently claimed Ca II H & K are the strongest lines in the visible and of the dom- this). inant ionization stage so that the two-level pairs in Fig.3 suf- This wide confusing literature was reviewed by Rutten and fice. In the ALC7 star they have ε ≈ 10−4 [SSX (pdf 88)], as do Uitenbroek (1991) concluding here (pdf 51) that K2V and H2V Mg II h & k [SSX (pdf 89)] while Lyα reaches ε ≈ 10−6 [SSX “cell grains” are an acoustic interference phenomenon along the (pdf 90)]. The other strong lines in the visible are the Balmer lines of Liu and Skumanich (1974) – then brilliantly proven with lines, the Ca II infrared triplet, the Mg I b triplet, and the Na I D observation-driven RADYN simulations by Carlsson and Stein and K I resonance doublets. These are all scattering lines with (1997)( endnote 49). ε ≤ 10−2 but from minority stages60 with sizable contributions However, the Ca II H & K core reversals also govern the fa- by the detour pairs in Fig.3 , for Hα ionization-recombination mous cool-star relation in Fig. 1 (pdf 9) of Wilson and Vainu loops as at right in Fig. 2 (pdf 4) of Rutten (2019) contributing Bappu (1957) which to my knowledge remains unexplained. up to 10% [SSX (pdf 96)]. These detour contributions have no Doing so will make you famous! Since the main contribution coherency. In addition, the Balmer lines have wide Doppler re- to the peaks in solar irradiance is from activity-measuring cal- distribution from the small hydrogen mass and their outer wings cium network and plage on classic spectroheliograms I sug- are formed deeply were collisional Stark redistribution is large gest an approach along the lines of Scharmer (1976) and Ayres (endnote 48). For Na I D2 Uitenbroek and Bruls (1992) found (1979) but applied to fluxtube atmospheres with attention to the PRD signatures only near the limb, none for K I 7699 Å. The slow network waves in the center panels of Fig. 3 (pdf 4) of much stronger Mg I 2852 Å resonance line has significant PRD Lites et al. (1993) and fluxtube shocks as in Kato et al. (2011) effects (Canfield and Cram 1977). while subtracting or including the basal flux contribution (end- Cross redistribution occurs for the Ca II triplet lines. In the note 61). [Back][ Main call] detour sketch at left in Fig. 2 (pdf 4) of Rutten (2019) this concerns the consecutive loop 3933 up, 8542 down and back 12. Ba II 4554 Å and Mn I lines up, 3933 down. It was analyzed by Uitenbroek (1989) who Ba II 4554 Å showed PRD signature near the limb in my eclipse showed in Fig. 2 (pdf 7) that PRD frequency dependence of the observation in Fig. 9 (pdf 17) of Rutten (1978). It was confirmed Ca II K source function affects the Ca II 8542 Å source func- in Figs. 1 – 3 (pdf 4) of Rutten and Milkey (1979) and Fig. 6 tion strongly through the Wien-reversing sensitivity transcrip- (pdf 5) and Fig. 11 (pdf 7) of Uitenbroek and Bruls (1992). The tion (endnote 25). However, since H & K are much stronger these line is a valuable spectropolarimetry diagnostic (Belluzzi et al. large PRD effects do not affect emergent Ca II 8542 Å profiles 2007) and also a good Doppler diagnostic of the upper photo- from plane-parallel atmospheres as modeled there because in sphere (Shchukina et al. 2009, Kostik et al. 2009) because it these the infrared lines are formed well below the heights of PRD combines large atomic mass, hence small thermal broadening, source function splitting for H & K. This is likely not the case with considerable core broadening by isotope splitting and hy- when both H & K and the triplet lines form together in chromo- perfine structure making it suited for filter instruments and so spheric features without such height separation (e.g., spicules-II yielding the astounding Dopplergram in Fig. 5 (pdf 4) of Sütter- and return fibrils, endnote 62). [Back][ Main call] lin et al. (2001). Other lines with large hyperfine broadening are the Mn I ones 14. Dark EUV features recommended by G. Elste to W.C. Livingston as candidates for long-term spectral irradiance monitoring because they are in- When you see the same black features on bright backgrounds sensitive to infamous microturbulence. Livingston found larger (on the disk but outside no-emissivity coronal holes) in multi- ple EUV wavelengths then out-of-the-passband scattering is the cycle-dependent variation in Mn I 5394.7 Å than for other atomic most probable mechanism. Cartoon in [SSX (pdf 10)] (where lines (Fig. 16 (pdf 9) of Livingston et al. 2007). This enhanced “volume blocking” is a bad misnomer); I drew it as Fig. 10 sensitivity was attributed by Doyle et al. (2001) to pumping of (pdf 12) of Rutten (1999) to answer a TRACE question by Mn I 5394.7 Å by Mg II k but erroneously assuming CRD for C.J. Schrijver. It shows that when black features are also present Mg II k in the spectral synthesis of their demonstration. Vitas in He I 584 Å images (not in AIA alas) the agent is photoioniza- et al. (2009) showed that the actual reason is that hyperfine- tion in the H I Lyman continuum. broadened lines indeed lack the “microturbulent” thermal and This may occur in filaments seen in EUV lines as in the granular Doppler smearing through which all narrower photo- IBIS active-region blinker of endnote 71. Other such features spheric lines lose such sensitivity (MURaM demonstration in are “condensation” (meaning recombination) clouds, e.g., in their Fig. 6 (pdf 8) and Fig. 7 (pdf 9)). I reviewed this issue in Rutten (2011, not on ADS but posted 60 For Balmer lines regard only the hydrogen top with n=2 as “ground here) where Fig. 1 (pdf 2) summarizes Livingston’s results and level” because its population is fully controlled by Lyα doing its own Fig. 2 (pdf 3) shows that in plage the Mn I blend (at -0.7 Å from thing (Sect.7 ), making this top a low-abundance severe-minority stage −7 Mg II k center) is not pumped at all because the blend is at the since n2/nH < 10 [SSX (pdf 91)] while the proton fraction np/nH > − wavelength where the PRD source function departs most from 10 4 where the Hα core forms. In this fixed-bottom top Hα is a reso- nance line comparable to the alkali resonance lines (endnote 26). Compendium solar spectrum formation 21 Notes coronal rain (Antolin and Rouppe van der Voort 2012, An- comparable resulting diagrams for both references in my IDL tolin et al. 2015) visible in Hα via non-E recombination cas- edulib. An overview of total extinction including lines in the cade giving large opacity on the disk and large emissivity off- FALC atmosphere is shown in Fig. 6 (pdf 13) of Rutten (2019). limb as in the temporal cycling of Froment et al. (2020)(Fig. 7 The minima are treated in endnote 19.[Back][ Main call] (pdf 10)). [Back][ Main call] 17. Emission measure, DEM 15. White-light corona Colleagues studying the coronal EUV spectrum usually do not The reason why we must travel far to see the corona61 is evident worry about photons after their creation when not ending in their from the visibility diagram in [SSX (pdf 5)] after Fig. 6 (pdf 51) telescope.64 Unhindered travel to the telescope, no extinction, in the classic review by van de Hulst (1953). only emissivity counts in their production per pair d. This must The pearly white color of the K corona suggests photospheric be evaluated for all lines within the passband65 at every location light that is bent towards the viewer without wavelength depen- along the line of sight and integrated with confusion an issue. Lo- dence (unlike blue-favoring Rayleigh scattering, hence Thomson cally the emissivity scales with the product nl Ne of lower-level scattering, Mie scattering, or reflection without refraction). The population density and electron density, the latter representing absence of photospheric Fraunhofer lines in the coronal spec- the required collider for excitation per d which is most likely trum, unlike their presence in the spectrum scattered by your an electron. Regard a resonance line from an ion ground state. nose, was a large disappointment for the first observers bringing Its population enters as (nl/NE) AE NH with element E abun- spectroscopes to totality, played in a role in the debate whether dance AE = NE/NH, hence local emissivity ∼ g f AE Gi j(T) Ne NH the corona is real or an illusion, and remained enigmatic during where Gi j is a conglomerate function of temperature describ- many decades. Thomson and Mie scattering are coherent so re- ing the element partitioning defining nl/NE for this line i j and tain frequencies and hence should retain Fraunhofer lines. Then needing all sophistication of CHIANTI (endnote 18) for evalua- Grotrian (1931) noticed that Ca II H & K leave a washed-out tion within the CE assumption. For a review see Del Zanna and dip and proposed here (pdf 22) that their 100 Å smearing cor- Mason (2018). If the line is from an excited state function Gi j responds to 7500 km s−1 Doppler redistribution [SSX (pdf 6)]. depends also on electron density. Integration over the volume of an unresolved emitting structure (say a loop) yields emission However, he excluded thermal broadening with a reminder that R the solar-atmosphere temperature is at most 6000 K, the deep- ∼ g f AE Gi j(T) Ne NH dV where NH is often replaced by Ne photosphere value – everybody knew that heat flows from hot to taking NH = Np ≈ Ne (1 − 2 B) with B the helium/hydrogen cool so further out may only be cooler. abundance ratio (Sect. 5 (pdf 20) of Pottasch 1964). The density His second-law-challenging attribution to million-K thermal dependencies Ne NH may be removed using the differential emis- −5 −1 electron motion came after Edlén’s high-ionization line inven- sion measure DEM(T) ≡ NeNH (ds/dT) in cm K , a com- tory enabling identification of the enigmatic bright “coronium” bined local gradient measure that defines the emission measure R lines [SSX (pdf 8)]. Detection of synchrotron radio emission integrated over s along the line of sight as EM ≡ N N ds = R e H was the clincher for the greatest solar physics discovery so far. DEM(T) dT in cm−5 (Sect. 7 (pdf 113) of Del Zanna and Ma- Explaining million-K coronal temperatures beyond incantation son 2018). With this density removal the emission in the line “Parker nanoflare heating” remains a challenge while dielec- R ∼ tronic line emission helps keeping the corona cool (endnote 10). from a feature becomes f AE DEM(T) g Gi j(T) dT where For smooth-disk illumination the scattered brightness is lo- the Gaunt factor g is under the integral because it depends on temperature (Eq. 2 (pdf 3) of Withbroe 1978). [Back][ Main cally proportional to Ne, but as pointed out by Minnaert (1930) limb darkening (and perhaps spots) must be accounted for in the call] local scattering source function S = J , and also the integra- ν ν 18. CHIANTI tion along the entire line of sight with problematic confusion between different structures, as evident in the marvelous Druck- CHIANTI website. (Name suggestion: Coronal Holism In Ap- müller eclipse images. preciating Non-Thick Ionization.) Community effort by these The outer-corona F component in [SSX (pdf 5)] is from colleagues. It handles coronal spectroscopy with a large tran- Mie scattering by dust particles that move slower, maintaining sition inventory for evaluating CE ionization equilibria, radia- Fraunhofer-line signatures. [Back][ Main call] tive losses, DEMs and more. The need to include numerous lines contributed by dielectronic processes (endnote 10) makes it ever- 66 16. Extinction diagrams growing. [Back][ Main call] Diagrams of bound-free and free-free He and He extinction for 64 Solar physicists used to be either thin or thick. Thin ones excelled in different stellar atmospheres assuming SB are shown in the four atomic and later in plasma physics, thick ones in 2-level scattering and Figs. 8.562 of Gray (2005) and in [SSF (pdf 53)]. Three ear- later in MHD. They now grow together (but heed my scorn at the end of lier versions of Gray (1992) are copied into [RTSA Fig. 8.14 endnote 34). The Sun doesn’t care and magnanimously sends photons to (pdf 211)]. The classic “confusograms” in Figs. I–XV (pdf 17) of both types adhering to her kind “Principle of Solar Communicativity” Vitense (1951) were redrawn in English in Novotny (1973), with (pdf 9) formulated at the only solar IAU symposium ever in the USSR the solar one shown in [SSX (pdf 54)] and reprinted in [RTSA (Rutten 1990). Fig. 8.5 (pdf 199)] and four stellar ones in [RTSA Figs. 8.15 – 65 For spectral-passband imagers as AIA usually a complex multi- 8.18 (pdf 212)]. Their not-so-trivial interpretation63 is treated in temperature combination from registering multiple lines. Spectrometers [RTSA Problem 10 (pdf 252)]. I have emulation programs and as Hinode’s EIS register spectrally dispersed bands but do not inform precisely where from – slitjaw imaging as by IRIS is highly desirable 61 Personal history: after starting my career with spectroscopy at three – or miss the more interesting happening just off the slit. It is a pity eclipses my wife and I became eclipse tourists. Our traveling in 2020 to that Fabry-Pérot imaging spectroscopy is not feasible at short wave- my 15th, completing three Saros 142 periods, was covid-canceled. lengths; SST/CHROMIS is the shortest (Scharmer et al. 2019). My par- 62 Erratum: Fig. 8.5a has (in Gray’s words) a small error of 10−26. adise dream is slitless 2D imaging spectrometry resolving wavelength 63 across Lyα (continued in endnote 75). When E. Vitense wrote her textbook Böhm-Vitense (1989) she didn’t 66 copy her thesis graphs but simplified them in Fig. 7.8 = [SSI (pdf 7)]. Example from P.R.Young for He-like O VII: currently CHIANTI specifies 577 levels, 127 below the ionization threshold and autoion- Compendium solar spectrum formation 22 Notes

(relatively increasing) electron density defining the collision 19. Deepest wavelengths frequency handling destruction pair a. The gas density drops The FALC extinction spectra in Fig. 6 (pdf 13) of Rutten (2019) exponentially outward [SSX (pdf 72)] but electrons dominate − show the canonical minimum at the 1.6 µm threshold of H the collision frequency by moving much√ faster than atoms. bound-free extinction where one sees deepest into the Sun (e.g., Avrett’s canonical isothermal + constant ε demonstrations Lamy et al. 1979). A nicer display is Avrett’s holistic sketch in [SSF (pdf 80)] therefore apply well to strong lines from ALC7. [SSX (pdf 36)] from Avrett (1990) plotting brightness temper- In addition, strong lines experience an atmosphere as more ature (endnote6 ) through the entire VALIIIC spectrum.67 Both isothermal than the continuum because their total τ scales are left and right of the double-hump maximum at 1 µm the extinc- compressed in height with respect to the continuum τ scale due tion and formation height increase, sampling the photospheric to their steeper buildup [SSF (pdf 38)], flattening their B(τ) gra- temperature decrease until the VALIIIC minimum above which dient [SSF (pdf 85)]. Where the temperature decays obeying RE lines and edges reverse into emission, with line crowding in as in the upper photosphere (endnote 44), lines see yet shallower the blue and UV sketched schematically. The shallow dip be- Bν gradients. √ tween the humps (continuum enlarged in [SSI (pdf 5)] and [SSF Hence line core darkening by ε scattering applies to many (pdf 55)]) spans the visible.68 The second hump is the opacity strong lines. [Back][ Main call] minimum at 4000 Å. At shorter wavelengths metal ionization edges steeply increase the summed continuous extinction [SSI 22. Menzel departure coefficients (pdf 7)]. Be aware that Avrett and Fontenla used the Menzel definition of The 4000 Å extinction minimum provides nominally less population departure coefficients that applies the next-ion pop- deep sampling than the 1.6 µm dip [SSF (pdf 55)], but effec- ulation as normalization and is described carefully in VALIII tively better because contribution functions steepen and inten- Eqs. 13 – 18 (pdf 29). I use the more intuitive Zwaan definition sity contrasts increase with the Wien response to temperature [RTSA Eq. 2.104 (pdf 53)] [SSF (pdf 73)] using the total element variations (Ayres 1989). In addition, the diffraction limit is also population as normalization.69 Between the two the values re- 6/5 four times better (but not seeing: Fried aperture r0 ∼ λ ). verse for majority species as H I, which led Fontenla et al. (2009) Hence, SST/CHROMIS 4000 Å imaging is currently the best in to misinterpret their own results (Rutten and Uitenbroek 2012). observing deepest granulation (Scharmer et al. 2019). For deep Various other authors (not to be named) reversed the meaning of imaging of intergranular MCs the SST 3950 Å wide-band filter VALIII Fig. 30 (pdf 29) into its b1 curve implying chromospheric is yet better because it samples the overlapping outer wings of hydrogen ground-state population ten times the local gas density Ca II H & K between these lines where reduction of collisional – ignoring equation (18) above the figure, haha. [Back][ Main damping by lower density produces wing opacity reduction, call] hence hole deepening in MCs and bright points with increased contrast (example with comparison to HMI in Fig. 2 (pdf 4) of 23. Epsilon is small Rutten 2020). [Back][ Main call] Another T-shirt slogan70 for this course. Its value sinks from unity in the deep photosphere to ε ≈ 10−2 − 10−6 for Na I [SSX 20. Square root epsilon law (pdf 80)], Ca II K [SSX (pdf 88)], Mg II k [SSX (pdf 89)], Lyα √ [SSX (pdf 90)] making the line source function mostly J hence This T-shirt-qualifying law says Sν(0) = εν Bν(T) at the top of an isothermal constant-ε 2-level-atom plane-parallel non- NLS. For estimating emergent intensity with the Eddington- irradiated atmosphere [RTSA Eq. 4.81 (pdf 117)], as shown in Barbier approximation you then need to know the local intensity Avrett’s [RTSA Fig. 4.21 (pdf 128)]. The emergent intensity in all directions defined by source function values elsewhere. If the η S d contribution is significant you also need to know the [RTSA Eq. 4.84 (pdf 117)] is similarly low: scattering lines get ν0 ν0 dark. intensity in all directions in other transitions with other wave- It seems a straightforward result from simple linear-looking lengths. Worse, the same may hold for evaluating the τν(h) scale Eq.2 and Eq.3 for such a simple atmosphere – but Hubený needed in your Eddington-Barbier interpretation or your line (1987b, 1987a) renamed it the Rybicki equation and analyzed synthesis results or in your “inversion” code. RT is easy only why and how it works (which I did not understand and there- when LTE holds everywhere at all participating wavelengths for fore handed out as student presentation topic in my RTSA all S and α processes or when there is no RT as in CE. course). I added two pages explanation myself [RTSA discus- The smallness of ε is also numerically challenging. The sion (pdf 117)] (which Zwaan said he did not understand). source term ε B may seem neglectable but it isn’t because with- I also taught the definition of ε wrongly in IART (now cor- out it there are no photons to scatter. In the early days of rect in [IART Eq. 7.10 (pdf 101)]) until I finally grasped that electronic computing it was tempting to stop Λ iteration [SSF the cancelation of c and g even applies to Thomson scatter- (pdf 101)] when the convergence reached tiny steps at the ma- ing. [Back][ Main call] chine precision – but they were not negligible because still far from the answer needing many more [SSF (pdf 103)]. This was √ 21. Validity of ε law remedied by operator splitting in approximate Λ iteration (ALI) [SSF (pdf 102)] [RTSA Sect. 5.3.2 (pdf 145)] and convergence Avrett’s ALC7 “chromosphere” is near-isothermal [SSX acceleration [RTSA (pdf 150)] [SSF (pdf 103)]. (pdf 72)] with near-constant ε in lines as Na I D [SSX (pdf 80)] 1 Also, at small ε the up and down radiative rates may be large from increasing hydrogen ionization and hence near-constant but nearly cancel; then using their “net radiative bracket” nor- ization rates specified from 128 up. Since CHIANTI v9 the 25 levels of malized difference can be numerically better and useful (Thomas O VIII are added also. 1960, Eq. 6.75 ff (pdf 151) of Jefferies 1968, Eq. 10 (pdf 6) of 67 I used this graph as RTSA exam by white-pasting most labels and Carlsson and Stein 2002). [Back][ Main call] ask explanation of every curve (my teaching is always graph-oriented). 69 68 We see the Sun at her smallest and sometimes just covered by our I used to write these as βl and βu following Wijbenga and Zwaan moon (endnote 15). A true solar physicist (pdf 9) observing the Earth (1972) but stopped because colleagues use b wrongly anyhow. 70 also sees deepest in our optical. Often hinted to students but no success yet. Compendium solar spectrum formation 23 Notes

the ultraviolet (Fig. 1 (pdf 5)). At smaller temperature sensitivity 24. No Thomas photoelectric control larger Texc excess is needed for the imposed bu/bl ratio. Hα has a sizable ηS d contribution from bound-free detours as The lefthand detour loop in Fig. 2 (pdf 4) of Rutten (2019) the loop shown in Fig. 2 (pdf 4) of Rutten (2019). It was fa- is for Ca II. Similar transcription of the Wien sensitivity of the mously called “photoelectric control” by Thomas (1957) with PRD-split source function of Ca II K to Ca II 8542 Å is evident a schematic source function diagram in Fig. 3 (pdf 5) of Jefferies for the wing source functions (label g) of Ca II K in panel c and and Thomas (1959) that was reprinted in Fig. 12-9 of Mihalas of Ca II 8542 Å in panel g of Fig. 2 (pdf 7) of Uitenbroek (1989). (1970) and Fig. 11-11 of Mihalas (1978). It looks similar to the Back Main call l [ ][ ] Hα source function in [SSX (pdf 76)] with Sν leveling out in the photosphere higher than Bν in the temperature minimum, indeed 26. Photon loss, no nanoflame heating, photon suction similarly to the Balmer continuum [SSX (pdf 76)]. “Photon loss” means losing the energy content of a photon alto- However, this “photoelectric control” designation was wrong gether from the object (star, atmosphere, feature), as “neutrino because even Hα is mostly scattering as shown in Fig. 12 (pdf 9) losses” means energy losses from a stellar core. In both the LTE of Rutten and Uitenbroek (2012); the upper-photosphere excess and CE extremes the energy loss is direct because the emer- is due to backscattering from the chromosphere. Such excess gent photon was locally created collisionally (bound-bound pair does not occur in the RE atmosphere at left in Fig. 7 (pdf 6) d or in LTE e, their bound-free equivalents, free-free). Where of Rutten and Uitenbroek (2012) whereas with “photoelectric d the Maxwell distribution holds this is direct loss of local thermal control” by ηS it should since the ultraviolet continua in Fig. 5 energy. ≈ (pdf 5) there have similar S J > B excess with or without a For LTE it should be a negligible leak (otherwise LTE is not chromosphere. It was a red herring. There are some lines as the Mg II triplet around Mg II h & k valid); an example is forbidden Mg I 4571 Å from ALC7 with S l ≈ B τ ≈ − for which source function enhancement in cooling recombining out to log( ) 4 in [SSX (pdf 82)] (but not its opacity gas (possibly non-E) producing emission is likely, but the many which is heavily depleted by ultraviolet overionization). Thomas-citing typecastings of Hα as “photoelectric” were mis- For CE summing all emergent photons from an EUV-bright feature quantifies radiative losses that measure required heating led. Hα is instead sensitive to non-E opacity boosting in cooling√ recombining gas (endnote 67) that controls the location of its ε for a stable structure that obeys SE (otherwise CE is not valid). source function decay (endnote 55). [Back][ Main call] Hence CHIANTI-computed radiative loss curves following clas- sic Fig. 2 (pdf 6) of Cook et al. (1989) specify required coronal 25. Photon pumping, lasering, Wien transcription heating. It is a mistake to count summed brightness of short-lived bright features as contribution to local heating: to the contrary, Photon pumping is multi-level NLW interaction where one tran- this energy is radiated away, only the remainder of the feature sition may feed another photons and Wien sensitivity. Classic heating may count. Observing small flaring flames (or “camp- nighttime cases are planetary-nebulae and nebulium lines ([SSX fires” endnote 66) spoils their nanoflare heating candidacy. (pdf 142)] ff). The two detour loops in Fig. 2 (pdf 4) of Rutten For a two-level scattering line (initial d or e followed by (2019) are good solar examples. many b or c steps with small chance of collisional destruction) At right the Balmer continuum overionizes hydrogen from the actual photon losses occur near the mathematical τ = 0 sur- level n = 2 in the upper photosphere because it escapes much face, in the last scattering around depth τ = 1, but the losses deeper and at short wavelengths where the large Wien sensitiv- start being felt already near the thermalization depth. There the ity produces Jν > Bν as in [SSX (pdf 76)] for the ALC7 star. radiation is nearly isotropic (isotropic in the isothermal simpli- Hence the ALC7 proton density is “pumped up” with respect to fication) and becomes more and more anisotropic further out, the local SB value (the n=2 level has b ≈1 because Lyα is ther- 2 ending up purely outward√ directed at the surface. The scatter- malized (lefthand plot [SSX (pdf 76)]). This photon excess con- ing decline towards the ε value there may be seen as inward- tributes together with Brackett, Paschen and Balmer-line photon diffusing photon losses. The same holds within effectively thick losses (suction, endnote 26) into downward overpopulation flow features rather than a plane-parallel atmosphere. through the H I Rydberg levels (Fig. 3 (pdf 6) with explanation For two-level atoms the losses deplete the upper level with in Rutten and Carlsson 1994). Hotter features in the deep photo- larger fractional bu decrease than the corresponding fractional bl sphere cause larger hydrogen overionization higher up and also increase of the lower level due to the Boltzmann ratio between sideways in cooler surroundings. l them: in two-level atoms the lower level is the main population The S /B ≈ bu/bl ratio may so be controlled in the Wien reservoir with n1  n2 and therefore close to LTE because there domain and transcribe to the Rayleigh-Jeans domain where is nowhere else to go: nearly the whole abundance sits in n = 1. given departure ratio from unity has larger effect as plotted The b curves in the first plot of [SSX (pdf 89)] for Mg II k in in [SSF (pdf 74)] from [RTSA Fig. 2.8 (pdf 54)]. Blow-up and ALC7 are the textbook example. These ions are not two-level then negative lasering is reached at right (line extinction and atoms but close to it because Mg is predominantly once-ionized source function both become negative for sufficient bu/bl [SSF throughout the ALC7 photosphere and chromosphere, as shown (pdf 74)] [RTSA Eq, 2.116 (pdf 55)] but emissivity remains pos- by the fractional lower-level population curve in the same plot. itive). Such blow-up (but not yet lasering) is reached in the “Photon suction” is a multi-level detour phenomenon, done Mg I 12 micron lines (endnote 34). l by strong lines when the population reservoir is in the next stage, For Fe II lines long-wavelength S /B blow-up from pump- as for most neutral metals in the photosphere. It is the reverse of ing in the ultraviolet with transcription of the exponential Wien photon pumping in which incoming photons overexcite higher sensitivity holding there to linear Rayleigh-Jeans sensitivity is levels from a lower population reservoir (endnote 25); in photon illustrated in Fig. 5 (pdf 6) pf Cram et al. (1980) and in Fig. 3 suction outgoing photons draw population from a higher reser- (pdf 10) of Rutten (1988). In the latter the high peak of the ex- voir down into lower levels. The alkali resonance lines are the citation temperature for Fe II lines at 7000 Å arises because it prime example.71 Figure 6 (pdf 8) of Bruls et al. (1992) shows must accommodate the imposed bu/bl ratio, here bu overpopu- 71 lation excess because bl ≈ 1 for dominant-stage ground states. Solar low-l history: the alkali resonance lines and some others The bu overpopulation results from radiative Jν > Bν feeding in are suited to build gas-filled resonance cell detectors in which atoms Compendium solar spectrum formation 24 Notes ground-level departure curves b1 (still called β following Wij- an early precursor of Carlsson and Stein (2002) indeed hold- benga and Zwaan 1972) of K I and Na I increasing slightly until ing for spicules (endnote 62); the height where ultraviolet underionization causes steep rises. – Chae (2014): second cloud model for the incident profile; This slight increase was surprising because neutral metals suf- – Heinzel et al. (2015): multi-thread RT with mutual irradia- fer ultraviolet overionization in the upper photosphere by sam- tion. pling the hot Balmer continuum coming from below, evident in the ALC7 population panels for Mg I 4571 Å in [SSX (pdf 82)] For the incident profile often a field average or a nearby sub- and Fe I 6301.5 Å in [SSX (pdf 83)]. Such pumping may happen field average is used but this is wrong in the case of backscat- from excited levels, but the alkali ground-state ionization edges tering as observed for Hα filaments as sub-filament brightening towards the limb by Kostik and Orlova (1975), proposed for Hα at 2856 Å for K I and 2411 Å for Na I suffer it also. For them fibrils by Al et al. (2004), and occurring in 1D standard mod- the overionization is fully compensated by population replen- els (endnote 24). The near-equality of the dashed profiles in the ishment from the low-lying next-ion population reservoir, drawn lower panels of their Fig. 8 (pdf 7) made Rutten and Uitenbroek from there in a net recombination flow along high levels that is 72 (2012) suggest as recipe (pdf 9) to use the outward intensity pro- driven by the photon losses in the resonance lines . If one starts file in a RE model at optical depth equal to the cloud’s thickness a NLTE modeling iteration with LTE populations the ultraviolet as impinging background profile [SSX (pdf 3)]. This recipe im- Jν > Bν excesses deplete the ground state but the NLTE scat- proves the finding of Bostanci and Al Erdogan (2010) that a RE tering losses diffuse inward and draw recombination population emergent-intensity estimate works best. [Back][ Main call] until SE is reached with full suction compensation for the alkali’s with their low ionization energies. A two-level-plus-continuum 29. Na I D, Mg I b, Ca II 8542 Å heights demonstration is shown in Fig. 2 (pdf 4) of Rutten and Carlsson (1994), with explanation. It was made by Bruls(1992) to adorn [Spoiler: million difference in line identification [SSF (pdf 41)].] the cover of his PhD thesis (not on ADS). [Back][ Main call] The same Fig. 4 (pdf 5) of Uitenbroek and Bruls (1992) is also the subject of less simple but recommended [RTSA Problem 7 27. Coronal clouds (pdf 250)] on Na I D2 line formation in the FALC star (Sect.6 ) [Spoiler: my answers to most RTSA problems]. The Na I D lines The simplest clouds are fully transparent as CE assumes for are the deepest in the visible (Fig. 4 (pdf 17) of Rutten and van coronal structures, e.g., loops. No RT whatsoever. In optically- der Zalm 1984b, Fig. 10 (pdf 18) of Rutten 2019)73. The Na I D thick RT terms: the homogeneous thin cloud solution [ISSF lines are often called chromospheric. Your Eddington-Barbier Eq. 27 (pdf 13)] has no impinging intensity from behind (I (0)≈ ν estimate h ≈ 500 km with [SSF (pdf 40)] puts the Na I D for- 0) in CE hence I (D) ≈ S τ (D) = j D with D the geometrical 2 ν ν ν ν mation just above the temperature minimum of the FALC star, cloud thickness along the line of sight. This holds for Thom- √ but the line source function shows a ε scattering decline as son scattering producing the white-light off-limb corona (with if the FALC star has no non-RE heating whatsoever. The pho- S ≈ J , endnote 15) and also for locally thermally created EUV ν ν tons emerging there were made much deeper, near thermaliza- radiation. In terms of Eq.2 for the latter (applicable if we count d d tion depth in the photosphere, and scattered out from there ig- escaping photons into αν with ην ≈ 1 for CE) both Jν ≈ 0 and d l l a noring the temperature rise and then escape. Intensity varia- Sν ≈ 0 so that Sν ≈ εν Bν(T) or jν ≈ α Bν(T) with exceedingly a tions of the core therefore correspond to photospheric temper- small α and εν due to absence of collisional deexcitation. ature variations, not at the τ ≈ 1 Eddington-Barbier height. This Even though all line photons are thermally created the inten- scattering is quantified for Na I D1 in the similar ALC7 star in sity stays far below saturation to the Planck function for lack of [SSX (pdf 85)] where the S, B, J graph shows optical depth lo- local enclosure. The corona is far hotter than the photosphere cations, destruction probability, thermalization length and de- and its radiation is thermal and defined by its temperature just as coupling height ([SSX (pdf 80)], endnote 45). The Na I D lines for the optical continuum from the photosphere [SSI (pdf 4)]– are the deepest because they are the strongest scattering lines of but in contrast to the photosphere the corona radiates far below which the emergent intensity is fully controlled within the photo- its thermal capacity because it lets all new photons depart instead sphere. However, in measuring Doppler or Zeeman modulation of boxing them in. [Back][ Main call] the Eddington-Barbier estimate for the inner wings may be ap- propriate because these signals are encoded at the last scattering. 28. Sophisticated clouds Similar FALC results in Fig. 13 (pdf 11) of Rutten Cloud complexity was added by e.g., et al. (2011) compare Na I D1 formation with Mg I b2 and – Molowny-Horas et al. (1999): model look-up tables for “in- Ca II 8542 Å formation. Mg I b2 forms similar to Na I D1, decou- version”; pling slightly higher, but the source function of Ca II 8542 Å de- – Heinzel et al. (1999) and Mein et al. (1996): source function couples above the FALC temperature minimum so that this line variation along the line of sight. This was already done by does convey FALC-chromosphere signature. The same is seen Beckers (1968) in a famous spicule review by using NLTE for the ALC7 star in [SSX (pdf 27)] and in comparing [SSX population tables for hydrogen from Giovanelli (1967) – al- (pdf 85)] [SSX (pdf 84)] [SSX (pdf 87)]. though warning against assuming SE below Eq. 10 (pdf 32), Rutten et al. (2011) used high-quality SST observations to compare the three lines in intensity in Fig. 4 (pdf 5) in compar- resonance-scatter optical sunlight, either irradiance for full-disk profile ison with simulation results in Fig. 5 (pdf 6) and those smeared sampling following the initial application for Sr I 4607.3 Å of Blamont and Roddier (1961) or with imaging using magneto-optical encoding 73 History: the first was drawn with a Calcomp plotter, a mechanical following Öhman (1960, not on ADS) and Cimino et al. (1968). Us- device moving a pen with an ink cartridge transversely across a rotat- ing resonance-cell Dopplershift gadgets worldwide became a large in- ing drum with paper chart from a large roll, both motions computer- dustry in g-mode-searching low-l helioseismology but a romantic over- controlled in 0.01 inch steps. Earlier the entire university computer (one promise/overclaim saga rather than a success story as endnote 53. for all) was required for such plotting. For job scheduling (one at a time) 72 The compensation was earlier noted in Fig. 1 (pdf 3) of Gomez et al. I had to count all steps. The white-coated operators liked plotting these (1991) but there incomplete from lacking high levels in the model atom. crowded spectral charts for giving them long breaks. Each strip was a meter long and photographically reduced to photo-offset print size. Compendium solar spectrum formation 25 Notes to observation-resolution in Fig. 6 (pdf 7), and then observed and Fig. 12 (pdf 14) of Carlsson et al. (1992) varying many transi- simulated line-core Dopplergrams in Fig. 7 (pdf 8). These figures tions in their Mg I model atom one by one (endnote 34), the confirm non-chromospheric sampling by Na I D1 and Mg I b2 successive step-function temperature perturbation to quantify whereas Ca II 8542 Å does show chromospheric signatures (end- UV intensity response of Fossum and Carlsson (2005b) and note 67). The main difference between the first two is that the the NLTE-nonE versus LTE simulation-run comparison in Fig. 4 Na I D lines have higher Doppler sensitivity than the Mg I b lines (pdf 7) of Leenaarts et al. (2007)( Sect.7 ). by having narrower cores with steeper flanks. This explains the Scorn-avoiding moral: even if your nice simulation repro- remarkable difference between their wings in sampling upper- duces actual observations you still have to delve per sensitivity photosphere reversed granulation in Fig. 1 (pdf 2). analysis into your simulation to diagnose how and why, as diffi- The steep Na I D flanks also make Na I D Dopplergrams a cult as diagnosing observations but with the advantages that you sharp proxy for localizing upper-photosphere shocks in MCs, as can check below τ = 1 and know the inputs. Carlsson and Stein summarized in [SSX (pdf 54)].74 [Back][ Main call] (1994) (more formally Carlsson and Stein 1997) remains the canonical example: intricate observed patterns, convincing re- 30. Contribution, response, sensitivity production with data-driven simulation (endnote 49, four-panel When using “the formation height” you may earn my scorn (like breakdowns (endnote 50). [Back][ Main call] the thinnies of endnote 34). Only thin clouds at specific height (as in endnote 27) share their height between different formation 31. MULTI program quantities. For thick objects a common characterizing escape The NLTE spectral synthesis code of Carlsson (1986) is obvi- height may be applicable when LTE holds for all contributing ously named for going beyond the 2-level approximation. It ini- processes, but for scattering lines the intensity is set by temper- tially used linearization following Scharmer and Carlsson (1985) ature much deeper and more widely within the atmosphere or with the brilliant Scharmer operator [RTSA Eq. 5.52 (pdf 148)] within an effectively thick structure than the Doppler and Zee- (yes, the same as the brilliant SST builder) to obtain Iν rather man response to line-of-sight velocities and magnetism (as in the than Jν from Sν with a local two-way along-the-beam Eddington- Na I D example of endnote 29). Also, for weak lines the emer- Barbier approximation. It was later extended with the faster com- gent intensity is dominated by background continuum processes bination of the diagonal [RTSA OAB operator (pdf 150)] of Ol- rather than line processes (exemplified by their factor-million son et al. (1986) and [RTSA convergence acceleration (pdf 150)] height ratio in [SSF (pdf 41)] with Schuster-Schwarzschild cloud of Ng (1974).[Back][ Main call] modeling appropriate for the telluric blends). In plane-parallel atmospheres the vertically emerging inten- 32. MULTI3D program R sity is Iν = Sν exp(−τν) dτν with integrand or “contribution Extension of MULTI to 3D parallel computation by Leenaarts function” Sν exp(−τν) sketched in [SSF (pdf 29)], an extended and Carlsson (2009). It uses short characteristics with further- plume. The photons do not “escape from τν = 1” but through- away contributions corrected by iteration. [Back][ Main call] out this plume. Such contribution plumes are plotted against height, i.e., dI/dh = jν(h) exp(−τν(h)), for many wavelengths 33. Classic abundance determination in the VALIIIC spectrum in fabulous many-page Fig. 36 of Ver- Holweger used the ALGOL “Kieler program” of Baschek et al. nazza et al. (1981) with excerpts in SSX. The first one in [SSX (1966) for LTE spectral synthesis in a classical approach: impose (pdf 61)] for the 500 nm continuum shows the expected con- a T(h) “model” relation, evaluate density stratifications for the tribution shape with initial jν(−h) rise set primarily by near- given element mix by assuming hydrostatic equilibrium [RTSA exponential inward density increase until stopped by the steep Sect. 7.2.3 (pdf 166)], compute continuum extinction from H− exp(−τ) decrease and with I(1, 0) ≈ S (log τ = 0). The last one in and extinction of a given line per SB (bu =bl =1 in [RTSA Eq. 9.6 [SSX (pdf 67)] for λ = 40 nm shows weird dI/dh shape without (pdf 224)]), use the integral solution of the transport equation peak near log τ = 0 and two orders-of-magnitude deviation from [RTSA Eq. 2.43 (pdf 38)] to obtain the emergent profile. He used the Eddington-Barbier value. line areas (“equivalent width” [RTSA Sect. 9.1.2 (pdf 224)]) as Earlier, sophistication beyond this simple integrand was quality gauge, repeated this for many optical lines, and derived based on weighting functions (e.g., Pecker 1951) and de- a best-fit temperature and corresponding density stratification by riving such for line depression rather than emergent inten- trial-and-error manual T(h) adjustment. sity (Gurtovenko et al. 1974). The most informative format His many followers in classic abundance determination then became Carlsson’s integrand breakdown rewriting dI/dh into took these stratifications (“model atmosphere”) for granted and [Sν][τν exp(−τν)] [d(ln(τν))/dz] in Eq. 6 (pdf 15) of Carlsson did similar line fitting for different lines of a given element and Stein (1994) and plotting these three terms and their product to establish its abundance, with much debate about the re- as function of height and wavelength through the line of interest liability of oscillator strength (transition probability) g f val- (Fig. 8 ff (pdf 16), endnote 50). ues. [Back][ Main call] In the meantime “response functions” were formulated that estimate the effect of perturbations of a given quantity on an 34. Mg I 12-micron emission lines emergent line profile (classics Mein 1971, Beckers and Milkey Assuming LTE implies predicting self-reversing chromospheric 1975, Magain 1986) but nowadays this is done by numerical cores in e.g., the Na I D lines treated in practical SSB 3 (pdf 22) sensitivity analysis, for example the “multi-Multi” analysis in but also in the strong Fe I lines that Holweger employed. Their absence was attributed correctly to NLTE scattering by oth- 74 Repeating the corollary on page 15 (pdf 15)) of Rutten et al. (2011) that Fourier analysis of full-disk Doppler sampling with the proposed ers but Holweger insisted tenaciously that NLTE computer pro- GOLF-NG space instrument (Turck-Chièze et al. 2012) will not achieve grams fatally overestimate NLTE departures by not including chromospheric helioseismology as advertised, whereas hoped-for g- sufficient hydrogen-atom collisions [SSX (pdf 39)] and that the mode detection will be hampered by upper-photosphere reversed gran- actual chromosphere starts higher. When he refereed Carlsson ulation, fluxtube shocks and “shutter” modulation by overlying features et al. (1992) explaining the enigmatic Mg I emission lines at as in Figs 5 and 6 (pdf 6) of De Wijn et al. (2007). Compendium solar spectrum formation 26 Notes

12 micron this was also his main complaint, but the Rydberg gradient gets too shallow and the same as for locations with recombination modeled there was already collision-dominated. less steep temperature decline. Different surface structures with The formation of these striking features is interesting. Nat- different gradients (as granules and lanes) can so be fitted erro- urally they were first attributed75 to chromospheric formation neously by a single shallow gradient. (Brault and Noyes 1983) but they are photospheric, combining Additionally, the neglect of source function NLTE for higher- NLW photon pumping in photo-ionization (endnote 25), NLW up scattering in stronger lines results in undoing an actual chro- photon suction (endnote 26) and a diffusive radiative-collisional mospheric temperature rise, masking that also. Rydberg population departure flow. In the Rydberg top this is Fig. 8 (pdf 8) of Rutten and Kostik (1982) shows how dominated by collisional decays, the opposite of the coronal di- a HOLMUL-like model so results from the steeper-gradient electronic radiation-down case of endnote 10. See the cartoon HSRA model of Gingerich et al. (1971), summary in [SSX explanations in Rutten and Carlsson (1994) including Mats’ (pdf 40)]. A cartoon is shown in Fig. 4 (pdf 14) of Rutten “Rydberg flows for kayakers” in Fig. 4 (pdf 7)(better quality). (1988).[Back][ Main call] Carlsson et al. (1992) obtained excellent fits in Fig. 6 (pdf 10) with a 1D RE model atmosphere and likewise in Fig. 14 (pdf 15) 37. Milne-Eddington approximation with the non-RE MACKKL model. These fits are uncommonly Inverters often apply the further simplification of assuming good because these lines form in the upper photosphere where Milne-Eddington line formation, i.e., constancy of the line-to- RE reigns and 1D models apply best (endnote 44). continuum extinction ratio with height. This assumption was Actually, Lemke and Holweger (1987) had come close to tested with a MURaM simulation and LTE line synthesis for model these emission lines properly but their Mg I model atom comparable weak Fe I lines at different excitation energy in did not reach high enough up the collision-dominated Rydberg Fig. 8 (pdf 10) of Vitas et al. (2009); summary in [SSX (pdf 41)]. ladder. Their paper was wildly and wrongly attacked by Zirin It is bad at low excitation, better at higher excitation as al- and Popp (1989), in turn severely rebutted by me in Carlsson ready evident in Fig. 3 (pdf 16) of Rutten and van der Zalm et al. (1992). This impressed Zirin so much that he saw me (1984b).[Back][ Main call] as potential successor. I saw him as thin because in the ab- stract he stated that these Mg I lines are thin. They are not in 38. Interlocking Fe I, Ce II, Fe II lines his sense, but that itself is a no-no because lines cannot be −1 Grotrian diagrams of e.g., Fe I and Fe II are rich in bound-bound thin (define a halfwidth threshold in Å or km s ?). A solar- transitions. Optical multiplets are mostly high up but connected atmosphere structure can be thin or thick or effectively thin or by ultraviolet transitions to levels lower down (ground state or thick in a line (maybe the Sun is thin in neutrino lines) but no low-lying metastable levels), comparable to the Ca II triangle line is thin or thick. Never say “this line is thin” or earn my in Fig. 2 (pdf 4) of Rutten (2019) but in large number. At small scorn! [Back][ Main call] lower-level excitation energy these ultraviolet transitions tend to be strong, maintaining S l ≈ B further out as in the righthand car- 35. Plane-parallel atmosphere files toon of [SSF (pdf 111)] and forcing bu ≈ bl on their upper levels I copied files with plane-parallel model atmospheres from the at the deeper heights where their weaker optical siblings reach 1980s that I found I still have into my solar file archive started τ ≈ 1. When the latter lines so share strong-line bl in their upper while writing endnote8 . These include HOLMUL, the HSRA and lower levels they have LTE source functions (but their opac- model of Gingerich et al. (1971), the thesis model of Lites ity scales with bl and senses ultraviolet overionization as in the (1972), the VALIIM model of Vernazza et al. (1976), the lefthand cartoon of [SSF (pdf 111)]). radiative-equilibrium model of Bell et al. (1976), the models Other interlocking lines occur in lanthanide rare earths. They VAL3A– VAL3F of Vernazza et al. (1981) and the MACKKL have multiplets with many members that overlap closely in model of Maltby et al. (1986). The FALC model of Fontenla energy. Ce II has most lines, peaking at 10 lines/Å around et al. (1993) is supplied and analyzed in practical SSB1. All Ca II H & K. These are weak lines but scatter strongly with much four FAL models of Fontenla et al. (1993) are available in RH crosstalk. They show up in emission inside the limb both out- (endnote 47) input format in the RH github Atmos directory, but side and within the extended wings of Ca II H & K which was my new solar file archive adds some more including VAL3C, explained in the classic analysis of Canfield (1971b) following MACKKL, models of Fontenla et al. (2009) and ALC7 of Avrett on Canfield (1971a) and his earlier thesis studies. They appear and Loeser (2008).[Back][ Main call] bland from wide multi-wavelength scattering producing their J, whereas similar Fe II blends in the wings of H & K (Rutten and 36. NLTE masking Stencel 1980) show much spatial variation due to pumping by For a steep outward temperature decline the outward increas- particular ultraviolet transitions (Cram et al. 1980) and so gain- ing ultraviolet J > B excess from scattering and Λ results ing ultraviolet Wien sensitivity (endnote 25). [Back][ Main call] in NLW-type NLTE overionization of minority-species metals [SSF (pdf 111)]. Instead adopting Saha ionization then puts the 39. Line haze modeling height assigned to τ = 1 too high so that the fitted temperature Photospheric species as Fe I and Fe II 76 are immensely rich in 75 History: being unusual they were blamed on some instrumental de- bound-bound transitions with dense line crowding into the “veil- fect and white-pasted by Goldman et al. (1980) in their infrared spec- ing line haze” of Labs and Neckel (1972) and Greve and Zwaan trum atlas. Brault and Testerman had already spotted these intriguing (1980), densest in the violet and ultraviolet (first extinction curve features (Fig. 1 (pdf 2) of Brault and Noyes 1983) at Kitt Peak; Brault in Fig. 3 (pdf 13) of Rutten (2019). It is discussed as major ob- therefore phoned them to ask why their atlas showed gaps precisely at stacle in Sect. 6.1 (pdf 17) of Rutten (2019) together with various these locations. Answer: “everybody knows that Fraunhofer lines are recipes to treat it in spectral synthesis, varying from the sim- in absorption”. Murcray et al. (1981) (including two more Murcrays) ple Zwaan-inspired Fig. 2 fudge (pdf 3) of Bruls et al. (1992) to then sensed their gaffe, dewhitened them and claimed discovery (with- out mentioning Brault) with fancy attribution to Si VIII and Mg VIII but 76 The classic Grotrian diagrams of Moore and Merrill (1968) are only not daring to exclude telluric origin. They were identified as solar Mg I partial, with their Fe I (pdf 66) and Fe II (pdf 65) charts showing only n=7–6 transitions by Chang and Noyes (1983). major systems. Complete ones would be black all over. Compendium solar spectrum formation 27 Notes brute-force solving all NLTE population equations for all transi- the “International Study week” at the Bilderberg Hotel (famous tions in Fontenla et al. (2015). The RH code (endnote 47) offers from more famous Bilderberg conferences) that De Jager orga- both these and also LTE and equivalent two-level scattering syn- nized in order to jointly derive the Bilderberg Continuum Atmo- thesis for all blend lines not treated in NLTE detail. The latter sphere of Gingerich and de Jager (1968) that later got improved two are compared with observations in Fig. 10 (pdf 18) of Rut- to the HSRA of Gingerich et al. (1971). As student I operated ten (2019), showing that both fail. The Pandora code uses an the slide projector and later described this truly expert meet- l l imposed gradual transition from S ν = Bν to S ν = Jν with SB ing in Rutten (2002) with a group photo (pdf 5) sporting Hol- opacities (endnote 43). At the end of the section (p. 19 (pdf 19)) weger (more photographs at the bottom of my astronomershots). I suggested a recipe solving for NLTE populations in a represen- E.A. Müller commented at its conclusion that she wouldn’t use tative simple atom and applying the resulting departures to all the new model with NLTE but stick to Holweger’s model and blends. [Back][ Main call] LTE, then did so with HOLMUL (Holweger and Müller 1974). Nevertheless, the Auer–Mihalas NLTE results did upset Un- 40. Fe I 6302 Å doublet, Mg I 4571 Å söld’s LTE bible “Physik der Sternatmosphären” (Unsold 1955, + A good example is the Fe I 6302 Å polarimetry doublet. The 900 citations). The translation by H. Panofsky and A.K. Pierce ALC7 star [SSX (pdf 83)] has a near-LTE source function from was never printed; I recollect that Unsöld insisted that any as- ultraviolet interlocking (endnote 38) and large opacity deficit trophysicist could and should learn German to read his book. I from ultraviolet overionization. Their NLTE “inversion” fitting have it (and can read it) but haven’t opened it in years. So it should use RH (endnote 47) with a large Fe model atom, sizable goes. [Back][ Main call] electron-donor atoms and a suited line-haze recipe or it will un- derestimate their source functions from erroneous scattering and 43. Pandora program overestimate their opacities from erroneous ionization. NLTE spectral synthesis code of Avrett and Loeser (1992); more Another good example is Mg I 4571 Å which is a forbidden here and here and here. Similarly to Holweger’s method it ap- transition and therefore has S = B up to large height but suf- plies hydrostatic plane-parallel equilibrium to a trial manually- fers considerable extinction NLTE even in 1D modeling [SSX fitted T(h) temperature stratification but uses SE NLTE for radi- (pdf 82)]. “Inversion” fitting needs a large Mg model atom plus ation evaluation in a giant “equivalent two-level” iteration loop the other electron donors and a line-haze recipe. [Back][ Main treating each bound-bound or bound-free transition that is ex- call] plicitly taken into account (not added as background opacity) with two-level rate equations, letting iteration take care of multi- 41. Granulation contrast level crosstalk. The equation system is large enough to earn the At the time Å. Nordlund had a habit of presenting his revolu- ominous name. tionary work in non-page-limited proceedings (e.g., non-scanned Whereas Holweger concentrated on fitting observed optical Nordlund 1985a, 1985b). In scanned ADS-available Nordlund lines, mostly from Fe I, Avrett concentrated on fitting observed (1984b) he argued that the large difference in granulation con- disk-center continua, especially in the ultraviolet. For the line trast between his numerical simulation and observations was to haze (endnote 39) Avrett added increasing numbers of lines from be blamed on the latter, becoming the prototype of the numeri- the ever-growing list of Kurucz (2009) and treated them with an cal simulator habit to blame observed data if these do not match imposed ad-hoc gradual source function transition from S = B the wonderful and necessarily correct computation77 – but ac- in the model photosphere to the continuum S = J in the model tually Nordlund was right: Scharmer et al. (2019) showed that chromosphere to avoid non-observed core reversals that other- telescopes that do not reach 10% rms contrast around 5000 Å wise show up above the temperature minimum in LTE sampling are not doing their job while the actual contrast is higher. (p. 243 ff (pdf 15) of Avrett and Loeser 2008). Seems reason- Nordlund stems from B. Gustafsson at Uppsala who also sug- able, but using LTE opacities does not for over-ionized minority gested RT theory as venue to G.B. Scharmer at Stockholm who species. then became thesis adviser to Gustafsson’s student M. Carlsson Fontenla diverged from Avrett after Fontenla et al. (1993), (endnote 31). Nordlund and Carlsson each began fruitful col- writing a Pandora alternative and adding tricks to enhance out- laborations with R.F. Stein at East Lansing who has Scan- put similarity to the non-1D non-static non-SE spectra of the dinavian roots. This illustrious Nordic quartet eventually co- actual Sun. First he invoked the Farley–Buneman instability to authored Carlsson et al. (2004). M. Asplund is also a Gustafsson warrant adding ad-hoc artificial pressure (the Avrett stars use ob- pupil.78 [Back][ Main call] served microturbulence (pdf 19) for turbulent pressure); later he mucked gravity ad-hoc to extend the model photosphere to re- 42. LTE – NLTE schism produce dark CO lines that probably have 3D(t) NLT formation (Uitenbroek 2000). History: I wonder what Unsöld, who at Kiel had told Holweger I prefer the ALC7 star of Avrett and Loeser (2008) for my to produce a thesis proving that LTE suits the solar spectrum, spectrum-formation demonstrations ([SSX (pdf 72)] ff). It isn’t made of the Auer – Mihalas results, representing extremes of the the Sun but it is self-consistent within its RTSA physics and then rampant LTE-versus-NLTE debates. In 1967 Unsöld wrote well-suited to showcase this physics. to Minnaert that Holweger had definitely proven the LTE ve- Pandora is web-available but documented only in many racity of the solar spectrum, in a recommendation for attending pages of coding requests and specifications to R. Loeser that 77 Generally I prefer simulations to not reproduce observations so that reside in Avrett’s CfA office. P.Heinzel holds copies in On- there is evidently something to learn without temptation to claim cor- drejov. [Back][ Main call] rectness from happenstance data reproduction (endnote 55). 78 History: Gustafsson’s astronomer production reminds me of 44. RE upper photosphere C. Zwaan (e.g., me, H.C. Spruit, A.A. van Ballegooijen, C.J. Schrijver, The reason for the near-equality of the upper photosphere in all R.G.M. Rutten, L.H. Strous, also PhD theses of A. Greve, K.L. Harvey- 1D models since MACKKL to the older HOLMUL and RE mod- Angle, K.F. Tapping). Utrecht University has a great past but no present els is that this is the most homogeneous domain of all. The gran- in solar physics. So it goes. ular convection has stopped. Acoustic waves do not yet shock. Compendium solar spectrum formation 28 Notes

Magnetic fields are still mostly confined to slender fluxtubes at I supply some plane-parallel 1D “standard” model atmo- small density in the modeled “quiet” areas. Gravity waves do not spheres in RH format beyond the github-supplied FAL models yet couple into them (endnote 51). of Fontenla et al. (1993) in my solar file archive( endnote 35). In addition, the bulk of the escaping solar radiation is in For the ALC7 spectrum-formation graphs starting at [SSX the optical and forces radiative equilibrium (RE). In “grey” (pdf 73)] I used RH 1D version-2 with H, He, Si, Al, Mg, Fe, RE [RTSA Sect. 7.3.2 (pdf 173)] the continuum gets gradient Ca, Na, and Ba active, C, N, O, S, and Ni passive, and with S (τ)∼(1+1.5 τ) [RTSA Eq. 7.43 (pdf 176)] because for that gra- 20 mÅ sampling of 343 000 lines between 1000 and 8000 Å dient Λ(S ) ≈ S [SSF (pdf 45)]. The continuum in the optical is in the atomic and molecular line list of Kurucz (2009). Fig. 6 − dominated by H which has near-LTE Sν ≈ Bν in the bound-free (pdf 13) of Rutten (2019) gives an overview of their extinction range and strict LTE Sν = Bν in the free-free range. Most lines in the FALC star. The line haze is prominent between log λ = 3 in the optical also have Sν ≈ Bν (endnote 38). Hence the tem- and 4. [Back][ Main call] perature gradient is set LTE-wise by the RE condition [RTSA Eq. 7.33 (pdf 174)]. 48. Quiet-Sun Hα scenes Indeed, more modern Bifrost (endnote 59) simulations ob- Note the deep population chasm in the dashed curve in the first tain about the same average temperature decay with relatively panel of [SSX (pdf 91)], first described by Schoolman (1972) small spread at these heights [SSX (pdf 59)]. Profile fitting with and compared with Ca II H in [SSF (pdf 61)] from Leenaarts 1D models is least suspect for lines formed there, as the Mg I et al. (2006). Observationally the gap is evident when you shift 12 micron lines of endnote 34.[Back][ Main call] the SST/CRISP Hα passband from line center to a wing for a quiet-Sun target. Up to about ∆λ = ±0.5 Å the scene usually ALC7 line formation plots 45. shows fibrils extending from network and then you suddenly The 17 triple-graph displays [SSX (pdf 78)] ff show populations, drop into deep-photosphere granulation, skipping the middle and source functions, emergent profiles for 13 representative lines of upper photosphere. Only in rare utterly quiet areas you may see the ALC7 star with common intercomparable graph formats. somewhat higher-up reversed granulation near ∆λ=±0.5 Å. The first graph [SSX (pdf 79)] shows lower and upper level I offer two SST examples. The first is a very quiet tar- population departure coefficients bl and bu. The other curves are get corner in the data of Rouppe van der Voort et al. (2007) NLTE and LTE fractional populations. Their divergence follows in this Ca II H + Hα blue-wing movie where reversed Hα at the log(b ) departure from zero (Eq.4 ). l ∆λ = −0.45 Å shows reversed granulation as grey underlying The second graph [SSX (pdf 80)] shows B, J and S l (these pancakes. must always be shown all three). They are on equivalent temper- The second is a similarly very quiet SST target (in a disk- ature scales, deWiened to enable comparisons between different center coronal hole) that you may inspect yourself by doing my lines across the spectrum (endnote6 ). On logarithmic intensity alignment practical( endnote 64). Its showex commands for its scales the S l separations from B would equal the logarithmic SST data enable Hα wing blinking against reversed granulation separations of bu from bl in the first graph. The thick solid curve usually lies between the other two, identifying it as S l obeying in the Ca II 8542 Å wings to identify it in Hα.[Back][ Main S l ≈(1−ε) J+ε B. For PRD lines it is split between monofrequent call] samplings. The log(τ) marks are defined by the lower-level frac- tional population, abundance and density drop. The other two 49. Carlsson-Stein shocks curves measure scattering governing S l − B split. In Carlsson and Stein (1994) they first described their mar- The third graph [SSX (pdf 81)] is the computed emergent velous numerical reproduction of observed internetwork shocks profile on the same equivalent-temperature scale as the sec- in Fig. 2 (pdf 4) of Lites et al. (1993) using a data-driven RADYN ond graph to enable checking the Eddington-Barbier approx- simulation. Summed over space and time these shocks contribute imation. I did not add observed disk-center atlas profiles in (as “H2V and K2V cell grains”) to the asymmetry between the vi- view of the Ha–Ha Erlebnis in endnote 55: you might be olet and red Ca II H & K core reversals (endnote 11). The charac- tempted to overclaim a good fit as observation-proven model teristic intricate shock development pattern had been established truth. [Back][ Main call] in Fig. 4 (pdf 9) of Cram and Damé (1983) and was sketched in Fig. 2 (pdf 7) of Rutten and Uitenbroek (1991). The resulting 46. My IDL programs excellent match of the observed and simulated Ca II H2V grain Under my Recipes for IDL with this inventory. RH plot pro- sequences in Fig. 14 (pdf 23) is highly convincing (and adorns grams in rhlib, LTE line formation programs in ltelib, SDO stuff the proceedings cover). These well-known results were then for- in sdolib, SST data-handling programs in sstlib, image sequence mally published in Carlsson and Stein (1995) and Carlsson and browser showex in imagelib, etc. Installation is described in this Stein (1997). More in endnote 58.[Back][ Main call] SDO manual. I even offer a beginner IDL manual. For my SDO pipeline (endnote 64) do my alignment practical.[Back][ Main 50. Carlsson four-panel breakdown call] Figs 8 – 12 (pdf 16) of Carlsson and Stein (1994) also introduced the informative Carlsson four-panel spectral-feature breakdown 47. RH program diagrams. For Ca II H one is shown in [SSX (pdf 156)], four in The RH spectrum synthesis code of Uitenbroek (2001) is named Figs. 4 (pdf 5)– 7 (pdf 8) of Carlsson and Stein (1997), three for after Rybicki and Hummer (1992) and follows their scheme Na I D1 and one for Ca II 8542 Å in Fig. 9 (pdf 8) of Leenaarts of multi-level approximate lambda iteration [RTSA Sect. 5.3.2 et al. (2010). I recommend studying these – by understanding (pdf 145)] not iterating Λ for Jν but the Ψ operator for jν. It per- them you join a select club! mits overlapping lines, includes PRD and full-Stokes options, In my courses they became examination material. An astute exists in 1D, 2D, 3D, spherical, and Cartesian versions, and more student noted that the second panel of [SSX (pdf 156)] exhibits recently also in parallel multi-column “1.5D” (Pereira and Uiten- the wrong assumption of CRD. Carlsson and Stein argued that broek 2015). It has various options to treat the worrisome ultra- their neglect of Ca II H PRD was compensated by not including violet line haze (endnote 39). Mg II h & k. Compendium solar spectrum formation 29 Notes

These diagrams clearly demonstrate the importance of NLT (endnote 49) appear as fast-moving erratically-distributed com- effects even though the spectral synthesis assumed SE. The ob- ing and going brightness wisps in internetwork regions. served bright Ca II H2V grains (endnote 11) which the simulation Both the photospheric and chromosphere-pistoning patterns reproduced so admirably exist only because previous shocks de- appear random and chaotic. White and Cha (1973) studied ap- fine the opacity distribution along the line of sight. The bright parent phase patterns of the 5-minute oscillation and concluded H2V grain represents an upward traveling acoustic wave at too that they are totally random: “any description of the observed low height to show shock heating at its wavelength if the higher motion in terms of simple deterministic functions will be inad- atmosphere were stationary; it gains its bright visibility from equate”. Pure chaos – so much for future helioseismology! I large higher-up downdraft after the earlier passage of a wave was the referee, found it convincing, and gave a Utrecht col- that shocked higher up. The downdraft pulls the line core red- loquium reviewing 5-minute oscillations but only the discovery ward and uncovers the grain. [Back][ Main call] part (showing 35-mm slides of the ingenious pre-IDL 10K×10K image processing plates of Leighton et al. 1962) and the confu- 51. Acoustic and gravity wave heating sion part (speculating about gravity waves turning into Alfvénic Fossum and Carlsson (2005a) also laid the long-debated attri- waves after Lighthill 1967) itemized in this later display (pdf 1) bution of coronal heating to acoustic waves (e.g., Ulmschnei- but yet without enlightenment because I was not aware of the 82 der and Kalkofen 2003) to rest. An earlier candidate for coro- marvelous prediction of Ulrich (1970) . nal heating were gravity waves proposed by Whitaker (1963)79 At the time the same held for F.-L. Deubner who was then who introduced the “diagnostic” k − ω diagram (pdf 22) in so- contesting Frazier’s finding of substructure in the (k, ω) dia- 83 lar physics following geophysics Hines (1960)( endnote 53). The gram (Deubner 1972) . However, Deubner was asked to ref- k − ω diagrams in Fig. 3 (pdf 3) of Rutten and Krijger (2003) eree Wolff (1973) and so became aware of Ulrich (1970) and exhibit gravity-wave power and phase difference signatures. I also Wolff (1972) stating that the Sun is spherical and that Ul- doubt that gravity waves heat the corona (but see Straus et al. rich’s modes must be spherical harmonics. Deubner then recog- 2008); however, since gravity waves couple relatively easy into nized the parabolic mode-location curves in Ulrich’s predictive Alfvénic waves (Lighthill 1967) they may play a role in spicule- (k, ω) diagram in Fig. 2 (pdf 6) as the fish to catch. He knew from II heating of the quiet chromosphere (endnote 62). [Back][ Main Fig. 6 (pdf 17) of the authoritative (in beautiful French) (k, ω) 84 call] description of Mein (1966) that long duration and large ex- tent were the way to (k, ω) resolution and fast cadence and fine 52. MC shifts 1600–1700 Å sampling to (k, ω) extent. He found the ridges, best in his hand- drawn note-in-press diagram (pdf 4) in Deubner (1975) adding Away from disk-center magnetic bright points in AIA 1600 Å the detailed spherical-harmonic predictions of Ando and Osaki images shift limbward from their 1700 Å counterparts. These (1975). He showed it first at a conference in Nice, apologizing shifts do not result from higher formation but from increasing to Frazier (Deubner and Mc Intyre 1976)85, where it became the hole transparency. They are shown, quantified, undone and ex- proceedings cover (Cayrel and Steinberg 1976) and made him plained in Figs. 70–73 (pdf 88) of Rutten (2020). My SDO align- famous.86 Upshot: the photospheric p-mode pattern pistoning ment pipeline (endnote 64) has an option to unshift and subtract chromospheric shocks appears random but is a most regular and them to bring out flaring active-region fibrils (Vissers et al. 2015, beautiful harmonic chord (pdf 50).87 [Back][ Main call] Rutten 2016) that are brightened independently in 1600 Å im- ages by emission in C IV lines. [Back][ Main call] 54. Lockyer chromosphere, He I D3, flash color I typed Lockyer (1868) into ADS, remain its main citer so 53. Solar p-mode history far, and look forward to write “less refrangible” in publications “Clapotis” on sea charts marks locations with “wild” waves. In ([SSX (pdf 12)] ff). Lockyer’s low total of ADS citations (but regular wind-driven ocean waves material motion is mostly ver- steepening reads) shows how citation ranking utterly fails for the tical; the same holds for the 5-minute oscillation of the solar sur- first professor in astrophysics and the founding editor (during 50 face. Sea clapotis results from pattern interference by wave re- years!) of Nature. I recommend his delightful “The Spectroscope flection off capes, harbor quays, or in meeting of different swell and Its Applications”. systems growing into “wild” nonlinearity. The chromospheric In this report he discovered He I D3 (as did “Indian observer” 3-min shocks are vertically NLT interfering in complex fashion, Janssen shortly before him). It contributed its splendid yellow wild enough to need Carlsson breakdown diagrams (endnote 50). color (familiar to us from sodium street lights) to the chromatic On the other hand, the photospheric 5-minute oscillations that beauty of the emission lines he knew as prominence lines but saw drive the shocks are also complex interference patterns but out- standingly regular, linear, and precisely predictable. I summarize 82 The acknowledgment thanks adviser R.F. Christy for suggesting that this solar success story80 here because both types of omnipresent the self-excited oscillations might be more than numerical instability of oscillation affect solar lines (more in my decades-old 5-minute Christy’s stellar-structure code (the Nobel for starting helioseismology lecture (pdf 1)). should go to that). E.N. Frazier is acknowledged for suggesting that the In HMI Dopplergram movies the 5-minute oscillation discov- observed 5-minute oscillation might be the surface manifestation. 83 ered by Leighton et al. (1962) appears everywhere as patches of The telltale markings are by R.W. Noyes whose Solar Physics vol- 81 umes were scanned by ADS. apparently coherent oscillation somewhat larger than granules . 84 The patches are the surface interference pattern of many solar p- Deubner and Mein both completed conservatory schooling. I made them perform the Debussy cello sonata at the Zwaan retirement work- modes and contribute “radial macroturbulence” to the wash-out shop (Rutten and Schrijver 1994). of narrower lines than Mn I ones (endnote 12). In AIA 1700 Å 85 McIntyre was not involved but wrote an ADS-misleading superfluous movies the 3-min components pistoning Carlsson-Stein shocks comment on terrestrial seismology. Bahcall (1999) did better. 86 Scooping Ulrich and Rhodes who also had detected the ridges, while 79 Parker’s first PhD student and with his blessing (before nanoflares). Ulrich (1970) may have scooped Leibacher and Stein (1971). 80 Without the low-l part, see footnote 71. 87 I take my flute to public lectures “The Singing Sun” to demonstrate 81 A DKIST show-off granulation movie will need (k, ω) filtering as in overtone harmonics. P.Sütterlin’s conefilt.pro for the DOT. Compendium solar spectrum formation 30 Notes everywhere just outside the limb with his new (but “incomplete” wavelengths through NLT H I-top population boosting (end- and delayed beyond Janssen) spectroscope [SSX (pdf 12)], con- note 77). [Back][ Main call] cluding that “prominences are merely local aggregations of a gaseous medium which entirely envelopes the sun”. 57. Dynamic fibrils This was well before he saw himself the beautiful purple-pink Field-guided shocked waves in slanted active-region MCs that of the flash chromosphere at eclipses, but he knew the visually- (in my opinion) became the first success story in the century-old integrated color from prominences and chose the name to distin- struggle to understand the rich zoo of solar Hα features (e.g., guish it from “the white light-giving photosphere” and the “cool Hansteen et al. 2006, De Pontieu et al. 2007, Heggland et al. absorbing atmosphere” – a reversing-layer Fraunhofer-line in- 2007 [SSX (pdf 22)]). Their slant contributes cut-off frequency terpretation prior to Schuster-Schwarzschild modeling treated in lowering (Eq. 4 (pdf 6) of Michalitsanos 1973; Fig. 1 (pdf 2) of the SSA 3.3 (pdf 31) practical, the opposite to Milne-Eddington Bel and Leroy 1977). Their signature is repetitive parabolic ex- modeling. tension and (non-ballistic) retraction along their length in x − t The chromospheric flash color at eclipses during a few sec- timeslices (Fig. 7 (pdf 7) of De Pontieu et al. 2007). Many are onds after second and before third contact (longer near the edge visible as dark up-and-down dancing stalks in phased rows in of the totality strip where Houtgast went on purpose) sums this magnificent 5-s cadence SST/SOUP Hα movie (courtesy Lockyer’s lines dominated by Hα, plus white electron scatter- L.H.M. Rouppe van der Voort)88. The rows map MC arrange- ing (Sect. 7 on non-opening pdf page 8 of Jejciˇ cˇ and Heinzel ments in active-region network and plage; the phasing along 2009). The beautiful color represents H I, He I and freed H and rows is set by the interference-pattern undulations of the com- He electrons: enjoy it yet more by appreciating it as the stuff the bined p-modes (endnote 53). Their tops may heat to high tem- universe is made of. [Back][ Main call] perature and show up as bright grains in the Si IV lines in IRIS 1400 Å slitjaw images (Skogsrud et al. 2016). Short ones are ob- 55. Ha–Ha scattering served in sunspots (Rouppe van der Voort and de la Cruz Ro- Avrett stars yield Hα disk-center profiles that correspond rea- dríguez 2013). [Back][ Main call] sonably well to the observed disk-center atlas profile. Naively this may be seen as vindication of the model “chromosphere”. 58. Internetwork shocks in Ca II H and Hα However, about the same profile is also produced in a NLTE RH When Mats Carlsson hunted in the data of Lites et al. (1993) computation of Hα from an 1D RE atmosphere having no chro- he found only few pixels with multi-cycle spectrum-versus-time mosphere whatsoever, only an RE photospheric temperature de- behavior matching his RADYN output for the same pixel and se- cline ([SSX (pdf 3)], Figs 7–8 (pdf 6) of Rutten and Uitenbroek lected the best for Fig. 16 (pdf 19) of Carlsson and Stein (1997), 2012). √ showing two less-well matching samplings in Fig. 17 (pdf 20). The reason is that Hα makes its own ε scattering decline, The first two are columns labeled 110 and 132 in the observed more or less the same in Hα optical depth scale irrespective of time slices in Fig. 2 (pdf 4) of Rutten et al. (2008), the latter two geometrical height scale and presence or absence of an upper- columns 149 and 30. The latter shocks may have been too slanted photosphere opacity gap. I call this “Ha–Ha” line formation. It to produce vertical-column NLT signature as in the 1D simula- similarly implies that actual Hα from higher-located optically tion. and effectively thick fibrils samples similar scattering declines. The very quiet-target SST movie linked in endnote 48 and Furthermore,√ in the ALC7 atmosphere Ca II 8542 Å has here again shows fast-varying threads constituting Hα “mush- nearly the same ε decline with closely the same τ sampling room” shells following on and around Ca II H2V grains, indeed (comparison in [SSX (pdf 29)]). The solar-atlas profiles of Hα not radial along the line of sight: internetwork shocks clearly are and Ca II 8542 Å in Fig. 10 (pdf 10) of Rutten et al. (2011) differ 3D(t) rather than 1D or 2D phenomena even if the p-mode wave substantially, but mostly from the difference in atomic mass af- pistoning is primarily radial. Per-pixel Fourier analysis to estab- fecting thermal broadening and Hα’s abnormal opacity gap, both lish vertical phase relations is invalid for the resulting shocks and known to ALC7 so that it reproduces both profiles reasonably their mushroom aftermaths. The cadence of the movie is 7 sec- well – but so does a Kurucz RE atmosphere without chromo- onds but faster is clearly desirable, confirming the need for 1 s sphere and so do actual higher-up chromospheric fibrils (ha–ha). spectral Hα imaging with DKIST (endnote5 ). In quiet-Sun resolving observations the actual nonthermal The showex SST displays in my alignment practical show motions yield different scenes in the cores of these lines (Fig. 3 Hα mushroom threads well in the upper-left quiet corner of the (pdf 4) of Cauzzi et al. 2009) except for their time-averaged core SST field but at slower cadence. [Back][ Main call] widths which both sense temperature (Fig. 4 (pdf 5) and Fig. 6 (pdf 7) of Cauzzi et al. 2009). In the wings both lines show 59. Bifrost program spicules-II but Hα more and longer (Fig. 2 (pdf 4) of Rouppe This MHD simulation code of Gudiksen et al. (2011) is named van der Voort et al. 2009, Fig. 1 (pdf 3) of Sekse et al. 2012). after old-Norse Bifröst in the Icelandic Edda: burning rainbow You can see all this well with the last showex SST-data in- bridge between Earth and the realm of gods (slightly more invit- spection of my alignment practical.[Back][ Main call] ing than Pandora). Reference lists (until 2018 at the time of writ- ing) are given in [SSX (pdf 160)] ff. Many studies employ the 56. Rydberg H I lines public snapshot of Carlsson et al. (2016). The H I Rydberg ladder produces emission lines in the in- ADS abstract search "Bifrost" and "solar" gives publication frared by radiative-collisional recombination departure diffusion count and Hirsch N =85, h=19, solid enough for a solid perma- (Carlsson and Rutten 1992) similar to the Mg I 12-micron lines nent position. "MURaM" and "solar" N =96, h=20 idem for the in endnote 34. At that wavelength the H I and Mg I lines are MHD simulation program of Vögler et al. (2005). Maybe you still apart but they creep together for higher levels and longer should simulate. [Back][ Main call] wavelengths. The Rydberg levels producing lines near 12 micron 88 have larger population for Mg I than for H I (Fig. 1 (pdf 2)) of For Acrobat the 19 MB size may exceed its limit, needing man- Rutten and Carlsson 1994) but this likely reverses at ALMA ual download from https:robrutten.nl/rrweb/rjr-pubstuff/ lar-2/halpha-soup-half.mpg. Compendium solar spectrum formation 31 Notes

Since this contribution should operate similarly in any 60. Transition region versus heated chromosphere star with convective granulation and supergranulation, even I skip the so-called “transition region” because I don’t think such stars without any internal dynamo, it may contribute cycle- exists as a separable domain in the solar atmosphere. independently to building network and hence to cool-star basal The shocks in the Fig 1 movie of Leenaarts et al. (2007) lift flux.90 In solar quiet-Sun ultraviolet images (TRACE 1550, the corona (let’s define that as fully ionized hydrogen) above 1600, 1700 Å; SDO 1600, 1700 Å; IRIS 1400 Å slitjaws) net- them tremendously up and down, doing so effortlessly with work (also quiet network) usually shows larger brightness con- a steep space-time-varying temperature jump as instantaneous tribution (from MC hole deepening) than internetwork acoustics. highly-warped envelope, not as the global shell of standard mod- When network produced by the surface dynamo outweighs the els. The simulated corona there reaches lower in the MCs. This acoustic contribution by internetwork shocks it may dominate simulation describes quiet chromosphere; in active regions the the basal flux observed at these wavelengths (Table 2.7 on p. 68 outbursts are yet further from 1D shell description. of Schrijver and Zwaan 2008). Schrijver et al. (1989) found that The transition region in Avrett atmospheres is sampled by basal flux in Ca II H & K from solar-type stars equals the so- He II 304 Å but in closed-field regions the AIA images in this lar flux in supergranular cell centers dominated by internetwork line show quiet-Sun chromosphere. It resembles the fibrilar Hα shocks, but Schrijver (1992) found additional basal contribution chromosphere so well that one may co-align reversed-intensity in weak network and plage. With a cool-star survey Ayres (2021) quiet-Sun Hα and 304 Å images by pattern correlation, as shown also finds evidence that basal flux measures the presence of net- in blinkable Figs. 11 – 12 (pdf 29) and discussed in Appendix B work not made by an activity dynamo. [Back][ Main call] (pdf 11) of Rutten (2020). This pattern similarity extends even to AIA 131 Å images. 62. Spicules-II The 304 and 131 Å emitting gas around quiet network is In my (not so humble) opinion these are the main agents causing not a thin transitional envelope around Hα fibrils but repre- the fibrilar Hα and heated He II 304 Å (endnote 60) quiet-Sun sents hot ionization phases, i.e., hot tips of spicules-II and their chromospheres (while not heating the corona) being the only still hot aftermaths whereas Hα fibrils show cooling recom- candidate agent sufficiently ubiquitous even in quiet network. bining return gas NLT–darkest near the spicule-II ejecting net- I use the name generically for off-limb and on-disk manifesta- work (endnote 62). Both outline canopies with similar surface tions. patterns because each samples the quiet chromosphere around I noted them as “straws” near the limb in DOT Ca II H network in successive manifestations. The chromosphere ain’t movies(Fig. 1 (pdf 2)) of Rutten 2006). They were then de- stacked in layers but is dynamically structured and unstuck in tected and measured using Hinode off-limb Ca II H sequences, time.89 [Back][ Main call] attributed to Alvénic waves and named “spicules type-II” by de Pontieu et al. (2007) and De Pontieu et al. (2007)91. 61. Internetwork fields and basal flux Spicules-I reach less high above the limb, crowd into the Internetwork MCs riding supergranular outflows to cell bound- “spicule forest” in Hα limb images, and were attributed to aries as the “persistent flasher” of Brandt et al. (1994) are active-region dynamic fibrils (endnote 57) by Pereira et al. traditionally identified as disjoint members of newly emerged (2012). tiny bipolar ephemeral regions (Harvey and Martin 1973) on Spicules-II are ubiquitously ejected from quiet-Sun network their way to sustain quiet network (Schrijver et al. 1997, re- including monopolar network and also in coronal holes. On the view Sect. 5.1.2 of Schrijver and Zwaan 2008). Such quiet net- disk these ejections where likely described earlier as “fine mot- work tends to be small-scale bipolar whereas plage in active re- tles” (cf. de Jager 1957), as upflows by Title (1966) in his thesis gions and plage emerging in smaller failed active regions without (now on ADS) and as outer-wing Hα “jets” by Suematsu et al. sunspots tends to be unipolar obeying the large-scale dipole na- (1995) and in later Big Bear reports. More recently they became ture of active region emergence. Subsequent shredding into more known as “rapid blue-wing excursion” (RBE) (Rouppe van der open plage and active network then adds bipolarity. Voort et al. 2009) and “rapid red-wing excursion” (RRE) (Sekse More recently another type of internetwork field became at- et al. 2013), first in the wings of Ca II 8542 Å and then in the tributed to local field production by a granular-scale near-surface outer wings of Hα, with heating up to EUV visibility (Henriques turbulent dynamo. Its action was observationally suggested from et al. 2016). Hanle depolarization (Trujillo Bueno et al. 2004) and by the de- Spicule-II aftermaths of cooling recombining return flows tection of abundant primarily horizontal weak fields at granu- cause dark Hα-core fibrils around network (known as “course lar scales in sensitive full-Stokes spectropolarimetry with Hin- mottles” in spectroheliogram days) minutes later with their Hα ode (Lites et al. 2008, 2017), theoretically simulated by Pietarila opacities probably much enhanced by NLT cool-down retar- Graham et al. (2010). The primarily horizontal signature implies 90 that opposite-polarity closing loops occur on granular scales History: defining the minimal or “basal” level of chromospheric ac- within the photosphere, a lower and denser “magnetic carpet” tivity in solar-type stars was a major Zwaan-initiated activity at Utrecht than the traditional network-scale one made by ephemeral re- (e.g., Mewe et al. 1981, Schrijver et al. 1982, Zwaan 1986, Schrijver 1987b, Rutten et al. 1991). See Duncan et al. (1991) and reviews by gions and invoked to supply coronal heating (Simon et al. 2001, Schrijver (1987a, 1995) and in Sect. 2.7 of Schrijver and Zwaan (2008). Priest et al. 2002). No direct chromospheric interest therefore, 91 Author spelling: no bibtex a, b distinction here because the first spells but subsequent supergranular-flow migration giving frequent en- Dutch so that bibtex classifies a split personality (but wrongly: Dutch counters causing cancelation or enhancement up to fluxtube col- would be B. de Pontieu but De Pontieu and De Pontieu, B. as for De lapse led Gošic´ et al. (2014) to suggest that this superficial dy- Jager, C. for Kees de Jager – misspelled as de Jager, C. by ADS in 600+ namo action competes with ephemeral active region emergence bibitems, as many other Dutch). The proper spellings B. De Pontieu in quiet-network replenishment and L.H.M. Rouppe van der Voort, respectively, indicate where they originated. Presently ADS offers 4 solar publications by “van der Voort, 89 Billy Pilgrim syndrome. Multi-wavelength chromosphere observa- L.”, 102 for Kostyk, R.I. and 61 for Kostik, R.I. while half the responses tion is as Tralfamadorian multi-time viewing; Hα fibril canopies show for “Rutten, R.” are not mine. What’s in a name? Using my ADS author what happened before (Rutten et al. 2019). strings may help. With these I maintain my solar abstract collection. Compendium solar spectrum formation 32 Notes dance as in the post-shock clouds of Leenaarts et al. (2007). Ev- throughout the AIA database. They are the smallest cyan fea- idence are the darkest-darkest time-delay correlations in Rutten tures in the detector images in the ten triple blinkers starting at et al. (2019) and the downflowing RREs described by Bose et al. Fig 17 (pdf 35). Counting the many CBPs in these snapshots sug- (2021). More in endnote 67. gests that at any moment smaller CBPs are more numerous than The spicule-II driver mechanism is not known (Pereira et al. larger CBPs but the latter live much longer (see fifteen cutout 2012) but because spicules-II exhibit torsion modes near and at sequences starting at Fig. 48 (pdf 66)). the limb (De Pontieu et al. 2012, Rutten 2013) vorticity in gran- The cyan detector CBP feet are at chromospheric heights, ular convection affecting embedded network fluxtubes seems a come from above as in Fig. 12 (pdf 14) of Frogner et al. (2020), candidate, perhaps with fly-by reconnection as in Meyer et al. and mark (not cause) higher-up CBP loop reconnection caused (2012) kicking up Alvénic waves and perhaps with coupling to by bipolar network as evident in blinking with the triple magne- gravity waves (endnote 51). [Back][ Main call] tograms. It is a mistake to sum observed CBP brightnesses as con- 63. SDO quiet-chromosphere images tribution to coronal heating since this energy is radiated away I recommend page-blinking the three members of each set in (endnote 26); bright-flame statistics contain no information on Fig. 17 – 46 (pdf 35) of Rutten (2020). I find it uncanny how non-lost heating but flames do represent telltales on where it oc- well the detector images separate localized chromospheric heat- curs. The actual coronal heating patterns in the AIA 193 Å triple ing (ubiquitous grey patches around any network outside coro- images suggest negligible heating from the smallest (“campfire”) nal holes) from sparser localized coronal heating (cyan CBP feet CBPs and most from the largest. [Back][ Main call] in denser bipolar network also in coronal holes). This separa- tion implies that atmospheric heating is not traditionally “via 67. RBE–RRE–fibril ionization-recombination the chromosphere to the corona” (with larger energy need for RBEs and RREs reach higher degrees of hydrogen ionization the denser chromosphere) but that it is done by unrelated mech- than the shocks in the Fig 2 movie of Leenaarts et al. (2007). anisms for quiet chromosphere (Alfvénic-wave spicules-II in They extend further in Hα than in Ca II 8542 Å (lower-left panel my opinion) and quiet corona (reconnective CBPs in my opin- of Fig. 1 (pdf 3) of Sekse et al. 2012) because Ca II ionizes before ion). [Back][ Main call] H I and is much less abundant. Where their tips reach AIA EUV visibility (Henriques et al. 2016) hydrogen is nearly or fully ion- 64. SDO alignment programs ized. The resulting recombining fibrils are more numerous and A collection of IDL programs to request and download SDO also much longer in Hα than in Ca II 8542 Å (Fig. 3 (pdf 4) of sequences for any field, location, time, duration, typically Cauzzi et al. 2009). The time-averaged core-widths correlation for co-alignment with image sequences from another smaller- in the bottom row of Fig. 6 (pdf 7) of Cauzzi et al. (2009) sug- field higher-resolution telescope such as the SST. The AIA gests that both lines sample recombining-fibril temperatures but sequences get cross-aligned to HMI by projection-correcting with much larger fibril opacity in Hα, as in the return contrail in cross-correlation of large parallel cutouts around solar center, Fig. 5 (pdf 6) of Rutten and Rouppe van der Voort (2017). usually ten times better than by aia_prep.pro. The SDO se- You should check these findings yourself with the high- quences are then co-aligned with 0.1 arcsec precision to the other quality SST quiet-Sun data in my alignment practical. Its last telescope (or vice versa). See program manual and alignment showex command browses SST Hα and Ca II 8542 Å image se- practical. Alignment results are shown in Fig. 75 (pdf 93) of Rut- quences sampling many wavelengths and also adding the three ten (2020). core measures of Cauzzi et al. (2009) for both lines. Blink these For small fields this pipeline is fast by not handling full-disk while activating scatter diagrams and also hit the button for se- images but only JSOC im_patch cutouts. All programs use IDL quence time-averaging. For example, the well-correlated width- assoc to avoid memory loads. There is a versatile sequence width blink shows how much further the Hα spicules-II reach, browser showex.pro which can play, blink, manipulate, cross- and time-delay blinking Hα-0.6 Å against Hα+0.2 Å confirms correlate many movies in parallel with zoom-in magnification. Rutten et al. (2019) with yet better data showing striking darkest- The alignment practical contains examples browsing the darkest correlation between RBEs and subsequent redshifted SDO-aligned quiet-Sun SST data of endnote 48. The disk-center core fibrils around 5-min delay, and also that these reach close to coronal hole sampled there also sports a few “SolO campfires” the spicule-II firing network – I think because in the cooling re- (endnote 66). [Back][ Main call] √ turn Hα gains largest NLT overopacity there keeping it ε dark. 65. CBP foot visibility The minima and shift blinks show that in Ca II 8542 Å these re- turn fibrils get visible only close to network where also Ca II CBP feet are bright in AIA 304 Å and its “fire detector” (pdf 87) recombines from Ca III. Thus, large spicule-II H ionization is product with 131 Å but sharper in IRIS 1400 Å slitjaw images followed by cooling return recombination producing core fibrils because Si IV has similar temperature-defined presence as He II that are much more prominent in the Hα core. but 103 lower abundance. This is evident for the complex CBP The few-percent ionization in the shocks in the Fig 2 movie in the 1400 Å and 304 Å columns of Fig. 2 (pdf 3) of Rutten of Leenaarts et al. (2007) with its two-orders of magnitude elec- and Rouppe van der Voort (2017), discussed at the end of Ap- tron density increase already suffices to reach near-SB proton pendix B (pdf 14) of Rutten (2020). The SDO cutouts for a densities (3rd movie row). Hence large electron densities and SB dozen CBPs in Fig. 48 ff (pdf 66) there suggest additional 304 Å partitioning must also occur in spicules-II, warranting straight- brightness around bright CBP feet due to EUV irradiation from forward SB extinction comparisons as in Fig. 7 (pdf 7) of Rut- above. [Back][ Main call] ten and Rouppe van der Voort (2017) for the hot spicule-II phases while recognizing that these high extinction values per- 66. SolO campfires sist NLT-wise afterwards in the recombining return fibrils con- Press-release hullabalooed as new coronal heating agent trig- stituting opaque canopies while instantaneous SB values drop gering Rutten (2020) finding them already reported as smallest Boltzmann-steep with their temperature. CBPs by Falconer et al. (1998) and non-hullabalooed present Compendium solar spectrum formation 33 Notes

Similar SB curves imply the same for the mm wavelengths of Rutten et al. 2013). They are therefore similarly shielded and ALMA in Fig. 1 (pdf 4) of Rutten (2017b) where the H I free-free invisible for ALMA (prediction 9 (pdf 9) of Rutten 2017b con- extinction curves rise as steeply as for Hα in the Boltzmann in- firmed by da Silva Santos et al. 2020). crease and then saturate at increasingly higher values for longer EBs are not visible in the Na I D and Mg I b lines (Eller- wavelengths due to the Rayleigh-Jeans λ2 in free-free extinc- man 1917, Rutten et al. 2015) because they are hot enough to ff −3/2 2 tion αλ ∼Ne Np T λ [RTSA Eq. 2.79 (pdf 47)]. Hence, return fully ionize neutral metals. Similarly EBs appear bright in UV fibrils that gain sufficient non-E opacity in Hα to cause opaque continua from Balmer continuum brightening with neutral-atom canopies will be similarly opaque for ALMA. They appear dif- bound-free opacity ionized away (as in flare feet showing optical ferently because the fibril contrasts√ in ALMA canopies show dif- flares). Vissers et al. (2019) have defined and calibrated an EB ferent temperatures, not different ε darkening as for Hα, but detector using AIA 1700 Å images which opens the entire SDO apparent fibril widths are controlled by lateral scattering in Lyα database for studying EB occurrence. and similar. In AIA 1600 Å images EBs gain additional brightness from Summary: in heating gas the hydrogen-top populations and the C IV lines in the passband, similarly to the large Si IV en- extinction slide fast up along the SB curves to high values at hancements in Fig. 4 (pdf 7) and Fig. 8 (pdf 9) of Vissers et al. large hydrogen ionization, but in subsequent cooling they hang (2015). The latter figure also shows Mn I, Fe II and Ni II absorp- multiple minutes near these high values before finally dropping tion blends on high Si IV, C II and Mg II emission peaks. Fig- back to the SB value for the actual temperature if not already ure 2 (pdf 5) of Rutten (2016) added a schematic EB into the refreshed before then. The resulting temporary overpopulations last panel of Fig. 1 of the simulation of Leenaarts et al. (2007) 12 can be enormous since the b2 = 10 overpopulations in post- to show how these can arise from undisturbed EB surroundings shock cooling in the Fig 1 movie of Leenaarts et al. (2007) (not part of the phenomenon as suggested for IRIS bombs by are reached already after only few percent hydrogen ionization Peter et al. 2014). (Sect.7 ). [Back][ Main call] Because EBs are hotter and denser than internetwork and dynamic-fibril shocks in which Lyα already reaches collisional 68. ALMA–GONG–SDO alignment equilibrium their Hα extinction obeys SB partitioning, permit- Another confirmation that ALMA sees quiet-Sun chromospheric ting the straightforward SB opacity comparisons in Fig. 5 (pdf 9) canopies is that quiet-Sun ALMA images can be co-aligned with of Rutten (2016). In this figure Hα reaches similar extinction Hα images through cross-correlation. These are almost always maxima as the Mg II and Si IV IRIS lines, so that the actual EB available from GONG but with large quality variations with time Hα profile under the canopy is probably similar to Mg II k with and between stations, fortunately with ALMA’s longitude sam- a high peak and a small central dip as in Fig. 4 (pdf 7) of Vissers pling the usually better ones (Big Bear and Cerro Tololo). Since et al. (2015), hence with much larger profile-summed brightness the GONG Hα’s are full-disk they may contain a suited large than observed. It is a mistake to count the observed Hα EB wing quiet closed-field area somewhere for precise co-alignment (with enhancements as contribution to coronal heating since this en- Hα intensity reversed) to simultaneous AIA 304 Å images as in ergy is radiated away (endnote 26) while the Hα core photons the blink for the first SolO/EUI target in Figs. 11 – 12 (pdf 29) of are boxed in by overlying non-E-opaque canopies but appear in- Rutten (2020) and discussed in Appendix B (pdf 11) there. The sufficient to brighten these. GONG Hα part covering the ALMA target field may then serve Small hot photospheric structures as EBs containing intense (also reversedly) as ALMA-to-SDO intermediary for wherever Lyα irradiate their cool surroundings causing aureole boosting ALMA pointed. of the hydrogen top including ionization [SSX (pdf 171)] and I succeeded in this for quiet-Sun ALMA data with a good making the feature visible even if the actual reconnection site GONG sequence, found the ALMA-to-GONG matching eas- is too small. Non-E recombination is slow so that such boosts ier than GONG-to-304 Å matching, and indeed found good pat- linger longer, in a wake if the reconnection progresses upward. tern correlation between the ALMA, GONG and AIA chromo- The Zanstra mechanism for detour photon conversion in plan- spheres. [Back][ Main call] etary nebulae [SSX (pdf 144)] [IART Sect. 8.3.3.1 (pdf 137)] likely also operates in their cool surroundings. 69. Ellerman bombs and flaring active-region fibrils Presently EB interest has moved to non-active quiet-Sun QSEBs (“quiet-Sun Ellerman-like brightening”) marking ubiq- EBs are active-region MC reconnection events in the low photo- uitous but hard-to-detect small-scale photospheric reconnection sphere. After growing interest in the past decade92 and success- (Rouppe van der Voort et al. 2016, Danilovic 2017, Joshi et al. ful simulation by Hansteen et al. (2017) and Ortiz et al. (2020) 2020). they now represent another Hα feature that is largely understood. FAFs were called “transient loops” by Pariat et al. (2009). Here they are of chromospheric interest because the spectral Vissers et al. (2015) used “flaring arch filaments” but I renamed signature of these events is complex and instructive. EBs appear them to “flaring active-region fibril” in Rutten (2016) to avoid bright only in the Hα wings because their Hα line cores remain confusion with arch filament systems. In active regions with shielded by overlying dark active-region fibrils (Fig. 3 (pdf 4) of much flux emergence they appear in AIA 1600 Å image se- 92 I participated in Watanabe et al. (2011), Rutten et al. (2013), Vissers quences with shorter duration and more abrupt changes than et al. (2013), Rutten et al. (2015), Vissers et al. (2015), Rutten (2016), comparatively stable and roundish EBs, have elongated mor- Young et al. (2018), Vissers et al. (2019) and think that I helped boost phology, and show fast apparent brightness motion along fila- interest in Ellerman (1917): after half a century of disregard it shot up mentary strands. Further description including EB and FAF dis- in ADS citations reaching 100 in its centennial year. Moral: if your nice tinction in Appendix A (pdf 9) of Rutten (2020).[Back][ Main discovery goes unnoticed just relax waiting for ADS to revive it. But call] note more moral: this revival got my laptop confiscated by border po- lice spotting the b* word and taking months to classify me as harmless 70. Moss Ellerman b*er – therefore using b=burst in Young et al. (2018). But also note that a solar EB releasing 20 Gt (Rezaei and Beck 2015) is a bigger “Spongy” active-region plage showing many short loops and blast than the champion terrestrial 50 Mt “Tsar Bomba” hydrogen b*. bright footpoints interspersed with darker patches at small scales Compendium solar spectrum formation 34 Notes in 171 Å and other EUV lines (e.g., Berger et al. 1999b, 1999a, namic fibrils observed down-the-throat. There are long slender de Pontieu et al. 1999, Martens et al. 2000, De Pontieu et al. dark active-region filaments above polarity inversion lines (blink 2003). Polito et al. (2020) differentiate moss and outflow regions either against p2). The spot flare shows its concurrent loop and in active-region plage where the latter likely contribute to the feet brightening similarly in both line cores. slow wind and non-CME FIP generation (Fig. 8.8 (pdf 170) of Hα core I – minimum shift (p6 – p7): the second is the shift of Reames 2021; endnote 73). [Back][ Main call] the Hα profile minimum per pixel following Cauzzi et al. (2009), with blueshift dark. The hearts of the plage areas have small- 71. IBIS–SDO active-region blinker scale grainy appearance with sizable patches of large blueshift. Download here. Most of these images have enhanced grey scales The fibrils extending away show lower-amplitude fibril-to-fibril by using the square root or the square of the measured intensity modulation, a signature of their intermittent repetitivity in ex- plus clipping. All have a small area near center magnified four tending out and shrinking back (called “breathing” by a col- times in the upper-right corner for more detailed blinking. league). Blinking the inset shows short dark blueshifted heat- Here are my inventory and tentative interpretations (better ing onsets; the longer grey fibrils are return flows. The latter blink yourself! – how in footnote 42): show similar but finer structure in Ca II 8542 Å I (blink p5 – HMI Stokes I – V (p1 – p2): Stokes I shows only granules, p7); showex zoom-in shows two-branch scatter diagrams there pores and the spot, not the network and plage seen in Stokes V because Doppler-shifting the narrow Ca II 8542 Å core bright- (higher resolution would show abnormal granulation in plage). ens its intensity for both blue- and redshift. In Hα the umbra Plage brightens in such continuum sampling only near the limb, shows flashes on their way to running penumbral waves. Only from MC transparency in slanted viewing (righthand cartoon in the brightest are visible in Ca II 8542 Å but they are well-known Fig. 73 (pdf 91) of Rutten 2020). The flare is along connections from Ca II H & K as in Figs. 2 – 4 (pdf 4) of Rouppe van der from the spot to opposite polarity in the sunspot moat (blink p2 – Voort et al. (2003). The spot flare shows a bright core-shift fea- p6). ture in the loop and also a dark shift companion feature further HMI Stokes V – AIA 1700 Å (p2 – p3): ultraviolet MC bright- out, at image edge. Diametrically opposite on the other side of ening is not from higher-up heating (as commonly thought) but the spot there is a conspicuous bright-dark feature pair that also from magnetic hole deepening (lefthand cartoon and caption seems connected. of Fig. 73 (pdf 91) of Rutten 2020). Blinking these disk-center Hα core width (p8): width of the Hα core per pixel following images shows co-spatiality but towards the limb 1700 Å bright Cauzzi et al. (2009). The bottom panels of their Fig. 6 (pdf 7) points shift relatively limbward, 1600 Å yet more, also not due with quantification in their Fig. 9 (pdf 9) showed that this is a to offset in formation height but to increasing MC transparency proxy for temperature thanks to the small H atom mass mak- (endnote 52). A few bright 1700 Å points in internetwork mark ing thermal broadening dominate over other broadening. This is isolated MCs (footnote 61). Elsewhere the grey fine structure the most telltale scene in this blinker with respect to heating. in 1700 Å internetwork areas represents clapotispheric shocks The plage hearts are grainy and mostly hot. The bushes extend- (endnote 49). ing away from there are shorter than the fibrils in core intensity AIA 1700 Å – Ca II 8542 Å wing I (p3 – p4): in the (p6) because these are RBE-type heating onsets. Network RBEs reach much further out while heating and accelerating, but their Ca II 8542 Å wing MC brightening in plage and network is profile signature moves out from the core to the far wing (Fig. 5 hole brightening from less collisional damping and additional (pdf 6) of Rouppe van der Voort et al. 2009). The same likely downdrafts, but the inset suggests much shielding within plage. happens for these active-plage fibrils, as for the active-network Enhanced bright points in the sunspot moat are likely not fibrils in the data of Rutten and Rouppe van der Voort (2017) EBs since they lack corresponding excess 1700 Å brightness analyzed in Rutten et al. (2019). Blinking the inset against min- (Vissers et al. 2019) and may mark MC downdrafts common in imum shift (p7) confirms that the short dark blueshifted features moats (Beck et al. 2007). Internetwork areas show predominant there are part of these. The Dopplershift features of the spot flare reversed granulation (inset) and also sporadic shock brightening show heating here, also on the other side of the spot. in grains (Sect. 5.5 (pdf 13) of Rutten et al. 2011). Slender dark He I D3 – He II 304 Å (p9 – p10): the He I D3 equivalent- fibrils extending from plage and network represent the plage width image is naturally noisy but nevertheless intriguing. The counterpart of network RBEs combining blueshift and line two scenes show fair bright-bright correspondence, better fur- broadening in this blue-wing sampling as in Fig. 1 (pdf 3) of ther away from the spot and not 1:1. The proposed mechanism Rouppe van der Voort et al. (2009). that He I D3 obtains opacity by He I photoionization from irradi- Ca II 8542 Å wing I – core I (p4 – p5): dramatic difference. ation by hotter overlying gas, often invoked to explain coronal- The core shows the chromosphere with bright-grained network hole visibility in optical He I lines (particularly in more opaque and plage hearts signifying onset of heating. Magnified blinking He I 10830 Å), is therefore not the whole story. Furthermore, this shows that many of these dense bright grains are roughly cospa- blink shows significant feature shifts, also in the inset, for the tial with wing MCs, suggesting upward heating, mostly in unipo- bright beginnings of dark fibrils; they occur further out from lar plage (p2). Blinking to 171 Å (p11) shows a few areas appear- network in 304 Å whereas the He I D -opaquest (bright) parts ing as EUV moss (endnote 70) over bipolar plage (one patch is 3 remain co-located with corresponding Hα features. These off- above the cutout square but HMI does not show all opposite- sets suggest outwardly increasing heating as occurs for network polarity field). Dark fibrils extend far over internetwork but do RBEs in the bottom panels of fig. 6 (pdf 8) of Sekse et al. (2012), not cover all of that, some internetwork still shows up in the qui- here for fatter fibrils originating from plage. This blink also in- etest areas with shock brightening in grains of which some are also visible in the wing image. dicates that the 304 Å dark filaments are also present in He I D3, barely outlined by showing slightly less noise. Ca II 8542 Å core I – Hα core I (p5 – p6): similar scenes but all fibrils are much fatter and longer in Hα which I attribute 304 Å (p10): blinking against Hα core I (p6), Hα core to large non-E opacity (endnote 67). All internetwork is cov- Doppler (p7), Hα core width (p8) and He I D3 equivalent width ered by fibrils in Hα. The bright plage grains often correspond (p9) are all interesting. The first (p10 – p6) shows coarse sim- well between the two (inset). Some may be heated tips of dy- ilarity between Hα and 304 Å: active chromosphere roughly Compendium solar spectrum formation 35 Notes corresponds between these diagnostics, opposite in sign to the tween more horizontal and more vertical. On-disk RBEs and reversed-Hα – 304 Å correlation for quiet closed-field chromo- RREs are easiest detected in the Hα wings when they reach far sphere (endnote 60). The correlation improves with spatial and out from their launching network, therefore preferentially show- temporal averaging (options in showex). I have used this like- ing wide-spread around heating, but at the limb spicules-II are ness for direct co-alignment of Hα to AIA 304 Å for SST targets easiest noted and isolated amid their projective confusion when with much activity and many active-region filaments and it be- they are more upward, as shown in the inclination histograms in came an option in my alignment pipeline (endnote 64), but usu- Fig. 4 (pdf 5) of Pereira et al. (2012) and seen near the limb e.g., ally alignment of quiet-area AIA 304 Å to HMI magnetograms in this DOT straw movie and in Fig. 5 (pdf 5) of De Pontieu et al. plus SST alignment to HMI per granulation works better. The (2012). Hence distinction between around versus upward mor- phology and heating seems gradual, suggesting similar driving second (p10 – p7) blink shows that the brightest 304 Å patches but difference in appearance and detectability set by variation in and correspondingly the brightest He I D opacity patches match 3 dominating field topography defined by MC clustering density best with the plage hearts where Hα Doppler is fine-scale black- and polarity. and-white speckled indicating upward heating (endnote 72). The recent CLASP2 – IRIS co-observation of Ishikawa et al. 171 Å (p11): obviously mostly corona, no chromosphere. (2021) indicates that fields keep longitudinal strength up to Blinking to HMI Stokes I (p1) shows that all long loops connect 300 Gauss high in the chromospheric above dense plage, sug- spot areas. The mossy patches near image center correspond to gesting preponderance of straight-up fields. The precise pattern bright areas in both He diagnostics (p9 and p10); in Hα (p6) correlation between field strength and Mg II k core intensity and these may harbor dynamic fibrils but their identification needs the good correspondence with AIA He II 304 Å intensity in their a time sequence. Otherwise, the grainy areas of chromospheric Fig. S1 (pdf 44) confirm straight-up heating. plage heating in Hα width (p8) show no coronal counterpart. For He I D the sparse and minor on-disk presence in the The separation between chromospheric and coronal heating in 3 IBIS blinker (pdf 9) versus its bright Balmer-like prominence quiet closed-field areas is already striking (endnote 63); here it in the flash spectrum, also away from active regions, suggests extends to an active region. Some thick 304 Å filament parts are sensitive response to upward network and plage heating (end- also dark in 171 Å, suggesting bound-free scattering as in end- note 74). [Back][ Main call] note 14.[Back][ Main call] 73. FIP effect 72. Around and upward heating Significant offset between fast-and slow-wind composition In quiet areas dark fibrilar Hα canopies and similar-pattern but for the photospheric electron-donor elements Fe, Si and Mg bright He II 304 Å canopies are made by spicules-II extending (neutral-atom ionization “potential” below 10 eV) which are pre- from network (endnote 62). The active region in the IBIS blinker dominantly ionized where H and hence the gas are near-neutral shows similar but denser, fatter, longerH α fibril canopies (pdf 6) (e.g., Meyer 1985, 1991, Fig. 7 (pdf 11) of von Steiger and Geiss that are relatively dark in 304 Å (pdf 10). It seems likely that such 1989, Fig. 8.5 (pdf 164) of Reames 2021). It is generally at- active canopies stem from similar mechanisms as quiet canopies; tributed to ion-neutral fractionation involving the degree of line- I call these “around heating”. However, the hearts of plage and tying to open field and hence attributed to low-atmosphere ori- network show dense grainy heating inH α core width (pdf 8), gins (e.g., von Steiger and Geiss 1989, cf. Laming et al. 2019, especially where mostly unipolar, that I call “upward heating” – Reames 2021). meaning at chromospheric heights, not high-reaching as coronal The IRIS study of Polito et al. (2020) points to non-mossy plumes in coronal holes as in Wang et al. 1997. This heating upward active-region heating (endnote 72) in slow-wind gener- follows fields that do bend back to the surface but connect along ation. The match in the IBIS blinker of “brightest” plage-heart long coronal loops to faraway activity. patches in He I D3 opacity with grainiest Hα core Dopplershift Hα core-width (pdf 8) shows remarkable difference: grainy patches in plage (p7 versus p9 of the blinker) marking upward bright upward heating in the plage hearts versus around heating heating suggests to test the He I lines (endnote 74) as poten- as extended bushes pointing away from network and plage that tial proxies for locating upward-heating FIP-offset source ar- seem comparable to spicule-II onsets (endnote 71). I attribute eas. [Back][ Main call] this difference to MC clustering and upward field forcing. Ac- tive network and plage contain discrete MCs in larger number 74. Optical He I lines and at larger density than quiet network. In classic magnetostatic The formation of optical He I lines remains an outstanding fluxtube cartoons [SSX (pdf 44)] the field flares out with height; enigma. The main ones are He I 10830 Å, weaker He I D3 at in classic “Zürich wineglass” models (Fig. 4 (pdf 6) of Bünte 5876 Å, and He I 6678 Å “occasionally making its appearance” et al. 1993) neighboring fluxtubes inhibit flaring wide. Denser near-unipolar collections so have more upward field, as in the (Lockyer 1868). He I 10830 Å reaches largest extinction because simulation-assisted inversion map in Fig. 7 (pdf 9) of Jafarzadeh its lower level is metastable and collects population. Its observed et al. (2017). depth is claimed as coronal activity mapper in the optical (e.g., I think that the resulting dominance of upward heating over Fig. 1 (pdf 3) of Penn et al. 1998) – but blinking He I D3 against AIA EUV in the IBIS blinker shows better He I matching with around heating causes larger brightness in 304 Å, outshining the surrounding fibril canopies that then appear dark in byte-scaled chromospheric 304 Å than with coronal 171 Å. The line serves as comparison. Where plage and network are more bipolar they coronal hole mapper in the McIntosh Archive (e.g., Webb et al. have small-scale short-loop connectivity leading to moss appear- 2017), perhaps better than coronal activity mapping since the endnote 63 ance (endnote 70). AIA fire detectors ( ) show less heated quiet chromo- For both around and upward heating the driving mechanisms sphere in holes, hence less He II and He I-top population. are not known, as is the case for spicules-II producing canopies Lockyer proposed element “helios” for He I D3 which be- around quiet network in closed-field areas (fewer in coronal came helium and remains the only extra-terrestrially discovered holes, e.g., Fig. 42 (pdf 60) of Rutten 2020). However, also for element since coronium ([SSX (pdf 8)]) and nebulium ([SSX quiet canopies there is observational inclination dichotomy be- (pdf 142)] ff) went away. I am not aware of a sound explana- Compendium solar spectrum formation 36 Notes tion why his He I D3 is so extraordinary strong off-limb, equal [SSX (pdf 165)], the n =2 level gains near-Boltzmann popula- l to Hα in his engraved Fig 42 (pdf 95) and ranking between Hγ tion, Lyα has Sν ≈ Bν(T) and Hα has SB line extinction. and Hβ in brightness and off-limb extent in pages 94–95 of Ta- Hα does not show such features in its core when they are ble 3A (pdf 94) of Dunn et al. (1968), much brighter and higher underneath the opaque Hα fibril canopies spreading from net- than its Na I D neighbors whereas it is nearly invisible on the work and then samples them only in its outer wings, bright for disk and hence not even mentioned on page 258 (pdf 292) of the EBs but Doppler-shifted thermally-broadened dark for spicule-II line table of Moore et al. (1966) from the Utrecht Atlas of Min- RBEs and RREs. Hα also skips internetwork shocks by its opac- naert et al. (1940) in which the Na I D lines are the darkest in the ity chasm in the upper photosphere (endnote 48, [SSX (pdf 91)]). optical and the Balmer lines are prominent wide dips (Fig. 10 Lyα does not have such gap so that shifting my narrow dream (pdf 18) of Rutten 2019). Similarly, He I 10830 Å is absent in the passband from core to wing means sampling the whole non- ALC7 spectrum [SSX (pdf 94)]. How can the Balmer lines in coronal atmosphere down to the low (dark) photosphere (first the top of the hydrogen atom and the He I lines in the top of the plot of [SSX (pdf 90)]). helium atom differ so enormously between transverse and top- In addition, Lyα is the worst PRD line in the spectrum (sec- down viewing? 93 ond plot of [SSX (pdf 90)]) but this is actually a boon because Off-limb (no continuum) large emissivity can go with small it gives each part of the line its own source function and re- extinction through large εB and/or large ηS d contributions in sponse, as in the [SSF (pdf 111)] cartoon, and hence specific Eq.3 . Because He I D3 shows Balmer-like extent above the limb surface-sampling signature in the dream-resolved emergent pro- while the Balmer extent is dominated by quiet-Sun spicules- file. At each wavelength the emergent intensity does not sam- II and He II 304 Å shows quiet-chromosphere canopies pattern- ple the LTE source function√ within hot and dense features but matching the spicule-II-made reversed-Hα ones outside coro- the outer monofrequent εν scattering decline towards the fea- nal holes (endnote 68) I suggest that in quiet-Sun around heat- ture surface.96 As in endnote 29 the emergent intensity is de- ing (endnote 72) the He I top shares in the spicule-II emissivity fined around τeff = 1, polarization and Doppler signature around phase giving large εB off-limb brightness, but less in the subse- τ = 1. Outer-wing passbands that form in the upper photo- quent recombination phase giving Hα sufficient non-E√ extinc- sphere should be well modelable with MURaM and Bifrost (end- tion to show opaque canopies on the disk with ε darkness. note 44). [Back][ Main call] Spicule-II RBEs and RREs and return fibrils in Hα occur ubiq- uitously also in coronal holes but these show less heated 304 Å 76. Lyα features chromosphere94. Maybe He I outside coronal holes remains con- The Lyα image in Fig. 2 (pdf 6) of Rutten (2017b) is from the siderably photoionized by blanket-like EUV active-corona irra- VAULT-II flight (Korendyke et al. 2001, Vourlidas et al. 2010, diation from wide-around above so that the He I and He II 304 Å Patsourakos et al. 2007) and is one of 17 images near disk center. lines map coronal holes by their absence while Hα shows on- A second sequence of four limb images was assembled by Koza disk canopies there also. et al. (2009) into a limb movie and a centerward movie. They In upward heating in active plage non-E retardance in cool- show dense bright grains in plage and active region of which a after-hot cycling may affect the He I lines yet more strongly than few change already during the 51-s movies. Koza et al. (2009) Hα because the He I top is at doubled excitation energy. The cor- suggested that the latter are Lyα counterparts to Hα dynamic fibrils (endnote 57). respondence between He I D3 absorption and He II 304 Å bright- ness in the IBIS – SDO blinker (endnote 71) may be due to The dense more stable grains seem tips of longer-lived up- large retarded population of the He I top (cf. Golding et al. ward heating jets that may also produce the numerous bright 2017). [Back][ Main call] grains in non-mossy hearts of network and plage in the IBIS Hα core-width scene (pdf 8) and contribute the upward chromo- 75. Lyα dream sphere heating of endnote 72. Away from the bright-grain plage areas in the Lyα panel of Let me dream further re imaging spectroscopy in Lyα (foot- Fig. 2 (pdf 6) the emission blobs are fuzzier and appear remark- note 65; paradise = 5 × 5 arcmin field, 0.1 arcsec resolution, 20 ably as the short bright bushes extending towards internetwork wavelengths all at 1 s cadence)95 and compare that with what in the IBIS display (pdf 8), suggesting that they sample similar we see (awake) with SST/CRISP in Hα, recapitulating results spicule-II-like roundabout heating onsets. above. The lines are both H I α and even share the n=2 level but The dense non-E Hα fibril canopies in the lefthand panel differ tremendously. of Fig. 2 (pdf 6) of Rutten (2017b) and in the IBIS Hα image Any feature visible in Hα has much higher opacity in Lyα (pdf 6) appear dark not from sampling temperature but from that will only show its outer surface. Within hot and dense fea- √ Hα’s self-made ε scattering decline ([SSX (pdf 91)], end- tures as internetwork shocks, dynamic fibrils, spicules-II, EBs note 55) with large non-E overopacity for cooling recombining and EB-like lower-atmosphere reconnection sites Lyα has fast hydrogen. These canopies must be yet more opaque in Lyα with collision-up rates [SSX (pdf 166)] while remaining boxed-in S l ≈ b B(T) superthermal internal source function boosting but even though being champion resonance scatterer (many random 2 √ similarly sampling the monochromatic ε Lyα declines towards steps but of small length [SSX (pdf 90)]). The high radiative 1–2 the canopy surface, as suggested by dark quieter areas in Fig. 2 up and down rates then nearly cancel to net radiative rate zero (pdf 6) of Rutten (2017b) where they are not out-radiated by 93 The 400+-page book by Athay and Thomas (1961) (wrong order at brighter features in profile-summed brightness. The same applies ADS) addressed this question without clear answer. It was grandly titled to filaments. [Back][ Main call]

“Physics of the solar chromosphere” earning a scathing review from De 96 Jager ([SSX (pdf 15)] ff) for not mentioning magnetism nor dynamism. Ca II H & K are richer in sampling chromospheric canopies than Hα, De Jager added that Thomas’ “photoelectric control” producing S > B as in Fig. 2 (pdf 5) of Jafarzadeh et al. (2017) and in Fig. 3 (pdf 4) of was the main item of interest but even that was misleading (endnote 24). Bose et al. (2019) partly thanks to such PRD source function splitting 94 High-resolution showcase: blink the AIA fire detector and SST Hα but also because thermal broadening is less wide for heavier atoms. core-width images in my alignment practical. 95 Yet rosier dreaming: plus imaging spectroscopy in He I 584 Å. Compendium solar spectrum formation 37 References

References 77. Rydberg H I candidate for ALMA Al, N., Bendlin, C., Hirzberger, J., Kneer, F., & Trujillo Bueno, J., “Dynamics There may exist chromospheric features with sufficient Rydberg- of an enhanced network region observed in Hα”, 2004A&A...418.1131A ladder (endnote 56) NLT recombination emissivity for visibility ADS in the H I 30-α line at 231.901 GHz. The predictions of Carls- Allende Prieto, C. & García López, R. J., “A catalogue of accurate wavelengths son and Rutten (1992) together with α-line observations in Fig. 6 in the optical spectrum of the Sun”, 1998A&AS..131..431A ADS Ando, H. & Osaki, Y., “Nonadiabatic nonradial oscillations: an application to (pdf 5) of Clark et al. (2000) suggested no visibility even at the the five-minute oscillation of the sun.”, 1975PASJ...27..581A ADS limb for H I-α lines above n≈25 but our predictions were for the Antolin, P. & Rouppe van der Voort, L., “Observing the Fine Structure MACKKL star assuming SE. The expected NLT overextinction of Loops through High-resolution Spectroscopic Observations of Coro- of the H I top may make 30-α both higher-formed and visible nal Rain with the CRISP Instrument at the Swedish Solar Telescope”, (Sect. 6 (pdf 13) of Rutten 2017a). 2012ApJ...745..152A ADS Antolin, P., Vissers, G., Pereira, T. M. D., Rouppe van der Voort, L., & Scul- If so detection enables exciting measurement of profile- lion, E., “The Multithermal and Multi-stranded Nature of Coronal Rain”, 2 resolving Zeeman splitting which scales ∼ λ while competing 2015ApJ...806...81A ADS Dopplerwidth scales ∼λ [RTSA Eq. 3.66 (pdf 79)]97. Their ratio Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P., “The Chemical Composi- produces significant profile signature already at 12 µm for 300 – tion of the Sun”, 2009ARA&A..47..481A ADS Athay, R. G. & Skumanich, A., “Broadening of H and K Emission Cores and 500 Gauss in Fig. 9 (pdf 12) of Rutten and Carlsson (1994); at the Wilson-Bappu Effect.”, 1967AJ.....72..784A ADS hundredfold wavelength one Gauss chromospheric sensitivity is Athay, R. G. & Thomas, R. N.: 1961, Physics of the solar chromosphere in sight. [Back][ Main call] 1961psc..book.....A ADS Athay, R. G., Low, B. C., & White, O. R., “The Solar Interior-Atmospheric 78. Filament blobs System”, 2008ASPC..383..315A ADS Auer, L. H. & Mihalas, D., “Non-Lte Model Atmospheres. I. Radiative Equilib- Filaments and prominences may be seen as “merely local aggre- rium Models with - Alpha”, 1969ApJ...156..157A ADS gations of a gaseous medium which entirely envelopes the sun” Auer, L. H. & Mihalas, D., “Non-Lte Model Atmospheres. II. Effects of Balmer (Lockyer 1868) but yet they present the toughest features to un- α”, 1969ApJ...156..681A ADS Avrett, E. H. & Loeser, R., “The PANDORA Atmosphere Program (Invited Re- derstand in the rich solar Hα zoo. How can cool chromospheric view)”, 1992ASPC...26..489A ADS gas show up endlessly in the hot corona in these weird, exceed- Avrett, E. H. & Loeser, R., “Models of the Solar Chromosphere and ingly long, slender, high-reaching features rich in fine structure Transition Region from SUMER and HRTS Observations: Formation of and dynamics? A nightmare for plane-parallel colleagues. Who the Extreme-Ultraviolet Spectrum of Hydrogen, Carbon, and Oxygen”, ordered these? 2008ApJS..175..229A ADS Avrett, E., Tian, H., Landi, E., Curdt, W., & Wülser, J. P., “Modeling the Chro- The SST/SOUP Hα movie of endnote 57 shows an active- mosphere of a Sunspot and the Quiet Sun”, 2015ApJ...811...87A ADS region filament where small bright bullet-like blobs run inter- Avrett, E. H., “Solutions of the Two-Level Line Transfer Problem with Complete mittently along its length (first panel of Fig. 1 (pdf 2) in the anal- Redistribution”, 1965SAOSR.174..101A ADS ysis by Lin et al. 2012). Since then I have noted many such fast Avrett, E. H., “Models of the Solar Outer Photosphere”, 1990IAUS..138....3A ADS disturbances running through filaments in AIA 304 Å movies, Ayres, T. R., “Chromospheric scaling laws, width-luminosity correlations, and suggesting frequent ionization-recombination cycling. I suspect the Wilson-Bappu effect.”, 1979ApJ...228..509A ADS that they are not disturbances of a mean state but that they actu- Ayres, T. R., “How deep can one see into the Sun?”, 1989SoPh..124...15A ally make and remake filaments continuously with NLT cooling- ADS Ayres, T. R., “In the Trenches of the Solar-Stellar Connection. III. The HST/COS gas opacity enhancement causing the extraordinary He I D3-like Ecliptic-poles Stellar Survey (EclipSS)”, 2021ApJ...910...71A ADS off-limb visibility of prominences in Hα and He II 304 Å and Bahcall, J., “Ulrich’s Explanation for the Solar Five Minute Oscillations”, that this cannot be modeled assuming static SE. Let me re- 1999ApJ...525C1199B ADS peat (Sect.7 ) that modest hydrogen ionization in simple inter- Barklem, P. S., Piskunov, N., & O’Mara, B. J., “VizieR Online Data Cat- alog: Broadening of metallic lines by H collisions (Barklem+ 2000)”, network shocks already produces post-shock hydrogen-top re- 2000yCat..41420467B ADS tardance with 12 dex Hα overopacity (last panel of Fig 1 movie Baschek, B., Holweger, H., & Traving, G., “Ein ALGOL-Programm fur die of Leenaarts et al. 2007), indeed in the “extraordinary” ballpark. quantitative Analyse von Sternspektren”, 1966AAHam...8...26B ADS That filaments may reform in place as happened after the ex- Beck, C., Bellot Rubio, L. R., Schlichenmaier, R., & Sütterlin, P., “Magnetic properties of G-band bright points in a sunspot moat”, emplary failed Rosetta eruption of Mason et al. (2021) also sug- 2007A&A...472..607B ADS gests dynamic refurbishment. [Back][ Main call] Beckers, J. M. & Milkey, R. W., “The Line Response Function of Stellar Atmospheres and the Effective Depth of Line Formation”, 1975SoPh...43..289B ADS Beckers, J. M., Bridges, C. A., & Gilliam, L. B.: 1976, High resolu- tion spectral atlas of the solar irradiance from 380 to 700 nanometers 1976hrsa.book.....B ADS Beckers, J. M.: 1964, “A Study of the Fine Structures in the Solar Chromo- sphere”„ Ph.D. thesis, Sacramento Peak Observatory, Air Force Cambridge Research Laboratories, Mass., USA 1964PhDT...... 83B ADS Beckers, J. M., “Solar Spicules (Invited Review Paper)”, 1968SoPh....3..367B ADS Bel, N. & Leroy, B., “Analytical Study of Magnetoacoustic Gravity Waves”, 1977A&A....55..239B ADS Bell, R. A., Eriksson, K., Gustafsson, B., & Nordlund, A., “A grid of model atmospheres for metal-deficient giant stars. II.”, 1976A&AS...23...37B ADS Belluzzi, L., Trujillo Bueno, J., & Landi Degl’Innocenti, E., “The Magnetic Sensitivity of the Ba II D1 and D2 Lines of the Fraunhofer Spectrum”, 2007ApJ...666..588B ADS Berger, T. E., de Pontieu, B., Fletcher, L., et al., “What is Moss?”, 1999SoPh..190..409B ADS √ Berger, T. E., De Pontieu, B., Schrijver, C. J., & Title, A. M., “High- 97 Erratum under [RTSA Eq. 3.66 (pdf 79)]: FWHM = 2 ln(2)∆λD ≈ resolution Imaging of the Solar Chromosphere/Corona Transition Region”, 1.665∆λD 1999ApJ...519L..97B ADS Compendium solar spectrum formation 38 References

Blamont, J. E. & Roddier, F., “Precise Observation of the Profile of the Fraun- Clark, T. A., Naylor, D. A., & Davis, G. R., “Detection and limb brightening of hofer Strontium Resonance Line. Evidence for the Gravitational Red Shift on the H bt I n=20-19 Rydberg line in the submillimetre spectrum of the Sun”, the Sun”, 1961PhRvL...7..437B ADS 2000A&A...357..757C ADS Böhm-Vitense, E.: 1989, Introduction to stellar astrophysics. Vol. 2. Stellar at- Cook, J. W., Cheng, C. C., Jacobs, V. L., & Antiochos, S. K., “Ef- mospheres., Vol. 2 1989isa2.book.....B ADS fect of Coronal Elemental Abundances on the Radiative Loss Function”, Bonnet, R. M., “C. de Jager, Europe, Space, and the Sun”, 1989ApJ...338.1176C ADS 1996SoPh..169..233B ADS Cram, L. E. & Damé, L., “High spatial and temporal resolution observations of Bose, S., Henriques, V. M. J., Joshi, J., & Rouppe van der Voort, L., “Character- the solar CA II H line”, 1983ApJ...272..355C ADS ization and formation of on-disk spicules in the Ca II K and Mg II k spectral Cram, L. E., Rutten, R. J., & Lites, B. W., “On the formation of Fe II lines”, 2019A&A...631L...5B ADS lines in stellar spectra. I - Solar spatial intensity variation of λ 3969.4”, Bose, S., Joshi, J., Henriques, V. M. J., & Rouppe van der Voort, L., “Spicules 1980ApJ...241..374C ADS and downflows in the solar chromosphere”, 2021A&A...647A.147B ADS Curdt, W., Brekke, P., Feldman, U., et al., “The SUMER spectral atlas of solar- Bostanci, Z. F. & Al Erdogan, N., “Cloud modeling of a quiet solar region in disk features”, 2001A&A...375..591C ADS Halpha .”, 2010MmSAI..81..769B ADS da Silva Santos, J. M., de la Cruz Rodríguez, J., White, S. M., et al., Brandt, P. N., Rutten, R. J., Shine, R. A., & Trujillo Bueno, J., “On photospheric “ALMA observations of transient heating in a solar active region”, flows and chromospheric corks”, 1994ASIC..433..251B ADS 2020A&A...643A..41D ADS Brault, J. & Noyes, R., “Solar emission lines near 12 microns”, Danilovic, S., “Simulating Ellerman bomb-like events”, 1983ApJ...269L..61B ADS 2017A&A...601A.122D ADS Bruls, J. H. M. J., Rutten, R. J., & Shchukina, N. G., “The formation of helioseis- de Boer, T., Busse, R., Kappes, A., Klasen, M., & Zeinstra, S., “Neu- mology lines. I. NLTE effects in alkali spectra.”, 1992A&A...265..237B trino constraints to scotogenic dark matter interacting in the Sun”, ADS 2021arXiv210505613D ADS Bruls, J. H. M. J.: 1992, Formation of diagnostic lines in the solar spectrum, de Jager, C., “The hydrogen spectrum of the sun”, 1952RAOU....1.....D PhD thesis, Utrecht University Bruls1992b ADS ADS Bünte, M., Solanki, S. K., & Steiner, O., “Centre-to-limb variation of the Stokes de Jager, C., “The interpretation of hydrogen spectroheliograms”, V asymmetry in solar magnetic flux tubes”, 1993A&A...268..736B ADS 1957BAN....13..133D ADS Burgess, A., “Delectronic Recombination and the Temperature of the Solar de Jager, C., “The temperature minimum in the upper solar photosphere”, Corona.”, 1964ApJ...139..776B ADS 1963BAN....17..209D ADS Canfield, R. C. & Cram, L. E., “Analysis of the solar magnesium I spec- de Pontieu, B., Berger, T. E., Schrijver, C. J., & Title, A. M., “Dynamics of trum. II. Sensitivity of lambda 2852 to partial redistribution effects.”, Transition Region ‘Moss’ at high time resolution”, 1999SoPh..190..419D 1977ApJ...216..654C ADS ADS Canfield, R. C., “Deviations from Local Thermodynamic Equilibrium in Weak De Pontieu, B., Tarbell, T., & Erdélyi, R., “Correlations on Arcsecond Scales Complex Spectra”, 1971A&A....10...54C ADS between Chromospheric and Transition Region Emission in Active Regions”, Canfield, R. C., “Formation of Solar Rare-Earth Lines inside and outside H and 2003ApJ...590..502D ADS K”, 1971A&A....10...64C ADS De Pontieu, B., Hansteen, V. H., Rouppe van der Voort, L., van Noort, M., & Cannon, C. J., “Frequency-Quadrature Perturbations in Radiative-Transfer The- Carlsson, M., “High-Resolution Observations and Modeling of Dynamic Fib- ory”, 1973ApJ...185..621C ADS rils”, 2007ApJ...655..624D ADS Carlsson, M. & Rutten, R. J., “Solar hydrogen lines in the infrared”, de Pontieu, B., McIntosh, S., Hansteen, V. H., et al., “A Tale of 1992A&A...259L..53C ADS Two Spicules: The Impact of Spicules on the Magnetic Chromosphere”, Carlsson, M. & Stein, R. F., “Radiation shock dynamics in the solar chromo- 2007PASJ...59S.655D ADS sphere - results of numerical simulations”, 1994chdy.conf...47C ADS De Pontieu, B., McIntosh, S. W., Carlsson, M., et al., “Chromospheric Alfvénic Carlsson, M. & Stein, R. F., “Does a Nonmagnetic Solar Chromosphere Exist?”, Waves Strong Enough to Power the Solar Wind”, 2007Sci...318.1574D 1995ApJ...440L..29C ADS ADS Carlsson, M. & Stein, R. F., “Formation of Solar Calcium H and K Bright De Pontieu, B., Carlsson, M., Rouppe van der Voort, L. H. M., et al., “Ubiquitous Grains”, 1997ApJ...481..500C ADS torsional motions in type-II spicules”, 2012ApJ...752L..12D ADS Carlsson, M. & Stein, R. F., “Dynamic Hydrogen Ionization”, De Wijn, A. G., De Pontieu, B., & Rutten, R. J., “Fourier Analysis of Active- 2007ApJ...654.1128D 2002ApJ...572..626C ADS Region Plage”, ADS Introduction to Spectropolarimetry Carlsson, M., Rutten, R. J., & Shchukina, N. G., “The formation of the MG I del Toro Iniesta, J. C.: 2007, 2007insp.book.....D emission features near 12 microns”, 1992A&A...253..567C ADS ADS Del Zanna, G. & Mason, H. E., “Solar UV and X-ray spectral diagnostics”, Carlsson, M., Rutten, R. J., Bruls, J. H. M. J., & Shchukina, N. G., “The non-LTE 2018LRSP...15....5D ADS formation of Li I lines in cool stars”, 1994A&A...288..860C ADS Delbouille, L., Roland, G., & Neven, L.: 1973, Atlas photometrique du spectre Carlsson, M., Stein, R. F., Nordlund, Å., & Scharmer, G. B., “Observa- solaire de [lambda] 3000 a [lambda] 10000 1973apds.book.....D ADS tional Manifestations of Solar Magnetoconvection: Center-to-Limb Varia- Dere, K. P., “Ionization rate coefficients for the elements hydrogen through zinc”, tion”, 2004ApJ...610L.137C ADS 2007A&A...466..771D ADS Carlsson, M., Hansteen, V. H., Gudiksen, B. V., Leenaarts, J., & De Pontieu, B., Deubner, F. L. & Mc Intyre, “Observations of low wavenumber non radial “A publicly available simulation of an enhanced network region of the Sun”, eigenmodes of the sun”, 1976pmas.conf..259D ADS 2016A&A...585A...4C ADS Deubner, F.-L., “Some Properties of Velocity Fields in the Solar Photosphere. IV: Carlsson, M., De Pontieu, B., & Hansteen, V. H., “New View of the Solar Chro- Long Periods, Five Minute Oscillations, and the Supergranulation at Lower mosphere 2019ARA&A..57..189C ”, ADS Layers”, 1972SoPh...22..263D ADS Carlsson, M., “A computer program for solving multi-level non- Deubner, F. L., “Observations of low wavenumber nonradial eigenmodes of the LTE radiative transferproblems in moving or static atmospheres.”, sun.”, 1975A&A....44..371D ADS 1986UppOR..33.....C ADS Doyle, J. G., Jevremovic,´ D., Short, C. I., et al., “Solar Mn I 5432/5395 Å line Castor, J. I.: 2007, Radiation Hydrodynamics 2007rahy.book.....C ADS formation explained”, 2001A&A...369L..13D ADS Cauzzi, G. & Reardon, K., “The IBIS Mosaic”, 2012IAUSS...6E.511C ADS Duncan, D. K., Vaughan, A. H., Wilson, O. C., et al., “CA II H Cauzzi, G., Reardon, K., Rutten, R. J., Tritschler, A., & Uitenbroek, H., “The and K Measurements Made at Mount Wilson Observatory, 1966–1983”, solar chromosphere at high resolution with IBIS. IV. Dual-line evidence of 1991ApJS...76..383D ADS heating in chromospheric network”, 2009A&A...503..577C ADS Dunn, R. B., Evans, J. W., Jefferies, J. T., et al., “The Chromospheric Spectrum Cayrel, R. & Steinberg, M. (Eds.): 1976, Physique des Mouvements dans at the 1962 Eclipse”, 1968ApJS...15..275D ADS les Atmospheres Stellaires, Proc. CNRS Coll. Nice, September 1975 Eddington, A. S., “The formation of absorption lines”, 1976pmas.conf.....C ADS 1929MNRAS..89..620E ADS Chae, J., “Spectral Inversion of the Hα Line for a Plasma Feature in the Upper Eklund, H., Wedemeyer, S., Snow, B., et al., “Characterization of shock Chromosphere of the Quiet Sun”, 2014ApJ...780..109C ADS wave signatures at millimetre wavelengths from Bifrost simulations”, Chandrasekhar, S. & Breen, F. H., “On the Continuous Absorption Coefficient 2021RSPTA.37900185E ADS of the Negative Hydrogen Ion. III.”, 1946ApJ...104..430C ADS Ellerman, F., “Solar Hydrogen “bombs””, 1917ApJ....46..298E ADS Chang, E. S. & Noyes, R. W., “Identification of the solar emission lines near 12 Falconer, D. A., Moore, R. L., Porter, J. G., & Hathaway, D. H., “Network microns”, 1983ApJ...275L..11C ADS Coronal Bright Points: Coronal Heating Concentrations Found in the Solar Cimino, M., Cacciani, A., & Sopranzi, N., “An Instrument to measure Solar Magnetic Network”, 1998ApJ...501..386F ADS Magnetic Fields by an Atomic-Beam Method”, 1968SoPh....3..618C Feautrier, P., “Sur la resolution numerique de l’equation de transfert.”, ADS 1964CR....258.3189F ADS Compendium solar spectrum formation 39 References

Fokker, A. D. & Rutten, R. J., “Coronal scattering of radiation from an Heinzel, P., Gunár, S., & Anzer, U., “Fast approximate radiative transfer method anisotropically radiating solar radio source.”, 1967BAN....19..254F for visualizing the fine structure of prominences in the hydrogen Hα line”, ADS 2015A&A...579A..16H ADS Fontenla, J. M., Avrett, E. H., & Loeser, R., “Energy Balance in the Solar Transi- Henriques, V. M. J., Kuridze, D., Mathioudakis, M., & Keenan, F. P., “Quiet-Sun tion Region. III. Helium Emission in Hydrostatic, Constant-Abundance Mod- Hα Transients and Corresponding Small-scale Transition Region and Coro- els with Diffusion”, 1993ApJ...406..319F ADS nal Heating”, 2016ApJ...820..124H ADS Fontenla, J. M., Curdt, W., Haberreiter, M., Harder, J., & Tian, H., Hines, C. O., “Internal atmospheric gravity waves at ionospheric heights”, “Semiempirical Models of the Solar Atmosphere. III. Set of Non-LTE 1960CaJPh..38.1441H ADS Models for Far-Ultraviolet/Extreme-Ultraviolet Irradiance Computation”, Holweger, H. & Müller, E. A., “The Photospheric Barium Spectrum: So- 2009ApJ...707..482F ADS lar Abundance and Collision Broadening of Ba II Lines by Hydrogen”, Fontenla, J. M., Stancil, P. C., & Landi, E., “Solar Spectral Irradiance, Solar 1974SoPh...39...19H ADS Activity, and the Near-Ultra-Violet”, 2015ApJ...809..157F ADS Holweger, H., “Ein empirisches Modell der Sonnenatmosphäre mit lokalem Fossum, A. & Carlsson, M., “High-frequency acoustic waves are not sufficient thermodynamischem Gleichgewicht”, 1967ZA.....65..365H ADS to heat the solar chromosphere”, 2005Natur.435..919F ADS House, L. L., “Ionization Equilibrium of the Elements from H to Fe.”, Fossum, A. & Carlsson, M., “Response Functions of the Ultraviolet Fil- 1964ApJS....8..307H ADS ters of TRACE and the Detectability of High-Frequency Acoustic Waves”, Houtgast, J.: 1942, “The variations in the profiles of strong Fraunhofer lines 2005ApJ...625..556F ADS along a radius of the solar disc”„ Ph.D. thesis,- 1942PhDT...... 12H Frogner, L., Gudiksen, B. V., & Bakke, H., “Accelerated particle beams in a 3D ADS simulation of the quiet Sun”, 2020A&A...643A..27F ADS Hubený, I. & Mihalas, D.: 2014, Theory of Stellar Atmospheres, Princeton Froment, C., Antolin, P., Henriques, V. M. J., Kohutova, P., & Rouppe van der University Press, Princeton 2014tsa..book.....H ADS Voort, L. H. M., “Multi-scale observations of thermal non-equilibrium cycles Hubený, I., “Probabilistic Interpretation of Radiative Transfer - Part Two - Ry- in coronal loops”, 2020A&A...633A..11F ADS bicki Equation”, 1987A&A...185..336H ADS Gary, D. E. & Hurford, G. J.: 2004, Radio Spectral Diagnostics, Vol. 314, 71 Hubený, I., “Probabilistic interpretation of radiative transfer. I - The square root 2004ASSL..314...71G ADS of epsilon law. II - Rybicki equation”, 1987A&A...185..332H ADS Gingerich, O. & de Jager, C., “The Bilderberg Model of the Photosphere and Hummer, D. G., “Non-coherent scattering: I. The redistribution function with Low Chromosphere”, 1968SoPh....3....5G ADS Doppler broadening”, 1962MNRAS.125...21H ADS Gingerich, O., Noyes, R. W., Kalkofen, W., & Cuny, Y., “The Harvard- Ishikawa, R., Bueno, J. T., del Pino Alemán, T., et al., “Mapping so- Smithsonian reference atmosphere”, 1971SoPh...18..347G ADS lar magnetic fields from the photosphere to the base of the corona”, Giovanelli, R. G., “Excitation of hydrogen and CaII under chromospheric con- 2021SciA....7.8406I ADS ditions”, 1967AuJPh..20...81G ADS Jafarzadeh, S., Rutten, R. J., Solanki, S. K., et al., “Slender Ca II Goldberg, L. & Noyes, R. W., “Origin of Emission Cores in Lines of Ionized H Fibrils Mapping Magnetic Fields in the Low Solar Chromosphere”, Calcium and Magnesium.”, 1964AJ.....69R.542G ADS 2017ApJS..229...11J ADS Goldberg, L., “The Origin of the Emission Reversals in the Fraunhofer h- and Jefferies, J. T. & Thomas, R. N., “Source Function in a Non-Equilibrium Atmo- K-Lines”, 1964ApJ...140..384G ADS sphere. III. The Influence of a Chromosphere.”, 1959ApJ...129..401J Golding, T. P., Leenaarts, J., & Carlsson, M., “Formation of the helium extreme- ADS UV resonance lines”, 2017A&A...597A.102G ADS Jefferies, J. T.: 1968, Spectral line formation, Blaisdell, Waltham Goldman, A., Blatherwick, R. D., Murcray, F. H., et al.: 1980, New atlas of IR so- 1968slf..book.....J ADS lar spectra, NASA STI/Recon Technical Report N 1980STIN...8031298G Jejciˇ c,ˇ S. & Heinzel, P., “Electron Densities in Quiescent Prominences Derived ADS from Eclipse Observations”, 2009SoPh..254...89J ADS Gomez, M. T., Severino, G., & Rutten, R. J., “Photospheric dynamics and the Jordan, C., “The ionization equilibrium of elements between carbon and nickel”, 1969MNRAS.142..501J NLTE formation of the solar K I 769.9 NM line”, 1991A&A...244..501G ADS ADS Jordan, S. D., “Obituary: Richard Nelson Thomas, 1921-1996”, 1996BAAS...28.1465J Gošic,´ M., Bellot Rubio, L. R., Orozco Suárez, D., Katsukawa, Y., & del Toro ADS Joshi, J., Rouppe van der Voort, L. H. M., & de la Cruz Rodríguez, J., “Sig- Iniesta, J. C., “The Solar Internetwork. I. Contribution to the Network Mag- natures of ubiquitous magnetic reconnection in the lower solar atmosphere”, netic Flux”, 2014ApJ...797...49G ADS 2020A&A...641L...5J ADS Gray, D. F., “A test of the micro-macroturbulence model on the solar flux spec- Kato, Y., Steiner, O., Steffen, M., & Suematsu, Y., “Excitation of Slow trum.”, 1977ApJ...218..530G ADS Modes in Network Magnetic Elements Through Magnetic Pumping”, Gray, D. F.: 1992, The Observation and Analysis of Stellar Photospheres., 2nd 2011ApJ...730L..24K ADS edition, Cambridge University Press 1992oasp.book.....G ADS Keller, C. U. & von der Lühe, O., “Solar speckle polarimetry”, Gray, D. F.: 2005, The Observation and Analysis of Stellar Photospheres, 3rd 1992A&A...261..321K ADS edition, Cambridge University Press 2005oasp.book.....G ADS Korendyke, C. M., Vourlidas, A., Cook, J. W., et al., “High-resolution Imaging Greve, A. & Zwaan, C., “Methods for the analysis of stellar spectra veiled by of the Upper Solar Chromosphere: First Light Performance of the Very-high- lines (III).”, 1980A&A....90..239G ADS Resolution Advanced ULtraviolet Telescope”, 2001SoPh..200...63K Grotrian, W., “Ergebnisse der Potsdamer Expedition zur Beobachtung der ADS Sonnenfinsternis am 9. Mai 1929 in Takengon (Nordsumatra). 6. Mit- Kostik, R. I. & Orlova, T. V., “On the Possible Mechanism of Formation of teilung. Über die Intensitätsverteilung des kontinuierlichen Spektrums der Emission Rim in Hydrogen Filaments”, 1975SoPh...45..119K ADS inneren Korona. Mit 8 Abbildungen. (Eingegangen am 27. Juni 1931)”, Kostik, R. I., Shchukina, N. G., & Rutten, R. J., “The solar iron abundance: not 1931ZA...... 3..199G ADS the last word.”, 1996A&A...305..325K ADS Gudiksen, B. V., Carlsson, M., Hansteen, V. H., et al., “The stellar Kostik, R., Khomenko, E., & Shchukina, N., “Solar granulation from atmosphere simulation code Bifrost. Code description and validation”, photosphere to low chromosphere observed in Ba II 4554 Å line”, 2011A&A...531A.154G ADS 2009A&A...506.1405K ADS Gurtovenko, E., Ratnikova, V., & de Jager, C., “On the Average Optical Depth Kourganoff, V.: 1952, Basic methods in transfer problems; radiative equilibrium of Formation of Weak Fraunhofer Lines”, 1974SoPh...37...43G ADS and neutron diffusion, Clarendon Press, Oxford 1952bmtp.book.....K Hansteen, V. H., De Pontieu, B., Rouppe van der Voort, L., van Noort, M., & ADS Carlsson, M., “Dynamic Fibrils Are Driven by Magnetoacoustic Shocks”, Koza, J., Rutten, R. J., & Vourlidas, A., “Dynamic Lyα jets”, 2006ApJ...647L..73H ADS 2009A&A...499..917K ADS Hansteen, V. H., Archontis, V., Pereira, T. M. D., et al., “Bombs and Flares at Krijger, J. M., Rutten, R. J., Lites, B. W., et al., “Dynamics of the so- the Surface and Lower Atmosphere of the Sun”, 2017ApJ...839...22H lar chromosphere. III. Ultraviolet brightness oscillations from TRACE”, ADS 2001A&A...379.1052K ADS Harvey, K. L. & Martin, S. F., “Ephemeral Active Regions”, Kurucz, R. L., “Including All the Lines”, 2009AIPC.1171...43K ADS 1973SoPh...32..389H ADS Labs, D. & Neckel, H., “Remarks on the Convergency of Photospheric Model Heeck, J., “How Stable is the Photon?”, 2013PhRvL.111b1801H ADS Conceptions and the Solar Quasi Continuum”, 1972SoPh...22...64L Heggland, L., De Pontieu, B., & Hansteen, V. H., “Numerical Simulations of ADS Shock Wave-driven Chromospheric Jets”, 2007ApJ...666.1277H ADS Lambert, J., Paletou, F., Josselin, E., & Glorian, J.-M., “Numerical Heinzel, P., Mein, N., & Mein, P., “Cloud model with variable source function radiative transfer with state-of-the-art iterative methods made easy”, for solar Hα structures. II. Dynamical models”, 1999A&A...346..322H 2016EJPh...37a5603L ADS ADS Laming, J. M., Vourlidas, A., Korendyke, C., et al., “Element Abundances: A New Diagnostic for the Solar Wind”, 2019ApJ...879..124L ADS Compendium solar spectrum formation 40 References

Lamy, P. L., Nguyen-Trong, T., Adjabschirzadeh, A., & Koutchmy, Mein, P., “Inhomogeneities in the Solar Atmosphere from the Ca II Infra-Red S., “Astronomical applications of infrared television imaging.”, Lines”, 1971SoPh...20....3M ADS 1979A&A....77..257L ADS Mewe, R., Schrijver, C. J., & Zwaan, C., “Coronal Activity an F-Type G-Type Landi Degl’Innocenti, E. & Landolfi, M.: 2004, Polarization in Spectral Lines, and K-Type Stars”, 1981SSRv...30..191M ADS Vol. 307 2004ASSL..307.....L ADS Meyer, K. A., Mackay, D. H., & van Ballegooijen, A. A., “Solar Mag- Leenaarts, J. & Carlsson, M., “MULTI3D: A Domain-Decomposed 3D Radiative netic Carpet II: Coronal Interactions of Small-Scale Magnetic Fields”, Transfer Code”, 2009ASPC..415...87L ADS 2012SoPh..278..149M ADS Leenaarts, J., Rutten, R. J., Sütterlin, P., Carlsson, M., & Uitenbroek, H., “DOT Meyer, J. P., “Solar-stellar outer atmospheres and energetic particles, and tomography of the solar atmosphere. VI. Magnetic elements as bright points galactic cosmic rays”, 1985ApJS...57..173M ADS in the blue wing of Hα”, 2006A&A...449.1209L ADS Meyer, J.-P., “Diagnostic methods for coronal abundances”, Leenaarts, J., Carlsson, M., Hansteen, V., & Rutten, R. J., “Non- 1991AdSpR..11a.269M ADS equilibrium hydrogen ionization in 2D simulations of the solar atmosphere”, Michalitsanos, A. G., “The Five Minute Period Oscillation in Magnetically Ac- 2007A&A...473..625L ADS tive Regions”, 1973SoPh...30...47M ADS Leenaarts, J., Rutten, R. J., Reardon, K., Carlsson, M., & Hansteen, V., Mihalas, D.: 1970, Stellar atmospheres, W. H. Freeman and Co, San Francisco “The Quiet Solar Atmosphere Observed and Simulated in Na I D1”, 1970stat.book.....M ADS 2010ApJ...709.1362L ADS Mihalas, D.: 1978, Stellar atmospheres, 2nd edition, W. H. Freeman and Co, Leenaarts, J., Carlsson, M., & Rouppe van der Voort, L., “The Formation of the San Francisco 1978stat.book.....M ADS Hα Line in the Solar Chromosphere”, 2012ApJ...749..136L ADS Minnaert, M., Houtgast, J., & Mulders, G. F. W.: 1940, Photomet- Leenaarts, J., Pereira, T., & Uitenbroek, H., “Fast approximation ric atlas of the solar spectrum, Sterrewacht “Sonnenborgh”, Utrecht of angle-dependent partial redistribution in moving atmospheres”, 1940pass.book.....M ADS 2012A&A...543A.109L ADS Minnaert, M., “On the continuous spectrum of the corona and its polarisation. Leenaarts, J., “Radiation hydrodynamics in simulations of the solar atmo- With 3 figures. (Received July 30, 1930)”, 1930ZA...... 1..209M ADS sphere”, 2020LRSP...17....3L ADS Molnar, M. E., Reardon, K. P., Chai, Y., et al., “Solar Chromospheric Temper- Leibacher, J. W. & Stein, R. F., “A New Description of the Solar Five-Minute ature Diagnostics: A Joint ALMA-Hα Analysis”, 2019ApJ...881...99M Oscillation”, 1971ApL.....7..191L ADS ADS Leighton, R. B., Noyes, R. W., & Simon, G. W., “Velocity Fields in the Solar Molowny-Horas, R., Heinzel, P., Mein, P., & Mein, N., “A non-LTE inversion Atmosphere. I. Preliminary Report.”, 1962ApJ...135..474L ADS procedure for chromospheric cloud-like features”, 1999A&A...345..618M Lemke, M. & Holweger, H., “A non-LTE study of the solar emission lines near ADS 12 microns”, 1987A&A...173..375L ADS Moore, C. E. & Merrill, P. W.: 1968, Partial Grotrian Diagrams of Astrophysical Lighthill, M. J., “Predictions on the Velocity Field Coming from Acoustic Noise Interest 1968pgda.book.....M ADS and a Generalized Turbulence in a Layer Overlying a Convectively Unstable Moore, C. E., Minnaert, M. G. J., & Houtgast, J.: 1966, The solar spectrum 2935 Atmospheric Region”, 1967IAUS...28..429L ADS Å to 8770 Å, NBS 1966sst..book.....M ADS Lin, Y., Engvold, O., & Rouppe van der Voort, L. H. M., “Small- Moore, C. E.: 1959, A multiplet table of astrophysical interest. Part 1 scale, Dynamic Bright Blobs in Solar Filaments and Active Regions”, 1959mtai.book.....M ADS 2012ApJ...747..129L ADS Murcray, F. J., Goldman, A., Murcray, F. H., et al., “Observation of new emission Linsky, J. L. & Avrett, E. H., “The Solar H and K Lines”, lines in the infrared solar spectrum near 12.33, 12.22, and 7.38 microns”, 1970PASP...82..169L ADS 1981ApJ...247L..97M ADS Lites, B. W., Rutten, R. J., & Kalkofen, W., “Dynamics of the Solar Chro- Neckel, H. & Labs, D., “The solar radiation between 3300 and 12500 Å”, mosphere. I. Long-Period Network Oscillations”, 1993ApJ...414..345L 1984SoPh...90..205N ADS ADS Neckel, H., “Announcement”, 1999SoPh..184..421N ADS Lites, B. W., Kubo, M., Socas-Navarro, H., et al., “The Horizontal Magnetic Neckel, H., “On the sun’s absolute disk-center and mean disk intensi- Flux of the Quiet-Sun Internetwork as Observed with the Hinode Spectro- ties, its limb darkening, and its ‘limb temperature’ (λλ330 to 1099 nm)”, Polarimeter”, 2008ApJ...672.1237L ADS 2003SoPh..212..239N ADS Lites, B. W., Rempel, M., Borrero, J. M., & Danilovic, S., “Are Internet- Neckel, H., “Analytical Reference Functions F(λ) for the Sun’s Limb work Magnetic Fields in the Solar Photosphere Horizontal or Vertical?”, Darkening and Its Absolute Continuum Intensities (λλ 300 to 1100 m)”, 2017ApJ...835...14L ADS 2005SoPh..229...13N ADS Observation and Analysis of the Solar Neutral Iron Lites, B. W.: 1972, “ Ng, K. C., “Hypernetted chain solutions for the classical one-component plasma Spectrum. Ph.D. thesis ”„ , UNIVERSITY OF COLORADO AT BOULDER. up to Γ=7000”, 1974JChPh..61.2680N ADS 1972PhDT...... 7L ADS Nordlund, Å., Stein, R. F., & Asplund, M., “Solar Surface Convection”, Lites, B. (Ed.): 1985, Chromospheric Diagnostics and Modelling, Sacramento 2009LRSP....6....2N ADS Peak Workshop, Sunspot, August 1984 1985cdm..proc.....L ADS Nordlund, A., “Modelling of Small-Scale Dynamical Processes: Convection and Liu, S. Y. & Skumanich, A., “An empirical interpretation for the time evolution Wave Generation (Keynote)”, 1984ssdp.conf..181N ADS of the Ca ii K line”, 1974SoPh...38..109L ADS Nordlund, A., “A Re-evaluation of the Granular ∆ Irms”, Livingston, W., Wallace, L., White, O. R., & Giampapa, M. S., “Sun-as-a-Star 1984ssdp.conf..174N ADS Spectrum Variations 1974-2006”, 2007ApJ...657.1137L ADS Nordlund, Å., “The 3-D structure of the magnetic field, and its interaction with Lockyer, J. N., “Spectroscopic Observation of the Sun, No. II.”, granulation.”, 1985MPARp.212..101N ADS 1868RSPS...17..131L ADS , “The dynamics of granulation, and its interaction with the radia- Madjarska, M. S., Chae, J., Moreno-Insertis, F., et al., “The chromospheric Nordlund, Å. 1985MPARp.212....1N component of coronal bright points. Coronal and chromospheric responses tion field.”, ADS Introduction to stellar atmospheres and interiors to magnetic-flux emergence”, 2021A&A...646A.107M ADS Novotny, E.: 1973, , Oxford 1973itsa.book.....N Madjarska, M. S., “Coronal bright points”, 2019LRSP...16....2M ADS University Press ADS Magain, P., “Contribution functions and the depths of formation of spectral Noyes, R. W., “Dielectronic Recombination and the Solar H and K Lines”, lines”, 1986A&A...163..135M ADS 1965SAOSR.174..405N ADS Malanushenko, O., Jones, H. P., & Livingston, W., “Imaging the Chromosphere Olson, G. L., Auer, L. H., & Buchler, J. R., “A rapidly convergent iterative solu- using Photospheric Mn 539.4 nm”, 2004IAUS..223..645M ADS tion of the non-LTE radiation transfer problem.”, 1986JQSRT..35..431O Maltby, P., Avrett, E. H., Carlsson, M., et al., “A New Sunspot Umbral Model ADS and Its Variation with the Solar Cycle”, 1986ApJ...306..284M ADS Ortiz, A., Hansteen, V. H., Nóbrega-Siverio, D., & Rouppe van der Voort, L., Martens, P. C. H., Kankelborg, C. C., & Berger, T. E., “On the Nature of the “Ellerman bombs and UV bursts: reconnection at different atmospheric lay- “Moss” Observed by TRACE”, 2000ApJ...537..471M ADS ers”, 2020A&A...633A..58O ADS Martínez-Sykora, J., De Pontieu, B., de la Cruz Rodriguéz, J., & Chintzoglou, Paletou, F., “On Milne-Barbier-Unsöld relationships”, G., “The Formation Height of Millimeter-wavelength Emission in the Solar 2018OAst...27...76P ADS Chromosphere”, 2020ApJ...891L...8M ADS Pariat, E., Masson, S., & Aulanier, G., “Current Buildup in Emerging Serpentine Mason, E. I., Antiochos, S. K., & Vourlidas, A., “An Observational Study of a Flux Tubes”, 2009ApJ...701.1911P ADS “Rosetta Stone” Solar Eruption”, 2021ApJ...914L...8M ADS Patsourakos, S., Gouttebroze, P., & Vourlidas, A., “The Quiet Sun Network at Mein, N., Mein, P., Heinzel, P., et al., “Cloud model with variable source func- Subarcsecond Resolution: VAULT Observations and Radiative Transfer Mod- tion for solar Hα structures.”, 1996A&A...309..275M ADS eling of Cool Loops”, 2007ApJ...664.1214P ADS Mein, P., “Champ macroscopique des vitesses dans l’atmosphère so- Pecker, J.-C., “Contribution à la théorie du type spectral: IV. La formation des laire d’après les mesures de déplacements des raies de Fraunhofer”, raies dans les spectres stellaires”, 1951AnAp...14..115P ADS 1966AnAp...29..153M ADS Compendium solar spectrum formation 41 References

Penn, M., Altrock, R. C., Henry, T., & Guhathakurta, M., “Synoptic Rutten, R. J. & van der Zalm, E. B. J., “Revision of solar equiva- Coronal Temperature; Magnetic Field and He 11083 NM Observations”, lent widths, Fe I oscillator strengths and the solar iron abundance.”, 1998ASPC..140..325P ADS 1984A&AS...55..143R ADS Pereira, T. M. D. & Uitenbroek, H., “RH 1.5D: a massively parallel code for Rutten, R. J. & Zwaan, C., “Empirical NLTE analyses of solar spectral lines. IV multi-level radiative transfer with partial frequency redistribution and Zee- - The Fe I curve of growth”, 1983A&A...117...21R ADS man polarisation”, 2015A&A...574A...3P ADS Rutten, R. G. M., Schrijver, C. J., Lemmens, A. F. P., & Zwaan, C., “Magnetic Pereira, T. M. D., De Pontieu, B., & Carlsson, M., “Quantifying Spicules”, structure in cool stars. XVII. Minimum radiative losses from theouter atmo- 2012ApJ...759...18P ADS sphere.”, 1991A&A...252..203R ADS Pereira, T. M. D., De Pontieu, B., Carlsson, M., et al., “An In- Rutten, R. J., van Veelen, B., & Sütterlin, P., “DOT Tomography of terface Region Imaging Spectrograph First View on Solar Spicules”, the Solar Atmosphere VII. Chromospheric Response to Acoustic Events”, 2014ApJ...792L..15P ADS 2008SoPh..251..533R ADS Pereira, T. M. D., “The dynamic chromosphere: Pushing the boundaries of ob- Rutten, R. J., Leenaarts, J., Rouppe van der Voort, L. H. M., et al., “Quiet-Sun servations and models”, 2019AdSpR..63.1434P ADS imaging asymmetries in Na I D_1 compared with other strong Fraunhofer Peter, H., Tian, H., Curdt, W., et al., “Hot explosions in the cool atmosphere of lines”, 2011A&A...531A..17R ADS the Sun”, 2014Sci...346C.315P ADS Rutten, R. J., Vissers, G. J. M., Rouppe van der Voort, L. H. M., Süt- Pierce, A. K. & Breckenridge, J. B.: 1974, The Kitt Peak table of photographic terlin, P., & Vitas, N., “Ellerman bombs: fallacies, fads, usage”, solar spectrum wavelengths 1974kptp.book.....P ADS 2013JPhCS.440a2007R ADS Pietarila Graham, J., Cameron, R., & Schüssler, M., “Turbulent Small-Scale Rutten, R. J., Rouppe van der Voort, L. H. M., & Vissers, G. J. M., Dynamo Action in Solar Surface Simulations”, 2010ApJ...714.1606P “Ellerman Bombs at High Resolution. IV. Visibility in Na I and Mg I”, ADS 2015ApJ...808..133R ADS Polito, V., De Pontieu, B., Testa, P., Brooks, D. H., & Hansteen, V., “IRIS Ob- Rutten, R. J., Rouppe van der Voort, L. H. M., & De Pontieu, B., “Solar Hα fea- servations of the Low-atmosphere Counterparts of Active Region Outflows”, tures with hot onsets. IV. Network fibrils”, 2019A&A...632A..96R ADS 2020ApJ...903...68P ADS Rutten, R. J., “Empirical NLTE analyses of solar spectral lines. II - The forma- Pottasch, S. R., “On the Interpretation of the Solar Ultraviolet Emission Line tion of the Ba II λ 4554 resonance line”, 1978SoPh...56..237R ADS Spectrum”, 1964SSRv....3..816P ADS Rutten, R. J., “The NLTE formation of iron lines in the solar photosphere”, Priest, E. R., Heyvaerts, J. F., & Title, A. M., “A Flux-Tube Tecton- 1988ASSL..138..185R ADS ics Model for Solar Coronal Heating Driven by the Magnetic Carpet”, Rutten, R. J., “Summary Lecture”, 1990IAUS..138..501R ADS 2002ApJ...576..533P ADS Rutten, R. J., “Chromospheric Oscillations”, 1995ESASP.376a.151R ADS Reames, D. V.: 2021, Solar Energetic Particles. A Modern Primer Rutten, R. J., “(Inter-),Network Structure and Dynamics”, on Understanding Sources, Acceleration and Propagation, Vol. 978 1999ASPC..184..181R ADS 2021LNP...978.....R ADS Rutten, R. J., “Solar Atmosphere Models”, 2002JAD.....8....8R ADS Reardon, K. & Cauzzi, G., “New insight on the coupling of the solar atmosphere Rutten, R. J., “On the Nature of the Solar Chromosphere”, from imaging spectroscopy”, 2012decs.confE..20R ADS 2006ASPC..354..276R ADS Reardon, K. P. & Cauzzi, G., “The IBIS Mosaic - A Broad View Of The Solar Rutten, R. J., “Observing the Solar Chromosphere”, 2007ASPC..368...27R Atmosphere”, 2012AAS...22020111R ADS ADS Rezaei, R. & Beck, C., “Multiwavelength spectropolarimetric observations of Rutten, R. J., “The solar activity sensitivity of Mn I lines”, 2011bandung ADS an Ellerman bomb”, 2015A&A...582A.104R ADS Rutten, R. J., “Twists to Solar Spicules”, 2013ASPC..470...49R ADS Rosner, R., Tucker, W. H., & Vaiana, G. S., “Dynamics of the quiescent solar Rutten, R. J., “Hα features with hot onsets. I. Ellerman bombs”, corona.”, 1978ApJ...220..643R ADS 2016A&A...590A.124R ADS Rouppe van der Voort, L. & de la Cruz Rodríguez, J., “Short Dynamic Fibrils in Rutten, R. J., “Solar ALMA predictions: tutorial”, 2017IAUS..327....1R Sunspot Chromospheres”, 2013ApJ...776...56R ADS ADS Rouppe van der Voort, L. H. M., Rutten, R. J., Sütterlin, P., Sloover, Rutten, R. J., “Solar H-alpha features with hot onsets. III. Long fibrils in Lyman- P. J., & Krijger, J. M., “La Palma observations of umbral flashes”, alpha and with ALMA”, 2017A&A...598A..89R ADS 2003A&A...403..277R ADS Rutten, R., “A year at Sunspot. Reminiscences of 1977-1978 at the Sacramento Rouppe van der Voort, L. H. M., De Pontieu, B., Hansteen, V. H., Carlsson, M., & Peak Observatory”, 2018arXiv180408709R ADS van Noort, M., “Magnetoacoustic Shocks as a Driver of Quiet-Sun Mottles”, Rutten, R. J., “Non-Equilibrium Spectrum Formation Affecting Solar Irradi- 2007ApJ...660L.169R ADS ance”, 2019SoPh..294..165R ADS Rouppe van der Voort, L., Leenaarts, J., de Pontieu, B., Carlsson, M., & Vissers, Rutten, R. J., “Small-scale solar surface magnetism”, 2021arXiv210514533R G., “On-disk Counterparts of Type II Spicules in the Ca II 854.2 nm and Hα ADS Lines”, 2009ApJ...705..272R ADS Rutten, R. J., “SolO campfires in SDO images”, 2020arXiv200900376R Rouppe van der Voort, L. H. M., Rutten, R. J., & Vissers, G. J. M., “Reconnec- ADS tion brightenings in the quiet solar photosphere”, 2016A&A...592A.100R Rybicki, G. B. & Hummer, D. G., “An accelerated lambda iteration method ADS for multilevel radiative transfer. II. Overlapping transitions with full contin- Rutten, R. J. & Carlsson, M., “The Formation of Infrared Rydberg Lines”, uum.”, 1992A&A...262..209R ADS 1994IAUS..154..309R ADS Rybicki, G. B. & Lightman, A. P.: 1986, Radiative Processes in Astrophysics, Rutten, R. J. & Kostik, R. I., “Empirical NLTE analyses of solar spec- Wiley, New York 1986rpa..book.....R ADS tral lines. III - Iron lines versus LTE models of the photosphere”, Samain, D., “Solar continuum data on absolute intensities, center to 1982A&A...115..104R ADS limb variations and Laplace inversion between 1400 and 2100 Å.”, Rutten, R. J. & Krijger, J. M., “Dynamics of the solar chromosphere IV. Evidence 1979A&A....74..225S ADS for atmospheric gravity waves from TRACE”, 2003A&A...407..735R Scharmer, G. B. & Carlsson, M., “A new approach to multi-level non-LTE ra- ADS diative transfer problems.”, 1985JCoPh..59...56S ADS Rutten, R. J. & Milkey, R. W., “Partial redistribution in the solar photospheric Scharmer, G. B., Löfdahl, M. G., Sliepen, G., & de la Cruz Rodríguez, J., “Is Ba II spectrum”, 1979ApJ...231..277R ADS the sky the limit?. Performance of the revamped Swedish 1-m Solar Telescope Rutten, R. J. & Rouppe van der Voort, L. H. M., “Solar Hα features with hot and its blue- and red-beam reimaging systems”, 2019A&A...626A..55S onsets. II. A contrail fibril”, 2017A&A...597A.138R ADS ADS Rutten, R. J. & Schrijver, C. J. (Eds.): 1994, Solar Surface Magnetism, Procs. Scharmer, G. B., “The Wilson-Bappu effect as a result of supersonic turbu- Soesterberg workshop, NATO ASI Series C 433 1994ASIC..433.....R lence.”, 1976A&A....53..341S ADS ADS Schoolman, S. A., “Formation of the Solar Hα Profile”, Rutten, R. J. & Stencel, R. E., “Solar limb emission lines near CA II H & K and 1972SoPh...22..344S ADS their spatial intensity variations”, 1980A&AS...39..415R ADS Schrijver, C. J. & Zwaan, C.: 2008, Solar and Stellar Magnetic Activity Rutten, R. J. & Uitenbroek, H., “Ca II H2v and K2v Cell Grains”, 2008ssma.book.....S ADS 1991SoPh..134...15R ADS Schrijver, C. J., Mewe, R., & Zwaan, C., “Coronal activity in F-, G-, and K-type Rutten, R. J. & Uitenbroek, H., “Chromospheric backradiation in ultraviolet stars; relations between parameters characterizing stellar structures and X- continua and Hα”, 2012A&A...540A..86R ADS ray emission”, 1982AdSpR...2i.243S ADS Rutten, R. J. & van Amerongen, H. J., “A very simple digital microdensitometer- Schrijver, C. J., Dobson, A. K., & Radick, R. R., “The Magnetic, Basal, and comparator”, 1975ASSL...54..261R ADS Radiative Equilibrium Components in Mount Wilson CA II H+K Fluxes”, Rutten, R. J. & van der Zalm, E. B. J., “Clean lines in the solar flux spectrum”, 1989ApJ...341.1035S ADS 1984A&AS...55..171R ADS Compendium solar spectrum formation 42 References

Schrijver, C. J., Title, A. M., van Ballegooijen, A. A., Hagenaar, H. J., & Uitenbroek, H. & Criscuoli, S., “Why One-dimensional Models Fail in the Shine, R. A., “Sustaining the Quiet Photospheric Network: The Bal- Diagnosis of Average Spectra from Inhomogeneous Stellar Atmospheres”, ance of Flux Emergence, Fragmentation, Merging, and Cancellation”, 2011ApJ...736...69U ADS 1997ApJ...487..424S ADS Uitenbroek, H., “Operator perturbation method for multi-level line transfer with Schrijver, C. J.: 1987a, Heating of Stellar Chromospheres and Coronae: partial redistribution”, 1989A&A...213..360U ADS Evidence for Non-Magnetic Heating (Invited review), Vol. 291, 135–145 Uitenbroek, H., “The CO Fundamental Vibration-Rotation Lines in the Solar 1987LNP...291..135S ADS Spectrum. I. Imaging Spectroscopy and Multidimensional LTE Modeling”, 2000ApJ...531..571U ADS Schrijver, C. J., “Magnetic structure in cool stars. XI. Relations between ra- Uitenbroek, H., “Multilevel Radiative Transfer with Partial Frequency Redistri- diative fluxes mesuring stellar activity, and evidence for two components in bution”, 2001ApJ...557..389U ADS stellar chromospheres.”, 1987A&A...172..111S ADS Ulmschneider, P. & Kalkofen, W.: 2003, Heating of the solar chromosphere, Schrijver, C. J., “The basal and strong-field components of the solar atmo- 181–195 2003dysu.book..181U ADS sphere”, 1992A&A...258..507S ADS Ulrich, R. K., “The Five-Minute Oscillations on the Solar Surface”, Schrijver, C. J., “Basal heating in the atmospheres of cool stars”, 1970ApJ...162..993U ADS 1995A&ARv...6..181S ADS Unruh, Y. C., Solanki, S. K., & Fligge, M., “The spectral dependence of facular Sekse, D. H., Rouppe van der Voort, L., & De Pontieu, B., “Statistical contrast and solar irradiance variations”, 1999A&A...345..635U ADS Properties of the Disk Counterparts of Type II Spicules from Simultane- Unsold, A.: 1955, Physik der Sternatmospharen, MIT besonderer Berucksichti- gung der Sonne. 1955psmb.book.....U ADS ous Observations of Rapid Blueshifted Excursions in Ca II 8542 and Hα”, van de Hulst, H. C.: 1953, The Chromosphere and the Corona, 207 2012ApJ...752..108S ADS 1953sun..book..207V ADS Sekse, D. H., Rouppe van der Voort, L., De Pontieu, B., & Scullion, E., “Inter- van Noort, M., Rouppe van der Voort, L., & Löfdahl, M. G., “Solar Image play of Three Kinds of Motion in the Disk Counterpart of Type II Spicules: Restoration By Use Of Multi-frame Blind De-convolution With Multiple Ob- Upflow, Transversal, and Torsional Motions”, 2013ApJ...769...44S jects And Phase Diversity”, 2005SoPh..228..191V ADS ADS Vernazza, J. E., Avrett, E. H., & Loeser, R., “Structure of the solar chromo- Shchukina, N. & Trujillo Bueno, J., “The Iron Line Formation Problem sphere. II. The underlying photosphere and temperature-minimum region.”, in Three-dimensional Hydrodynamic Models of Solar-like Photospheres”, 1976ApJS...30....1V ADS 2001ApJ...550..970S ADS Vernazza, J. E., Avrett, E. H., & Loeser, R., “Structure of the solar chromo- Shchukina, N. G., Olshevsky, V. L., & Khomenko, E. V., “The solar BaII 4554 Å sphere. III. Models of the EUV brightness components of the quiet sun.”, line as a Doppler diagnostic: NLTE analysis in 3D hydrodynamical model”, 1981ApJS...45..635V ADS , “Eller- 2009A&A...506.1393S ADS Vissers, G. J. M., Rouppe van der Voort, L. H. M., & Rutten, R. J. Sheminova, V. A., Rutten, R. J., & Rouppe van der Voort, L. H. M., “The wings man Bombs at High Resolution. II. Triggering, Visibility, and Effect on Upper Atmosphere”, 2013ApJ...774...32V ADS of Ca II H and K as solar fluxtube diagnostics”, 2005A&A...437.1069S Vissers, G. J. M., Rouppe van der Voort, L. H. M., Rutten, R. J., Carlsson, M., ADS & De Pontieu, B., “Ellerman Bombs at High Resolution. III. Simultaneous Shine, R. A., Milkey, R. W., & Mihalas, D., “Resonance Line Transfer with Observations with IRIS and SST”, 2015ApJ...812...11V ADS Partial Redistribution. V. The Solar CA II Lines”, 1975ApJ...199..724S Vissers, G. J. M., Rouppe van der Voort, L. H. M., & Rutten, R. J., “Automating ADS Ellerman bomb detection in ultraviolet continua”, 2019A&A...626A...4V Simões, P. J. A., Reid, H. A. S., Milligan, R. O., & Fletcher, L., “The Spectral ADS Content of SDO/AIA 1600 and 1700 Å Filters from Flare and Plage Observa- Vitas, N., Viticchiè, B., Rutten, R. J., & Vögler, A., “Explanation of the activity tions”, 2019ApJ...870..114S ADS sensitivity of Mn I 5394.7 Å”, 2009A&A...499..301V ADS Simon, G. W., Title, A. M., & Weiss, N. O., “Sustaining the Sun’s Magnetic Vitense, E., “Der Aufbau der Sternatmosphären. IV. Teil. Kontinuierliche Ab- Network with Emerging Bipoles”, 2001ApJ...561..427S ADS sorption und Streuung als Funktion von Druck und Temperatur. Mit 16 Skogsrud, H., Rouppe van der Voort, L., & De Pontieu, B., “On the Active Textabbildungen”, 1951ZA.....28...81V ADS Region Bright Grains Observed in the Transition Region Imaging Channels Vögler, A., Shelyag, S., Schüssler, M., et al., “Simulations of magneto- of IRIS”, 2016ApJ...817..124S ADS convection in the solar photosphere. Equations, methods, and results of the MURaM code”, 2005A&A...429..335V ADS Smitha, H. N., Holzreuter, R., van Noort, M., & Solanki, S. K., “The influence von Steiger, R. & Geiss, J., “Supply of fractionated gases to the corona”, of NLTE effects in Fe I lines on an inverted atmosphere. I. 6301 Å and 6302 1989A&A...225..222V ADS Å lines formed in 1D NLTE”, 2020A&A...633A.157S ADS Vourlidas, A., Sánchez Andrade-Nuño, B., Landi, E., et al., “The Struc- Smitha, H. N., Holzreuter, R., van Noort, M., & Solanki, S. K., “The influence ture and Dynamics of the Upper Chromosphere and Lower Transi- of NLTE effects in Fe I lines on an inverted atmosphere. II. 6301 Å and 6302 tion Region as Revealed by the Subarcsecond VAULT Observations”, Å lines formed in 3D NLTE”, 2021A&A...647A..46S ADS 2010SoPh..261...53V ADS Straus, T., Fleck, B., Jefferies, S. M., et al., “The Energy Flux of Internal Gravity Wang, Y. M., Sheeley, N. R., Dere, K. P., et al., “Association of Extreme- Waves in the Lower Solar Atmosphere”, 2008ApJ...681L.125S ADS Ultraviolet Imaging Telescope (EIT) Polar Plumes with Mixed-Polarity Mag- Suematsu, Y., Wang, H., & Zirin, H., “High-Resolution Observation of Disk netic Network”, 1997ApJ...484L..75W ADS Spicules. I. Evolution and Kinematics of Spicules in the Enhanced Network”, Watanabe, H., Vissers, G., Kitai, R., Rouppe van der Voort, L., & Rutten, R. J., 1995ApJ...450..411S ADS “Ellerman Bombs at High Resolution. I. Morphological Evidence for Photo- Sukhorukov, A. V. & Leenaarts, J., “Partial redistribution in 3D non-LTE radia- spheric Reconnection”, 2011ApJ...736...71W ADS tive transfer in solar-atmosphere models”, 2017A&A...597A..46S ADS Webb, D. F., Gibson, S. E., Hewins, I. M., et al., “Preserving a Unique Archive for Long-Term Solar Variability Studies”, 2017SpWea..15.1442W ADS Sütterlin, P., Rutten, R. J., & Skomorovsky, V. I., “Ba II 4554 Å speckle imaging Whitaker, W. A., “Heating of the Solar Corona by Gravity Waves.”, 2001A&A...378..251S as solar Doppler diagnostic”, ADS 1963ApJ...137..914W ADS Thomas, R. N., “The Source Function in a Non-Equilibrium Atmosphere. I. The White, O. R. & Cha, M. Y., “Analysis of the 5 min Oscillatory Photospheric Resonance Lines.”, 1957ApJ...125..260T ADS Motion. I: A Problem in Waveform Classification”, 1973SoPh...31...23W Thomas, R. N., “The Source Function in a Non-Equilibrium Atmo- ADS sphere.IV. Evaluation and Applications of the Net Radiative Bracket.”, Wiersma, J., Rutten, R. J., & Lanz, T., “NLTE in a Hot Hydrogen Star: Auer & 1960ApJ...131..429T ADS Mihalas Revisited”, 2003ASPC..288..130W ADS Thomas, R. N.: 1983, Stellar atmospheric structural patterns, Vol. 471 Wijbenga, J. W. & Zwaan, C., “Empirical NLTE Analyses of Solar Spec- 1983NASSP.471.....T ADS tral Lines. I: A Method and Some Applications to Earlier Analyses”, Title, A. M.: 1966, Selected Spectroheliograms, Mount Wilson and Palomar 1972SoPh...23..265W ADS Observatories, Pasadena Title1966a ADS Wilson, O. C. & Vainu Bappu, M. K., “H and K Emission in Late-Type Title, A. M.: 1966, “A Study of Velocity Fields in the Hydrogen-Alpha Chromo- Stars: Dependence of Line Width on Luminosity and Related Topics.”, sphere by Means of Time-Lapse Doppler Movies. Ph.D. thesis 1957ApJ...125..661W ADS ”„ , CALIFOR- Withbroe, G. L., “The thermal phase of a large solar flare.”, 1966PhDT...... 1T NIA INSTITUTE OF TECHNOLOGY. ADS 1978ApJ...225..641W ADS Trujillo Bueno, J., Shchukina, N., & Asensio Ramos, A., “A substantial amount Wolff, C. L., “The Five-Minute Oscillations as Nonradial Pulsations of the En- of hidden magnetic energy in the quiet Sun”, 2004Natur.430..326T ADS tire Sun”, 1972ApJ...177L..87W ADS Turck-Chièze, S., Carton, P. H., Barrière, J. C., et al., “Solar Global Oscillations Wolff, C. L., “What is the Horizontal Scale of the 5-min Oscillations?”, of Low-Degree modes (GOLD): The Status of the Multi-channel Resonance 1973SoPh...32...31W ADS Spectrometer GOLF-NG”, 2012ASPC..462..240T ADS Young, P. R., Tian, H., Peter, H., et al., “Solar Ultraviolet Bursts”, Tziotziou, K., “Chromospheric Cloud-Model Inversion Techniques”, 2018SSRv..214..120Y ADS 2007ASPC..368..217T ADS Zirin, H. & Popp, B., “Observations of the 12 Micron MG i Lines in Various , “The formation of helioseismology Solar Features”, 1989ApJ...340..571Z ADS Uitenbroek, H. & Bruls, J. H. M. J. Zwaan, C.: 1986, Relations Between Magnetic Activity and Stellar Properties, lines. III. Partial redistribution effects in weak solar resonance lines. ”, Vol. 254, 19 1986LNP...254...19Z ADS 1992A&A...265..268U ADS