## MHDlensing in inhomogeneous ISM for qualitative understanding of the morphology of supernova remnants

## Yoshiaki Sofue

Institute of Astronomy, The University of Tokyo, Mitaka, Tokyo 186-0015, Japan

Accepted for publication in Astrophysics &amp; Space Sciences (ApSS)

## ABSTRACT

Morphological evolution of expanding shells of fast-mode magnetohydrodynamic (MHD) waves through an inhomogeneous ISM is investigated in order to qualitatively understand the complicated morphology of shell-type supernova remnants (SNR). Interstellar clouds with high Alfvén velocity act as concave lenses to diverge the MHD waves, while those with slow Alfvén velocity act as convex lenses to converge the waves to the focal points. By combination of various types of clouds and fluctuations with different Alfvén velocities, sizes, or wavelengths, the MHD-waveshells attain various morphological structures, exhibiting filaments, arcs, loops, holes, and focal strings, mimicking old and deformed SNRs.

Key words: ISM: magnetic fields - ISM: clouds - ISM: supernova remnant - magnetohydrodynamics (MHD)

## 1 INTRODUCTION

The interaction of expanding blast waves of supernova remnants (SNR) with the interstellar medium is one of the most fundamental subjects of interstellar physics (McKee &amp; Cowie 1975; TenorioTagle et al. 1985; Stone &amp; Norman 1992; Shelton et al. 1999; Chevalier 1999). The shock compression of gas clouds and magnetic fields by SNRs is also an important issue in the high-energy astrophysics for investigating the acceleration mechanisms of cosmic-rays. Extensive observations of the interacting clouds near SNRs have been obtained from 𝛾 -ray to radio wavelengths, on which several models for cosmic-ray accelerations have been proposed (Fujita et al. 2009; Acciari et al. 2009; Inoue et al. 2012; Cardillo et al. 2014; Archambault et al. 2017; Sano et al. 2019, 2020; Abeysekara et al. 2020; Giuffrida et al. 2022; Sakemi et al. 2023).

the SNR shell. The interaction of an SNR with the ISM takes place most directly at the contact surface of the outermost blast wave, which is observed as a thin shell with optical emission due to the recombination process between the blast and cold ISM as well as the non-thermal radio emission due to the Fermi acceleration of cosmic-rays at the highest compression surface.

Multi-wavelengths observations from 𝛾 to radio have been obtained also for investigating the morphological evolution of SNRs interacting with the ISM, which are used at the same time to diagnose the physical conditions and magnetic fields of the ISM (Bocchino et al. 2000, 2011; Lopez et al. 2011; Miceli et al. 2006, 2017; Bykov et al. 2008; West et al. 2016; Sano et al. 2021; Orlando et al. 2022; Ohmura et al. 2021). In order to explain the observed properties and morphologies of SNRs, a number of numerical simulations of hydro- and magneto-hydrdoyanamics have been obtained about the interaction of shock waves with the inhomogeneous/turbulent ISM and stellar winds (Velázquez et al. 2002, 2003; Schneiter et al. 2006; Fang et al. 2014; Toledo-Roy et al. 2014; Moranchel-Basurto et al. 2017; Zhang &amp; Chevalier 2019; Ustamujic et al. 2021; Kupilas et al. 2022).

The current studies of the SNR-cloud interaction have been devoted mainly to study the effects of SNRs as the active source on the shocked medium to change the physical conditions interior to

At the same time, the morphology of SNRs also manifests the past paths of propagation of the shell front through the interstellar medium, because the blast wave propagation is strongly influenced by the density and magnetic fields in the ISM. This means that the morphology of SNRs may be used to diagnose the gas density and magnetic field distributions in the ISM before the passage of SNRshells (Sofue 1978). Therefore, morphological study of images of aged SNRs in optical (van den Bergh et al. 1973) and radiocontinuum emissions (Green et al. 2014; Kothes et al. 2006; Helfand et al. 2006), would provide with a tool to diagnose the interstellar medium by analysing the shapes that have been deformed during the propagation.

In this paper we aim at understanding the morphology of old shell-type SNRs by the deformation of expanding shells in the inhomogeneous ISM. We show that various types of morphology of old SNRs are qualitatively understood by the lensing effect due to the variation of propagation speed of the waves in the ISM. Since the treatment of the wave in this paper is linear, we discuss only the geometrical shapes of the front, or the morphology, of old SNRs. Non-linear physical quantities such as the density, magnetic strength, cosmic-ray acceleration, emission rates, etc.. cannot be calculated in the present linear approximation. However, the relative distribution of the brightness of wave front, which defines the morphology, can be mapped by the distribution of column density of wave packets projected on the sky. Thereby, the brightness of a sheet of MHD wave is highest, when it is seen edge on and the

tangential sheet shines brightest, whereas a face-on sheet is almost invisible.

In section 2 we formulate the method using fast-mode MHD waves. In section 3 we describe result of simulation for large-scale deformation of shell fronts. Section 3 deals with more complicated morphology, and section 4 compares the model with the well-known SNRs having unique morphology. In section 6 we present some speculative considerations towards understanding the morphology of SNRs in more general cases based on the MHD model, and discuss the difficulty and limitation of the model. Section 7 summarizes the results.

## 2 MHDLENSING

## 2.1 Fast-mode MHD wave and Alfvén velocity in the ISM

The Alfvén velocity in a medium with magnetic-field strength 𝐵 and gas density 𝜌 is given by

$$V & = B / \sqrt { 4 \pi \rho }. & ( 1 ) \quad \ \exp \text{lo}$$

Its variation with the cloud density depends on the coupling property of the magnetic field and gas. There are three types of coupling:

- · (a) Frozen-in field depends on the gas density as 𝐵 ∝ 𝜌 2 3 / , which leads to Alfvén velocity varying as 𝑉 ∝ 𝜌 1 6 / . This happens when the magnetized gas cloud is highly turbulent.
- · (b) Cylindrical magnetic field uniformly filled with the gas varies as 𝐵 ∝ 𝜌 , so that we have 𝑉 ∝ 𝜌 1 2 / . This occurs in stretched clouds and filaments extending along the magnetic lines of force.
- · (c) Gas flowing along magnetic lines of force with the field strength kept consant yields variation of the Alfvén velocity as 𝑉 ∝ 𝜌 - / 1 2 . Such a flow is found in a magnetic tube inflating from the Galactic plane by the Parker instability with a condensed cloud at the bottom of the tube.

Interstellar clouds in cases (a) and (b) yield higher 𝑉 than the ambient medium, so that they act as concave lenses to diverge the MHD waves. Clouds in case (c) act as convex lenses to converge the waves. In this paper we represent the property of a cloud by its 𝑉 variation relative to the mean value, and represent the variation by equation 8.

MHD disturbances excited by explosive events such as due to SNe in the interstellar space propagate as spherical shock waves in the initial expansion phase. In the fully expanded phase when the expanding velocity becomes comparable to the local Alfvén velocity or sound velocity, they propagate as the fast-mode MHD compression wave (hereafter, MHD wave), Alfvén wave, and the sound wave.

