Magnetohydrodynamic Waves in an Asymmetric Magnetic Slab

Analytical models of solar atmospheric magnetic structures have been crucial for our understanding of magnetohydrodynamic (MHD) wave behaviour and in the development of the field of solar magneto-seismology. Here, an analytical approach is used to derive the dispersion relation for MHD waves in a magnetic slab of homogeneous plasma enclosed on its two sides by non-magnetic, semi-infinite plasma with different densities and temperatures. This generalises the classic magnetic slab model, which is symmetric about the slab. The dispersion relation, unlike that governing a symmetric slab, cannot be decoupled into the well-known sausage and kink modes, i.e. the modes have mixed properties. The eigenmodes of an asymmetric magnetic slab are better labelled as quasi-sausage and quasi-kink modes. Given that the solar atmosphere is highly inhomogeneous, this has implications for MHD mode identification in a range of solar structures. A parametric analysis of how the mode properties (in particular the phase speed, eigenfrequencies, and amplitudes) vary in terms of the introduced asymmetry is conducted. In particular, avoided crossings occur between quasi-sausage and quasi-kink surface modes, allowing modes to adopt different properties for different parameters in the external region.


Introduction
Dynamic solar events have been widely observed to induce perturbations in the magnetically dominated coronal plasma (Banerjee et al., 2007;McLaughlin, Hood, and de Moortel, 2011;Arregui, Oliver, and Ballester, 2012;Mathioudakis, Jess, and Erdélyi, 2013;Komm et al., 2015). The inhomogeneous structuring of the plasma parameters determines the characteristics and speed of the waves propagating in these structures. The detection and analysis B R. Erdélyi robertus@sheffield.ac.uk 1 Solar Physics and Space Plasma Research Centre, School of Mathematics and Statistics, University of Sheffield, Hicks Building, Hounsfield Road, Sheffield, S3 7RH, UK of waves provides an indirect way to gain information about the structure of the magnetic waveguide environment through the techniques of solar magneto-seismology, see e.g. reviews by Nakariakov and Verwichte (2005), Andries et al. (2009), Erdélyi (2009), De Moortel andNakariakov (2012). Analytical models of waves in solar magnetic environments are crucial in the accurate employment of this technique.
Linear MHD waves that propagate along magnetic field and density stratifications have been widely researched in Cartesian geometry to gain first insight into the often complex magnetic structures, including e.g. wave propagation in sunspots, coronal loops, prominences, or coronal hole boundaries, to name a few applications. One of the simplest approaches is when a single planar magnetic and/or density interface is considered and the properties of the available surface waves that propagate along this interface are deduced (Roberts, 1981a). Adding a second planar interface to this model, in order to form a simple way of mimicking structuring, gives a magnetic slab where the external environment is not magnetised and has a uniform but different density from the internal slab density (Roberts, 1981b). Next, a homogeneous magnetic field may then be added to the external plasma (Edwin and Roberts, 1982), and the eigenvalues, i.e. allowed wave modes of a magnetised slab system, can be determined. After this initial surge in analytical linear models of solar atmospheric magnetic structures, Cartesian geometry was replaced in favour of cylindrical geometry, for which significant research has led to the development and applications of solar magneto-seismology (Rosenberg, 1970;Uchida, 1970;Zajtsev and Stepanov, 1975;Roberts, Edwin, and Benz, 1984;Goossens, Andries, and Aschwanden, 2002; see a review by Ruderman and Erdélyi, 2009) and hypothesised mechanisms for solar plasma heating including mechanisms such as phase mixing (Heyvaerts and Priest, 1983) and resonant absorption (see a review by Goossens, Erdélyi, and Ruderman, 2011, with plenty of references), to name a few. This trend is largely a result of ubiquitous observations of solar atmospheric magnetic flux tube oscillations, which lend themselves more naturally to cylindrical geometry.
The purpose of the present work is to generalise the isolated magnetic slab model by investigating the linear wave physics that arises when the density and temperature of one external plasma is different from that of the other, and external densities and temperatures are both different from those inside the slab. The interest in this generalisation comes from the asymmetry of the system, which gives rise to asymmetric quasi-sausage and quasi-kink eigenmodes that demonstrate mixed properties.
Magnetic structures in the corona where solar magneto-seismology has previously been employed (Nakariakov and Ofman, 2001;Nakariakov and Verwichte, 2005) are often better modelled by cylindrical than Cartesian geometry. Instead, structures closer to the photosphere, such as the magnetic canopy, provide an asymmetric slab structure for application of this work. The magnetic canopy is a region of dominant magnetic field parallel to the surface of the Sun between the much less magnetised photosphere and the chromosphere that connects the magnetic field lines between active regions. Another application may be to oscillations in magnetic bright points (MBPs). MBPs are often elongated vertical magnetic structures between granular cells of different densities and temperatures. Further information on lower atmospheric solar magneto-seismology can be found in e.g. de Pontieu and Erdélyi (2006).
First, a derivation of the dispersion relation is presented in Section 2. The eigenmode solutions of the dispersion relation are compared to those of a symmetric slab in Section 2.1. The properties of the waves in an incompressible and low-beta plasma are discussed in Sections 3.2 and 3.3. Wide and thin slab approximations are made in Sections 3.4 and 3.5. A numerical procedure is employed in Section 4 to investigate the effect that varying the external densities has on wave dispersion.

