Showing 1–50 of 53 results for author: Török, T

Global Coronal Plasma Diagnostics Based on Multi-slit EUV Spectroscopy

Authors: Lami Chan, Hui Tian, Xianyu Liu, Tibor Török, Xianyong Bai, Yufei Feng, Dipankar Banerjee

Abstract: Full-disk spectroscopic observations of the solar corona are highly desired to forecast solar eruptions and their impact on planets and to uncover the origin of solar wind. In this paper, we introduce a new multi-slit design (5 slits) to obtain extreme ultraviolet (EUV) spectra simultaneously. The selected spectrometer wavelength range (184-197 Å) contains several bright EUV lines that can be used… ▽ More Full-disk spectroscopic observations of the solar corona are highly desired to forecast solar eruptions and their impact on planets and to uncover the origin of solar wind. In this paper, we introduce a new multi-slit design (5 slits) to obtain extreme ultraviolet (EUV) spectra simultaneously. The selected spectrometer wavelength range (184-197 Å) contains several bright EUV lines that can be used for spectral diagnostics. The multi-slit approach offers an unprecedented way to efficiently obtain the global spectral data but the ambiguity from different slits should be resolved. Using a numerical simulation of the global corona, we primarily concentrate on the optimization of the disambiguation process, with the objective of extracting decomposed spectral information of six primary lines. This subsequently facilitates a comprehensive series of plasma diagnostics, including density (Fe XII 195.12/186.89 Å), Doppler velocity (Fe XII 193.51 Å), line width (Fe XII 193.51 Å) and temperature diagnostics (Fe VIII 185.21 Å, Fe X 184.54 Å, Fe XI 188.22 Å, Fe XII 193.51 Å). We find a good agreement between the forward modeling parameters and the inverted results at the initial eruption stage of a coronal mass ejection, indicating the robustness of the decomposition method and its immense potential for global monitoring of the solar corona. △ Less

Submitted 10 June, 2024; v1 submitted 19 April, 2024; originally announced April 2024.

Comments: 14 pages, 5 figures, 2 tables. Accepted on 2024 April 18 for publication in ApJ. Published on 2024 May 30. The name of first author changed once from Linyi Chen (simplified Chinese) to Lami Chan (traditional Chinese) but for the same person

Journal ref: ApJ 967 162 (2024)

Solar Eruptions Triggered by Flux Emergence Below or Near a Coronal Flux Rope

Authors: T. Török, M. G. Linton, J. E. Leake, Z. Mikić, R. Lionello, V. S. Titov, C. Downs

Abstract: 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. He… ▽ More 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. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 18 pages, 8 figures, accepted for publication by The Astrophysical Journal

CORHEL-CME: An Interactive Tool For Modeling Solar Eruptions

Authors: Jon Linker, Tibor Torok, Cooper Downs, Ronald Caplan, Viacheslav Titov, Andres Reyes, Roberto Lionello, Pete Riley

Abstract: Coronal Mass Ejections (CMEs) are immense eruptions of plasma and magnetic fields that are propelled outward from the Sun, sometimes with velocities greater than 2000 km/s. They are responsible for some of the most severe space weather at Earth, including geomagnetic storms and solar energetic particle (SEP) events. We have developed CORHEL-CME, an interactive tool that allows non-expert users to… ▽ More Coronal Mass Ejections (CMEs) are immense eruptions of plasma and magnetic fields that are propelled outward from the Sun, sometimes with velocities greater than 2000 km/s. They are responsible for some of the most severe space weather at Earth, including geomagnetic storms and solar energetic particle (SEP) events. We have developed CORHEL-CME, an interactive tool that allows non-expert users to routinely model multiple CMEs in a realistic coronal and heliospheric environment. The tool features a web-based user interface that allows the user to select a time period of interest, and employs RBSL flux ropes to create stable and unstable pre-eruptive configurations within a background global magnetic field. The properties of these configurations can first be explored in a zero-beta magnetohydrodynamic (MHD) model, followed by complete CME simulations in thermodynamic MHD, with propagation out to 1 AU. We describe design features of the interface and computations, including the innovations required to efficiently compute results on practical timescales with moderate computational resources. CORHEL-CME is now implemented at NASA's Community Coordinated Modeling Center (CCMC) using NASA Amazon Web Services (AWS). It will be available to the public by the time this paper is published. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 10 pages, 6 figures, submitted to Journal of Physics Conference Series

Deflection of Coronal Mass Ejections in Unipolar Ambient Magnetic Fields

Authors: Michal Ben-Nun, Tibor Török, Erika Palmerio, Cooper Downs, Viacheslav S. Titov, Mark G. Linton, Ronald M. Caplan, Roberto Lionello

Abstract: The trajectories of coronal mass ejections (CMEs) are often seen to substantially deviate from a purely radial propagation direction. Such deviations occur predominantly in the corona and have been attributed to "channeling" or deflection of the eruptive flux by asymmetric ambient magnetic fields. Here, we investigate an additional mechanism that does not require any asymmetry of the pre-eruptive… ▽ More The trajectories of coronal mass ejections (CMEs) are often seen to substantially deviate from a purely radial propagation direction. Such deviations occur predominantly in the corona and have been attributed to "channeling" or deflection of the eruptive flux by asymmetric ambient magnetic fields. Here, we investigate an additional mechanism that does not require any asymmetry of the pre-eruptive ambient field. Using magnetohydrodynamic numerical simulations, we show that the trajectory of CMEs through the solar corona can significantly deviate from a radial direction when propagation takes place in a unipolar radial field. We demonstrate that the deviation is most prominent below ~15 solar radii and can be attributed to an "effective IxB force" that arises from the intrusion of a magnetic flux rope with a net axial electric current into a unipolar background field. These results are important for predictions of CME trajectories in the context of space weather forecasts, as well as for reaching a deeper understanding of the fundamental physics underlying CME interactions with the ambient fields in the extended solar corona. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 15 pages, 12 figures

New Observations Needed to Advance Our Understanding of Coronal Mass Ejections

Authors: Erika Palmerio, Benjamin J. Lynch, Christina O. Lee, Lan K. Jian, Teresa Nieves-Chinchilla, Emma E. Davies, Brian E. Wood, Noé Lugaz, Réka M. Winslow, Tibor Török, Nada Al-Haddad, Florian Regnault, Meng Jin, Camilla Scolini, Fernando Carcaboso, Charles J. Farrugia, Vincent E. Ledvina, Cooper Downs, Christina Kay, Sanchita Pal, Tarik M. Salman, Robert C. Allen

Abstract: Coronal mass ejections (CMEs) are large eruptions from the Sun that propagate through the heliosphere after launch. Observational studies of these transient phenomena are usually based on 2D images of the Sun, corona, and heliosphere (remote-sensing data), as well as magnetic field, plasma, and particle samples along a 1D spacecraft trajectory (in-situ data). Given the large scales involved and th… ▽ More Coronal mass ejections (CMEs) are large eruptions from the Sun that propagate through the heliosphere after launch. Observational studies of these transient phenomena are usually based on 2D images of the Sun, corona, and heliosphere (remote-sensing data), as well as magnetic field, plasma, and particle samples along a 1D spacecraft trajectory (in-situ data). Given the large scales involved and the 3D nature of CMEs, such measurements are generally insufficient to build a comprehensive picture, especially in terms of local variations and overall geometry of the whole structure. This White Paper aims to address this issue by identifying the data sets and observational priorities that are needed to effectively advance our current understanding of the structure and evolution of CMEs, in both the remote-sensing and in-situ regimes. It also provides an outlook of possible missions and instruments that may yield significant improvements into the subject. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: White Paper submitted to the Heliophysics 2024-2033 Decadal Survey, 9 pages, 4 figures

Journal ref: Bulletin of the AAS, 55(3), 307, 2023

Work-Life Balance Starts with Proper Deadlines and Exemplary Agencies

Authors: Noé Lugaz, Réka M. Winslow, Nada Al-Haddad, Christina O. Lee, Sarah K. Vines, Katharine Reeves, Amir Caspi, Daniel Seaton, Cooper Downs, Lindsay Glesener, Angelos Vourlidas, Camilla Scolini, Tibor Török, Robert Allen, Erika Palmerio

