Resolved magnetohydrodynamic wave lensing in the solar corona

Electromagnetic wave lensing, a common physical phenomenon recognized in visible light for centuries, finds extensive applications in manipulating light in optical systems such as telescopes and cameras. Magnetohydrodynamic wave is a common perturbation phenomenon in the corona. By using high spatio-temporal resolution observations from the Solar Dynamics Observatory, here, we report the observation of a magnetohydrodynamic wave lensing in the highly ionized and magnetized coronal plasma, where quasi-periodic wavefronts emanated from a flare converged at a specific point after traversing a coronal hole. The entire process resembles an electromagnetic wave lensing from the source to the focus. Meanwhile, the magnetohydrodynamic wave lensing is well reproduced through a magnetohydrodynamic numerical simulation with full spatio-temporal resolution. We further investigate potential applications for coronal seismology, as the lensing process encodes information on the Alfvén speed, in conjunction with favorable geometric and density variations.


Introduction
Lensing of various wave phenomena, such as electromagnetic (EM) (light) or sound waves, is ubiquitous.Similar to a convex glass lens that focuses light, acoustic lenses can focus and manipulate acoustic waves and have a wide range of valuable applications, including biomedical engineering, imaging, cancer treatment, and surgery [1][2][3][4] .Acoustic lensing is responsible for the Rotunda effect (sometimes called whispering galleries) 5 , where a conversation can be held on opposite sides of a round-shaped auditorium.An intuitive lensing effect in daily life also occurs in water waves: shallower water areas in the ocean can act as lenses.Their amplitude and energy will be enhanced near the focal point when they pass through these areas 6,7 .Light maneuvering with optical lenses has been practiced for many centuries, which has evolved into a mature industry for lens manufacturing, such as cameras, telescopes, microscopes, and lasers.In general relativity, a ray of light follows the curvature of space-time, and its path would be bent around a massive object.This effect, known as gravitational lensing, could be used for detecting black holes and dark matter in the Universe 8,9 .
Magnetohydrodynamic (MHD) waves play essential roles in various fundamental processes in the corona.For example, they might heat the corona plasma 10,11 and accelerate the solar wind 12 .Moreover, through a technique known as coronal seismology 13 , physical information carried by MHD waves can be utilized to decode coronal physical parameters (such as the strength of the magnetic field) that can not be measured directly by our current technology.Fast-mode magnetoacoustic waves are a type of MHD waves that perturb coronal magnetic field and plasma density simultaneously, enabling its detection in the extreme ultraviolet (EUV) bandpass.Although typical wave phenomena such as refraction, transmissions, and reflection in the corona have been reported in previous studies [14][15][16] , the lensing effect of fast-mode MHD waves has not been directly observed.
Here, we present comprehensive spatio-temporal resolved observations of a globally propagating fast-mode MHD wave that exhibits an unambiguous lensing process upon its passage through a coronal hole (CH) that acts as a convex lens.The entire lensing process of the wave, including the transmission through the lens and the focusing at the focal point, is comprehensively captured by the Atmospheric Imaging Assembly (AIA) 17 instrument onboard the Solar Dynamics Observatory (SDO) 18 with full spatio-temporal resolution.Through numerical simulations, we not only reproduce the observed lensing effect but also conduct a comprehensive investigation of the essential physical parameters that govern the occurrence of lensing effects on MHD waves in structured plasma environments.