Figure 1
Equilibrium state inside the slab (|x| ≤ x 0 ) and outside the slab (x < −x 0 and x > x 0 ). The red arrows illustrate the vertical magnetic field, B(x) z, and the dashed black lines indicate the boundaries of the slab.

Derivation of the Dispersion Relation
Consider an unbounded inviscid static plasma under an equilibrium magnetic field, B(x) z, where z is the unit vector in the vertical direction, and where Within each region, the plasma is uniform and the equilibrium plasma pressure, density, and temperature are denoted by p i , ρ i , and T i , respectively, for i = 0, 1, 2 (see Figure 1). The effects of gravity are ignored throughout; it is important to note, however, that equilibrium density stratification itself in the solar atmosphere can be a consequence of gravity, but gravity can be ignored if the gravity scale height is large compared to the wavelength and the thickness of the magnetic slab, which is safe to assume for many small-scale solar atmospheric structures. For the equilibrium conditions to be stable, there must be equilibrium pressure balance across each interface: Here, μ 0 is the permeability of free space. The sound speed in each region is denoted by c i = √ γ p i /ρ i for i = 0, 1, 2, where γ is the adiabatic index, which can be taken as constant across the whole domain if we take a single-fluid approximation in each region. Equilibrium pressure balance, Equation (2), gives the relationship between the external sound speeds and densities: That is, to satisfy the equilibrium structure, the external densities and sound speeds must obey this relation. The governing equations for the disturbance within the magnetic slab are the ideal MHD equations: d dt where the variables v = (v x , v y , v z ), B, p, and ρ are the velocity, magnetic field, pressure, and density, at time t . After linearising about a static basic state, we can combine the governing equations and seek solutions of the form where k is the component of the wavenumber vector in the z-direction and ω is the angular frequency. This restricts the investigation to waves propagating parallel to the equilibrium magnetic field, with velocity perturbation amplitude v x (x) in the x-direction, and v z (x) in the z-direction. With this ansatz, the linearised equations reduce to where v A = B 0 / √ ρ 0 μ is the Alfvén speed and = d/dx. These equations can be combined to give an ordinary differential equation where which is valid for perturbations inside the slab. This is identical to the corresponding equation for a symmetric slab derived by Roberts (1981b). By setting v A = 0, a change of subscripts in Equation (11) is all that is required to derive the ordinary differential equation governing linear perturbations outside the slab, namely where The solutions of Equations (11) and (13) are a linear combination of the hyperbolic functions. We restrict our model to waves trapped by the slab by ensuring v x → 0 as |x| → ∞. The general solution for the velocity perturbation in the x-direction is where A, B, C, and D are arbitrary constants (with respect to x). The plasma pressure perturbation can be assumed to be of the form p(x, t) = p(x)e i (kz−ωt) , then the total pressure perturbation amplitude across the whole domain is where The boundary conditions are continuity of velocity and total pressure across the slab boundaries at x = ±x 0 , which gives the four coupled homogeneous algebraic equations ⎛ where c i = cosh m i x 0 and s i = sinh m i x 0 for i = 0, 1, 2, for brevity. The condition for the existence of non-trivial solutions to this system of equations is that the determinant of the above matrix of coefficients is zero. Applying this condition gives us the dispersion relation for an asymmetric slab, namely Using the notation introduced by Equation (14), we arrive at