Abstract: Diversity, equity and inclusion (DEI) programs can only be implemented successfully if proper work-life balance is possible in Heliophysics (and in STEM field in general). One of the core issues stems from the culture of "work-above-life" associated with mission concepts, development, and implementation but also the expectations that seem to originate from numerous announcements from NASA (and oth… ▽ More Diversity, equity and inclusion (DEI) programs can only be implemented successfully if proper work-life balance is possible in Heliophysics (and in STEM field in general). One of the core issues stems from the culture of "work-above-life" associated with mission concepts, development, and implementation but also the expectations that seem to originate from numerous announcements from NASA (and other agencies). The benefits of work-life balance are well documented; however, the entire system surrounding research in Heliophysics hinders or discourages proper work-life balance. For example, there does not seem to be attention paid by NASA Headquarters (HQ) on the timing of their announcements regarding how it will be perceived by researchers, and how the timing may promote a culture where work trumps personal life. The same is true for remarks by NASA HQ program officers during panels or informal discussions, where seemingly innocuous comments may give a perception that work is expected after "normal" work hours. In addition, we are calling for work-life balance plans and implementation to be one of the criteria used for down-selection and confirmation of missions (Key Decision Points: KDP-B, KDP-C). △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 6 pages

Journal ref: Bulletin of the AAS, Vol. 55, Issue 3, Whitepaper #250 (6pp); 2023 July 31

Connecting Solar and Stellar Flares/CMEs: Expanding Heliophysics to Encompass Exoplanetary Space Weather

Authors: B. J. Lynch, B. E. Wood, M. Jin, T. Török, X. Sun, E. Palmerio, R. A. Osten, A. A. Vidotto, O. Cohen, J. D. Alvarado-Gómez, J. J. Drake, V. S. Airapetian, Y. Notsu, A. Veronig, K. Namekata, R. M. Winslow, L. K. Jian, A. Vourlidas, N. Lugaz, N. Al-Haddad, W. B. Manchester, C. Scolini, C. J. Farrugia, E. E. Davies, T. Nieves-Chinchilla , et al. (3 additional authors not shown)

Abstract: The aim of this white paper is to briefly summarize some of the outstanding gaps in the observations and modeling of stellar flares, CMEs, and exoplanetary space weather, and to discuss how the theoretical and computational tools and methods that have been developed in heliophysics can play a critical role in meeting these challenges. The maturity of data-inspired and data-constrained modeling of… ▽ More The aim of this white paper is to briefly summarize some of the outstanding gaps in the observations and modeling of stellar flares, CMEs, and exoplanetary space weather, and to discuss how the theoretical and computational tools and methods that have been developed in heliophysics can play a critical role in meeting these challenges. The maturity of data-inspired and data-constrained modeling of the Sun-to-Earth space weather chain provides a natural starting point for the development of new, multidisciplinary research and applications to other stars and their exoplanetary systems. Here we present recommendations for future solar CME research to further advance stellar flare and CME studies. These recommendations will require institutional and funding agency support for both fundamental research (e.g. theoretical considerations and idealized eruptive flare/CME numerical modeling) and applied research (e.g. data inspired/constrained modeling and estimating exoplanetary space weather impacts). In short, we recommend continued and expanded support for: (1.) Theoretical and numerical studies of CME initiation and low coronal evolution, including confinement of "failed" eruptions; (2.) Systematic analyses of Sun-as-a-star observations to develop and improve stellar CME detection techniques and alternatives; (3.) Improvements in data-inspired and data-constrained MHD modeling of solar CMEs and their application to stellar systems; and (4.) Encouraging comprehensive solar--stellar research collaborations and conferences through new interdisciplinary and multi-agency/division funding mechanisms. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: 9 pages, 5 figures, white paper submitted to the Heliophysics 2024--2033 Decadal Survey

Journal ref: Bulletin of the American Astronomical Society, 55(3), 254, 2023

A Magnetogram-matching Method for Energizing Magnetic Flux Ropes Toward Eruption

Authors: Viacheslav S. Titov, Cooper Downs, Tibor Török, Jon A. Linker

Abstract: We propose a new ``helicity-pumping'' method for energizing coronal equilibria that contain a magnetic flux rope (MFR) toward an eruption. We achieve this in a sequence of magnetohydrodynamics relaxations of small line-tied pulses of magnetic helicity, each of which is simulated by a suitable rescaling of the current-carrying part of the field. The whole procedure is ``magnetogram-matching'' becau… ▽ More We propose a new ``helicity-pumping'' method for energizing coronal equilibria that contain a magnetic flux rope (MFR) toward an eruption. We achieve this in a sequence of magnetohydrodynamics relaxations of small line-tied pulses of magnetic helicity, each of which is simulated by a suitable rescaling of the current-carrying part of the field. The whole procedure is ``magnetogram-matching'' because it involves no changes to the normal component of the field at the photospheric boundary. The method is illustrated by applying it to an observed force-free configuration whose MFR is modeled with our regularized Biot--Savart law method. We find that, in spite of the bipolar character of the external field, the MFR eruption is sustained by two reconnection processes. The first, which we refer to as breakthrough reconnection, is analogous to breakout reconnection in quadrupolar configurations. It occurs at a quasi-separator inside a current layer that wraps around the erupting MFR and is caused by the photospheric line-tying effect. The second process is the classical flare reconnection, which develops at the second quasi-separator inside a vertical current layer that is formed below the erupting MFR. Both reconnection processes work in tandem with the magnetic forces of the unstable MFR to propel it through the overlying ambient field, and their interplay may also be relevant for the thermal processes occurring in the plasma of solar flares. The considered example suggests that our method will be beneficial for both the modeling of observed eruptive events and theoretical studies of eruptions in idealized magnetic configurations. △ Less

Submitted 18 August, 2022; v1 submitted 8 May, 2022; originally announced May 2022.

Comments: 9 pages, 3 figures, new comments and references added, relative magnetic helicity calculated, fig. 2 modified, typos corrected, accepted to ApJ

A Model of Homologous Confined and Ejective Eruptions Involving Kink Instability and Flux Cancellation

Authors: Alshaimaa Hassanin, Bernhard Kliem, Norbert Seehafer, Tibor Török

Abstract: In this study, we model a sequence of a confined and a full eruption, employing the relaxed end state of the confined eruption of a kink-unstable flux rope as the initial condition for the ejective one. The full eruption, a model of a coronal mass ejection, develops as a result of converging motions imposed at the photospheric boundary, which drive flux cancellation. In this process, parts of the… ▽ More In this study, we model a sequence of a confined and a full eruption, employing the relaxed end state of the confined eruption of a kink-unstable flux rope as the initial condition for the ejective one. The full eruption, a model of a coronal mass ejection, develops as a result of converging motions imposed at the photospheric boundary, which drive flux cancellation. In this process, parts of the positive and negative external flux converge toward the polarity inversion line, reconnect, and cancel each other. Flux of the same amount as the canceled flux transfers to a flux rope, increasing the free magnetic energy of the coronal field. With sustained flux cancellation and the associated progressive weakening of the magnetic tension of the overlying flux, we find that a flux reduction of $\approx\!11\%$ initiates the torus instability of the flux rope, which leads to a full eruption. These results demonstrate that a homologous full eruption, following a confined one, can be driven by flux cancellation. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 7 pages, 5 figures

Journal ref: The Astrophysical Journal Letters, 929:L23 (7pp), 2022 April 20

Torus-Stable Zone Above Starspots

Authors: Xudong Sun, Tibor Török, Marc L. DeRosa

Abstract: Whilst intense solar flares are almost always accompanied by a coronal mass ejection (CME), reports on stellar CMEs are rare, despite the frequent detection of stellar 'super flares'. The torus instability of magnetic flux ropes is believed to be one of the main driving mechanisms of solar CMEs. Suppression of the torus instability, due to a confining background coronal magnetic field that decreas… ▽ More Whilst intense solar flares are almost always accompanied by a coronal mass ejection (CME), reports on stellar CMEs are rare, despite the frequent detection of stellar 'super flares'. The torus instability of magnetic flux ropes is believed to be one of the main driving mechanisms of solar CMEs. Suppression of the torus instability, due to a confining background coronal magnetic field that decreases sufficiently slowly with height, may contribute to the lack of stellar CME detection. Here we use the solar magnetic field as a template to estimate the vertical extent of this 'torus-stable zone' (TSZ) above a stellar active region. For an idealised potential field model comprising the fields of a local bipole (mimicking a pair of starspots) and a global dipole, we show that the upper bound of the TSZ increases with the bipole size, the dipole strength, and the source surface radius where the coronal field becomes radial. The boundaries of the TSZ depend on the interplay between the spots' and the dipole's magnetic fields, which provide the local- and global-scale confinement, respectively. They range from about half the bipole size to a significant fraction of the stellar radius. For smaller spots and an intermediate dipole field, a secondary TSZ arises at a higher altitude, which may increase the likelihood of 'failed eruptions'. Our results suggest that the low apparent CME occurrence rate on cool stars is, at least partially, due to the presence of extended TSZs. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 12 pages, 10 figures. Data and 3 interactive figures available via Zenodo. Accepted to MNRAS

Optimization of Magnetic Flux Ropes Modeled with the RBSL method

Authors: V. S. Titov, C. Downs, T. Török, J. A. Linker, R. M. Caplan, R. Lionello

Abstract: The so-called regularized Biot-Savart laws (RBSLs) provide an efficient and flexible method for modeling pre-eruptive magnetic configurations of coronal mass ejections (CMEs) whose characteristics are constrained by observational images and magnetic-field data. This method allows one to calculate the field of magnetic flux ropes (MFRs) with small circular cross-sections and an arbitrary axis shape… ▽ More The so-called regularized Biot-Savart laws (RBSLs) provide an efficient and flexible method for modeling pre-eruptive magnetic configurations of coronal mass ejections (CMEs) whose characteristics are constrained by observational images and magnetic-field data. This method allows one to calculate the field of magnetic flux ropes (MFRs) with small circular cross-sections and an arbitrary axis shape. The field of the whole configuration is constructed as a superposition of (1) such a flux-rope field and (2) an ambient potential field derived, for example, from an observed magnetogram. The RBSL kernels are determined from the requirement that the MFR field for a straight cylinder must be exactly force-free. For a curved MFR, however, the magnetic forces are generally unbalanced over the whole path of the MFR. To minimize these forces, we apply a modified Gauss-Newton method to find optimal MFR parameters. This is done by iteratively adjusting the MFR axis path and axial current. We then try to relax the resulting optimized configuration in a subsequent line-tied zero-beta magnetohydrodynamic simulation toward a force-free equilibrium. By considering two models of the sigmoidal pre-eruption configuration for the 2009 February 13 CME, we demonstrate how this approach works and what it is capable of. We show, in particular, that the building blocks of the core magnetic structure described by these models match to morphological features typically observed in such type of configurations. Our method will be useful for both the modeling of particular eruptive events and theoretical studies of idealized pre-eruptive MFR configurations. △ Less

Submitted 5 June, 2021; originally announced June 2021.

Comments: 21 pages, 11 figures, accepted to ApJS

Energetic Proton Propagation and Acceleration Simulated for the Bastille Day Event of July 14, 2000

Authors: Matthew A. Young, Nathan A. Schwadron, Matthew Gorby, Jon Linker, Ronald M. Caplan, Cooper Downs, Tibor Török, Pete Riley, Roberto Lionello, Viacheslav Titov, Richard A. Mewaldt, Christina M. S. Cohen

Abstract: This work presents results from simulations of the 14 July 2000 ("Bastille Day") solar proton event. We used the Energetic Particle Radiation Environment Model (EPREM) and the CORona-HELiosphere (CORHEL) software suite within the SPE Threat Assessment Tool (STAT) framework to model proton acceleration to GeV energies due to the passage of a CME through the low solar corona, and compared the model… ▽ More This work presents results from simulations of the 14 July 2000 ("Bastille Day") solar proton event. We used the Energetic Particle Radiation Environment Model (EPREM) and the CORona-HELiosphere (CORHEL) software suite within the SPE Threat Assessment Tool (STAT) framework to model proton acceleration to GeV energies due to the passage of a CME through the low solar corona, and compared the model results to GOES-08 observations. The coupled simulation models particle acceleration from 1 to 20 $R_\odot$, after which it models only particle transport. The simulation roughly reproduces the peak event fluxes, and timing and spatial location of the energetic particle event. While peak fluxes and overall variation within the first few hours of the simulation agree well with observations, the modeled CME moves beyond the inner simulation boundary after several hours. The model therefore accurately describes the acceleration processes in the low corona and resolves the sites of most rapid acceleration close to the Sun. Plots of integral flux envelopes from multiple simulated observers near Earth further improve the comparison to observations and increase potential for predicting solar particle events. Broken-power-law fits to fluence spectra agree with diffusive acceleration theory over the low energy range. Over the high energy range, they demonstrate the variability in acceleration rate and mirror the inter-event variability observed solar-cycle 23 GLEs. We discuss ways to improve STAT predictions, including using corrected GOES energy bins and computing fits to the seed spectrum. This paper demonstrates a predictive tool for simulating low-coronal SEP acceleration. △ Less

Submitted 25 January, 2021; v1 submitted 16 December, 2020; originally announced December 2020.

Comments: 18 pages; 9 figures

Decoding the Pre-Eruptive Magnetic Field Configurations of Coronal Mass Ejections

Authors: S. Patsourakos, A. Vourlidas, T. Török, B. Kliem, S. K. Antiochos, V. Archontis, G. Aulanier, X. Cheng, G. Chintzoglou, M. K. Georgoulis, L. M. Green, J. E. Leake, R. Moore, A. Nindos, P. Syntelis, S. L. Yardley, V. Yurchyshyn, J. Zhang

Abstract: A clear understanding of the nature of the pre-eruptive magnetic field configurations of Coronal Mass Ejections (CMEs) is required for understanding and eventually predicting solar eruptions. Only two, but seemingly disparate, magnetic configurations are considered viable; namely, sheared magnetic arcades (SMA) and magnetic flux ropes (MFR). They can form via three physical mechanisms (flux emerge… ▽ More A clear understanding of the nature of the pre-eruptive magnetic field configurations of Coronal Mass Ejections (CMEs) is required for understanding and eventually predicting solar eruptions. Only two, but seemingly disparate, magnetic configurations are considered viable; namely, sheared magnetic arcades (SMA) and magnetic flux ropes (MFR). They can form via three physical mechanisms (flux emergence, flux cancellation, helicity condensation) . Whether the CME culprit is an SMA or an MFR, however, has been strongly debated for thirty years. We formed an International Space Science Institute (ISSI) team to address and resolve this issue and report the outcome here. We review the status of the field across modeling and observations, identify the open and closed issues, compile lists of SMA and MFR observables to be tested against observations and outline research activities to close the gaps in our current understanding. We propose that the combination of multi-viewpoint multi-thermal coronal observations and multi-height vector magnetic field measurements is the optimal approach for resolving the issue conclusively. We demonstrate the approach using MHD simulations and synthetic coronal images. Our key conclusion is that the differentiation of pre-eruptive configurations in terms of SMAs and MFRs seems artificial. Both observations and modeling can be made consistent if the pre-eruptive configuration exists in a hybrid state that is continuously evolving from an SMA to an MFR. Thus, the 'dominant' nature of a given configuration will largely depend on its evolutionary stage (SMA-like early-on, MFR-like near the eruption). △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: Space Science Reviews, accepted for publication

Exploring Plasma Heating in the Current Sheet Region in a Three-Dimensional Coronal Mass Ejection Simulation

Authors: Katharine K. Reeves, Tibor Török, Zoran Mikić, Jon Linker, Nicholas A. Murphy

Abstract: We simulate a coronal mass ejection (CME) using a three-dimensional magnetohydrodynamic (MHD) code that includes coronal heating, thermal conduction, and radiative cooling in the energy equation. The magnetic flux distribution at 1 R$_s$ is produced by a localized subsurface dipole superimposed on a global dipole field, mimicking the presence of an active region within the global corona. Transvers… ▽ More We simulate a coronal mass ejection (CME) using a three-dimensional magnetohydrodynamic (MHD) code that includes coronal heating, thermal conduction, and radiative cooling in the energy equation. The magnetic flux distribution at 1 R$_s$ is produced by a localized subsurface dipole superimposed on a global dipole field, mimicking the presence of an active region within the global corona. Transverse electric fields are applied near the polarity inversion line to introduce a transverse magnetic field, followed by the imposition of a converging flow to form and destabilize a flux rope, producing an eruption. We examine the quantities responsible for plasma heating and cooling during the eruption, including thermal conduction, radiation, adiabatic effects, coronal heating, and ohmic heating. We find that ohmic heating is an important contributor to hot temperatures in the current sheet region early in the eruption, but in the late phase adiabatic compression plays an important role in heating the plasma there. Thermal conduction also plays an important role in the transport of thermal energy away from the current sheet region throughout the reconnection process, producing a ``thermal halo'' and widening the region of high temperatures. We simulate emission from solar telescopes for this eruption and find that there is evidence for emission from heated plasma above the flare loops late in the eruption, when the adiabatic heating is the dominant heating term. These results provide an explanation for hot supra-arcade plasma sheets that are often observed in X-rays and extreme ultraviolet wavelengths during the decay phase of large flares. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 22 pages, 12 figures

Coupled MHD-Focused Transport Simulations for Modeling Solar Particle Events

Authors: Jon A. Linker, Ronald M. Caplan, Nathan Schwadron, Matthew Gorby, Cooper Downs, Tibor Torok, Roberto Lionello, Janvier Wijaya

Abstract: We describe the initial version of the Solar Particle Event (SPE) Threat Assessment Tool or STAT. STAT relies on elements of Corona-Heliosphere (CORHEL) and the Earth-Moon-Mars Radiation Environment Module (EMMREM), and allows users to investigate coronal mass ejection (CME) driven SPEs using coupled magnetohydrodynamic (MHD) and focused transport solutions. At the present time STAT focuses on mod… ▽ More We describe the initial version of the Solar Particle Event (SPE) Threat Assessment Tool or STAT. STAT relies on elements of Corona-Heliosphere (CORHEL) and the Earth-Moon-Mars Radiation Environment Module (EMMREM), and allows users to investigate coronal mass ejection (CME) driven SPEs using coupled magnetohydrodynamic (MHD) and focused transport solutions. At the present time STAT focuses on modeling solar energetic particle (SEP) acceleration in and transport from the low corona, where the highest energy SEP events are generated. We illustrate STAT's capabilities with a model of the July 14, 2000 "Bastille Day" event, including innovative diagnostics for understanding the three-dimensional distribution of particle fluxes and their relation to the structure of the underlying CME driver. A preliminary comparison with NOAA GOES measurements is shown. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Comments: 13 pages, 6 figures, accepted for publication in Journal of Physics Conf. Ser. ASTRONUM 2018

GPU Acceleration of an Established Solar MHD Code using OpenACC

Authors: R. M. Caplan, J. A. Linker, Z. Mikić, C. Downs, T. Török, V. S. Titov

Abstract: 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… ▽ More 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. △ Less

Submitted 6 November, 2018; originally announced November 2018.

Comments: 13 pages, 9 figures

Global Non-Potential Magnetic Models of the Solar Corona During the March 2015 Eclipse

Authors: A. R. Yeates, T. Amari, I. Contopoulos, X. Feng, D. H. Mackay, Z. Mikić, T. Wiegelmann, J. Hutton, C. A. Lowder, H. Morgan, G. Petrie, L. A. Rachmeler, L. A. Upton, A. Canou, P. Chopin, C. Downs, M. Druckmüller, J. A. Linker, D. B. Seaton, T. Török

Abstract: 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 a… ▽ More 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. △ Less

Submitted 2 August, 2018; originally announced August 2018.

Comments: 29 pages, 11 figures, accepted for publication in Space Science Reviews

Sequential eruptions triggered by flux emergence - observations and modeling

Authors: Sally Dacie, Tibor Torok, Pascal Demoulin, Mark Linton, Cooper Downs, Lidia van Driel-Gesztelyi, David Long, James Leake

Abstract: We describe and analyze observations by the Solar Dynamics Observatory of the emergence of a small, bipolar active region within an area of unipolar magnetic flux that was surrounded by a circular, quiescent filament. Within only eight hours of the start of the emergence, a partial splitting of the filament and two consecutive coronal mass ejections took place. We argue that all three dynamic even… ▽ More We describe and analyze observations by the Solar Dynamics Observatory of the emergence of a small, bipolar active region within an area of unipolar magnetic flux that was surrounded by a circular, quiescent filament. Within only eight hours of the start of the emergence, a partial splitting of the filament and two consecutive coronal mass ejections took place. We argue that all three dynamic events occurred as a result of particular magnetic-reconnection episodes between the emerging bipole and the pre-existing coronal magnetic field. In order to substantiate our interpretation, we consider three-dimensional magnetohydrodynamic simulations that model the emergence of magnetic flux in the vicinity of a large-scale coronal flux rope. The simulations qualitatively reproduce most of the reconnection episodes suggested by the observations; as well as the filament-splitting, the first eruption, and the formation of sheared/twisted fields that may have played a role in the second eruption. Our results suggest that the position of emerging flux with respect to the background magnetic configuration is a crucial factor for the resulting evolution, while previous results suggest that parameters such as the orientation or the amount of emerging flux are important as well. This poses a challenge for predicting the onset of eruptions that are triggered by flux emergence, and it calls for a detailed survey of the relevant parameter space by means of numerical simulations. △ Less

Submitted 29 June, 2018; originally announced July 2018.

Comments: 17 pages, 8 figures, accepted for publication by the Astrophysical Journal

Sun-to-Earth MHD Simulation of the 14 July 2000 "Bastille Day" Eruption

Authors: Tibor Török, Cooper Downs, Jon A. Linker, Roberto Lionello, Viacheslav S. Titov, Zoran Mikić, Pete Riley, Ron M. Caplan, Janvier Wijaya

Abstract: Solar eruptions are the main driver of space-weather disturbances at the Earth. Extreme events are of particular interest, not only because of the scientific challenges they pose, but also because of their possible societal consequences. Here we present a magnetohydrodynamic (MHD) simulation of the 14 July 2000 Bastille Day eruption, which produced a very strong geomagnetic storm. After constructi… ▽ More Solar eruptions are the main driver of space-weather disturbances at the Earth. Extreme events are of particular interest, not only because of the scientific challenges they pose, but also because of their possible societal consequences. Here we present a magnetohydrodynamic (MHD) simulation of the 14 July 2000 Bastille Day eruption, which produced a very strong geomagnetic storm. After constructing a thermodynamic MHD model of the corona and solar wind, we insert a magnetically stable flux rope along the polarity inversion line of the eruption's source region and initiate the eruption by boundary flows. More than 10^33 ergs of magnetic energy are released in the eruption within a few minutes, driving a flare, an EUV wave, and a coronal mass ejection (CME) that travels in the outer corona at about 1500 km/s, close to the observed speed. We then propagate the CME to Earth, using a heliospheric MHD code. Our simulation thus provides the opportunity to test how well in situ observations of extreme events are matched if the eruption is initiated from a stable magnetic-equilibrium state. We find that the flux-rope center is very similar in character to the observed magnetic cloud, but arrives about 8.5 hours later and about 15 degrees too far to the North, with field strengths that are too weak by a factor of about 1.6. The front of the flux rope is highly distorted, exhibiting localized magnetic-field concentrations as it passes 1 AU. We discuss these properties with regard to the development of space-weather predictions based on MHD simulations of solar eruptions. △ Less

Submitted 17 January, 2018; originally announced January 2018.

Comments: 27 pages, 13 figures, under revision for publication in the Astrophysical Journal

The origin, early evolution and predictability of solar eruptions

Authors: Lucie Green, Tibor Torok, Bojan Vrsnak, Ward Manchester IV, Astrid Veronig

Abstract: Coronal mass ejections (CMEs) were discovered in the early 1970s when space-borne coronagraphs revealed that eruptions of plasma are ejected from the Sun. Today, it is known that the Sun produces eruptive flares, filament eruptions, coronal mass ejections and failed eruptions; all thought to be due to a release of energy stored in the coronal magnetic field during its drastic reconfiguration. This… ▽ More Coronal mass ejections (CMEs) were discovered in the early 1970s when space-borne coronagraphs revealed that eruptions of plasma are ejected from the Sun. Today, it is known that the Sun produces eruptive flares, filament eruptions, coronal mass ejections and failed eruptions; all thought to be due to a release of energy stored in the coronal magnetic field during its drastic reconfiguration. This review discusses the observations and physical mechanisms behind this eruptive activity, with a view to making an assessment of the current capability of forecasting these events for space weather risk and impact mitigation. Whilst a wealth of observations exist, and detailed models have been developed, there still exists a need to draw these approaches together. In particular more realistic models are encouraged in order to asses the full range of complexity of the solar atmosphere and the criteria for which an eruption is formed. From the observational side, a more detailed understanding of the role of photospheric flows and reconnection is needed in order to identify the evolutionary path that ultimately means a magnetic structure will erupt. △ Less

Submitted 14 January, 2018; originally announced January 2018.

Comments: 65 pages, 8 figures, accepted by Space Science Reviews

Regularized Biot-Savart Laws for Modeling Magnetic Flux Ropes

Authors: Viacheslav S. Titov, Cooper Downs, Zoran Mikić, Tibor Török, Jon A. Linker, Ronald M. Caplan

Abstract: Many existing models assume that magnetic flux ropes play a key role in solar flares and coronal mass ejections (CMEs). It is therefore important to develop efficient methods for constructing flux-rope configurations constrained by observed magnetic data and the morphology of the pre-eruptive source region. For this purpose, we have derived and implemented a compact analytical form that represents… ▽ More Many existing models assume that magnetic flux ropes play a key role in solar flares and coronal mass ejections (CMEs). It is therefore important to develop efficient methods for constructing flux-rope configurations constrained by observed magnetic data and the morphology of the pre-eruptive source region. For this purpose, we have derived and implemented a compact analytical form that represents the magnetic field of a thin flux rope with an axis of arbitrary shape and circular cross-sections. This form implies that the flux rope carries axial current $I$ and axial flux $F$, so that the respective magnetic field is the curl of the sum of axial and azimuthal vector potentials proportional to $I$ and $F$, respectively. We expressed the vector potentials in terms of modified Biot-Savart laws whose kernels are regularized at the axis in such a way that, when the axis is straight, these laws define a cylindrical force-free flux rope with a parabolic profile for the axial current density. For the cases we have studied so far, we determined the shape of the rope axis by following the polarity inversion line of the eruptions' source region, using observed magnetograms. The height variation along the axis and other flux-rope parameters are estimated by means of potential field extrapolations. Using this heuristic approach, we were able to construct pre-eruption configurations for the 2009 February 13 and 2011 October 1 CME events. These applications demonstrate the flexibility and efficiency of our new method for energizing pre-eruptive configurations in simulations of CMEs. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: 9 pages with 5 figures

Electric-Current Neutralization, Magnetic Shear, and Eruptive Activity in Solar Active Regions

Authors: Yang Liu, Xudong Sun, Tibor Török, Viacheslav S. Titov, James E. Leake

Abstract: The physical conditions that determine whether or not solar active regions (ARs) produce strong flares and coronal mass ejections (CMEs) are not yet well understood. Here we investigate the association between electric-current neutralization, magnetic shear along polarity inversion lines (PILs), and eruptive activity in four ARs; two emerging and two well-developed ones. We find that the CME-produ… ▽ More The physical conditions that determine whether or not solar active regions (ARs) produce strong flares and coronal mass ejections (CMEs) are not yet well understood. Here we investigate the association between electric-current neutralization, magnetic shear along polarity inversion lines (PILs), and eruptive activity in four ARs; two emerging and two well-developed ones. We find that the CME-producing ARs are characterized by a strongly non-neutralized total current, while the total current in the ARs that did not produce CMEs is almost perfectly neutralized. The difference in the PIL-shear between these two groups is much less pronounced, which suggests that the degree of current-neutralization may serve as a better proxy for assessing the ability of ARs to produce CMEs. △ Less

Submitted 15 August, 2017; originally announced August 2017.

Comments: ApJL in press

2010 August 1-2 sympathetic eruptions: II. Magnetic topology of the MHD background field

Authors: Viacheslav S. Titov, Zoran Mikic, Tibor Torok, Jon A. Linker, Olga Panasenco

Abstract: Using a potential field source surface (PFSS) model, we recently analyzed the global topology of the background coronal magnetic field for a sequence of coronal mass ejections (CMEs) that occurred on 2010 August 1-2. Here we repeat this analysis for the background field reproduced by a magnetohydrodynamic (MHD) model that incorporates plasma thermodynamics. As for the PFSS model, we find that all… ▽ More Using a potential field source surface (PFSS) model, we recently analyzed the global topology of the background coronal magnetic field for a sequence of coronal mass ejections (CMEs) that occurred on 2010 August 1-2. Here we repeat this analysis for the background field reproduced by a magnetohydrodynamic (MHD) model that incorporates plasma thermodynamics. As for the PFSS model, we find that all three CME source regions contain a coronal hole that is separated from neighboring coronal holes by topologically very similar pseudo-streamer structures. However, the two models yield very different results for the size, shape, and flux of the coronal holes. We find that the helmet-streamer cusp line, which corresponds to a source-surface null line in the PFSS model, is structurally unstable and does not form in the MHD model. Our analysis indicates that generally, in MHD configurations, this line rather consists of a multiple-null separator passing along the edge of disconnected flux regions. Some of these regions are transient and may be the origin of so-called streamer blobs. We show that the core topological structure of such blobs is a three-dimensional "plasmoid", consisting of two conjoined flux ropes of opposite handedness, which connect at a spiral null point of the magnetic field. Our analysis reveals that such plasmoids appear also in pseudo-streamers on much smaller scales. These new insights into the coronal magnetic topology provide some intriguing implications for solar energetic particle events and for the properties of the slow solar wind. △ Less

Submitted 24 July, 2017; originally announced July 2017.

Comments: 29 pages, 7 figures

The Contribution of Coronal Jets To The Solar Wind

Authors: R. Lionello, T. TÖrÖk, V. S. Titov, J. E. Leake, Z. MikiĆ, J. A. Linker, M. G. Linton

Abstract: 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… ▽ More 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. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Solar Coronal Jets: Observations, Theory, and Modeling

Authors: N. E. Raouafi, S. Patsourakos, E. Pariat, P. R. Young, A. C. Sterling, A. Savcheva, M. Shimojo, F. Moreno-Insertis, C. R. DeVore, V. Archontis, T. Török, H. Mason, W. Curdt, K. Meyer, K. Dalmasse, Y. Matsui

Abstract: Coronal jets represent important manifestations of ubiquitous solar transients, which may be the source of significant mass and energy input to the upper solar atmosphere and the solar wind. While the energy involved in a jet-like event is smaller than that of "nominal" solar flares and coronal mass ejections (CMEs), jets share many common properties with these phenomena, in particular, the explos… ▽ More Coronal jets represent important manifestations of ubiquitous solar transients, which may be the source of significant mass and energy input to the upper solar atmosphere and the solar wind. While the energy involved in a jet-like event is smaller than that of "nominal" solar flares and coronal mass ejections (CMEs), jets share many common properties with these phenomena, in particular, the explosive magnetically driven dynamics. Studies of jets could, therefore, provide critical insight for understanding the larger, more complex drivers of the solar activity. On the other side of the size-spectrum, the study of jets could also supply important clues on the physics of transients close or at the limit of the current spatial resolution such as spicules. Furthermore, jet phenomena may hint to basic process for heating the corona and accelerating the solar wind; consequently their study gives us the opportunity to attack a broad range of solar-heliospheric problems. △ Less

Submitted 7 July, 2016; originally announced July 2016.

Comments: 53 pages, 24 figures

Journal ref: Space Science Reviews, 04 July 2016

Modeling Jets in the Corona and Solar Wind

Authors: T. Torok, R. Lionello, V. S. Titov, J. E. Leake, Z. Mikic, J. A. Linker, M. G. Linton

Abstract: 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 bri… ▽ More 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. △ Less

Submitted 30 November, 2015; originally announced November 2015.

Comments: 12 pages, 4 figures, to appear in ASP proceedings

The origin of net electric currents in solar active regions

Authors: K. Dalmasse, G. Aulanier, P. Démoulin, B. Kliem, T. Török, E. Pariat

Abstract: There is a recurring question in solar physics about whether or not electric currents are neutralized in active regions (ARs). This question was recently revisited using three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetic flux emergence into the solar atmosphere. Such simulations showed that flux emergence can generate a substantial net current in ARs. Another source… ▽ More There is a recurring question in solar physics about whether or not electric currents are neutralized in active regions (ARs). This question was recently revisited using three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetic flux emergence into the solar atmosphere. Such simulations showed that flux emergence can generate a substantial net current in ARs. Another source of AR currents are photospheric horizontal flows. Our aim is to determine the conditions for the occurrence of net vs. neutralized currents with this second mechanism. Using 3D MHD simulations, we systematically impose line-tied, quasi-static, photospheric twisting and shearing motions to a bipolar potential magnetic field. We find that such flows: (1) produce both {\it direct} and {\it return} currents, (2) induce very weak compression currents - not observed in 2.5D - in the ambient field present in the close vicinity of the current-carrying field, and (3) can generate force-free magnetic fields with a net current. We demonstrate that neutralized currents are in general produced only in the absence of magnetic shear at the photospheric polarity inversion line - a special condition rarely observed. We conclude that, as magnetic flux emergence, photospheric flows can build up net currents in the solar atmosphere, in agreement with recent observations. These results thus provide support for eruption models based on pre-eruption magnetic fields possessing a net coronal current. △ Less

Submitted 17 July, 2015; originally announced July 2015.

Comments: 14 pages and 11 figures (Accepted in The Astrophysical Journal)

Slow Rise and Partial Eruption of a Double-Decker Filament. II. Modeling by a Double Flux Rope Equilibrium

Authors: Bernhard Kliem, Tibor Török, Viacheslav S. Titov, Roberto Lionello, Jon A. Linker, Rui Liu, Chang Liu, Haimin Wang

Abstract: Force-free equilibria containing two vertically arranged magnetic flux ropes of like chirality and current direction are considered as a model for split filaments/prominences and filament-sigmoid systems. Such equilibria are constructed analytically through an extension of the methods developed in Titov & Démoulin (1999) and numerically through an evolutionary sequence including shear flows, flux… ▽ More Force-free equilibria containing two vertically arranged magnetic flux ropes of like chirality and current direction are considered as a model for split filaments/prominences and filament-sigmoid systems. Such equilibria are constructed analytically through an extension of the methods developed in Titov & Démoulin (1999) and numerically through an evolutionary sequence including shear flows, flux emergence, and flux cancellation in the photospheric boundary. It is demonstrated that the analytical equilibria are stable if an external toroidal (shear) field component exceeding a threshold value is included. If this component decreases sufficiently, then both flux ropes turn unstable for conditions typical of solar active regions, with the lower rope typically being unstable first. Either both flux ropes erupt upward, or only the upper rope erupts while the lower rope reconnects with the ambient flux low in the corona and is destroyed. However, for shear field strengths staying somewhat above the threshold value, the configuration also admits evolutions which lead to partial eruptions with only the upper flux rope becoming unstable and the lower one remaining in place. This can be triggered by a transfer of flux and current from the lower to the upper rope, as suggested by the observations of a split filament in Paper~I (Liu et al. 2012). It can also result from tether-cutting reconnection with the ambient flux at the X-type structure between the flux ropes, which similarly influences their stability properties in opposite ways. This is demonstrated for the numerically constructed equilibrium. △ Less

Submitted 8 July, 2014; originally announced July 2014.

Comments: To appear in the Astrophys. J

Coronal magnetic reconnection driven by CME expansion -- the 2011 June 7 event

Authors: L. van Driel-Gesztelyi, D. Baker, T. Torok, E. Pariat, L. M. Green, D. R. Williams, J. Carlyle, G. Valori, P. Demoulin, B. Kliem, D. M. Long, S. A. Matthews, J. -M. Malherbe

Abstract: Coronal mass ejections (CMEs) erupt and expand in a magnetically structured solar corona. Various indirect observational pieces of evidence have shown that the magnetic field of CMEs reconnects with surrounding magnetic fields, forming, e.g., dimming regions distant from the CME source regions. Analyzing Solar Dynamics Observatory (SDO) observations of the eruption from AR 11226 on 2011 June 7, we… ▽ More Coronal mass ejections (CMEs) erupt and expand in a magnetically structured solar corona. Various indirect observational pieces of evidence have shown that the magnetic field of CMEs reconnects with surrounding magnetic fields, forming, e.g., dimming regions distant from the CME source regions. Analyzing Solar Dynamics Observatory (SDO) observations of the eruption from AR 11226 on 2011 June 7, we present the first direct evidence of coronal magnetic reconnection between the fields of two adjacent ARs during a CME. The observations are presented jointly with a data-constrained numerical simulation, demonstrating the formation/intensification of current sheets along a hyperbolic flux tube (HFT) at the interface between the CME and the neighbouring AR 11227. Reconnection resulted in the formation of new magnetic connections between the erupting magnetic structure from AR 11226 and the neighboring active region AR 11227 about 200 Mm from the eruption site. The onset of reconnection first becomes apparent in the SDO/AIA images when filament plasma, originally contained within the erupting flux rope, is re-directed towards remote areas in AR 11227, tracing the change of large-scale magnetic connectivity. The location of the coronal reconnection region becomes bright and directly observable at SDO/AIA wavelengths, owing to the presence of down-flowing cool, dense (10^{10} cm^{-3}) filament plasma in its vicinity. The high-density plasma around the reconnection region is heated to coronal temperatures, presumably by slow-mode shocks and Coulomb collisions. These results provide the first direct observational evidence that CMEs reconnect with surrounding magnetic structures, leading to a large-scale re-configuration of the coronal magnetic field. △ Less

Submitted 12 June, 2014; originally announced June 2014.

Comments: 12 pages, 12 figures

Journal ref: The Astrophysical Journal, Volume 788, Issue 1, 85, 12 pp. (2014)

Catastrophe versus instability for the eruption of a toroidal solar magnetic flux rope

Authors: B. Kliem, J. Lin, T. G. Forbes, E. R. Priest, T. Török

Abstract: The onset of a solar eruption is formulated here as either a magnetic catastrophe or as an instability. Both start with the same equation of force balance governing the underlying equilibria. Using a toroidal flux rope in an external bipolar or quadrupolar field as a model for the current-carrying flux, we demonstrate the occurrence of a fold catastrophe by loss of equilibrium for several represen… ▽ More The onset of a solar eruption is formulated here as either a magnetic catastrophe or as an instability. Both start with the same equation of force balance governing the underlying equilibria. Using a toroidal flux rope in an external bipolar or quadrupolar field as a model for the current-carrying flux, we demonstrate the occurrence of a fold catastrophe by loss of equilibrium for several representative evolutionary sequences in the stable domain of parameter space. We verify that this catastrophe and the torus instability occur at the same point; they are thus equivalent descriptions for the onset condition of solar eruptions. △ Less

Submitted 8 July, 2014; v1 submitted 23 April, 2014; originally announced April 2014.

Comments: V2: update to conform to the published article; new choice for internal inductance of torus; updated Fig. 2; new Figs. 3, 5, and 8

Journal ref: Astrophys. J. 789, 46, 2014

The Evolution of Writhe in Kink-Unstable Flux Ropes and Erupting Filaments

Authors: Tibor Torok, Bernhard Kliem, Mitchell A. Berger, Mark G. Linton, Pascal Demoulin, Lidia van Driel-Gesztelyi

Abstract: The helical kink instability of a twisted magnetic flux tube has been suggested as a trigger mechanism for solar filament eruptions and coronal mass ejections (CMEs). In order to investigate if estimations of the pre-eruptive twist can be obtained from observations of writhe in such events, we quantitatively analyze the conversion of twist into writhe in the course of the instability, using numeri… ▽ More The helical kink instability of a twisted magnetic flux tube has been suggested as a trigger mechanism for solar filament eruptions and coronal mass ejections (CMEs). In order to investigate if estimations of the pre-eruptive twist can be obtained from observations of writhe in such events, we quantitatively analyze the conversion of twist into writhe in the course of the instability, using numerical simulations. We consider the line tied, cylindrically symmetric Gold-Hoyle flux rope model and measure the writhe using the formulae by Berger and Prior which express the quantity as a single integral in space. We find that the amount of twist converted into writhe does not simply scale with the initial flux rope twist, but depends mainly on the growth rates of the instability eigenmodes of higher longitudinal order than the basic mode. The saturation levels of the writhe, as well as the shapes of the kinked flux ropes, are very similar for considerable ranges of initial flux rope twists, which essentially precludes estimations of pre-eruptive twist from measurements of writhe. However, our simulations suggest an upper twist limit of $\sim 6π$ for the majority of filaments prior to their eruption. △ Less

Submitted 6 March, 2014; originally announced March 2014.

Comments: 16 pages, 4 figures, to appear in Plasma Physics and Controlled Fusion

Distribution of Electric Currents in Solar Active Regions

Authors: Tibor Török, James E. Leake, Viacheslav S. Titov, Vasilis Archontis, Zoran Mikić, Mark G. Linton, Kévin Dalmasse, Guillaume Aulanier, Bernhard Kliem

Abstract: There has been a long-lasting debate on the question of whether or not electric currents in solar active regions are neutralized. That is, whether or not the main (or direct) coronal currents connecting the active region polarities are surrounded by shielding (or return) currents of equal total value and opposite direction. Both theory and observations are not yet fully conclusive regarding this q… ▽ More There has been a long-lasting debate on the question of whether or not electric currents in solar active regions are neutralized. That is, whether or not the main (or direct) coronal currents connecting the active region polarities are surrounded by shielding (or return) currents of equal total value and opposite direction. Both theory and observations are not yet fully conclusive regarding this question, and numerical simulations have, surprisingly, barely been used to address it. Here we quantify the evolution of electric currents during the formation of a bipolar active region by considering a three-dimensional magnetohydrodynamic simulation of the emergence of a sub-photospheric, current-neutralized magnetic flux rope into the solar atmosphere. We find that a strong deviation from current neutralization develops simultaneously with the onset of significant flux emergence into the corona, accompanied by the development of substantial magnetic shear along the active region's polarity inversion line. After the region has formed and flux emergence has ceased, the strong magnetic fields in the region's center are connected solely by direct currents, and the total direct current is several times larger than the total return current. These results suggest that active regions, the main sources of coronal mass ejections and flares, are born with substantial net currents, in agreement with recent observations. Furthermore, they support eruption models that employ pre-eruption magnetic fields containing such currents. △ Less

Submitted 13 January, 2014; originally announced January 2014.

Comments: 6 pages, 5 figures, to appear in Astrophysical Journal Letters

Initiation of Coronal Mass Ejections by Sunspot Rotation

Authors: Tibor Török, Manuela Temmer, Gherardo Valori, Astrid Veronig, Lidia van Driel-Gesztelyi, Bojan Vršnak

Abstract: We study a filament eruption, two-ribbon flare, and coronal mass ejection (CME) that occurred in Active Region NOAA 10898 on 6 July 2006. The filament was located South of a strong sunspot that dominated the region. In the evolution leading up to the eruption, and for some time after it, a counter-clockwise rotation of the sunspot of about 30 degrees was observed. We suggest that the rotation trig… ▽ More We study a filament eruption, two-ribbon flare, and coronal mass ejection (CME) that occurred in Active Region NOAA 10898 on 6 July 2006. The filament was located South of a strong sunspot that dominated the region. In the evolution leading up to the eruption, and for some time after it, a counter-clockwise rotation of the sunspot of about 30 degrees was observed. We suggest that the rotation triggered the eruption by progressively expanding the magnetic field above the filament. To test this scenario, we study the effect of twisting the initially potential field overlying a pre-existing flux-rope, using three-dimensional zero-$β$ MHD simulations. We first consider a relatively simple and symmetric system, and then study a more complex and asymmetric magnetic configuration, whose photospheric flux distribution and coronal structure are guided by the observations and a potential field extrapolation. In both cases, we find that the twisting leads to the expansion of the overlying field. As a consequence of the progressively reduced magnetic tension, the flux-rope quasi-statically adapts to the changed environmental field, rising slowly. Once the tension is sufficiently reduced, a distinct second phase of evolution occurs where the flux-rope enters an unstable regime characterized by a strong acceleration. Our simulations thus suggest a new mechanism for the triggering of eruptions in the vicinity of rotating sunspots. △ Less

Submitted 13 January, 2014; originally announced January 2014.

Comments: 28 pages, 10 figures

Journal ref: Solar Physics, Volume 286, Issue 2, pp.453-477 (2013)

Simulations of Emerging Magnetic Flux. I: The Formation of Stable Coronal Flux Ropes

Authors: James E. Leake, Mark G. Linton, Tibor Török

Abstract: We present results from 3D visco-resistive magnetohydrodynamic (MHD) simulations of the emergence of a convection zone magnetic flux tube into a solar atmosphere containing a pre-existing dipole coronal field, which is orientated to minimize reconnection with the emerging field. We observe that the emergence process is capable of producing a coronal flux rope by the transfer of twist from the conv… ▽ More We present results from 3D visco-resistive magnetohydrodynamic (MHD) simulations of the emergence of a convection zone magnetic flux tube into a solar atmosphere containing a pre-existing dipole coronal field, which is orientated to minimize reconnection with the emerging field. We observe that the emergence process is capable of producing a coronal flux rope by the transfer of twist from the convection zone as found in previous simulations. We find that this flux rope is stable, with no evidence of a fast rise, and that its ultimate height in the corona is determined by the strength of the pre-existing dipole field. We also find that although the electric currents in the initial convection zone flux tube are almost perfectly neutralized, the resultant coronal flux rope carries a significant net current. These results suggest that flux tube emergence is capable of creating non-current-neutralized stable flux ropes in the corona, tethered by overlying potential fields, a magnetic configuration that is believed to be the source of coronal mass ejections. △ Less

Submitted 28 August, 2013; originally announced August 2013.

Comments: 14 figures

2010 August 1-2 sympathetic eruptions: I. Magnetic topology of the source-surface background field

Authors: V. S. Titov, Z. Mikic, T. Török, J. A. Linker, O. Panasenco

Abstract: 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… ▽ More 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. △ Less

Submitted 25 September, 2012; originally announced September 2012.

Comments: 34 pages, 11 figures, to appear in ApJ in November 2012

Contracting and Erupting Components of Sigmoidal Active Regions

Authors: Rui Liu, Chang Liu, Tibor Török, Yuming Wang, Haimin Wang

Abstract: It is recently noted that solar eruptions can be associated with the contraction of coronal loops that are not involved in magnetic reconnection processes. In this paper, we investigate five coronal eruptions originating from four sigmoidal active regions, using high-cadence, high-resolution narrowband EUV images obtained by the Solar Dynamic Observatory (SDO}). The magnitudes of the flares associ… ▽ More It is recently noted that solar eruptions can be associated with the contraction of coronal loops that are not involved in magnetic reconnection processes. In this paper, we investigate five coronal eruptions originating from four sigmoidal active regions, using high-cadence, high-resolution narrowband EUV images obtained by the Solar Dynamic Observatory (SDO}). The magnitudes of the flares associated with the eruptions range from the GOES-class B to X. Owing to the high-sensitivity and broad temperature coverage of the Atmospheric Imaging Assembly (AIA) onboard SDO, we are able to identify both the contracting and erupting components of the eruptions: the former is observed in cold AIA channels as the contracting coronal loops overlying the elbows of the sigmoid, and the latter is preferentially observed in warm/hot AIA channels as an expanding bubble originating from the center of the sigmoid. The initiation of eruption always precedes the contraction, and in the energetically mild events (B and C flares), it also precedes the increase in GOES soft X-ray fluxes. In the more energetic events, the eruption is simultaneous with the impulsive phase of the nonthermal hard X-ray emission. These observations confirm the loop contraction as an integrated process in eruptions with partially opened arcades. The consequence of contraction is a new equilibrium with reduced magnetic energy, as the contracting loops never regain their original positions. The contracting process is a direct consequence of flare energy release, as evidenced by the strong correlation of the maximal contracting speed, and strong anti-correlation of the time delay of contraction relative to expansion, with the peak soft X-ray flux. This is also implied by the relationship between contraction and expansion, i.e., their timing and speed. △ Less

