Small-scale magnetic flux ropes (SFRs) are a type of structures in the solar wind that possess helical magnetic field lines. In a recent report (Chen & Hu 2020), we presented the radial variations of the properties of SFR from 0.29 to 8 au using in situ measurements from the Helios, ACE/Wind, Ulysses, and Voyager spacecraft. With the launch of the Parker Solar Probe (PSP), we extend our previous investigation further into the inner heliosphere. We apply a Grad-Shafranov-based algorithm to identify SFRs during the first two PSP encounters. We find that the number of SFRs detected near the Sun is much less than that at larger radial distances, where magnetohydrodynamic (MHD) turbulence may act as the local source to produce these structures. The prevalence of Alfvenic structures significantly suppresses the detection of SFRs at closer distances. We compare the SFR event list with other event identification methods, yielding a dozen well-matched events. The cross-section maps of two selected events confirm the cylindrical magnetic flux rope configuration. The power-law relation between the SFR magnetic field and heliocentric distances seems to hold down to 0.16 au.

Gravitational wave (GW) data can be used to test the parity symmetry of gravity by investigating the difference between left-hand and right-hand circular polarization modes. In this article, we develop a method to decompose the circular polarizations of GWs produced during the inspiralling stage of compact binaries, with the help of stationary phase approximation. The foremost advantage is that this method is simple, clean, independent of GW waveform, and is applicable to the existing detector network. Applying it to the mock data, we test the parity symmetry of gravity by constraining the velocity birefringence of GWs. If a nearly edge-on binary neutron-stars with observed electromagnetic counterparts at 40 Mpc is detected by the second-generation detector network, one could derive the model-independent test on the parity symmetry in gravity: the lower limit of the energy scale of parity violation can be constrained within $\mathcal{O}(10^4{\rm eV})$.

The slow solar wind is typically characterized as having low Alfv\'enicity. However, Parker Solar Probe (PSP) observed predominately Alfv\'enic slow solar wind during several of its initial encounters. From its first encounter observations, about 55.3\% of the slow solar wind inside 0.25 au is highly Alfv\'enic ($|\sigma_C| > 0.7$) at current solar minimum, which is much higher than the fraction of quiet-Sun-associated highly Alfv\'enic slow wind observed at solar maximum at 1 au. Intervals of slow solar wind with different Alfv\'enicities seem to show similar plasma characteristics and temperature anisotropy distributions. Some low Alfv\'enicity slow wind intervals even show high temperature anisotropies, because the slow wind may experience perpendicular heating as fast wind does when close to the Sun. This signature is confirmed by Wind spacecraft measurements as we track PSP observations to 1 au. Further, with nearly 15 years of Wind measurements, we find that the distributions of plasma characteristics, temperature anisotropy and helium abundance ratio ($N_\alpha/N_p$) are similar in slow winds with different Alfv\'enicities, but the distributions are different from those in the fast solar wind. Highly Alfv\'enic slow solar wind contains both helium-rich ($N_\alpha/N_p\sim0.045$) and helium-poor ($N_\alpha/N_p\sim0.015$) populations, implying it may originate from multiple source regions. These results suggest that highly Alfv\'enic slow solar wind shares similar temperature anisotropy and helium abundance properties with regular slow solar winds, and they thus should have multiple origins.

With recent advances in the field of machine learning, the use of deep neural networks for time series forecasting has become more prevalent. The quasi-periodic nature of the solar cycle makes it a good candidate for applying time series forecasting methods. We employ a combination of WaveNet and LSTM neural networks to forecast the sunspot number using the years 1749 to 2019 and total sunspot area using the years 1874 to 2019 time series data for the upcoming Solar Cycle 25. Three other models involving the use of LSTMs and 1D ConvNets are also compared with our best model. Our analysis shows that the WaveNet and LSTM model is able to better capture the overall trend and learn the inherent long and short term dependencies in time series data. Using this method we forecast 11 years of monthly averaged data for Solar Cycle 25. Our forecasts show that the upcoming Solar Cycle 25 will have a maximum sunspot number around 106 $\pm$ 19.75 and maximum total sunspot area around 1771 $\pm$ 381.17. This indicates that the cycle would be slightly weaker than Solar Cycle 24.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Theoretically, CME kinematics are related to magnetic reconnection processes in the solar corona. However, the current quantitative understanding of this relationship is based on the analysis of only a handful of events. Here we report a statistical study of 60 CME-flare events from August 2010 to December 2013. We investigate kinematic properties of CMEs and magnetic reconnection in the low corona during the early phase of the eruptions, by combining limb observations from STEREO with simultaneous on-disk views from SDO. For a subset of 42 events with reconnection rate evaluated by the magnetic fluxes swept by the flare ribbons on the solar disk observed from SDO, we find a strong correlation between the peak CME acceleration and the peak reconnection rate. Also, the maximum velocities of relatively fast CMEs (> 600 km/s) are positively correlated with the reconnection flux, but no such correlation is found for slow CMEs. A time-lagged correlation analysis suggests that the distribution of the time lag of CME acceleration relative to reconnection rate exhibits three peaks, approximately 10 minutes apart, and on average, acceleration-lead events have smaller reconnection rates. We further compare the CME total mechanical energy with the estimated energy in the current sheet. The comparison suggests that, for small-flare events, reconnection in the current sheet alone is insufficient to fuel CMEs. Results from this study suggest that flare reconnection may dominate the acceleration of fast CMEs, but for events of slow CMEs and weak reconnection, other mechanisms may be more important.

Small-scale magnetic flux ropes (SFRs), in the solar wind, have been studied for decades. Statistical analysis utilizing various in situ spacecraft measurements is the main observational approach which helps investigate the generation and evolution of these small-scale structures. Based on the Grad-Shafranov (GS) reconstruction technique, we use the automated detection algorithm to build the databases of these small-scale structures via various spacecraft measurements at different heliocentric distances. We present the SFR properties including the magnetic field and plasma parameters at different radial distances from the sun near the ecliptic plane. It is found that the event occurrence rate is still in the order of a few hundreds per month, the duration and scale size distributions follow power laws, and the flux rope axis orientations are approximately centered around the local Parker spiral directions. In general, most SFR properties exhibit radial decays. In addition, with various databases established, we derive scaling laws for the changes of average field magnitude, event counts, and SFR scale sizes, with respect to the radial distances, ranging from $\sim$ 0.3 au for Helios to $\sim$ 7 au for the Voyager spacecraft. The implications of our results for comparisons with the relevant theoretical works and for the application to the Parker Solar Probe (PSP) mission are discussed.

