Issue 
A&A
Volume 626, June 2019



Article Number  A46  
Number of page(s)  9  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201935326  
Published online  10 June 2019 
Intermediate shock substructures within a slowmode shock occurring in partially ionised plasma
University of Exeter, Exeter EX4 4QF, UK
email: b.snow@exeter.ac.uk
Received:
21
February
2019
Accepted:
28
April
2019
Context. Slowmode shocks are important in understanding fast magnetic reconnection, jet formation and heating in the solar atmosphere, and other astrophysical systems. The atmospheric conditions in the solar chromosphere allow both ionised and neutral particles to exist and interact. Under such conditions, fine substructures exist within slowmode shocks due to the decoupling and recoupling of the plasma and neutral species.
Aims. We study numerically the fine substructure within slowmode shocks in a partially ionised plasma, in particular, analysing the formation of an intermediate transition within the slowmode shock.
Methods. Highresolution 1D numerical simulations were performed using the (PIP) code using a twofluid approach.
Results. We discover that longlived intermediate (Alfvén) shocks can form within the slowmode shock, where there is a shock transition from above to below the Alfvén speed and a reversal of the magnetic field across the shock front. The collisional coupling provides frictional heating to the neutral fluid, resulting in a SedovTaylorlike expansion with overshoots in the neutral velocity and neutral density. The increase in density results in a decrease of the Alfvén speed and with this the plasma inflow is accelerated to above the Alfvén speed within the finite width of the shock leading to the intermediate transition. This process occurs for a wide range of physical parameters and an intermediate shock is present for all investigated values of plasmaβ, neutral fraction, and magnetic angle. As time advances the magnitude of the magnetic field reversal decreases since the neutral pressure cannot balance the Lorentz force. The intermediate shock is longlived enough to be considered a physical structure, independent of the initial conditions.
Conclusions. Intermediate shocks are a physical feature that can exist as shock substructure for long periods of time in partially ionised plasma due to collisional coupling between species.
Key words: magnetohydrodynamics (MHD) / shock waves / Sun: chromosphere
© ESO 2019
1. Introduction
Shocks occur readily in the lower solar atmosphere, driven by wave steepening, for example, umbral flashes (Beckers & Tallant 1969; Houston et al. 2018), or magnetic reconnection driven events, such as Ellerman bombs (Ellerman 1917; Nelson et al. 2013). The listed phenomena occur in regions of the sun where partially ionised effects are thought to play a key role in the underlying physics. As such, to understand phenomena occurring in the lower solar atmosphere, we must also understand the role of partial ionisation in shocks.
In magnetohydrodynamics (MHD), there are three characteristic wave speeds (slow, Alfvén and fast) which leads to a multitude of potential shock transitions (see for example Delmont & Keppens 2011). The slowmode shock (transition from superslow to subslow flow speeds) is of particular importance due to its role in magnetic reconnection, where a change in connectivity of the magnetic field results in a release of stored magnetic energy. The energy can be released instantaneously via Joule heating, or postreconnection by the influence of the magnetic field on the plasma. The classical 2D schematic for fast reconnection features slowmode shocks (Petschek 1964). Recent work suggests more complicated reconnection configurations occur in the solar atmosphere, for example, plasmoid instability, 3D topology, extended MHD (see reviews by Yamada et al. 2010; Pontin 2011; Loureiro & Uzdensky 2016; Cassak et al. 2017). In the Petschek model, the presence of slowmode shocks in the outflow region allow efficient transport of energy away from the reconnection region, increasing the reconnection rate. Slowmode shocks are also found to form in alternative models for fast reconnection (Liu et al. 2012; Innocenti et al. 2015; Shibayama et al. 2015).
The intermediate shock (transition from above to below the Alfvén velocity) has been shown to exist in singlefluid resistive MHD representative of the Earth’s magnetopause (Karimabadi 1995). The physical mechanism being that the nonideal region around the shock allows a separation of the magnetic field and fluid, and the formation of an intermediate shock between the two ideal regions. The stability of intermediate shocks was proved by Wu & Hada (1991) who studied numerically the formation of these shocks due to wave steepening. The relationship between intermediate shocks and slowmode shocks was shown analytically by Hau & Sonnerup (1989) for resistive MHD shocks, finding a relation between the downstream slowmode speed and the strength of the intermediate transition.
In partially ionised plasma, such as in the solar chromosphere or prominances, shocks become more complex. Hillier et al. (2016) studied the formation of slowmode shocks in partially ionised plasmas. Around the shock front, the plasma and neutral species decouple and recouple resulting in a finite width slowmode shock. The Lorentz force indirectly affects neutrals through collisions and hence the drift between ion and neutral species becomes an important parameter in partiallyionised plasma. However, observing the ionneutral drift has been found to be difficult due to lineofsight effects (Anan et al. 2017). A review of partiallyionised modelling can be found in Khomenko (2017).
In this paper, we have used highresolution twofluid numerical simulations to investigate the substructure within slowmode shocks. We discover that longlived intermediate shocks can form within the finite width of the shock as a result of the fluid coupling and decoupling around the shock front. This result has applications to the solar chromosphere, as well as interplanetary and interstellar partially ionised plasma. We studied the formation of compound intermediate shocks in partially ionised plasma, and the conditions in which such shock structures can form.
The outline of this paper is as follows. First we define the shock transitions and classifications to be used in this paper, and the analytical solution to the magnetohydrodynamic (MHD) and partially ionised plasma (PIP) equations in the Hoffman–Teller shock frame. Next, the numerical methods and initial conditions under consideration are introduced. Following this, a reference case is presented where we compare an MHD and a PIP simulation, and identify the key differences, investigating the physical evolution that produces the intermediate shock. Finally, we consider a parameter study and investigate the effect of the plasma properties on the lifetime and magnitude of the shock. We find that for our initial conditions, an intermediate shock will always form in a partiallyionised plasma. However the lifetime and magnitude depend heavily on the plasma properties.
2. Methodology
2.1. Shock classifications
MHD waves have three characteristic speeds: Alfvén (V_{A}), slow (V_{s}), and fast (V_{f}). As such, multiple shock transitions are possible, depending on the magnitude of the velocity, relative to the characteristic speeds in the pre and postshock regions.
Following the approach of Delmont & Keppens (2011), we classify shock transitions using the relationship between the normal flow velocity v_{⊥} and the characteristic speeds:

(1) superfast: V_{f} < v_{⊥},

(2) subfast: V_{A} < v_{⊥}< V_{f},

(3) superslow: V_{s} < v_{⊥}< V_{A},

(4) subslow: 0 < v_{⊥}< V_{s},

(∞) static: v_{⊥} = 0.
Defining the upstream condition u and downstream condition d, several shocks of the form u → d are possible:

1 → 2 fast shocks,

3 → 4 slow shocks,

1 → 2 = 3 switchon,

2 = 3 → 4 switchoff,

1 → 3, 1 → 4, 2 → 3, 2 → 4 intermediate shocks.
2.2. Hoffman–Teller equations
The Hoffman–Teller frame allows for jump relations to be derived from the MHD equations. In this choice of rest frame, the velocity and magnetic field vectors are in the same plane either side of the shock, that is, the electric field across the shock is zero. This reduces the MHD equations to a twodimensional problem for variables perpendicular (⊥) and parallel (∥) to the shock front.
2.2.1. MHD solution
In the Hoffmann–Teller frame, the MHD equations can be integrated and the upstream (^{u}) and downstream (^{d}) conditions can be equated:
We can also define the following: Alfvén Mach number as , plasmaβ as , and the angle between the velocity and magnetic field θ.
A solution to Eqs. (1)–(6) can be found in Hau & Sonnerup (1989) relating the upstream Alfvén velocity to the upstream plasmaβ, upstream θ, and downstream Alfvén velocity, that is,
This is shown graphically in Fig. 1 for a choice of plasmaβ^{u} = 0.1 and θ^{u} = π/4. The trivial solution is that the upstream and downstream velocities are identical, that is, no shock transition. The shaded region shows impossible solutions, where the velocity is higher downstream than upstream. The nontrivial solutions show the possible shock transitions for given plasmaβ and θ values. This curve intersects the A^{u} = A^{d} line at three points. The intersect labelled s^{u} denotes the upstream slow speed. At the point , a rotational discontinuity is possible, where there is a change in the angle of the magnetic field but plasma properties remain the same. Slowmode transitions (blue line) are bounded by . The downstream slow speed (point labelled on Fig. 1 as s^{d}) corresponds to the critical value that separates strong (2 → 4) and weak (2 → 3) intermediate transitions. The strong transition (2 → 4) is closely linked to the slowmode transition (3 → 4) through its relation to the downstream slow speed.
Fig. 1. Hau–Sonnerup shock solution for plasmaβ^{u} = 0.1 and θ^{u} = 0.4. Shaded region shows impossible transitions. Possible transitions include: slowmode shock (blue), intermediate shock (red), rotational discontinuity (circle), and fastmode shocks (green). 