Comparison with a Symmetric Slab
There is an intrinsic difference between perturbations along symmetric and asymmetric magnetic slabs. The dispersion relation governing an asymmetric slab is a single equation, whereas the dispersion relation governing a symmetric slab (Roberts, 1981a) consists of two independent equations, corresponding to the sausage and kink eigenmodes. Under the approximation that the densities and temperatures of the external plasma are of the same order with respect to each other, the dispersion relation, Equation (20), can be factorised to give the approximate dispersion relation The expressions for the variables Λ i for i = 0, 1, 2 in Equations (17) can be employed to yield the approximately symmetric dispersion relation This equation is now in an analogous form to the dispersion relation corresponding to MHD waves along a symmetric magnetic slab (Roberts, 1981b), namely where external parameters are denoted by subscript e.

Asymmetric Eigenmodes
There is a rich spectrum of MHD waves supported by a magnetic slab. In a symmetric slab, the dispersion relation, Equation (23), consists of two decoupled equations that correspond to the two types of fundamental wave supported by the slab: the "sausage" and "kink" MHD waves. Sausage and kink modes have been observed, for example, in chromospheric fibrils, where the ubiquity of these waves has been linked to coronal heating (Morton et al., 2012). Sausage and kink modes can be further categorised into "surface" and "body" modes. Surface modes are waves more enhanced at the slab boundaries, whereas body waves are characterised by oscillations permeating spatially throughout the slab, having their maximum amplitude within the slab. Mathematically, surface waves correspond to exponential solutions of Equation (11), that is, they exist when m 2 0 > 0, which occurs when the phase speed ω/k satisfies Body waves correspond to spatially oscillatory solutions (i.e. they can have nodes inside the slab), that is, they exist when m 2 0 < 0, which occurs when For an asymmetric slab, the sausage and kink modes are modified by the external density difference, causing an asymmetry of the oscillation amplitude on each side of the slab (for visualisation see Figures 2a and 2b). We call these asymmetric eigenmodes quasi-sausage and quasi-kink modes. In a symmetric slab, sausage modes are characterised by a line of zero perturbation at the centre of the slab. In an asymmetric slab, this line is shifted towards the side of greatest external density. For symmetric kink modes, the width of the perturbed slab remains constant along the slab, but this characteristic is dropped in an asymmetric slab. This highlights the mixed nature of these modes.
The surface and body properties also take a modified form in an asymmetric slab. For a surface mode (for visualisation see Figures 3a and 3b), the wave power distribution across the slab has a single minimum. The displacement of this minimum from the centre of the slab is a consequence of the asymmetry in the external plasma. The intensity of the maximum amplitudes on the left and right boundaries of the slab (x = ±x 0 ) is different, reflecting the asymmetry in the external plasma.
Body modes are also affected by the asymmetric external environment (for visualisation see Figures 3c and 3d). Local maxima and minima in wave power are shifted towards the external plasma of higher density for a quasi-kink body mode and towards the external plasma of lower density for a quasi-sausage mode. However, body modes depend only weakly on the external plasma parameters and are therefore less affected than surface modes. This is shown analytically in Sections 3.4 and 3.5.

