Analogous Hawking Effect: S -Matrix and Thermoﬁeld Dynamics

work. Abstract: We consider the full S -matrix in the scattering giving rise to analogous Hawking radiation in dispersive media. We show the general structure of the scattering in the weak dispersion approximation and discuss some unnoticed features of the primary process, with a possible generalization of the phenomenology of the Hawking effect. In particular, we stress that the Hawking particle and its antiparticle partner a priori could also be produced with different rates. We provide a general parameterization of the S -matrix, adopting the Iwasawa decomposition for the matrix itself. Then, we assume that a perturbative structure in a suitable sense is allowed and display the corresponding expansion. In connection with the general structure of the S -matrix at the leading order, we also consider the thermoﬁeld dynamics (TFD) framework and show that the TFD picture is still available, with a doubling of the degrees of freedom emerging in a natural way, as for the astrophysical black hole case. Furthermore, we show that particles on the thermal vacuum can be identiﬁed with real particles appearing in the scattering.


Introduction
In [1,2] we proposed a fourth-order ordinary differential equation as a master equation allowing to deal with the analogous Hawking effect in condensed matter systems in a systematic way, in the approximation of weak dispersive effects. Systems between the most studied and relevant ones, such as Bose-Einstein condensates (BEC), dielectric media, and water are encompassed in our picture. The only change is in the identification of the weak dispersion parameter appearing in the master equation.
Herein, we take into consideration the full structure of the S-matrix, which is associated with the analogous Hawking effect in condensed matter. The S-matrix has been discussed previously in the literature, see, e.g., [3][4][5][6][7], for BEC, water, and for the Corley model [8]. In particular, in a general picture encompassing the previously discussed ones, we show that the S-matrix, at the leading order in , is very simplified with respect to the original complexity of the complete scattering process, due to the fact that two modes participating in the scattering substantially decouple, at least as far as the primary process, i.e., the pair-creation process, is concerned (they still can interact at the level of subsequent scattering with respect to the primary process). The manifest advantage of the proposed framework [1,2] consists in the fact that a unified picture for many interesting cases occurring in the physical literature can be provided, and this is true also at the level of the structure of the S-matrix. Moreover, by our analysis it is also clear that, whenever a direct participation of the aforementioned modes in the primary process would occur, one could maintain thermality for the spectrum of the particles detected at infinity. However, there would occur a difference in the number of created particles, to be identified with the Hawking particles, and the would be Hawking partners (would-be anti-particles, where "would be" stresses their possible different production ratio w.r.t. particles), and then a generalization of the Hawking process would occur, with a thermal spectrum but separate behavior between the Hawking particle and its would-be partner. This is also a feature that could be a very interesting subject of experimental verification. In particular, a joint measurement of the Hawking particle could represent in any case a corroboration of the general picture for the analog Hawking effect.
Further, we take into account a series expansion of the S-matrix, as we are interested in looking for corrections to the leading structure S-matrix due to subleading processes. After a naive approach, we adopt a rigorous picture involving the Iwasawa decomposition of the unitary group associated with the S-matrix itself, by means of which we can corroborate the aforementioned naive approach. The proposed Iwasawa decomposition represents a novel tool in the analysis of the S-matrix and will be the subject of further analysis.
As an interesting consequence of the above-mentioned picture involving decoupling of the two modes, we show that also in the analog Hawking effect there exists still the possibility to associate a thermofield dynamics framework for the primary process. One can again show that the doubling of the number of degrees of freedom, with the introduction of a (would-be fictitious) Hilbert space, arises in a natural way, as in the case of astrophysical black holes, as Israel pointed out [9]. Furthermore, we show that particles over the thermal vacuum can also be identified with real particles appearing in the scattering and detectable in the final scattering states, as far as the white hole situation is concerned. For the black hole case, they still appear as ingoing states, and we point out that experiments involving stimulated Hawking radiation in the black hole case have been performed in [10].