Open with DEXTER 
2.2.2. PIP solution
In the PIP equations, collisional terms are taken into account and therefore, the equations have a nonzero righthand side, see Eqs. (8)–(17). As such, when the equations are integrated in the Hoffmann–Teller frame, there are integral terms on the right hand side, see Eqs. (A.1)–(A.13), as opposed to constants in the MHD Equations. However, by adding the neutral and plasma species together, the integral terms can be eliminated, and by choosing a point upstream and downstream such that the drift velocity equals zero, one recovers the MHD equations in the Hoffman–Teller frame, where the Alfvén speed and plasmaβ values depend on the total plasma density (ρ_{n} + ρ_{p}), see Appendix A. These equations are also independent from the neutral fraction and hence, in the steady state solution, the values either side of the shock are governed by the MHD solution, when points are chosen such that the drift velocity is zero. However, within the shock, the species decouple and recouple hence in partiallyionised plasmas, it is of interest to study the substructure that occurs within the finitewidth of the shock.
2.3. Numerical methods
Twofluid numerical simulations are performed using the (PIP) code (Hillier et al. 2016) which solves the interactions of a neutral fluid, and a coupled ion electron plasma. The simulations are 1D and use high resolution to resolve the substructure present within the shock. A first order HLLD (HartenLaxvan LeerDiscontiunities) scheme is used to prevent spurious oscillations from occurring around the shock interface. The normalised PIP equations are given below:
for neutral (subscript n) and plasma (subscript p) species. The neutral Eqs. (8)–(11) are independent of the magnetic field. The plasma Eqs. (12)–(17) are similar to the MHD equations however include collisional terms that couple the plasma to the neutral fluid. The collisional coefficient α_{c} is defined as
where T_{0} is the initial temperature and, in the normalised form, α_{0} = 1. The magnetic field is normalised using . Details of the equations and their implementation in the (PIP) code can be found in Hillier et al. (2016).
Our simulation data is translated into the shock frame by calculating the shock propagation speed in the MHD solution (v_{s}). The position and velocity variables are remapped to the shock frame, that is, v_{xns, xps} = v_{xn, xp} − v_{s}, x_{s} = x/t − v_{s}.
2.4. Initial conditions
The initial conditions in this paper are an extension of the model for slowmode shock created in reconnection proposed by Petschek (1964). Initially, a discontinuous magnetic field is specified across the boundary. The plasma and neutral fluids are assumed to be in thermal equilibrium at the start of the simulation.
where ξ_{n} and ξ_{i} are the neutral and ion fractions respectively.
The x = 0 boundary is treated as reflective, but such that magnetic field can penetrate the boundary. The system is normalised to have an Alfvén velocity of unity. The collisional time as determined by the bulk fluid density is calculated as τ = 1/(α(T_{0})ρ_{t}) = 1.
These initial and boundary conditions were used in the work of Hillier et al. (2016) to investigate the slowmode shock. Here we investigate the formation and lifetimes of a feature that was present in Figs. 4a and 5a of Hillier et al. (2016) but not discussed, namely the intermediate shock transition, which features a reversal of magnetic field across the shock interface.
3. Results
In this section, we present an MHD simulation and a PIP simulation and highlight the key differences. A snapshot of the results is shown in Fig. 2. Both simulations use 256 000 grid cells to resolve the spatial dimension. These two simulations use the following parameters: β = 1.0, B_{x} = 0.1, ξ_{n} = 0.9. We note that the neutral fraction is only used in the PIP case and the MHD simulation is fully ionised. The effect of different parameters is investigated in Sect. 4.2.
Fig. 2. MHD (black) and PIP (red neutral, blue plasma) reference solution for β = 1.0, B_{x} = 0.1, ξ_{n} = 0.9. Panel a: v_{x} velocity. Panel b: B_{y} magnetic field. Panel c: v_{y} velocity. Panel d: pressure. Green line indicates total pressure (P_{n} + P_{p}). The PIP solution is plotted after 2500 collisional times. 