The typical Alfvén velocity is 𝑉 ∼ 10 km s -1 with magnetic strength 𝐵 ∼ 5 𝜇 Gandgasdensity 𝜌 ∼ 1Hcm -3 , whereas the sound velocity is 𝑐 s ∼ 0 1 and 1.3 km s . -1 for molecular and HI gases with temperatures 𝑇 mol ∼ 20 K and 𝑇 HI ∼ 100 K, respectively. Therefore, among the three modes, the fast MHD wave is most efficient to convey the energy into the ISM, which propagates at the group velocity approximately equal to the Alfvén velocity, 𝑣 g ≃ 𝑉 , where the Alfvén velocity is higher than the sound velocity, or the magnetic pressure is higher than the thermal pressure (Porter et al. 1994).

## 2.2 Alfvén vs turbulent velocities

Oneconcernthatdisturbs the above assumption is that the Alfvén velocity is comparable to the turbulent velocity as the consequence of the energy balance (equilibrium) between the magnetic and thermal energy densities in the Galactic disc. Because the Eikonal equations hold locally, the wave propagation is superposed by the local turbulent motion. This effect is significant on the ray paths slowly bent (refracted) by low𝑉 clouds. On the other hand, it does not affect so much the ray paths against strongly disturbed (reflected) ray paths by high𝑉 lensing. In this paper, we assume that the turbulent motion is negligible, which, even if existed, disturbs the result on the same order as that due to low𝑉 lensing. The amount of ray deflection is much smaller than that by high𝑉 lensing, and the entire morphological shape is more effectively determined by the high𝑉 lensing.

## 2.3 Propagation of MHD waves by Eikonal equations

The propagation of MHD waves is traced by solving the Eikonal equations developed for the Solar coronal Moreton waves (Uchida 1970, 1974; Zhou et al. 2024). The method has been applied to the SNRs (Sofue 1978, 2023) as well as to the Galactic Centre explosions (Sofue 1978, 2023). Gas clouds with higher density with constant and weaker magnetic strength act as a convex lens (low𝑉 lens) gathering the waves. On the other hand, gas clouds with frozen-in magnetic field and gaseous cavity (hole) with constant and higher magnetic strength act as a concave lens (high𝑉 lens) reflecting the waves.

The Eikonal equations describing the propagation of an MHD wave packet are given as follows (Uchida 1970, 1974).

$$\underset { \dots } { \text{ched} } \quad \frac { d r } { d t } = V \frac { p _ { r } } { p },$$

$$\begin{smallmatrix} \cdot _ { s } \infty \\ \text{from} \\ \cdot \text{load} \end{smallmatrix} \quad \frac { d \theta } { d t } = V \frac { p _ { \theta } } { r p },$$

$$\text{bient} _ { \text{HD} } \quad \frac { d \phi } { d t } = V \frac { p _ { \phi } } { r p \sin \theta },$$

$$\Psi ^ { y \text{ its} } _ { \text{ n bv} } \quad \frac { d p _ { r } } { d t } = - p \frac { \partial V } { \partial r } + \frac { V } { r p } ( p _ { \theta } ^ { 2 } + p _ { \phi } ^ { 2 } ),$$

$$\underset { \text{waves} } { \text{ue to} } \quad \frac { d p _ { \theta } } { d t } = - \frac { p } { r } \frac { \partial V } { \partial \theta } - \frac { V } { r p } ( p _ { \theta } p _ { r } - p _ { \phi } ^ { 2 } \cot \theta ),$$

$$\underset { d \text{ the} } { \text{view} } \underset { d \text{ the} } \frac { d p _ { \phi } } { \text{d} } = - \frac { p } { \sin \theta } \frac { \partial V } { \partial \phi } - \frac { V } { r p } ( p _ { \phi } p _ { r } + p _ { \phi } p _ { \theta } \cot \theta ),$$

where 𝑉 = 𝑉 𝑟, 𝜃, 𝜙 ( ) = 𝑉 𝑥, 𝑦, 𝑧 ( ) is the Alfvén velocity, the vector p = ( 𝑝 𝑟 , 𝑝 𝜃 , 𝑝 𝜙 ) = grad Φ is defined by the gradient of the eikonal Φ , 𝑝 = ( 𝑝 2 𝑟 + 𝑝 2 𝜙 + 𝑝 2 𝜃 ) 1 2 / , ( 𝑟, 𝜃, 𝜙 ) and ( 𝑥, 𝑦, 𝑧 ) are the polar and Cartesian coordinates.

An interstellar cloud with lower Alfvén velocity deflects the path toward its center, acting as a convex lens, and the converging waves are focused around its focal point to produce amplified wave front. On the other hand, a cloud with higher Alfvén velocity acts as a concave lens, and diverges the waves to produce a hole with the deflected waves making a ring.

Figure 1 shows the propagation of MHD wave packets produced at the center intervened by clouds with higher and lower Alfvén velocities than the ambient velocity, where the 𝑉 variation is expressed by Gaussian function with an 𝑒 -folding scale radius of 0.5.

In the following simulations, we solve the Eikonal equations

## (A) Low𝑉 cloud (convex lens)

(B) High𝑉 cloud (concave lens)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-77bcb844-09e5-4fc3-902b-542f7cfc64cc.jpg)

Figure 1. MHDray path from an explosion center encountered by spherical clouds located at 𝑥 = 1 with (A) lower and (B) higher 𝑉 , acting as convex and concave lens, respectively (Sofue 1978).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-b0ef88d3-c720-4fe7-aa7e-f308409cf4de.jpg)

assuming the functional form of the 𝑉 variation in the space, which is expressed in the following form:

$$V & = V _ { 0 } [ 1 + \epsilon \, f ( x, y, z ) ], & ( 8 ) & \quad V = 1$$

where 𝑉 0 is the mean 𝑉 in the ISM, 𝑓 ( 𝑥, 𝑦, 𝑧 ) is a function expressing the shape of the turbulence, and 𝜖 represents the relative amplitude of variation. Here, 𝜖 ≥ -1, and low- and high𝑉 clouds have negative and positive 𝜖 , respectively.

## 3 LARGE-SCALE DEFORMATION

## 3.1 Lopsidedness by 𝑉 gradient

The gradient of Alfvén velocity produces a lopsided elongation of the front shape, so that the lower𝑉 side is stopped, and the shell expands faster in the higher𝑉 side. Figure 2 (A) shows a case when the wave expansion takes place in an ISM with one-directional

𝑉 variation expressed by

$$V = V _ { 0 } \, [ 1 + \exp ( z / h ) ],$$

