S. K. Antiochos

The Dynamic Formation of Prominence Condensations

We present simulations of a model for the formation of a prominence condensation in a coronal loop. The key idea behind the model is that the spatial localization of loop heating near the chromosphere leads to a catastrophic cooling in the corona. Using a new adaptive grid code, we simulate the complete growth of a condensation and find that after ~5000 s it reaches a quasi-steady state. We show that the size and growth time of the condensation are in good agreement with data and discuss the implications of the model for coronal heating and for observations of prominences and the surrounding corona.

Dynamical Formation and Stability of Helical Prominence Magnetic Fields

Abstract : We numerically simulated an initially bipolar magnetic field subjected to shear motions concentrated near and parallel to the photospheric polarity inversion line. The simulations yield three principal results: (1) For footprint displacements comparable to the bipoles depth, the sheared core field acquires a dipped geometry that can support cool prominence material against gravity. This confirms previous force-free equilibrium models for forming dipped prominence fields by differential shear, and extends them to much larger applied shears and time-dependent dynamics with dissipation. (2) At larger shears, we discover a new mechanism for forming the helical magnetic fields of prominences. It entails a two-step process of magnetic reconnection in the corona. First, flux in the sheared core reconnects with flux in the unsheared, restraining arcade, producing new pairs of interlinked field lines. Second, as these interlinked fields continue to be sheared, they are brought together and reconnect again, producing helical field threading and enveloping the body of the prominence. This mechanism can account for the twist that is often observed in both quiescent and erupting prominences. (3) Even for very large shears, the dipped, helical structure settles into an apparently stable equilibrium, despite the substantial amount of reconnection and twist in the magnetic field. We conclude that neither a kink instability of the helical core field, nor a tether-cutting instability of the restraining arcade, is operating in our low-lying model prominence. This concurs with both observations and a theoretical model for prominence stability.

A Numerical Study of the Breakout Model for Coronal Mass Ejection Initiation

A leading theory for the initiation of coronal mass ejections (CMEs) is the breakout model, in which magnetic reconnection above a filament channel is responsible for disrupting the coronal magnetic field. We present the first simulations of the complete breakout process including the initiation, the plasmoid formation and ejection, and the eventual relaxation of the coronal field to a more potential state. These simulations were performed using a new numerical code that solves the numerical cavitation problems that prevented previous simulations from calculating a complete ejection. Furthermore, the position of the outer boundary in the new simulations is increased out to 30 R☉, which enables determination of the full structure and dynamics of the ejected plasmoid. Our results show that the ejection occurs at a speed on the order of the coronal Alfven speed and hence that the breakout model can produce fast CMEs. Another key result is that the ejection speed is not sensitive to the refinement level of the grid used in the calculations, which implies that, at least for the numerical resistivity of these simulations, the speed is not sensitive to the Lundquist number. We also calculate, in detail, the helicity of the system and show that the helicity is well conserved during the breakout process. Most of the helicity is ejected from the Sun with the escaping plasmoid, but a significant fraction (of order 10%) remains in the corona. The implications of these results for observation and prediction of CMEs and eruptive flares is discussed.

Reconnection of Twisted Flux Tubes as a Function of Contact Angle

The collision and reconnection of magnetic flux tubes in the solar corona has been proposed as a mechanism for solar flares and in some cases as a model for coronal mass ejections. We study this process by simulating the collision of pairs of twisted flux tubes with a massively parallel, collocation, viscoresistive, magnetohydrodynamic code using up to 256 × 256 × 256 Fourier modes. Our aim is to investigate the energy release and possible global topological changes that can occur in flux-tube reconnection. We have performed a number of simulations for different angles between the colliding flux tubes and for either co- or counterhelicity flux tubes. We find the following four classes of interaction: (1) bounce (no appreciable reconnection), (2) merge, (3) slingshot (the most efficient reconnection), and (4) tunnel (a double reconnection). We will describe these four classes of flux-tube reconnection and discuss in what range of parameter space each class occurs and the implications our results have for models of flares and coronal mass ejections.

SOLAR PROMINENCE INTERACTIONS

We report numerical simulations of the formation, interaction, and magnetic reconnection between pairs of solar prominences within the sheared-arcade model. Our experiments consider the four possible basic combinations of chiralities (identical or opposite) and axial magnetic fields (aligned or opposed) between the participating prominences. When the topology of the global flux system comprising the prominences and arcades is bipolar, so that a single polarity inversion line is shared by the two structures, then identical chiralities necessarily imply aligned axial fields, while opposite chiralities imply opposed axial fields. In the former case, external magnetic reconnections forming field lines linking the two prominences occur; in the latter, such reconnections are disfavored, and no linkage takes place. These results concur with empirical rules for prominence interactions. When the topology instead is quadrupolar, so that a second polarity inversion line crossing the first lies between the prominences, then the converse relation holds between chirality and axial-field alignment. External reconnections forming linking field lines now occur between prominences with opposite chiralities; they also occur, but result only in footpoint exchanges, between prominences with identical chiralities. These findings conflict with the accepted empirical rules but may not have been tested in observations to date. All of our model prominences, especially those that undergo linking reconnections, contain substantial magnetic shear and twist. Nevertheless, none exhibits any sign of onset of instability or loss of equilibrium that might culminate in an eruption.

The Origin of High-Speed Motions and Threads in Prominences

