Novel triple barrier potential for axial gravitational perturbations of a family of Lorentzian wormholes

We study the behavior of a specific Lorentzian wormhole family under gravitational perturbations. In earlier work (EPJC 80:850, 2020), we have proved the stability of a test scalar field in the background of the wormhole family, where the effective potential was that of a double barrier. Continuing with the stability analysis, here we focus on the more physically relevant scenario, that of axial gravitational perturbations. Interestingly, we find that the effective potential is a triple barrier for lower angular momentum modes. This raises important questions on the ringdown of the corresponding wormhole geometry as well as the gravitational wave echo profile that we try to answer through our work. We study in detail how the geometry of each member wormhole affects the quasinormal modes, the time evolution of the signal as well as echoes which are, in general, very feeble in comparison to the main signal. Different ‘cleaning’ techniques have been used to obtain the echo profile in the time evolution of the signal. Lastly, we dwell on the possibility of our wormhole family as a candidate black hole mimicker, as long as its stability is proven under all kinds of perturbations. We briefly present a comparison of the ringdown characteristics of these wormholes with that of a black hole, in support of this speculation.


Introduction
Wormhole spacetimes have gripped the imagination of both the scientific community and the general public for decades ever since the term was coined by Wheeler [1,2] to describe geometries that act as bridges between two universes. They were first discussed as early as 1916 in the pioneering work of Flamm [3] and then extended by Einstein and Rosen [4]. Their use as passageways for interstellar travel made wormholes a staple for science fiction. The physics community simultaneously took up the challenge of constructing a e-mail: poulamiphysics@iitkgp.ac.in (corresponding author) an actual wormhole geometry and probing the possibility of their existence in nature. It turned out, in the seminal work of Misner and Thorne [5,6], that such spacetimes need exotic matter to remain stable and traversable within the confines of General Relativity [7][8][9] . If one is willing to explore alternate theories of gravity, then we get many scenarios where indeed wormholes can exist without the violation of energy conditions [10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26].
Another hindrance one encounters while investigating wormholes is the issue of radial instability observed in some models. There have been number of stability studies on wormholes supported by scalar fields with the simplest being the Ellis-Bronnikov geometry. Unfortunately, these models have been found to possess a growing mode that perturbs the throat to rapidly collapse or expand within a finite time thus stopping any form of passage through it [27][28][29][30]. Specifically, in [27], the dynamical instability of the Ellis-Bronnikov wormhole, which arises as a special case (n = 2) in our wormhole family, has been studied. However, the question of stability of the other wormholes in the family considered here (i.e. those with n = 2) is yet to be addressed. These results have motivated research in finding the possibility of stable wormhole configurations through different avenues. Some examples include using alternate theories of gravity [18,31,32], introducing certain number of Klein-Gordon scalar fields [33], considering source fluid with specific equation of state [34,35], using specific kind of exotic matter with thin shell [36] -all providing interesting but very specific examples of wormhole geometries. Thus a robust wormhole model stable under all perturbation scenarios seems yet to be discovered. This highlights the significance and need for exploring the stability of various wormhole models under different perturbation schemes.
In spite of these issues plaguing wormholes, they are taken as serious contenders in black hole mimicker analysis. Recently, after the detection of gravitational waves [37,38], wormholes are being considered as a possible black hole mimicker candidate along with a host of other exotic compact objects (ECO) [39][40][41][42]. The huge array of ECOs which are horizonless include gravastars [43,44], boson stars [45][46][47], quasiblack holes [41,48,49], fuzzballs [50,51] and also wormholes [39,52,53]. For wormholes, there has been numerous studies where different aspects of such spacetimes have been looked at for their ability to mimic the behavior of a black hole [39,[54][55][56][57]. We will be particularly interested in the ringdown behavior of the wormholes involving their quasi-normal mode spectrum [39,53,[58][59][60]. The characteristic complex frequencies, dominating the ringdown stage, depend only on the source's parameters which makes them an excellent probe to deduce the stability and nature of the associated object. The ringdown spectrum can also be used as a tool to distinguish the wormholes from black holes which we will discuss in our work as well. Apart from stability analysis, the QNMs associated with a wormhole geometry has been used to determine its shape [59] and have been connected to the shadow radius produced by wormhole [61].
A further distinguishing feature associated with ECOs are the presence of echoes in the ringdown signal [62][63][64][65][66][67]. There have been multiple studies on the formation of echoes in various spacetimes especially for wormholes and wormholeblackhole transitions [68,69]. Even black holes can have echoes too due to presence of additional structure near the horizon, hence making echoes a probe for the strong gravity regime [70][71][72][73][74][75]. Therefore, detecting echoes in gravitational wave observations will always be associated with new physics either at the near-horizon region or near a compact object's surface. In fact recently, many studies have focused on finding traces of these echoes from the LIGO as well as Virgo data [76][77][78][79]. We will encounter these echoes in our work as well. Other detection channels for wormholes include gravitational lensing analysis [80,81] and also looking for orbital perturbations caused by objects on other side of the wormhole throat, as formulated recently in [82,83].
With these motivations in mind, we study the behavior of a two parameter family of wormhole spacetimes under gravitational perturbations. Investigating the relationship between the QNMs and the wormhole's 'shape' (geometry) is the prime objective of this work. A similar analysis of the wormhole family was done for scalar perturbations in our earlier work [84] which also included a detailed discussion of its matter content and energy conditions satisfied. In our present work, we perturb the metric itself, which can be split into axial and polar components [85,86] and is physically more relevant in the context of gravitational wave detection. We focus only on the axial part because of its simplicity as well as the fact that matter remains unperturbed for our wormhole geometries. In the absence of an exact form of wellestablished matter sourcing the wormhole geometries, we use Einstein's equations to get the form of effective energymomentum tensor which is interpreted as an anisotropic fluid. Following [87], we then find the matter sector remaining unperturbed indicated by a vanishing perturbation of the energy-momentum tensor under axial metric perturbations. Recently, there have been attempts to formulate matter sources that generate our wormhole family. For example, wormhole solutions can be realized in the context of asymptotically safe gravity without any need for exotic matter, and these solutions are part of our wormhole family as shown in [88]. Our wormhole solutions can also be realized in f (R) extended gravity where anisotropic dark matter sources them [89]. Both methods have found results that further motivate studies to investigate the nature of matter sourcing our wormhole family conclusively. For now continuing with the stability analysis, we calculate the corresponding quasi-normal modes associated with each member of the wormhole family. Such analysis will contribute towards developing proper templates for possible detection of wormholes through future gravitational wave observation.
A novel feature that emerged during our studies is the triple barrier effective potential associated with the perturbation equation for our wormholes. To this author's knowledge, such a scenario has not been observed in this context. This is true only for the lower angular momentum modes. For higher modes though, the potential becomes a single barrier. It is intriguing to study the echo profiles generated by the triple barrier because they are feeble and not observed directly in the time evolution spectrum. Special 'cleaning' methods [84,90] have to be employed in order to visualise the echoes generated.
The paper is organised as follows. In Sect. 2, we review the essential characteristics of our wormhole family. An indepth analysis of these features has been done in our previous work. We then move on with the stability analysis under axial perturbations in Sect. 3. The dependence of the triple barrier effective potential on different parameter values of the metric is discussed as well. Then, in Sect. 4, we calculate the quasinormal modes and gravitational wave echoes as produced by the effective potential barrier. Finally, in Sect. 5 we study the possibility of this wormhole family being a viable black hole mimicker candidate by comparing the ringdown profiles with that of a black hole. We end with a brief discussion of our observations in Sect. 6.

Review of the wormhole family and its properties
We will begin with a brief overview of the wormhole spacetime we are interested in. The well-known reflectionsymmetric Ellis-Bronnikov wormhole, that was proposed independently by Ellis [91] and Bronnikov [92] in their 1973 papers, is of the form where b 0 denotes the throat radius of the 'drainhole'. Such a spacetime is geodesically complete, spherically symmetric, horizonless and is formed by a massless scalar field with a negative kinetic energy term, thereby making the matter threading the wormhole 'exotic'. This is one of the simplest example of an ultra-static form of Morris-Thorne wormhole [5]. A similar wormhole spacetime but with a more general form of the shape function was suggested in [93] where the authors studied the transmission resonances associated with a line element where b(r ), the shape function, is given by The form of the shape function can be attributed to the relation assumed between the tortoise coordinate and the radial coordinate r , Using one can rewrite the metric in Eq. (2) in a simpler form as Thus, we end up with a geometry that denotes a two parameter family of wormholes where the parameter n is assumed to take only even values to ensure the smooth behavior of r ( ) over the entire region of (−∞ < < ∞). The geometry (shape) of the wormhole is controlled by n while the other parameter is simply the throat radius b 0 . Substituting n = 2 in the above metric gives back the Ellis-Bronnikov geometry while n > 2 values correspond to distinct new wormholes that can be visualized through their respective embedding diagrams ( Fig. 1). In our earlier work, which from now on, we will refer to as Paper-I [84], we had given an in-depth analysis of this wormhole family including its geometry, matter, embedding diagrams, the behavior of r ( ) and the scalar quasi-normal modes associated with each member wormhole. For completeness we will review here, some of the salient features of the n > 2 wormhole geometries in this work as well.
While we know the n = 2 wormhole i.e. the Ellis-Bronnikov case is supported by a massless phantom scalar field, we are yet to understand fully the nature of matter threading the n > 2 geometries. It is known that in GR we need exotic matter, atleast at the throat, to support a traversable wormhole [8]. This necessary violation of energy conditions become a major drawback for the wormholes to exist within the confines of GR. Even so, Visser argues Fig. 1 Embedding diagram of a t = constant, θ = π/2 spatial slice in 3-D Euclidean space for different values of n and b 0 = 1. The cylindrical coordinates z( ) and ζ( ) are plotted with being the parameter that violation of WEC does not necessarily imply the nonexistence of wormholes in nature but simply a possibility of new physics that needs to be explored [94][95][96]. In our wormhole family, all n > 2 geometries show violation of WEC not at the throat but a small distance away from it.
This behaviour can be understood by writing the energymomentum tensor components as a sum of contribution from the phantom scalar field and an additional matter term as shown in Eqs. (10), (11) and (12) of Paper-I. The nature of this extra matter term is not yet known but it was found to satisfy averaged null energy condition (ANEC). As a consequence, the flaring out of the geometry for n > 2 wormholes happens a bit away from the throat. It can be contrasted with that of the n = 2 wormhole where the flaring happens just from the throat as shown in Fig. 1. The n > 2 geometries can thus be called as 'long-necked' wormholes.
We also discussed the quasi-normal modes and echoes associated with each member of the family under the propagation of a massless scalar field. The stability of all wormholes of this family were proven under scalar wave propagation as observed in the damped signal of time domain evolution. After confirming its stability for the scalar case, we are now in a position to continue the stability analysis of the wormhole family under the more physically relevant gravitational perturbations. tion gets complicated in presence of matter which is generally the case for wormhole spacetimes. Recently in [98], the authors have calculated the radial master equation for a general wormhole spacetime using the Newman-Penrose method and the Regge-Wheeler gauge. This analysis is applicable for spacetimes where the matter remains unaffected under axial perturbations thus excluding spacetimes that are solutions to Einstein-Maxwell equations. We will see that the final form of the radial equation (Eq. (32)) and effective potential (Eq. (33)) in [98] matches exactly with our calculation, hence verifying the validity of our analysis.
Going back to our channel of analysis, we have the background metric of Eq. (2) which is spherically symmetric and static but in the perturbed state it might not remain so. Following Chandrasekhar [85], we use a non-stationary but axisymmetric metric to denote the perturbed state of the wormholes. The line element denoting such a geometry is of the form, The components of the metric are functions of (t, r, θ) but not φ, due to axial symmetry. Also, as is evident, since the metric components depend on t, the perturbed object will evolve with time as it should. In the unperturbed state, it reduces to the background wormhole geometry, with the metric functions given as, We notice that a perturbed state of our wormhole geometry will correspond to non-zero values of the quantities σ, q r and q θ while ν, μ r , μ θ , ψ will have small increments δν, δμ r , δμ θ , δψ. The non-zero value of the first set of metric quantities lead to presence of cross-terms which are of odd parity. This corresponds to a dragging of the frame of the perturbed object thus denoting rotation and is called axial perturbation. We will see that for our wormhole family the matter sector remains unaffected (shown in appendix). On the other hand, the small increments δν, δμ r , δμ θ , δψ do not impart any rotational effects and correspond to even parity. The matter content of the spacetime is indeed affected by such polar perturbation modes and while dealing with them, perturbations of the matter fields need to be taken into consideration. In this work, we will focus on the axial perturbation modes mainly because of their simplicity. The perturbation equations for any general Morris-Thorne wormhole geometry has been derived in multiple works in literature [30,99,100]. We will not repeat the entire calculations here but mention the important equations and results. For the axial case, we havē where the Ricci tensorR μν corresponds to that of the perturbed metric shown in Eq. (6). A detailed calculation of the derivation of Eq. (9) is shown in the appendix. The perturbation equation mentioned above holds true only for the components relevant for axial case as in general the spacetime has matter and the RHS will not be 0. Considering the r φ and θφ components of Eq. (9) and simplifying them using the ansatz Q(t, r, θ) = Q r θ sin 3 θ r r 2 − b(r )r with e −iωt as time dependence we get where Q r θ = q r , θ −q θ , r . Thus, the quantity Q r θ involves the derivatives of the metric elements q r and q θ . In order to be able to separate the above equation into its radial and angular parts, The parameter m denotes the angular momentum mode arising due to the separation of variables. We can make a comment regarding the allowed values of m by remembering the relation between the Gegenbauer polynomial and the Legendre polynomial P m (θ ) One notices that for m = 0, 1 the Gegenbauer polynomial is identically 0. So, only for m ≥ 2 one gets non-zero values and as the above equation must hold for all values of angular momentum modes, the radial master equation for axial perturbation will be Thus we have finally arrived at the master equation for the radial component of axial perturbation for our wormhole family. We can further simplify this by implementing the relation R(r ) = r Z(r ), and using the tortoise coordinate we get the familiar Schrödinger-like form of the radial equation, from which we can read off the effective potential [99] 3.1 Effective potential: single/triple barriers The effective potential for our wormhole family can be obtained by substituting b(r ) of Eq. (3) into Eq. (14) so that we get the potential as (15) while in terms of the tortoise coordinates it becomes We can reformulate the potential also in terms of dimension- Once the form of the effective potential has been obtained, it can be plotted for different geometries and angular momentum modes for observing its behavior. We can study the dependence of V(x) on x or alternatively we can take the throat radius as b 0 = 1 and plot V ( ) as a function of .
From the figures of the potential for different parameter values we can make certain observations which are summarised in the following three points.
• In Fig. 2 we find the potential to be a single barrier for wormholes corresponding to both n = 2 and n = 4 irrespective of the angular momentum mode. While the Ellis-Bronnikov wormhole always had a single barrier even for the scalar case (see Paper-I), the n = 4 geometry, unlike its behavior under axial perturbation, had a double barrier potential just like its sister wormholes under scalar wave propagation. • Moving on to n > 4 geometries, we observe the effective potential to be characterized by triple barriers for lower angular momentum modes as can be seen in Fig. 3a. To the knowledge of this author, this is a novel behavior unique to our wormhole family. The height of the two symmetric peaks about the throat increases with increasing n.
The effective potential eventually becomes a single barrier for all geometries as we go to higher m values. The two symmetric peaks about the central one merge to form a strict single barrier as shown in Fig. 3b. The exact value of 'm' from which the potential becomes a single barrier goes on increasing as we go to higher 'n' geometries. In contrast, for the scalar wave propagation case, it apparently seemed that the potential was a single barrier for higher m values but in reality, on zooming in, the double barrier nature was found to be still prevalent as discussed in Paper-I. • Finally, going to even higher values of n (see Fig. 4), we find the triple barrier to be still present for lower m modes but the depth of the well has now increased and is negative. It is known that negative potential wells may harbour bound states and hence have instabilities [101]. So extra caution is required while studying the stability of the higher n wormhole geometries. Also, as mentioned above, with increasing n, the heights of the symmetric peaks about the throat increases. We find that for n ≥ 26 wormholes, the two peaks become higher than the central one at the throat. Hence we choose n = 26 for the plot shown in Fig. 4. For higher modes, the height of the two peaks gradually decreases, later merging into a single barrier.
In an attempt to derive more information about the extrema of the potential, especially for lower angular momentum modes, we again go back to dimensionless variable x = b 0 so that the points of maxima and minima for V ( ) will be a We see that x = = 0 is always a point of extrema irrespective of geometry and m value. For m = 2, 3, i.e. lower angular momentum modes, we can check whether the bracketed quantity gives any real roots for n = 4 as they will correspond to points of extrema other than = 0. We find that for n = 4, m = 2 the quantity in bracket has only imaginary roots indicating that there are no points of extrema. Therefore, n = 4 wormhole always has a single barrier similar to n = 2 case. If one calculates for n = 6, there will be two solutions other than = 0 indicating a maxima and a minima and hence a triple barrier for m = 2. Unfortunately we cannot solve the equation for any general m and n value and hence can comment on each geometry (n value) separately.

Stability analysis: quasi-normal modes and gravitational wave echoes
The ringdown of a perturbed object is a well studied phenomenon using linear perturbation theory. Due to the recent detections of gravitational waves there has been a huge interest in the study of the ringdown profile, which is dominated by the quasi-normal modes, for various objects like black holes, ECOs as well as neutron stars. Quasi-normal mode (QNM) frequencies are complex characteristic frequencies through which a perturbed object looses its energy and settles into an equilibrium state [102][103][104]. These modes are very significant because of their dependence only on the parameters of the final object and not on the cause producing them. As per our assumption, we have the time dependence as e −iωt with ω = ω r + iω i where a negative ω i will indicate damping and hence stability of the system over time. The QNMs apart from proving the stability of an object can also be used as a test for GR and other theories of gravity [105,106].

Numerical computation of QNMs
We will calculate the QNMs for our wormhole family by solving the radial differential equation Eq. (13) numerically. The boundary conditions for a wormhole are similar to that of a black hole, but with outgoing waves at spatial infinity and at the throat. We use the direct integration method and the Prony extraction technique to obtain the dominant fundamental QNM frequencies from the time domain profiles. In direct integration, the radial equation is integrated numerically by imposing proper boundary conditions at the throat and at spatial infinity. Since our wormhole has reflection symmetry about the throat we can separate the solutions into even and odd cases. The method is discussed in detail in [84,107]. On the other hand, for implementing the Prony extraction technique we need to obtain the time evo-  ψ(0, 0). For details on the discretization scheme, see [108] lution plots for different wormhole geometries. We begin by recasting the master radial wave equation (keeping the time dependence) as with ψ denoting the perturbation. Writing the above equation in light cone coordinates (du = dt − d and dv = dt + d ) and using a Gaussian pulse as initial condition along the u and v grid lines we numerically integrate to obtain the timedomain profiles as shown in Fig. 5. The damped ringing in time, exhibiting the decay of the axial gravitational perturbation is clearly visible in the plot. Once we obtain it, we can extract the most dominant frequency by fitting an exponentially damped signal to it. Reference [108] summarises the process in great detail along with the discretization scheme used for generating the plot. Using the numerical methods discussed above we calculate the QNM frequencies associated with different wormhole geometries and for different angular momentum modes. All computations are done using Mathematica 12. The Tables 1, 2 and 3 lists QNM frequencies calculated for n = 2, 4 and n = 6. The values calculated from both methods show better matching for lower angular momentum modes. As m value increases, DI becomes less stable for wide range of matching points. So, for further analysis we will prefer the results as obtained using the Prony method.
We now study how the geometry of the wormholes influence the QNM values i.e. how the QNMs depend on parameter n.
• The behavior of the real part of the QNM frequency for different modes and geometries are similar. As m increases, the magnitude of the QNM increases indicating a higher value of frequency. This is because of the fact that with increasing m, height of the potential barrier also increases and so only waves with enough energy i.e. high frequency are able to cross the barrier. • The behavior of the imaginary part is rather different for n = 2. As m increases, the magnitude of ω i increases, indicating that m = 2 is the least damped mode. For n > 2 geometries, the value of ω i first increases and then decreases with m. So ω i for m = 3 generally has higher magnitude than all other m. Also we cannot predict the most dominant mode for n > 2 geometries as the least damping time will correspond to m → ∞. The Fig. 7 shows the effect of the long-lived modes resulting in the slow damping of the signal over time. The red curve corresponds to m = 20 which clearly decays very slowly than the m = 8 mode for n = 10 wormhole geometry.  A similar result is reported in [109] for a completely different wormhole scenario. Hence, the above analysis highlights the fact that we can easily distinguish the wormholes of different 'shapes' just from their fundamental axial quasinormal mode frequency which was one of the goals of this work.

Are higher n wormhole geometries unstable?
As we have seen in the potential plots, the minima for some cases reach negative values (large n, small m) (see Fig. 4), which might harbour bound states and may have instabilities. In all the previous time domain plots, we have used a grid spacing of h = 0.1 with the convergence of the integration scheme being checked for smaller h values. But for high n geometries choosing h = 0.1 shows instability in the time domain profile which is a numerical artifact as they disappear and we get damped signal for smaller h values like 0.02 or 0.01. This happens because the well width is so narrow that the variation in the potential cannot be detected during the numerical analysis. Hence we need to reduce the grid spacing sufficiently so that it can trace the variations of the potential near the wells. No instability has been observed for any parameter values in the time evolution profiles of our

Observing echoes in different wormhole geometries
In this section we will have a look at the echo profiles in the time domain signal generated by the triple potential barrier. Whenever a potential is characterized by a second peak or a reflecting surface, the gravitational wave signal not only decays over time but gets reflected. This leads to repetitive bursts of signal, after the initial ringdown, with decreasing amplitudes called 'echoes'. Under scalar wave propagation our wormhole geometries had a double barrier and hence had distinct echo signatures in their time domain profiles. These echo profiles became more prominent for larger n geometries because of the distinct reflection occurring from the sharp peaks in the potential. The separation between two echo signals is equal to the time taken by the signal to travel to-and-fro between the two peaks.
We expect to see similar echo signatures in the case of axial perturbation as well because here too we have triple potential barriers. It is needless to say that the generation of the echo signal for this case will be much more rich and complicated because the signal has to get reflected between multiple barriers. But the first question that we might ask is where are the echoes? In Fig. 5 we see that even for m = 2 case there are no echoes. Even if the time domain profile is calculated for longer time, echoes are not visible. One might think that the situation is similar to the scalar case (as discussed in Paper-I) where low n wormholes did not show echoes due to the wide peaks causing back scattering and hence going to large n solved the issue as the potentials had sharp and narrow peaks. But it is not so simple in the case of axial perturbations. Here too the back scattering from the wide tails of the triple barrier suppresses the echoes. To add to the problem, the separation between two consecutive peaks is very small and there is very less chance of observing an echo. The situation does not get better even if we go to higher n geometries because for axial perturbation, with increasing n, the depth of the well (minima) increases but the height of the peaks remain almost unchanged and so does the width of the barriers.
One possible way out of this scenario is to 'clean' the profile. Such a technique was used by us in Paper-I to observe echoes for small n geometries following [90]. The main essence of the process was to remove, from the full spectrum, the effect of the scattering of signal from the single barrier (see Fig. 9). After subtraction, one should ideally be left with only the scattering happening due to the double barrier i.e. echoes. But implementing such a method for the axial perturbation case becomes quite tricky. We first need to ascertain the effect of which potential peak needs to be subtracted from the full spectrum since now we have three distinct single barriers. We begin by taking a look at the three possible scenarios at hand for a specific wormhole geometry, say n = 10, m = 2, We take the single barrier to be the right most peak among the three barriers and so the potential is that of the wormhole from ≥ 0.955 which is the point of minima. For region < 0.955, the potential is kept constant and equal to the value of the potential at = 0.955 i.e. 1.6. Similarly, the double barrier potential is considered by keeping the potential peaks at = 0 and the one to its right so that the potential is that of the wormhole for ≥ −0.955 while for < −0.955 it is again constant with value 1.6. All these values correspond to the n = 10, m = 2 geometry (Fig. 9). The well potential simply denotes the first well a signal encounters on the positive side of the throat so that for ≥ 0 we have the wormhole potential but on the other side the potential becomes constant and equal to the value at throat. Before we move on with the subtraction procedure let us see how the spectrum for each of these potential cases behave.
To make sure that after cleaning if we observe the echo, that is not the effect of solely the potentials that are subtracted, we plot the scattering from these potentials alone. It can be seen from Fig. 10 that none of the potential scenarios possess echoes but just a damped ringdown signal. We now subtract the effect of scattering from each of these potentials from the full spectrum. The 'cleaned' profile corresponding to each of the above cases is shown in Fig. 11. The echo signal can now clearly be seen after cleaning for all three cases. In the third figure (Fig. 11c), we subtract the effect of the first potential minima or the well and observe an echo. The difference of this signal with subtracting the double barrier potential is that we can observe the scattering effect occurring from the tail of the central peak on the left. Finally, this procedure helps us clearly visualize the effect of the triple barrier on the echo signal. For such a potential the signal will not only be reflected between the first two peaks but some part will be transmitted into the other potential gap. There again it will undergo multiple reflection and some part will be transmitted to infinity while some will again get transferred to the other potential gap. In this way the echo signal for the triple barrier is a superposition of multiple echoes produced via reflections between the two potential gaps.
From the figures showing echo profiles, we see that the height of the signal increases with time and then decreases, indicating the profile to be an 'echo'. Also, there is no initial ringdown because it has been subtracted and we are left only with the effect of reflections from the potential barriers.
We now continue with some interesting observations which can be made from the echo profiles.
• Which peak, out of the three, should be chosen as the single barrier?
The single barrier that we used for 'cleaning' the profile and observe echoes is taken as the first peak of the potential on the right of = 0 that a signal from infinity will encounter first. We have also checked with the cen- The echo profile obtained by removing the effect of the single barrier for n = 6 (blue curve) and n = 10 (red curve) geometries with m = 2 tral peak but the echoes are not distinct after cleaning. This may be due to the fact that we are taking constant potentials on both sides of the peak which changes the behavior of the potential at both infinities.

• Why is only one echo signal observed?
It is very difficult to observe multiple echo packets because the signal damps very rapidly and the amplitude of the signal is too low to observe echoes even with cleaning procedure. Thus all the reflections between the different potential peaks superimpose to form the dominant echo signal with significant amplitude, which we observe.

• Comparison of signal from different cleaning procedures
In Fig. 12 we observe the echo profiles obtained by subtracting different sections of the potential for n = 10 wormhole. The blue curve denotes the echo observed while subtracting just the effect of the first single peak. Hence the signal starts at an earlier time as well. Now, when we subtract the effect of the well potential we get the red curve which denotes that the echo starts a bit late and is of weaker strength. Finally, we get the green curve corresponding to the echo obtained by subtracting the effect of the double barrier that starts even later and is more weak than the other two. • Why echo amplitude decreases with increasing n ?
With increasing n we find the echo amplitude to be decreasing as seen in Fig. 13. This occurs because of the increasing depth of the well for higher geometries which makes it more difficult for a trapped wave to escape from the well. So even though the peak height is small and the potential peaks have wide tails for small n, we still get higher echo amplitude after 'cleaning' because of the small well depth. • How does echo signal change with increasing m ? As m increases for a particular n, the symmetric peaks on both sides of = 0 will be reducing in height and finally merging with the central peak to give a single barrier for high m. As a result the echo amplitude also decreases. Figure 14 shows how small the echo signal is for m = 3 case in comparison to the m = 2 when we just subtract the single barrier. These behaviours are as expected because of the small height of potential peaks.
From all the above observations regarding the echo structures we can come to the conclusion that the time evolution of a signal in the background of our wormhole geometry is much more diverse and complex because of the triple potential barrier.

Possible black hole mimicker?
One of the main motivations behind the study of our wormhole geometry is to develop the wormhole family as a potential black hole mimicker model. After the detection of gravitational waves, interest in studying observable features of black holes have increased tremendously -in particular, with reference to post-merger ringdown behavior via quasinormal modes. Parallely, one notices an increase in studies on black hole mimickers which include a variety of compact objects sans the event horizon. These objects are expected to mimic some features of a traditional black hole and thus question the very existence and detection of the black hole through observations. It is therefore of utmost importance to study all possible mimicker models, not only to verify the existence of black holes without doubt but in this process also investigate any new physics that may arise. Having said that, any spacetime to exist in nature and mimic a black hole signature must pass through all possible stability analysis tests. The instability of some wormhole geometries under radial perturbation is well studied in literature [27][28][29][30]. So apart from the scalar wave propagation [84] and axial perturbation already performed that confirm the stability of our wormhole family, other perturbation tests need to be explored in future studies to ensure the viability of the mimicker model.
Wormholes are one of the contenders which have been studied as a mimicker in many works [39,42,53,110]. Through our wormhole family metric we can simultaneously study the possibility of all the member wormholes being mimickers. We focus on the quasinormal ringing of the wormholes and try to find parameter values for which the ringing of the wormholes will be closest to that of a black hole in a similar vein, as done in [53]. For the stability analysis, we had kept the throat radius parameter b 0 as unity. Now we will tweak it as well as 'n'. We take the case of a Schwarzschild black hole with 2M = 1 so that for m = 3 the fundamental QNM frequency is 1.19882 − i0.185822. We aim to find wormhole geometries in our family that might have similar QNM frequencies. Following are some scenarios where we find the closest mimicking behavior of the wormholes where the listed fundamental QNM values have been calculated using the Prony extraction scheme, For the above mentioned wormhole geometries with specific throat radius, the dominant damping rate is very similar to that of the black hole but the corresponding oscillation frequency i.e. the real part of the QNM frequency is quite different. So the frequency of the QNM ringing can be used to distinguish any black hole mimicker from an actual black hole provided we observe the m = 3 mode. As value of n increases, it becomes easier to distinguish the wormholes from a black hole as the real part of the QNM increases. Figure 15 shows the time domain signal for the above parameter values of the wormhole geometries in comparison to the behavior of the black hole. The damping rate is similar but the ringing structure determined by the real part of QNM frequency is quite different. The above mimicking properties of the wormholes hold true for the m = 3 angular momentum mode. For mode m = 2, the period over which QNM ringing occurs in the time domain evolution plots for Schwarzschild BH is much less than that of wormholes even though the dominant QNM frequencies can be obtained to be similar for proper choice of throat radius (see Fig. 16). Thus, while for Schwarzschild black hole with 2M = 1 we have 0.747−i0.1779 as the dominant QNM frequency; for wormholes it becomes 0.976232− i0.171425 corresponding to n = 2, m = 2, b 0 = 1.78 and for n = 4, m = 2, b 0 = 1.54, it is 1.20934 − i0.174385.
From the above analysis we come to the conclusion that our wormhole geometries can mimic only the damping rate of the ringdown profile of a black hole under certain parameter choices. The real part of the QNM will be quite distinct for the wormhole giving away their identity. On the other hand, higher n wormhole geometries are poor mimickers and can easily be distinguished from their axial quasi-normal mode frequencies.

Discussion
Our work comprises of a complete analysis of the behaviour of a two parameter Lorentzian wormhole family under axial gravitational perturbation. Following Chandrasekhar we calculated the master radial equation where an axially symmetric non-stationary spacetime represents the perturbed metric. The first significant result of this study is the observation of a triple potential barrier corresponding to lower angular momentum modes in n > 4 geometries. As the mode number increases, the peaks come closer and finally merge into a single barrier. For n = 2 and 4 wormholes, the potential is always a single barrier for all angular momentum modes. This is followed by the calculation of QNM frequencies for different wormhole geometries. The signal frequency was increasing as we went to higher angular momentum modes for all values of n. However, the behavior of the imaginary part of the QNMs that controls the damping rate is rather interesting. We observe long-lived modes which correspond to low damping rates as we go to higher angular momentum values in the n > 4 wormholes. For n = 2 case the imaginary part of QNM increases with increasing m value making m = 2 the lowest damped mode. Such long-lived modes might lead to insta-  [111] which needs to be explored in future. These modes can arise in different wormhole spacetimes as well which have been studied extensively in the literature [112,113]. One of our aim of distinguishing the 'shape' (geometry) of the wormholes just from their fundamental QNM frequency is also achieved as each geometry has a distinct frequency evident from the plots shown in Fig. 6 corresponding to a particular throat radius. The multi-peak nature of the potential for our wormhole family makes it appropriate for studying the generation of gravitational wave echoes. The formation of echoes is much more involved and rich for a triple barrier potential as compared to a double barrier. The absence of sharpness in the potential peaks and the peaks being very closely spaced makes the echoes feeble and hard to see directly in the time domain profile. As a remedy, we try to 'clean' the spectrum of the effects of scattering of signal from each barrier peak. First, we remove the effect due to the single barrier, then move on to the double barrier and finally remove the effect of scattering from the first well. This subtraction procedure leaves us with the net scattering off the entire triple barrier a.k.a the Fig. 16 Schwarzschild profile (red), n = 2 (green), n = 4 (blue) with the b 0 values mentioned above echoes. The echo profiles are observed and compared after all the subtraction schemes mentioned above. Even after such a cleaning procedure we are able to observe only a single echo signal because the amplitude of the signal gets damped very quickly making observation of multiple echoes rather impossible. We compare the strength of the echo profile for different parameter values as well as angular momentum modes.
As expected, the strength of the echo signal decreases with an increasing value of m because the potential peaks are slowly merging to form single barrier. Another interesting aspect regarding the stability of the spacetime is the presence of well in the potential whose depth increases as we go to higher n geometry wormholes. The presence of well which may support bound states indicates possible instability. But for all wormhole geometries belonging to our family, we get a damped signal in the time domain evolution hence suggesting stability under linear axial perturbation.
We also briefly study the possibility of our wormhole family being a black hole mimicker by comparing the ringdown signal with that of a black hole. However, before exploring the mimicking features of our wormholes, it must be borne in mind that many wormhole geometries studied in literature do show instability under radial perturbation. Our wormhole family has narrowly escaped the instability scenario under axial perturbation despite the presence of 'potential well' and 'long-lived' modes. Thus, to ensure the robustness of our family, we need a complete stability analysis including radial perturbations, which we intend to pursue and complete in future. Keeping aside the above comments for now and focusing on our present work, we note that some parameter values can yield QNM frequencies that are very close to that of a black hole. However even though the damping rate is nearly identical, the real part of the QNM is distinct and distinguishable from that of the black hole. So even though we get identical damping rates, the frequency would give away the identity of the wormholes when compared to a black hole. The larger n wormhole geometries are even poor mimickers and can easily be distinguished using the axial QNMs. Future works can focus on the polar perturbations of the wormhole family and an even more accurate study would be to include rotation in the metric, since, in nature, all astrophysical objects are known to be rotating.
In summary, we performed a detailed analysis of the QNM spectrum and the echo structure for our wormhole family. We have also speculated on how such wormholes can be of use as black hole mimicker candidates. gesting the problem, carefully reading the manuscript and for his valuable comments in improving it. She also thanks Indian Institute of Technology, Kharagpur, India for support and for allowing her to use all available facilities there.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There is no data involved in the work.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.

Appendix A: Perturbation of energy-momentum tensor
Axial perturbation of the metric is associated with inducing rotation of the perturbed object due to presence of non-zero cross-terms in the metric as shown in Eq. (6). Such a perturbation will affect both the metric as well as the matter content of the spacetime. While deriving the perturbation equation we follow the notation used in [85] with the coordinates being t (= x 0 ), φ(= x 1 ), r (= x 2 ) and θ(= x 3 ). The Einstein's equation for the perturbed metric of Eq. (6) will bē where the total energy-momentum tensor is taken as a sum of the background contribution T μν and a small perturbation to it. The metric of Eq. (6) can be written in tetrad basis as ds 2 = η ab e a e b with η ab being the Minkowskian metric and a, b denoting the tetrad indices. Hence the corresponding tetrads will be of the form By setting q r = q θ = σ = 0 the tetrads for the background unperturbed metric will be obtained. For axial perturbation, we will be interested in φt, φr and φθ components of the perturbation equation which will correspond to 10, 12 and 13 components in frame basis. Thus we get R ab = δT ab (A6) as η 10 = η 12 = η 13 = 0 and T 10 = T 12 = T 13 = 0. To complete the analysis, the perturbation in the energy-momentum tensor needs to be calculated. In case of a known matter field sourcing the geometry, the energy-momentum tensor can be directly perturbed to find the RHS of the perturbation equation. It can be seen that for a spacetime sourced by a scalar field, δT ab = 0 however for solutions of Einstein-Maxwell systems δT ab = 0 like in the Reissner-Nordstrom black hole [85] (see also [30] for example in a wormhole spacetime). To continue a similar analysis in absence of an underlying matter theory sourcing our n > 2 wormhole geometries, we follow [87] and consider the effective energy-momentum tensor as obtained from Einstein's equation with r ( ) = ( n + b n 0 ) 1/n . Note that the Ellis-Bronnikov wormhole corresponding to n = 2 case is supported by a phantom massless scalar field of the form φ( /b 0 ) = √ 2 tan −1 ( /b 0 ). Again, for a phantom scalar field, ρ φ = τ φ = −p φ = −(1/2)φ 2 where the derivative is with respect to . Thus the above equations with r ( ) = 2 + b 2 0 are indeed satisfied by the scalar field giving, These have been discussed in detail in Paper-I along with properties of energy-momentum tensor satisfied by n > 2 wormholes. Now, the energy-momentum tensor, as shown in Eqs. A7, A8 and A9, is interpreted as an anisotropic fluid of the form where g μν corresponds to the background metric and ρ, τ, p are the energy density, radial and tangential pressure respectively. u μ (= u t , 0, 0, 0) and x μ (= 0, x r , 0, 0) are respectively timelike four-velocity and spacelike unit vector orthogonal to u μ . Converting to frame basis and taking small perturbation we arrive at δT ab = δ(ρ + p)u a u b + (ρ + p)δ(u a u b ) It can now be easily visualized that δT 10 = δT 12 = δT 13 = 0 which leads to the final perturbation equation being the one shown in Eq. (9) with the RHS being 0.