Figure 2
Quasi-kink and quasi-sausage modes with external density ordering ρ 1 > ρ 2 . The red lines illustrate the perturbed magnetic field, the thick solid black lines illustrate the perturbed slab boundaries, and the dashed lines illustrate the future position of the slab boundaries after half a period.

Figure 3
The transverse velocity perturbation amplitude, v x as a function of the transverse spatial coordinate, x, for quasi-sausage and quasi-kink modes in an isolated magnetic slab with external density ordering ρ 1 > ρ 2 .

Analytical Solutions
In this section, simplifications are made to the dispersion relation, Equation (20), and approximate dispersion relation, Equation (22). Thin slab, wide slab, low-beta, and incompressible approximations are made.

Spurious Solutions
The solution to the exact and approximate dispersion relations, Equations (20) and (22), given by ω = kv A , is spurious and does not correspond to an eigenmode. To see this, observe that m 0 = 0 for this solution, which leads to a linear (rather than oscillatory or exponential) solution to the governing differential equation, Equation (2.20). The same can be said for the solutions ω = kc 0 and ω = kc T . This rules out the possibility of the existence of pure sound waves and pure Alfvén waves.

Incompressible Approximation
Compressibility is essential for the propagation of sound waves. Consider the dispersion relation, Equation (20), in the limit of incompressibility, that is, when γ → ∞, where γ is the adiabatic index. In this limit, the sound speeds become unbounded, therefore the tube speed behaves like c T → v A . This means that m j → k for j = 0, 1, 2, so that the dispersion relation reduces to This is a special case of the dispersion relation previously derived by Ruderman (1992), who found solitons propagating on a system of N tangential discontinuities. This is a quadratic equation in ω 2 , with solutions given by where σ = 1 2 These solutions hold for all kx 0 and describe surface modes with sub-Alfvénic phase speed. Figures 4a -4d demonstrate that in a thin incompressible slab, the phase speeds of these modes approaches zero or the Alfvén speed. In a symmetric wide incompressible slab, the phase speeds converge to the same speed (Figure 4a), whereas in an asymmetric slab, the phase speeds converge to different speeds (Figures 4b -4d) that depend upon the values of the external densities. This observation is mirrored by both fast and slow surface modes in the more general solutions of a compressible slab solved numerically to give Figure 7a.

Low-Beta Approximation
The following section concerns the case when the magnetic pressure strongly dominates the gas pressure within the slab, i.e. β := 2μ 0 p 0 /B 2 0 1. This is known as the low-beta approximation and corresponds to the Alfvén speed dominating the sound speed in the slab; this provides a good approximation of the solar coronal environment. However, for application here an external magnetic field needs to be considered; this is the subject of a follow-up article.
Under this speed ordering, m 2 0 ≈ k 2 − ω 2 /v 2 A . After a numerical investigation, it is clear that the frequency of waves in this approximation satisfies ω 2 k 2 v 2 A , in which case m 2 0 ≈ k 2 provides a valid approximation. This means that m 2 0 > 0 and the solutions are surface modes. For a symmetric slab of low-beta plasma (e.g. Roberts 1981b), the dispersion relation reduces to a quadratic expression in ω 2 whose solutions are the fast sausage and kink surface modes given by where c e is the external sound speed, along with a spurious solution.
Unfortunately, for the more general case of an asymmetric slab of low-beta plasma, the dispersion relation does not reduce to an analytically solvable equation. However, we find numerically that there are two fast surface modes. The quasi-sausage surface mode is not present for small kx 0 , but becomes a solution at an intermediate value of kx 0 with phase speed ω 2 /k 2 = min (c 2 1 , c 2 2 ). The quasi-kink surface mode is present for all values of kx 0 . Qualitatively, the solutions for a low-beta plasma are analogous to the fast quasi-sausage and quasi-kink mode solutions, discussed later in Sections 3.4 and 3.5.
In the following section, thin and wide slab approximations are made to the approximate dispersion relation, Equation (22), rather than the exact dispersion relation, Equation (20). Here, the aim is to retrieve the variety of wave modes with future applications in mind.