The Master Equation: An Orr-Sommerfeld Type Fourth-Order Equation
Three significant cases of wave equations in dispersive analogue gravity can be reconduced to the equation where the upper sign occurs in the case of subluminal dispersion and the lower one in the case of superluminal dispersion. The latter case is considered in Nishimoto's works (see, e.g., [11] and references therein). Furthermore, is assumed. As → 0, one finds the reduced equation Solutions of p 30 (x) = 0 (4) define the turning points (TPs) of the equation, and in the analysis of the reduced equation, the behavior of solutions in the neighborhood of the TPs is of utmost relevance for the scattering problem we mean to delve into. In the following, we limit ourselves to the case of a single TP, which can be identified with x = 0 without loss of generality. In [11], it is assumed that the reduced equation displays a Fuchsian singularity at the TP (nothing actually prevents the general equation in itself to admit a regular behavior).
It can be shown [11] that near the TP x = a the original equation is replaced by the fourth-order differential equation where, as discussed at length in [1], w(z) represents the (rescaled [1]) wave-function in the near horizon approximation, to be then matched in the linear region, with WKB solutions that hold far from the turning point. We refer the reader to [1] for further details. In (5), the upper sign is for the subluminal case, and the lower one is for the superluminal one, and and Solutions to Equation (5) are found by means of Laplace integrals with a suitable choice for the paths C j in the complex t-plane. Of course, there is also a constant solution, which appears to be a Dirac delta in the space of the Laplace transform. This further solution is fundamental in order to obtain a complete basis for the scattering problem also in the region near the TP (i.e., near the horizon). As we discussed in [1,2], Equation (5) is universal in form and is at the root of the analogous Hawking effect in dispersive media, as far as they are governed by the above Orr-Sommerfeld-like equation.
In the following, we shall limit ourselves to the explicit discussion of the subluminal case. For the superluminal one, we refer the reader to [2] for further details. See also the following section. For the specific case of dielectric media, see also [12].

The S-Matrix in Presence of Dispersion
As known, the spectator mode v does not participate in the process at the turning point, at least as far as the weak dispersion limit 0 < 1 is assumed, where is the parameter involved in the weak dispersion approximation, with = 0 indicating the absence of dispersion. One may construct the straddling mode, which is obtained by analytic continuation of the Hawking mode in the black hole region, to be identified with the l-mode. The analytic continuation is unique and involves also the choice of the branch, i.e., of the sheet of the Riemannian surface related to the logarithmic branch point associated with the solutions of the differential equation governing the scattering process.
To be more precise: let us first individuate the OUT modes of the scattering, and in particular let us call the u-mode the Hawking mode in the external region, the l-mode its partner in the black hole region, and d-mode the (regular) mode again in the black hole region. The l-mode has a negative norm (antiparticle). As far as the IN modes are concerned, let us call the p-mode and the n-mode the high wave-number modes entering the black hole, where the n-mode is a negative norm one. Last, but not least, we call v-mode the (regular) mode entering the black hole. See Figure 1 for a pictorial representation.
As we have discussed in [1,2], at least at the leading order in the weak dispersion approximation, the v-mode substantially decouples from the other ones, leaving a very simple structure for the S-matrix, which will be discussed in the following subsection.

The Complete S-Matrix
The following scattering operator is taken into account in literature (see, e.g., [6] for the subluminal case): if one considers as I N-modes the ones moving towards the TP and as OUT-modes the ones moving towards infinity (i.e., x → ±∞), one may write where and One has In order to allow a better visualization of the overall process, we represent schematically in Figure 1 the modes participating to the scattering. In particular, we display the subluminal case for the black hole scattering. The IN modes are traveling towards the horizon, whereas the OUT ones are traveling towards the asymptotic regions (x → ∞ in the asymptotic external region and x → −∞ in the deep black hole region). In the superluminal case, the only change occurring is that the p-mode and the n-mode travel towards the horizon from the black hole region.  One obtains, by taking into account that S ∈ U(1, −1, 1), The above formulas amount to current conservations. As it is evident, there are three current conservations occurring for the processes at hand. The first formula amounts to the current conservation (13) for the Hawking process. As discussed in [1,2], in the near-horizon approximation one can individuate for the Hawking process a diagram in the Laplace transform space where the particles involved in the aforementioned current conservation appear. One should expect that also for the other processes involved in the other two current conservations where analogous diagrams can be drawn. As we noticed in [1], for the Hawking process the v-mode is not enabled to participate in the process directly, i.e., at the level of the near-horizon process. For the sake of completeness, we also present in Figure 2 the diagram involved in the Hawking process in the subluminal case, which was called Corley's diagram in [1]. See [8] and also [4,13]. For the superluminal case, there are only little changes. See [2,4,8,13].   [8]. Given the universality of (5), this diagram holds true for all the subluminal cases, which can be rewritten in terms of the Orr-Sommerfeld-like equation discussed in [1]. The paths for the dispersive modes, called the p-mode and the n-mode in the text, are indicated, as well as the one for the Hawking mode (u-mode) and for the decaying mode. The last mode is the one in the black hole region x < 0. As remarked by Corley [8], the path for the decaying mode can be deformed in the paths for the Hawking particle and for the dispersive modes. For the superluminal case, we refer the reader to Figure 1 displayed in [2].
For the Hawking mode at frequency ω, one obtains (cf. [1]) where κ corresponds to the surface gravity for the analogous black hole one is considering (see [1,2]), whereas for the n-mode and the p-mode one has [1]