Submitted 2 August, 2012; originally announced August 2012.

Comments: Accepted for publication in ApJ

Slow Rise and Partial Eruption of a Double-Decker Filament. I Observations and Interpretation

Authors: Rui Liu, Bernhard Kliem, Tibor Török, Chang Liu, Viacheslav S. Titov, Roberto Lionello, Jon A. Linker, Haimin Wang

Abstract: We study an active-region dextral filament which was composed of two branches separated in height by about 13 Mm. This "double-decker" configuration sustained for days before the upper branch erupted with a GOES-class M1.0 flare on 2010 August 7. Analyzing this evolution, we obtain the following main results. 1) During hours before the eruption, filament threads within the lower branch were observ… ▽ More We study an active-region dextral filament which was composed of two branches separated in height by about 13 Mm. This "double-decker" configuration sustained for days before the upper branch erupted with a GOES-class M1.0 flare on 2010 August 7. Analyzing this evolution, we obtain the following main results. 1) During hours before the eruption, filament threads within the lower branch were observed to intermittently brighten up, lift upward, and then merge with the upper branch. The merging process contributed magnetic flux and current to the upper branch, resulting in its quasi-static ascent. 2) This transfer might serve as the key mechanism for the upper branch to lose equilibrium by reaching the limiting flux that can be stably held down by the overlying field or by reaching the threshold of the torus instability. 3) The erupting branch first straightened from a reverse S shape that followed the polarity inversion line and then writhed into a forward S shape. This shows a transfer of left-handed helicity in a sequence of writhe-twist-writhe. The fact that the initial writhe is converted into the twist of the flux rope excludes the helical kink instability as the trigger process of the eruption, but supports the occurrence of the instability in the main phase, which is indeed indicated by the very strong writhing motion. 4) A hard X-ray sigmoid, likely of coronal origin, formed in the gap between the two original filament branches in the impulsive phase of the associated flare. This supports a model of transient sigmoids forming in the vertical flare current sheet. 5) Left-handed magnetic helicity is inferred for both branches of the dextral filament. 6) Two types of force-free magnetic configurations are compatible with the data, a double flux rope equilibrium and a single flux rope situated above a loop arcade. △ Less

