Magnetoacoustic and Alfve'nic Black Holes with Hawking Radiation at Horizons Made of Magnephonons and Alphonons

TL;DR

This paper proposes a new type of analogue black holes based on magnetohydrodynamics, which trap magnetoacoustic and Alfvén waves, emitting Hawking radiation composed of quantized vibrations called magnephonons and Alphonons, with potential laboratory creation.

Contribution

It introduces magnetoacoustic and Alfvénic black holes as novel analogue models, detailing their theoretical properties and potential experimental realization using magnetic fields and fluid flows.

Findings

Hawking temperature depends on magnetic field, density, and tube radius.

Estimated Hawking temperature is approximately 0.0266 K.

Proposes laboratory setup for creating these analogue black holes.

Abstract

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\'{e}nic BHs, the magnetoacoustic and Alfve'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 Alfve'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-Alfve'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.

Equations

$$\frac{\partial \rho }{\partial t}+\nabla .(\rho \mathbf{v})=0,$$

$$\rho \left(\frac{\partial \mathbf{v}}{\partial t}+\frac{1}{2}\nabla {\text{v}}^2\right)=-\nabla p+\frac{1}{\mu }(\nabla \times \mathbf{B})\times \mathbf{B},$$

$$\frac{\partial \mathbf{B}}{\partial t}=\nabla \times (\mathbf{v}\times \mathbf{B}),$$

$$\nabla .\mathbf{B}=0,$$

$$p=k{\rho }^{\gamma },$$

$$\frac{{\partial }^2\mathbf{v}}{\partial t^2}=c_0^2\nabla (\nabla .\mathbf{v})+{\text{v}}_A^2(\nabla \times (\nabla \times (\mathbf{v}\times {\hat{B}}_0)))\times {\hat{B}}_0.$$

$${\omega }^2\hat{v}=c_0^2\mathbf{k}(\mathbf{k}.\hat{v})+{\text{v}}_A^2(\mathbf{k}\times (\mathbf{k}\times (\hat{v}\times {\hat{B}}_0)))\times {\hat{B}}_0.$$

$${\omega }^2\hat{v}=c_0^2{k}^2\hat{v}+{\text{v}}_A^2(\hat{v}-{\hat{B}}_0(\hat{v}.{\hat{B}}_0))k^2.$$

$$v_{ph}^2=\frac{{\omega }^2}{k^2}=c_0^2-{\text{v}}_A^2(\cos \theta -1).$$

$${\omega }^2\hat{v}={\text{v}}_A^2(\mathbf{k}\times (\mathbf{k}\times (\hat{v}\times {\hat{B}}_0)))\times {\hat{B}}_0.$$

$$\rho (\mathbf{r},t)$$

$$\mathbf{v}(\mathbf{r},t)$$

$$\mathbf{B}(\mathbf{r},t)$$

$$\frac{\partial {\rho }_0}{\partial t}+{\text{v}}_0\frac{\partial {\rho }_0}{\partial x}+{\rho }_0\frac{\partial {\text{v}}_0}{\partial x}=0,$$

$${\rho }_0{\text{v}}_0\frac{\partial {\text{v}}_0}{\partial x}=-c_0^2\frac{\partial {\rho }_0}{\partial x}$$

$$\nabla \times ({\text{v}}_0\times {\mathbf{B}}_0)=0,$$

$$\nabla .{\mathbf{B}}_0=0,$$

$$p_0=k{\rho }_0^{\gamma }.$$

$$\frac{\partial \psi }{\partial t}+\nabla \ln {\rho }_0.\nabla \phi +{\mathbf{v}}_0.\nabla \psi +{\nabla }^2\phi =0,\text{hfill}$$

$${\rho }_0\nabla \left(\frac{\partial \phi }{\partial t}+{\mathbf{v}}_0.\nabla \phi +c_0^2\psi \right)=\frac{1}{\mu }\left(\nabla \times {\mathbf{B}}_1\right)\times {\mathbf{B}}_0,\text{hfill}$$

$$\frac{\partial {\mathbf{B}}_1}{\partial t}=\nabla \times \left({\mathbf{v}}_0\times {\mathbf{B}}_1\right)+\nabla \times \left({\mathbf{v}}_1\times {\mathbf{B}}_0\right),\text{hfill}$$

$$\nabla .{\mathbf{B}}_1=0,$$

$$\frac{\partial }{\partial x}\left(\frac{\partial \phi }{\partial t}+{\mathbf{v}}_0.\nabla \phi +c_0^2\psi \right)=0.$$