The Hawking Partner
One may obtain the Hawking partner by choosing a suitable analytical continuation for x < 0. It turns out that, by choosing the branch where −1 = e −iπ , for the further solution one obtains As to the latter mode, the second line (14) represents the process involving the Hawking partner, i.e., the l-mode, which is specular with respect to the one for the Hawking particle (u-mode). As to the associated diagram, there are some subtleties that can be suitably pointed out. Of course, as the number of Hawking modes can be obtained by considering the ratio for the l-mode, to be identified with the Hawking partner, one must get |β partner ω and this expression must coincide with the former one, if the identification is correct.
According to the discussion in [13], in order to maintain the same ratio between the two currents in the two latest equations, it is necessary to exchange the roles of the n-mode and the p-mode, and this can be obtained by shifting the cut on the positive real axis. See  We follow a different strategy, as we maintain the identification of the Hawking partner with the mode that is obtained by analytic continuation at z < 0 from the Hawking mode, as indicated above [1] (where we chose −1 = e −iπ ). Being an antiparticle, its relation is with the conjugate solutions corresponding to the p-mode and to the n-mode (so that the former is in its antiparticle state and the latter in the particle one). Then, one obtains As a consequence, one is able to maintain the same ratio for (20) and for (21). We can also obtain the same result with a different choice. It is straightforward to show that the l-mode can be obtained simply by taking into account a different leaf of the logarithm log(z) (or, equivalently of the power function z a , with a being non-rational): The principal branch is given by m = 0. Let us consider the branch m = −1. We check if the same Corley diagram as in Figure 2, but this time with the path around the cut representing the Hawking particle with x < 0, can represent the right physical situation, now on the branch m = −1. This time, we choose −1 = e iπ , and it is easy to show that we get again (23) for the Hawking partner. What happens in the case of the p-mode and n-mode? They coincide with the ones indicated in (22) and (23), respectively, and then the ratio in (20) and (21) is the same. From the physical point of view, one would wonder if there is the necessity to pass to the different leaf as we did. In a certain sense, given that in the non-dispersive Hawking process the Hawking particle and its partner are pair-created, one would expect the same diagram for both the modes in the dispersive case, at least as far as u-modes and l-modes can be considered as the partners of a pair. This gives a nice support to the picture involving the m = −1 branch. At the same time, one can also interpret the particles appearing in that branch as belonging to the would-be fictitious Hilbert space involved in the TFD picture of the phenomenon. See also the following section.

Generalized Hawking Effect in Analogue Gravity
There is also a nice consequence to be immediately pointed out: let us write the number of created particles for the two aforementioned cases, by taking into account that the signal of thermality for the Hawking process is in general assumed to be given by where β h := 2π/κ is the inverse of the analogous Hawking temperature, and, analogously, also for the Hawking partner one must expect One obtains respectively. A first interesting consequence is that, at the level of particle creation, the identification of the l-mode as the Hawking partner, which is then created at the same rate as the Hawking mode, requires that at least for what concerns the primary process, i.e., we mean the process of pair-creation regardless of any subsequent scattering depleting the flux of the particles arriving at infinity: It might happen that subsequent scattering is able to deplete the flux of Hawking particles at infinity, still leaving unaltered the above balance for the primary pair-creation process. This is a very interesting consequence of the overall picture concerning the analogous Hawking effect in dispersive media. Indeed, the above considerations are modelindependent, as they hold true in general for a dispersion relation like the one emerging for any model that is governed by a fourth-order ordinary differential equation, and, as far as (25) and (26) are assumed to be implemented, as it is commonly done at least in the weak dispersion approximation, any deviation from pure thermality at the level of the primary process would clearly be of hindrance to the identification of the l-mode as the Hawking partner, as a different rate of creation would occur. To some extent, this is not a tricky problem, as the thermality one is interested in would still be ensured. However, it is remarkable that one should also stress that a process where the aforementioned difference in the rate of particle creation occurs would be a generalized Hawking effect. It is remarkable also that the Hawking effect discussed, e.g., in [14,15] is (probably?) to be ascribed to the latter generalization, as a grey-body factor arises at the level of the primary process.