Submitted 14 August, 2012; v1 submitted 7 July, 2012; originally announced July 2012.

Journal ref: 2012, ApJ, 756, 59

A Parametric Study of Erupting Flux Rope Rotation. Modeling the "Cartwheel CME" on 9 April 2008

Authors: B. Kliem, T. Török, W. T. Thompson

Abstract: The rotation of erupting filaments in the solar corona is addressed through a parametric simulation study of unstable, rotating flux ropes in bipolar force-free initial equilibrium. The Lorentz force due to the external shear field component and the relaxation of tension in the twisted field are the major contributors to the rotation in this model, while reconnection with the ambient field is of m… ▽ More The rotation of erupting filaments in the solar corona is addressed through a parametric simulation study of unstable, rotating flux ropes in bipolar force-free initial equilibrium. The Lorentz force due to the external shear field component and the relaxation of tension in the twisted field are the major contributors to the rotation in this model, while reconnection with the ambient field is of minor importance. Both major mechanisms writhe the flux rope axis, converting part of the initial twist helicity, and produce rotation profiles which, to a large part, are very similar in a range of shear-twist combinations. A difference lies in the tendency of twist-driven rotation to saturate at lower heights than shear-driven rotation. For parameters characteristic of the source regions of erupting filaments and coronal mass ejections, the shear field is found to be the dominant origin of rotations in the corona and to be required if the rotation reaches angles of order 90 degrees and higher; it dominates even if the twist exceeds the threshold of the helical kink instability. The contributions by shear and twist to the total rotation can be disentangled in the analysis of observations if the rotation and rise profiles are simultaneously compared with model calculations. The resulting twist estimate allows one to judge whether the helical kink instability occurred. This is demonstrated for the erupting prominence in the "Cartwheel CME" on 9 April 2008, which has shown a rotation of \approx 115 degrees up to a height of 1.5 R_sun above the photosphere. Out of a range of initial equilibria which include strongly kink-unstable (twist Phi=5pi), weakly kink-unstable (Phi=3.5pi), and kink-stable (Phi=2.5pi) configurations, only the evolution of the weakly kink-unstable flux rope matches the observations in their entirety. △ Less

