Magnetic flux generated and intensified by the solar dynamo emerges into the solar atmosphere, forming active regions (ARs) including sunspots. Existing theories of flux emergence suggest that the magnetic flux can rise buoyantly through the convection zone but is trapped at the photosphere, while its further rising into the atmosphere resorts to the Parker buoyancy instability. To trigger such an instability, the Lorentz force in the photosphere needs to be as large as the gas pressure gradient to hold up an extra amount of mass against gravity. This naturally results in a strongly non-force-free photosphere, which is indeed shown in typical idealized numerical simulations of flux tube buoyancy from below the photosphere into the corona. Here we conduct a statistical study of the extents of normalized Lorentz forces and torques in the emerging photospheric magnetic field with a substantially large sample of SDO/HMI vector magnetograms. We found that the photospheric field has a rather small Lorentz force and torque on average, and thus is very close to a force-free state, which is not consistent with theories as well as idealized simulations of flux emergence. Furthermore, the small extents of forces and torques seem not to be influenced by the emerging AR's size, the emergence rate, or the non-potentiality of the field. This result puts an important constraint on future development of theories and simulations of flux emergence.

Solar eruptions are the most powerful drivers of space weather. To understand their cause and nature, it is crucial to know how the coronal magnetic field evolves before eruption. Here we study the formation process of a relatively large-scale magnetic flux rope (MFR) in active region NOAA~12371 that erupts with a major flare and coronal mass ejection on 2015 June 21. A data-driven numerical magnetohydrodynamic model is employed to simulate three-dimensional coronal magnetic field evolution of one-day duration before the eruption. Comparison between the observed features and our modeled magnetic field discloses how the pre-eruption MFR forms. Initially, the magnetic field lines were weakly twisted as being simple sheared arcades. Then a long MFR was formed along the polarity inversion line due to the complex photospheric motion, which is mainly shearing rather than twisting. The presence of the MFR is evidenced by a coherent set of magnetic field lines with twist number above unity. Below the MFR a current sheet is shown in the model, suggesting that tether-cutting reconnection plays a key role in the MFR formation. The MFR's flux grows as more and more field lines are twisted due to continuous injection of magnetic helicity by the photospheric motions. Meanwhile, the height of the MFR's axis increases monotonely from its formation. By an analysis of the decay index of its overlying field, we suggest that it is because the MFR runs into the torus instability regime and becomes unstable that finally triggers the eruption.

Magnetic flux ropes (MFRs) are thought to be the central structure of solar eruptions, and their ideal MHD instabilities can trigger the eruption. Here we performed a study of all the MFR configurations that lead to major solar flares, either eruptive or confined, from 2011 to 2017 near the solar disk center. The coronal magnetic field is reconstructed from observed magnetograms, and based on magnetic twist distribution, we identified the MFR, which is defined as a coherent group of magnetic field lines winding an axis with more than one turn. It is found that 90% of the events possess pre-flare MFRs, and their three-dimensional structures are much more complex in details than theoretical MFR models. We further constructed a diagram based on two parameters, the magnetic twist number which controls the kink instability (KI), and the decay index which controls the torus instability (TI). It clearly shows lower limits for TI and KI thresholds, which are $n_{\rm crit} = 1.3$ and $|T_w|_{\rm crit} = 2$, respectively, as all the events above $n_{\rm crit}$ and nearly 90% of the events above $|T_w|_{\rm crit}$ erupted. Furthermore, by such criterion, over 70% of the events can be discriminated between eruptive and confined flares, and KI seems to play a nearly equally important role as TI in discriminating between the two types of flare. There are more than half of events with both parameters below the lower limits, and 29% are eruptive. These events might be triggered by magnetic reconnection rather than MHD instabilities.