The Third Process
One has further interesting consequences for the third process for the third current. Additionally, the d-mode cannot participate in the process (it corresponds to the constant solution in the near-horizon approximation [1]), as well as the v-mode itself. There is no chance to draw a diagram involving also the modes p, n involved in the third current conservation. As a consequence, one must also find at the leading order which implies also |σ dv | = 1.

The S-Matrix
As a consequence of the weak dispersion ansatz discussed in [1], one obtains at the leading order in a trivial decoupling of the S-matrix: It is easy to show that, and then, at the leading order, Analogously, one finds and then, at the leading order, As already noticed in the literature, the above S-matrix satisfies the equation where This condition arises from the following considerations: we are actually considering the Bogoliubov transformation B from the initial scattering state to the final one. As is well-known (see, e.g., [16,17]) in the case of a scalar field one has where one defines the operator q : and where P ± are the projectors on the positive/negative norm solutions, respectively. This condition is easily seen to be equivalent to (37). It may be noted that, from the point of view of quantum field theory (QFT), the Smatrix discussed above is different from the usual S-matrix of the standard QFT, albeit there is a direct relation between them. In QFT, there is usually another definition of S-matrix: where we have indicated with S the scattering operator of QFT. As well-known, it has to be an unitary operator: S † S = SS † = I, such that |0, out = S|0, in holds, which is ensured if a finite number of particles is created in the transition between the two aforementioned vacua [16,18,19].

A Perturbative Approach to the S-Matrix
For simplicity of notation, henceforth we shift to a more standard writing of the scattering matrix. The S-matrix has entries s ij . Condition (37) leads to the following six equations: |s 31 | 2 − |s 32 | 2 + |s 33 | 2 = 1, One may explore the next-to-leading-order correction in by assuming an analytic behavior of the S-matrix in the parameter . This might be a non-trivial assumption, as our analysis is involved with quantum field theory in external fields, where enters in a non-trivial way. Then, we still maintain an expansion parameter, say , leaving open the problem to find out the relation between and . Indeed, in our case, our fourthorder equation of the Orr-Sommerfeld kind is involved in a singular perturbation theory problem, as higher-order derivatives appear in the perturbation itself, and this represents a nontrivial mathematical problem. It is, e.g., not easy to determine in a complete way the first correction to the wave functions in the near-horizon region. As to general difficulties occurring in the case of singular perturbation theory, see, e.g., Examples 1.19 and 1.20, pp. 435-436 in [20]. We proceed as in [21], with the difference that our expansion parameter is : We put: It is supposed that s 13 , s 23 are at least order of , i.e., are small (cf. [21]). With simple algebraic manipulations of Equations (46) and (47) we get the following equations: and then we obtain, as in [21], i.e., s A less-naive approach to the parameterization of the S-matrix is presented in the following subsection. As to the meaning of the parameters s 23 , it is better to go back to the previous notation: one has s (1) uv , which corresponds to the amplitude, at the first order, of conversion of an initial v-mode to a final u-mode (Hawking particle), whereas s (1) lv corresponds to the amplitude, at the first order, of conversion of an initial v-mode to a final l-mode (Hawking partner). We recall that we indicated the l-mode as would-be Hawking partner because of the possibility to find a different amplitude of production of the l-mode with respect to the amplitude of production of the u-mode, as previously discussed. This difference does not arise in our framework of -expansion. This is obviously also true in the above expansion in . A relation (yet-unknown) = h( ), where h is a function to be determined, is expected.