Results
On 2011 February 24, an energetic GOES M3.5 solar flare occurred in NOAA active region AR11163.From the viewpoint of Solar TErrestrial RElations Observatory Behind (STEREO-B) 19 , we can see a crescent-shaped CH adjacent to the active region, as shown in 195 Å image captured by the Extreme UltraViolet Imager (EUVI) 20 onboard the STEREO (Figure 1(a)).The flare excited large-scale, quasi-periodic intensity perturbations propagating along the solar surface.These largescale perturbations, a form of MHD waves [21][22][23][24] , showed a chain of arc-shaped wavefronts, which remained approximately concentric with the flare center during the initial phase (Figure 1(d)).Subsquently, the wave train propagated towards the center of the solar disk (Figure 1(e)) and transmitted through a low latitude CH (the dark region in Figure 1(a)-(b)).Following, the transmission, the wave train entered a quiet-Sun region with a speed of about 350 km s −1 .Such a speed is consis-tent with the typical Alfvén speed of the quiet corona 25 .Interestingly, as the wavefronts propagated westward after crossing the CH's far-side boundary, the original arc-shaped wavefronts gradually became anti-arc-shaped and converged towards the focal point far from the CH (Figure 1(f)).The complete evolution of the MHD wave is display in Figure 1(c).More details can be found in Supplementary Video I.These observational characteristics strongly suggest that the observed wave train is a propagating quasi-periodic fast-mode MHD wave in nature, and its physical property is similar to the globally propagating EUV waves [26][27][28][29] .In addition, the whole evolution process of the wave train from its source to the focusing resembles an EM wave passing through a converging lens; or mimics acoustic wave refocusing, but this time dictated by geometric variation of the Alfvén speed.In the present case, the flare center, the CH, and the final focal point make up the key components of the MHD lensing system.The lensing effect of the wave train could be detected clearly in AIA 211 Å (2.0 MK), 193 Å (1.6 MK), and 171 Å (0.6 MK) channels.Here, we only present the observations of the AIA 193 Å channel.
A CH is characterized by dark regions in X-ray and EUV images, where a single magnetic polarity dominates 30 .At the same time, the densities and temperatures are lower than in the CH surroundings [31][32][33] .Therefore, a CH is often a region of high Alfvén speed with respect to the quiet-Sun.Theoretical investigation suggests that large-scale coronal MHD waves will be refracted away from regions of high Alfvén speed towards the low Alfvén speed areas.Meanwhile, as a linear fast-mode MHD wave enters the low Alfvén speed region, it may steepen into a shock wave due to refraction and the deceleration of the wave 26,34 .For an optical lensing system, the light speed is usually reduced significantly inside the lens with respect to that outside the lens.However, for the MHD lensing system presented here, the fast-mode MHD wave propagated faster inside the CH than in the quiet-Sun region; it also showed an obvious focusing effect after its transmission through the CH due to the special shape of the CH that acts as the MHD lens.The behavior of the observed wave train is hence consistent with theoretical prediction 34 .The particular crescent geometry and enhanced Alfvén speed of the CH are important for the occurrence of the observed lensing effect.
To quantify the MHD lensing effect, we measure the intensity amplitude and energy flux density of the wave train along its propagation path.Although the STEREO-B had an advantageous viewing angle to observe the eruption source region, its spatial and temporal resolutions are too low to resolve the lensing process.Luckily, the SDO/AIA images are sufficient to determine the detailed evolution process of the wave, though its viewing angle was not ideal for measuring the wave signal before and during the transmission of the wave train through the CH.As shown in Figure 1(b), the eruption source region was on the backside of the Sun, and the wave signal inside the CH was too weak to measure.Thus, we only measure the intensity amplitude variation after the wave train passed through the CH using the AIA images, i.e., during the focusing phase.Figure 2 shows the intensity amplitude variations along three selected paths (Figure 1(c)).As the wave approached the focal point, the intensity amplitude increases by about 2 to 6 times along all three paths.The intensity amplitudes reached their peaks around the focal point, and then they started dropping off quickly to the level before the focusing.We calculate the energy flux density of the wave to estimate the energy-focusing effect (see Methods section).The results indicate that energy flux densities of the wave were respectively about 5.2 × 10 4 erg • cm −2 • s −1 and 3.5 × 10 5 erg • cm −2 • s −1 before the focusing and at the focal point, indicating that the wave energy flux density at the focal point increased about seven times compared to the pre-focusing level.Therefore, the results suggest that an MHD lens can effectively focus wave energy on the focal point.Since magnetized plasma and magnetic field pervade the whole Universe and the ideal MHD description is scale-invariant, such MHD lensing is expected to be duplicated in magnetized planetary, stellar, and galactic counterparts and may play a significant role in the convergence of MHD wave energy.
As mentioned above, the observation suggests that a particular structuring of magnetized plasma could act as an MHD lens that maneuvers an MHD wave and the associated energy flux.We model this scenario in a 2.5-dimensional magnetized plasma to further demonstrate the MHD lensing effect.The essential aspects of the lensing effect of a fast-mode MHD wave can be investigated in a simple idealized setup, where the influence of the model parameters can be shown clearly.This part is done in Methods section, where the numerical approach and setup details are also given.Here, we discuss a more realistic scenario that uses the observed geometric shape of the CH.
As shown in Figure 3, the simulation ignores the solar surface's curvature and takes the CH's shape from the STEREO-B/EUVI 195 Å image as shown in Figure 1(a).The location of active region AR11163 relative to the CH is also determined by this image.We launch a point-source induced fast-mode MHD wave from the center of the active region, which shows as a circular wavefront and travels at the local Alfvén speed away from the source region.The rightward edge of the induced wavefront gets reflected and transmitted at the left edge of the CH.Eventually, some of the transmitted waves converge to a particular point.More details can be found in Supplementary Video II.Transmitted wave patterns traverse the CH at the locally higher Alfvén speed, in direct correspondence to the adopted density (and temperature) contrast used in the model.The background magnetic field was set to be orthogonal to the simulation plane (i.e., a vertical, open field is assumed, as appropriate for the CH region).The sum of gas and magnetic pressures in the CH region was set in balance with the surrounding coronal plasma, where input parameters control the external plasma beta (here β = 0.05), and a prescribed density and temperature contrast reflects the observed lower T and ρ conditions in the CH.After using reasonable density and temperature ratios, wave amplitudes, and source duration, the CH shape dictates the wave behavior throughout the simulated coronal region.We find that the complex shape of the CH causes a single circular wavefront to become reflected/transmitted in an intricate pattern of constructively and destructively interfering waves that sense the changes in background quantities at the CH edges.After the first transmitted wave has wholly traversed the CH, the concave part found at the opposing side of the wave source region causes an intricate pattern of constructively interfering waves in a location similar to the observations that find the lensing to occur.We can then isolate the main ingredients of an Alfvén-speed-lensing in a different idealized setup by simplifying the CH shape to be composed of a straight left edge with a concave parabolic right edge.This different idealized model is then parametrically explored in Methods section (this has the advantage that the complex wave interference patterns due to the complex CH shape boundary are omitted).There, we find that the density ratio determines the precise shape of the lensing at the focal point, and the overall plasma beta sets the time of traversal from source to focus.In the idealized setup, after passing the artificial MHD lens (or CH), the amplitude increases as it approaches the focal point and then becomes weaker as it moves away from the focal point.The simulation's wave propagation and focusing phase are consistent with our observation, implying that observing a lensing event allows us to deduce local Alfvén speed, density variations, and plasma beta conditions.At the same time, the needed geometric information can be directly seen in EUV images.This is of immediate use for coronal seismology.

