## Universal scaling behavior of resonant absorption

## D. J. Yu ∗

School of Space Research, Kyung Hee Uiversity 1732, Deogyeong-daero, Yongin, Gyeonggi, 17104, Korea

## T. Van Doorsselaere

Centre for mathematical Plasma Astrophysics, Department of Mathematics, KU Leuven Celestijnenlaan 200B bus 2400, B-3001 Leuven, Belgium

Mode conversion and resonant absorption are crucial mechanisms for wave transport and absorption. Scaling behavior of mode conversion or resonant absorption is well-known for electromagnetic and MHD waves in planar geometry. Our recent study showed that such a scaling behavior of resonant absorption could also exist for coronal loop oscillations with cylindrical geometry, but it was only tested for one density profile. Here we generalise our previous study on the scaling behavior of resonant absorption by considering multiple density profiles. Applying an invariant imbedding method to the ideal MHD wave equations, we show that the scaling behavior also exists for these density models. We thus generalise our earlier results and show that such a universal scaling exists in cylindrical geometry, too. Given these results and the earlier results in planar geometry, we formulate a hypothesis that a universal scaling behavior exists regardless of the type of mode conversion or resonant absorption.

Mode conversion and resonant absorption can often occur in inhomogeneous plasmas such that two different waves are coupled with an exchange of energy. This mechanism has been considered as crucial for plasma heating and transport [4, 5, 14-17, 23, 28-31, 35]. The previous literature is mainly concerned with mode conversion between electromagnetic waves and electrostatic waves or between magnetoacoustic waves and Alfv´ en waves [e.g., 23, 44]. From theoretical and numerical studies, it is well established that a scaling behavior exists for the mode conversion coefficient in a planar geometry [8, 10-13, 17, 19, 20, 23, 27, 33, 38, 43]. The scaling parameter has in general a certain relation with the density gradient, wave vector, and wave frequency [see, e.g., 17, 19, 20, 33, 38]. The mode conversion coefficient, which denotes the ratio of energy transformation from one mode to another mode, has a maximum value of about 0.5 for the linear or parabolic density profile. When the density profile is more complex, the scaling behavior is broken and the conversion coefficient can vary greatly [39-42]. While this scaling behavior is well-known in planar geometry, its behavior in cylindrical geometry was unknown until recently. In our recent study, we considered kink modes in coronal loops, where the coronal loop was assumed as a infinitely long, straight, and axisymmetric cylinder with a radial inhomogeneity. We found that resonant absorption of the kink modes in the Alfv´ en resonance also show a certain scaling behavior. We called this behavior sub-universal therein. However, we only considered a sinusoidal density profile in the transitional layer [44]. Given that it was studied only for a particular density profile, it was unclear whether this behavior would be generally valid, as it is for planar geometry. This motivates the question if the results from our

∗ djyu79@gmail.com

previous paper can be generalised to other density profiles. If they can be generalised, it would imply that the scaling behavior is universal in inhomogeneous plasma, regardless of the geometry, density profile, or mode conversion.

In this work, to answer this question if the scaling law also occurs for other density profiles, we extend our previous study on resonant absorption by comparing linear, sinusoidal, and parabolic density profiles. Adding to our previous results on the sinusoidal profile, we show that this kind of scaling behavior exists also for linear and parabolic density profiles. Thus, we show that the scaling behavior also exists in cylindrical geometry, for monotonic density profiles. Based on these results, we speculate that a scaling feature of mode conversion (or resonant absorption) can generally occur regardless of the geometry, by synthesizing the newly-found scaling behaviors in cylindrical geometry with the previous ones in planar geometry.

The governing equation for the wave amplitude that describes the mode conversion usually contains a singularity in it, which makes it difficult to get exact solutions. The invariant imbedding method (IIM) transforms the concerned differential equations with boundary conditions into differential equations with initial conditions. It is an efficient and powerful tool for this kind of problems when one considers a one-dimensional inhomogeneity in the wave equations [2, 7, 18, 19, 22, 24, 25, 32]. We apply the recently used invariant imbedding method [see 44, and references therein] to resonant absorption of externally-driven standing kink modes in the Alfv´n cone tinuum.

