Similar to the compactness parameter ($\beta=M/R$), the gravitational binding energy (GBE) is also a characteristic parameter which can reflect the internal structure of a neutron star and thus can be used to expressing the universal relations. Scaling by the stellar mass, this investigation demonstrates a perfect universal relation between the GBE and the moment of inertia, where both of the normal neutron stars and the quark stars satisfy the same universal relation. Moreover, a fine empirical relation between the GBE and the tidal deformability is proposed, where the difference of the relations can be used to distinguish wether a pulsar is a normal neutron star or a quark star if the stellar mass and the tidal deformability can be observed or estimated rather accurately. These universal relations provide a potential way to estimate the GBE if the stellar mass and the moment of inertia/the tidal deformability are precisely measured.

A new code, named MAP, is written in Fortran language for magnetohydrodynamics (MHD) calculation with the adaptive mesh refinement (AMR) and Message Passing Interface (MPI) parallelization. There are several optional numerical schemes for computing the MHD part, namely, modified Mac Cormack Scheme (MMC), Lax-Friedrichs scheme (LF) and weighted essentially non-oscillatory (WENO) scheme. All of them are second order, two-step, component-wise schemes for hyperbolic conservative equations. The total variation diminishing (TVD) limiters and approximate Riemann solvers are also equipped. A high resolution can be achieved by the hierarchical block-structured AMR mesh. We use the extended generalized Lagrange multiplier (EGLM) MHD equations to reduce the non-divergence free error produced by the scheme in the magnetic induction equation. The numerical algorithms for the non-ideal terms, e.g., the resistivity and the thermal conduction, are also equipped in the MAP code. The details of the AMR and MPI algorithms are described in the paper.

Microflares are small activities in solar low atmosphere, some are in the low corona, and others in the chromosphere. Observations show that some of the microflares are triggered by magnetic reconnection between emerging flux and a pre-existing background magnetic field. We perform 2.5D compressible resistive MHD simulations of magnetic reconnection with gravity considered. The background magnetic field is a canopy-type configuration which is rooted at the boundary of the solar supergranule. By changing the bottom boundary conditions in the simulation, new magnetic flux emerges up at the center of the supergranule and reconnects with the canopy-type magnetic field. We successfully simulate the coronal and chromospheric microflares, whose current sheets are located at the corona and the chromosphere, respectively. The microflare of coronal origin has a bigger size and a higher temperature enhancement than that of chromospheric origin. In the microflares of coronal origin, we also found a hot jet ($\sim$$1.8 \times 10^6$ K), which is probably related to the observational EUV/SXR jets, and a cold jet ($\sim$$10^4$ K), which is similar to the observational H$\alpha$/Ca surges, whereas there is only an H$\alpha$/Ca bright point in the microflares of chromospheric origin. The study of parameter dependence shows that the size and strength of the emerging magnetic flux are the key parameters which determine the height of the reconnection location, and further determine the different observational features of the microflares.

We performed three dimensional resistive magnetohydrodynamic simulations to study the magnetic reconnection using an initially shearing magnetic field configuration (force free field with a current sheet in the middle of the computational box). It is shown that there are two types of reconnection jets: the ordinary reconnection jets and fan-shaped jets, which are formed along the guide magnetic field. The fan-shaped jets are much different from the ordinary reconnection jets which are ejected by magnetic tension force. There are two driving forces for accelerating the fan-shaped jets. The one is the Lorentz force which dominates the motion of fluid elements at first and then the gas pressure gradient force accelerates the fluid elements in the later stage. The dependence on magnetic reconnection angle and resistivity value has also been studied. The formation and evolution of these jets provide a new understanding of dynamic magnetohydrodynamic jets.

Magnetic reconnection is a fundamental process in space and astrophysical plasmas in which oppositely directed magnetic fields changes its connectivity and eventually converts its energy into kinetic and thermal energy of the plasma. Recently, ubiquitous jets (for example, chromospheric anemone jets, penumbral microjets, umbral light bridge jets) have been observed by Solar Optical Telescope on board the satellite Hinode. These tiny and frequently occurring jets are considered to be a possible evidence of small-scale ubiquitous reconnection in the solar atmosphere. However, the details of three dimensional magnetic configuration are still not very clear. Here we propose a new model based on three dimensional simulations of magnetic reconnection using a typical current sheet magnetic configuration with a strong guide field. The most interesting feature is that the jets produced by the reconnection eventually move along the guide field lines. This model provides a fresh understanding of newly discovered ubiquitous jets and moreover a new observational basis for the theory of astrophysical magnetic reconnection.

With gravity, ionization, and radiation being considered, we perform 2.5D compressible resistive MHD simulations of chromospheric magnetic reconnection using the CIP-MOCCT scheme. The temperature distribution of the quiet-Sun atmospheric model VALC and the helium abundance (10%) are adopted. Our 2.5D MHD simulation reproduces qualitatively the temperature enhancement observed in chromospheric microflares. The temperature enhancement $\Delta T$ is demonstrated to be sensitive to the background magnetic field, whereas the total evolution time $\Delta t$ is sensitive to the magnitude of the anomalous resistivity. Moveover, we found a scaling law, which is described as $\Delta T/\Delta t \sim {n_H}^{-1.5} B^{2.1} {\eta_0}^{0.88}$. Our results also indicate that the velocity of the upward jet is much greater than that of the downward jet and the X-point may move up or down.