Thin Slab Approximation
Consider the case where the wavelength, λ, of the waves propagating in the system is much greater than the width of the slab, 2x 0 , i.e x 0 /λ = kx 0 /2π 1. First, consider the quasi-sausage surface modes, which are governed by the tanh version of Equation (22), for m 2 0 > 0. In the thin slab limit, this equation reduces to Clearly, ω 2 = k 2 v 2 A is a solution, but as previously noted, this is spurious. The other solution for ω 2 behaves like ω 2 → k 2 c 2 T as kx 0 → 0. To first order in kx 0 , this solution is a slow quasi-sausage surface mode given by which is less than k 2 c 2 T and exists only when c 1 > c T and c 2 > c T . It is interesting to note that if c 1 = c 2 = c e (and therefore ρ 1 = ρ 2 = ρ e by Equation (3)), then there exists a second solution to Equation (30). By letting ω 2 = k 2 c 2 e (1 + ν) for some ν 1 in Equation (30), we find the solution in the thin slab limit. This is a fast sausage surface mode, and it degenerates (as a solution in the thin slab limit) as c 1 and c 2 become distinct, but this mode can still exist with a phase speed below the cut-off at min(c 1 , c 2 ). By solving the dispersion relation numerically, it can be shown that although the fast quasi-sausage surface mode does not exist in the thin slab limit for an asymmetric slab, it can exist in the wide slab limit (i.e. as a solution to Equation (36)) when c 0 < v A , c 1 < v A , and c 2 < v A . For example, in Figure 5a, the minimum of c 1 and c 2 becomes a new cut-off, causing the fast quasi-sausage surface mode to degenerate for small kx 0 . Next, consider quasi-kink surface mode solutions in the thin slab limit, which are governed by the coth version of Equation (22), for m 2 0 > 0. As kx 0 → 0, we have m 0 x 0 → 0, so coth m 0 x 0 → 1/m 0 x 0 . This simplifies the dispersion relation to This is a slow quasi-kink surface mode that behaves like ω/k → 0 in the thin slab limit. For body waves in the thin slab approximation, following the same procedure as for surface waves turns out to be fruitless, so we must reconsider our assumptions. Unfortunately, letting m 0 x 0 → 0 as kx 0 → 0, whilst valid for surface modes, does not exhaust all of the modes of the slab. Instead, we must consider the scenario where m 0 x 0 remains finite as kx 0 → 0. This can occur only if |m 2 0 | → ∞ as kx 0 → 0. To ensure |m 2 0 | → ∞, we are restricted to solutions that behave like ω 2 → k 2 c 2 T as kx 0 → 0. Upon consideration of Equation (22), this can only be the case when m 2 0 < 0, i.e. only for body modes. To find these solutions, set ω 2 = k 2 c 2 T (1 + ν(kx 0 ) 2 ) for some ν > 0 that is to be determined. To see why this form has been chosen, a substitution into the definition of m 2 0 demonstrates that |m 2 0 | → ∞ and m 0 x 0 remains bounded as kx 0 → 0, as required. Using this ansatz, Equation (22) has a countably infinite set of quasi-sausage body solutions, which in the thin slab limit, behave like There are also quasi-kink body solutions that in the thin slab limit behave like Equations (34) and (35) show us that to quadratic order in kx 0 , the quasi-sausage and quasikink body modes do not depend on the external environment parameters. The effects of external density and temperature are felt in the higher order terms, which explains why Equations (34) and (34) are identical to the corresponding solutions in a thin symmetric slab derived by Roberts (1981b).