where 𝑉 0 and ℎ = 0 5 are constants. . Hereafter, the coordinates are ( 𝑥, 𝑦, 𝑧 ) measured in arbitrary unit. The time is measured in unit of unity scale divided by 𝑉 , 𝑡 unit = 1 / 𝑉 0 . In this diagram, the wave front at 𝑡 = 0 25 0 5 0 75 and 1 . , . , . 𝑡 unit are plotted. The SN explosion takes place at ( 𝑥, 𝑦, 𝑧 ) = ( 0 0 0 , , ) spherically, and the wave front is spherical in the initial stage. As it expands, the front is elongated in the 𝑧 direction of increasing 𝑉 , and attains an open cone structure. The front finally forms a one-sided curved sheet open toward positive higher𝑉 direction.

If the wave expansion takes place in a disc with increasing 𝑉 gradient toward both sides expressed by a cosh function,

$$V = V _ { 0 } [ 1 + 0. 5 \cosh ( z / 0. 5 ) ],$$

the front shape becomes bi-conical, as shown in figure 2 (B) for 𝑧 0 = 0. If the disc is thin enough compared with the expanding shell radius, the front becomes bipolar, exhibiting a dumbbellshaped bubbles for negative-Gaussian disc (figure 3 (C)). If the disc has higher𝑉 than the surrounding ISM having positive-Gaussian distribution, the front attains bipolar,or bilateral caps (figure 2 (D)), somehow mimicking such an axisymmetric bilateral SNR as G296.5+10.0 (Giacani et al. 2000).

## 3.2 High- and low𝑉 walls

If the MHD-wave front encounters a wall with sufficiently high 𝑉 with finite thickness, the front is reflected and the shell is deformed to become a half-spherical shape. Figure 3 (A) shows an example of such an encounter, where a horizontal wall is located at 𝑧 0 = 0 5 . expressed by

$$V = V _ { 0 } [ 1 + \exp ( ( z - z _ { 0 } ) / h ) ^ { 2 } ]$$

with 𝑧 = 0 5. .

If the wall has lower 𝑉 , the front is strongly refracted toward the wall and absorbed there, resulting in a half shell, and dense, flat front along the wall is formed, as shown in panel (B) of figure 3.

## 3.3 Clouds with low𝑉 vs high𝑉 : Convex vs concave lenses

A low𝑉 cloud acts as a convex lens for the ray paths of the MHD wave packets. Figure 4 shows front propagation of a spherical MHD wave front encountering a low𝑉 cloud, where the Alfvén velocity is expressed by a negative Gaussian function as

$$( 8 ) \quad V = V _ { 0 } [ 1 + \epsilon \exp ( - ( r / r _ { 0 } ) ^ { 2 } ) ].$$

Here, 𝑟 = [( 𝑥 -𝑥 0 ) 2 + ( 𝑦 -𝑦 0 ) 2 + ( 𝑧 -𝑧 0 ) 2 ] 1 2 / is the distance from the cloud center at ( 𝑥, 𝑦, 𝑧 ) = ( 𝑥 0 , 𝑦 0 , 𝑧 0 ) ∼ ( 0 5 . , 0 5 0 5 , , , , ) 𝑟 0 ∼ 1 is the radius of the cloud, and 𝜖 ∼ -0 5 is the negative amplitude of . 𝑉 .

The front is reflected by the lens and bent toward the focal point of the cloud. Accordingly, the wave forms a concave surface shape and the amplitude increases.

If 𝑉 of the cloud is small enough, the front evolves into an elephant-trunk shaped elongated lobe. However, the focusing effect to the trunk's head amplifies the magnetic field together with wave density, which works to increase the Alfvén velocity and suppresses the growth of the trunk.

On the other hand, when the Alfvén velocity of the cloud is higher than the ambient medium, it acts as a concave lens. The ray

Figure 2. Explosion in a medium having (A) exponential-gradient in 𝑧 direction, in a (B) hyperbolic-cosine 𝑧 disc, (C) flat negative Gaussian disc (low𝑉 disc), and (D) flat positive Gaussian disc (high𝑉 disc).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-7a765c50-3da5-4c8c-888c-0a8dc37cc57c.jpg)

Figure 3. Encounter with widely extended (A) high- and (B) low𝑉 plate (wall).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-4337d7de-6968-4246-9edf-0fbe86692b78.jpg)

paths are diverged by the lens, and the front expands locally, making a vacant hole. Thereby, reflected rays are locally converged on the edge of the hole, and form a focal ring (loop) of amplified waves. Figure 4 shows an example of evolution of the front through such a high𝑉 cloud with 𝜖 ∼ + 0 5. .

## 3.4 Large-scale deformation by large clouds

## 3.4.1 IC443 and Sh 2

Large-scale deformation such as observed for the SNRs IC 443 (figure 5) can be explained by interaction of a larger-sized clouds and density fluctuations than the size of the expanding front. When a

(B) High𝑉 cloud

Figure 4. (A) Encounter with a low 𝑉 cloud. ) (B) Encounter with a high-V cloud, making a hole surrounded by a ring of amplified wave density.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-ac56db7a-74b8-4edd-97d2-9f507553dace.jpg)

spherical front encounters a wall or a sheet with higher𝑉 , the front suddenly expands faster and the shape becomes lopsided, so that one side of the shell is more extended. Our MHD-wave simulation of IC 443 is shown in figure 5, where the large-scale deformation of the shell is reasonably mimicked in coincidence with that obtained by the 3D numerical simulations (Ustamujic et al. 2021).

( 𝑥, 𝑧 ) plane ( 𝑦 ∼ 0) across the explosion center, yielding wavy front around the equator. When the shell is seen from the pole direction ( 𝑦 axis) axis, the projected shell exhibits wavy loop and arcs due to the tangential view of the corrugated surface of figure 7) (A). When the shell is viewed from the 𝑥 direction, filament structure is observed as equatorial filaments across the center.

## 3.5 Manatee (W50) SNR

When the SN explosion takes place in a high𝑉 tube, the front expands in the direction of the axis of the tube, attaining a spindle shape, as shown in figure 6. Here, the magnetic tube is further superposed with oblique sheets of variable 𝑉 , which produce the oblique stripes and wavy deformation of the front. The simulated shape mimics the Manatee nebula (W50), whose large-scale morphology (Farnes et al. 2017) can thus be possibly explained by the interaction of the SNR exploded in a magnetic tube, although more sophisticated model has been required to explain its extremely peculiar properties related to the central active source (Ohmura et al. 2021; Sakemi et al. 2023). The off-plane location and large angle of the major axis to the Galactic plane of the Manatee nebula (Dubner et al. 1998) suggests that the SNR exploded in a vertical magnetic tube inflating from the Galactic disc due to the Parker instability. The overlaid oblique stripes are also explained by corrugation of the front by sheer inhomogeneity of the medium, where the entire region is embedded in sheets obliquely tilted by 45 ◦ from the 𝑥 axis having sinusoidal variation of 𝑉 with 𝛿 V ∼ 0 01 and . 𝜆 ∼ 0 2. .

## 3.6 Corrugated filaments by sheet fluctuation

## 3.6.1 Filamentary loop

If ISM has sheet-like fluctuations in one direction, the MHD front is disturbed to form wavy fluctuations near the equator parallel to the sheets' plane. Figure 7 shows a case when the 𝑉 variation is sinusoidal in the 𝑦 direction by ∼ 1 percent of the ambient value,

