Search

Observations have shown a clear association of filament/prominence eruptions with the emergence of magnetic flux in or near filament channels. Magnetohydrodynamic (MHD) simulations have been employed to systematically study the conditions under which such eruptions occur. These simulations to date have modeled filament channels as two-dimensional (2D) flux ropes or 3D uniformly sheared arcades. Here we present MHD simulations of flux emergence into a more realistic configuration consisting of a bipolar active region containing a line-tied 3D flux rope. We use the coronal flux-rope model of Titov et al. (2014) as the initial condition and drive our simulations by imposing boundary conditions extracted from a flux-emergence simulation by Leake et al. (2013). We identify three mechanisms that determine the evolution of the system: (i) reconnection displacing foot points of field lines overlying the coronal flux rope, (ii) changes of the ambient field due to the intrusion of new flux at the boundary, and (iii) interaction of the (axial) electric currents in the pre-existing and newly emerging flux systems. The relative contributions and effects of these mechanisms depend on the properties of the pre-existing and emerging flux systems. Here we focus on the location and orientation of the emerging flux relative to the coronal flux rope. Varying these parameters, we investigate under which conditions an eruption of the latter is triggered., Comment: 18 pages, 8 figures, accepted for publication by The Astrophysical Journal

The topology of the coronal magnetic field produces a strong impact on the properties of the solar corona and presumably on the origin of the slow solar wind. To advance our understanding of this impact, we revisit the concept of so-called slip-back mapping (Titov et al. 2009) and adapt it for determining open, closed, and disconnected flux systems that are formed in the solar corona by magnetic reconnection during a given time interval. The developed method allows us, in particular, to describe the magnetic flux transfer between open and closed flux regions via so-called interchange reconnection with an unprecedented level of details. We illustrate the application of this method to the analysis of the global MHD evolution of the solar corona that is driven by an idealized differential rotation of the photospheric plasma., Comment: To appear in Ap.J

GPU accelerators have had a notable impact on high-performance computing across many disciplines. They provide high performance with low cost/power, and therefore have become a primary compute resource on many of the largest supercomputers. Here, we implement multi-GPU acceleration into our Solar MHD code (MAS) using OpenACC in a fully portable, single-source manner. Our preliminary implementation is focused on MAS running in a reduced physics "zero-beta" mode. While valuable on its own, our main goal is to pave the way for a full physics, thermodynamic MHD implementation. We describe the OpenACC implementation methodology and challenges. "Time-to-solution" performance results of a production-level flux rope eruption simulation on multi-CPU and multi-GPU systems are shown. We find that the GPU-accelerated MAS code has the ability to run "zero-beta" simulations on a single multi-GPU server at speeds previously requiring multiple CPU server-nodes of a supercomputer., Comment: 13 pages, 9 figures

Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models., Comment: 29 pages, 11 figures, accepted for publication in Space Science Reviews

Long-period EUV pulsations, recently discovered to be common in active regions, are understood to be the coronal manifestation of thermal non-equilibrium (TNE). The active regions previously studied with EIT/SOHO and AIA/SDO indicated that long-period intensity pulsations are localized in only one or two loop bundles. The basic idea of this study is to understand why. For this purpose, we tested the response of different loop systems, using different magnetic configurations, to different stratifications and strengths of the heating. We present an extensive parameter-space study using 1D hydrodynamic simulations (1,020 in total) and conclude that the occurrence of TNE requires specific combinations of parameters. Our study shows that the TNE cycles are confined to specific ranges in parameter space. This naturally explains why only some loops undergo constant periodic pulsations over several days: since the loop geometry and the heating properties generally vary from one loop to another in an active region, only the ones in which these parameters are compatible exhibits TNE cycles. Furthermore, these parameters (heating and geometry) are likely to vary significantly over the duration of a cycle, which potentially limits the possibilities of periodic behavior. This study also confirms that long-period intensity pulsations and coronal rain are two aspects of the same phenomenon: both phenomena can occur for similar heating conditions and can appear simultaneously in the simulations., Comment: 26 pages, 17 figures, accepted for publication in The Astrophysical journal