Discussion
This paper reports an observation of the MHD lensing effect with full spatial and temporal resolution.A quasi-periodic fast-mode MHD magnetoacoustic wave was excited by a solar flare and transmitted through a structured CH with low density.Due to the particular shape and the higher Alfvén speed of the CH, it acted as an MHD converging lens and focused the wave energy towards a small region beyond the CH, although the focussing is not total.Thanks to the full disk and high spatio-temporal resolution observations taken by SDO/AIA, we successfully observed the dynamical lensing effect of MHD waves in solar, magnetized plasma conditions.We demonstrated that the lensing of a fast-mode magnetoacoustic wave conducted by a CH follows our numerical experiment, where the main parameters influencing the lensing process can be identified in a simple idealized setup.We note that deflection, reflection, and transmission effects of globally propagating EUV waves have been observed and confirmed by many numerical simulations during their interaction with coronal structures such as active regions and CHs 14-16, 28, 35-39 .However, none of these works considered the influence of the shape of CHs on wave propagation.In the present case, the appearance of the lensing effect is believed to be due to the sharp gradients of the temperature, background plasma density, and magnetic field strength at the CH boundary, as well as the particular shape of the CH, which could be explained using the method of geometrical acoustics 26 .Another focus of coronal fast-mode MHD waves appears at magnetic null points that commonly exist in breakout 40,41 and fan-spine 42,43 magnetic systems.Since the magnetic field at a coronal null point approaches zero, the low Alfvén speed property around the null point will cause a fast-mode MHD wave to focus towards and wrap around the null point when the wave approaches the null 44 .This convergence of an MHD wave at a coronal null point differs from the focusing of an MHD wave at a non-zero Alfvén speed point, as is happening in the MHD lensing effect discussed here, though they are both due to refraction effects.In the first case, the wave and its energy accumulate at the zero point, resulting in a convergence effect due to the inability of the wave to pass through the magnetic zero point 45 .The focusing effect reported in the present paper is a manifestation of magnetoacoustic focusing in the solar coronal plasma environment and can be explained in terms of classical geometric acoustics 26 .
The MHD lensing effect documented in this study must have counterparts across various length scales since the MHD description of plasma dynamics is scale invariant 46 .A localized and specially shaped MHD lensing setup could focus MHD wave energy to a specific destination.This process can be reproduced and experimented with in laboratory plasmas.The MHD lens could be designed to either converge or diverge MHD wave energy.For example, the magnetosphere of the Earth diverts solar storms from the dayside to the magnetotail.Therefore, it can be viewed as a natural diverging MHD lens that effectively protects life on the Earth 47 .A converging MHD lens is observed and modeled, showing clear diagnostic information on the background plasma.We illustrated how the focus relies on low plasma beta conditions and specific shapes of CHs and encodes their density contrast and Alfvén speed variations.Since fast-mode MHD waves can be easily excited in space and laboratory plasma, the MHD lensing effect could be used for diagnosing various plasma properties, as initially proposed for solar coronal Moreton waves 26 .