A General Parametrization of the S-Matrix
As we said, the S-matrix constructed above is an element of the group U(1, −1, 1), a real form of U(3, C) preserving the quadratic Hermitian form η :=diag(1, −1, 1) in C 3 . Its general element depends on nine real parameters and can be realized by means of the Iwasawa parametrization [22,23]. It means that the generic form of S can be put in the following form: where K is a generic element of the maximal compact subgroup; A is a maximal noncompact abelian subgroup; and N is the exponentiation of the nilpotent generated by the positive roots associated to A. Let us do it in detail, starting from the associated Lie algebra u(1, −1, 1). It consists of the 3 × 3 matrices A that satisfy Accordingly, we choose the following basis for the Lie algebra: Notice that they are trace orthogonal but not uniformly normalized. This choice simplifies the expressions we need. The first 5 matrices are compact (their squares have negative trace) and generate the maximal compact subgroup K. It is a U(2) × U(1): the first four matrices generate U(2), of which the first three generate SU (2), and the fifth matrix generates the U(1). The quotient U(1, −1, 1)/(U(2) × U(1)) has rank 1, so there is only one generator for A. We choose µ 6 . Finally, we have to determine the eigenmatrices of the adjoint action of µ 6 corresponding to positive roots. These are easily obtained by computing the commutators [µ 6 , µ j ], j = 7, 8, 9. We get Parametrizing the subgroup SU(2) à la Euler, we get with More explicitly The range of the parameter to cover the whole group is Exploiting all the products would provide us the most general parametrization for coefficients of S. However, for several analyses it may be convenient to keep the terms factorized. As an example, let us consider again the case of a perturbative deviation from (32). To this end, we can write with sinh y ω = β ω . We could compute the most general first-order correction, but we want now to limit ourselves to consider only real corrections, so we keep and keep the first order in . This gives S = S (0) + S (1) + · · · , with This expression is more general than (55). While the first 2 × 2 block coincides with the previous one, the last row and column are more general, since they are not forced to be symmetric, and so we have three independent parameters in place of just two. Notice that imposing S (1) 32 in (71) automatically gives symmetry and reproduces exactly (55). We stress that, while the comparison we have just done seems to suggest a direct identification = , such a conclusion would be too naive as one has to determine the final identification in the full quantum field theory framework.
Our parametrization allows to include easily also the corrections due to the non-real terms (by switching on x 1 , x 3 , x 4 , x 5 , z 1 , z 3 ) or to extend the perturbation at any given order or to provide exact general expressions, after relating the generalized Iwasawa parameters with the physical ones.

Thermofield Dynamics and Thermality
From the origins of the Hawking effect and of the thermofield dynamics (TFD) framework for describing thermal equilibrium, there exists a fine relation between thermal emission by black holes [9] and thermofield formalism [24][25][26][27]. As discovered by Werner Israel, the thermal vacuum of TFD coincides with the Hartle-Hawking state, sometimes also indicated as a Hartle-Hawking-Israel state in tribute of Israel's finding [9]. See also [13,28,29]. Our possibility to show that the p-modes and n-modes correspond to the particle/antiparticle states of the thermal vacuum, i.e., as particles of the Hartle-Hawking-Israel vacuum, is a noticeable consequence of the aforementioned scheme.

A Summary of TFD
Thermofield dynamics are introduced with the aim of describing quantum field theory at a finite temperature. The key idea is to set up a formalism where the ensemble average of any operator at thermal equilibrium is replaced by an expectation value of the same operator with respect to a thermal vacuum: Tr(exp(−βH)A) =: 0(β)|A|0(β) .
In the previous equation, Z(β) is the partition function, and H is the Hamiltonian operator for the system at hand; |0(β) is the thermal vacuum, and β is the inverse temperature, as usual. It is not possible to obtain that |0(β) stays in the same Hilbert space H as the bosonic theory at hand, but, in order to implement the above mapping, it is necessary to duplicate the degrees of freedom by introducing a fictitious copyH of the original Hilbert space, with associated creation and annihilation operators. The thermal state |0(β) can be arranged for in the tensor product H ⊗H. Let us introduce in the vacuum states |0 and |0 , with corresponding annihilator operators a l ,ã l and b l ,b l for particles and antiparticles, respectively. A physical (fictitious) n-particle state whose energy is E n is denoted by |n (|ñ ) and |n,ñ = |n |ñ is the n-particle tensor product state.
The vacuum state in the extended Hilbert state is Z(β) is a normalization.
In the calculation of (72), one has to sum over all the fictitious statesñ; the result is One can also introduce thermal-state operators, according to the standard constructions in thermofield dynamics, and introduce the thermal state |0(β) and thermal-state annihilation operators a l (β),ã l (β), b l (β),b l (β), which are such that Let us recall that, for bosons, we have E n = nω. We are mainly interested in the following Bogoliubov relations: and the analogous ones for b-operators (which correspond to operators for antiparticles, i.e., for negative frequency states; cf. [24]), with