Submitted 14 December, 2011; originally announced December 2011.

Comments: Solar Physics, submitted

3D Reconstruction of a Rotating Erupting Prominence

Authors: W. T. Thompson, B. Kliem, T. Török

Abstract: A bright prominence associated with a coronal mass ejection (CME) was seen erupting from the Sun on 9 April 2008. This prominence was tracked by both the Solar Terrestrial Relations Observatory (STEREO) EUVI and COR1 telescopes, and was seen to rotate about the line of sight as it erupted; therefore, the event has been nicknamed the "Cartwheel CME." The threads of the prominence in the core of the… ▽ More A bright prominence associated with a coronal mass ejection (CME) was seen erupting from the Sun on 9 April 2008. This prominence was tracked by both the Solar Terrestrial Relations Observatory (STEREO) EUVI and COR1 telescopes, and was seen to rotate about the line of sight as it erupted; therefore, the event has been nicknamed the "Cartwheel CME." The threads of the prominence in the core of the CME quite clearly indicate the structure of a weakly to moderately twisted flux rope throughout the field of view, up to heliocentric heights of 4 solar radii. Although the STEREO separation was 48 degrees, it was possible to match some sharp features in the later part of the eruption as seen in the 304 Å line in EUVI and in the Hα-sensitive bandpass of COR1 by both STEREO Ahead and Behind. These features could then be traced out in three-dimensional space, and reprojected into a view in which the eruption is directed towards the observer. The reconstructed view shows that the alignment of the prominence to the vertical axis rotates as it rises up to a leading-edge height of \approx 2.5 solar radii, and then remains approximately constant. The alignment at 2.5 solar radii differs by about 115 degrees from the original filament orientation inferred from Hα and EUV data, and the height profile of the rotation, obtained here for the first time, shows that two thirds of the total rotation is reached within \approx 0.5 solar radii above the photosphere. These features are well reproduced by numerical simulations of an unstable moderately twisted flux rope embedded in external flux with a relatively strong shear field component. △ Less