Open with DEXTER 
3.1. Reference MHD solution
In the MHD case (Fig. 2, black lines), the initial conditions produce a rarefaction wave that drives fluid at the local Alfvén speed towards a slowmode shock. The solution is expanding in time, however, is steadystate when plotted on a timenormalised axis (x/t), hence can be described as pseudosteady.
It is important to note that in the MHD solution, the slowmode shock possesses no complex substructure, has a discontinuous jump (no finite width), and is purely a 3 → 4 transition. There is no reversal in magnetic field across the shock (see Fig. 2). These features do not hold when the model is extended to include additional effects (e.g. partial ionisation) where substructure can form in the shock front affecting the shock dynamics and resultant heating and energy transport.
3.2. Reference PIP solution
The red and blue lines in Fig. 2 show the neutral and ion fluid respectively for the PIP simulation, using the same parameters as the MHD case. The snapshot is taken after 2500 collisional times when the system can be assumed to be in a quasiselfsimilar state. In the PIP case, the ion and neutral species decouple and recouple around the shock front, resulting in a finite width slowmode shock where shock substructure can occur (compared to the discontinuous jump in the MHD case).
The initial conditions are a discontinuity in the magnetic field only, therefore directly affect the plasma only. The neutral fluid can only be influenced indirectly by collisions with the ionised plasma. As such, the system tends towards a pseudosteadystate as time tends towards infinity. Hillier et al. (2016) presented the dynamic evolution for a similar case:
Initialisation. The initial conditions are a discontinuity in the magnetic field and hence directly affect the plasma only, and the neutral fluid is affected via collisions. At τ = 1 the system is highly decoupled. A SedovTaylorlike expansion occurs in the neutral fluid due to frictional heating between the two species.
Weak coupling. A fastmode rarefaction wave forms and drives fluid towards a slowmode shock. At this time, the rarefaction wave and slowmode shock front are coupled and interact.
Strong coupling. In this phase, the rarefaction wave and slowmode shock decouple. The slowmode shock can be considered to be independently evolving.
Quasiselfsimilar state. System tends towards a selfsimilar state as time tend towards infinity. This is plotted in Fig. 2. The rarefaction wave is separated from the slowmode shock, Fig. 2a. The two fluids are reasonably well coupled except within the slowmode shock where large drift velocities are present, Fig. 2a. The neutral overshoot is present at late times indicating that it is a stable feature of the system (Fig. 2a).
A feature present but not discussed in the work of Hillier et al. (2016) is the reversal of the magnetic field across the shock front (Fig. 2b). This is a longlived (but transient) feature of the system that exists long after the rarefaction wave and slowmode shock have separated. The magnetic field reversal is a signature that an intermediate transition exists within the slowmode shock.
4. Analysis
4.1. Intermediate shock
The reversal in magnetic field in Fig. 2b is a key indicator of an intermediate shock, where the plasma transitions from above the local Alfvén speed, to below it. An intermediate shock is present in all partially ionised simulations as a transient feature of the system, however is not present in the MHD cases. The feature arises due to the interactions of ion and neutral particles. To analyse this feature, a new frame of reference is chosen such that the shock is stationary, that is, x_{s} = x/t − v_{s} where v_{s} is the propagation speed of the MHD shock. In this frame, the velocity is also adjusted to account for the shock propagation speed, that is, v_{x} − v_{s}. The shock frame allows us to correctly compute the transitions across the shock front and compare with the analytical results. The evolution of this structure is as follows:
1. SedovTaylorlike expansion of the neutral fluid. The initial conditions are a discontinuity in the magnetic field that affects the plasma only, resulting in a separation of neutral and plasma species in the shock front (Fig. 3a). The large drift velocity creates a rapid heating of the neutral fluid. Subsequently, a hydrodynamic, SedovTaylorlike expansion occurs in the neutral fluid, resulting in an overshoot of neutral velocity (Fig. 3a), and a sudden increase in neutral density (Fig. 3b).
Fig. 3. MHD (black dashed) and PIP (red neutral, blue plasma) solutions within the shock. The figures are in the shock frame (x_{s} = x/t − v_{s}, where v_{s} is the propagation speed of the shock in the MHD solution). Panel a: perpendicular velocity in the shock frame v_{x} − v_{s}. Panel b: density. Green shows the total (ρ_{n} + ρ_{p}) density. Panel c: Alfvén Mach number using the plasma density (magenta) and total density (green). Panel d: Alfvén speed using the total density. Panel e: shock transitions. Panel f: pressure. Green line shows total pressure (P_{n} + P_{p}). 