Analogue Hawking Effect and TFD
In the framework of analogue gravity, one may show that particle of the thermal vacuum can be identified with particles appearing in the spectrum of the analogous system displaying the analogous Hawking effect. In agreement with the leading term in the Smatrix we discussed above, we have Γ ω = 1, where Γ ω stays for the grey-body factor (see (85) below), and we have restored its standard notation for the Hawking effect itself. The decoupling of the modes v, d at the leading order, as discussed the previous section, is fundamental in order to allow a thermofield dynamics picture also in analogue gravity, as we show below. Indeed, as a consequence of such a decoupling, we can write a u ω = α ω a p ω + β ω a n † ω , to be compared with the thermofield formula where The fictitious Hilbert space of TFD is instead populated by the Hawking partner: to be compared with the thermofield formulã As a consequence, we can claim that the p-modes and the n-modes correspond to particle states associated with the thermal vacuum, to be compared with the Hawking mode (u-mode) and the Hawking partner (l-mode), which are states of the "standard" (i.e., non-thermal) vacuum, and belong to the Hilbert space and to the would-be fictitious Hilbert space, respectively. Further, we note that, in the case t of a white hole, modes p, n are outgoing scattering states; in the black hole case, modes p, n are ingoing scattering states, which can be associated with vacuum fluctuations (as far as the spontaneous process is concerned) and real particles in the stimulated case. Then, it is worth mentioning that, in both the white hole case and the black hole one, p, n may be experimentally measured, and so we have also a physical measure of the aforementioned particles excited over the thermal vacuum. This happens only in the case of analogue-gravity models, where a nontrivial dispersion appears. Furthermore, we can also interpret the fictitious part of the Hilbert space involving the Hawking partner as associated with the m = −1 branch of the analytic function we discussed in the previous section.
We also think that the above picture would be of strong interest also for theoreticians of TFD, as the particles of the thermal vacuum appear in the present case as real particles, which can actually be measured, as discussed above (e.g., in the case of a white hole configuration of a nonlinear dielectric medium, one may detect photons corresponding to the p, n states in a suitable photon detector). See also [10].
One may wonder if it is possible to extend our considerations also to the case where Γ ω < 1 because of a depletion process of the actually created particles. This means that σ uv = 0, and then Γ ω = 1 − |σ uv | 2 < 1.
As discussed in [1], one maintains the relation (33), and then one obtains A very naive ansatz would be the following: the picture concerning states u, l and states p, n proceeds as above, with the difference that the depletion of Hawking particles (u-modes) at infinity because of the backscattering (v-modes) is such that and v ω (β h ) = Γ ω exp(βω) − 1 .
Still, one has to remark that the effects of depletion not only can but even must in principle act in a different way in the case of the Hawking mode and in the case of the Hawking partner, as Equations (27) and (28) show, disrupting the standard TFD picture and maybe leaving room for a generalization of the standard formalism. We do not delve into this problem herein.

Conclusions
We took into account the general structure of the Hawking scattering process in analogue gravity with weak dispersion. We discussed in detail the process involving the Hawking partner, and we showed that, in general, one could expect different creation rates for the Hawking particle and the (would-be) Hawking partner, unless a weak dispersion takes place, as in the framework of the generalized Orr-Sommerfeld equation we discussed in [1,2], encompassing the framework introduced in [4,8], at least at the leading order in the weak expansion parameter . Including higher-order correction in in general may be tricky. Therefore, in order to show which kind of general form the further corrections one must expect at the subleading order for the S-matrix operator, we proposed a general method based on a rigorous Iwasawa decomposition of the group U(1, −1, 1) in order to prescribe the general form of the S-matrix. In principle, this method allows to describe any possible modification of the S-matrix, perturbative or not, and it makes easily controllable any expansion with respect to an infinitesimal parameter , once the latter is related to the physical perturbation parameter . We did not try to propose here a direct relation between and , since in general it is made tricky by technical subtleties; quantum field theory does not allow to just make naive identifications. A direct perturbative expansion in requires in general further assumptions (as a reasonable smallness of some non-diagonal entries of the S-matrix, see [21]). We leave this delicate question to a future work devoted to a full investigation of the potentiality of the Iwasawa parameterization approach.
We think that the analysis we presented here, together with the ones in [1,2], provides the possibility to improve our knowledge of the analog Hawking effect on the theoretical side, but also on the experimental one, as it is made possible in line with principles a moreaccurate black hole spectroscopy (albeit probably only in a semi-analytical way for what concerns actual calculations of higher-order corrections to the S-matrix). Furthermore, we took into account the thermofield dynamics framework, and we have showed that, at least at the leading order, it maintains its performability also in the dispersive analoguegravity framework, with the (apparently unnoticed before) further bonus that particles on the thermal vacuum can be identified with the modes p and n, and then it becomes experimentally measurable, as discussed in the previous section. The TFD framework is instead jeopardized by the presence of grey-body factors for the Hawking particle and the would-be Hawking partner at the leading order, as different rates no more allow to define a thermal vacuum in the standard way.