In this paper, we report our multiwavelength observations of subarcsecond blobs in coronal jets. In AR 12149, a C5.5 circular-ribbon flare occurred on 2014 August 24. Two jets (jet1 and jet2) were related to the flare. Jet1 appeared first. Jet2 appeared 6 minutes later. During its initial phase, a big plasmoid was ejected out of jet2. After the flare peak time, multiple bright and compact blobs appeared in the lower part of jet2, which were observed by IRIS/SJI. The blobs observed by SJI in 1330 {\AA} have sizes of 0$\farcs$45$-$1$\farcs$35, nearly 84\% of which are subarcsecond ($<$1\arcsec). The mean value and standard deviation of the sizes are 0$\farcs$78 and 0$\farcs$19, respectively. The velocities of the blobs range from 10 to more than 220 km s$^{-1}$. Three of the blobs had their counterparts in EUV wavelengths observed by SDO/AIA. The velocities are almost identical in UV and EUV wavelengths. We propose that the blobs observed in 1330 {\AA} are the cool component ($\sim$0.025 MK), while the blobs observed in EUV are the hot component of several MK. In jet1, only one blob was present. We conclude that the blobs are created by the tearing-mode instability of the current sheet at the base or inside the coronal jets. Our results have important implication for uncovering the fine structures of coronal jets and understanding the relationship between the blobs observed in UV and EUV wavelengths.

Magnetic reconnection is the most likely mechanism responsible for the high temperature events that are observed in strongly magnetized locations around the temperature minimum in the low solar chromosphere. This work improves upon our previous work ["Magnetic Reconnection in Strongly Magnetized Regions of the Low Solar Chromosphere", The Astrophysical Journal 852, 95 (2018) ] by using a more realistic radiative cooling model computed from the OPACITY project and the CHIANTI database. We find that the rate of ionization of the neutral component of the plasma is still faster than recombination within the current sheet region. For low $\beta$ plasmas, the ionized and neutral fluid flows are well-coupled throughout the reconnection region resembling the single-fluid Sweet-Parker model dynamics. Decoupling of the ion and neutral inflows appears in the higher $\beta$ case with $\beta_0=1.46$, which leads to a reconnection rate about three times faster than the rate predicted by the Sweet-Parker model. In all cases, the plasma temperature increases with time inside the current sheet, and the maximum value is above $2\times10^4$~K when the reconnection magnetic field strength is greater than $500$~G. While the more realistic radiative cooling model does not result in qualitative changes of the characteristics of magnetic reconnection, it is necessary for studying the variations of the plasma temperature and ionization fraction inside current sheets in strongly magnetized regions of the low solar atmosphere. It is also important for studying energy conversion during the magnetic reconnection process when the hydrogen-dominated plasma approaches full ionization.

We numerically study magnetic reconnection on different spatial scales and at different heights in the weakly ionized plasma of the low solar atmosphere (around $300-800$~km above the solar surface) within a reactive 2.5 D multi-fluid plasma-neutral model. We consider a strongly magnetized plasma ($\beta \sim 6\% $) evolving from a force-free magnetic configuration and perturbed to initialize formation of a reconnection current sheet. On large scales, the resulting current sheets are observed to undergo a secondary 'plasmoid' instability. A series of simulations at different scales demonstrate a cascading current sheet formation process that terminates for current sheets with width of ~2m and length of $\sim100$~m, corresponding to the critical current sheet aspect ratio of $\sim50$. We also observe that the plasmoid instability is the primary physical mechanism accelerating the magnetic reconnection in this plasma parameter regime. After plasmoid instabilities appear, the reconnection rate sharply increases to a value of $\sim$ 0.035, observed to be independent of the Lundquist number. These characteristics are very similar to magnetic reconnection in fully ionized plasmas. In this low $\beta$ guide field reconnection regime, both the recombination and collisionless effects are observed to have a small contribution to the reconnection rate. The simulations show that it is difficult to heat the dense weakly ionized photospheric plasmas to above $2\times10^4$~K during the magnetic reconnection process. However, the plasmas in the low solar chromosphere can be heated above $3\times10^4$~K with reconnection magnetic fields of $500$~G or stronger.

The Kelvin-Hemholtz (K-H) instability in the corona EUV jet is studied via 2.5 MHD numerical simulations. The jet results from magnetic reconnection due to the interation of the new emerging magnetic field and the pre-existing magnetic field in the corona. Our results show that the Alfven Mach number along the jet is about 5-14 just before the instability occurs, and it is even higher than 14 at some local areas. During the K-H instability process, several vortex-like plasma blobs of high temperature and high density appear along the jet, and magnetic fields have also been rolled up and the magnetic configuration including anti-parallel magnetic fields forms, which leads to magnetic reconnection at many X-points and current sheet fragments inside the vortex-like blob. After magnetic islands appear inside the main current sheet, the total kinetic energy of the reconnection outflows decreases, and cannot support the formation of the vortex-like blob along the jet any longer, then the K-H instability eventually disappears. We also present the results about how the guide field and the flux emerging speed affect the K-H instability. We find that the strong guide field inhibits the shock formation in the reconnection upward outflow regions but helps secondary magnetic islands appear earlier in the main current sheet, and then apparently suppresses the K-H instability. As the speed of the emerging magnetic field decreases, the K-H instability appears later, the highest temperature inside the vortex blob gets lower and the vortex structure gets smaller.