$${\rho }_0\frac{\partial }{\partial z}\left(\frac{\partial \phi }{\partial t}+{\mathbf{v}}_0.\nabla \phi +c_0^2\psi \right)=\frac{B_0}{\mu }\left(\frac{\partial B_z}{\partial x}-\frac{\partial B_x}{\partial z}\right).$$

$$-\frac{\partial }{\partial t}\left(\frac{\partial \phi }{\partial t}+{\text{v}}_0\frac{\partial \phi }{\partial x}\right)+\frac{\partial }{\partial x}\left(c_0^2\frac{\partial \phi }{\partial x}-{\text{v}}_0\left(\frac{\partial \phi }{\partial t}+{\text{v}}_0\frac{\partial \phi }{\partial x}\right)\right)+c_0^2{\nabla }^2\phi =0.$$

$${\text{v}}_0\frac{\partial }{\partial x}\left(\frac{\partial \phi }{\partial t}+{\mathbf{v}}_0.\nabla \phi \right)+{\text{v}}_0\psi \frac{\partial c_0^2}{\partial x}$$

$$-c_0^2\left(\frac{\partial \psi }{\partial t}+\nabla \ln {\rho }_0.\nabla \phi +{\nabla }^2\phi \right)=0.$$

$$\frac{\partial }{\partial z}\left(\frac{1}{c_0^2}\frac{\partial }{\partial x}\frac{d\phi }{dt}\right)=\frac{1}{{\text{v}}_0}{\nabla }^2\frac{\partial \phi }{\partial z}+\frac{1}{{\text{v}}_0}\nabla \ln {\rho }_0.\nabla \frac{\partial \phi }{\partial z}\text{hfill}$$

$$+\frac{1}{{\text{v}}_0}\frac{\partial }{\partial t}\left(-\frac{B_0}{\mu {\rho }_0{c}_0^2}{\left(\nabla \times {\mathbf{B}}_1\right)}_y-\frac{1}{c_0^2}\frac{\partial }{\partial z}\frac{d\phi }{dt}\right)\text{hfill}$$

$$+\frac{B_0}{\mu {\rho }_0{c}_0^4}\frac{\partial c_0^2}{\partial x}{\left(\nabla \times {\mathbf{B}}_1\right)}_y+\frac{\partial }{\partial z}\left(\frac{1}{c_0^4}\frac{\partial c_0^2}{\partial x}\frac{d\phi }{dt}\right).\text{hfill}$$

$$\frac{\partial }{\partial t}\left(\frac{1}{N}\left(\frac{1}{c_0^2}\frac{\partial }{\partial x}\frac{d\phi }{dt}-\frac{1}{c_0^4}\frac{\partial c_0^2}{\partial x}\frac{d\phi }{dt}-\frac{{\nabla }^2\phi }{{\text{v}}_0}-\frac{\nabla \text{ln}{\rho }_0.\nabla \phi }{{\text{v}}_0}\right)\right)\text{hfill}$$

$$+\frac{\partial }{\partial t}\left(\frac{1}{N}\left(\frac{1}{{\text{v}}_0}\frac{\partial }{\partial t}\left(\frac{1}{c_0^2}\frac{d\phi }{dt}\right)-\frac{{\text{v}}_{\text{A}}^2}{c_0^2{\text{v}}_0}{\nabla }^2\phi \right)\right)=B_0{\nabla }^2\phi \text{, }\text{hfill}$$

$$\begin{array}{c}N=\frac{B_0}{\mu {\text{v}}_0}\frac{\partial }{\partial t}\left(\frac{1}{{\rho }_0{c}_0^2}\right)-\frac{B_0}{\mu {\rho }_0{c}_0^4}\frac{\partial c_0^2}{\partial x}\text{hfill}\\ =\frac{{\text{v}}_A^2}{B_0{c}_0^2}\left(\frac{\gamma }{{\text{v}}_0}+\frac{{\text{v}}_0}{c_0^2}\right)\frac{d{\text{v}}_0}{dx}.\text{hfill}\\ \text{hfill}\end{array}$$

$$-\frac{\partial }{\partial t}\left(\frac{\partial \phi }{\partial t}+{\text{v}}_0\frac{\partial \phi }{\partial x}\right)+\frac{\partial }{\partial x}\left(c_0^2\frac{\partial \phi }{\partial x}-{\text{v}}_0\left(\frac{\partial \phi }{\partial t}+{\text{v}}_0\frac{\partial \phi }{\partial x}\right)\right)\text{hfill}$$