The acceleration of electrons in 3C 279 is investigated through analyzing the injected electron energy distribution (EED) in a time-dependent synchrotron self-Compton + external Compton emission model. In this model, it is assumed that relativistic electrons are continuously injected into the emission region, and the injected EED [$Q_e^\prime(\gamma^\prime)$] follows a single power-law form with low- and high-energy cutoffs $\rm \gamma_{min}'$ and $\rm \gamma_{max}'$, respectively, and the spectral index $n$, i.e, $Q_e^\prime(\gamma^\prime)\propto\gamma^{\prime-n}$. This model is applied to 14 quasi-simultaneous spectral energy distributions (SEDs) of 3C 279. The Markov Chain Monte Carlo fitting technique is performed to obtain the best-fitting parameters and the uncertainties on the parameters. The results show that the injected EED is well constrained in each state. The value of $n$ is in the range of 2.5 to 3.8, which is larger than that expected by the classic non-relativistic shock acceleration. However, the large value of $n$ can be explained by the relativistic oblique shock acceleration. The flaring activity seems to be related to an increased acceleration efficiency, reflected in an increased $\gamma'_{\rm min}$ and electron injection power.

In past decades, much progress has been achieved on the origin and evolution of coronal mass ejections (CMEs). In-situ observations of the counterparts of CMEs, especially magnetic clouds (MCs) near the Earth, have provided measurements of the structure and total flux of CME flux ropes. However, it has been difficult to measure these properties in the erupting CME flux rope, in particular in the pre-existing flux rope. In this work, we propose a model to estimate the toroidal flux of the pre-existing flux rope by subtracting the flux contributed by magnetic reconnection during the eruption from the flux measured in the MC. The flux by the reconnection is derived from geometric properties of two-ribbon flares based on a quasi-2D reconnection model. We then apply the model to four CME/flare events and find that the ratio of toroidal flux in the pre-existing flux rope to that of the associated MC lies in the range of 0.40--0.88. It indicates that the toroidal flux of the pre-existing flux rope has an important contribution to that of the CME flux rope and is usually at least as large as the flux arising from the eruption process for the selected events.

The Parker Solar Probe (PSP) observed an interplanetary coronal mass ejection (ICME) event during its first orbit around the sun, among many other events. This event is analyzed by applying a wavelet analysis technique to obtain the reduced magnetic helicity, cross helicity, and residual energy, the first two of which are magnetohydrodynamics (MHD) invariants. Our results show that the ICME, as a large scale magnetic flux rope, possesses high magnetic helicity, very low cross helicity, and highly negative residual energy, thus pointing to a magnetic fluctuation dominated structure. Using the same technique, we also search for small-scale coherent magnetic flux rope structures during the period from 2018/10/22--2018/11/21, which are intrinsic to quasi-2D MHD turbulence in the solar wind. Multiple structures with duration between 8 and 300 minutes are identified from PSP in-situ spacecraft measurements. The location and scales of these structures are characterized by wavelet spectrograms of the normalized reduced magnetic helicity, normalized cross helicity and normalized residual energy. Transport theory suggests that these small-scale magnetic flux ropes may contribute to the acceleration of charged particles through magnetic reconnection processes, and the dissipation of these structures may be important for understanding the coronal heating processes.

The Godbillon-Vey invariant occurs in homology theory, and algebraic topology, when conditions for a co-dimension 1, foliation of a 3D manifold are satisfied. The magnetic Godbillon-Vey helicity invariant in magnetohydrodynamics (MHD) is a higher order helicity invariant that occurs for flows, in which the magnetic helicity density $h_m={\bf A}{\bf\cdot}{\bf B}={\bf A}{\bf\cdot}(\nabla\times{\bf A})=0$, where ${\bf A}$ is the magnetic vector potential and ${\bf B}$ is the magnetic induction. This paper obtains evolution equations for the magnetic Godbillon-Vey field $\boldsymbol{\eta}={\bf A}\times{\bf B}/|{\bf A}|^2$ and the Godbillon-Vey helicity density $h_{gv}=\boldsymbol{\eta}{\bf\cdot}(\nabla\times{\boldsymbol\eta})$ in general MHD flows in which either $h_m=0$ or $h_m\neq 0$. A conservation law for $h_{gv}$ occurs in flows for which $h_m=0$. For $h_m\neq 0$ the evolution equation for $h_{gv}$ contains a source term in which $h_m$ is coupled to $h_{gv}$ via the shear tensor of the background flow. The transport equation for $h_{gv}$ also depends on the electric field potential $\psi$, which is related to the gauge for ${\bf A}$, which takes its simplest form for the advected ${\bf A}$ gauge in which $\psi={\bf A\cdot u}$ where ${\bf u}$ is the fluid velocity. An application of the Godbillon-Vey magnetic helicity to nonlinear force-free magnetic fields used in solar physics is investigated. The possible uses of the Godbillon-Vey helicity in zero helicity flows in ideal fluid mechanics, and in zero helicity Lagrangian kinematics of three-dimensional advection are discussed.

A novel technique is presented for describing and visualizing the local topology of the magnetic field using single-spacecraft data in the solar wind. The approach merges two established techniques: the Grad-Shafranov (GS) reconstruction method, which provides a plausible regional two-dimensional magnetic field surrounding the spacecraft trajectory, and the Partial Variance of Increments (PVI) technique that identifies coherent magnetic structures, such as current sheets. When applied to one month of Wind magnetic field data at 1-minute resolution, we find that the quasi-two-dimensional turbulence emerges as a sea of magnetic islands and current sheets. Statistical analysis confirms that current sheets associated with high values of PVI are mostly located between and within the GS magnetic islands, corresponding to X-points and internal boundaries. The method shows great promise for visualizing and analyzing single-spacecraft data from missions such as Parker Solar Probe and Solar Orbiter, as well as 1 AU Space Weather monitors such as ACE, Wind and IMAP.