Submitted 14 December, 2011; originally announced December 2011.

Comments: published in Solar Physics (Online First)

A model for magnetically coupled sympathetic eruptions

Authors: T. Torok, O. Panasenco, V. S. Titov, Z. Mikic, K. K. Reeves, M. Velli, J. A. Linker, G. De Toma

Abstract: 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… ▽ More 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. △ Less

Submitted 9 August, 2011; originally announced August 2011.

Comments: 6 pages, 4 figures, accepted in ApJ Letters

Homologous Flares and Magnetic Field Topology in Active Region NOAA 10501 on 20 November 2003

Authors: R. Chandra, B. Schmieder, C. H. Mandrini, P. Démoulin, E. Pariat, T. Torok, W. Uddin

Abstract: We present and interpret observations of two morphologically homologous flares that occurred in active region (AR) NOAA 10501 on 20 November 2003. Both flares displayed four homologous H-alpha ribbons and were both accompanied by coronal mass ejections (CMEs). The central flare ribbons were located at the site of an emerging bipole in the center of the active region. The negative polarity of this… ▽ More We present and interpret observations of two morphologically homologous flares that occurred in active region (AR) NOAA 10501 on 20 November 2003. Both flares displayed four homologous H-alpha ribbons and were both accompanied by coronal mass ejections (CMEs). The central flare ribbons were located at the site of an emerging bipole in the center of the active region. The negative polarity of this bipole fragmented in two main pieces, one rotating around the positive polarity by ~ 110 deg within 32 hours. We model the coronal magnetic field and compute its topology, using as boundary condition the magnetogram closest in time to each flare. In particular, we calculate the location of quasiseparatrix layers (QSLs) in order to understand the connectivity between the flare ribbons. Though several polarities were present in AR 10501, the global magnetic field topology corresponds to a quadrupolar magnetic field distribution without magnetic null points. For both flares, the photospheric traces of QSLs are similar and match well the locations of the four H-alpha ribbons. This globally unchanged topology and the continuous shearing by the rotating bipole are two key factors responsible for the flare homology. However, our analyses also indicate that different magnetic connectivity domains of the quadrupolar configuration become unstable during each flare, so that magnetic reconnection proceeds differently in both events. △ Less