Wide Slab Approximation
We now turn our attention to the behaviour of solutions to the dispersion relation in the wide slab limit, kx 0 1. As the width of the slab increases with respect to the wavelength, the slab boundaries have diminishing effect on the modes and each other.
The supposition that m 0 x 0 1 when kx 0 1 (which is verified by Roberts, 1981b) leads to the approximation tanh m 0 x 0 ≈ 1, which reduces the dispersion relation, Equation (22), to A comparison with Roberts (1981a) demonstrates that in the wide slab limit, the dispersion relation reduces to the single interface problem, but with a modified density ratio and wavenumber. This does not qualitatively change the solutions. Unfortunately, the body waves have no parallel in the single interface model because body waves in a slab that are externally evanescent owe their existence to the two interfaces. For body modes in the wide slab limit, there exists a solution that behaves like ω 2 → k 2 c 2 0 as kx 0 → ∞. To see this, substitute the ansatz ω 2 = k 2 c 2 0 (1 + ν/(kx 0 ) 2 ) into the dispersion relation to retrieve the family of quasi-sausage body modes given by in the wide slab limit. Similarly, there exist quasi-kink body mode solutions given by in the wide slab limit. These are valid solutions for the characteristic speed ordering v A > c 0 . This analysis may then be conducted for v A < c 0 to find that in the wide slab limit, there exist quasi-sausage body mode solutions given by and quasi-kink body mode solutions of the form These solutions to the approximate asymmetric slab dispersion relation demonstrate that to inverse quadratic order of kx 0 , the wide slab body modes are independent of the external plasma parameters. Therefore, Equations (37), (38), (39), and (40) are identical to the body mode solutions in a wide symmetric slab (Roberts, 1981b). Equations (31) -(35), (37) and (38) also appear in Li, Habbal, and Chen (2013) for a symmetric magnetic slab with shear flow when the shear flow speed is set to zero.

Numerical Solutions
In this section, the asymmetric dispersion relation is solved numerically, with particular interest placed on the effect of changing the ratio of the external densities.