$$V & = V _ { 0 } [ 1 + \Sigma _ { j } \epsilon _ { j } \cos ( 2 \pi y / \lambda _ { i } ) ], & ( 1 3 ) \quad \text{ is m} \, \\ & \text{cated}$$

where 𝜖 𝑖 ∼ 0 01 and . 𝜆 𝑖 ∼ 0 1. .

The perturbation influences the rays that propagate near the

000 , 000-000 (2024)

## 3.6.2 Thin filaments in Cygnus Loop

Filamentary structure in the tangential projection of the outer shell is more ubiquitous in older SNRs. The very thin filaments in the Cygnus Loop as shown in the right panel of figure 7 (B) (Blair et al. 2005) can be naturally reproduced by a tangential view of corrugated front, as simulated by the MHD-wave model shown in figure 7 (B).

## 4 TURBULENT AND PECULIAR MORPHOLOGY

## 4.1 Randomly distributed high- and low𝑉 clouds

When the MHD wave front expands into a turbulent medium, the front is deformed to produce clumps, filaments, arcs, loops and holes. We simulate three cases with randomly distributed 100 Gaussian clouds having random sizes, amplitudes, and positions in the plotted frame.

Figure 8 compares the time evolution of MHD fronts at 𝑡 = 1 and 2 unit , expanding into turbulent medium filled with randomly 𝑡 distributed high-V (left) and low-V (right) Gaussian clouds. High-V clouds deform the shell more strongly than low-V clouds.

If the medium is filled with high𝑉 clouds, the front is largely deformed to exhibit filaments, arcs, and holes as shown in figure 8 (A). When the interstellar medium is composed of low𝑉 random clouds, the front is deformed to make randomly focusing, concave, andlocally amplified clumps, as shown in figure 8 (B). If the medium is mixed with low- and high-V clouds, the front results in a complicated deformed shell with the mixture of concave clumps, holes, and filaments as shown in figure 8 (D). Here, the 𝑉 distribution in the ISM is given superposition of random Gaussian 𝑉 clouds expressed

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-f3079095-4ab7-4e6d-90fe-ff5e1774a4f8.jpg)

𝑉

(B) High- plate: Sh2

Figure 5. (A) Dual half shells with different sizes divided by the encounter with a high𝑉 plate, mimicking the SNR IC 443 ( https://apod.nasa.gov/apod/ap231226.html) (B) Dual half shells divided by a high𝑉 plate simulating SNR Sh 2-224 (wing-and-shell nebula, VRO 42.05.01) (Arias et al. 2024) (https://apod.astronomia.com/wp-content/uploads/2022/01/sh2-224.jpg).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-fa609b9d-9399-4dd7-92b2-3856ba324feb.jpg)

by

$$V = 1 + \Sigma _ { i } \delta V _ { i }$$

$$( 1 4 ) \quad v _ { i } = 2 ( \mathcal { R } ( 5, i ) - 0. 5 ),$$

with

$$\delta V _ { i } = v _ { i } \exp \left ( - \frac { ( x - x _ { i } ) ^ { 2 } + ( y - y _ { i } ) - 2 + ( z - z _ { i } ) ^ { 2 } } { w _ { i } ^ { 2 } } \right ), \quad \ \ ( 1 5 ) \quad \ ] \, \text{result}.$$

$$( 2 0 )$$

and R( 𝑗 , 𝑖 ) represents random number between 0 and 1, as shown in figure 8 (C).

where

𝑤 𝑖 = . + . × R( , 𝑖 ) ,

$$w _ { i } & = 0. 0 5 + 0. 1 \times \mathcal { R } ( 4, i ), & ( 1 6 ) \quad \text{ an e}.$$

$$x _ { i } & = 2 ( \mathcal { R } ( 1, i ) - 0. 5 ),$$

$$y _ { i } & = 2 ( \mathcal { R } ( 2, i ) - 0. 5 ), & ( 1 8 ) & \quad \text{strinc}$$

$$z _ { i } = 2 ( \mathcal { R } ( 3, i ) - 0. 5 ),$$

(19)

Expansion into a turbulent ISM filled with high-V clumps results in a front filled with cell-like clumps of arc-shaped filaments. As the wave front expands further, the turbulent ISM with high-V clumps causes more disturbed structure. The Vela nebula would be an example of such deformed front, where a number of filamentary arcs and loops are observed as shown in figure 9 (B).

## 4.2 Expansion into sheets and tubes: oblique stripes

Plane parallel sheets with sinusoidal 𝑉 variation creates longitudinal stripes. If it is superposed by sinusoidal variation in the 𝑦 and 𝑧 directions, mimicking high- and low𝑉 tubes, the longitudinal stripes appear to cover the surface as shown in figure 11, where the wavelengths are taken as 𝜆 𝑦 = 𝜆 𝑧 = 1.

## (A) Manatee: Explosion in a high𝑉 tube

Figure 6. Explosion inside high𝑉 cigar-shaped cloud, mimicking the SNR, Manatee nebula (W50) in radio continuum. The image was reproduced from Dubner et al. (1998), https:// www.nrao.edu/pr/ 2013/w50/ .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-5d9008e2-2b62-4cbb-9b90-f63fdf3a34dd.jpg)

## (A) Filaments as corrugated surface: SNR 0509-67.5

(B) Filaments near the edge of Cygnus Loop

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-16030913-5291-4ac4-839e-aaffb3bb594a.jpg)

Figure 7. (A) Sheet filaments by tangential view of a corrugated surface near the edge of an almost perfect shell. (B) Close up of the filaments in Cygnus Loop (Blair et al. 2005) (https://apod.nasa.gov/apod/ap200928.html ) compared with the simulation.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-09bad7e0-9bae-4f73-bbd8-606cc9eddb8d.jpg)

Figure 8. (A) MHD front in random (turbulent) Gaussian clouds: (B) Same, but in low𝑉 clouds. (C) Randomly distributed 100 low- and high𝑉 clouds. Blue clouds have low𝑉 , and red high𝑉 . Note, enlarged scale in this panel. (D) MHD front in the medium given by panel C.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-f648c192-3864-427d-9546-c7f92f4f1ac7.jpg)

(A) Filaments and loops by high- and low𝑉 clouds (B) The Vela nebula

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-2b6e6b62-4047-4d20-beea-f86be33a2552.jpg)

Figure 9. (A) MHD-wave front expanding in randomly distributed highand low𝑉 clouds. (B) The Vela nebula reproduced from https://apod.nasa.gov/apod/ap100910.html.

MHDwave front expanding into ordered ISM with directional filaments with high𝑉 results in a shell covered by filaments and stripes parallel to the perturbing ISM. Low𝑉 filaments result in less deformed filamentary surface. Figure 10 (A) shows the simulated results.

However, if the deformed filaments are viewed from the direction parallel to the filament axes, the deformed front resembles the result for random high𝑉 clouds. Figure 10 (B) shows the same result for filamentary ISM turbulence viewed from two different directions: The left panel is a view perpendicular to the filaments, and the right is parallel. The projection parallel to the filaments resembles to random-cloud cases.