Solid evidence of magnetic reconnection is rarely reported within sunspots, the darkest regions with the strongest magnetic fields and lowest temperatures in the solar atmosphere. Using the world's largest solar telescope, the 1.6-meter Goode Solar Telescope, we detect prevalent reconnection through frequently occurring fine-scale jets in the H${\alpha}$ line wings at light bridges, the bright lanes that may divide the dark sunspot core into multiple parts. Many jets have an inverted Y-shape, shown by models to be typical of reconnection in a unipolar field environment. Simultaneous spectral imaging data from the Interface Region Imaging Spectrograph show that the reconnection drives bidirectional flows up to 200~km~s$^{-1}$, and that the weakly ionized plasma is heated by at least an order of magnitude up to $\sim$80,000 K. Such highly dynamic reconnection jets and efficient heating should be properly accounted for in future modeling efforts of sunspots. Our observations also reveal that the surge-like activity previously reported above light bridges in some chromospheric passbands such as the H${\alpha}$ core has two components: the ever-present short surges likely to be related to the upward leakage of magnetoacoustic waves from the photosphere, and the occasionally occurring long and fast surges that are obviously caused by the intermittent reconnection jets.

Magnetic reconnection in strongly magnetized regions around the temperature minimum region of the low solar atmosphere is studied by employing MHD-based simulations of a partially ionized plasma within a reactive 2.5D multi-fluid model. It is shown that in the absence of magnetic nulls in a low $\beta$ plasma the ionized and neutral fluid flows are well-coupled throughout the reconnection region. However, non-equilibrium ionization-recombination dynamics play a critical role in determining the structure of the reconnection region, lead to much lower temperature increases and a faster magnetic reconnection rate as compared to simulations that assume plasma to be in ionization-recombination equilibrium. The rate of ionization of the neutral component of the plasma is always faster than recombination within the current sheet region even when the initial plasma $\beta$ is as high as $\beta_0=1.46$. When the reconnecting magnetic field is in excess of a kilogauss and the plasma $\beta$ is lower than 0.0145, the initially weakly ionized plasmas can become fully ionized within the reconnection region and the current sheet can be strongly heated to above $2.5\times10^4$~K, even as most of the collisionally dissipated magnetic energy is radiated away. The Hall effect increases the reconnection rate slightly, but in the absence of magnetic nulls it does not result in significant asymmetries or change the characteristics of the reconnection current sheet down to meter scales.

We perform two-dimensional resistive magnetohydrodynamic simulations of coronal jets driven by flux emergence along the lower boundary. The reconnection layers are susceptible the formation of blobs that are ejected in the jet. Our simulation with low plasma $\beta$ (Case I) shows that magnetic islands form easily and propagate upwards in the jet. These islands are multithermal and thus are predicted to show up in hot channels (335 \AA\, and 211 \AA) and the cool channel (304 \AA) in observations by the Atmospheric Imaging Assembly (AIA) on the \emph{Solar Dynamics Observatory}. The islands have maximum temperatures of 8 MK, lifetimes of 120 s, diameters of 6 Mm, and velocities of 200 km s$^{-1}$. These parameters are similar to the properties of blobs observed in EUV jets by AIA\@. The Kelvin-Helmholtz instability develops in our simulation with moderately high plasma $\beta$ (Case II), and leads to the formation of bright vortex-like blobs above the multiple high magneto-sonic Mach number regions that appear along the jet. These vortex-like blobs can also be identified in the AIA channels. However, they eventually move downward and disappear after the high magneto-sonic Mach number regions disappear. In the lower plasma $\beta$ case, the lifetime for the jet is shorter, the jet and magnetic islands are formed with higher velocities and temperatures, the current sheet fragments are more chaotic, and more magnetic islands are generated. Our results show that the plasmoid instability and Kelvin-Helmholtz instability along the jet are both possible causes of the formation of blobs observed at extreme ultraviolet (EUV) wavelengths.