Small-scale magnetic flux ropes, in the solar wind, have been studied for decades via the approach of both simulation and observation. Statistical analysis utilizing various in-situ spacecraft measurements is the main observational approach. In this study, we extend the automated detection of small-scale flux ropes based on the Grad-Shafranov reconstruction to the complete dataset of \emph{Ulysses} spacecraft in-situ measurements. We first discuss the temporal variation of the bulk properties of 22,719 flux ropes found through our approach, namely, the average magnetic field and plasma parameters, etc., as functions of the heliographical latitudes and heliocentric radial distances. We then categorize all identified events into three groups based on event distributions in different latitudes separated at 30$^{\circ}$, at different radial distances, and under different solar activities, respectively. By the detailed statistical analysis, we conclude as follows. (1) The properties of flux ropes, such as the duration, scale size, etc., follow the power-law distributions, but with different slope indices, especially for distributions at different radial distances. (2) Also, they are affected by the solar wind speed which has different distributions under the different solar activities, which is manifested as the latitudinal effect. (3) The main difference in flux rope properties between the low and high latitudes is attributed to possible Alfv\'enic structures or waves. (4) Flux ropes with longer duration and larger scale sizes occur more often at larger radial distances. (5) With more strict Wal\'en slope threshold, more events are excluded at higher latitudes. The entire database is published online at \url{http://www.fluxrope.info}.

It is well accepted that a magnetic flux rope (MFR) is a critical component of many coronal mass ejections (CMEs), yet how it evolves toward eruption remains unclear. Here we investigate the continuous evolution of a pre-existing MFR, which is rooted in strong photospheric magnetic fields and electric currents. The evolution of the MFR is observed by the Solar Terrestrial Relations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) from multiple viewpoints. From STEREO's perspective, the MFR starts to rise slowly above the limb five hours before it erupts as a halo CME on 2012 June 14. In SDO observations, conjugate dimmings develop on the disk, simultaneously with the gradual expansion of the MFR, suggesting that the dimmings map the MFR's feet. The evolution comprises a two-stage gradual expansion followed by another stage of rapid acceleration/eruption. Quantitative measurements indicate that magnetic twist of the MFR increases from 1.0 +/- 0.5 to 2.0 +/- 0.5 turns during the five-hour expansion, and further increases to about 4.0 turns per AU when detected as a magnetic cloud at 1 AU two day later. In addition, each stage is preceded by flare(s), implying reconnection is actively involved in the evolution and eruption of the MFR. The implications of these measurements on the CME initiation mechanisms are discussed.

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.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

