@article{DalmasseAulanierDemoulinetal.2015, author = {Dalmasse, Kevin and Aulanier, Guillaume and Demoulin, Pascal and Kliem, Bernhard and T{\"o}r{\"o}k, Tibor and Pariat, E.}, title = {The origin of net electric currents in solar active regions}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {810}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.1088/0004-637X/810/1/17}, pages = {14}, year = {2015}, abstract = {There is a recurring question in solar physics regarding 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. Other sources of AR currents are photospheric horizontal flows. Our aim is to determine the conditions for the occurrence of net versus 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 direct and 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 that is rarely observed. We conclude that. photospheric flows,. as magnetic flux emergence, 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 that possess a net coronal current.}, language = {en} } @article{ToeroekKliemBergeretal.2014, author = {Toeroek, T. and Kliem, Bernhard and Berger, M. A. and Linton, M. G. and Demoulin, Pascal and van Driel-Gesztelyi, L.}, title = {The evolution of writhe in kink-unstable flux ropes and erupting filaments}, series = {Plasma physics and controlled fusion}, volume = {56}, journal = {Plasma physics and controlled fusion}, number = {6}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0741-3335}, doi = {10.1088/0741-3335/56/6/064012}, pages = {7}, year = {2014}, 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-emptive 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 similar to 6 pi for the majority of filaments prior to their eruption.}, language = {en} } @article{vanDrielGesztelyiBakerToeroeketal.2014, author = {van Driel-Gesztelyi, L. and Baker, Daniel N. and Toeroek, T. and Pariat, E. and Green, L. M. and Williams, D. R. and Carlyle, J. and Valori, G. and Demoulin, Pascal and Kliem, Bernhard and Long, D. M. and Matthews, S. A. and Malherbe, J. -M.}, title = {Coronal magnetic reconnection driven by CME expansion-the 2011 June 7 event}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {788}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.1088/0004-637X/788/1/85}, pages = {12}, year = {2014}, 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 present the first direct evidence of coronal magnetic reconnection between the fields of two adjacent active regions 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 at the interface between the CME and the neighboring 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 redirected toward 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 (1010 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 reconfiguration of the coronal magnetic field.}, language = {en} } @article{ValoriGreenDemoulinetal.2012, author = {Valori, Gherardo and Green, Lucie M. and Demoulin, Pascal and Vargas Dominguez, S. and van Driel-Gesztelyi, L. and Wallace, A. and Baker, Daniel N. and Fuhrmann, Marcel}, title = {Nonlinear force-free extrapolation of emerging flux with a global twist and serpentine fine structures}, series = {Solar physics : a journal for solar and solar-stellar research and the study of solar terrestrial physics}, volume = {278}, journal = {Solar physics : a journal for solar and solar-stellar research and the study of solar terrestrial physics}, number = {1}, publisher = {Springer}, address = {Dordrecht}, issn = {0038-0938}, doi = {10.1007/s11207-011-9865-8}, pages = {73 -- 97}, year = {2012}, abstract = {We study the flux emergence process in NOAA active region 11024, between 29 June and 7 July 2009, by means of multi-wavelength observations and nonlinear force-free extrapolation. The main aim is to extend previous investigations by combining, as much as possible, high spatial resolution observations to test our present understanding of small-scale (undulatory) flux emergence, whilst putting these small-scale events in the context of the global evolution of the active region. The combination of these techniques allows us to follow the whole process, from the first appearance of the bipolar axial field on the east limb, until the buoyancy instability could set in and raise the main body of the twisted flux tube through the photosphere, forming magnetic tongues and signatures of serpentine field, until the simplification of the magnetic structure into a main bipole by the time the active region reaches the west limb. At the crucial time of the main emergence phase high spatial resolution spectropolarimetric measurements of the photospheric field are employed to reconstruct the three-dimensional structure of the nonlinear force-free coronal field, which is then used to test the current understanding of flux emergence processes. In particular, knowledge of the coronal connectivity confirms the identity of the magnetic tongues as seen in their photospheric signatures, and it exemplifies how the twisted flux, which is emerging on small scales in the form of a sea-serpent, is subsequently rearranged by reconnection into the large-scale field of the active region. In this way, the multi-wavelength observations combined with a nonlinear force-free extrapolation provide a coherent picture of the emergence process of small-scale magnetic bipoles, which subsequently reconnect to form a large-scale structure in the corona.}, language = {en} } @article{ValoriDemoulinPariat2012, author = {Valori, Gherardo and Demoulin, Pascal and Pariat, E.}, title = {Comparing values of the relative magnetic helicity in finite volumes}, series = {Solar physics : a journal for solar and solar-stellar research and the study of solar terrestrial physics}, volume = {278}, journal = {Solar physics : a journal for solar and solar-stellar research and the study of solar terrestrial physics}, number = {2}, publisher = {Springer}, address = {Dordrecht}, issn = {0038-0938}, doi = {10.1007/s11207-012-9951-6}, pages = {347 -- 366}, year = {2012}, abstract = {Relative magnetic helicity, as a conserved quantity of ideal magnetohydrodynamics, has been highlighted as an important quantity to study in plasma physics. Due to its nonlocal nature, its estimation is not straightforward in both observational and numerical data. In this study we derive expressions for the practical computation of the gauge-independent relative magnetic helicity in three-dimensional finite domains. The derived expressions are easy to implement and rapid to compute. They are derived in Cartesian coordinates, but can be easily written in other coordinate systems. We apply our method to a numerical model of a force-free equilibrium containing a flux rope, and compare the results with those obtained employing known half-space equations. We find that our method requires a much smaller volume than half-space expressions to derive the full helicity content. We also prove that values of relative magnetic helicity of different magnetic fields can be compared with each other in the same sense as free-energy values can. Therefore, relative magnetic helicity can be meaningfully and directly compared between different datasets, such as those from different active regions, but also within the same dataset at different times. Typical applications of our formulae include the helicity computation in three-dimensional models of the solar atmosphere, e.g., coronal-field reconstructions by force-free extrapolation and discretized magnetic fields of numerical simulations.}, language = {en} } @article{ValoriDemoulinPariatetal.2013, author = {Valori, Gherarod and Demoulin, Pascal and Pariat, E. and Masson, S.}, title = {Accuracy of magnetic energy computations}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {553}, journal = {Astronomy and astrophysics : an international weekly journal}, number = {2}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201220982}, pages = {14}, year = {2013}, abstract = {Context. For magnetically driven events, the magnetic energy of the system is the prime energy reservoir that fuels the dynamical evolution. In the solar context, the free energy (i.e., the energy in excess of the potential field energy) is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. A trustworthy estimation of the magnetic energy is therefore needed in three-dimensional (3D) models of the solar atmosphere, e. g., in coronal fields reconstructions or numerical simulations. Aims. The expression of the energy of a system as the sum of its potential energy and its free energy (Thomson's theorem) is strictly valid when the magnetic field is exactly solenoidal. For numerical realizations on a discrete grid, this property may be only approximately fulfilled. We show that the imperfect solenoidality induces terms in the energy that can lead to misinterpreting the amount of free energy present in a magnetic configuration. Methods. We consider a decomposition of the energy in solenoidal and nonsolenoidal parts which allows the unambiguous estimation of the nonsolenoidal contribution to the energy. We apply this decomposition to six typical cases broadly used in solar physics. We quantify to what extent the Thomson theorem is not satisfied when approximately solenoidal fields are used. Results. The quantified errors on energy vary from negligible to significant errors, depending on the extent of the nonsolenoidal component of the field. We identify the main source of errors and analyze the implications of adding a variable amount of divergence to various solenoidal fields. Finally, we present pathological unphysical situations where the estimated free energy would appear to be negative, as found in some previous works, and we identify the source of this error to be the presence of a finite divergence. Conclusions. We provide a method of quantifying the effect of a finite divergence in numerical fields, together with detailed diagnostics of its sources. We also compare the efficiency of two divergence-cleaning techniques. These results are applicable to a broad range of numerical realizations of magnetic fields.}, language = {en} }