## 4.3 Morphological evolution

Wehave so far presented snap shots of the simulated result. Besides snap shots, it is also interesting to see how the morphology evolves during the expansion. In figure 12 and 13 we show the time variation of an expanding MHD wave front in a turbulent medium filled with sinusoidal and randomly distributed high- and low𝑉 clouds.

The shell expands spherically in the initial stage, until the size becomes comparable to the scale size of the turbulence. After a while, the front is deformed to exhibit asymmetric and elongated shell due to the spatial gradients in the 𝑉 distribution. As it further expands, the deformation is amplified to attain a turbulent morphology covered by randomly distributed filaments, arcs, holes, and concave focal clumps.

## 5 QUALITATIVE UNDERSTANDING OF PECULIAR MORPHOLOGY OF SNRS

## 5.1 Cygnus Loop: Encounter with low𝑉 bars

The Cygnus Loop is known as the most typical shell-type SNR, but it also shows peculiarly open and complicated morphology divided by a bunch of long straight filaments crossing the center (Fesen et al. 2018). It consists of largely deformed arcs and filaments: The Veil Nebulae (NGC 6992) composes the eastern limb of a bunch of bright arcs and filaments, and the western limb also composes a bunch of veil nebulae and extends to the south. The entire SNR is largely deformed and is opened toward the south, making a long protrusion. The most spectacular structure is the straightly extending bunch of vertical (north-south) filaments, which separate the whole shell into eastern and western halves.

These striking features have been well simulated by a lopsided outburst of the progenitor star in a uniform ISM (Fang et al. 2017). We here try to simulate such peculiar morphology qualitatively by the lensing effect, where an initially spherical shell expands into an inhomogeneous medium. We assume that the front encounters a vertical (NS-oriented) long bar-shaped cloud with low𝑉 and a high𝑉 cloud near the southern edge. Figure 14 shows the results compared with a photograph of the SNR, where blue and red dots represent the low and high 𝑉 clouds, respectively. In the top panel we show a stereoscopic view of the simulated results, one panel rotated 10 ◦ from the other perpendicular to the line of sight to explore the 3D structure of the SNR.

## 5.2 The 𝜏 nebula: two orthogonal low𝑉 bars

The 𝜏 nebula (N49/DEM-L-90) in the LMC is an SNR with the most strange morphology, exhibiting two orthogonal bunches of straight filaments. In our simulation, such morphology is understood as a result of an encounter of the MHD front with two orthogonal barshaped clouds with low𝑉 . Figure 15 shows the simulation result and assumed low𝑉 clouds, as compared with the HST optical image.

## 5.3 SNR S147: loops, arcs, and converging filaments

The old shell-type SNR S147 has been simulated in our early paper (Sofue 1978) based on the optical image (van den Bergh et al. 1973). In figure 16 we reproduce our early work by the present simulation. We further performed higher-resolution simulation, and show the result in the lower panels of the figure in comparison with an optical image of S147 in the H 𝛼 line.

Low𝑉 clouds 1 and 2 causes the concave filaments due to the focusing effect. High𝑉 clouds 3 and 4 leads to the convex deformation in the eastern and western edges.

Figure 17 enlarges the NE filamentary region compared with the simulation. Due to the focusing effect by the low𝑉 cloud, wave front is strongly deformed and many converging waves produce a condensation of overlapped filaments, leading to a concave bright clump.

## 6 DISCUSSION

## 6.1 Limitation of the method and speculative application of the model to SNRs

By definition, the present MHD wave simulation cannot be applied to the early phase of SNR expansion, when the mass of plowed gas is less than the mass of the progenitor. Also, the method is not applicable to SNRs with expansion velocity much faster than the Alfvén velocity with 𝑉 expa &gt;&gt; 10 km s -1 as well. However, we argue that the result would be applied to shell-type SNRs in order to qualitatively understand their morphology

Strictly speaking, the Eikonal equation with linear wave treatment is not applicable to waves of finite amplitude, but such waves do not exist in nature, and this is true for any linear theory of astrophysics. Nevertheless, linear theories are widely used to understand essential aspects of phenomena. Here we discuss the possibility of applying the MHD wave method to more general SNR cases. We consider whether shock waves are subject to lensing effects, or reflection and refraction by an inhomogeneous ISM.

There have been extensive studies of the interaction of expanding shells of SNRs with interstellar clouds (McKee &amp; Cowie 1975; Stone &amp; Norman 1992; Inoue et al. 2012; Sofue et al. 2021) in the premise that the interaction may trigger star formation and acceleration of cosmic rays by the compression of the magnetized molecular gas. There have been also a number of hydrodynamical and MHD numerical simulations of the interaction with the clouds (Mac Low et al. 1994; Klein et al. 1994; Bedogni &amp; Woodward 1990; Nakamura et al. 2006; Velázquez et al. 2002, 2003; Schneiter et al. 2006; Fang et al. 2014; Toledo-Roy et al. 2014; Moranchel-Basurto et al. 2017).

Commonly seen in these simulations is the implosion of the shock front converging toward the cloud's center axis by refraction of propagation direction according to the deceleration of shock velocity due to momentum conservation. As a result, the shock wavesuffers from a similar lensing effect to the MHD-wave focusing around the low𝑉 cloud's center, as shown in figures 4 (A) and 17. On the other hand, no numerical simulation has been obtained yet

Figure 10. (A) Formation of filamentary MHD-wve fronts by encounter with elongated low- and high𝑉 strings of clouds. (B) Encounter with elongated clouds, but shows the difference of the appearances depending on the viewing angle. The left panel shows a perpendicular view, while the right is an end-on view parallel to the filaments.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-7c38b307-bb58-4f0a-adf7-c2601e945cc3.jpg)

(A) Sinusoidal sheets

(B) Sinusoidal filaments

Figure 11. MHDfronts interacting with sinusoidal 𝑉 sheets and tubes.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-84fe3d58-9f51-4655-85e7-e5a1f254afab.jpg)

Figure 12. Morphological evolution of MHD-wave front expanding into a medium with sinusoidal 𝑉 variation, plotted at 𝑡 = 0 5, 1, 1.5 and 2. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-c9b180ce-13b4-4dbc-808c-3debb19b4d98.jpg)

Figure 13. Morphological evolution of an expanding MHD-wave front into a turbulent medium filled with randomly distributed high- and low𝑉 clouds plotted at 𝑡 = 0 5, 1, 1.5 and 2. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-c5b4ca3d-f894-47ad-bb62-a44de9ecce2b.jpg)