The heliospheric magnetic field is of pivotal importance in solar and space physics. The field is rooted in the Sun's photosphere, where it has been observed for many years. Global maps of the solar magnetic field based on full disk magnetograms are commonly used as boundary conditions for coronal and solar wind models. Two primary observational constraints on the models are (1) the open field regions in the model should approximately correspond to coronal holes observed in emission, and (2) the magnitude of the open magnetic flux in the model should match that inferred from in situ spacecraft measurements. In this study, we calculate both MHD and PFSS solutions using fourteen different magnetic maps produced from five different types of observatory magnetograms, for the time period surrounding July, 2010. We have found that for all of the model/map combinations, models that have coronal hole areas close to observations underestimate the interplanetary magnetic flux, or, conversely, for models to match the interplanetary flux, the modeled open field regions are larger than coronal holes observed in EUV emission. In an alternative approach, we estimate the open magnetic flux entirely from solar observations by combining automatically detected coronal holes for Carrington rotation 2098 with observatory synoptic magnetic maps. This approach also underestimates the interplanetary magnetic flux. Our results imply that either typical observatory maps underestimate the Sun's magnetic flux, or a significant portion of the open magnetic flux is not rooted in regions that are obviously dark in EUV and X-ray emission.

Transient collimated plasma eruptions in the solar corona, commonly known as coronal (or X-ray) jets, are among the most interesting manifestations of solar activity. It has been suggested that these events contribute to the mass and energy content of the corona and solar wind, but the extent of these contributions remains uncertain. We have recently modeled the formation and evolution of coronal jets using a three-dimensional (3D) magnetohydrodynamic (MHD) code with thermodynamics in a large spherical domain that includes the solar wind. Our model is coupled to 3D MHD flux-emergence simulations, i.e, we use boundary conditions provided by such simulations to drive a time-dependent coronal evolution. The model includes parametric coronal heating, radiative losses, and thermal conduction, which enables us to simulate the dynamics and plasma properties of coronal jets in a more realistic manner than done so far. Here we employ these simulations to calculate the amount of mass and energy transported by coronal jets into the outer corona and inner heliosphere. Based on observed jet-occurrence rates, we then estimate the total contribution of coronal jets to the mass and energy content of the solar wind to (0.4-3.0) % and (0.3-1.0) %, respectively. Our results are largely consistent with the few previous rough estimates obtained from observations, supporting the conjecture that coronal jets provide only a small amount of mass and energy to the solar wind. We emphasize, however, that more advanced observations and simulations are needed to substantiate this conjecture.

Coronal jets are transient, collimated eruptions that occur in regions of predominantly open magnetic field in the solar corona. Our understanding of these events has greatly evolved in recent years but several open questions, such as the contribution of coronal jets to the solar wind, remain. Here we present an overview of the observations and numerical modeling of coronal jets, followed by a brief description of "next-generation" simulations that include an advanced description of the energy transfer in the corona ("thermodynamic MHD"), large spherical computational domains, and the solar wind. These new models will allow us to address some of the open questions., Comment: 12 pages, 4 figures, to appear in ASP proceedings

Measurement of the coronal magnetic field is a crucial ingredient in understanding the nature of solar coronal phenomena at all scales. We employed STEREO/COR1 data obtained during a deep minimum of solar activity in February 2008 (Carrington rotation CR 2066) to retrieve and analyze the three-dimensional (3D) coronal electron density in the range of heights from 1.5 to 4 Rsun using a tomography method. With this, we qualitatively deduced structures of the coronal magnetic field. The 3D electron density analysis is complemented by the 3D STEREO/EUVI emissivity in the 195 A band obtained by tomography for the same CR. A global 3D MHD model of the solar corona was used to relate the reconstructed 3D density and emissivity to open/closed magnetic field structures. We show that the density maximum locations can serve as an indicator of current sheet position, while the locations of the density gradient maximum can be a reliable indicator of coronal hole boundaries. We find that the magnetic field configuration during CR 2066 has a tendency to become radially open at heliocentric distances greater than 2.5 Rsun. We also find that the potential field model with a fixed source surface (PFSS) is inconsistent with the boundaries between the regions with open and closed magnetic field structures. This indicates that the assumption of the potential nature of the coronal global magnetic field is not satisfied even during the deep solar minimum. Results of our 3D density reconstruction will help to constrain solar coronal field models and test the accuracy of the magnetic field approximations for coronal modeling., Comment: Published in "Solar Physics"

Knowledge about the background solar wind plays a crucial role in the framework of space weather forecasting. In-situ measurements of the background solar wind are only available for a few points in the heliosphere where spacecraft are located, therefore we have to rely on heliospheric models to derive the distribution of solar wind parameters in interplanetary space. We test the performance of different solar wind models, namely Magnetohydrodynamic Algorithm outside a Sphere/ENLIL (MAS/ENLIL), Wang-Sheeley-Arge/ENLIL (WSA/ENLIL), and MAS/MAS, by comparing model results with in-situ measurements from spacecraft located at 1 AU distance to the Sun (ACE, Wind). To exclude the influence of interplanetary coronal mass ejections (ICMEs), we chose the year 2007 as a time period with low solar activity for our comparison. We found that the general structure of the background solar wind is well reproduced by all models. The best model results were obtained for the parameter solar wind speed. However, the predicted arrival times of high-speed solar wind streams have typical uncertainties of the order of about one day. Comparison of model runs with synoptic magnetic maps from different observatories revealed that the choice of the synoptic map significantly affects the model performance., Comment: published on-line, Solar Physics, November 2013