$$+{\text{v}}_{\text{A}}^2\frac{{\partial }^2\phi }{\partial x^2}+{\text{(v}}_{\text{A}}^2+c_0^2)\left(\frac{{\partial }^2\phi }{\partial y^2}+\frac{{\partial }^2\phi }{\partial z^2}\right)=0.\text{hfill}$$

$$\frac{1}{\sqrt{-g}}{\partial }_{\mu }\left(\sqrt{-g}{g}^{\mu \nu }{\partial }_{\nu }\phi \right)=0,$$

$$\displaystyle{g^{\mu\nu}}=\frac{1}{{{\rho_{0}}c_{f}^{3}}}\left({\begin{array}[]{*{20}{c}}{-1}&{-{\text{}}{{\text{v}}_{0}}}&0&0\\ {-{\text{}}{{\text{v}}_{0}}}&{c_{f}^{2}-{\text{v}}_{0{\text{}}}^{2}}&0&0\\ 0&0&{c_{f}^{2}}&0\\ 0&0&0&{c_{f}^{2}}\\ \end{array}}\right),\,\,\,\,\,\hfill{g_{\mu\nu}}=-{\rho_{0}}{c_{f}}\left({\begin{array}[]{*{20}{c}}{c_{f}^{2}-{\text{v}}_{0{\text{}}}^{2}}&{{\text{}}{{\text{v}}_{0}}}&0&0\\ {{\text{}}{{\text{v}}_{0}}}&{-1}&0&0\\ 0&0&{-1}&0\\ 0&0&0&{-1}\\ \end{array}}\right),\,\,\,\,\,\,\,\,\,\hfill$$

$$d{s}^2=g_{\mu \nu }d{x}^{\mu }d{x}^{\nu }\text{hfill}=-{\rho }_0{c}_f\left((c_f^2-{\text{v}}_0^2)d{t}^2+2{\text{v}}_{0}dtdx-d{x}^2-d{y}^2-d{z}^2\right).\text{hfill}$$

$$d{s}^2=\text{hfill}{c}_f{\rho }_0\left(-\left(1-\frac{{\text{v}}_0^2}{c_f^2}\right)c_f^{2}d{\tau }^2+\frac{d{x}^2}{\left(1-\frac{{\text{v}}_0^2}{c_f^2}\right)}+d{y}^2+d{z}^2\right).\text{hfill}$$

Full text

black hole physics—waves

1 Introduction

