Magnetoacoustic and Alfvénic black holes

We introduce analogue black holes (BHs) based on ideal magnetohydrodynamic equations. Similar to acoustic BHs, which trap phonons and emit Hawking radiation (HR) at the sonic horizon where the flow speed changes from super- to sub-sonic, in the horizon of magnetoacoustic and Alfvénic BHs, the magnetoacoustic and Alfvén waves will be trapped and emit HR made of quantized vibrations similar to phonons which we call magnephonons and Alphonons. We proposed that magnetoacoustic and Alfvénic BHs may be created in the laboratory using a tube with variable cross section embedded in a uniform magnetic field, and a super-magnetoacoustic or a super-Alfvénic flow. We show that the Hawking temperature for both BHs is a function of the background magnetic field, number density of fluid, and radius of the tube. For a typical setup, the temperature is estimated to be about 0.0266 K.


Introduction
In 1916 Schwarzschild gave a metric as a solution of the Einstein field equation. Singularity of such a metric predicted a gravitational BH with event horizon at Schwarzschild radius [1]. Although, based on classical physics everything, even light, is absorbed by BHs and cannot escape, in the context of quantum field theory in curved space, Hawking showed that BHs should emit black body radiation [2,3]. HR in the universe has not been observed yet, but numerous attempts have been done to simulate the interesting phenomena in the laboratory. Unruh showed that HR is not only a characteristic of gravitational BHs, it is also a characteristic of the acoustic analogue BH [4,5]. After 1981, most of attempts are proposed based on Bose-Einstein condensates of quantum fluid, [6][7][8][9][10] quasi particles in superfluid [11], ultra-cold fermions [12], in plasmas and ion rings, [13][14][15][16][17] slow light in an atomic vapor, [18][19][20][21] in water [22,23], etc. Recently, a e-mail: safari@znu.ac.ir observation of self-amplifying HR in an analogue BH laser suggested a very promising experiment method for probing the inside of a BH [24]. From a theoretical point of view, the acoustic analogue BH models are developed in geometrical acoustics and physical acoustics [25]. Using the linearized hydrodynamic equations in the presence of initial material flow, a wave equation for velocity potential was obtained. Tensorial form of the wave equation results in an acoustic metric. The acoustic metric is singular at a point where the local sound speed is equal to the flow speed [26]. This was interpreted as characteristic of a sonic BH [25,26].The effect of magnetic field on the acoustic BH and HR has not been studied, yet.
Alternatively, the idea for definition of acoustic BH can be applied to introduce new analogue magnetoacoustic and Alfvénic BHs in the magnetohydrodynamics (MHD) framework. Magnetohydrodynamics (MHD) is an useful approach to analyze characteristics (flow, wave and dissipation etc) of the laboratory and astrophysical plasma [33]. After Alfvén [34], MHD waves (Alfvén and magnetoacoustic) have been detected, using the laboratory experiments Alfvén wave (e.g. [35][36][37], in the earth-ionosphere and magnetic field [38,39], a variety normal modes of the solar corona [40][41][42][43][44]. Here, first we treat the magnetoacoustic, Alfvén, and acoustic waves based on Helmholtz theorem for a uniform and stationary medium with constant background magnetic field. Second, mimicking the definition of acoustic BH in a nozzle, we introduce magnetoacoustic and Alfvénic BHs with using a slightly variable cross section tube. We conclude that at the horizon of magnetoacoustic and Alfvénic BHs should emit radiations made of the Magnephonon and Alphonon, respectively. We define two quasi-particles Magnephonon and Alphonon correspond to quantum of the magnetoacoustic and Alfvén waves, respectively. This paper is organized as follows: Sect. 2 gives the properties of magnetoacoustic waves (fast, slow, and Alfvén waves) using the Helmholtz decomposition and explains a derivation of magnetoacoustic metric in the basis of linearized ideal magnetohydrodynamic equations. Section 4 introduces the acoustic, magnetoacoustic, and Alfvénic black holes, respectively. Section 5 calculate the Hawking temperature for the acoustic, magnetoacoustic, and Alfvénic black holes. Section 6 describes the conclusions.

Magnetoacoustic waves and metric
Here, we give the conditions for a definition of a magnetoacoustic metric, in the non-relativistic magnetohydrodynamic (MHD) framework. In the MHD approach, the behavior of continuous plasma is governed by a non-relativistic form of Maxwell's equations, together with Ohm's law, a gas law, equation of mass continuity, motion and energy equations [28]. The ideal MHD equations for an adiabatic process and irrotational flow ∇ × v = 0 are given by: where, ρ, p, v, B, µ, γ , and k are density, pressure, flow velocity, magnetic field, magnetic permeability, atomicity coefficient, and a constant, respectively. For the derivation of a metric for magnetoacoustic wave we need to linearised the MHD equations with choosing the velocity disturbances as a velocity potential in similar analytical process for derivation of acoustic metric from wave equation. To do this, first we treat the propagation of magnetoacoustic waves in homogenous unbounded medium choosing velocity disturbance as a velocity potential and a vector potential. Second, the magnetoacoustic and Alfvénic metrics are calculated.

Helmholtz theorem and magnetoacoustic waves
Properties of magnetoacoustic waves (fast, slow, and Alfvén waves) in an unbounded homogenous and stationary medium with constant density (ρ 0 = cte), constant pressure ( p 0 = cte), and uniform background magnetic field (B 0 ), were investigated in the literature [28]. The linearised ideal MHD equations can be reduced to a single equation for disturbed velocity as [28], where c 0 = γ p 0 ρ 0 is the sound speed, v A = B 0 √ μρ 0 is speed of the Alfvén wave, andB 0 = B 0 /B 0 is a unit vector. Here, we focus on studying of the characteristics of the above mentioned waves using a fundamental theorem of calculus (Helmholtz's theorem). Based on the Helmholtz decomposition, a vector field with sufficient smoothness and decay conditions [29], can be decomposed to an irrotational part (∇φ where φ is a scaler) and a solenoidal part (∇ × A where A is the vector potential and satisfy ∇.A = 0). The irrotational (gradient) and divergence-free solutions can be used for treating the longitudinal and transversal waves, respectively. Suppose where k, ω,φ, andv = ikφ, are wave vector, angular frequency of oscillations, a constant, and wave amplitude, respectively. Furthermore Eq. (6) yields, Equation (7) simplifies as We note that the wave vector (k) and wave amplitude (v) are parallel (k v). Considering the wave propagation parallel to the background magnetic field (k B 0 ), Eq. (8) reduces to the dispersion relation ω 2 = c 2 0 k 2 . This indicates that, in the case of an irrotational solution parallel to the magnetic field, only the acoustic wave can propagate. In the case of propagation perpendicular to the background magnetic field (k⊥B 0 ), Eq. (8) gives the dispersion relation ω 2 = ( v 2 A + c 2 0 )k 2 . This is the well-known characteristic of the fast magnetoacoustic wave with the phase speed c 2 For the oblique propagation k.B 0 = k B 0 cos θ (θ is the angle between B 0 and k), the phase speed of the slow magnetoacoustic wave is given by From the above analysis we see that, choosing the disturbed velocity as an irrotational velocity field v = ∇ϕ, the Alfvén waves cannot propagate along the background magnetic field We see that the velocity amplitude is perpendicular to the wave vector (k.v = 0), and only the transversal Alfvén wave with phase speed v ph = v A can propagate.
In the remainder of this section, the metrics for acoustic, magnetoacoustic, and Alfvén longitudinal waves (irrotational solutions) are derived.

Magnetoacoustic metric
Here, in the presence of initial material flow the magnetoacoustic metric using the magnetoacoustic wave is derived. For an irrotational flow the velocity is satisfied by a scalar field v = ∇ϕ. We consider small perturbations from equilibrium as where, equilibrium quantities indicated by subscript "0" are function of position x and time t, ψ(r, t), ϕ(r, t), and B 1 (r, t) are perturbed quantities [30,31]. Equilibrium quantities (ρ 0 , v 0 , B 0 , and p 0 ) are satisfied by Linearization of Eqs. (1)-(5) (products and squares of the small perturbations are neglected) and after some mathematical manipulations, give As we explained in the previous section, by choosing v = ∇ϕ the parallel propagation of Alfvén waves along the magnetic field (B 0 = B 0x ) is absent. Therefore, the first term on right hand of Eq. (21) is set to zero, ∇ × (v 0 × B 1 ) = 0. A combination of this assumption and the solenoidal condition for magnetic field (∇.B 0 = 0), gives ∂ B 1x ∂ x = 0. Thus the propagation of Alfvén waves along the background magnetic field is removed from our analysis. We assume the irrotational part of the vector v 0 × B 1 = ∇δ(x, y, z, t) in which δ(x, y, z, t) is a function. The x and z components of Eq. (20) are Here after, first, we focus our analysis on the derivation of acoustic waves with propagation along the background magnetic field directions (x) and second, magnetoacoustic wave propagation in all directions except the magnetic field direction. First, combining Eqs. (19) and (23) gives Equation (25) is the well-known Klein-Gordon equation for acoustic waves. Eliminating ∂ψ/∂x between Eqs. (19) and Briefly, by choosing the irrotational solution for the velocity disturbance, the transversal Alfvén wave is absent, and the propagation of the longitudinal Alfvén wave parallel with the background magnetic field is also absent as expected. Our analysis shows that, along the magnetic field only the acoustic wave can propagate. Because our goal is to analyse the magnetoacoustic black hole, hereafter, we focus our analysis in all directions except the magnetic field direction. Second, by differentiating Eq. (26) with respect to z and substituting ∂ψ/∂z from Eq. (24) one finds in which, ∂z from Eq. (21), into Eq. (27) and differentiating with respect to t using ∂ ∂t (∇ × B 1 ) y = −B 0 ∇ 2 ∂ϕ ∂z we derived a single equation for velocity potential (ϕ) where, For a flow having a slight change in the speed (dv 0 /dx 1), for high frequency waves (short period t 1) the term t (dv 0 /dx) becomes too small. In this case, the right hand side term B 0 ∇ 2 ϕ of Eq. (28) Equation (31) describes the propagation of acoustic, Alfvén, and magnetoacoustic waves in laboratory and astrophysical plasma. This equation is in the form of the well-known Klein-Gordon equation. As expected, in the case of unmagnetized fluid (B 0 = 0), Eq. (31) reduces to the acoustic wave equation for velocity potential. Usually, a d'Alembertian equation (for a minimally coupled massless scalar field) of motion was derived for velocity potential in a barotropic, inviscid, and rotational free flow [25]. Equation (31) can be reformulate in a tensorial form where, g μν and g are inverse metric tensor and its determinant, μ and ν runs from 0 (indicates the time coordinate) to 3 (denotes the spatial coordinate). The magnetoacoustic inverse metric tensor g μν and metric tensor g μν are obtained as where, c 2 f = v 2 A + c 2 0 . Using metric tensor, Eq. (33), the magnetoacoustic interval can be defined as Inserting the specific time interval dτ = dt (34), one obtains In the following section, the properties of acoustic, magnetoacoustic, and Alfvénic BHs are investigated.

Acoustic black hole
The theory of gravitational BHs has been developed into the supersonic flow by Unruh [26]. For a moving fluid medium at the horizon where speed of medium is closed to propagation speed of the acoustic signals " then nothing can fight its way back upstream and signals are trapped" [27]. In the case of unmagnetised gas(B 0 = 0), Eq. (28) reduces to the Klein-Gordon equation, for propagation of acoustic waves in the presence of material flow. Using the resultant equation the acoustic metric can be derived [26]. The acoustic metric can be obtained by setting v A = 0 in Eq. (33) Equivalently, the acoustic interval can be expressed as Combination of continuity and momentum equations (Eqs. 14 and 15) in stationary state, the relation between cross section S and velocity v is given by This relation shows that for (d S < 0) a subsonic flow (v 0 < c 0 ) will be accelerated and a supersonic flow (v 0 > c 0 ) will be decelerated. If the nozzle is sufficiently narrow and with a slightly variable cross section the speed of flow exceeds to sound speed at the throat (sonic horizon). This shows the acoustic interval Eq. (37) interpreted acoustic BH which has a sonic horizon and trapped phonon in the surface gravity of acoustic BH. This means that, when the acoustic waves cross from upstream to downstream, the acoustic wave quanta (phonons) are captured in the horizon of the BH where they emit Hawking radiation made by phonons. In this regard, in many papers a Laval nozzle setup ( Fig. 1) has been proposed to discuss the above mentioned acoustic BH. This setup uses, an axisymmetric sufficiently thin tube with slightly decreasing cross section (S(x)) that reaches its minimum cross section at the throat and then slightly increases that. An initial material flow (v 0 = v 0 (x)x) along the tube axis is considered.

Magnetoacoustic black hole
The magnetoacoustic metric Eq. (35) is singular at the magnetoacoustic point, where v 0 = c f , determines a magnetoacoustic horizon. The speed of super magnetoacoustic plasma flow reduces to local propagation speed of magnetoacoustic wave at horizon; then signal of magnetoacoustic wave is trapped and therefore it can be called magnetoacoustic BH. Similar to the HR emitted from acoustic and gravitational BHs, the magnetoacoustic BH also should emit HR.
In this regard, we propose a setup (Fig. 2) to discuss the above mentioned BH. The setup consists of an axisymmetric sufficiently thin tube with slightly variable cross section (S(x)), a uniform force free magnetic field B 0 = B 0x and an initial material flow (v 0 = v 0 (x)x) along the tube axis (Fig. 2). A similar treatment of sub and supersonic flow in tube configuration can be explained for sub and supermagnetoacoustic flow based on Eq. (38). In other words, the super-magnetoacoustic flow (v 0 > c f ) will be decelerated along the tube where its cross section slightly decreases. It is possible to release a super-magnetoacoustic flow in the tube, which its speed tends to the speed of magnetoacoustic wave (v 0 = c f ) at the horizon. The boundary between sub-magnetoacoustic and super-magnetoacoustic flow could be called the magnetoacoustic horizon, analogous to the sonic horizon in acoustic BHs. Phononic version of HR is an inevitable result of trapping acoustic wave at the acoustic horizon [24,32]. Under a likely scenario, the magnetoacoustic wave cannot escape from the magnetoacoustic horizon, therefore should emit HR made of magnephonon. Indeed, a magnephonon will be a quantum for magnetoacoustic wave, analogous to the phonon which is a quantum for acoustic wave.

Alfvénic black hole
In the limit of c 0 << v A (zero β plasma condition), the magnetoacoustic wave (Eq. 31) reduces to an Alfvén wave in the presence of initial material flow. The resultant Alfvén wave equation is in the form of Klein-Gordon equation. Immediately, the Alfvénic metric can be derived from magnetoacoustic metric (Eqs. 33 and 35) by setting c 0 tends to 0, Equivalent Alfvénic interval is given by The Alfvénic metric (Eq. 40) is singular at the location of Alfvénic point where v 0 = v A . This singular behaviour of Alfvénic metric leads to an Alfvénic BH. We propose a setups that is illustrated in Fig. 3, the condition for occurrence of an Alfvénic BH can be discussed. For compressional (longitudinal) Alfvén wave, the plasma density is nearly constant. As a result of mass continuity, the flux, S(x)v 0 (x), is constant in the tube cross section. Consider a super-Alfvénic flow in a tube with increasing cross section (Fig. 3), speed of the flow decreases along the tube axis and reaches to Alfvén velocity v 0 = v A at Alfvénic point and then transformed to sub-Alfvénic flow. Therefore, we call Alfvénic horizon to be the interface between super-Alfvénic flow and sub-Alfvénic flow. In horizon of the Alfvénic BH, the Alfvén wave is trapped and it is expected to be radiated by Alphonon. Alphonon is introduced as a quantum particle for Alfvén wave packet.

Hawking temperature
Hawking temperature is an important characteristic of BHs. Unruh showed that fluid flows mimic BHs. Hawking temperature (T H ) for acoustic BH was obtained where,h = h/2π , h is the plank constant and k is the Boltzmann constant. Since the Hawking temperature is independent of metric conformal factor. It will therefore be as following for the magnetoacoustic BH, In the limit of v A tends to zero, Eq. (42) then reduces to Hawking temperature for acoustic BH. Although, in the zero β plasma condition the above mentioned equation can describe the Hawking temperature for Alfvénic BH, At the horizon where the tube cross section radius is equal to R, Eq. (42) can be simplified as where, the term 1

Conclusion
In this study, we introduced the magnetoacoustic and Alfvénic analogue BHs. In the horizon of magnetoacoustic and Alfvénic BHs, the magnetoacoustic and Alfvén waves are trapped, respectively, and should emit magnephonons and Alphonons version of HR at Hawking temperature. The next logical step is to investigate the physical properties of both magnephonon and Alphonon quasi-particles based on quantum approach. As stated in the literature, for acoustic BH in a natural fluid, the Hawking temperature is a function of sound speed and geometry of nozzle setup at the horizon. However, Hawking temperature for magnetoacoustic BH is related to the sound speed with additional positive terms that depends on the magnetic field and density. Magnephonons and Alphonons are particles (quanta) corresponding to magnetoacoustic and Alfvén waves, respectively. The idea for definition of acoustic, magnetoacoustic, and Alfvénic BHs can be applied to theoretical and/or experimental prediction of new non-gravitational BHs. Perhaps, the study of new BHs could help us observe HR.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .