Waves and shocks traveling through the solar chromospheric plasma are influenced by its partial ionization and weak collisional coupling, and may become susceptible to multi-fluid effects, similar to interstellar shock waves. In this study, we consider fast magneto-acoustic shock wave formation and propagation in a stratified medium, that is permeated by a horizontal magnetic field, with properties similar to that of the solar chromosphere. The evolution of plasma and neutrals is modeled using a two-fluid code that evolves a set of coupled equations for two separate fluids. We observed that waves in neutrals and plasma, initially coupled at the upper photosphere, become uncoupled at higher heights in the chromosphere. This decoupling can be a consequence of either the characteristic spatial scale at the shock front, that becomes similar to the collisional scale, or the change in the relation between the wave frequency, ion cyclotron frequency, and the collisional frequency with height. The decoupling height is a sensitive function of the wave frequency, wave amplitude, and the magnetic field strength. We observed that decoupling causes damping of waves and an increase in the background temperature due to the frictional heating. The comparison between analytical and numerical results allows us to separate the role of the nonlinear effects from the linear ones on the decoupling and damping of waves.

Solar chromosphere consists of a partially ionized plasma, which makes modeling the solar chromosphere a particularly challenging numerical task. Here we numerically model chromospheric waves using a two-fluid approach with a newly developed numerical code. The code solves two-fluid equations of conservation of mass, momentum and energy, together with the induction equation, for the case of the purely hydrogen plasma with collisional coupling between the charged and neutral fluid components. The implementation of a semi-implicit algorithm allows us to overcome the numerical stability constraints due to the stiff collisional terms. We test the code against analytical solutions of acoustic and Alfv\'en wave propagation in uniform medium in several regimes of collisional coupling.The results of our simulations are consistent with the analytical estimates, and with other results described in the literature. In the limit of a large collisional frequency, the waves propagate with a common speed of a single fluid. In the other limit of a vanishingly small collisional frequency, the Alfv\'en waves propagate with an Alfv\'en speed of the charged fluid only, while the perturbation in neutral fluid is very small. The acoustic waves in these limits propagate with the sound speed corresponding to either the charges or the neutrals, while the perturbation in the other fluid component is very small. Otherwise, when the collision frequency is similar to the real part of the wave frequency, the interaction between charges and neutrals through momentum transfer collisions cause alterations of the waves frequencies and damping of the wave amplitudes.

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.

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.

Magnetic reconnection in the solar corona is thought to be unstable to the formation of multiple interacting plasmoids, and previous studies have shown that plasmoid dynamics can trigger MHD waves of different modes propagating outward from the reconnection site. However, variations in plasma parameters and magnetic field strength in the vicinity of a coronal reconnection site may lead to wave reflection and mode conversion. In this paper we investigate the reflection and refraction of fast magnetoacoustic waves near a reconnection site. Under a justified assumption of an analytically specified Alfv\'{e}n speed profile, we derive and solve analytically the full wave equation governing propagation of fast mode waves in a non-uniform background plasma without recourse to the small-wavelength approximation. We show that the waves undergo reflection near the reconnection current sheet due to the Alfv\'en speed gradient and that the reflection efficiently depends on the plasma-$\beta$ parameter as well as on the wave frequency. In particular, we find that waves are reflected more efficiently near reconnection sites in a low-$\beta$ plasma which is typical for the solar coronal conditions. Also, the reflection is larger for lower frequency waves while high frequency waves propagate outward from the reconnection region almost without the reflection. We discuss the implications of efficient wave reflection near magnetic reconnection sites in strongly magnetized coronal plasma for particle acceleration, and also the effect this might have on First Ionization Potential (FIP) fractionation by the ponderomotive force of these waves in the chromosphere.

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 present the results of 3D-hydrodynamical simulations of accretion flow in the eclipsing dwarf nova V1239 Her in quiescence. The model includes the optical star filling its Roche lobe, a gas stream emanating from the inner Lagrangian point of the binary system, and the accretion disc structure. A cold hydrogen gas stream is initially emitted towards a point-like gravitational centre. A stationary accretion disc is formed in about 15 orbital periods after the beginning of accretion. The model takes into account partial ionization of hydrogen and uses realistic cooling function for hydrogen. The light curve of the system is calculated as the volume emission of optically thin layers along the line of sight up to the optical depth \tau =2/3 calculated using Planck-averaged opacities. The calculated eclipse light curves show good agreement with observations, with the changing shape of pre-eclipse and post-eclipse light curves being explained entirely due to ~ 50% variations in the mass accretion rate through the gas stream.

We review our understanding of ionized plasma and neutral gas coupling in the weakly ionized, stratified, electromagnetically-permeated regions of the Sun's chromosphere and Earth's ionosphere/thermosphere. Using representative models for each environment we derive fundamental descriptions of the coupling of the constituent parts to each other and to the electric and magnetic fields, and we examine the variation in magnetization of the ionized component. Using these descriptions we compare related phenomena in the two environments, and discuss electric currents, energy transfer and dissipation. We present a coupled theoretical and numerical study of plasma instabilities in the two environments that serves as an example of how the chromospheric and ionospheric communities can further collaborate. We also suggest future collaborative studies that will help improve our understanding of these two different atmospheres which share many similarities, but have large disparities in key quantities.

It has been proposed that particles bouncing between magnetized flows converging in a reconnection region can be accelerated by the first order Fermi mechanism. Analytical considerations of this mechanism have shown that the spectral index of accelerated particles is related to the total plasma compression within the reconnection region similarly to the case of diffusive shock acceleration mechanism. As a first step to investigate the efficiency of Fermi acceleration in reconnection regions in producing hard energy spectra of particles in the solar corona, we explore the degree of plasma compression that can be achieved at reconnection sites. In particular, we aim to determine the conditions for the strong compressions to form. Using a two-dimensional resistive MHD numerical model we consider a set of magnetic field configurations where magnetic reconnection can occur including a Harris current sheet, a force-free current sheet, and two merging flux ropes. Plasma parameters are taken to be characteristic of the solar corona. Numerical simulations show that strong plasma compressions ($\geq 4$) in the reconnection regions can form when the plasma heating due to reconnection is efficiently removed by fast thermal conduction or radiative cooling process. The radiative cooling process which is negligible in the typical 1 MK corona can play an important role in the low corona/transition region. It is found that plasma compression is expected to be strongest in low-beta plasma $\beta \sim 0.01-0.07$ at reconnection magnetic nulls.

Using RXTE/ASM archival data, we investigate the behaviour of the spectral hardness ratio as a function of X-ray luminosity in a sample of six transient X-ray pulsars (EXO 2030+375, GX 304-1, 4U 0115+63, V 0332+63, A 0535+26 and MXB 0656-072). In all sources we find that the spectral hardness ratio defined as $F_{5-12\mathrm{keV}}/ F_{1.33-3\mathrm{keV}}$ increases with the ASM flux (1.33--12 keV) at low luminosities and then saturates or even slightly decreases above some critical X-ray luminosity falling into the range $\sim(3-7)\times10^{37}$~erg~s$^{-1}$. Two-dimensional structure of accretion columns in the radiation-diffusion limit is calculated for two possible geometries (filled and hollow cylinder) for mass accretion rates $\dot M$ ranging from $10^{17}$ to 1.2$\times 10^{18}$~g s$^{-1}$. The observed spectral behaviour in the transient X-ray pulsars with increasing $\dot M$ can be reproduced by a Compton saturated sidewall emission from optically thick magnetized accretion columns with taking into account the emission reflected from the neutron star atmosphere. At $\dot M$ above some critical value $\dot M_{cr}\sim (6-8)\times 10^{17}$~g~s$^{-1}$, the hight of the column becomes such that the contribution of the reflected component to the total emission starts decreasing, which leads to the saturation and even slight decrease of the spectral hardness. Hollow-cylinder columns have a smaller height than the filled-cylinder ones, and the contribution of the reflected component in the total emission does not virtually change with $\dot M$ (and hence the hardness of the continuum monotonically increases) up to higher mass accretion rates than $\dot M_{cr}$ for the filled columns.

Realistic models of magnetic reconnection in the solar chromosphere must take into account that the plasma is partially ionized and that plasma conditions within any two magnetic flux bundles undergoing reconnection may not be the same. Asymmetric reconnection in the chromosphere may occur when newly emerged flux interacts with pre-existing, overlying flux. We present 2.5D simulations of asymmetric reconnection in weakly ionized, reacting plasmas where the magnetic field strengths, ion and neutral densities, and temperatures are different in each upstream region. The plasma and neutral components are evolved separately to allow non-equilibrium ionization. As in previous simulations of chromospheric reconnection, the current sheet thins to the scale of the neutral-ion mean free path and the ion and neutral outflows are strongly coupled. However, the ion and neutral inflows are asymmetrically decoupled. In cases with magnetic asymmetry, a net flow of neutrals through the current sheet from the weak field (high density) upstream region into the strong field upstream region results from a neutral pressure gradient. Consequently, neutrals dragged along with the outflow are more likely to originate from the weak field region. The Hall effect leads to the development of a characteristic quadrupole magnetic field modified by asymmetry, but the X-point geometry expected during Hall reconnection does not occur. All simulations show the development of plasmoids after an initial laminar phase.