Density Ratio Variation
First, consider a magnetised slab with symmetric non-magnetic external plasma, as described by Roberts (1981b). Figures 6a -6c illustrate how varying the ratio of external to internal density affects the propagation speeds of the slow kink and sausage surface modes. An increase in the density ratio, ρ e /ρ 0 , causes a decrease in the propagation speed of the slow modes. The fast surface modes demonstrate an identical behaviour (not shown). The body modes are weakly dependent on the external density, so that the propagation speed decreases only negligibly as ρ e /ρ 0 increases. More generally, consider an asymmetric slab whose equilibrium conditions are given by Figure 1. In Section 2.1 it was shown that the dispersion relation does not decouple into separate sausage and kink mode equations. However, their characteristics remain by and large, therefore the labels of "quasi-sausage" and "quasi-kink" may be used by referring to the anti-phase and in-phase oscillatory behaviour of the slab boundaries. However, we note that the quasi-kink mode now does not retain the width of the perturbed slab like the symmetric kink mode does. Figures 7a -7e illustrate the behaviour of the slow surface modes as the external density on one side of the slab is varied while holding fixed the other external density. The slice where ρ 1 /ρ 0 = 2 corresponds to a symmetric slab, where the usual behaviour is observed: the phase speed of the two slow surface modes converges to a speed that is slower than the tube speed, c T , as the slab width increases. However, as the external densities become distinct, the phase speeds of these modes become distinct for a wide slab. This can also be seen in Figures 5a and 5b.
For a wide slab width, kx 0 1, Figure 7e illustrates that the eigencurves of the slow surface modes demonstrate a wave phenomenon known as "avoided crossing". Avoided crossings occur when the phase speeds of two wave modes avoid crossing when a parameter of the system is varied due to constraints preventing them from being equal; it demonstrates a transferral of properties between the two modes and can be used to give insight into the modal structure. There is rich literature regarding avoided crossings for the eigensolutions of a wide range of physical processes including energy level repulsion in quantum physics (Naqvi and Brown, 1972) and coupled spring oscillations in classical mechanics (Novotny, 2010); in MHD the subject has been covered only briefly, for example between fast and slow magneto-acoustic gravity waves in a magnetically stratified plasma by Abdelatif (1990) and Mather and Erdélyi (2016).
In the present study, the avoided crossing occurs between quasi-kink and quasi-sausage surface solutions to the asymmetric slab. This explains why the dispersion relation does not decouple into a sausage and a kink equation (Section 2.1). Figures 8a and 8b demonstrate that during the transition across the avoided crossing, the quasi-sausage and quasi-kink modes exchange the slab boundary upon which the largest perturbation occurs. For example, the left plots of Figure 8b show that the quasi-sausage mode has its highest amplitude on the interface of highest local phase-speed (equivalently, lowest external density). The quasi-kink mode demonstrates the opposite behaviour. The central plots show the special case of a symmetric slab, where ρ 1 = ρ 2 , demonstrating the spatial antisymmetry and symmetry in the symmetric sausage and kink mode, respectively. As the left external density, ρ 1 , dominates the right external density, ρ 2 , the right plots of Figure 8b show that again the quasi-sausage mode has its higher amplitude on the interface of higher local phase-speed, but this is now on the other interface. By the term local phase speed we are referring to the phase-speed of a slow surface mode travelling parallel to the magnetic field at an isolated density and magnetic interface between a magnetised and non-magnetised plasma, i.e. half of the magnetic slab model. Each interface of an asymmetric magnetic slab has a distinct "local" phase-speed, and according to Roberts (1981a), these phase speeds are inversely proportional to the density in the non-magnetic region.
When they exist, the fast quasi-sausage and quasi-kink surface modes demonstrate an identical behaviour (not shown). As demonstrated analytically in Equations (34) -(35) and (37) -(40), the body modes are not dependent on internal or external densities to leading order in kx 0 . This means that body modes demonstrate only a weak dependence on the external densities, and an avoided crossing does not occur between these modes as the external densities are varied.   (20), are plotted showing the variation of the dispersion as the ratio of one external density to the internal density is changed. The other density ratio is held fixed at ρ 2 /ρ 0 = 2 and the non-dimensionalised slab width kx 0 = 1.5. The characteristic speed orderings are c 2 = 0.7c 0 , v A = 0.4c 0 , and c 1 varies to satisfy equilibrium pressure balance, given by Equation (3). The parameters at each blue and red dot in panel (a) are used to plot the spatial variation of the transverse displacement perturbation, ξ x , given by panel (b). The upper (lower) plots in panel (b) correspond to the quasi-sausage (quasi-kink) mode solutions.