Prominences are among the most spectacular manifestations of both quiescent and eruptive solar activity, yet the origins of their magnetic-field and plasma structures remain poorly understood. We have made steady progress toward a comprehensive model of prominence formation and evolution with our sheared three-dimensional arcade model for the magnetic field and our thermal nonequilibrium model for the cool, dense material suspended in the corona. According to the thermal nonequilibrium model, condensations form readily along long, low-lying magnetic field lines when the heating is localized near the chromosphere. In most cases this process yields a dynamic cycle in which condensations repetitively form, stream along the field, and ultimately disappear by falling onto the nearest footpoint. Two key observed features were not adequately explained by our earlier simulations of thermal nonequilibrium, however: the threadlike (i.e., elongated) horizontal structure and high-speed motions of many condensations. In this paper we discuss how simple modifications to the radiative loss function, the heating scale, and the geometry of our model largely eliminate these discrepancies. In particular, condensations in nearly horizontal flux tubes are most likely to develop both transient high-speed motions and elongated threads. These results strengthen the case for thermal nonequilibrium as the origin of prominence condensations and support low-twist models of prominence magnetic structure.

Constraints on Active Region Coronal Heating

We derive constraints on the time variability of coronal heating from observations of the so-called active region moss by the Transition Region and Coronal Explorer (TRACE). The moss is believed to be due to million-degree emission from the transition regions at the footpoints of coronal loops whose maximum temperatures are several million degrees. The two key results from the TRACE observations discussed in this paper are that in the moss regions one generally sees only moss, not million-degree loops, and that the moss emission exhibits only weak intensity variations, � 10% over periods of hours. TRACE movies showing these results are presented. We demonstrate, using both analytic and numerical calculations, that the lack of observable million-degree loops in the moss regions places severe constraints on the possible time variability of coronal heating in the loops overlying the moss. In particular, the heating in the hot moss loops cannot be truly flarelike with a sharp cutoff, but instead must be quasi-steady to an excellent approximation. Furthermore, the lack of significant variations in the moss intensity implies that the heating magnitude is only weakly varying. The implications of these conclusions for coronal heating models will be discussed. Subject headings: Sun: corona — Sun: transition region — Sun: UV radiation On-line material: mpg animation

Structure and Dynamics of the Sun’s Open Magnetic Field

The solar magnetic field is the primary agent that drives solar activity and couples the Sun to the heliosphere. Although the details of this coupling depend on the quantitative properties of the field, many important aspects of the corona-solar wind connection can be understood by considering only the general topological properties of those regions on the Sun where the field extends from the photosphere out to interplanetary space, the so-called open field regions that are usually observed as coronal holes. From the simple assumptions that underlie the standard quasi-steady corona-wind theoretical models, and that are likely to hold for the Sun as well, we derive two conjectures as to the possible structure and dynamics of coronal holes: (1) coronal holes are unique in that every unipolar region on the photosphere can contain at most one coronal hole, and (2) coronal holes of nested polarity regions must themselves be nested. Magnetic reconnection plays the central role in enforcing these constraints on the field topology. From these conjectures we derive additional properties for the topology of open field regions, and propose several observational predictions for both the slowly varying and transient corona/solar wind.

A Transient Heating Model for Coronal Structure and Dynamics

A wealth of observational evidence for flows and intensity variations in nonflaring coronal loops leads to the conclusion that coronal heating is intrinsically unsteady and concentrated near the chromosphere. We have investigated the hydrodynamic behavior of coronal loops undergoing transient heating with one-dimensional numerical simulations in which the timescale assumed for the heating variations (3000 s) is comparable to the coronal radiative cooling time and the assumed heating location and scale height (10 Mm) are consistent with the values derived from TRACE studies. The model loops represent typical active region loops: 40-80 Mm in length, reaching peak temperatures up to 6 MK. We use ARGOS, our state-of-the-art numerical code with adaptive mesh refinement, in order to resolve adequately the dynamic chromospheric-coronal transition region sections of the loop. The major new results from our work are the following: (1) During much of the cooling phase, the loops exhibit densities significantly larger than those predicted by the well-known loop scaling laws, thus potentially explaining recent TRACE observations of overdense loops. (2) Throughout the transient heating interval, downflows appear in the lower transition region (T ~ 0.1 MK) whose key signature would be persistent, redshifted UV and EUV line emission, as have long been observed. (3) Strongly unequal heating in the two legs of the loop drives siphon flows from the more strongly heated footpoint to the other end, thus explaining the substantial bulk flows in loops recently observed by the Coronal Diagnostic Spectrometer and the Solar Ultraviolet Measurement of Emission Radiation instrument. We discuss the implications of our studies for the physical origins of coronal heating and related dynamic phenomena.

Condensation Formation by Impulsive Heating in Prominences

Our thermal nonequilibrium model for prominence formation provides an explanation for the well-observed presence of predominantly dynamic, cool, dense material suspended in the corona above filament channels. According to this model, condensations form readily along long, low-lying magnetic field lines when heating is localized near the chromosphere. Often this process yields a dynamic cycle in which condensations repeatedly form, stream along the field, and ultimately disappear by falling onto the nearest footpoint. Our previous studies employed only steady heating, as is consistent with some coronal observations, but many coronal heating models predict transient episodes of localized energy release (e.g., nanoflares). Here we present the results of a numerical investigation of impulsive heating in a model prominence flux tube and compare the outcome with previous steady-heating simulations. We find that condensations form readily when the average interval between heating events is less than the coronal radiative cooling time (~2000 s). As the average interval between pulses decreases, the plasma evolution more closely resembles the steady-heating case. The heating scale and presence or absence of background heating also determine whether or not condensations form and how they evolve. Our results place important constraints on coronal heating in filament channels and strengthen the case for thermal nonequilibrium as the process responsible for the plasma structure in prominences.