Open with DEXTER 
2. Acceleration of plasma. The rarefaction wave drives velocity towards the slowmode shock at the Alfvén speed (Fig. 3c). This is also true in the MHD case and is a feature of the shock problem being studied. Within the slowmode shock front, the bulk Alfvén speed decreases due to the increase in neutral density (Fig. 3d) and the reversal in the magnetic field (Fig. 2b). The plasma inside the shock is accelerated to above the Alfvén speed (Fig. 3c), resulting in an intermediate transition (Fig. 3e).
3. Decay with time. The intermediate shock is sustained by the neutral pressure (Fig. 3f), which, because of the finite coupling, cannot balance the Lorentz force. As time tends to infinity, the neutral pressure equalises and the magnetic field reversal becomes decreases in magnitude, see Fig. 4. This is an indicator that the intermediate shock disappears with time. For these parameters, the magnetic field reversal is still prominent after 10 000 collisional times, indicating that whilst this is a transient feature of the system, it is sufficiently longlived to be physical and independent of the initial conditions.
Fig. 4. Maximum neutral pressure (dashed) and maximum magnetic field reversal (solid) through time. Both quantities are maximum near the intermediate transition. 

Open with DEXTER 
The frictional heating (α_{c}(T_{n}, T_{p})ρ_{n}ρ_{p}(v_{n} − v_{p})^{2}) and thermal damping (3α_{c}(T_{n}, T_{p})ρ_{n}ρ_{p}(P_{n}/ρ_{n} − P_{p}/(2ρ_{p}))) between the two species are shown in Fig. 5 through the finitewidth of the shock. Both the frictional heating and thermal damping reach their maximum values at the intermediate transition. Here, there is the largest drift velocity (Fig. 3a). Either side of the shock, the species are reasonably well coupled and hence there is minimal frictional heating.
Fig. 5. Frictional heating (black line) and thermal damping (red and blue lines). The red line indicates ions losing heat to neutrals, and vice versa for the blue line. Dashed lines indicate the finite width of the shock. 

Open with DEXTER 
4.2. Parameter study
4.2.1. Changing ionisation fraction
The simulations in this section use β = 0.1, B_{x} = 0.3 with the ionisation fraction ξ_{i} varying. The results can be scaled by a rough estimate of the time scale changes by multiplying by the ionisation fraction, that is, tξ_{i}. As the neutral fraction increases, the influence of the neutral pressure becomes more important. The initial evolution stages generate an overshoot in the neutral pressure due to the collisional coupling. There is then an interplay between the neutral pressure and the magnetic tension. When ξ_{n} ≤ 0.9, the neutral pressure increases with time to balance the magnetic tension. When ξ_{n} > 0.9 the neutral pressure overshoot decreases with time. This is shown in Fig. 6 with the ξ_{n} = 0.9 (black dashed) and the ξ_{n} = 0.99999 (magenta dashed) lines. Whilst the two results look very different, the same mechanism is occurring whereby the system seeks an equilibrium; the main difference is the magnitude of the neutral pressure overshoot due to the initial heating. All cases tend towards a constant maximum neutral pressure.
Fig. 6. Maximum neutral pressure (dashed) and magnetic field reversal (solid) for different neutral fractions ξ_{n} = 0.9 (black), 0.99 (red), 0.999 (blue), 0.9999 (green), 0.99999 (magenta). 

Open with DEXTER 
All neutral fractions investigated feature the intermediate transition as a substructure within the slow mode shock. The normalised time scale does not map the behaviour perfectly and there are differences in magnitude and gradients of the magnetic field across the different neutral fractions. The neutral pressure is the key variable due to its role in the equalising the Lorentz force. As the neutral fraction tends to unity, the maximum neutral pressure is reasonably similar across simulations.
4.2.2. Changing beta
Figure 7 shows the neutral pressure and magnetic field reversal for the lowbeta regime. All three cases plotted contain roughly the same gradient in magnetic field over time but have different equilibrium neutral pressures. As the plasmaβ increases, the maximum equilibrium neutral pressure increases drastically. This results in a larger restoring force and the magnetic field decreases more rapidly.
Fig. 7. Neutral pressure (dashed) and magnetic field reversal (solid) for different beta values β = 0.1 (black), 0.01 (red), and 0.001 (blue). 