Flux emergence is widely recognized to play an important role in the initiation of coronal mass ejections. The Chen-Shibata (2000) model, which addresses the connection between emerging flux and flux rope eruptions, can be implemented numerically to study how emerging flux through the photosphere can impact the eruption of a pre-existing coronal flux rope. The model's sensitivity to the initial conditions and reconnection micro-physics is investigated with a parameter study. In particular, we aim to understand the stability of the coronal flux rope in the context of X-point collapse and the effects of boundary driving in both unstratified and stratified atmospheres. In the absence of driving, we assess the behavior of waves in the vicinity of the X-point. With boundary driving applied, we study the effects of reconnection micro-physics and atmospheric stratification on the eruption. We find that the Chen-Shibata equilibrium can be unstable to an X-point collapse even in the absence of driving due to wave accumulation at the X-point. However, the equilibrium can be stabilized by reducing the compressibility of the plasma, which allows small-amplitude waves to pass through the X-point without accumulation. Simulations with the photospheric boundary driving evaluate the impact of reconnection micro-physics and atmospheric stratification on the resulting dynamics: we show the evolution of the system to be determined primarily by the structure of the global magnetic fields with little sensitivity to the micro-physics of magnetic reconnection; and in a stratified atmosphere, we identify a novel mechanism for producing quasi-periodic behavior at the reconnection site behind a rising flux rope as a possible explanation of similar phenomena observed in solar and stellar flares.

In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment. In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities. In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.

Magnetic reconnection in partially ionized plasmas is a ubiquitous phenomenon spanning the range from laboratory to intergalactic scales, yet it remains poorly understood and relatively little studied. Here, we present results from a self-consistent multi-fluid simulation of magnetic reconnection in a weakly ionized reacting plasma with a particular focus on the parameter regime of the solar chromosphere. The numerical model includes collisional transport, interaction and reactions between the species, and optically thin radiative losses. This model improves upon our previous work in Leake et al. 2012 \cite{Leake2012} by considering realistic chromospheric transport coefficients, and by solving a generalized Ohm's law that accounts for finite ion-inertia and electron-neutral drag. We find that during the two dimensional reconnection of a Harris current sheet with an initial width larger than the neutral-ion collisional coupling scale, the current sheet thins until its width becomes less than this coupling scale, and the neutral and ion fluids decouple upstream from the reconnection site. During this process of decoupling, we observe reconnection faster than the single-fluid Sweet-Parker prediction, with recombination and plasma outflow both playing a role in determining the reconnection rate. As the current sheet thins further and elongates it becomes unstable to the secondary tearing instability, and plasmoids are seen. The reconnection rate, outflows and plasmoids observed in this simulation provide evidence that magnetic reconnection in the chromosphere could be responsible for jet-like transient phenomena such as spicules and chromospheric jets.

We present results from the first self-consistent multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma. We simulate two dimensional magnetic reconnection in a Harris current sheet with a numerical model which includes ion-neutral scattering collisions, ionization, recombination, optically thin radiative loss, collisional heating, and thermal conduction. In the resulting tearing mode reconnection the neutral and ion fluids become decoupled upstream from the reconnection site, creating an excess of ions in the reconnection region and therefore an ionization imbalance. Ion recombination in the reconnection region, combined with Alfv\'{e}nic outflows, quickly removes ions from the reconnection site, leading to a fast reconnection rate independent of Lundquist number. In addition to allowing fast reconnection, we find that these non-equilibria partial ionization effects lead to the onset of the nonlinear secondary tearing instability at lower values of the Lundquist number than has been found in fully ionized plasmas.These simulations provide evidence that magnetic reconnection in the chromosphere could be responsible for jet-like transient phenomena such as spicules and chromospheric jets.

We present the science case for a broadband X-ray imager with high-resolution spectroscopy, including simulations of X-ray spectral diagnostics of both active regions and solar flares. This is part of a trilogy of white papers discussing science, instrument (Bandler et al. 2010), and missions (Bookbinder et al. 2010) to exploit major advances recently made in transition-edge sensor (TES) detector technology that enable resolution better than 2 eV in an array that can handle high count rates. Combined with a modest X-ray mirror, this instrument would combine arcsecondscale imaging with high-resolution spectra over a field of view sufficiently large for the study of active regions and flares, enabling a wide range of studies such as the detection of microheating in active regions, ion-resolved velocity flows, and the presence of non-thermal electrons in hot plasmas. It would also enable more direct comparisons between solar and stellar soft X-ray spectra, a waveband in which (unusually) we currently have much better stellar data than we do of the Sun.