Submitted 4 November, 2010; originally announced November 2010.

Comments: 24 pages, 10 figures, Solar Physics (accepted)

Reconnection of a kinking flux rope triggering the ejection of a microwave and hard X-ray source. II. Numerical Modeling

Authors: B. Kliem, M. G. Linton, T. Török, M. Karlický

Abstract: Numerical simulations of the helical ($m\!=\!1$) kink instability of an arched, line-tied flux rope demonstrate that the helical deformation enforces reconnection between the legs of the rope if modes with two helical turns are dominant as a result of high initial twist in the range $Φ\gtrsim6π$. Such reconnection is complex, involving also the ambient field. In addition to breaking up the origina… ▽ More Numerical simulations of the helical ($m\!=\!1$) kink instability of an arched, line-tied flux rope demonstrate that the helical deformation enforces reconnection between the legs of the rope if modes with two helical turns are dominant as a result of high initial twist in the range $Φ\gtrsim6π$. Such reconnection is complex, involving also the ambient field. In addition to breaking up the original rope, it can form a new, low-lying, less twisted flux rope. The new flux rope is pushed downward by the reconnection outflow, which typically forces it to break as well by reconnecting with the ambient field. The top part of the original rope, largely rooted in the sources of the ambient flux after the break-up, can fully erupt or be halted at low heights, producing a "failed eruption." The helical current sheet associated with the instability is squeezed between the approaching legs, temporarily forming a double current sheet. The leg-leg reconnection proceeds at a high rate, producing sufficiently strong electric fields that it would be able to accelerate particles. It may also form plasmoids, or plasmoid-like structures, which trap energetic particles and propagate out of the reconnection region up to the top of the erupting flux rope along the helical current sheet. The kinking of a highly twisted flux rope involving leg-leg reconnection can explain key features of an eruptive but partially occulted solar flare on 18 April 2001, which ejected a relatively compact hard X-ray and microwave source and was associated with a fast coronal mass ejection. △ Less

Submitted 13 July, 2010; originally announced July 2010.

Comments: Solar Physics, in press

Testing magnetofrictional extrapolation with the Titov-Démoulin model of solar active regions

Authors: G. Valori, B. Kliem, T. Török, V. S. Titov

Abstract: We examine the nonlinear magnetofrictional extrapolation scheme using the solar active region model by Titov and Démoulin as test field. This model consists of an arched, line-tied current channel held in force-free equilibrium by the potential field of a bipolar flux distribution in the bottom boundary. A modified version, having a parabolic current density profile, is employed here. We find that… ▽ More We examine the nonlinear magnetofrictional extrapolation scheme using the solar active region model by Titov and Démoulin as test field. This model consists of an arched, line-tied current channel held in force-free equilibrium by the potential field of a bipolar flux distribution in the bottom boundary. A modified version, having a parabolic current density profile, is employed here. We find that the equilibrium is reconstructed with very high accuracy in a representative range of parameter space, using only the vector field in the bottom boundary as input. Structural features formed in the interface between the flux rope and the surrounding arcade-"hyperbolic flux tube" and "bald patch separatrix surface"-are reliably reproduced, as are the flux rope twist and the energy and helicity of the configuration. This demonstrates that force-free fields containing these basic structural elements of solar active regions can be obtained by extrapolation. The influence of the chosen initial condition on the accuracy of reconstruction is also addressed, confirming that the initial field that best matches the external potential field of the model quite naturally leads to the best reconstruction. Extrapolating the magnetogram of a Titov-Démoulin equilibrium in the unstable range of parameter space yields a sequence of two opposing evolutionary phases which clearly indicate the unstable nature of the configuration: a partial buildup of the flux rope with rising free energy is followed by destruction of the rope, losing most of the free energy. △ Less

Submitted 3 May, 2010; originally announced May 2010.

Comments: 14 pages, 10 figures

The writhe of helical structures in the solar corona

Authors: Tibor Toeroek, Mitchell A. Berger, Bernhard Kliem