Open with DEXTER 
For large plasmaβ values, a sonic shock can occur in the neutral fluid, as seen in Hillier et al. (2016). In these simulations, there can be two shocks occurring within the finitewidth shock region: intermediate shock in the plasma, and sonic shock in the neutrals.
4.2.3. Changing magnetic angle
The initial B_{x} component of magnetic field is modified to alter the angle of the magnetic field. Note that the initial B_{y} component remains the same and hence the total B strength varies across simulations. Interestingly, the preshock region has B≈1 in all simulations as a consequence of the investigated system.
Changing the B_{x} value changes the propagation speeds of waves and hence phenomena occurs on different time scales. From Fig. 8 it is clear that a larger B_{x} results in faster equalisation of the neutral pressure and hence the magnetic field reversal and the intermediate shock decay faster, compare to low B_{x} values. Magnetic tension is the equalising force for the intermediate shock. The magnetic tension increases for larger B_{x} values and hence the intermediate shock decays faster.
Fig. 8. Neutral pressure (dashed) and magnetic field reversal (solid) for B_{x} = 0.1 (black), 0.2 (red), 0.3 (blue), 0.6 (green). 

Open with DEXTER 
5. Comments on potential observations
In the MHD case, the fluid and the magnetic field are frozen together, so in the rest frame of the fluid the electric field becomes E ∝ v × B = 0. However, for the two fluid case, there is a separation of the species within the finitewidth of the shock. Within the finitewidth, even though there is no electric field felt by the plasma fluid there is therefore an electric field felt by the neutrals. This can be calculated by looking at how fast the neutrals move across the magnetic field, which is the drift velocity, that is, E_{n} ∝ (v_{n} − v_{p})×B. This metric is zero either side of the shock since the species are fully coupled. Inside the shock, the ions and neutrals decouple and hence there is a localised increase in the neutral electric field inside the shock substructure. The electric field felt by the neutrals is a potential observable for measuring the neutralion drift (e.g. Anan et al. 2014).
The maximum electric field felt by the neutrals within the shock is plotted through time in Fig. 9 for the reference PIP simulation (β = 1, B_{x} = 0.1, ξ = 0.9). As time advances, this tends towards a constant value indicating that there will always be a localised increase in neutral electric filed within the shock. The presence of this enhanced neutral electric field could be a potential observable for the effects of partial ionisation in shock waves.
Fig. 9. Peak neutral electric field increase within the shock through time. Note that the electric field outside the shock in this frame is zero. 