Observations and Data Processing
This study used imaging observations provided by AIA.The AIA takes full-disk corona images at seven EUV wavelength bands that are sensitive to plasma emissions with a wide temperature coverage ranging from about 500,000 K to over 20,000,000 K. Its time resolution is 12 seconds, and each pixel corresponds to an angular width of 0.6 arc second or about 420 km on the Sun.In this study, we focus on the 193 Å bandpass (Fe XII line) that is mainly from plasma emission at about 1,000,000 K.We also used images taken by the EUVI onboard the STEREO.EUVI observes the corona with an angle of about 95 • apart from AIA with respect to the center of the Sun.In the EUVI's viewpoint, AR11163 and the CH were observed with negligible projection effects (Figure 1(a)-(b)).Both AIA and EUVI images were calibrated with standard routines in the Solar Software (SSW).This process calibrates the images by removing the CCD bias and dark frames, correcting for the flat field, and then normalizing them with the exposure time.

Kinematic Study and Energy Estimation
In order to highlight the wavefronts, we calculated the forward running difference of the image sequence in the time dimension, as shown in Figure 1(d)-(f).The wavefronts were tracked manually to measure wave parameters, as illustrated in Figure 1(c).The emission intensities along these slits were stacked in time to form time-distance plots.Figure 4 shows the wavefronts propagation in the form of base-difference time-distance maps, in which the slopes of the ridges represent the speeds.In Figure 4, we can see that the speed of the wave along Slice 2 (about 312 km s −1 ) is smaller than Slice 1 (about 368 km s −1 ) and Slice 3 (about 385 km s −1 ) after its passage through the CH, since the paths of the latter two are longer.Additionally, the wavefront showed a visible intensity enhancement at the focal point as shown in Figure 4(b).The average phase speed is 350 km s −1 (with projection effect), which is in agreement with the typical fast-mode MHD wave speed in the solar corona 25 .The black horizontal lines in Figure 4 labeled p1-p10 indicate the positions for extracting the amplitudes (emission intensities) of the wave.After obtaining the value of the wave crest at each position, we get the amplitude variations along the three paths, as shown in Figure 2. The uncertainties on the wave amplitude were estimated mainly by accounting for the photon noise, as suggested by Yuan et al. 48nsidering that the group velocity of the wave is more difficult to obtain, we do not use the equation for calculating the energy flow using the group velocity suggested elsewhere 49 but employed the following equation for estimating the energy flow using the phase velocity: where ρ is the plasma density, δv is the perturbation speed of the local plasma, and v ph is the phase speed of the wave 50,51 .The relative perturbation speed could be estimated by δv/v ph ≥ δρ/ρ = δI/(2I), where δρ and δI are the perturbations to density and emission intensity.We assumed that the emission intensity of the AIA 193 Å bandpass is proportional to the square of density (and hence the square of electron number density 52,53 ), I ∝ ρ 2 ∝ n 2 e .Then Equation 1 could be rearranged for practical use as done in Ref 54 (2) δI/I was measured to be about 12% at the right boundary of the CH, while it was about 31% at the focal point.If we use the average phase speed of v ph = 350 km s −1 and number density n e = 3.4 × 10 8 cm −3 (an empirical value 55 ), then we obtain that the energy flux of the transmitted wave through the CH was about 5.2 × 10 4 erg • cm −2 • s −1 , while the energy flux at the focal point was about 3.5 × 10 5 erg • cm −2 • s −1 .This means that the energy flux of the wave increased by about 6.7 times after the focusing.