Figure 14. (A) Stereographic presentation of the simulated MHD front for the Cygnus Loop. See the left and right panels individually with the left and right eyes, respectively. (B) Encounter with vertical bar-shaped clouds of low 𝑉 , resulting in a long, extended bunch of filaments, and with a high𝑉 cloud in the southern edge, mimicking the SNR, Cygnus loop. (C) Same, but a view from the top (north). (D) Optical image of the Cygnus Loop (North to the top: https://apod.nasa.gov/apod/ap191031.html ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-794fdfdb-50cb-46a0-b02b-86ae0d19070d.jpg)

(A) MHD wave front

(B) Two orthogonal bar clouds (C) The 𝜏 nebula.

Figure 15. Deformed MHD-wave front (A) by an encounter with (B) two orthogonal bar-shaped clouds having low 𝑉 , qualitatively mimicking (C) the 𝜏 nebula, N49 in LMC ( https://hubblesite.org/contents/media/images/ ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-f090fb96-235e-4dc9-976b-64ce4bb00fae.jpg)

about shock waves interacting with a gaseous cavity or a high𝑉 cloud containing compressed frozen-in magnetic fields. However, we may speculate that a shock wave encountering such a cavity would behave similarly to MHD-waves that expand inside the cavity and make a hole in the front surrounded by a focal ring, as shown in figure 4 (B). Also, a shock wave inside high𝑉 cloud would become sub-Alfvénic due to the compression of frozen-in fields, which will propagate as an MHD-wave, so that the behavior would become more similar to that of MHD waves.

In summary, shock waves in SNRs behave similarly to MHD waves and result in similar morphological evolution. The difference is that shock waves disturb/destroy the clouds they pass through, while MHD waves are passive and do not disturb the environment. The deformation of the shock front may be weaker than that of the MHD wave front, but the sense of deformation is similar to each other. We therefore believe that the morphological peculiarities observed in shell-type SNRs can be qualitatively explained by the lensing effect of the inhomogeneous ISM. Figure 18 shows some examples of such SNRs and the corresponding MHD wave models.

Through the simulation, we learnt that the shape of MHD wave front suffers from significant deformation due to the lensing effect by the inhomogeneous medium. This means in turn that a simplified analysis to find association of shell-shaped molecular gas around a round SNR edge, as often obtained about molecular-gas association (Sofue et al. 2021), may not be sufficient to study the past interaction with clouds.

## 6.2 Comment on the 'mixed morphology'

The term, 'mixed morphology' SNR, is often used to represent an SNRwith different morphologies observed at different wavelengths such as in radio and X rays, which originate from the outer shell and interior hot gas, respectively (Rho &amp; Petre 1998). Since the physical properties of the two components are quite different, the difference in their distributions is trivial. That means, the term may better be read as 'mixed physics', because the existence of such multiple conditions with different spatial distribution is trivial as the general property of a shock-wave driven shell by point explosion.

On the other hand, we have shown that even the radio (or optical) shell, which is the outermost surface of the SNR in contact with the ISM, shows partially different morphologies within the same SNR. Such "mixed morphologies" with the same physical condition within a single SNR shell are explained by interactions with interstellar clouds with different Alfvén velocities.

## 6.3 Different kinds of 'filament'

We have thus far shown that the filaments in SNRs are classified into two categories, either strings or sheets.

- · (a) Sheet filaments: Many of the observed 'filaments' are tangential view of corrugated (curved) sheets projected on the sky, but not real strings. The corrugation occurs by the reflection and refraction of the wave front by the lens effect of high- and low𝑉 interstellar fluctuations acting as concave and convex lenses, respectively. Traditionally, most of the shell-type SNRs have been recognized by their loop shapes of tangential-sheet filaments near their edges.
- · (b) String filaments A real string (line)-structured filament is formed in such a case that a low𝑉 tube or an elongated low𝑉 cloud acts as a cylindrical lens and the waves are converged onto a focal string (line). The straight filaments crossing the Cygnus Loop and the 𝜏 nebula are shown to be filaments of this type. It must be noted that the string filament exhibits higher amplification of the wave due to the focusing effect compared to the corrugated-sheet filaments. Anextreme case of wave amplification of this type is the focusing by a spherical low𝑉 cloud to the focal point, where the wave exhibits a small and concave spherical front with strongly amplified wave.
- · (c) Ring around hole The third type of filament is a focusing ring/loop around a hole on the shell surface produced by the diverging reflection of the wave front by a high𝑉 cloud. The tangled rings and arcs of filaments in the Vela nebula and S147 are of this type, which were produced by the diverging-lens effect through many high𝑉 clouds.
- · (d) Focus of converging front: Another type of filament is a dense cluster of converging waves due to focusing effect by a low𝑉 cloud. The north-eastern clump of high density condensation of filaments in S147 is a typical case of such focusing filaments.

These types of filaments have different formation mechanisms, and hence, different physical conditions.

Sheet filaments are the tangential view of wave fronts whose large line-of-sight depth makes the apparent brightness much higher than the surroundings.

On the other hand, the string filaments, as well as the converging clumps by convex lenses, are formed by focusing of the

(A) Cloud distribution around S147

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-d75f4891-ab07-4717-abfb-88eccfdd498a.jpg)

(B) Approximate result.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-73e56c0a-aeb3-41c2-89ab-95c2aea7820e.jpg)

(C) MHD simulation of S147

(D) Same, but in fits image

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-2684fae7-c9b1-4d3a-a756-9ed5a3252a54.jpg)

(E) S147 (North to top-left)

Figure 16. (A) Approximate cloud distribution and (B) an early simulation result for S147 (Sofue 1978). (C) Filamentary loops, arcs, and concave focusing clumps as the result of encounter with low- and high𝑉 clouds as marked by the bule and red dots, mimicking the SNR S147. (D) shows the same simulation, but by a fits-formated image. (E) Photo of S147 reproduced from https://apod.nasa.gov/apod/ap121009.html

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-280dd284-a66c-44e4-9e5c-6860c0134eaf.jpg)

(B) NE clump of S147

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-b2c48503-d28b-4e9c-8518-95cbc96793ff.jpg)