From the ideal MHD equations, the wave equation for the perturbed total pressure P under coronal conditions is written as, by assuming a dependence exp[ ( i mφ + k z z -

ωt )] in cylindrical coordinates, [e.g., 34, 44]

$$\frac { d ^ { 2 } P } { d r ^ { 2 } } + \stackrel { \mu } { \left [ \frac { 1 } { r } - \frac { \omega ^ { 2 } } { \rho _ { 0 } ( \omega ^ { 2 } - \omega _ { A } ^ { 2 } ) } } \frac { d \rho _ { 0 } } { d r } \right ] \frac { u } { d r } + \stackrel { \mu } { \left [ \frac { \omega ^ { 2 } } { v _ { A } ^ { 2 } } - k _ { z } ^ { 2 } - \frac { m ^ { 2 } } { r ^ { 2 } } \right ] } P = 0, \quad \begin{array} { c } \omega / c _ { 0 }, \\ \text{in vao} \\ r > I \\ \text{and } F \end{array}$$

where P ≈ B · b /µ 0 = B b /µ 0 z 0 , b z is the perturbed magnetic field, the background magnetic field is B = (0 , 0 , B 0 ), ω 2 A ( r ) = k v 2 z 2 A = k B /ρ 2 z 2 0 0 ( r µ ) 0 , k z is the longitudinal wave number, m is the azimuthal wavenumber, and µ 0 is the magnetic permeability.

We assume that the density varies from ρ i inside the loop to ρ e outside the loop through the transitional layer of thickness l ( ρ t ):

$$\rho _ { 0 } ( r ) = \begin{cases} \stackrel { \bullet } { \rho _ { i } } & \text{if $r\leq R_{0}-l/2$} & \text{$\omega=$} \\ \rho _ { t } ( r ) & \text{if $R_{0}-l/2<r<R_{0}+l/2$} \, \quad ( 2 ) & \text{We d} \\ \rho _ { e } & \text{if $r\geq R_{0}+l/2$} & \text{as} \ [ 9 ] } \end{cases}$$

 where R 0 represents the radius of the loop. For the density ρ t ( r ) of the transitional layer we consider linear ( li ), parabolic ( p ), and sinusoidal ( s ) profiles [34]:

$$\rho _ { l i } ( r ) & = ( \rho _ { e } - \rho _ { i } ) ( r - R _ { 0 } ) / l + ( \rho _ { e } + \rho _ { i } ) / 2, \quad \ \ ( 3 ) \quad \text{where} \\ \lambda \, \widehat { \nu } \, = \, \lambda \, - \, \widehat { \nu } \, - \, \widehat { \nu } \, - \, \widehat { \nu } \, - \, \widehat { \nu } \, + \, 1 / \mathfrak { M } ^ { 2 } / 1 ^ { 2 } \quad \ \ \ ( 4 \, \mathfrak { M } \, \quad \ \ \mathfrak { T }.$$

$$\rho _ { s } ( r ) & = [ ( \rho _ { e } + \rho _ { i } ) + ( \rho _ { e } - \rho _ { i } ) \sin ( \pi ( r - R _ { 0 } ) / l ) ] / 2, \ \ ( 5 ) \quad \text{versu} \\ \mathfrak { 1 } \cdot \mathfrak { 1 } \quad \mathfrak { 1 } \quad \mathfrak { 1 } \quad \mathfrak { 1 }$$

$$\mu _ { p } ( r ) & = \mu _ { i } - ( \rho _ { i } - \rho _ { e } ) ( r - R _ { 0 } + l / 2 ) ^ { 2 } / l ^ { 2 }, \quad \ ( 4 ) \quad \ \text{In} \\ \alpha _ { l } ( r ) & = \left ( l \alpha _ { l } + \alpha _ { l } \right \rangle + \alpha _ { l } - \alpha _ { l } \right \rangle \sin ( \pi ( r - R _ { \cap } \lambda / l ) ) / 9 \quad l 5 \quad v e s i n$$

which are shown in Fig. 1.

For mode conversion and the related resonant absorption to occur, the wave frequency should reside in the Aflv´n continuum ( e ω Ae &lt; ω &lt; ω Ai ) of the loop. We consider that an external wave is incident on a coronal loop with the frequency of the fundamental standing kink mode, ω k = k B z 0 √ 2 /µ 0 ( ρ i + ρ e ) with k z = π/L 0 , where L 0 is the loop length. Thus mode conversion occurs where ω A = ω k or, in other words, ρ r ( ) = ( ρ e + ρ i ) / 2. We consider the frequency dependence of mode conversion afterward. The wave outside the loop can be described with an incident wave and a scattered wave, which is

$$P ( r, \phi ) = & \sum _ { m } a _ { m } e ^ { i m \phi } \{ J _ { m } [ k ( r - R ) + c ] \quad \begin{array} { c } \text{in} \\ & + H _ { m } ^ { ( 1 ) } [ k ( r - R ) + c ] r _ { m } ( R ) \}, \quad r > R, \quad ( 6 ) \end{array} \text{ \quad \begin{array} { c } \text{in} \\ \text{effic} \\ & \text{$l/R_{0}$} \\ \text{where $J_{m}(H_{m}^{(1)} )$ is the Bessel $(Hankel)$ function of the} \quad \text{$the$ a fin} \\ \text{$i_m_=1$} \end{array}$$

where J m ( H (1) m ) is the Bessel (Hankel) function of the first kind, r m the scattering coefficient, c a constant equivalent to kR , and k (= √ ( ω /v 2 2 A ) -k 2 z ) the radial wave number for r &gt; R . A constant c is used both to apply the invariant imbedding method and for Eq. (6) to be eigenfunctions for r &gt; R [44]. The value a m depends on the form of the incident wave. For instance, when a plane wave is incident to the x direction, a m = i m [36]. Applying IIM to Eq. (1) with Eq. (6) and using the thickness R as a new variable (imbedding parameter), we obtain for r m [44]

$$\frac { \infty \ a \text{ new values } \max \max \longrightarrow \gamma _ { m }, \, \text{ we } \text{voxam} \quad \overline { \overline { d r ^ { 2 } } } } { \overline { d r _ { m } } } = \frac { k } { H _ { m } } \frac { \mu ( R ) } { \mu _ { 1 } } [ J _ { m } ^ { \prime } + H _ { m } ^ { ( 1 ) ^ { \prime } } r _ { m } ( R ) ] & ( 7 ) \\ + k \left [ \frac { \mu ( R ) } { \mu _ { 1 } } \frac { H _ { m } ^ { ( 1 ) ^ { \prime } } } { H _ { m } ^ { ( 1 ) } } + \frac { 1 } { k R } \right ] & \text{where} \quad \kappa = ( \cdot \\ \text{ only } \\ \times \frac { [ J _ { m } ^ { \prime } + H _ { m } ^ { ( 1 ) ^ { \prime } } r _ { m } ( R ) ] [ J _ { m } + H _ { m } ^ { ( 1 ) } r _ { m } ( R ) ] } { J _ { m } ^ { \prime } H _ { m } ^ { ( 1 ) } - H _ { m } ^ { ( 1 ) ^ { \prime } } J _ { m } } & \text{ Thus} \\ \text{ files.} \\ + k \left [ \frac { \epsilon ( R ) } { \epsilon _ { 1 } } - \frac { \mu _ { 1 } } { \mu ( R ) } \frac { m ^ { 2 } } { k ^ { 2 } R ^ { 2 } } \right ] \frac { [ J _ { m } + H _ { m } ^ { ( m ) } r _ { m } ( R ) ] ^ { 2 } } { J _ { m } ^ { \prime } H _ { m } ^ { ( 1 ) } - H _ { m } ^ { ( 1 ) ^ { \prime } } J _ { m } }, & \text{ profil.} \\ 1 \right ] ^ { 1 / 3 }$$

where µ r ( ) = ρ 0 ( r )( ω 2 -ω 2 A ( r )), /epsilon1 ( r ) = µ /k B 0 2 0 2 0 , k 0 = ω/c 0 , and c 0 is the speed of the electromagnetic wave in vacuum, and µ 1 and /epsilon1 1 are the values of µ , /epsilon1 for r &gt; R , prime denotes df ( x /dx ) for f ( x ), J m = J m ( c ), and H (1) m = H (1) m ( c ). Eq. (7) is a first order ordinary differential equation with R as an integration variable. We integrate Eq. (7) from 0 to R using initial conditions r m (0) = 0 to obtain r m ( r = R ) [44], where the 1 /kR term gives rise to a singularity at R = 0 which we avoid by letting the values of parameters for R &lt; δ equal those for R = δ /lessmuch 1. This singularity is due to the cylindrical geometry we consider here. The value of δ needs to be sufficiently small to not affect the results. We also set ω = ω k + iγ to avoid the singularity at the resonance position ( µ = 0) where γ /lessmuch 1, not affecting the results. We define the mode conversion (absorption) coefficient A as [9, 44]

$$A = - \sum _ { m = \pm 1 } \text{Re} ( r _ { m } + \left | r _ { m } \right | ^ { 2 } ), \text{ \quad \ \ } ( 8 )$$

where two kink modes ( m = ± 1) are considered.

In Fig. 2, we show the mode conversion coefficient A versus l/R 0 for each density profile when ∆ ρ (= ρ /ρ e i ) = 2 and 10, and L /R 0 0 = 50, 100, and 150, while other fixed parameters are ω = ω k , k z = π/L 0 , R = 3 R 0 , ρ i = 1 67353 . × 10 -12 kg/m 3 , B 0 = 10 -3 T , R 0 = 2 × 10 6 m , γ = 10 -8 s -1 , and δ = 10 -6 . We find that the mode conversion is strong when l/R 0 is small, depending on the loop length and density contrast. The peak positions shift towards a lower value of l/R 0 as ∆ ρ decreases and L /R 0 0 increases. The figure shows that mode conversion for the parabolic profile is less effective than the other two density profiles, which is due to the difference of the resonance point.

In Fig. 3 (a)-(c), we plot the mode conversion coefficient A versus L /R 0 0 for each density profile when l/R 0 = 0.002, 0.2 and 2.0, and ∆ ρ = 2, 5, and 20. For each l/R 0 , the peak value of A for each density profile has the same value regardless of ∆ ρ , which draws us to a finding below. The mode conversion coefficient can be plotted as a single curve when expressed with a scaling factor for some simple density (or Alfv´n speed) proe files [8, 12, 13, 19, 23]. Approximating the density profile near the resonance (mode conversion) position r = R ∗ as a linear function ρ t ( r ) = a r ( -R ∗ ) + b ( a and b are constants) and including ω = ω k + iγ , the wave equation reduces to [44]

$$\overset { \mu } { \text{max} } \, \underset { \text{$f$} \, \} } { \text{tain} } \, \quad \frac { d ^ { 2 } P } { d \bar { r } ^ { 2 } } + \left ( \frac { 1 } { \bar { r } } - \frac { 1 } { \bar { r } - q + i \bar { \gamma } } \right ) \frac { d P } { d \bar { r } } + \left ( \bar { r } - q + i \bar { \gamma } - \frac { m ^ { 2 } } { \bar { r } ^ { 2 } } \right ) P = 0,$$

where ¯ = r κr , q = κR ∗ , ¯ = 2( γ a/b ) -2 / 3 k z 2 / 3 ( η/ω k ) and κ = ( a/b ) 1 / 3 ( k z ) 2 / 3 . In the limit ¯ ( γ γ ) → 0, this equation only depends on the value of q for a given m [19, 21, 44]. Thus one can expect a single curve of the mode conversion coefficient in terms of q for simple density profiles. For the linear, sinusoidal, and parabolic density profiles, we obtain that q equals [( πR /l 0 )(∆ ρ -1) / (∆ ρ + 1)] 1 / 3 ( k R z 0 ) 2 / 3 , [(2 R /l 0 )(∆ ρ -1) / (∆ ρ +1)] 1 / 3 ( k R z 0 ) 2 / 3 ,

and [(2 √ 2 R /l 0 )(∆ ρ -1) / (∆ ρ +1)] 1 / 3 ( k R z 0 ) 2 / 3 , respectively [cf., 44]. With those scaling factors, the mode conversion coefficient A in Fig. 3 (a)-(c) reduces to a single curve in Fig. 3 (d)-(f) for each density profile. A difference from the results in planar geometry is that A has some dependence on l/R 0 despite q is a function of l/R 0 . However, Fig. 3 points out that the scaling behavior is a general property if the density profile can approximately be linearized at the resonance point, reinforcing and extending the results from planar geometry. It is straightforward to apply these results to the mode conversion of electromagnetic waves into electrostatic modes by using the mathematical analogy of the two governing wave equations [19, 23, 44]. Comparing the scaling behavior in planar and cylindrical geometries, we propose that a scaling behavior can exist regardless of the geometry for density configurations that are approximately linear around the resonant point.

In contrast with previous results, we here consider the frequency dependence of the mode conversion. Fig. 4 presents the dependence of the mode conversion on the wave frequency. We draw A versus ω/ω k for ∆ =5, ρ l/R 0 =0.002, 0.2 and 2.0, and L /R 0 0 =20 and 100. Mode conversion is strongest when ω ≈ ω k and l/R 0 is small, but the peak position shifts downwards as l/R 0 increases. Comparing the top and bottom rows, the shape of the mode conversion coefficient is almost independent of L /R 0 0 .

In this paper, we have generalised our earlier results on universal scaling in resonant absorption to systems with different density profiles. We have done this using the invariant embedding approach. We found, irrespective of the particular density profile, that mode conversion (resonant absorption) is efficient when the thickness of the transitional layers ( l/R 0 ) is small and the loop length ( L /R 0 0 ) is relatively large. From this dependence irrespective of the density profile (as long as it is approximately linear near the resonant point), we thus infer that a universal scaling behavior of mode conversion is found for coronal loops exhibiting kink oscillations. In this particular case, the scaling formulas is useful, because it could be possible to infer the wave energy absorption

- [1] Arregui, I., Advances in Space Research 61 , 655 (2018).
- [2] Babkin, G. I. and Klyatskin, V. I., Soviet Journal of Experimental and Theoretical Physics 52 , 416 (1980).
- [3] Brodin, G. and Marklund, M., Phys. Rev. Lett. 82 , 3012 (1999).
- [4] Chen, L. and Hasegawa, A., Physics of Fluids 17 , 1399 (1974).
- [5] Davila, J. M., Astrophys. J. 317 , 514 (1987).
- [6] Ding, Y. S., Chan, C. T., and Wang, R. P., Scientific Reports 3 , 2954 (2013).
- [7] Doucot, B. and Rammal, R., Journal de Physique 48 , 527 (1987).
- [8] Forslund, D. W., Kindel, J. M., Lee, K., Lindman, E. L., and Morse, R. L., Phys. Rev. A 11 , 679 (1975).

due to resonant absorption from the observed wave parameters [1]. We expect that a similar scaling behavior is expected for resonant absorption of usual coronal loop oscillations for which the density contrast ∆ ρ &lt; 1.

The results of our study also have implications beyond the field of solar physics. Our results show that the universal scaling behaviour for mode conversion in resonant absorption extends to cylindrical geometry, adding to the previous knowledge on planar geometry. The universal scaling is thus truely universal, and happens irrespective of geometry. It is generally applicable to other forms of mode conversion as well. In particular, we can infer that the universal scaling also applies to mode conversion of electromagnetic waves into electrostatics waves in unmagnetized plasmas or in metamaterials [6, 21, 37] in cylindrical (or spherical) geometry, because that is also described with the same mathematical analogy of the wave equations [44] as here. Furthermore, going beyond mode conversion of electromagnetic or magnetohydrodynamic waves, other kind of mode conversion like e.g. between gravitational waves and electromagnetic waves are also possible [3, 26]. Our results allow us to postulate that, even for these exotic cases, a scaling behavior also exists. We can thus conclude the paper by formulating the hypothesis: 'a scaling behavior of the mode conversion coefficient is possible for monotonic density profiles regardless of the geometry and the type of mode conversion.'

## ACKNOWLEDGMENTS

This work was based on the presentation in 43rd European Physical Society (EPS) Conference on Plasma Physics (2016). D.J.Y. thanks the support from the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea. T.V.D. thanks the support from GOA-2015-014 (KU Leuven), and the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation programme (grant agreement No. 724326).

- [9] Hanson, C. S. and Cally, P. S., Astrophys. J. 791 , 129 (2014).
- [10] Hinkel-Lipsker, D. E., Fried, B. D., and Morales, G. J., Phys. Rev. Lett. 62 , 2680 (1989).
- [11] Hinkel-Lipsker, D. E., Fried, B. D., and Morales, G. J., Phys. Rev. Lett. 66 , 1862 (1991).
- [12] Hinkel-Lipsker, D. E., Fried, B. D., and Morales, G. J., Physics of Fluids B 4 , 559 (1992).
- [13] Hinkel-Lipsker, D. E., Fried, B. D., and Morales, G. J., Physics of Fluids B 5 , 1746 (1993).
- [14] Ionson, J. A., Astrophys. J. 226 , 650 (1978).
- [15] Johnson, J. R., Cheng, C. Z., and Song, P., Geophysical Research Letters 28 , 227 (2001).

- [16] Kim, E.-H., Cairns, I. H., and Robinson, P. A., Phys. Rev. Lett. 99 , 015003 (2007).
- [17] Kim, E.-H., Cairns, I. H., and Robinson, P. A., Physics of Plasmas 15 , 102110 (2008).
- [18] Kim, K., Phys. Rev. B 58 , 6153 (1998).
- [19] Kim, K. and Lee, D.-H., Physics of Plasmas 12 , 062101 (2005).
- [20] Kim, K. and Lee, D.-H., Physics of Plasmas 13 , 042103 (2006).
- [21] Kim, K., Lee, D. H., and Lim, H., Optics Express 16 , 18505 (2008).
- [22] Kim, K., Lim, H., and Lee, D.-H., Journal of Korean Physical Society 39 , L956 (2001).
- [23] Kivelson, M. G. and Southwood, D. J., Journal of Geophysical Research 91 , 4345 (1986).
- [24] Klyatskin, V. I., Progress in optics 33 , 1 (1994).
- [25] Klyatskin, V. I., Stochastic equations through the eye of the physicist: Basic concepts, exact results and asymptotic approximations (Elsevier, Amsterdam, 2005).
- [26] Marklund, M., Brodin, G., and Dunsby, P. K. S., Astrophys. J. 536 , 875 (2000).
- [27] Mjølhus, E., Radio Science 25 , 1321 (1990).
- [28] Nakariakov, V. M., Ofman, L., Deluca, E. E., Roberts, B., and Davila, J. M., Science 285 , 862 (1999).
- [29] Ofman, L., Davila, J. M., and Steinolfson, R. S., Astrophys. J. 444 , 471 (1995).
- [30] Okamoto, T. J., Antolin, P., De Pontieu, B., Uitenbroek, H., Van Doorsselaere, T., and Yokoyama, T., Astrophys. J. 809 , 71 (2015).
- [31] Poedts, S., Goossens, M., and Kerner, W., Astrophys. J. 360 , 279 (1990).
- [32] Rammal, R. and Doucot, B., Journal de Physique 48 , 509 (1987).
- [33]
- Schleyer, F., Cairns, I. H., and Kim, E.-H., Journal of Geophysical Research (Space Physics) 119 , 3392 (2014).
- [34] Soler, R., Goossens, M., Terradas, J., and Oliver, R., Astrophys. J. 777 , 158 (2013).
- [35] Steinolfson, R. S. and Davila, J. M., Astrophys. J. 415 , 354 (1993).
- [36] Stratton, J. A., Electromagnetic theory (John Wiley &amp; Sons, 2007).
- [37] Sun, J., Liu, X., Zhou, J., Kudyshev, Z., and Litchinitser, N. M., Scientific Reports 5 , 16154 (2015).
- [38] Willes, A. J. and Cairns, I. H., Publ. Astron. Soc. Aust. 18 , 355 (2001).
- [39] Willes, A. J. and Cairns, I. H., Physics of Plasmas 10 , 4072 (2003).
- [40] Yu, D. J. and Kim, K., Physics of Plasmas 20 , 122104 (2013).
- [41] Yu, D. J. and Kim, K., Physics of Plasmas 23 , 032112 (2016).
- [42] Yu, D. J., Kim, K., and Lee, D.-H., Physics of Plasmas 17 , 102110 (2010).
- [43] Yu, D. J., Kim, K., and Lee, D.-H., Physics of Plasmas 20 , 062109 (2013).
- [44] Yu, D. J. and Van Doorsselaere, T., Astrophys. J. 831 , 30 (2016).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00214v1-figure-06900a4f-43ba-4e89-b834-ff83fcf09380.jpg)

FIG. 1. Density profiles in the transitional layer when l/R 0 = 0 2 and ∆ . ρ (= ρ /ρ e i ) = 5: linear (red), sinusoidal (blue), and parabolic (green) profiles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00214v1-figure-b24bfe20-8b9b-4973-9b35-0c81eb8b5e06.jpg)