Three-dimensional magnetic topology is crucial to understanding the explosive release of magnetic energy in the corona during solar flares. Much attention has been given to the pre-flare magnetic topology to identify candidate sites of magnetic reconnection, yet it is unclear how the magnetic reconnection and its attendant topological changes shape the eruptive structure and how the topology evolves during the eruption. Here we employed a realistic, data-constrained magnetohydrodynamic simulation to study the evolving magnetic topology for an X9.3 eruptive flare that occurred on 2017 September 6. The simulation successfully reproduces the eruptive features and processes in unprecedented detail. The numerical results reveal that the pre-flare corona contains multiple twisted flux systems with different connections, and during the eruption, these twisted fluxes form a coherent flux rope through tether-cutting-like magnetic reconnection below the rope. Topological analysis shows that the rising flux rope is wrapped by a quasi-separatrix layer, which intersects itself below the rope, forming a topological structure known as hyperbolic flux tube, where a current sheet develops, triggering the reconnection. By mapping footpoints of the newly-reconnected field lines, we are able to reproduce both the spatial location and, for the first time, the temporal separation of the observed flare ribbons, as well as the dynamic boundary of the flux rope's feet. Futhermore, the temporal profile of the total reconnection flux is comparable to the soft X-ray light curve. Such a sophisticated characterization of the evolving magnetic topology provides important insight into the eventual understanding and forecast of solar eruptions.

Three-dimensional magnetic topology of solar flare plays a crucial role in understanding its explosive release of magnetic energy in the corona. However, such three-dimensional coronal magnetic field is still elusive in direct observation. Here we realistically simulate the magnetic evolution during the eruptive process of a great flare, using a numerical magnetohydrodynamic model constrained by observed solar vector magnetogram. The numerical results reveal that the pre-flare corona contains multi-set twisted magnetic flux, which forms a coherent rope during the eruption. The rising flux rope is wrapped by a quasi-separatrix layer, which intersects itself below the rope, forming a hyperbolic flux tube and magnetic reconnection is triggered there. By tracing the footprint of the newly-reconnected field lines, we reproduce both the spatial location and its temporal evolution of flare ribbons with an expected accuracy in comparison of observed images. This scenario strongly confirms the three-dimensional version of standard flare model.

As a fundamental magnetic structure in the solar corona, electric current sheets (CSs) can form either prior to or during solar flare, and they are essential for magnetic energy dissipation in the solar corona by enabling magnetic reconnection. However static reconstruction of CS is rare, possibly due to limitation inherent in available coronal field extrapolation codes. Here we present the reconstruction of a large-scale pre-flare CS in solar active region 11967 using an MHD-relaxation model constrained by SDO/HMI vector magnetogram. The CS is found to be associated with a set of peculiar homologous flares that exhibit unique X-shaped ribbons and loops occurring in a quadrupolar magnetic configuration. This is evidenced by that the field lines traced from the CS to the photosphere form an X shape which nearly precisely reproduces the shape of the observed flare ribbons, suggesting that the flare is a product of the dissipation of the CS through reconnection. The CS forms in a hyperbolic flux tube, which is an intersection of two quasi-separatrix layers. The recurrence of the X-shaped flares might be attributed to the repetitive formation and dissipation of the CS, as driven by the photospheric footpoint motions. These results demonstrate the power of data-constrained MHD model in reproducing CS in the corona as well as providing insight into the magnetic mechanism of solar flares.

Magnetic field extrapolation is an important tool to study the three-dimensional (3D) solar coronal magnetic field which is difficult to directly measure. Various analytic models and numerical codes exist but their results often drastically differ. Thus a critical comparison of the modeled magnetic field lines with the observed coronal loops is strongly required to establish the credibility of the model. Here we compare two different non-potential extrapolation codes, a non-linear force-free field code (CESE-MHD-NLFFF) and a non-force-free field (NFFF) code in modeling a solar active region (AR) that has a sigmoidal configuration just before a major flare erupted from the region. A 2D coronal-loop tracing and fitting method is employed to study the 3D misalignment angles between the extrapolated magnetic field lines and the EUV loops as imaged by SDO/AIA. It is found that the CESE-MHD-NLFFF code with preprocessed magnetogram performs the best, outputting a field which matches the coronal loops in the AR core imaged in AIA 94 {\AA} with a misalignment angle of ~10 degree. This suggests that the CESE-MHD-NLFFF code, even without using the information of coronal loops in constraining the magnetic field, performs as good as some coronal-loop forward-fitting models. For the loops as imaged by AIA 171 {\AA} in the outskirts of the AR, all the codes including the potential-field give comparable results of mean misalignment angle (~30 degree). Thus further improvement of the codes is needed for a better reconstruction of the long loops enveloping the core region.