Flux rope, hyperbolic flux tube, and late EUV phases in a non-eruptive circular-ribbon flare

We present a detailed study of a confined circular flare dynamics associated with 3 UV late phases in order to understand more precisely which topological elements are present and how they constrain the dynamics of the flare. We perform a non-linear force free field extrapolation of the confined flare observed with the HMI and AIA instruments onboard SDO. From the 3D magnetic field we compute the squashing factor and we analyse its distribution. Conjointly, we analyse the AIA EUV light curves and images in order to identify the post-flare loops, their temporal and thermal evolution. By combining both analysis we are able to propose a detailed scenario that explains the dynamics of the flare. Our topological analysis shows that in addition to a null-point topology with the fan separatrix, the spine lines and its surrounding Quasi-Separatix Layers halo (typical for a circular flare), a flux rope and its hyperbolic flux tube (HFT) are enclosed below the null. By comparing the magnetic field topology and the EUV post-flare loops we obtain an almost perfect match 1) between the footpoints of the separatrices and the EUV 1600~\AA{} ribbons and 2) between the HFT's field line footpoints and bright spots observed inside the circular ribbons. We showed, for the first time in a confined flare, that magnetic reconnection occured initially at the HFT, below the flux rope. Reconnection at the null point between the flux rope and the overlying field is only initiated in a second phase. In addition, we showed that the EUV late phase observed after the main flare episode are caused by the cooling loops of different length which have all reconnected at the null point during the impulsive phase.Comment: Astronomy & Astrophysics, in pres

Magneto-frictional Modeling Of Coronal Nonlinear Force-free Fields. I. Testing With Analytic Solutions

We report our implementation of the magneto-frictional method in the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The method aims at applications where local adaptive mesh refinement (AMR) is essential to make follow-up dynamical modeling affordable. We quantify its performance in both domain-decomposed uniform grids and block-adaptive AMR computations, using all frequently employed force-free, divergence-free, and other vector comparison metrics. As test cases, we revisit the semi-analytic solution of Low and Lou in both Cartesian and spherical geometries, along with the topologically challenging Titov-Démoulin model. We compare different combinations of spatial and temporal discretizations, and find that the fourth-order central difference with a local Lax-Friedrichs dissipation term in a single-step marching scheme is an optimal combination. The initial condition is provided by the potential field, which is the potential field source surface model in spherical geometry. Various boundary conditions are adopted, ranging from fully prescribed cases where all boundaries are assigned with the semi-analytic models, to solar-like cases where only the magnetic field at the bottom is known. Our results demonstrate that all the metrics compare favorably to previous works in both Cartesian and spherical coordinates. Cases with several AMR levels perform in accordance with their effective resolutions. The magneto-frictional method in MPI-AMRVAC allows us to model a region of interest with high spatial resolution and large field of view simultaneously, as required by observation-constrained extrapolations using vector data provided with modern instruments. The applications of the magneto-frictional method to observations are shown in an accompanying paper

The Magnetic Environment of a Stealth Coronal Mass Ejection

Interest in stealth coronal mass ejections (CMEs) is increasing due to their relatively high occurrence rate and space weather impact. However, typical CME signatures such as extreme-ultraviolet dimmings and post-eruptive arcades are hard to identify and require extensive image processing techniques. These weak observational signatures mean that little is currently understood about the physics of these events. We present an extensive study of the magnetic field configuration in which the stealth CME of 2011 March 3 occurred. Three distinct episodes of flare ribbon formation are observed in the stealth CME source active region (AR). Two occurred prior to the eruption and suggest the occurrence of magnetic reconnection that builds the structure that will become eruptive. The third occurs in a time close to the eruption of a cavity that is observed in STEREO-B 171 Å data; this subsequently becomes part of the propagating CME observed in coronagraph data. We use both local (Cartesian) and global (spherical) models of the coronal magnetic field, which are complemented and verified by the observational analysis. We find evidence of a coronal null point, with field lines computed from its neighborhood connecting the stealth CME source region to two ARs in the northern hemisphere. We conclude that reconnection at the null point aids the eruption of the stealth CME by removing the field that acted to stabilize the preeruptive structure. This stealth CME, despite its weak signatures, has the main characteristics of other CMEs, and its eruption is driven by similar mechanisms