We simulate several magnetic reconnection processes in the low solar chromosphere/photosphere, the radiation cooling, heat conduction and ambipolar diffusion are all included. Our numerical results indicate that both the high temperature($ \gtrsim 8\times10^4$~K) and low temperature($\sim 10^4$~K) magnetic reconnection events can happen in the low solar atmosphere ($100\sim600$~km above the solar surface). The plasma $\beta$ controlled by plasma density and magnetic fields is one important factor to decide how much the plasma can be heated up. The low temperature event is formed in a high $\beta$ magnetic reconnection process, Joule heating is the main mechanism to heat plasma and the maximum temperature increase is only several thousand Kelvin. The high temperature explosions can be generated in a low $\beta$ magnetic reconnection process, slow and fast-mode shocks attached at the edges of the well developed plasmoids are the main physical mechanisms to heat the plasma from several thousand Kelvin to over $8\times10^4$~K. Gravity in the low chromosphere can strongly hinder the plasmoind instability and the formation of slow-mode shocks in a vertical current sheet. Only small secondary islands are formed; these islands, however, are not well developed as those in the horizontal current sheets. This work can be applied for understanding the heating mechanism in the low solar atmosphere and could possibly be extended to explain the formation of common low temperature EBs ($\sim10^4$~K) and the high tenperature IRIS bombs ($\gtrsim 8\times10^4$) in the future.

In this paper, we study the energy conversion and spectra in a corona current sheet by 2.5-dimensional MHD numerical simulations. Numerical results show that many Petschek-like fine structures with slow-mode shocks mediated by plasmoid instabilities develop during the magnetic reconnection process. The termination shocks can also be formed above the primary magnetic island and at the head of secondary islands. These shocks play important roles in generating thermal energy in a corona current sheet. For a numerical simulation with initial conditions close to the solar corona environment, the ratio of the generated thermal energy to the total dissipated magnetic energy is around $1/5$ before secondary islands appear. After secondary islands appear, the generated thermal energy starts to increase sharply and this ratio can reach a value about $3/5$. In an environment with a relatively lower plasma density and plasma $\beta$, the plasma can be heated to a much higher temperature. After secondary islands appear, the one dimensional energy spectra along the current sheet do not behave as a simple power law and the spectrum index increases with the wave number. The average spectrum index for the magnetic energy spectrum along the current sheet is about $1.8$. The two dimensional spectra intuitively show that part of the high energy is cascaded to large $kx$ and $ky$ space after secondary islands appear. The plasmoid distribution function calculated from numerical simulations behaves as a power law closer to $f(\psi) \sim \psi^{-1}$ in the intermediate $\psi$ regime. By using $\eta_{eff} = v_{inflow}\cdot L$, the effective magnetic diffusivity is estimated about $10^{11}\sim10^{12}$~m$^2$\,s$^{-1}$.

Magnetic reconnection in the partially ionized solar chromosphere is studied in 2.5-dimensional magnetohydrodynamic simulations including radiative cooling and ambipolar diffusion. A Harris current sheet with and without a guide field is considered. Characteristic values of the parameters in the middle chromosphere imply a high magnetic Reynolds number of $\sim10^{6}\mbox{--}10^7$ in the present simulations. Fast magnetic reconnection then develops as a consequence of the plasmoid instability without the need to invoke anomalous resistivity enhancements. Multiple levels of the instability are followed as it cascades to smaller scales, which approach the ion inertial length. The reconnection rate, normalized to the asymptotic values of magnetic field and Alfv\'en velocity in the inflow region, reaches values in the range $\sim0.01\mbox{--}0.03$ throughout the cascading plasmoid formation and for zero as well as for strong guide field. The out-flow velocity reaches $\approx40$~km\,s$^{-1}$. Slow-mode shocks extend from the $X$-points, heating the plasmoids up to $\sim 8\times10^4$~K. In the case of zero guide field, the inclusion of ambipolar diffusion and radiative cooling both cause a rapid thinning of the current sheet (down to $\sim30$~m) and early formation of secondary islands. Both of these processes have very little effect on the plasmoid instability for a strong guide field. The reconnection rates, temperature enhancements, and upward out-flow velocities from the vertical current sheet correspond well to their characteristic values in chromospheric jets.

In this paper we investigate, by means of two-dimensional magnetohydrodynamic simulations, the impact of temperature-dependent resistivity and thermal conduction on the development of plasmoid instabilities in reconnecting current sheets in the solar corona. We find that the plasma temperature in the current sheet region increases with time and it becomes greater than that in the inflow region. As secondary magnetic islands appear, the highest temperature is not always found at the reconnection $X$-points, but also inside the secondary islands. One of the effects of anisotropic thermal conduction is to decrease the temperature of the reconnecting $X-$points and transfer the heat into the $O-$points, the plasmoids, where it gets trapped. In the cases with temperature-dependent magnetic diffusivity, $\eta \sim T^{-3/2}$, the decrease in plasma temperature at the $X-$points leads to: (i) increase in the magnetic diffusivity until the characteristic time for magnetic diffusion becomes comparable to that of thermal conduction; (ii) increase in the reconnection rate; and, (iii) more efficient conversion of magnetic energy into thermal energy and kinetic energy of bulk motions. These results provide further explanation of the rapid release of magnetic energy into heat and kinetic energy seen during flares and coronal mass ejections. In this work, we demonstrate that the consideration of anisotropic thermal conduction and Spitzer-type, temperature-dependent magnetic diffusivity, as in the real solar corona, are crucially important for explaining the occurrence of fast reconnection during solar eruptions.