Discussion
For the first time, a mathematical model of an isolated magnetic slab in an asymmetric environment has been presented, generalising the classic symmetric magnetic slab model (Roberts, 1981b). Key analytical results demonstrate the fundamental differences between the behaviour of linear waves propagating along asymmetric and symmetric slabs.
Sausage and kink modes are traditionally thought of as the fundamental modes of symmetric magnetic slab and tube geometries. Unlike that of a symmetric slab, the dispersion relation governing linear waves along a magnetic slab in a non-magnetic asymmetric external environment does not decouple into two equations, which signifies that the eigenmodes of an asymmetric slab are not the pure sausage and kink modes that we are familiar with; instead, they are adjusted by the asymmetry in the external region and demonstrate mixed properties. For example, the quasi-kink mode does not have a spatially constant perturbed slab width like the symmetric kink mode.
To understand this further, we propose an analogy with a coupled mechanical simple harmonic oscillation system (more information can be found in Novotny, 2010). Consider a system of two identical masses, between two fixed walls, with springs connecting the left wall to the left mass, the masses together, and the right mass to the right wall (Figure 9a). The three springs correspond to the three regions of plasma, and the masses to the interfaces. If the left and right springs have the same spring constant, then the two normal modes of oscillation are the oscillations where the masses are in phase (Figure 9b) or in anti-phase ( Figure 9c). These are analogous to the kink and sausage surface modes, respectively, in a symmetric magnetic slab. Consider now the case when the spring constants in the left and right springs are distinct (Figure 9d). In this case, the normal modes are asymmetric in-phase ( Figure 9e) and in-anti-phase oscillations (Figure 9f), analogous to the asymmetric quasikink and quasi-sausage magnetic slab modes, respectively. In this mechanical analogy, it is clear that for the in-phase mode, the higher amplitude will be on the side with smaller spring constant, and vice versa for the in-anti-phase mode. A higher spring constant in the left or right springs in this analogy corresponds to a lower density outside the magnetic slab because a higher spring constant in an uncoupled spring gives a higher characteristic frequency. This gives us some motivation as to why the surface modes of the asymmetric magnetic slab have higher amplitudes on different sides for quasi-sausage and quasi-kink modes.
Mixed properties and mode coupling of asymmetric MHD waves are very interesting for solar physics because they signify interaction between different modes, and interaction can mean energy transfer. This indicates a potential plasma heating mechanism: waves with negligible dissipation are excited in the solar interior and propagate through the photosphere and chromosphere until they interact and transfer energy to modes with rapid dissipation in the upper solar atmosphere, i.e. corona, where the kinetic and magnetic energy of the wave is converted into thermal energy (Priest, 2014).
There is potential use of the ideas presented here as a diagnostic tool in the emerging field of solar magneto-seismology. Without loss of generality, consider a magnetic slab of plasma in an asymmetric environment with a higher density in the left external plasma than on the right (e.g. Figures 2a and 2b). The ratio between the amplitudes of oscillation on the left and right boundaries will be a function of the equilibrium and wave parameters of the slab. This ratio and several other parameters such as the slab width, the wave frequency, and the wavelength are directly observable, and parameters such as the density and temperature in each region can be determined through emission spectra, leaving only the magnetic field parameters unknown. Thus, in the context of solar observations, the cross-slab amplitude ratio could be used to calculate difficult-to-measure solar parameters such as the Alfvén speed and magnetic field strength when observing MHD waves in slab-like structures.
A similar solar magneto-seismological tool that shows its potential is the novel anti-node shift method, which uses the shift due to density structuring of the anti-nodes of standing modes to deduce the parameters of solar environments Verth et al., 2007;Erdélyi, Hague, and Nelson, 2014). The success of this method indicates the need in the solar physics community for novel solar magneto-seismology methods.
There are clear future generalisations to the asymmetric slab model presented here. The consideration of a constant, uniform shear flow in one or both external environments would generalise the symmetric slab with shear flow studied by Li, Habbal, and Chen (2013), and introduce the additional physics of the Kelvin Helmholtz instability and negative energy waves. This would allow for better application of the asymmetric slab model to dynamic solar environments such as the solar wind, where magnetic structures with flows have been modelled using tube and symmetric slab geometries (Parker, 1965;Nakariakov, Roberts, and Mann, 1996;Taroyan and Erdélyi, 2002;Ruderman, 2010;Terradas et al., 2011;Yu et al., 2016). The generalisation of including a uniform external magnetic field and gravity, either in the direction parallel or perpendicular to the slab, would improve the applicability of the model, at the cost of analytic difficulty. With the completion of the next generation of solar telescopes, including the eagerly awaited Daniel K. Inouye Solar Telescope (DKIST), observations of fine-scale solar magnetic structures look forward to significant improvements in both spatial and temporal resolution; the future looks bright for solar magneto-seismology.