We have developed a new automated small-scale magnetic flux rope (SSMFR) detection algorithm based on the Grad-Shafranov (GS) reconstruction technique. We have applied this detection algorithm to the Wind spacecraft in-situ measurements during 1996 - 2016, covering two solar cycles, and successfully detected a total number of 74,241 small-scale magnetic flux rope events with duration from 9 to 361 minutes. This large number of small-scale magnetic flux ropes has not been discovered by any other previous studies through this unique approach. We perform statistical analysis of the small-scale magnetic flux rope events based on our newly developed database, and summarize the main findings as follows. (1) The occurrence of small-scale flux ropes has strong solar cycle dependency with a rate of a few hundreds per month on average. (2) The small-scale magnetic flux ropes in the ecliptic plane tend to align along the Parker spiral. (3) In low speed ($<$ 400 km/s) solar wind, the flux ropes tend to have lower proton temperature and higher proton number density, while in high speed ($\ge$ 400 km/s) solar wind, they tend to have higher proton temperature and lower proton number density. (4) Both the duration and scale size distributions of the small-scale magnetic flux ropes obey a power law. (5) The waiting time distribution of small-scale magnetic flux ropes can be fitted by an exponential function (for shorter waiting times) and a power law function (for longer waiting times). (6) The wall-to-wall time distribution obeys double power laws with the break point at 60 minutes (corresponding to the correlation length). (7) The small-scale magnetic flux ropes tend to accumulate near the heliospheric current sheet (HCS). The entire database is available at \url{http://fluxrope.info} and in machine readable format in this article.

Magnetohydrodynamics of the solar corona is simulated numerically. The simulation is initialized with an extrapolated non-force-free magnetic field using the vector magnetogram of the active region (AR) NOAA 12192 obtained on the solar photosphere. Particularly, we focus on the magnetic reconnections occurring close to a magnetic null-point that resulted in appearance of circular chromospheric flare ribbons on October 24, 2014 around 21:21 UT, after peak of an X3.1 flare. The extrapolated field lines show the presence of the three-dimensional (3D) null near one of the polarity inversion lines---where the flare was observed. In the subsequent numerical simulation, we find magnetic reconnections occurring near the null point, where the magnetic field lines from the fan-plane of the 3D null form a X-type configuration with underlying arcade field lines. The footpoints of the dome-shaped field lines, inherent to the 3D null, show high gradients of the squashing factor. We find slipping reconnections at these quasi-separatrix layers, which are co-located with the post-flare circular brightening observed at the chromospheric heights. This demonstrates the viability of the initial non-force-free field along with the dynamics it initiates. Moreover, the initial field and its simulated evolution is found to be devoid of any flux rope, which is in congruence with the confined nature of the flare.

Magnetic reconnections (MRs) for various magnetic field line (MFL) topologies are believed to be the initiators of various solar eruptive events like flares and coronal mass ejections (CMEs). Consequently, important is a thorough understanding and quantification of the MFL topology and their evolution which leads to MRs. Contemporary standard is to extrapolate the coronal MFLs using equilibrium models where the Lorentz force on the coronal plasma is zero everywhere, because either there is no current or the current is parallel to the magnetic field. In tandem, a non-force-free-field (NFFF) extrapolation scheme has evolved and allows for a Lorentz force which is non-zero only at the photosphere but asymptotically vanishes with height. The paper reports magnetohydrodynamic (MHD)- simulations initiated by NFFF extrapolation of the coronal MFLs for the active region NOAA 11158. Interestingly, quasi-separatrix layers (QSLs) which facilitate MRs are detected in the extrapolated MFLs. The AR 11158 is flare producing and, the paper makes an attempt to asses the role of QSLs in the flare onsets.

The solar active region (AR) 12192 was one of the most flare productive region of solar cycle 24, which produced many X-class flares; the most energetic being an X3.1 flare on October 24, 2014 at 21:10 UT. Customarily, such events are believed to be triggered by magnetic reconnection in coronal magnetic fields. Here we use the vector magnetograms from solar photosphere, obtained from Heliospheric Magnetic Imager (HMI) to investigate the magnetic field topology prior to the X3.1 event, and ascertain the conditions that might have caused the flare. To infer the coronal magnetic field, a novel non-force-free field (NFFF) extrapolation technique of the photospheric field is used, which suitably mimics the Lorentz forces present in the photospheric plasma. We also highlight the presence of magnetic null points and quasi-separatrix layers (QSLs) in the magnetic field topology, which are preferred sites for magnetic reconnections and discuss the probable reconnection scenarios.

Magnetic flux rope (MFR) is the core structure of the greatest eruptions, i.e., the coronal mass ejections (CMEs), on the Sun, and magnetic clouds are post-eruption MFRs in interplanetary space. There is a strong debate about whether or not a MFR exists prior to a CME and how the MFR forms/grows through magnetic reconnection during the eruption. Here we report a rare event, in which a magnetic cloud was observed sequentially by four spacecraft near Mercury, Venus, Earth and Mars, respectively. With the aids of a uniform-twist flux rope model and a newly developed method that can recover a shock-compressed structure, we find that the axial magnetic flux and helicity of the magnetic cloud decreased when it propagated outward but the twist increased. Our analysis suggests that the `pancaking' effect and `erosion' effect may jointly cause such variations. The significance of the `pancaking' effect is difficult to be estimated, but the signature of the erosion can be found as the imbalance of the azimuthal flux of the cloud. The latter implies that the magnetic cloud was eroded significantly leaving its inner core exposed to the solar wind at far distance. The increase of the twist together with the presence of the erosion effect suggests that the post-eruption MFR may have a high-twist core enveloped by a less-twisted outer shell. These results pose a great challenge to the current understanding on the solar eruptions as well as the formation and instability of MFRs.

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.

GRB 170817A is the first short gamma-ray burst (GRB) with direct detection of the gravitational-wave radiation and also the spectroscopically identified macronova emission (i.e., AT 2017gfo). The prompt emission of this burst, however, is underluminous in comparison with the other short GRBs with known redshift. In this work, we examine whether GRB 170817A is indeed unique. We firstly show that GRB 130603B/macronova may be the on-axis "analogs" of GRB 170817A/AT 2017gfo, and the extremely dim { but long-lasting} afterglow emission of GRB 170817A may suggest a low number density ($\sim 10^{-5}~{\rm cm^{-3}}$) of its circumburst medium { and a structured outflow}. We then discuss whether GRB 070923, GRB 080121, GRB 090417A, GRB 111005A, and GRB 170817A form a new group of very nearby underluminous GRBs originated from neutron star mergers. If the short events GRB 070923, GRB 080121, and GRB 090417A are indeed at a redshift of $\sim 0.076,~0.046,~0.088$, respectively, their isotropic energies of the prompt emission are $\sim 10^{47}$ erg and thus comparable to the other two events. The non-detection of optical counterparts of GRB 070923, GRB 080121, GRB 090417A, and GRB 111005A, however, strongly suggests that the macronovae from neutron star mergers are significantly diverse in luminosities or, alternatively, there is the other origin channel (for instance, the white dwarf and black hole mergers). We finally suggest that GW170817/GRB 170817A are likely not alone and similar events will be detected by the upgraded/upcoming gravitational-wave detectors and the electromagnetic monitors.

We present unique {and additional} observational evidence for the self-generation of small-scale coherent magnetic flux rope structures in the solar wind. Such structures with durations between 9 and 361 minutes are identified from Wind in-situ spacecraft measurements through the Grad-Shafranov (GS) reconstruction approach. The event occurrence counts are on the order of 3,500 per year on average and have a clear solar cycle dependence. We build a database of small-scale magnetic flux ropes from twenty-year worth of Wind spacecraft data. We show a power-law distribution of the wall-to-wall time corresponding well to the inertial range turbulence, which agrees with relevant observations and numerical simulation results. We also provide the axial current density distribution from the GS-based observational analysis, which yields a non-Gaussian probability density function consistent with numerical simulation results.

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.

This article completes and extends a recent study of the Grad-Shafranov (GS) reconstruction in toroidal geometry, as applied to a two and a half dimensional configurations in space plasmas with rotational symmetry. A further application to the benchmark study of an analytic solution to the toroidal GS equation with added noise shows deviations in the reconstructed geometry of the flux rope configuration, characterized by the orientation of the rotation axis, the major radius, and the impact parameter. On the other hand, the physical properties of the flux rope, including the axial field strength, and the toroidal and poloidal magnetic flux, agree between the numerical and exact GS solutions. We also present a real event study of a magnetic cloud flux rope from \textit{in situ} spacecraft measurements. The devised procedures for toroidal GS reconstruction are successfully executed. Various geometrical and physical parameters are obtained with associated uncertainty estimates. The overall configuration of the flux rope from the GS reconstruction is compared with the corresponding morphological reconstruction based on white-light images. The results show overall consistency, but also discrepancy in that the inclination angle of the flux rope central axis with respect to the ecliptic plane differs by about 20-30 degrees in the plane of the sky. We also compare the results with the original straight-cylinder GS reconstruction and discuss our findings.

We develop an approach of Grad-Shafranov (GS) reconstruction for toroidal structures in space plasmas, based on in-situ spacecraft measurements. The underlying theory is the GS equation that describes two-dimensional magnetohydrostatic equilibrium as widely applied in fusion plasmas. The geometry is such that the arbitrary cross section of the torus has rotational symmetry about the rotation axis $Z$, with a major radius $r_0$. The magnetic field configuration is thus determined by a scalar flux function $\Psi$ and a functional $F$ that is a single-variable function of $\Psi$. The algorithm is implemented through a two-step approach: i) a trial-and-error process by minimizing the residue of the functional $F(\Psi)$ to determine an optimal $Z$ axis orientation, and ii) for the chosen $Z$, a $\chi^2$ minimization process resulting in the range of $r_0$. Benchmark studies of known analytic solutions to the toroidal GS equation with noise additions are presented to illustrate the two-step procedures and to demonstrate the performance of the numerical GS solver, separately. For the cases presented, the errors in $Z$ and $r_0$ are 9$^\circ$ and 22\%, respectively, and the relative percent error in the numerical GS solutions is less than 10\%. We also make public the computer codes for these implementations and benchmark studies.

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.