/s32

/s32

FIG. 2. Mode conversion coefficient A for linear (black), parabolic (red), and sinusoidal (blue) density profiles as a function of l/R 0 for ∆ ρ (= ρ /ρ e i ) = 2 (top), 10 (bottom) and for L/R 0 = 50 (left), 100 (center), and 150 (right), where ω = ω k , k z = π/L , R = 3 R 0 , ρ i = 1 67353 . × 10 -12 kg/m 3 , B 0 = 10 -3 T , R 0 = 2 × 10 6 m γ , = 10 -8 s -1 , and δ = 10 -6 . Mode conversion is strong when l/R 0 is small.

/s32

/s32

FIG. 3. Mode conversion coefficient A for linear (black), parabolic (red), and sinusoidal (blue) density profiles as a function of L /R 0 0 (a)-(c) and of q (d)-(f) for l/R 0 = 0.002 (left), 0.2 (center) and 2.0 (right), where ∆ ρ = 2, 5, and 20. In (a)-(c), as ∆ ρ increases the peak position of the curves shifts toward the higher values of L /R 0 0 . As a function of q , A reduces into a single curve for each density profiles, while the curves still depend on the value of l/R 0 . The scaling parameter q is defined as [ c 1 ( R /l 0 )(∆ ρ -1) / (∆ ρ +1)] 1 / 3 ( k R z 0 ) 2 / 3 where c 1 corresponds to π , 2 √ 2, and 2 for linear, parabolic, and sinusoidal profiles, respectively. The other parameters are the same as in Fig. 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00214v1-figure-34afc6da-c1e1-475c-be47-489656f31e87.jpg)