A sequence of apparently coupled eruptions was observed on 2010 August 1-2 by SDO and STEREO. The eruptions were closely synchronized with one another, even though some of them occurred at widely separated locations. In an attempt to identify a plausible reason for such synchronization, we study the large-scale structure of the background magnetic configuration. The coronal field was computed from the photospheric magnetic field observed at the appropriate time period by using the potential field source-surface model. We investigate the resulting field structure by analyzing the so-called squashing factor calculated at the photospheric and source-surface boundaries, as well as at different coronal cross-sections. Using this information as a guide, we determine the underlying structural skeleton of the configuration, including separatrix and quasi-separatrix surfaces. Our analysis reveals, in particular, several pseudo-streamers in the regions where the eruptions occurred. Of special interest to us are the magnetic null points and separators associated with the pseudo-streamers. We propose that magnetic reconnection triggered along these separators by the first eruption likely played a key role in establishing the assumed link between the sequential eruptions. The present work substantiates our recent simplified magnetohydrodynamic model of sympathetic eruptions and provides a guide for further deeper study of these phenomena. Several important implications of our results for the S-web model of the slow solar wind are also addressed., Comment: 34 pages, 11 figures, to appear in ApJ in November 2012

Sympathetic eruptions on the Sun have been observed for several decades, but the mechanisms by which one eruption can trigger another one remain poorly understood. We present a 3D MHD simulation that suggests two possible magnetic trigger mechanisms for sympathetic eruptions. We consider a configuration that contains two coronal flux ropes located within a pseudo-streamer and one rope located next to it. A sequence of eruptions is initiated by triggering the eruption of the flux rope next to the streamer. The expansion of the rope leads to two consecutive reconnection events, each of which triggers the eruption of a flux rope by removing a sufficient amount of overlying flux. The simulation qualitatively reproduces important aspects of the global sympathetic event on 2010 August 1 and provides a scenario for so-called twin filament eruptions. The suggested mechanisms are applicable also for sympathetic eruptions occurring in other magnetic configurations., Comment: 6 pages, 4 figures, accepted in ApJ Letters

Models for the origin of the slow solar wind must account for two seemingly contradictory observations: The slow wind has the composition of the closed field corona, implying that it originates from the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind also has large angular width, up to ~ 60{\circ}, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far from the heliospheric current sheet. We then use an MHD code and MDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind and magnetic field for a time period preceding the August 1, 2008 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere, and propose further tests of the model.

In recent work, Antiochos and coworkers argued that the boundary between the open and closed field regions on the Sun can be extremely complex with narrow corridors of open flux connecting seemingly disconnected coronal holes from the main polar holes, and that these corridors may be the sources of the slow solar wind. We examine, in detail, the topology of such magnetic configurations using an analytical source surface model that allows for analysis of the field with arbitrary resolution. Our analysis reveals three important new results: First, a coronal hole boundary can join stably to the separatrix boundary of a parasitic polarity region. Second, a single parasitic polarity region can produce multiple null points in the corona and, more important, separator lines connecting these points. It is known that such topologies are extremely favorable for magnetic reconnection, because they allow this process to occur over the entire length of the separators rather than being confined to a small region around the nulls. Finally, the coronal holes are not connected by an open-field corridor of finite width, but instead are linked by a singular line that coincides with the separatrix footprint of the parasitic polarity. We investigate how the topological features described above evolve in response to the motion of the parasitic polarity region. The implications of our results for the sources of the slow solar wind and for coronal and heliospheric observations are discussed., Comment: 32 pages, 10 figures, accepted to ApJ