With SDO observations and a data-constrained MHD model, we identify a confined multi-ribbon flare occurred on 2010 October 25 in solar active region 11117 as a magnetic bald patch (BP) flare with strong evidences. From the photospheric magnetic field observed by SDO/HMI, we find there is indeed magnetic BPs on the PILs which match parts of the flare ribbons. From the 3D coronal magnetic field derived from a MHD relaxation model constrained by the vector magnetograms, we find strikingly good agreement of the BP separatrix surface (BPSS) footpoints with the flare ribbons, and the BPSS itself with the hot flaring loop system. Moreover, the triggering of the BP flare can be attributed to a small flux emergence under the lobe of the BPSS, and the relevant change of the coronal magnetic field through the flare is well reproduced by the pre-flare and post-flare MHD solutions, which match the corresponding pre and post-flare AIA observations, respectively. Our work contributes to the study of non-typical flares that constitute the majority of solar flares but cannot be explained by the standard flare model.

Magnetic flux ropes (MFRs) are one kind of fundamental structures in the solar physics, and involved in various eruption phenomena. Twist, characterizing how the magnetic field lines wind around a main axis, is an intrinsic property of MFRs, closely related to the magnetic free energy and stableness. So far it is unclear how much amount of twist is carried by MFRs in the solar atmosphere and in heliosphere and what role the twist played in the eruptions of MFRs. Contrasting to the solar MFRs, there are lots of in-situ measurements of magnetic clouds (MCs), the large-scale MFRs in interplanetary space, providing some important information of the twist of MFRs. Thus, starting from MCs, we investigate the twist of interplanetary MFRs with the aid of a velocity-modified uniform-twist force-free flux rope model. It is found that most of MCs can be roughly fitted by the model and nearly half of them can be fitted fairly well though the derived twist is probably over-estimated by a factor of 2.5. By applying the model to 115 MCs observed at 1 AU, we find that (1) the twist angles of interplanetary MFRs generally follow a trend of about $0.6\frac{l}{R}$ radians, where $\frac{l}{R}$ is the aspect ratio of a MFR, with a cutoff at about $12\pi$ radians AU$^{-1}$, (2) most of them are significantly larger than $2.5\pi$ radians but well bounded by $2\frac{l}{R}$ radians, (3) strongly twisted magnetic field lines probably limit the expansion and size of MFRs, and (4) the magnetic field lines in the legs wind more tightly than those in the leading part of MFRs. These results not only advance our understanding of the properties and behavior of interplanetary MFRs, but also shed light on the formation and eruption of MFRs in the solar atmosphere. A discussion about the twist and stableness of solar MFRs are therefore given.

The largest geomagnetic storm so far in the solar cycle 24 was produced by a fast coronal mass ejection (CME) originating on 2015 March 15. It was an initially west-oriented CME and expected to only cause a weak geomagnetic disturbance. Why did this CME finally cause such a large geomagnetic storm? We try to find some clues by investigating its propagation from the Sun to 1 AU. First, we reconstruct the CME's kinematic properties in the corona from the SOHO and SDO imaging data with the aid of the graduated cylindrical shell (GCS) model. It is suggested that the CME propagated to the west $\sim$$33^\circ$$\pm$$10^\circ$ away from the Sun-Earth line with a speed of about 817 km s$^{-1}$ before leaving the field of view of the SOHO/LASCO C3 camera. A magnetic cloud (MC) corresponding to this CME was measured in-situ by the Wind spacecraft two days later. By applying two MC reconstruction methods, we infer the configuration of the MC as well as some kinematic information, which implies that the CME possibly experienced an eastward deflection on its way to 1 AU. However, due to the lack of observations from the STEREO spacecraft, the CME's kinematic evolution in interplanetary space is not clear. In order to fill this gap, we utilize numerical MHD simulation, drag-based CME propagation model (DBM) and the model for CME deflection in interplanetary space (DIPS) to recover the propagation process, especially the trajectory, of the CME from $30 R_S$ to 1 AU. It is suggested that the trajectory of the CME was deflected toward the Earth by about $12^\circ$, consistent with the implication from the MC reconstruction at 1 AU. This eastward deflection probably contributed to the CME's unexpected geoeffectiveness by pushing the center of the initially west-oriented CME closer to the Earth.

We study the physical mechanism of a major X-class solar flare that occurred in the super NOAA active region (AR) 12192 using a data-driven numerical magnetohydrodynamic (MHD) modeling complemented with observations. With the evolving magnetic fields observed at the solar surface as bottom boundary input, we drive an MHD system to evolve self-consistently in correspondence with the realistic coronal evolution. During a two-day time interval, the modeled coronal field has been slowly stressed by the photospheric field evolution,which gradually created a large-scale coronal current sheet, i.e., a narrow layer with intense current, in the core of the AR. The current layer was successively enhanced until it became so thin that a tether-cutting reconnection between the sheared magnetic arcades was set in, which led to a flare. The modeled reconnecting field lines and their footpoints match well the observed hot flaring loops and the flare ribbons, respectively, suggesting that the model has successfully "reproduced" the macroscopic magnetic process of the flare. In particular, with simulation, we explained why this event is a confined eruption-the consequent of the reconnection is the shared arcade instead of a newly formed flux rope. We also found much weaker magnetic implosion effect comparing to many other X-class flares