The Influence of Spatial Resolution on Nonlinear Force-Free Modeling

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.Comment: Accepted to ApJ; comments/corrections to this article are welcome via e-mail, even after publicatio

Solar Physics

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 non-local nature, its estimation is not straightforward in both observational and numerical data. In the present 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. Additionally, we 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 the same dataset at different times. Typical applications of our formulae include the helicity computation in threedimensional models of the solar atmosphere, e.g. coronal-field reconstructions by force-free extrapolation and discretized magnetic fields of numerical simulations

Experience of 170 cases of ductal carcinoma in ' (DCIS)

n/

Nonlinear force-free reconstruction of the global solar magnetic field: methodology

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday's equation, give rise to a respective normal field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modelled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition - the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.Comment: 18 pages, 5 figures, Solar Physic

Spectropolarimetric Fluctuations in a Sunspot Chromosphere

The instrumental advances made in this new era of 4-meter class solar telescopes with unmatched spectropolarimetric accuracy and sensitivity, will enable the study of chromospheric magnetic fields and their dynamics with unprecedented detail. In this regard, spectropolarimetric diagnostics can provide invaluable insight into magneto-hydrodynamic (MHD) wave processes. MHD waves and, in particular, Alfv\'enic fluctuations associated to particular wave modes, were recently recognized as important mechanisms not only for the heating of the outer layers of the Sun's atmosphere and the acceleration of the solar wind, but also for the elemental abundance anomaly observed in the corona of the Sun and other Sun-like stars (also known as first ionisation potential; FIP) effect. Here, we take advantage of state-of-the-art and unique spectropolarimetric IBIS observations to investigate the relation between intensity and circular polarisation (CP) fluctuations in a sunspot chromosphere. Our results show a clear link between the intensity and CP fluctuations in a patch which corresponds to a narrow range of magnetic field inclinations. This suggests the presence of Alfv\'enic perturbations in the sunspot.Comment: 15 pages, 5 figures. Accepted for publication in Philosophical Transactions of the Royal Society

A Helicity-Based Method to Infer the CME Magnetic Field Magnitude in Sun and Geospace: Generalization and Extension to Sun-Like and M-Dwarf Stars and Implications for Exoplanet Habitability

Patsourakos et al. (Astrophys. J. 817, 14, 2016) and Patsourakos and Georgoulis (Astron. Astrophys. 595, A121, 2016) introduced a method to infer the axial magnetic field in flux-rope coronal mass ejections (CMEs) in the solar corona and farther away in the interplanetary medium. The method, based on the conservation principle of magnetic helicity, uses the relative magnetic helicity of the solar source region as input estimates, along with the radius and length of the corresponding CME flux rope. The method was initially applied to cylindrical force-free flux ropes, with encouraging results. We hereby extend our framework along two distinct lines. First, we generalize our formalism to several possible flux-rope configurations (linear and nonlinear force-free, non-force-free, spheromak, and torus) to investigate the dependence of the resulting CME axial magnetic field on input parameters and the employed flux-rope configuration. Second, we generalize our framework to both Sun-like and active M-dwarf stars hosting superflares. In a qualitative sense, we find that Earth may not experience severe atmosphere-eroding magnetospheric compression even for eruptive solar superflares with energies ~ 10^4 times higher than those of the largest Geostationary Operational Environmental Satellite (GOES) X-class flares currently observed. In addition, the two recently discovered exoplanets with the highest Earth-similarity index, Kepler 438b and Proxima b, seem to lie in the prohibitive zone of atmospheric erosion due to interplanetary CMEs (ICMEs), except when they possess planetary magnetic fields that are much higher than that of Earth.Comment: http://adsabs.harvard.edu/abs/2017SoPh..292...89

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

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.Comment: Solar Physics, submitte