Low and mid-latitude coronal holes (CHs) observed on the Sun during the current solar activity minimum (from September 21, 2006, Carrington rotation (CR) 2048, until June 26, 2009 (CR 2084)) were analyzed using {\it SOHO}/EIT and STEREO-A SECCHI EUVI data. From both the observations and Potential Field Source Surface (PFSS) modeling, we find that the area occupied by CHs inside a belt of $\pm 40^\circ$ around the solar equator is larger in the current 2007 solar minimum relative to the similar phase of the previous 1996 solar minimum. The enhanced CH area is related to a recurrent appearance of five persistent CHs, which survived during 7-27 solar rotations. Three of the CHs are of positive magnetic polarity and two are negative. The most long-lived CH was being formed during 2 days and existed for 27 rotations. This CH was associated with fast solar wind at 1 AU of approximately 620$\pm 40$ km s$^{-1}$. The 3D MHD modeling for this time period shows an open field structure above this CH. We conclude that the global magnetic field of the Sun possessed a multi-pole structure during this time period. Calculation of the harmonic power spectrum of the solar magnetic field demonstrates a greater prevalence of multi-pole components over the dipole component in the 2007 solar minimum compared to the 1996 solar minimum. The unusual large separation between the dipole and multi-pole components is due to the very low magnitude of the dipole component, which is three times lower than that in the previous 1996 solar minimum., Comment: 14 pages, 7 figures

A general method for describing magnetic reconnection in arbitrary three-dimensional magnetic configurations is proposed. The method is based on the field-line mapping technique previously used only for the analysis of magnetic structure at a given time. This technique is extended here so as to analyze the evolution of magnetic structure. Such a generalization is made with the help of new dimensionless quantities called "slip-squashing factors". Their large values define the surfaces that border the reconnected or to-be-reconnected magnetic flux tubes for a given period of time during the magnetic evolution. The proposed method is universal, since it assumes only that the time sequence of evolving magnetic field and the tangential boundary flows are known. The application of the method is illustrated for simple examples, one of which was considered previously by Hesse and coworkers in the framework of the general magnetic reconnection theory. The examples help us to compare these two approaches; they reveal also that, just as for magnetic null points, hyperbolic and cusp minimum points of a magnetic field may serve as favorable sites for magnetic reconnection. The new method admits a straightforward numerical implementation and provides a powerful tool for the diagnostics of magnetic reconnection in numerical models of solar-flare-like phenomena in space and laboratory plasmas., Comment: 39 pages, 9 figures, corrected typos, to appear in ApJ, March 2009

Many stormy events in astrophysics occur due to the sudden magnetic energy release. This is possible if a magnetic configuration abruptly changes its topology, an event usually referred to as magnetic reconnection. It is known that pure Ohmic decay is inefficient, occurring during cosmological times (due to the huge characteristic scales $L$). It is recognized that the presence of current sheets speeds up the process, but still insufficiently$^{1,2,3,4,5}$. We show that, in highly compressible and substantially gravitational media, the reconnection is fast enough to account for stormy events. Thus, highly compressible situations offer exiting opportunities in explanations of violent events, although full-scale compressible and gravitational simulations proved to be quite challenging., Comment: 10 pages, 4 figures

A simple model of the coronal magnetic field prior to the CME eruption on May 12 1997 is developed. First, the magnetic field is constructed by superimposing a large-scale background field and a localized bipolar field to model the active region (AR) in the current-free approximation. Second, this potential configuration is quasi-statically sheared by photospheric vortex motions applied to two flux concentrations of the AR. Third, the resulting force-free field is then evolved by canceling the photospheric magnetic flux with the help of an appropriate tangential electric field applied to the central part of the AR. To understand the structure of the modeled configuration, we use the field line mapping technique by generalizing it to spherical geometry. It is demonstrated that the initial potential configuration contains a hyperbolic flux tube (HFT) which is a union of two intersecting quasi-separatrix layers. This HFT provides a partition of the closed magnetic flux between the AR and the global solar magnetic field. The vortex motions applied to the AR interlock the field lines in the coronal volume to form additionally two new HFTs pinched into thin current layers. Reconnection in these current layers helps to redistribute the magnetic flux and current within the AR in the flux-cancellation phase. In this phase, a magnetic flux rope is formed together with a bald patch separatrix surface wrapping around the rope. Other important implications of the identified structural features of the modeled configuration are also discussed., Comment: 25 pages, 11 figures, to appear in ApJ 2008

Numerical simulations of strongly compressible MHD corresponding to a stellar atmosphere with substantial gravity and near force-free magnetic fields show that the current sheet collapses (its width decreasing substantially). As a result, the reconnection rate increases dramatically., Comment: 24 pages 10 figures, typos corrected

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 on 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. (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., Comment: 26 pages, 6 figures

We use direct numerical simulations to study the evolution, or relaxation, of magnetic configurations to an equilibrium state. We use the full single-fluid equations of motion for a magnetized, non-resistive, but viscous fluid; and a Lagrangian approach is used to obtain exact solutions for the magnetic field. As a result, the topology of the magnetic field remains unchanged, which makes it possible to study the case of topological nonequilibrium. We find two cases for which such nonequilibrium appears, indicating that these configurations may develop singular current sheets., Comment: 10 pages, 5 figures