Magnetic clouds (MCs) are the interplanetary counterpart of coronal magnetic flux ropes. They can provide valuable information to reveal the flux rope characteristics at their eruption stage in the corona, which are unable to be explored in situ at present. In this paper, we make a comprehensive survey of the average iron charge state (<Q>Fe) distributions inside 96 MCs for solar cycle 23 using ACE (Advanced Composition Explorer) data. As the <Q>Fe in the solar wind are typically around 9+ to 11+, the Fe charge state is defined as high when the <Q>Fe is larger than 12+, which implies the existence of a considerable amount of Fe ions with high charge states (e.g., \geq 16+). The statistical results show that the <Q>Fe distributions of 92 (~ 96%) MCs can be classified into four groups with different characteristics. In group A (11 MCs), the <Q>Fe shows a bimodal distribution with both peaks higher than 12+. Group B (4 MCs) presents a unimodal distribution of <Q>Fe with its peak higher than 12+. In groups C (29 MCs) and D (48 MCs), the <Q>Fe remains higher and lower than 12+ throughout ACE passage through the MC, respectively. Possible explanations to these distributions are discussed.

Solar active region (AR) 11283 is a very magnetically complex region and it has produced many eruptions. However, there exists a non-eruptive filament in the plage region just next to an eruptive one in the AR, which gives us an opportunity to perform a comparison analysis of these two filaments. The coronal magnetic field extrapolated using a CESE-MHD-NLFFF code (Jiang & Feng 2013) reveals that two magnetic flux ropes (MFRs) exist in the same extrapolation box supporting these two filaments, respectively. Analysis of the magnetic field shows that the eruptive MFR contains a bald-patch separatrix surface (BPSS) co-spatial very well with a pre-eruptive EUV sigmoid, which is consistent with the BPSS model for coronal sigmoids. The magnetic dips of the non-eruptive MFRs match H{\alpha} observation of the non-eruptive filament strikingly well, which strongly supports the MFR-dip model for filaments. Compared with the non-eruptive MFR/filament (with a length of about 200 Mm), the eruptive MFR/filament is much smaller (with a length of about 20 Mm), but it contains most of the magnetic free energy in the extrapolation box and holds a much higher free energy density than the non-eruptive one. Both the MFRs are weakly twisted and cannot trigger kink instability. The AR eruptive MFR is unstable because its axis reaches above a critical height for torus instability, at which the overlying closed arcades can no longer confine the MFR stably. On the contrary, the quiescent MFR is very firmly held by its overlying field, as its axis apex is far below the torus-instability threshold height. Overall, this comparison investigation supports that MFR can exist prior to eruption and the ideal MHD instability can trigger MFR eruption.

Hot channels (HCs), high temperature erupting structures in the lower corona of the Sun, have been proposed as a proxy of magnetic flux ropes (MFRs) since their initial discovery. However, it is difficult to make definitive proof given the fact that there is no direct measurement of magnetic field in the corona. An alternative way is to use the magnetic field measurement in the solar wind from in-situ instruments. On 2012 July 12, an HC was observed prior to and during a coronal mass ejection (CME) by the AIA high-temperature images. The HC is invisible in the EUVI low-temperature images, which only show the cooler leading front (LF). However, both the LF and an ejecta can be observed in the coronagraphic images. These are consistent with the high temperature and high density of the HC and support that the ejecta is the erupted HC. In the meanwhile, the associated CME shock was identified ahead of the ejecta and the sheath through the COR2 images, and the corresponding ICME was detected by \textit{ACE}, showing the shock, sheath and magnetic cloud (MC) sequentially, which agrees with the coronagraphic observations. Further, the MC contained a low-ionization-state center and a high-ionization-state shell, consistent with the pre-existing HC observation and its growth through magnetic reconnection. All of these observations support that the MC detected near the Earth is the counterpart of the erupted HC in the corona for this event. Therefore, our study provides strong observational evidence of the HC as an MFR.

We report on the detailed and systematic study of field-line twist and length distributions within magnetic flux ropes embedded in Interplanetary Coronal Mass Ejections (ICMEs). The Grad-Shafranov reconstruction method is utilized together with a constant-twist nonlinear force-free (Gold-Hoyle) flux rope model to reveal the close relation between the field-line twist and length in cylindrical flux ropes, based on in-situ Wind spacecraft measurements. We show that the field-line twist distributions within interplanetary flux ropes are inconsistent with the Lundquist model. In particular we utilize the unique measurements of magnetic field-line lengths within selected ICME events as provided by Kahler et al. (2011) based on energetic electron burst observations at 1 AU and the associated type III radio emissions detected by the Wind spacecraft. These direct measurements are compared with our model calculations to help assess the flux-rope interpretation of the embedded magnetic structures. By using the different flux-rope models, we show that the in-situ direct measurements of field-line lengths are consistent with a flux-rope structure with spiral field lines of constant and low twist, largely different from that of the Lundquist model, especially for relatively large-scale flux ropes.

Hydromagnetic waves, especially those of frequencies in the range of a few milli-Hz to a few Hz observed in the Earth's magnetosphere, are categorized as Ultra Low Frequency (ULF) waves or pulsations. They have been extensively studied due to their importance in the interaction with radiation belt particles and in probing the structures of the magnetosphere. We developed an approach in examining the toroidal standing Aflv\'{e}n waves in a background magnetic field by recasting the wave equation into a Klein-Gordon (KG) form along individual field lines. The eigenvalue solutions to the system are characteristic of a propagation type when the corresponding eigen-frequency is greater than a cut-off frequency and an evanescent type otherwise. We apply the approach to a compressed dipole magnetic field model of the inner magnetosphere, and obtain the spatial profiles of relevant parameters and the spatial wave forms of harmonic oscillations. We further extend the approach to poloidal mode standing Alfv\'{e}n waves along field lines. In particular, we present a quantitative comparison with a recent spacecraft observation of a poloidal standing Alfv\'{e}n wave in the Earth's magnetosphere. Our analysis based on KG equation yields consistent results which agree with the spacecraft measurements of the wave period and the amplitude ratio between the magnetic field and electric field perturbations. We also present computational results of eigenvalue solutions to the compressional poloidal mode waves in the compressed dipole magnetic field geometry.