Abstract: Helicity is a fundamental property of magnetic fields, conserved in ideal MHD. In flux rope topology, it consists of twist and writhe helicity. Despite the common occurrence of helical structures in the solar atmosphere, little is known about how their shape relates to the writhe, which fraction of helicity is contained in writhe, and how much helicity is exchanged between twist and writhe when th… ▽ More Helicity is a fundamental property of magnetic fields, conserved in ideal MHD. In flux rope topology, it consists of twist and writhe helicity. Despite the common occurrence of helical structures in the solar atmosphere, little is known about how their shape relates to the writhe, which fraction of helicity is contained in writhe, and how much helicity is exchanged between twist and writhe when they erupt. Here we perform a quantitative investigation of these questions relevant for coronal flux ropes. The decomposition of the writhe of a curve into local and nonlocal components greatly facilitates its computation. We use it to study the relation between writhe and projected S shape of helical curves and to measure writhe and twist in numerical simulations of flux rope instabilities. The results are discussed with regard to filament eruptions and coronal mass ejections (CMEs). We conclude that the writhe is useful in interpreting S shaped coronal structures and in constraining models of eruptions. △ Less

Submitted 22 April, 2010; originally announced April 2010.

Comments: 11 pages, 12 figures, accepted by Astronomy and Astrophysics

Fan-spine topology formation through two-step reconnection driven by twisted flux emergence

Authors: T. Toeroek, G. Aulanier, B. Schmieder, K. K. Reeves, L. Golub

Abstract: We address the formation of 3D nullpoint topologies in the solar corona by combining Hinode/XRT observations of a small dynamic limb event, which occurred beside a non-erupting prominence cavity, with a 3D zero-beta MHD simulation. To this end, we model the boundary-driven kinematic emergence of a compact, intense, and uniformly twisted flux tube into a potential field arcade that overlies a wea… ▽ More We address the formation of 3D nullpoint topologies in the solar corona by combining Hinode/XRT observations of a small dynamic limb event, which occurred beside a non-erupting prominence cavity, with a 3D zero-beta MHD simulation. To this end, we model the boundary-driven kinematic emergence of a compact, intense, and uniformly twisted flux tube into a potential field arcade that overlies a weakly twisted coronal flux rope. The expansion of the emerging flux in the corona gives rise to the formation of a nullpoint at the interface of the emerging and the pre-existing fields. We unveil a two-step reconnection process at the nullpoint that eventually yields the formation of a broad 3D fan-spine configuration above the emerging bipole. The first reconnection involves emerging fields and a set of large-scale arcade field lines. It results in the launch of a torsional MHD wave that propagates along the arcades, and in the formation of a sheared loop system on one side of the emerging flux. The second reconnection occurs between these newly formed loops and remote arcade fields, and yields the formation of a second loop system on the opposite side of the emerging flux. The two loop systems collectively display an anenome pattern that is located below the fan surface. The flux that surrounds the inner spine field line of the nullpoint retains a fraction of the emerged twist, while the remaining twist is evacuated along the reconnected arcades. The nature and timing of the features which occur in the simulation do qualititatively reproduce those observed by XRT in the particular event studied in this paper. Moreover, the two-step reconnection process suggests a new consistent and generic model for the formation of anemone regions in the solar corona. △ Less

Submitted 11 September, 2009; originally announced September 2009.

Comments: Accepted for publication in ApJ, 11 pages and 5 figures

Journal ref: Astrophys.J.704:485-495,2009

Eruption of magnetic flux ropes during flux emergence

Authors: V. Archontis, T. Torok

Abstract: Aims: We investigate the formation of flux ropes in a flux emergence region and their rise into the outer atmosphere of the Sun. Methods: We perform 3D numerical experiments solving the time-dependent and resistive MHD equations. Results: A sub-photospheric twisted flux tube rises from the solar interior and expands into the corona. A flux rope is formed within the expanding field, due to sh… ▽ More Aims: We investigate the formation of flux ropes in a flux emergence region and their rise into the outer atmosphere of the Sun. Methods: We perform 3D numerical experiments solving the time-dependent and resistive MHD equations. Results: A sub-photospheric twisted flux tube rises from the solar interior and expands into the corona. A flux rope is formed within the expanding field, due to shearing and reconnection of field lines at low atmospheric heights. If the tube emerges into a non-magnetized atmosphere, the flux rope rises, but remains confined inside the expanding magnetized volume. On the contrary, if the expanding tube is allowed to reconnect with a preexisting coronal field, the flux rope experiences a full eruption with a rise profile which is in qualitative agreement with erupting filaments and Coronal Mass Ejections. △ Less

Submitted 7 November, 2008; originally announced November 2008.

Observations and modeling of the early acceleration phase of erupting filaments involved in coronal mass ejections

Authors: C. J. Schrijver, C. Elmore, B. Kliem, T. Toeroek, A. M. Title

Abstract: We examine the early phases of two near-limb filament destabilization involved in coronal mass ejections on 16 June and 27 July 2005, using high-resolution, high-cadence observations made with the Transition Region and Coronal Explorer (TRACE), complemented by coronagraphic observations by Mauna Loa and the SOlar and Heliospheric Observatory (SOHO). The filaments' heights above the solar limb in… ▽ More We examine the early phases of two near-limb filament destabilization involved in coronal mass ejections on 16 June and 27 July 2005, using high-resolution, high-cadence observations made with the Transition Region and Coronal Explorer (TRACE), complemented by coronagraphic observations by Mauna Loa and the SOlar and Heliospheric Observatory (SOHO). The filaments' heights above the solar limb in their rapid-acceleration phases are best characterized by a height dependence h(t) ~ t^m with m near, or slightly above, 3 for both events. Such profiles are incompatible with published results for breakout, MHD-instability, and catastrophe models. We show numerical simulations of the torus instability that approximate this height evolution in case a substantial initial velocity perturbation is applied to the developing instability. We argue that the sensitivity of magnetic instabilities to initial and boundary conditions requires higher fidelity modeling of all proposed mechanisms if observations of rise profiles are to be used to differentiate between them. The observations show no significant delays between the motions of the filament and of overlying loops: the filaments seem to move as part of the overall coronal field until several minutes after the onset of the rapid-acceleration phase. △ Less

Submitted 8 October, 2007; originally announced October 2007.

Comments: ApJ (2007, in press)

Numerical simulations of fast and slow coronal mass ejections

Authors: T. Toeroek, B. Kliem

Abstract: Solar coronal mass ejections (CMEs) show a large variety in their kinematic properties. CMEs originating in active regions and accompanied by strong flares are usually faster and accelerated more impulsively than CMEs associated with filament eruptions outside active regions and weak flares. It has been proposed more than two decades ago that there are two separate types of CMEs, fast (impulsive… ▽ More Solar coronal mass ejections (CMEs) show a large variety in their kinematic properties. CMEs originating in active regions and accompanied by strong flares are usually faster and accelerated more impulsively than CMEs associated with filament eruptions outside active regions and weak flares. It has been proposed more than two decades ago that there are two separate types of CMEs, fast (impulsive) CMEs and slow (gradual) CMEs. However, this concept may not be valid, since the large data sets acquired in recent years do not show two distinct peaks in the CME velocity distribution and reveal that both fast and slow CMEs can be accompanied by both weak and strong flares. We present numerical simulations which confirm our earlier analytical result that a flux-rope CME model permits describing fast and slow CMEs in a unified manner. We consider a force-free coronal magnetic flux rope embedded in the potential field of model bipolar and quadrupolar active regions. The eruption is driven by the torus instability which occurs if the field overlying the flux rope decreases sufficiently rapidly with height. The acceleration profile depends on the steepness of this field decrease, corresponding to fast CMEs for rapid decrease, as is typical of active regions, and to slow CMEs for gentle decrease, as is typical of the quiet Sun. Complex (quadrupolar) active regions lead to the fastest CMEs. △ Less

Submitted 11 October, 2007; v1 submitted 15 May, 2007; originally announced May 2007.

Comments: 4 pages, 3 figures, published in Astron. Nachr. 328, 743 (2007). Minor update to conform to published paper, with minor language improvements and citation updates included

Torus instability

Authors: B. Kliem, T. Toeroek

Abstract: The expansion instability of a toroidal current ring in low-beta magnetized plasma is investigated. Qualitative agreement is obtained with experiments on spheromak expansion and with essential properties of solar coronal mass ejections, unifying the two apparently disparate classes of fast and slow coronal mass ejections. The expansion instability of a toroidal current ring in low-beta magnetized plasma is investigated. Qualitative agreement is obtained with experiments on spheromak expansion and with essential properties of solar coronal mass ejections, unifying the two apparently disparate classes of fast and slow coronal mass ejections. △ Less

Submitted 4 August, 2006; v1 submitted 24 May, 2006; originally announced May 2006.

Comments: minor update to conform to printed version

Journal ref: Phys.Rev.Lett. 96 (2006) 255002

Confined and ejective eruptions of kink-unstable flux ropes

Authors: T. Toeroek, B. Kliem

Abstract: The ideal helical kink instability of a force-free coronal magnetic flux rope, anchored in the photosphere, is studied as a model for solar eruptions. Using the flux rope model of Titov & Demoulin (1999} as the initial condition in MHD simulations, both the development of helical shape and the rise profile of a confined (or failed) filament eruption (on 2002 May 27) are reproduced in very good a… ▽ More The ideal helical kink instability of a force-free coronal magnetic flux rope, anchored in the photosphere, is studied as a model for solar eruptions. Using the flux rope model of Titov & Demoulin (1999} as the initial condition in MHD simulations, both the development of helical shape and the rise profile of a confined (or failed) filament eruption (on 2002 May 27) are reproduced in very good agreement with the observations. By modifying the model such that the magnetic field decreases more rapidly with height above the flux rope, a full (or ejective) eruption of the rope is obtained in very good agreement with the developing helical shape and the exponential-to-linear rise profile of a fast coronal mass ejection (CME) (on 2001 May 15). This confirms that the helical kink instability of a twisted magnetic flux rope can be the mechanism of the initiation and the initial driver of solar eruptions. The agreement of the simulations with properties that are characteristic of many eruptions suggests that they are often triggered by the kink instability. The decrease of the overlying field with height is a main factor in deciding whether the instability leads to a confined event or to a CME. △ Less

Submitted 24 August, 2005; v1 submitted 28 July, 2005; originally announced July 2005.

Comments: minor update to conform to printed version; typo in table corrected

Journal ref: Astrophys.J.630:L97-L100,2005