Numerical model and analysis
To model a fully ionized and magnetized plasma in macroscopic scales, we used the Newtonian MHD equations, a scale-invariant model for plasma dynamics 46 .They combine the pre-Maxwell equations of electromagnetism and the equations of fluid dynamics as follows, where ρ, p, v, and B are the mass density, the gas pressure, the fluid velocity, and the magnetic field, respectively.Here H = [ρε + ρv 2 /2 + B 2 /8π] is the total energy density, where ε = p/(γρ − ρ) is the internal energy per unit mass.We adopted γ = 5/3 as the adiabatic index.This set of MHD equations is closed by the equation of state p = ρk B T /µm H , where m H is the hydrogen atom mass, and the mean molecular weight is taken as µ = 1.4/2.3 for a fully ionized plasma with a 10 : 1 abundance of hydrogen and helium.
We used the open-source MPI-AMRVAC code [56][57][58][59] to solve the MHD equations in a conservative form with a finite-volume scheme.In our setup, we used an HLLD solver with a Koren-type limiter for reconstruction and a three-step Runge-Kutta method in the time integration.Equations 3-7 were solved in a dimensionless form.The normalization uses the reference quantities L 0 , ρ 0 , and B 0 etc., as given in Table 1.
We modeled a CH surrounded by a uniform and static coronal plasma.The essential features of fast-mode lensing can be investigated both in a simple idealized setup and a more realistic scenario that uses the observed geometric shape of the CH.The latter is already shown in Figure 3. Here, we describe the numerical setups and present the simulation results of the idealized setup in detail.
For the idealized setup, the simulation domain includes a low-density (high Alfvén speed) region to act as an artificial MHD lens.This region is equivalent to a CH, and it is bounded by a straight line x/L 0 = −0.2 at the left edge and a parabolic curve x/L 0 = 4y 2 (see Figure 5) at the right edge.The background coronal plasma has a density of ρ bg /ρ 0 = 1 and a pressure of p bg /p 0 = 1/γ.The CH was set to have a cooler temperature and lower density than the background.The density and temperature ratios between the CH and background plasma are R ρ = 0.25 and R T = 0.1, respectively.The magnetic field was set normal to the 2D simulation plane; therefore, it only has a z-component initially.The background magnetic field strength was set to B bg /B 0 = 2p bg /(p 0 β bg ), where β bg = 0.05 is the plasma beta for the background coronal plasma, i.e. the ratio of the gas pressure to the magnetic pressure.The magnetic field of the CH was calculated by setting the total plasma pressure (p tot = p + B 2 /8π) in equilibrium over the simulation domain.
This simulation was done in a two-dimensional plasma region of −2 ≤ x/L 0 ≤ 2 and −2 ≤ y/L 0 ≤ 2, where L 0 is our normalization factor for length, see Table 1.The region has a base computational grid of 100 × 100 cells, but an adaptive mesh refinement (AMR) was set up to allow for five levels.Therefore, the finest grid reaches a physical cell size of 700 km.At the lateral boundaries, the z-component of the velocity v z and the x and y-components of the magnetic field B x and B y are fixed to zero while a continuous extrapolation applies to other variables.
A periodic fast-mode MHD wave is excited by setting a perturbation source located at (x s /L 0 , y s /L 0 ) = (−0.65,0) in the following form, where b .This, in essence, corresponds to the exact analytic eigenfunction variation for a fast (linear) magnetoacoustic wave 46 with frequency ω f .Geometric factors are centered on the source location as in cos θ = (x − x s )/r 0 , sin θ = (y − y s )/r 0 , where r 0 = (x − x s ) 2 + (y − y s ) 2 and r 0 < 0.05L 0 .
The wavefronts propagate outward radially in the form of concentric rings; partial reflection and transmission occur at the left boundary of the MHD lens.When the MHD wave passes throingugh the MHD lens, it starts to bend towards a focal point.The convergence process is similar to the observation, i.e., the focusing is not total, which is caused by the shape of the CH.More details can be found in Supplementary Video III.The focal length is clearly related to the geometry and structuring of the MHD lens.Along the selected path, we first extract the slices from the total energy density δH/p 0 images (here the δH is the difference between total energy density H and the background total energy density H 0 to remove the background and retain the wavefront signal.),and subsequently stack these slices according to the time series to obtain the time-distance plots, as shown in Figure 6(a), (b), and (c).Based on these obtained time-distance plots of the wavefront evolution, we then calculated the wave amplitude variation following three paths (Slice 4, Slice 5 and Slice 6 in Figure 5(b)) starting from the source to the focal point.As shown in Figure 6(d), the wave amplitude decreases as the wave departs from the source.After passing the artificial MHD lens (or CH), the amplitude increases as it approaches the focal point and then becomes weaker as it moves away from the focal point.The wave transmission and focusing phase in the simulation are consistent with our observation.
We also investigated the essential aspects of the MHD lensing by varying the various parameters, where the influence of the model parameters can be shown clearly.By doing so, we only adjusted one parameter in each case while keeping all other parameters unchanged.We first investigated how the MHD lensing effect changes with various density ratios R ρ between the CH and background plasma, as shown in Figure 7.The MHD lensing effect increases as R ρ decreases.When R ρ = 0.1, the transmitted wave is so strongly bent by the CH that a single focal point cannot be formed.The MHD lensing focal points are obvious to identify when 0.2 < R ρ < 0.5.When R ρ is increased towards 1, the MHD lensing effect gradually disappears because the refractive indices of the CH and the background are almost identical.Then, we investigated how the MHD lensing effect changes with various plasma beta β bg for the background coronal plasma as shown in Figure 8.The plasma beta is inversely proportional to the square of the Alfvén speed.Thus, the Alfvén wave propagation time is scaled with the square root of plasma beta as shown in Figure 8, where the time on the right is twice of the left and the plasma beta is 4 times.We also investigated how the position x s of the wave source changes the MHD lensing effect, as shown in Figure 9.The focal point of the transmitted wave shifts in the same direction to which the source position x s moves horizontally.We also discovered that changing R T does not have a major influence, nor does changing the amplitude P b .
In the main paper, we present a simulation that used the observed realistic geometric shape of the CH, where the numerical setup is the same as the idealized simple setup, except for the shape of the CH.The values of normalization factors p 0 , t 0 , and L 0 for the total energy density (H), time (t) and length (x, y, z) are listed in Table 1.The simulation domain includes a background coronal plasma (light blue) and an artificial MHD lens (blue) enclosed by a straight line and a parabolic curve.The values of normalization factors p 0 , t 0 , L 0 for the total energy density(H), time (t) and length (x, y, z) are listed in Table 1.  1.  1.