Magnetic reconnection is essential to release the flux rope during its ejection. The question remains: how does the magnetic reconnection change the flux rope structure? Following the original study of \citet{Qiu2007}, we compare properties of ICME/MC flux ropes measured at 1 AU and properties of associated solar progenitors including flares, filaments, and CMEs. In particular, the magnetic field-line twist distribution within interplanetary magnetic flux ropes is systematically derived and examined. Our analysis shows that for most of these events, the amount of twisted flux per AU in MCs is comparable with the total reconnection flux on the Sun, and the sign of the MC helicity is consistent with the sign of helicity of the solar source region judged from the geometry of post-flare loops. Remarkably, we find that about one half of the 18 magnetic flux ropes, most of them being associated with erupting filaments, have a nearly uniform and relatively low twist distribution from the axis to the edge, and the majority of the other flux ropes exhibit very high twist near the axis, of up to $\gtrsim 5$ turns per AU, which decreases toward the edge. The flux ropes are therefore not linear force free. We also conduct detailed case studies showing the contrast of two events with distinct twist distribution in MCs as well as different flare and dimming characteristics in solar source regions, and discuss how reconnection geometry reflected in flare morphology may be related to the structure of the flux rope formed on the Sun.

Solar filament are commonly thought to be supported in magnetic dips, in particular, of magnetic flux ropes (FRs). In this Letter, from the observed photospheric vector magnetogram, we implement a nonlinear force-free field (NLFFF) extrapolation of a coronal magnetic FR that supports a large-scale intermediate filament between an active region and a weak polarity region. This result is the first in that current NLFFF extrapolations with presence of FRs are limited to relatively small-scale filaments that are close to sunspots and along main polarity inversion line (PIL) with strong transverse field and magnetic shear, and the existence of a FR is usually predictable. In contrast, the present filament lies along the weak-field region (photospheric field strength $\lesssim 100$ G), where the PIL is very fragmented due to small parasitic polarities on both side of the PIL and the transverse field has a low value of signal-to-noise ratio. Thus it represents a far more difficult challenge to extrapolate a large-scale FR in such case. We demonstrate that our CESE--MHD--NLFFF code is competent for the challenge. The numerically reproduced magnetic dips of the extrapolated FR match observations of the filament and its barbs very well, which supports strongly the FR-dip model for filaments. The filament is stably sustained because the FR is weakly twisted and strongly confined by the overlying closed arcades.

The New Vacuum Solar Telescope (NVST) is a 1 meter vacuum solar telescope that aims to observe the fine structures on the Sun. The main tasks of NVST are high resolution imaging and spectral observations, including the measurements of solar magnetic field. NVST is the primary ground-based facility of Chinese solar community in this solar cycle. It is located by the Fuxian Lake of southwest China, where the seeing is good enough to perform high resolution observations. In this paper, we first introduce the general conditions of Fuxian Solar Observatory and the primary science cases of NVST. Then, the basic structures of this telescope and instruments are described in detail. Finally, some typical high resolution data of solar photosphere and chromosphere are also shown.

We present a magnetic analysis of the formation and eruption of an active region (AR) sigmoid in AR 11283 from 2011 September 4 to 6. To follow the quasi-static evolution of the coronal magnetic field, we reconstruct a time sequence of static fields using a recently developed nonlinear force-free field model constrained by the SDO/HMI vector magnetograms. A detailed analysis of the fields compared with the SDO/AIA observations suggests the following scenario for the evolution of the region. Initially, a new bipole emerges into the negative polarity of a pre-existing bipolar AR, forming a null point topology between the two flux systems. A weakly twisted flux rope (FR) is then built up slowly in the embedded core region, largely through flux-cancellation photospheric reconnections, forming a bald patch separatrix surface (BPSS) separating the FR from its ambient field. The FR grows gradually until its axis runs into a torus instability (TI) domain near the end of the third day, and the BPSS also develops a fully S-shape. Unlike in the case of standard TI, the FR does not erupt instantly since it is still attached at the photosphere along the bald patch (BP) portion of the polarity inversion line. The combined effects of the TI-driven expansion of the FR and the line-tying at the BP tear the FR into two parts with the upper portion freely expelled and the lower portion remaining behind the post-flare arcades. This process dynamically perturbs the BPSS and results in the transient enhanced brightening of the sigmoid. The accelerated expansion of the upper portion of the FR strongly pushes its envelope flux near the null point and triggers breakout reconnection at the null, as evidenced by a remarkable circular flare ribbon, which further facilitates the eruption. We discuss the important implications of these results for the formation and disruption of sigmoid region with FR.

Two major processes have been proposed to convert the coronal magnetic energy into the kinetic energy of a coronal mass ejection (CME): resistive magnetic reconnection and ideal macroscopic magnetohydrodynamic instability of magnetic flux rope. However, it remains elusive whether both processes play a comparable role or one of them prevails during a particular eruption. To shed light on this issue, we carefully studied energetic but flareless CMEs, \textit{i.e.}, fast CMEs not accompanied by any flares. Through searching the Coordinated Data Analysis Workshops (CDAW) database of CMEs observed in Solar Cycle 23, we found 13 such events with speeds larger than 1000 km s$^{-1}$. Other common observational features of these events are: (1) none of them originated in active regions; they were associated with eruptions of well-developed long filaments in quiet-Sun regions, (2) no apparent enhancement of flare emissions was present in soft X-ray, EUV and microwave data. Further studies of two events reveal that (1) the reconnection electric fields, as inferred from the product of the separation speed of post-eruption ribbons and the photospheric magnetic field measurement, were in general weak; (2) the period with a measurable reconnection electric field is considerably shorter than the total filament-CME acceleration time. These observations indicate that, for these fast CMEs, the magnetic energy was released mainly via the ideal flux rope instability through the work done by the large scale Lorentz force acting on the rope currents rather than via magnetic reconnections. We also suggest that reconnections play a less important role in accelerating CMEs in quiet Sun regions of weak magnetic field than those in active regions of strong magnetic field.