In 1916 Schwarzschild gave a metric as a solution of the Einstein field equation. Singularity of such a metric predicted a gravitational black hole (BH) with event horizon at Schwarzschild radius (see Misner et al., 1973). 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 (Hawking, 1974, 1975). Hawking radiation (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 (Unruh, 1981, 1995). After 1981, most of attempts are proposed based on Bose-Einstein condensates of quantum fluid Garay et al. (2000, 2001); Barceló et al. (2001); Recati et al. (2009); Zapata et al. (2011), quasi particles in superfluid (Volovik, 2000), ultra-cold fermions (Giovanazzi, 2005), in plasmas and ion rings Horstmann et al. (2010); Blatt & Roos (2012); Gerritsma et al. (2010, 2011); Casanova et al. (2011), slow light in an atomic vapor Brevik (2001); Leonhardt & Piwnicki (2000); Leonhardt (2002); Unruh & Schützhold (2003), in water (Rousseaux et al., 2008; Weinfurtner et al., 2011), etc. Recently, observation of self-amplifying HR in an analogue BH laser suggested a very promising experiment method for probing the inside of a BH (Steinhauer, 2014). From a theoretical point of view, the acoustic analogue BH models are developed in geometrical acoustics and physical acoustics (Barceló et al., 2011). 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 (Sakagami & Ohashi, 2002). This was interpreted as characteristic of a sonic BH (Barceló et al., 2011; Sakagami & Ohashi, 2002). 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 (e.g., (Goedbloe et al., 1973)). After 1942 (Alfvén, 1942), MHD (Alfvén and magnetoacoustic) waves have been detected, using the laboratory experiments (e.g., (Lundquist, 1949; Bostick & Levine, 1952; Lehnert, 1954)), in the earth-ionosphere and magnetic field (Berthold et al., 1960; Masahisa, 1961), a variety normal modes of the solar corona Aschwanden et al. (1999); Nakariakov et al. (1999); Wang et al. (2002); Aschwanden (2004); McIntosh et al. (2011). 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: Sec. 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. Sec. 3 introduces the acoustic, magnetoacoustic, and Alfvénic black holes, respectively. Sec. 4 calculate the Hawking temperature for the acoustic, magnetoacoustic, and Alfvénic black holes. Sec. 5 describes the conclusions.

2 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 (Priest, 2014). The ideal MHD equations for an adiabatic process and irrotational flow $\nabla\times{\bf v}=0$ are given by:

[TABLE]

where, $\rho$, $p$, ${\bf v}$, ${\bf B}$, $\mu$, $\gamma$, 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.

2.1 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 ($\rho_{0}=cte$), constant pressure ($p_{0}=cte$), and uniform background magnetic field (${\bf B}_{0}$), were investigated in the literature (e.g., (Priest, 2014)). The linearised ideal MHD equations can be reduced to a single equation for disturbed velocity as (Priest, 2014),

[TABLE]

where $c_{0}=\sqrt{\frac{\gamma p_{0}}{\rho_{0}}}$ is the sound speed, $\text{v}_{A}={\frac{B_{0}}{\sqrt{\mu\rho_{0}}}}$ is speed of the Alfvén wave, and ${\hat{B}_{0}}={\bf 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 (Gibbs & Wilson, 1901), can be decomposed to an irrotational part ($\nabla\phi$ where $\phi$ is a scaler) and a solenoidal part ($\nabla\times{\bf A}$ where ${\bf A}$ is the vector potential and satisfy $\nabla.{\bf A}=0$ ). The irrotational (gradient) and divergence-free solutions can be used for treating the longitudinal and transversal waves, respectively.

Suppose an irrotational solution, $\bf{v}=\nabla\varphi$ with $\varphi=\tilde{\varphi}\exp{(i{\bf k}.{\bf r}-i\omega t)}$ for Eq. (6), immediately we find ${\bf v}=\hat{v}\exp{(i{\bf k}.{\bf r}-i\omega t)}$ where $\bf k$, $\omega$, $\tilde{\varphi}$, and $\hat{v}=i{\bf k}\tilde{\varphi}$, are wave vector, angular frequency of oscillations, a constant, and wave amplitude, respectively. Furthermore Eq. (6) yields,

[TABLE]

Equation (7) simplifies as

[TABLE]

We note that the wave vector ($\bf k$) and wave amplitude ($\hat{v}$) are parallel (${\bf k}\|\hat{v}$). Considering the wave propagation parallel to the background magnetic field $({\bf k}\|{\bf B}_{0})$, Eq. (8) reduces to the dispersion relation $\omega^{2}=c_{0}^{2}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 (${\bf k}\bot{\bf B}_{0}$), Eq. (8) gives the dispersion relation $\omega^{2}=({\text{v}_{A}^{2}+c_{0}^{2}}){k^{2}}$. This is the well-known characteristic of the fast magnetoacoustic wave with the phase speed $c_{f}^{2}={\omega^{2}}/{k^{2}}=c_{0}^{2}+\text{v}_{A}^{2}$. For the oblique propagation ${\bf k}.{\bf B}_{0}=kB_{0}\cos\theta~{}(\theta$ is the angle between ${{\bf B}_{0}}$ and ${\bf k})$, the phase speed of the slow magnetoacoustic wave is given by

[TABLE]

From the above analysis we see that, choosing the disturbed velocity as an irrotational velocity field $\bf v=\nabla\varphi$, the Alfvén waves cannot propagate along the background magnetic field (${\bf B}_{0}$).

Suppose a divergence-free solution, $\bf{v}=\nabla\times{\bf A}$, with a planar wave solution ${\bf A}={\hat{A}}\exp{(i{\bf k}.{\bf r}-i\omega t)}$, $\hat{v}=i{\bf k}\times\hat{A}$ and $\hat{A}$ is a constant vector, Eq. (6) gives,

[TABLE]

We see that the velocity amplitude is perpendicular to the wave vector (${\bf k}.{\hat{v}}=0$), and only the transversal Alfvén wave with phase speed $v_{ph}=\text{v}_{A}$ can propagate.

In the remainder of this section, the metrics for acoustic, magnetoacoustic, and Alfvén longitudinal waves (irrotational solutions) are derived.

2.2 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 ${\bf v}=\nabla\varphi$. We consider small perturbations from equilibrium as

[TABLE]

where, equilibrium quantities indicated by subscript ”0” are function of position $x$ and time $t$, $\psi({\bf r},t)$, $\varphi({\bf r},t)$, and ${\bf B}_{1}({\bf r},t)$ are perturbed quantities (Safari et al., 2007; Edwin & Roberts, 1982). Equilibrium quantities ($\rho_{0},\text{v}_{0},{\bf B}_{0},$ and  $p_{0}$) are satisfied by

[TABLE]

Linearization of Eqs (1)-(5) (products and squares of the small perturbations are neglected) and after some mathematical manipulations, give

[TABLE]

As we explained in the previous section, by choosing ${\bf v}=\nabla\varphi$ the parallel propagation of Alfvén waves along the magnetic field $({\bf B}_{0}=B_{0}\hat{x})$ is absent. Therefore, the first term on right hand of Eq. (21) is set to zero, $\nabla\times\left({\mathbf{v}}_{0}\times{\mathbf{B}}_{1}\right)=0$. A combination of this assumption and the solenoidal condition for magnetic field $(\nabla.{\bf B}_{0}=0)$, gives $\frac{\partial B_{1x}}{\partial 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 ${{{\mathbf{v}}_{0}}\times{{\mathbf{B}}_{1}}}=\nabla\delta(x,y,z,t)$ in which $\delta(x,y,z,t)$ is a function.

The $x$ and $z$ component of Eq. (20) are

[TABLE]

[TABLE]

Here after, first, we focus our analysis on the derivation of acoustic waves with propagation along the background magnetic field directions $(\hat{x})$ and second, magnetoacoustic wave propagation in all directions except the magnetic field direction.

First, combining Eqs (19) and (23) gives

[TABLE]

Equation (25) is the well-known Klein-Gordon equation for acoustic waves. Eliminating $\partial\psi/\partial x$ between Eq. (19) and Eq. (23) one obtains

[TABLE]

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. (2.2) with respect to $z$ and substituting $\partial\psi/\partial z$ from Eq. (24) one finds

[TABLE]

in which, $d/dt=\partial/\partial t+{{\mathbf{v}}_{0}}.\nabla$. Substituting $-\frac{\partial}{{\partial t}}{\left({\nabla\times{{\mathbf{B}}_{1}}}\right)_{y}}={B_{0}}{\nabla^{2}}\frac{{\partial\varphi}}{{\partial z}}$ from Eq. (21), into Eq. (27) and differentiating with respect to $t$ using $\frac{\partial}{{\partial t}}{\left({\nabla\times{{\mathbf{B}}_{1}}}\right)_{y}}=-{B_{0}}{\nabla^{2}}\frac{{\partial\varphi}}{{\partial z}}$ we derived a single equation for velocity potential $(\varphi)$

[TABLE]

where,

[TABLE]

For a flow having a slight change in the speed $(d{{\text{v}}_{0}}/dx\ll 1)$, for high frequency waves (short period $\Delta t\ll 1$) the term $\Delta t(d{{\text{v}}_{0}}/dx)$ becomes too small. In this case, the right hand side term ${B_{0}}{\nabla^{2}}\varphi$ of Eq. (2.2) can be negligible compared to the last term $\frac{\partial}{{\partial t}}\left({\frac{{{\text{v}}_{{\text{A}}}^{2}}}{{Nc_{0}^{2}{{\text{v}}_{0}}}}{\nabla^{2}}\varphi}\right)$ of the left hand side. This leaves an equation for $\varphi$

[TABLE]

Equation (2.2) 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. (2.2) 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 (Barceló et al., 2011).

Equation (2.2) can be reformulate in a tensorial form

[TABLE]

where, $g^{\mu\nu}$ and $g$ are inverse metric tensor and its determinant, $\mu$ and $\nu$ runs from 0 (indicates the time coordinate) to 3 (denotes the spatial coordinate). The magnetoacoustic inverse metric tensor $g^{\mu\nu}$ and metric tensor $g_{\mu\nu}$ are obtained as

[TABLE]

where, $c_{f}^{2}={\text{v}}_{{\text{A}}}^{2}+c_{0}^{2}$. Using metric tensor, Eq. (39), the magnetoacoustic interval can be defined as

[TABLE]

Inserting the specific time interval $d\tau=dt+\frac{{{{\text{v}}_{\text{0}}}dx}}{{c_{f}^{2}-{\text{v}}_{0{\text{}}}^{2}}}$ into Eq. (40), one obtains

[TABLE]

In the following section, the properties of acoustic, magnetoacoustic, and Alfvénic BHs are investigated.

3 Acoustic, Magnetoacoustic, and Aflvénic black holes

3.1 Acoustic black hole

The theory of gravitational BHs has been developed into the supersonic flow by Unruh (Sakagami & Ohashi, 2002). 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” (Novello et al., 2002). In the case of unmagnetised gas($B_{0}=0$), Eq. (2.2) 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 (Sakagami & Ohashi, 2002).

The acoustic metric can be obtained by setting $\text{v}_{A}=0$ in Eq. (39)

[TABLE]

Equivalently, the acoustic interval can be expressed as

[TABLE]

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

[TABLE]

This relation shows that for $(dS<0)$ a subsonic flow $({{\text{v}}_{0}}<{c_{0}})$ will be accelerated and a supersonic flow $({{\text{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. (47) 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 Figure 1 has been proposed to discuss the above mentioned acoustic BH. This setup uses, an axisymmetric sufficiently thin tube with slightly decreasing cross section $\left({S(x)}\right)$ that reaches its minimum cross section at the throat and then slightly increases that. An initial material flow $({{{\bf v}}_{0}}={{\text{v}}_{0}}\left(x\right)\,\hat{x})$ along the tube axis is considered.

3.2 Magnetoacoustic black hole

The magnetoacoustic metric Eq. (41) is singular at the magnetoacoustic point, where ${\text{}}{{\text{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 Figure 2 to discuss the above mentioned BH. The setup consists of an axisymmetric sufficiently thin tube with slightly variable cross section $\left({S(x)}\right)$, a uniform force free magnetic field ${\bf B}_{0}={{\text{B}}_{0}}\hat{x}$ and an initial material flow $({{{\bf v}}_{0}}={{\text{v}}_{0}}\left(x\right)\,\hat{x})$ along the tube axis Figure 2.

A similar treatment of sub and supersonic flow in tube configuration can be explained for sub and super-magnetoacoustic flow based on Eq. (48). In other words, the super-magnetoacoustic flow $({{\text{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 (${{\text{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 (Steinhauer, 2014; Schley et al., 2013). 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.

3.3 Alfvénic black hole

In the limit of ${{\text{c}}_{0}}<<{{\text{v}}_{\text{A}}}$ (zero $\beta$ plasma condition), the magnetoacoustic wave (Eq. 2.2) 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 39 and 41) by setting $c_{0}$ tends to 0,

[TABLE]

Equivalent Alfvénic interval is given by

[TABLE]

The Alfvénic metric (Eq. 54) is singular at the location of Alfvénic point where ${{\text{v}}_{0}}={{\text{v}}_{\text{A}}}$. This singular behaviour of Alfvénic metric leads to an Alfvénic BH. We propose a setups that is illustrated in Figure 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){\text{v}_{0}}(x)$, is constant in the tube cross section. Consider a super-Alfvénic flow in a tube with increasing cross section Figure 3, speed of the flow decreases along the tube axis and reaches to Alfvén velocity $\text{v}_{0}=\text{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.

4 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

[TABLE]

where, $\hbar=h/2\pi$, $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,

[TABLE]

In the limit of $\text{v}_{A}$ tends to zero, Eq. (56) then reduces to Hawking temperature for acoustic BH. Although, in the zero $\beta$ plasma condition the above mentioned equation can describe the Hawking temperature for Alfvénic BH,

[TABLE]

At the horizon where the tube cross section radius is equal to $R$, Eq. (56) can be simplified as

[TABLE]

where, the term $\frac{1}{{{c_{f}}}}\frac{d}{{dx}}({c_{f}}-{{\text{v}}_{0}})$ is approximately equal to $\frac{1}{R}.$ For a plasma with a ratio of $\chi={c_{0}}/{{\text{v}}_{A}}$ in which $\chi$ is a positive number, Eq. (58) gives

[TABLE]

where, $n$ is the number density of plasma and all units are in SI. For a typical plasma with magnetic field strength $B_{0}$=1 Tesla, number density $n={10^{18}}\,{m^{-3}}$, $R=1$ mm, and ${\chi^{2}}\ll 1$, the Hawking temperature is estimated about 0.0266 K. For a typical natural fluid in a Lavel nozzle experiment, the Hawking temperature was estimated about $10^{-6}$ K (Novello et al., 2002).

5 conclusion

In this study, the behavior of MHD waves (Alfvén and magnetoacoustic) are are treated in the Helholtz decomposition frame-work, in the case of irrotational flow, 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.