Open with DEXTER 
6. Summary
This paper has demonstrated that intermediate shocks can occur as substructure inside slowmode shocks due to partial ionisation. Highresolution 1D numerical simulations were performed to fully resolve the substructure that occurs within the finite width of the slowmode shock. The physical process involved in forming and dissipating the intermediate shock substructure is as follows:
1. Collisional coupling results in overshoots in neutral velocity and density. The initial conditions drive the plasma only. As such, there is a large drift velocity between the two species. The neutral fluids response is to create a SedovTaylorlike expansion which creates localised increases in the neutral density and velocity at the interface.
2. Acceleration of the plasma. Our system has velocity driven towards the slowmode shock at the Alfvén speed. Within the finitewidth of the shock, the Alfvén velocity decreases due to the increased neutral pressure and reversal of the magnetic field, and hence the plasma velocity is accelerated inside the shock to above the Alfvén speed.
3. Decay with time. The neutral pressure is not sufficient to balance the Lorentz force and hence there is a gradual evolution whereby the magnetic field reversal tends to zero as time tends towards infinity. The intermediate shock is however sufficiently longlived to be independent of the initial conditions and considered a physical feature of partially ionised shocks.
The intermediate shock was present as substructure for all tested parameter regimes but the magnitude of the magnetic field reversal was parameter dependent. The larger the equilibrium neutral pressure, the larger the magnetic field reversal. As such, the intermediate shock is strongest for high plasmaβ values. We would therefore expect this feature to be most significant in the lower atmosphere, for example, Ellerman bombs.
The work in this paper is analogous to the magnetospheric intermediate shocks discussed in Karimabadi (1995) where the resistivity creates a dispersive region accelerating the plasma towards the shock region and resulting in an intermediate transition. Here we have a similar effect except it is driven by twofluid interaction in the absence of resistivity.
We have shown that an intermediate shock exists for a wide range of parameters, with varying degrees of magnitude and lifetime. Theoretically, this has the potential to be present across a wide range of phenomena in the solar atmosphere where partialionisation effects are important, from wavesteepening events (e.g. umbral flashes), to magnetic reconnection (Ellerman bombs, spicules). Future work will be to analyse the implications and observability of intermediate shocks in the lower solar atmosphere.
In summary, there are four main conclusions in this paper:
1. Ideal MHD heating across the shock. We have shown analytically that the twofluid equations reduce to the ideal MHD shock equations when the species are coupled either side of the shock (see Appendix A). The consequence of this is that one would expect to obtain ideal MHDlike heating across the shock, with no additional heating from the collisions.
2. Shocks as substructure in PIP case. Within the finite width of a partiallyionied shock, interactions between the two species can lead to the formation of intermediate transitions within the larger shock structure. Intermediate shocks are a transition from above to below the Alfvén speed (here superAlfvén to subslow) and feature a reversal in the magnetic field (see Fig. 2). This feature forms due to the collisional effects between the two species inside the finitewidth shock. An intermediate transition was present for all tested parameter regimes.
3. Potential for large currents. The formation of an intermediate shock features a sharp reversal in the magnetic field across the intermediate shock front, hence there is the potential for large currents to form inside the largescale slowmode transition. Large currents are known to play a role in particle acceleration. Hence the formation of intermediate shocks in partially ionised plasma may lead to an additional particle acceleration mechanism.
4. Localised electric field experienced by the neutrals. Within the finite width of the shock, there is a localised increase in neutral electric field. This tends towards a constant value as time tends to infinity, hence may be a potential observable of ionneutral interactions within shock fronts.
Acknowledgments
BS and AH are supported by STFC research grant ST/R000891/1. AH is also supported by STFC Ernest Rutherford Fellowship grant number ST/L00397X/2.
References
 Anan, T., Ichimoto, K., & Hillier, A. 2017, A&A, 601, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Anan, T., Casini, R., & Ichimoto, K. 2014, ApJ, 786, 94 [NASA ADS] [CrossRef] [Google Scholar]
 Beckers, J. M., & Tallant, P. E. 1969, Sol. Phys., 7, 351 [NASA ADS] [CrossRef] [Google Scholar]
 Cassak, P. A., Liu, Y.H., & Shay, M. A. 2017, J. Plasma Phys., 83, 715830501 [CrossRef] [Google Scholar]
 Delmont, P., & Keppens, R. 2011, J. Plasma Phys., 77, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Ellerman, F. 1917, ApJ, 46, 298 [NASA ADS] [CrossRef] [Google Scholar]
 Hau, L.N., & Sonnerup, B. U. O. 1989, J. Geophys. Res., 94, 6539 [NASA ADS] [CrossRef] [Google Scholar]
 Hillier, A., Takasao, S., & Nakamura, N. 2016, A&A, 591, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Houston, S. J., Jess, D. B., Asensio Ramos, A., et al. 2018, ApJ, 860, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Innocenti, M. E., Goldman, M., Newman, D., Markidis, S., & Lapenta, G. 2015, ApJ, 810, L19 [NASA ADS] [CrossRef] [Google Scholar]
 Karimabadi, H. 1995, Adv. Space Res., 15, 507 [NASA ADS] [CrossRef] [Google Scholar]
 Khomenko, E. 2017, Plasma Phys. Control. Fusion, 59, 014038 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, Y.H., Drake, J. F., & Swisdak, M. 2012, Phys. Plasmas, 19, 022110 [NASA ADS] [CrossRef] [Google Scholar]
 Loureiro, N. F., & Uzdensky, D. A. 2016, Plasma Phys. Control. Fusion, 58, 014021 [NASA ADS] [CrossRef] [Google Scholar]
 Nelson, C. J., Shelyag, S., Mathioudakis, M., et al. 2013, ApJ, 779, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Petschek, H. E. 1964, NASA Spec. Publ., 50, 425 [NASA ADS] [Google Scholar]
 Pontin, D. I. 2011, Adv. Space Res., 47, 1508 [NASA ADS] [CrossRef] [Google Scholar]
 Shibayama, T., Kusano, K., Miyoshi, T., Nakabou, T., & Vekstein, G. 2015, Phys. Plasmas, 22, 100706 [NASA ADS] [CrossRef] [Google Scholar]
 Wu, C. C., & Hada, T. 1991, J. Geophys. Res., 96, 3769 [NASA ADS] [CrossRef] [Google Scholar]
 Yamada, M., Kulsrud, R., & Ji, H. 2010, Rev. Mod. Phys., 82, 603 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Hoffman–Teller PIP equations
In the Hoffman–Teller frame, the twofluid PIP equations are:
The integral terms (I_{1}, I_{2}, I_{3}) can be removed by adding the neutral and ion equations together, that is,
The partial pressure and density can be expressed in terms of a total value using the neutral fraction ξ_{n}:
Substituting these into Eqs. (A.14)–(A.19) gives:
Furthermore, we can impose an additional constraint such that either side of the shock the drift velocity equals zero (v_{∥p} = v_{∥n} = v_{∥} and v_{⊥p} = v_{⊥n} = v_{⊥}):
Equations (A.30)–(A.35) are identical to the MHD equations except here the density is the total density ρ_{t} and pressure P_{t}. Therefore, the solution to these equations is identical to the Hau & Sonnerup (1989) solution for MHD (Eq. (7)), independent of the neutral fraction. It should be noted that this is only true over the larger shock structure and inside the shock there is substructure that is highly dependent on the collisional effects.
All Figures
Fig. 1. Hau–Sonnerup shock solution for plasmaβ^{u} = 0.1 and θ^{u} = 0.4. Shaded region shows impossible transitions. Possible transitions include: slowmode shock (blue), intermediate shock (red), rotational discontinuity (circle), and fastmode shocks (green). 

Open with DEXTER  
In the text 
Fig. 2. MHD (black) and PIP (red neutral, blue plasma) reference solution for β = 1.0, B_{x} = 0.1, ξ_{n} = 0.9. Panel a: v_{x} velocity. Panel b: B_{y} magnetic field. Panel c: v_{y} velocity. Panel d: pressure. Green line indicates total pressure (P_{n} + P_{p}). The PIP solution is plotted after 2500 collisional times. 

Open with DEXTER  
In the text 
Fig. 3. MHD (black dashed) and PIP (red neutral, blue plasma) solutions within the shock. The figures are in the shock frame (x_{s} = x/t − v_{s}, where v_{s} is the propagation speed of the shock in the MHD solution). Panel a: perpendicular velocity in the shock frame v_{x} − v_{s}. Panel b: density. Green shows the total (ρ_{n} + ρ_{p}) density. Panel c: Alfvén Mach number using the plasma density (magenta) and total density (green). Panel d: Alfvén speed using the total density. Panel e: shock transitions. Panel f: pressure. Green line shows total pressure (P_{n} + P_{p}). 

Open with DEXTER  
In the text 
Fig. 4. Maximum neutral pressure (dashed) and maximum magnetic field reversal (solid) through time. Both quantities are maximum near the intermediate transition. 

Open with DEXTER  
In the text 
Fig. 5. Frictional heating (black line) and thermal damping (red and blue lines). The red line indicates ions losing heat to neutrals, and vice versa for the blue line. Dashed lines indicate the finite width of the shock. 

Open with DEXTER  
In the text 
Fig. 6. Maximum neutral pressure (dashed) and magnetic field reversal (solid) for different neutral fractions ξ_{n} = 0.9 (black), 0.99 (red), 0.999 (blue), 0.9999 (green), 0.99999 (magenta). 

Open with DEXTER  
In the text 
Fig. 7. Neutral pressure (dashed) and magnetic field reversal (solid) for different beta values β = 0.1 (black), 0.01 (red), and 0.001 (blue). 

Open with DEXTER  
In the text 
Fig. 8. Neutral pressure (dashed) and magnetic field reversal (solid) for B_{x} = 0.1 (black), 0.2 (red), 0.3 (blue), 0.6 (green). 

Open with DEXTER  
In the text 
Fig. 9. Peak neutral electric field increase within the shock through time. Note that the electric field outside the shock in this frame is zero. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.