Current magnetohydrodynamic (MHD) simulations of the initiation of solar eruptions are still commonly carried out with idealized magnetic field models, whereas the realistic coronal field prior to eruptions can possibly be reconstructed from the observable photospheric field. Using a nonlinear force-free field extrapolation prior to a sigmoid eruption in AR 11283 as the initial condition in a MHD model, we successfully simulate the realistic initiation process of the eruption event, as is confirmed by a remarkable resemblance to the SDO/AIA observations. Analysis of the pre-eruption field reveals that the envelope flux of the sigmoidal core contains a coronal null and furthermore the flux rope is prone to a torus instability. Observations suggest that reconnection at the null cuts overlying tethers and likely triggers the torus instability of the flux rope, which results in the eruption. This kind of simulation demonstrates the capability of modeling the realistic solar eruptions to provide the initiation process.

We report on a heterodyne receiver designed to observe the astrophysically important neutral atomic oxygen [OI] line at 4.7448 THz. The local oscillator is a third-order distributed feedback Quantum Cascade Laser operating in continuous wave mode at 4.741 THz. A quasi-optical, superconducting NbN hot electron bolometer is used as the mixer. We recorded a double sideband receiver noise temperature (T^DSB_rec) of 815 K, which is ~7 times the quantum noise limit (h{\nu}/2k_B) and an Allan variance time of 15 s at an effective noise fluctuation bandwidth of 18 MHz. Heterodyne performance was confirmed by measuring a methanol line spectrum.

We apply a data-driven MHD model to investigate the three-dimensional (3D) magnetic field of NOAA active region (AR) 11117 around the time of a C-class confined flare occurred on 2010 October 25. The MHD model, based on the spacetime conservation-element and solution-element (CESE) scheme, is designed to focus on the magnetic-field evolution and to consider a simplified solar atomsphere with finite plasma $\beta$. Magnetic vector-field data derived from the observations at the photoshpere is inputted directly to constrain the model. Assuming that the dynamic evolution of the coronal magnetic field can be approximated by successive equilibria, we solve a time sequence of MHD equilibria basing on a set of vector magnetograms for AR 11117 taken by the Helioseismic and Magnetic Imager (HMI) on board the {\it Solar Dynamic Observatory (SDO)} around the time of flare. The model qualitatively reproduces the basic structures of the 3D magnetic field, as supported by the visual similarity between the field lines and the coronal loops observed by the Atmospheric Imaging Assembly (AIA), which shows that the coronal field can indeed be well characterized by the MHD equilibrium in most time. The magnetic configuration changes very limited during the studied time interval of two hours. A topological analysis reveals that the small flare is correlated with a bald patch (BP, where the magnetic field is tangent to the photoshpere), suggesting that the energy release of the flare can be understood by magnetic reconnection associated with the BP separatrices. The total magnetic flux and energy keep increasing slightly in spite of the flare, while the computed magnetic free energy drops during the flare with an amount of $\sim 10^{30}$ erg, which seems to be adequate to provide the energy budget of the minor C-class confined flare.

We present a newly developed approach to solar coronal magnetic field extrapolation from vector magnetograms, based on the Principle of Minimum Dissipation Rate (MDR). The MDR system was derived from a variational problem that is more suitable for an open and externally driven system, like the solar corona. The resulting magnetic field equation is more general than force-free. Its solution can be expressed as the superposition of two linear (constant-$\alpha$) force-free fields (LFFFs) with distinct $\alpha$ parameters, and one potential field. Thus the original extrapolation problem is decomposed into three LFFF extrapolations, utilizing boundary data. The full MDR-based approach requires two layers of vector magnetograph measurements on solar surface, while a slightly modified practical approach only requires one. We test both approaches against 3D MHD simulation data in a finite volume. Both yield quantitatively good results. The errors in the magnetic energy estimate are within a few percents. In particular, the main features of relatively strong perpendicular current density structures, representative of the non-force freeness of the solution, are well recovered.

Coronal mass ejections (CMEs) observed near the Sun via LASCO coronographic imaging are the most important solar drivers of geomagnetic storms. ICMEs, their interplanetary, near-Earth counterparts, can be detected in-situ, for example, by the Wind and ACE spacecraft. An ICME usually exhibits a complex structure that very often includes a magnetic cloud (MC). They can be commonly modelled as magnetic flux ropes and there is observational evidence to expect that the orientation of a halo CME elongation corresponds to the orientation of the flux rope. In this study, we compare orientations of elongated CME halos and the corresponding MCs, measured by Wind and ACE spacecraft. We characterize the MC structures by using the Grad-Shafranov reconstruction technique and three MC fitting methods to obtain their axis directions. The CME tilt angles and MC fitted axis angles were compared without taking into account handedness of the underlying flux rope field and the polarity of its axial field. We report that for about 64% of CME-MC events, we found a good correspondence between the orientation angles implying that for the majority of interplanetary ejecta their orientations do not change significantly (less than 45 deg rotation) while travelling from the Sun to the near Earth environment.

The SDSS Data Release 1 includes 1833 DA white dwarfs (WDs) and forms the largest homogeneous sample of WDs. This sample provides the best opportunity to study the statistical properties of WDs. We adopt a recently established theoretical model to calculate the mass and distance of each WD using the observational data. Then we adopt a bin-correction method to correct for selection effects and use the 1/V weight-factor method to calculate the luminosity function, the continuous mass function and the formation rate of these WDs. The SDSS DA WD sample is incomplete and suffers seriously from selection effects. After corrections for the selection effects, only 531 WDs remain. From this final sample we derive the most up-to-date luminosity function and mass function, in which we find a broad peak of WD masses centered around 0.58$M_{\odot}$. The DA WD space density is calculated as $8.81\times10^{-5}pc^{-3}$ and the formation rate is $2.579\times 10^{-13}pc^{-3}yr^{-1}$. The statistical properties of the SDSS DA WD sample are generally in good agreement with previous observational and theoretical studies, and provide us information on the formation and evolution of WDs. However, a larger and more complete all-sky WD sample is still needed to explain some subtle disagreements and unresolved issues.