Figure 1 :
Figure 1: Overview of the solar coronal structure and MHD wave propagation and focusing process.(a)-(b) Combined views of STEREO-B/EUVI 195 Å and SDO/AIA 193 Å images (about 95 degrees apart) showing the full disk of the Sun at 07:30 UT on 2011 February 24.The red circle and white closed curve in (a) marked the active region AR11163 and the contours of the interest CH, respectively.The contour of the CH as deduced from STEREO-B/EUVI 195 Å is also projected onto the SDO/AIA 193 Å image.The open green curves in (b) show the extrapolated magnetic field line using the potential-field source-surface (PFSS) model 60 , while the white box displays the field of view of (c)-(f).(c) Synoptic view showing the evolution of the MHD wavefront.The CH is marked by a gray-shade, and the wavefront is colored with time.The focal point at (x, y) = [−500 ′′ , 150 ′′ ] is labeled with a black plus sign.Three slits (Slice 1-3) are used to make time-distance plots as shown in Figure 4 to track the evolution of wave.(d)-(f) SDO/AIA 193 Å running-difference images showing the flaring core and the MHD wave evolution, where the evolving wavefronts are labeled with colored curves.(An animation Video I of this figure is available.)

Figure 2 :
Figure2: Intensity amplitude evolution of the MHD wave along three slits indicated by Slice 1, Slice 2, and Slice 3 in Figure1(c).The data sampling locations correspond to the positions marked by the white horizontal lines p1-p10 in the time-distance plots of Figure4.The colored shades give the errors obtained from the calculation method provided by Yuan et al48 for wave amplitude in units of data number (DN).The location of the CH was depicted in light brown color, while the brown area represents regions that were not measured due to overexposure.The vertical dashed line marks the position of the focal point.