Figure 17. (A) Same as figure 16 for S147, but enlarged for the NE complex, which is a focusing clump by a big low𝑉 cloud, comapred with (B) the enlarged image of S147 (https://apod.nasa.gov/apod/ap121009.html)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-d3909559-d288-41ef-b311-9cf3d9a6a196.jpg)

waves onto a line or a point, so that the wave amplitude is increased by a factor proportional to the surface-area ratio between before and after the focusing, often amounting to orders of magnitudes. Such difference of real amplitudes would result in different emission mechanisms and physical conditions in the filaments. However, discussing such difference quantitatively is beyond the present simulation on the linear wave approximation. It would be a subject for the future work and observations.

corrugated sheets of the wave fronts. Real string-shaped filaments are formed by cylindrical focusing by cigar-like low 𝑉 clouds, and arcs and rings around holes are formed by reflection around high𝑉 clouds. The various types of "filaments" on SNR surfaces are, thus, understood as the manifestation of ISM fluctuations, through which the fronts have propagated, suffering from the lensing effect. However, since the front is treated as a linear MHD wave, such nonlinear quantities as the gas density, temperature, magnetic strength, cosmic-ray acceleration, and electromagnetic emissions cannot be calculated in the present method.

## 6.4 SNRs with unique morphology

The 'Manatee' nebula (W50) is an SNR having the most strange morphology. While efforts have been made to explain its elongated shape by interaction of a round shell and a precessing jet from the central compact object, SS433, (Zavala et al. 2008), the oblique stripes superposed on the surface have been difficult to be understood. In the present MHD model, however, such unique morphology has been relatively easily reproduced by assuming interaction with the ISM with cigar-shaped high𝑉 clouds superposed by sinusoidal variation of 𝑉 .

The 𝜏 nebula (N49 in LMC) is also unique for its orthogonal bunches of numerous long and straight filaments (Bilikova et al. 2007). This morphology is relatively simply explained in the present MHD model as due to interaction of the SNR surface with two orthogonal bar-shaped clouds having low𝑉 . The central long filaments across the Cygnus Loop is classified into the same category due to interaction with an elongated cloud having low𝑉 .

## 7 SUMMARY

By the simulation we showed qualitatively that most of the morphologies observed in shell-type SNRs can be reasonably explained by the lensing effect due to highand low𝑉 interstellar clouds. Most of the "filaments" in SNRs are shown to be tangential views of

The SNR morphology can be used to constrain the ISM conditions such as the wavelength of fluctuations, cloud sizes, direction of elongated clouds, and variation of relative amplitudes of the 𝑉 velocity. However, the observed shape is a result of combination of the reflection and refraction paths of the waves, the filaments are significantly displaced from the original position of the lenses (clouds), which makes it not easy to quantitatively determine the parameters of the clouds. Deconvolution of the observed SNR filaments back into the real ray paths as well as the distribution of Alfvén velocities in the ISM would be a subject for the future.

## ACKNOWLEDGEMENTS

The computations were carried out at the Astronomy Data Center of the National Astronomical Observatory of Japan.

## DATA AVAILABILITY

The images of the supernova remnants were downloaded mainly from the home page of APOD (Astronomy Picture of the Day) of NASA at https://apod.nasa.gov/apod/, and from the Hubble site at https://hubblesite.org/images. Some were taken from the web sites described in the figure captions for the individual images.

Figure 18. Gallery of some more examples of peculiar morphology of SNRs possibly influenced by the lensing effect by the inhomogeneity of ISM as compared with the MHD-wave simulation. In each block, the left and right panels show the simulation and observed image(s), respectively. (a) SNR N132-d from the Hubble web site and https://apod.nasa.gov/apod/ compared with the MHD model encountering a high𝑉 cloud. (b) Bilateral axisymmetric SNR G296.5+10.0 in radio (Giacani et al. 2000) compaed with MHD model of explosion in a high𝑉 plate or disc. (c) SNR 0509-67.5 from the Hubble site (https://stsciopo.org/STScI-01EVT1GHC7N7YEC39WH79DAN1S.jpg) and SN 1006 (https://hubblesite.org/contents/media/images/) compared with MHD model with slight corrugation on the shell. (d) SNR G18.8+0.38 from VGPS survey (Stil et al. 2006) and Puppis A (Dubner et al. 2013) compared with MHD model encountering a low𝑉 plate. (e) SNR CTB 1 from Hubble web site, NASA, compared with MHD model with sinusoidal variation of 𝑉 . (f) SNR W49b from Hubble web site compared with MHD model with high- and low𝑉 filaments. (g) SNR L136 from Hubble web site (https://apod.nasa.gov/apod/image/0801/) compared with MHD model divided by low𝑉 plate.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00260v1-figure-25af994d-9d33-4707-9f77-fab688314544.jpg)

## CONFLICT OF INTEREST

The author has no conflict of interest.

## REFERENCES

Abeysekara, A. U., Archer, A., Benbow, W., et al. 2020, ApJ, 894, 51. doi:10.3847/1538-4357/ab8310

Acciari, V. A., Aliu, E., Arlen, T., et al. 2009, ApJ.L, 698, L133. doi:10.1088/0004-637X/698/2/L133

Archambault, S., Archer, A., Benbow, W., et al. 2017, ApJ, 836, 23. doi:10.3847/1538-4357/836/1/23

Arias, M., Zhou, P., Chiotellis, A., et al. 2024, A&amp;Ap, 684, A178. doi:10.1051/0004-6361/202348064

Bilikova, J., Williams, R. N. M., Chu, Y.-H., et al. 2007, AJ, 134, 2308. doi:10.1086/52230

Bedogni, R. &amp; Woodward, P. R. 1990, A&amp;Ap, 231, 481

Bocchino, F., Maggio, A., Sciortino, S., et al. 2000, A&amp;Ap, 359, 316

Bocchino, F., Orlando, S., Miceli, M., et al. 2011, A&amp;Ap, 531, A129. doi:10.1051/0004-6361/201016341

Bykov, A. M., Krassilchtchikov, A. M., Uvarov, Y. A., et al. 2008, ApJ, 676, 1050. doi:10.1086/529117

Blair, W. P., Sankrit, R., &amp; Raymond, J. C. 2005, AJ, 129, 2268. doi:10.1086/429381

Cardillo, M., Tavani, M., Giuliani, A., et al. 2014, A&amp;Ap, 565, A74. doi:10.1051/0004-6361/201322685

Chevalier, R. A. 1999, ApJ, 511, 798. doi:10.1086/306710

Dubner, G. M., Holdaway, M., Goss, W. M., et al. 1998, AJ, 116, 1842. doi:10.1086/300537

| Dubner, G., Loiseau, N., Rodríguez-Pascual, P., et al. 2013, A&Ap, 555, A9. doi:10.1051/0004-6361/201321401                                                                                                    |
|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Fang, J., Yu, H., & Zhang, L. 2014, MNRAS, 445, 2484. doi:10.1093/mnras/stu1950                                                                                                                                |
| Fang, J., Yu, H., & Zhang, L. 2017, MNRAS, 464, 940. doi:10.1093/mnras/stw2473                                                                                                                                 |
| Farnes, J. S., Gaensler, B. M., Purcell, C., et al. 2017, MNRAS, 467, 4777. doi:10.1093/mnras/stx338                                                                                                           |
| Fujita, Y., Ohira, Y., Tanaka, S. J., et al. 2009, ApJ.L, 707, L179. doi:10.1088/0004-637X/707/2/L179                                                                                                          |
| Fesen, R. A., Weil, K. E., Cisneros, I. A., et al. 2018, MNRAS, 481, 1786. doi:10.1093/mnras/sty2370 Giuffrida, R., Miceli, M., Caprioli, D., et al. 2022, Nature Communications,                              |
| 13, 5098. doi:10.1038/s41467-022-32781-4 Giacani, E. B., Dubner, G. M., Green, A. J., et al. 2000, AJ, 119, 281. doi:10.1086/301173                                                                            |
| Green, A. J., Reeves, S. N., & Murphy, T. 2014, PASA, 31, e042. doi:10.1017/pasa.2014.37                                                                                                                       |
| Helfand, D. J., Becker, R. H., White, R. L., et al. 2006, AJ, 131, 2525. doi:10.1086/503253                                                                                                                    |
| Inoue, T., Yamazaki, R., Inutsuka, S.-. ichiro ., et al. 2012, ApJ, 744, 71. doi:10.1088/0004-637X/744/1/71 Klein, R. I., McKee, C. F., & Colella, P. 1994, ApJ, 420,                                          |
| 213. doi:10.1086/173554 Kothes, R., Fedotov, K., Foster, T. J., et al. 2006, A&Ap, 457, 1081.                                                                                                                  |
| doi:10.1051/0004-6361:20065062 Kupilas, M. M., Pittard, J. M., Wareing, C. J., et al. 2022, MNRAS,                                                                                                             |
| 513, 3345. doi:10.1093/mnras/stac1104 Lopez, L. A., Ramirez-Ruiz, E., Huppenkothen, D., et al. 2011, ApJ,                                                                                                      |
| 732, 114. doi:10.1088/0004-637X/732/2/114 Mac Low, M.-M., McKee, C. F., Klein, R. I., et al. 1994, ApJ, 433,                                                                                                   |
| 757. doi:10.1086/1                                                                                                                                                                                             |
| McKee, C. F. & Cowie, L. L. 1975, ApJ, 195, 715. doi:10.1086/153373.74685                                                                                                                                      |
| Miceli, M., Reale, F., Orlando, S., et al. 2006, A&Ap, 458, 213. doi:10.1051/0004-6361:20065738                                                                                                                |
| Miceli, M., Bamba, A., Orlando, S., et al. 2017, A&Ap, 599, A45. doi:10.1051/0004-6361/201629842                                                                                                               |
| Moranchel-Basurto, A., Velázquez, P. F., Giacani, E., et al. 2017, MNRAS, 472, 2117. doi:10.1093/mnras/stx2086 Nakamura, F., McKee, C. F., Klein, R. I., et al. 2006, ApJ Suppl., 164,                         |
| 477. doi:10.1086/501530 Ohmura, T., Ono, K., Sakemi, H., et al. 2021, ApJ, 910, 149.                                                                                                                           |
| doi:10.3847/1538-4357/abe5a1 Orlando, S., Wongwathanarat, A., Janka, H.-T., et al. 2022, A&Ap, 666, A2. doi:10.1051/0004-6361/202243258 Porter, L. J., Klimchuk, J. A., & Sturrock, P. A. 1994, ApJ, 435, 482. |
| doi:10.1086/174830 Rho, J. &Petre, R. 1998, ApJ.L, 503, L167. doi:10.1086/311538 Sakemi, H., Machida, M., Yamamoto, H., et al. 2023, PASJ, 75, 338.                                                            |
| doi:10.1093/pasj/psad001 Sano, H., Rowell, G., Reynoso, E. M., et al. 2019, ApJ, 876, doi:10.3847/1538-4357/ab108f                                                                                             |
| 37. Sano, H., Plucinsky, P. P., Bamba, A., et al. 2020, ApJ, 902, 53.                                                                                                                                          |
| doi:10.3847/1538-4357/abb469                                                                                                                                                                                   |
| Sano, H., Yoshiike, S., Yamane, Y., et al. 2021, ApJ, 919, 123.                                                                                                                                                |
| doi:10.3847/1538-4357/ac0dba Schneiter, E. M., de La Fuente, E., &Velázquez, P. F. 2006, MNRAS, 371, 369. doi:10.1111/j.1365-2966.2006.10652.x                                                                 |
| Shelton, R. L., Cox, D. P., Maciejewski, W., et al. 1999, ApJ, 524, 192. doi:10.1086/307799 Sofue, Y. 1978, A&Ap, 67, 409                                                                                      |
| Sofue, Y. 2023, MNRAS, 518, 6273. doi:10.1093/mnras/stac3137 Sofue, Y., Kohno, M., & Umemoto, T. 2021, ApJ Suppl.,                                                                                             |
| 253, doi:10.3847/1538-4365/abd484                                                                                                                                                                              |
| J. M., Taylor, A. R., Dickey, J. M., et al. 2006, AJ, 132,                                                                                                                                                     |
| doi:10.1086/505940                                                                                                                                                                                             |
| Stil,                                                                                                                                                                                                          |
| 1158.                                                                                                                                                                                                          |
| Sofue, Y. 2020, PASJ, 72, L4. doi:10.1093/pasj/psaa011                                                                                                                                                         |
| 17.                                                                                                                                                                                                            |

Stone, J. M. &amp; Norman, M. L. 1992, ApJ.L, 390, L17. doi:10.1086/186361 Tenorio-Tagle, G., Bodenheimer, P., &amp; Yorke, H. W. 1985, A&amp;Ap, 145, 70 Toledo-Roy, J. C., Esquivel, A., Velázquez, P. F., et al. 2014, MNRAS, 442,

229. doi:10.1093/mnras/stu880

- Uchida, Y. 1970, PASJ, 22, 341
- Uchida, Y. 1974, Sol. Phys., 39, 431. doi:10.1007/BF00162436
- Ustamujic, S., Orlando, S., Greco, E., et al. 2021, A&amp;Ap, 649, A14. doi:10.1051/0004-6361/202039940
- Velázquez, P. F., Dubner, G. M., Goss, W. M., et al. 2002, AJ, 124, 2145. doi:10.1086/342936
- Velázquez, P. F., Koenigsberger, G., &amp; Raga, A. C. 2003, ApJ, 584, 284. doi:10.1086/345663
- van den Bergh, S., Marscher, A. P., &amp; Terzian, Y. 1973, ApJ Suppl., 26, 19. doi:10.1086/190278
- West, J. L., Safi-Harb, S., Jaffe, T., et al. 2016, A&amp;Ap, 587, A148. doi:10.1051/0004-6361/201527001
- Zavala, J., Velázquez, P. F., Cerqueira, A. H., et al. 2008, MNRAS, 387, 839. doi:10.1111/j.1365-2966.2008.13281.x
- Zhang, D. &amp; Chevalier, R. A. 2019, MNRAS, 482, 1602. doi:10.1093/mnras/sty2769
- Zhou, X., Shen, Y., Yuan, D., et al. 2024, Nature Communications, 15, 3281. doi:10.1038/s41467-024-46846-z
Yoshiaki Sofue
MHD lensing in inhomogeneous ISM for qualitative understanding of the morphology of supernova remnants
Morphological evolution of expanding shells of fast-mode magnetohydrodynamic
(MHD) waves through an inhomogeneous ISM is investigated in order to
qualitatively understand the complicated morphology of shell-type supernova
remnants (SNR). Interstellar clouds with high Alfv\'en velocity act as concave
lenses to diverge the MHD waves, while those with slow Alfv\'en velocity act as
convex lenses to converge the waves to the focal points. By combination of
various types of clouds and fluctuations with different Alfv\'en velocities,
sizes, or wavelengths, the MHD-wave shells attain various morphological
structures, exhibiting filaments, arcs, loops, holes, and focal strings,
mimicking old and deformed SNRs.