Revisiting a family of wormholes: geometry, matter, scalar quasinormal modes and echoes

We revisit a family of ultra-static Lorentzian wormholes which includes Ellis–Bronnikov spacetime as a special case. We first show how the required total matter stress energy (which violates the local energy conditions) may be split into a part due to a phantom scalar and another extra piece (which vanishes for Ellis–Bronnikov) satisfying the Averaged Null Energy Condition (ANEC) along radial null geodesics. Thereafter, we examine the effective potential for scalar wave propagation in a general setting. Conditions on the metric function, for which the effective potential may have double barrier features are written down and illustrated (using this class of wormholes). Subsequently, using numerous methods, we obtain the scalar quasinormal modes (QNMs). We note the behaviour of the QNMs as a function of n (the metric parameter) and b0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_0$$\end{document} (the wormhole throat radius). Thus, the shapes and sizes of the wormholes, governed by the metric parameter n and the throat radius b0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_0$$\end{document} are linked to the variation and the values of the QNMs. Finally, we demonstrate how, for large n, the time domain profiles exhibit, expectedly, the occurence of echoes. In summary, our results suggest that this family of wormholes may indeed be used as a template for further studies on the gravitational wave physics of exotic compact objects.


Introduction
Much of the interest today in traversable Lorentzian wormhole spacetimes (originally proposed in [1][2][3]) revolve around the question: do they exist? While the existence of black holes is no longer in doubt (more so after recent observations in M87 [4]), the wormhole story is far from complete. The same may be said about naked singularities and cosmic censorship too [5,6].
The existence question on wormholes is based on a couple of issues. Firstly, classical General Relativity (GR) along with the imposed energy conditions does not allow wormholes [7][8][9][10][11]. In other words, the shape of the spatial slice of a wormhole is such that a converging null geodesic congruence would have to be defocused, as long as a throat and a 'flare out to the other universe' (second asymptotically flat region) -both necessary geometric features of wormholeshave to be admitted [12][13][14]. To get away with this so-called 'defect' or 'problem' one may appeal to modified theories of gravity. In such theories, the energy conditions on matter may hold but the convergence condition is violated [15,16]. Examples of wormholes in modified gravity are numerous [17][18][19] and they largely have energy-condition-satisfying matter [20][21][22][23][24][25][26][27][28][29][30][31][32][33][34] and sometimes, non-phantom fields [35,36]. In addition to taking refuge in modified gravity, other ways of restricting the violation of energy conditions are known. These include dynamic wormholes [37][38][39][40][41], a proposal for limiting the amount of exotic matter [9] etc.
The second issue concerns possible signatures. Till recently, the most compelling suggested signature for wormholes appeared to be from gravitational lensing in such spacetimes [42][43][44]. However, with the advent of gravitational wave astronomy [45][46][47][48][49][50], one comes across the notion of black hole mimickers [51][52][53][54][55][56][57][58][59][60][61][62] (eg. wormholes, gravastars and other ultra-compact objects) which can ideally mimic the results found using black holes, in GW observations. In other words, one may, in some scenarios, be able to explain GW observations using such black hole mimickers as the end state of black hole and/or neutron star mergers [63,64]. Therefore, to improve upon viable templates for black hole mimickers, it is necessary to study various properties associated with them. One such property is a study of the quasinormal modes which can be used by observers at GW interferometers for verifying a wormhole proposal.
It is therefore of importance to study quasinormal modes of different types of wormhole spacetimes. This has been done to some extent in [55,[65][66][67][68][69][70][71][72]. One of the purposes of this article is to further this line of thought for another class of ultra-static wormholes.
The choice of the family of wormholes we make here is based on an earlier paper [73] where the well-known Ellis-Bronnikov spacetime [74][75][76][77] had been extended to provide a two-parameter family of spacetimes. We first show how this generalisation leads to a matter-stress energy which has (a) an energy condition violating part generated via a phantom scalar field and (b) another piece (of zero value for Ellis-Bronnikov) which satisfies the Averaged Null Energy Condition (ANEC) along radial null geodesics.
The scalar field which generates the wormhole can therefore be perturbed , with the ensuing perturbations satisfying φ = 0. Thus, scalar quasinormal modes exist and are of interest in the context of this family of wormholes. Surely, in future, gravitational perturbations will have to be looked at, in order to make direct contact with GW observations. One may also recall that scalar perturbations do appear to be relevant in modified theories of gravity, where they arise as a so-called extra 'breathing mode', which is absent in GR. Numerous ways of detecting such a mode using a network of detectors have been discussed in the literature [78,79].
Several types of wormhole spacetimes have been investigated recently, which may be thought of as templates for future studies especially with reference to GW observations. For most of these spacetimes, there is no clear understanding of the matter required, from the viewpoint of a Lagrangian based field theory coupled to gravity. One such well-known example is Damour-Solodukhin (DS) spacetime [54,67,70,[80][81][82][83]. Another interesting geometry is the black-bounce metric constructed in [84] where a continuous parameter in the metric allows a transition from a regular BH to a one-way wormhole with an extremal null throat (black bounce) and then to a traversable wormhole. The spacetime we work with here has distinctive features and may be viewed as another such wormhole template.
Our work reported in this article is organised as follows. In Sect. 2, we introduce the wormhole spacetime. We write down the Einstein tensor and equate it to the energymomentum tensor of the 'required matter'. We also show how the required matter can be modeled with a phantom scalar and another extra piece satisfying ANEC. The embedding diagram and the behaviour of the expansion of a geodesic congruence is also briefly discussed. In Sect. 3, we move on to studying scalar wave propagation. We show when the s-wave effective potential can behave as a double barrier by deriving the required conditions on the metric functions. Section 4 is devoted to the scalar quasinormal modes in this class of spacetimes The interesting possibility of the QNMs being used as a tool to determine the shape of a wormhole geometry in this family is presented here too. Further, the occurence of echoes in the time domain profile is shown and analysed in Sect. 5. Finally, in Sect. 6 we end with our conclusions and remarks.

The family of ultra-static spacetimes
In their 1973 papers [74,75], Ellis and Bronnikov, independently constructed a spacetime using a phantom (negative kinetic energy) scalar field source. Their work produced a static, spherically symmetric, geodesically complete, horizonless manifold with a throat (which Ellis called a 'drainhole') connecting two asymptotically flat regions. The line element of the spacetime constructed by Ellis and Bronnikov is given as, where b 0 is the throat radius of the 'drainhole'. This spacetime is known today as the Ellis-Bronnikov wormhole. If the Morris-Thorne conditions [7] necessary for the construction of a Lorentzian wormhole are considered then one finds that there is ample scope for various similar wormhole geometries to exist. Thus, a generalised version of the Ellis-Bronnikov wormhole geometry was suggested in [73] as a two-parameter (n and the throat radius b 0 ) family of Lorentzian wormholes. When the parameter n takes the value n = 2, we get back the Ellis-Bronnikov spacetime. The motivation behind such a construction was to study the geodesics and propagation of scalar fields for a wider class of wormhole spacetimes and note various differences as well as similarities. The work also included the observation of resonances in the transmission coefficient for n > 2 geometries thus indicating that n > 2 geometries are clearly different from the Ellis-Bronnikov geometry (n = 2). The line element of the generalised Ellis-Bronnikov spacetime involving the parameters n and b 0 is given as, The parameter n is allowed to take only even values to ensure the smooth behavior of r ( ) over the entire domain of the 'tortoise' or 'proper radial distance' coordinate (−∞ ≤ ≤ ∞). Note that the functional form of r ( ) has some curious features. At = 0 (wormhole throat), only the nth derivative of r ( ) is non-vanishing. Also, the function V ( ) = r r (V( ) is the effective potential discussed in Sect. 3.1) has a nonzero n-th derivative at = 0, which is negative in value for n = 2, but positive for all n > 2. These facts will be crucial while discussing the effective potential for scalar wave propagation later in this article.
The line element can also be written in an alternative form in terms of the usual radial coordinate r as, where r and are related through

Geometry
Given the metric functions in Eq. (2), we note that the spacetime is spherically symmetric and ultrastatic. All the metric components are independent of time and hence, all t = constant slices are identical. This will be useful when we embed a 2-D slice of the wormhole in flat space, in order to understand its shape. The embedding diagram encodes the shape of the wormhole through a variation of the shape function b(r ) or r ( ) (w.r.t the ' ' coordinate). It is easy to observe that the R 00 component of the Ricci tensor will always be zero for this family of wormholes, irrespective of the form of r ( ) or b(r ). The spacetime geometry owes such features since the metric is ultra-static. We will see the distinction between geometries for different n values explicitly when we plot the embedding diagrams for different n. One of the objectives of this paper is to distinguish the different wormhole geometries for different values of n using the corresponding scalar quasi-normal modes.

Matter, energy conditions
The energy conditions are a way of ensuring that a solution of Einstein's equations is physically viable. The creation and maintenance of any traversable wormhole was first studied by Morris and Thorne [7] through the energy conditions that they are supposed to satisfy. It was observed that for a traversable wormhole to exist in GR the Weak Energy Condition (WEC) must be violated atleast at the throat. This meant that one requires exotic matter (i.e. matter violating the energy conditions). Later studies showed that all classes of static wormholes in GR need exotic matter for stability [40,41,85]. This is a major drawback for wormholes and, as stated before, forbids their existence within the tenets of GR.
Let us first write down the energy momentum tensor for a general r ( ) using the Einstein equations and the Einstein tensor. The energy-momentum tensor defined by its diagonal components in the frame basis, i.e. T 00 = ρ( ), T 11 = τ ( ), T 22 = T 33 = p( ) has the following form with (8π G = c 2 = 1), It is clear from the above expressions that, for any r ( ), ρ + τ = −2 p = −2 r r . We will see later, how the single or double barrier character of the effective potential is related to the energy conditions and the nature of r ( ). In general for any wormhole, r ( ) is preferably an even function, r (0) = b 0 , and r ∼ as → ±∞. Thus r ( ) always has a minimum at = 0.
We now write down the energy-momentum tensor for our wormhole spacetime. Similar to the case of general r ( ), we calculate ρ, τ and p for our wormhole family using x = b 0 , We have written the ρ, τ, p as a sum of (i) a contribution from a phantom scalar field (ρ φ , τ φ , p φ , i.e. the first terms in the R. H.S. of (10), (11), (12)) and (ii) extra matter (ρ e , τ e , p e ), i.e. the second terms in the R. H. S. of (10), (11), (12). Recall that for n = 2, the well-known Ellis-Bronnikov spacetime can be obtained as an exact solution of the Einstein equations with a phantom scalar. The extra matter, which vanishes for n = 2, arises when we wish to consider generalisations for n = 2. One may note that the ρ φ , τ φ , p φ satisfy the relation where φ = dφ d and φ is a phantom scalar with a negative kinetic energy term in the action.
Thus, one may sum up and say that the the total action which leads to the equations of motion of which the n ≥ 2 line elements are solutions, is given as where S g is the standard Einstein-Hilbert gravity action and S ph represents the phantom scalar action (note the plus sign (our signature is -+++) which yields the negative kinetic energy for the phantom). The extra term in the action, S extra , is responsible for giving rise to the n > 2 wormhole geometries and it vanishes for the n = 2 Ellis-Bronnikov wormhole. Its presence is required because the scalar field alone cannot generate the entire set of n > 2 geometries. Though we do not know the exact form of S extra (i.e. in terms of a Lagrangian density), we do know, from the field equations, the energy-momentum tensor (the ρ e , τ e , p e mentioned above) resulting from S extra (see Eqs. (10), (11) and (12)). We will show below how this extra piece of the energy momentum tensor indeed satisfies the ANEC. The equation of motion for the phantom scalar φ is just φ = 0. We can solve for the scalar field for arbitrary n metrics to get If n = 2 one obtains φ(x) = √ 2 arctan x, which is the solution for Ellis-Bronnikov spacetime. The general solution with the hypergeometric function has a behaviour similar to the arctan x solution for n = 2. This feature can be noted in the graphs in Fig. 1, for different n and b 0 = 1. It is easily seen that the Euclidean trace ρ + τ + 2 p = 0 (follows from R 00 = 0, stated earlier). If we check the WEC inequalities: ρ ≥ 0, ρ + τ ≥ 0 and ρ + p ≥ 0, for the ρ, τ and p stated above, it can be seen that the second inequality is always violated for any value of parameters n and b 0 as well as . The other inequalities may not be violated for certain values of the parameters and over restricted domains of the coordinate. This can be observed more clearly when we plot ρ( ), ρ( )+ τ ( ) and ρ( ) + p( ) with respect to .
We find that for a range of we do get matter with positive energy density, as can be seen from Fig. 2a. The third energy condition is also partially satisfied (see Fig. 2c) but as already stated, the second energy condition is always violated (Fig. 2b) except at = 0 (for n > 2). Since the ρ( ) + τ ( ) and ρ( ) + p( ) inequalities are violated the Null Energy Condition (NEC) is violated as well.
It is clear from the expressions for ρ φ , τ φ and p φ that ρ φ ≤ 0, ρ φ +τ φ ≤ 0 and ρ φ + p φ = 0, where the equality in the first two inequalities happens only at the asymptotic infinities. Thus, the phantom scalar field violates the energy conditions everywhere. On the other hand, the additional matter given via ρ e ,τ e and p e has, as we shall see below, the property that it can satisfy the Averaged Null Energy Condition (ANEC) along radial null geodesics.
The ANEC integral, as is well-known is evaluated along null curves with tangent vector k i , and is given by: The ANEC is therefore stated as I ≥ 0. For our line element, it is easy to see that (ṫ,˙ ,θ,φ) = (1, 1, 0, 0) represents null geodesics. If we evaluate the ANEC integral along this set of null geodesics, we need to work out the integral (choosing = λ as the parameter labeling points on the null geodesics): Using the expressions for ρ e and τ e given in (10) and (11), we find that (for even n ≥ 2), This expression is manifestly positive (see Fig. 3) for all even n and has an asymptotic value equal to 4 3 . As is well-known, a phantom scalar is a source for the n = 2 Ellis-Bronnikov geometry. However, for n > 2, the phantom scalar alone cannot generate the geometries. Additional matter is required, which, interestingly satisfies the ANEC. This is actually the reason behind the fact that the violation at the throat does not occur for the n > 2 geometries. The fact that the 'flaring out' does not happen at the throat is related to this additional ANEC satisfying matter. This aspect which is related to the shape of the geometry will become clearer through the embedding discussed below.

Embedding
In order to visualise the shape of such generalised wormholes in 3+1 dimensions, for different values of n, the 2-D spatial slice (t = constant, θ = π/2) of the spacetime is embedded in 3-D Euclidean space with the line element written in cylindrical coordinates. To begin, we consider the t = constant, θ = π/2 slice given as, The metric has now been reduced to a 2-D geometry which can be visualised by embedding it in a 3-D flat space.
Since there is axial symmetry we use cylindrical coordinates (ζ, ψ, z) with the metric of the flat space being As the surface possesses axial symmetry, ψ = φ, z = z( ) and ζ = ζ( ). Thus we get, Comparing this with Eq. (19) we find, The Eq. (23) is numerically integrated with b 0 = 1 and z is plotted with respect to ζ in a parametric plot giving the embedding diagram for different values of parameter n.
As is evident in the plot of Fig. 4, the geometries for different n are quite distinct from each other. It is also observed that with increasing value of n the flaring out of the wormhole right from the throat becomes less prominent. For larger n values the embedding will resemble a uniform tunnel connecting two remote flat regions. One may also name these geometries as long necked wormholes.
The change in shape of the geometry with increasing n is also manifest in the behaviour of the expansion of a null geodesic congruence. It can be shown that the expansion of a null geodesic congruence (k i = (1, 1, 0, 0)) is equal to 2 r r . Evaluating for n = 2 one notices that it is zero at the throat but changes significantly for > 0 and < 0. In contrast, for n > 2, as shown in Fig. 5, there is a domain around = 0 where the expansion is almost zero (exactly so only at = 0). Thus, for large n, a converging null geodesic congruence entering from one universe, becomes almost parallel ( ∼ 0 ) over this extent of , near the throat, before it starts diverging to become zero again in the 'other universe'.
In the subsequent sections we will try to see, to what extent, the nature of the geometries for different n is reflected in the quasinormal modes and in echoes.

Scalar wave propagation and effective potentials
Any open system under perturbation emits radiation and looses its energy through a discrete set of complex frequencies of the form ω Q N M = ω r + iω i called the quasi-normal modes [55,66,68]. The QNMs are associated with specific boundary conditions, i.e. purely outgoing waves at spatial infinities. The real part of ω Q N M determines the oscillation frequency while the imaginary part gives the damping rate of the field over time. The QNMs can be excited due to scalar, vector (electromagnetic) or tensor perturbations of a given geometry. More details on the mathematical definition of QNMs and how to find them appears later in this section.
Let us now consider the propagation of a massless, minimally coupled scalar field in our wormhole spacetimes. This scalar field may be viewed as a perturbation of the phantom scalar which was part of the source for the metric. In other words, we may write φ + as the perturbed phantom scalar and note that will satisfy the same equation = 0 as for φ itself. We intend to calculate the corresponding scalar quasinormal modes (QNMs). The scalar QNMs will help us in understanding the stability of the spacetime under scalar perturbations and also give us an idea about whether we can distinguish the different geometries of the wormholes through the values of the scalar QNMs.
The Klein-Gordon equation for the scalar field is, As our background spacetime is spherically symmetric and static we use the following ansatz to decompose in terms of spherical harmonics, where the indices of Y (θ, φ) have been suppressed for simplicity.
Incorporating this in Eq. (24) we get the radial equation in the form of a Schrödinger-like equation in the tortoise coordinate , where and m is the azimuthal number arising from the separation of variables. For a general r ( ), it can easily be shown that Thus, for m = 0 (the s-wave) , V e f f = r r .
3.1 Effective potentials, energy conditions and single/double barriers An interesting fact to note about the m = 0 effective potential (for any r ( ) stated just above) is that it is linked to ρ( ) and τ ( ) via the Einstein equations for a general r ( ). We can easily note, Thus, to have an everywhere non-negative effective potential, we need to have ρ + τ ≤ 0 or p ≥ 0. The single or double barrier nature of the effective potential depends on the number of zeros of dp dr and their nature (maximum/minimum etc.). If dp dr has only one zero which is a maximum and since V e f f goes to zero asymptotically from the positive side if ρ + τ ≤ 0, we have a single barrier. On the other hand, if dp dr has three zeros (two maxima and one minimum) then we can have a double barrier, assuming r ( ) to behave as required for a wormhole. Thus, choosing a r ( ) and ensuring that dp dr has the right behaviour we can generate wormholes which will have a double barrier effective potential.
It is important to note that our n > 2 geometries necessarily have these properties. One can verify that the extrema for any value of n occurs at where = 0 is a minimum with r (0) = 0 and the other two symmetrically placed (about = 0) locations are maxima.
Many more examples can be worked out which exhibit double barriers. For example r ( ) = 2 + b 2 0 e −α 2 will yield a double barrier too for α = 0 and it will reduce to Ellis-Bronnikov for α = 0. It can be checked that none of the three standard wormhole spacetimes: spatial Schwarzschild where 1 n (maxima for m = 0) and . It is easy to note from the above expression that x 2 max < x 2 0 or |x max | < |x 0 |. This means the maxima shift towards the origin as m becomes larger and the value of V e f f at = 0 is gradually lifted upwards. For large m, the double barrier features almost (not fully) disappear and we end up being closer to a single barrier located at = 0. However, if n = 2, x 0 = 0 and the above equation has no solution apart from x max = 0, which is a maximum for all m.
We state below some features of the effective potential with reference to the associated figures.
• The plot in Fig. 6 shows the variation of V e f f ( ) w.r.t for different geometries at a particular value of m = 1  (Fig. 7) for higher n geometries, even for lower modes, shows that although the potential is still a double barrier, the potential curves for different geometries are nearly identical. • In Fig. 8, we find the potential being plotted once again for different geometries but with higher m value. It is interesting to observe that here the potential, for all geometries (i.e. for all n > 2), show an almost single barrier structure similar to the n = 2 case. We can see the two peaks getting flatter and merge into a 'nearly' single barrier when we zoom in very close to = 0 (see Fig. 9 for a zoomed plot of the effective potential-note the y-axis range and the scale here). • For all the potential plots we have taken the throat radius as unity. From Eq. (27) if we write the potential in terms of x = /b 0 , we observe that the V e f f varies as b −2 0 . So for higher values of b 0 the height of the potential peak will go on decreasing. Also, for any n > 2 geometry the width of the peak as well as their separation (for small m) increases for higher b 0 values.

Quasinormal modes
Given the features of the effective potentials, we now move on towards obtaining the scalar quasinormal modes in the background geometry of this family of wormholes. The existence of QNMs can be directly observed in the time evolution of the scalar field obtained by integrating the scalar wave equation following methods outlined in [86,87]. The wave equation is recast using light cone coordinates (du = dt − d and dv = dt + d ). Using appropriate initial conditions along the u and v lines we numerically integrate to obtain the timedomain profiles shown in Fig. 10(a, b). The damped ringing in time, exhibiting the decay of the scalar field is clearly visible in these plots. In order to find the quasinormal modes we use various available methods. Let us first discuss the methods briefly. Thereafter, we present our results and discuss their consequences in detail.

Methods for finding QNMs
As stated before, QNMs are complex frequencies associated with purely outgoing waves at spatial infinity. From Eq. (32) we see that ψ( , t) → e ∓iω e −iωt as → ∓∞ (note, for an asymptotically flat spacetime, V e f f → 0 at spatial infinity). Thus, if a QNM frequency has a negative imaginary part, the field will decay with time through these modes indicating a stable geometry.
We may find the scalar QNMs for different geometries using three different methods. The method of Prony fitting and direct integration are completely numerical techniques while the WKB method is semi-analytic. In the following, we discuss the direct integration method in detail and refer to well-known sources for the other two methods.
The time domain evolution of a signal can be used to extract the QNMs by fitting damped exponentials through the method of Prony fitting. This method has been applied multiple times in the literature (for details see [87]). We get a precision of upto 3 decimal places through this method. A possible source of error here is the lack of knowledge about the exact beginning of QNM ringing. However, this error does not really affect our results as we are interested only in the dominant mode obtained by fitting the late-time part of the signal with damped exponentials.  ψ(0, 0). For details on the discretization scheme, see [87] Another method for finding the QNMs is the semianalytical method or WKB approximation which was developed by Schutz and Will [88]. They calculated the QNM frequencies by taking the WKB solutions upto the eikonal limit which gives a simple analytical formula involving the parameters in the metric functions, where V 0 and V 0 denote the values of the effective potential and its second derivative at the maximum. p denotes the overtone number with p = 0 being the fundamental mode. In our work we have dealt with only the fundamental modes and compared some of the WKB values with the numerically obtained results. For a recent comprehensive review on WKB methods one may refer to [89]. We note that the WKB method is not perturbative. Hence, higher orders may not necessarily ensure better results. Also, one needs to keep in mind that the WKB formula used here is applicable for single barrier potentials i.e. when two turning points are involved. In our case though, we have double barriers which tend to almost single barriers for large m. Thus, WKB results will not be applicable in general (may be so, in an approximate sense, for large m) and the formula needs to be suitably modified for four turning points in order to obtain correct results. Finally, we discuss the direct integration method for finding QNMs, first given by Chandrasekhar and Detweiler [90]. Here, the differential equation (26) is numerically integrated using purely outgoing boundary conditions. Since the potential is symmetric about the throat of the wormhole at = 0, the solutions will be symmetric or anti-symmetric. Thus, we can impose additional conditions on the solutions namely: u(0) = 0 for anti-symmetric solutions and u (0) = 0 for symmetric solutions. This gives us two classes of QNMs corresponding to the condition imposed at = 0.
We observe that only e −iω wave will exist on the left side of the throat as there is no reflecting potential at = −∞. Hence the throat of the wormhole seems to play the role of the event horizon [68]. Ideally, we should integrate from = −∞ to ∞. However, since our potential is symmetric about = 0 we integrate in the positive half range of corresponding to r = b 0 ( = 0) to ∞ and then reflect the solution about the throat by using the conditions for obtaining symmetric or anti-symmetric solutions. In this way, we are taking into account the behavior of the solution for < 0 and we have the right boundary conditions (outgoing) at both the infinities in (i.e. → ±∞). A similar treatment can be found for another wormhole geometry in [71].
We begin with the series expansion of u at infinity given as, and obtain the coefficients A 1 , A 2 , . . . in terms of A 0 . We integrate the differential equation (26)    In our work, we will consider QNM frequencies related to the symmetric class of solutions, which have low damping. It may be noted that the error in this method arises due to the fact that the boundary condition imposed is not exactly at spatial infinity but at some finite point. This error may be minimised if we choose the matching point to be sufficiently far away from the throat. In our quoted values, we have a precision of six (6) digits for the imaginary part of the QNMs and five (5) digits for the real part.

Results from various methods
We have found the QNMs using each of the above numerical methods, assuming the throat radius b 0 = 1. The two tables (Tables 1, 2) list the fundamental QNM values calculated numerically for different modes, for the n = 4 and n = 10 geometries. We observe that the two numerical methods give nearly identical values and the results from both methods are suitable for use in further calculation. However, we use the QNM values from the Prony method in the next section, although using values from DI would lead to exactly the same observations.

Wormhole shapes from QNMs?
It is now reasonable to ask: is it possible, using QNMs to distinguish between geometries with different n values. The plot in Fig. 11 shows the variation of the real part of the fun- damental ω Q N M with the magnitude of its imaginary part, for different n. Each point in the plot for a particular n corresponds to its QNM frequency for a particular angular momentum mode. As we move from left to right, the value of m goes on increasing, so the left-most point corresponds to lowest m, while the rightmost point is for the highest m value. As mentioned before, the QNM values used in this plot have been calculated using the Prony method. From the plot of QNM frequencies in Fig. 11 we can draw the following conclusions: • For lower n, the geometries are distinguishable through their fundamental scalar QNMs. • When we move to higher n, the geometries begin to look nearly identical as evident from the QNMs (and also the effective potentials Fig. 7) for all modes. Hence it becomes difficult to distinguish different higher n geometries, solely from the QNMs. • Note that lower m modes are more suited for identifying the geometries. Higher m modes of all geometries have similar QNM values due to their nearly identical almost single barrier effective potentials (see Fig. 8).
• We also observe that as we go to higher n geometries the magnitude of the imaginary part of ω goes on decreasing. Thus the higher n wormhole families are likely to be relatively less stable owing to the fact that the perturbation for these wormholes takes a longer time to decay. • It is also clear that the n = 2 values are markedly different from those for other n. From a geometry standpoint, one is aware that the n = 2 geometry (Ellis-Bronnikov) is indeed special and has clear distinguishing features (see earlier discussion in Sect. 2).

Approximate Analytic Fit
In order to extract physically relevant information from our numerically obtained QNM values, we need an approximate model for ω which may be obtained by fitting the numerical data to analytical functions. We intend to study the variation of frequency (obtained from the real part of ω) with throat radius for different geometries and modes. We observe from Tables 1 and 2 that ω r increases with increasing m. The value of ω r also varies inversely with the throat radius. To imitate such a behavior we construct an approximate analytic model, where c is the speed of light, m is the corresponding mode which we want to fit, b 0 is in length units and the magnitudes of the coefficients a, g, h, k, p and the exponent d are obtained using NonLinearModel fit in Mathematica 10, which fits the above model with frequency corresponding to each mode, for each geometry, as obtained earlier using the Prony method ( Sect. 4). We have used six frequencies corresponding to six values of n for finding the coefficients and the exponent in the fit model. The analytical fits as obtained using Mathematica 10 for m = 1 and 10 has been discussed in detail in the Appendix. By looking at the explicit values of the coefficients as given in Eqs. (A.1) and (A.2) for modes m = 1 and 10 respectively, we can get an idea about the contribution made by each power of n, which is true for all modes. From Fig. 12(a, b) we note that the best-fit model matches the values of the frequencies quite well. If we refer to the values of the coefficients and the exponent in the fit model (given in the Appendix for m = 1 and 10), we find that, as expected, the higher powers of n contribute much less to the value of the frequency. We can now use the fit model to write the throat radius, in units of M (which is physically meaningful), corresponding to each geometry as with f (n, m) denoting the fitting function for each mode and geometry, G = 6.67 × 10 −11 m 3 kg −1 s −2 and M = 1.989 × 10 30 kg. The accuracy of the fit as given in Mathematica is upto 5 decimal places i.e. the results obtained from the fit exactly matches the ones obtained numerically, upto 5 decimal places. The plots shown in Fig. 12(a, b) have a precision of 5 decimal places. The fitting function (as given in the Appendix) has been used to generate the plots in Mathematica. Similarly, we can get a model fit for the imaginary part of the QNM frequency as well, to study the damping time. The approximate model can be taken as, ω i = c a n d + p n 2 + q n 3 + k n 4 + g n 5 b 0 m = cF(n, m) b 0 (38) where the magnitudes of the coefficients a, p, q, k and g and the exponent d are to be determined from the NonLinear-Model fit in Mathematica 10. As before, for the detailed fitting functions for modes m = 1 and 10, the reader is referred to the Appendix. The accuracy of this fitting model is the same as stated just above, for the real part. The damping time will be given by a relation The above relation can be used to compute the damping time for any geometry corresponding to any mode for a particular b 0 . We now try to find frequencies, using Eq. (37) (37)) and the frequency is plotted as a function of throat radius. In Fig. 13(a, b), we observe that for almost the same range of throat radius, the frequencies corresponding to various geometries have higher values for higher m modes. Hence, it is the lower modes which bear a closeness with the frequencies observable in the current generation of gravitational wave detectors. In the case of higher m values, the frequencies for geometries with different n are very similar (see Fig. 13b) over the entire range of the throat radius . Thus, even if we succeed in detecting the higher modes in some way, they will not help us in determining the corresponding geometry of the wormhole from the frequency. In contrast, for the lower m modes, the lower n geometries (n = 2, 4, 6) indeed have different frequencies for the same throat radius, thus making them distinguishable.
It is therefore fair to say that the QNMs can indeed be used as a tool to identify geometries with different n. It seems, for large n, in particular, difficulties appear. However, as we now illustrate, large n has another distinguishing featurethe occurence of echoes.

Observing echoes for large 'n' values
Earlier, we noted that for n > 2 and small m, the effective potential is a double barrier, as shown in Fig. 6. In the presence of such a potential, after the initial damped ringdown phase (i.e at later times), the transmitted wave  (37) emerges after getting successively reflected between the two potential peaks. Thus, the transmitted waves have a relative time delay along with a reduction in their peak amplitude. This is the well-known phenomenon of 'echoes in the timedomain profile' which is manifest in double barrier potentials. These echo signals appear as follow-up to the initial ringdown phase and dominate only at late times. They play an important role in distinguishing the spectrum of a black hole from exotic compact objects that are generally characterised by two potential surfaces, thus leading to the possibility of echoes. Echoes have been studied for a wide variety of cases, some of which can be found in [67,[91][92][93][94][95][96][97][98][99][100][101].
For our family of wormholes as we go to higher n geometries, we get distinct echo patterns. With an increasing n, the peaks of the effective potential become sharper resembling delta-functions, hence resulting in sufficient reflection of waves to generate prominent echoes. The time domain profiles for two geometries with different n values are shown in Fig. 14(a, b). Echoes are clearly visible in these profiles. Note that there is no particular reason for choosing n = 2000, 2600, except the fact that echoes are more prominent for large n.
For smaller n, echo patterns may be better observed with the help of a 'cleaning' procedure of the time domain profile as described in [96]. Ghersi et al. [96] have studied the scattering of wave packets from a Morris-Thorne [7] wormhole geometry which is constructed by joining two Schwarzschild geometries of equal mass at some r 0 > 2M. The potential barrier for such a wormhole is symmetric about the throat hence forming a cavity and harbouring echoes in the time evolution of a wave packet. Since the amplitude of the echoes is very small in comparison to the entire spectrum, one needs to subtract the effect of scattering due to the black hole (i.e. a single potential barrier) in order to clearly observe the echoes. Through this process the authors in [96] 'clean' the double barrier signal from the back-scattering occurring from the tail of the potential barrier, which, if not removed, suppresses the echo pattern. For scalar waves in our wormhole geometry we can perform a similar 'cleaning' procedure too. In small 'n' geometries, the potential barrier is smooth and hence has a wide tail which results in sufficient back scattering. Also the potential peaks do not have enough separation to provide strong reflections. So the echoes, even if they exist, are very weak and get damped easily. Thus, a similar 'cleaning' procedure might help us in observing the presence of echoes for smaller 'n' values which are normally not observable in the full spectrum. To proceed, we need a single barrier potential such that this single barrier, along with its mirror image, generates the double potential barrier. For this we consider only the peak present in the positive side. To create such a single barrier we have restricted the potential as given in Eq. (27) to ≥ 0. For < 0 we set the potential to a constant value 2 as the value of our double barrier potential around = 0 is always 2 for m = 1. Hence, the 'cleaned' profile will be given by: ψ clean = ψ doublebarrier − ψ singlebarrier . This will leave us with echoes produced by the double barrier. Even though such a cleaning process is non-unique as we can always multiply some constant factor to the ψ of the single barrier and still get the echo as a remnant of the subtraction process, it serves our purpose of demonstrating the existence of echoes for lower n geometries.
The significance of the cleaning procedure is evident in Fig. 15 where we observe that for n = 100 the echoes indeed exist and are observable after we clean the profile of the effect of the single barrier. We find that the echoes in Fig. 15b have an amplitude of the order of 10 −2 lower than the original spectrum.
Thus the existence of echoes in a full spectrum without 'cleaning' will be a tell-tale sign for a higher n geometry. For high values of n, the wormhole (a 2D slice embedded in 3D Euclidean space) resembles two flat sheets with holes which are connected by a cylinder. The wormhole throat is the radius of the cylinder which connects the flat regions. The effective potential for scalar waves looks close to a double delta-function potential, for which echoes are found. However, if we restrict ourselves to small n geometries, as is done in the QNM studies here, we can safely comment that the echoes do not appear to have a dominant presence or contribution (and can be found only after cleaning the spectrum).

Summary and concluding remarks
We summarize below the results obtained by us.
Working with a family of wormholes of which Ellis-Bronnikov spacetime is a special case, we first show how one may model the matter required for this generalised family. It turns out that it is possible to generate the members of this family of geometries, by using a massless phantom scalar field and extra matter (absent for EB spacetime) satisfying the ANEC.
Our next result is about the effective potentials for scalar perturbations. We obtain the general conditions on the metric function for which the effective potential could be a double barrier. Thereafter, we illustrate the conditions using the family of wormholes dealt with in this paper.
Subsequently, we focus on the scalar quasinormal modes for this family of Lorentzian wormholes. We obtain the QNM frequencies mainly using two methods: Prony fitting and direct integration and observe that both methods are equally suitable for our wormhole family.
Having found the fundamental QNMs, we attempt to distinguish between the members in this wormhole family, through their values and dependencies on metric parameters and azimuthal number m. The fact that the member wormholes of this family have different geometries for different  It turns that from the QNMs one may be able to distinguish between wormhole geometries, for lower n values. For higher n, we cannot unfortunately make such a distinction just from the QNMs, as the effective potentials become nearly identical, thereby making any such attempt, largely difficult. This inability to distinguish between different, higher n geometries through QNMs, is not a shortcoming of the methods used for finding them but is an innate feature of this wormhole family. Note also that the QNMs corresponding to lower m values are better suited for distinguishing between the geometries, as for higher m, the QNMs are similar in value due to similar effective potentials.
Finally, as another way to distinguish different geometries, we have studied the occurrence of 'echoes' through the time evolution of the field for different members ('n' values) of the family. We find that although all n > 2 geometries possess double barrier potentials, for low m modes distinct echoes appear only at very large values of n. This happens because of the fact that for small n wormholes the peaks of the potential are closely spaced and hence the reflections are not strong enough to produce clearly visible echoes. The echo signals get damped giving way to the dominant QNMs of the wormhole. Using a 'cleaning' procedure of subtracting the scattering of the single barrier from the full spectrum helps us in observing the echoes, if they exist, for such lower n geometries too. On the other hand for large n, the peaks are well separated and we get distinct echoes, even without cleaning.
Thus, in summary, it is fair to say that the QNMs can be used as a distinguishing tool for small n geometries while the presence of echoes in the ('uncleaned') time domain profile is a clear indication of a large n geometry.
We conclude with a remark. Our work, as reported here, provides a fairly thorough study of the geometry, matter, scalar QNMs and echoes in this family of ultrastatic wormholes. This is a first step in establishing this family of geometries as a viable template for wormholes. A useful future endeavour would be to analyse gravitational perturbations which will enable us to use GW observations directly to test our results and address the existence question of wormholes. We hope to pursue such investigations in future.

Expressions of fitting functions for real and imaginary parts of ω Q N M
In this brief appendix we provide the details of the approximate analytical fit for the QNMs, as obtained using Mathematica 10. The Appendix is a supplement for the Sect. 4.3.1. When the real part of QNM frequency is fitted to the approximate analytical model given in Eq. (36)