Figure 3 :Figure 4 :
Figure3: Numerical simulation of the MHD lensing process at t/t 0 = 0.185 based on the observed geometric shape of the CH.The simulation domain includes a background coronal plasma (light blue) and the CH that acts as an MHD lens (blue), as shown in the total energy density (H/p 0 ).A single-period fast-mode MHD wave is launched to the left of the which propagates with a circular wavefront initially.After being reflected/transmitted, the initial circular wavefront presents a complex pattern that senses the geometric shape at the CH edges.The initial wavefront, the reflected wavefront by the CH, the transmitted wave, and the focal point are annotated by black arrows.The values of normalization factors p 0 , t 0 , L 0 for the total energy density (H), time (t) and length (x, y, z) are listed in Table1.The red dot represents the location of the wave source, and it retains the same meaning in the subsequent Figures.(An animation Video II of this figure is available.)

Figure 5 :Figure 6 :
Figure 5: Evolution of perturbation to the total energy density (H/p 0 ) during the MHD lensing process in an idealized numerical simulation.(a) Numerical set-up and MHD wave excitation.The simulation domain includes a background coronal plasma (light blue) and an artificial MHD lens (blue) enclosed by a straight line and a parabolic curve.A fast-mode MHD wave train launched to the left of the MHD lens propagates outwards in the form of concentric rings.(b) Initial phase of MHD wave propagation and its interaction with the MHD lens.(c) and (d) MHD wave focusing process.The black arrows point out the initial wave, the reflected wave and the location of the focus.Three paths (Slice 4, Slice 5, and Slice 6) in (b) are used for making the time-distance plots as shown in Figure 6.The values of normalization factors p 0 , t 0 , L 0 for the total energy density (H), time (t) and length (x, y, z) are listed in Table 1.(An animation Video III of this figure is available.)

Figure 7 :
Figure 7: MHD lensing effects with various density ratios R ρ between the CH and background plasma (a)-(d) MHD lensing effect at R ρ = 0.1, R ρ = 0.25, R ρ = 0.5, and R ρ = 0.8, respectively.The simulation domain includes a background coronal plasma (light blue) and an artificial MHD lens (blue) enclosed by a straight line and a parabolic curve.The values of normalization factors p 0 , t 0 , L 0 for the total energy density(H), time (t) and length (x, y, z) are listed in Table1.

Figure 8 :
Figure 8: MHD lensing effects with various plasma beta β bg for the background coronal plasma.(a)-(b) The MHD lensing effect when the plasma beta β bg is 0.05 and 0.2.The simulation domain includes a background coronal plasma (light blue or light orange) and an artificial MHD lens (blue or deep blue) enclosed by a straight line and a parabolic curve.The values of normalization factors p 0 , t 0 , L 0 for the total energy density(H), time (t) and length (x, y, z) are listed in Table1.

Figure 9 :
Figure 9: MHD lensing effects with different positions x s of wave source.(a)-(b)The MHD lensing effect when the relative position x s /L 0 of the wave source is taken as -1.0 and as -0.65.The simulation domain includes a background coronal plasma (light blue) and an artificial MHD lens (blue) enclosed by a straight line and a parabolic curve.The values of normalization factors p 0 , t 0 , L 0 for the total energy density(H), time (t) and length (x, y, z) are listed in Table1.

Table 1 :
Quantities used for normalization