FIG. 4. Mode conversion coefficient A for linear (black), parabolic (red), and sinusoidal (blue) density profiles as a function of ω/ω k for ∆ ρ =5, l/R 0 =0.002 (left), 0.2 (center) and 1.0 (right), and L /R 0 0 =20 (top) and 100 (bottom). Mode conversion is strong when ω ≈ ω k and l/R 0 is small. The relative relations of A among the density profiles do not depend on the value of L /R 0 0 . The other parameters are the same as in Fig. 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00214v1-figure-e78fb226-de9c-41de-8261-821a1b90d271.jpg)
Dae Jung Yu Tom Van Doorsselaere
Universal scaling behavior of resonant absorption
Mode conversion and resonant absorption are crucial mechanisms for wave
transport and absorption. Scaling behavior of mode conversion or resonant
absorption is well-known for electromagnetic and MHD waves in planar geometry.
Our recent study showed that such a scaling behavior of resonant absorption
could also exist for coronal loop oscillations with cylindrical geometry, but
it was only tested for one density profile. Here we generalise our previous
study on the scaling behavior of resonant absorption by considering multiple
density profiles. Applying an invariant imbedding method to the ideal MHD wave
equations, we show that the scaling behavior also exists for these density
models. We thus generalise our earlier results and show that such a universal
scaling exists in cylindrical geometry, too. Given these results and the
earlier results in planar geometry, we formulate a hypothesis that a universal
scaling behavior exists regardless of the type of mode conversion or resonant
absorption.