## Two distinct types of echoes in compact objects

Shui-Fa Shen 1 2 , , Kai Lin 3 4 , , Tao Zhu 5 , Yu-Peng Yan 6 7 , , Cheng-Gang Shao 8 , and Wei-Liang Qian 4 9 , ∗

- 1 School of intelligent manufacturing, Zhejiang Guangsha Vocational and Technical University of Construction, 322100, Jinhua, Zhejiang, China

2 School of Nuclear Science and Engineering, East China University of Technology, 330013, Nanchang, Jiangxi, China

- 3 Universidade Federal de Campina Grande, Campina Grande, PB, Brasil
- 4 Universidade de São Paulo, São Paulo, Brasil

5 Institute for Theoretical Physics and Cosmology, Zhejiang University of Technology, 310032, Hangzhou, China

6 School of Physics, Suranaree University of Technology, 30000, Nakhon Ratchasima, Thailand

7

Thailand Center of Excellence in Physics, Commission on Higher Education, 10400, Ratchathewi, Bangkok, Thailand

8 MOE Key Laboratory of Fundamental Physical Quantities Measurement, Hubei Key Laboratory of Gravitation and Quantum Physics, PGMF, and School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China and

9 Center for Gravitation and Cosmology, School of Physical Science and Technology, Yangzhou University, 225002, Yangzhou, Jiangsu, China

(Dated: Sept. 12th, 2024)

In the black hole perturbation theory framework, two di ff erent physical pictures for echoes in compact objects have been proposed. The first mechanism interprets echoes as repeated reflections of gravitational waves within a potential well, where the echo period is defined by twice the distance related to the spatial displacement operator that separates two local maxima of the e ff ective potential. The second mechanism associates echoes with a discontinuity in the e ff ective potential, potentially associated with specific accretion processes, without necessarily introducing a second local maximum in the e ff ective potential. This discontinuity leads to echo signals that are typically attenuated over time more quickly, with their period dictated by the characteristics of the transfer amplitudes. In both scenarios, the echoes correspond to a new category of quasinormal modes with minor real parts, with their period connected to the spacing between successive modes in the frequency domain. This work elaborates on a unified framework in compact stars that encompasses both echo mechanisms. It suggests that these two types of echoes derive from di ff erent physical origins and can be independently triggered. The occurrence and interplay between these two types of echoes are demonstrated through numerical simulations.

## I. INTRODUCTION

The inception of gravitational wave (GW) observation, facilitated by direct waveform captures [1-4], has initiated a groundbreaking phase in observational astronomy. This novel research domain is dedicated to exploring the strong-field gravity sector, which has experienced rapid progress lately. Notably, current endeavors in space-based laser interferometry [5-7] have led to significant advancements aimed at achieving unmatched detector sensitivity [8-12].

Theoretical insights have been instrumental in assessing the feasibility of black hole spectroscopy [13-18]. Specifically, in the context of real-world phenomena, sources of gravitational radiation such as black holes and neutron stars are embedded within a medium, not existing in isolation. This consideration inevitably alters spacetime from a perfect symmetric metric, substantially modifying the GWs emitted from such an entity. Such observations have directed research focus towards 'dirty' black holes [19-22], thus broadening the scope of black hole perturbation theory.

Within this framework, there has been a concentrated e ff ort on constructing realistic models involving compact astrophysical bodies, including binaries of black holes or neutron stars. The investigation of black hole quasinormal modes (QNMs) [23-25] has been a primary focus, given their critical role during the ringdown phase post-merger. These modes, which are dissipative oscillations reflective of the core characteristics of the black hole spacetime, adhere to several no-hair theorems [26, 27]. Initial studies by Leung et al. [20] on scalar QNMs of non-rotating 'dirty' black holes, evaluating shifts in quasinormal frequencies through generalized logarithmic perturbation theory, laid the groundwork in this area. Further analysis by Barausse et al. on perturbations surrounding a central Schwarzschild black hole [22] demonstrated that the resulting QNMscould exhibit significant deviations from an isolated black hole's characteristics. Nevertheless, they posited that the influence of the astrophysical environment on black hole spectroscopy could be minimal if appropriate waveform templates were applied. Among the scenarios examined in [22], the thin shell model was identified as having a significant e ff ect on the QNM spectrum. The concept of the black hole pseudospectrum, introduced by Nollert and Price [28, 29], aligns with these insights, showing how minimal perturbations, depicted as step functions, could remarkably influence high-overtone QNM modes, thus revealing a sizable susceptibility of the QNM spectrum to 'ultraviolet' perturbations. Our prior work [30] contended that discontinuities could alter the QNMspectrum's asymptotic behavior non-perturbatively, causing high-overtone modes to transition along the real axis rather than the imaginary frequency axis [31, 32], a phenomenon observed irrespective of the discontinuity's proximity to the horizon or its order. Through the notion of structural stability, Jaramillo et al. [33-35] explored the e ff ects of metric randomizations, with their findings via the hyperboloidal coordinates [36] indicating a shift of the pseudospectrum's boundary towards the real frequency axis, thereby emphasizing the inherent instability of highovertone modes to ultraviolet perturbations. Recent findings by Cheung et al. [37] also suggest the vulnerability of even the fundamental mode to destabilization by generic perturbations.

Cardoso et al. introduced the concept of echoes [38], a novel phenomenon intersecting with late-stage ringing waveforms, potentially distinguishing di ff erent gravitational systems through unique near-horizon characteristics. This concept has spurred numerous investigations [39-41] into echoes across various systems, including exotic compact objects such as gravastars [42, 43], wormholes [44-50], among others.

On the analytic aspect, echoes' derivation from Green's function's properties, especially its asymptotic behavior, was furthered by Mark et al. [51] through frequency domain evaluation. Echo

∗ E-mail: wlqian@usp.br

phenomena in compact objects were then conceptualized by reinterpreting the response waveform as a series summation of reflection and transmission amplitudes alongside the convolution integration process in the inverse Fourier transform. From the viewpoint of a scattering process, investigations of echoes within the framework of Damour-Solodukhin type wormholes [48] by Bueno et al. [52] involved solving for specific frequencies at which the transition matrix turns singular, leading to identifiable quasinormal frequencies that give rise to echoes. While most scenarios of echoes occur for e ff ective potentials possessing two local maxima, some of us proposed an alternative mechanism for echoes [53] without necessarily introducing a second local maximum in the e ff ective potential. The Green's function's approach was adapted, and the echoes' emergence is linked to the asymptotic pole structure of the QNM spectrum [30, 54] due to some minor discontinuity planted into the potential. This suggests that discontinuity might serve an alternative ingredient for echo mechanism.

In principle, it can be argued that both approaches stand on an equivalent theoretical foundation concerning a peculiar property of the transit amplitudes. However, from a physical perspective, they can be viewed as two distinct mechanisms for generating echoes. In the first scenario, the e ff ective potential is characterized by two local maxima separated by a distance, creating a potential well [52, 55], with Damour-Solodukhin type wormholes [52] serving as a prime example. Echoes in this context can be intuitively understood as resulting from repeated reflections of GWs within such a potential well, with the echo period mathematically determined as twice the distance between the maxima of the e ff ective potential via the spatial displacement operator that separates the two local maxima. The second scenario entails a degree of discontinuity within the e ff ective potential [51, 53], potentially brought about by specific accretion processes giving rise to 'cuspy' [56-58] or 'splashback' [59, 60] profiles, without necessarily leading to a second local maximum in the e ff ective potential. This gives rise to echo signals typically attenuated over time, with their period dictated by the characteristics of the relevant transfer amplitudes. For both cases, the echoes correspond to a novel category of quasinormal modes with minor real parts, and the echo period is associated with the spacing, along the direction of the real frequency axis, between successive modes [53].

The present study aims to explore this topic further, proposing a unified model within compact stars that encompasses both echo mechanisms. We demonstrate that the two types of echoes, stemming from distinct physical origins, can be independently triggered within such a model, as evidenced through numerical simulations. The remainder of the work is structured as follows. In Sec. II, we delve into Green's function approach for the two distinct echo mechanisms. Subsequently, in Sec. III, we provide an overview of axial gravitational perturbations in compact stars, serving as a comprehensive theoretical setup for the emergence of echoes through both mechanisms. The master equation for the perturbations is formulated for compact stars of uniform density. Numerical results are presented in Sec. IV, where spatial-temporal evolutions are calculated using the finite di ff erence method, elucidating the principal characteristics of the echoes and the interplay between the two mechanisms. Additional discussions and concluding remarks are given in Sec. V.

## II. TWOPHYSICALPICTURESFOR ECHOESIN COMPACT OBJECTS

This section explores Green's function formalism within black hole perturbation theory and its application to echoes in compact objects. When variable separation is achievable, the evolution of perturbations simplifies, with the dynamics predominantly governed by the radial component of

the master equation [24, 25]:

$$\i i m _ { \lfloor \leftarrow \tau, \, \leftarrow$$

where x denotes the tortoise coordinate, and V e ff represents the e ff ective potential defined by the spacetime metric, spin ¯, s and angular momentum /lscript of the waveform. For example, the ReggeWheeler potential V RW for a Schwarzschild black hole metric is given by

$$\intertext { a r R W } & \text{for a $v$ in a $v$$

where

$$F = 1 - r _ { h } / r,$$

rh = 2 M represents the horizon, with M being the black hole mass, and the tortoise coordinate relates to the radial coordinate r through x = ∫ dr / F r ( ).

QNMsof black holes are determined by solving the eigenvalue problem presented in Eq. (1) in the frequency domain:

$$\begin{array} {$$

For asymptotically flat spacetimes, the boundary conditions are defined as:

$$\text{acetimes, the boundary conditions are defined as} \colon \\ \Psi \sim$$

 indicating an ingoing wave at the horizon and an outgoing wave at infinity, with n denoting the overtone number. Beyond the initial burst, the waveform Ψ is characterized by quasinormal oscillations with complex eigenvalues ω n , termed quasinormal frequencies. The dissipative nature of the system is reflected in the imaginary components of these frequencies.

The QNM properties are intricately linked to the analytic properties of the corresponding Green's function, satisfying

Following established procedures [61-63], Green's function is constructed as:

$$\text{on, saxifying} \\ \left [ \frac { d ^ { 2 } } { d x ^ { 2 } } + ( \omega _ { n } ^ { 2 } - V _ { \text{eff} } ( r ) ) \right ] \widetilde { G } ( \omega, x, y ) = \delta ( x - y ). \\ \text{blished procedures } [ 6 1 - 6 3 ], \text{Green's function is constructed as:}$$

where x &lt; ≡ min( x y , ), x &gt; ≡ max( x y , ), and

$$\widetilde { G } ( \omega, x, y ) = \frac { 1 } { W ( \omega ) } f ( \omega, x _ { < } ) g ( \omega, x _ { > } ), \\ y ), x _ { > } \equiv \max ( x, y ), \, \text{and}$$

$$W ( \omega ) \equiv W ( g, f ) = g f ^ { \prime } - f g ^ { \prime }$$

is the Wronskian of f and g , the two linearly independent solutions of the corresponding homogeneous equation satisfying the boundary conditions Eq. (5) at the horizon and infinity, respectively. To be specific, f and g demonstrate the following asymptotic behaviors

$$\text{ecific, $f$ and $g$ demonstrate the following asymptotic behaviors} \\ f ( \omega, x ) & \sim \begin{cases} e ^ { - i \omega x }, & x \to - \infty, \\ A _ { \text{out} } ( \omega ) e ^ { + i \omega x } + A _ { \text{in} } ( \omega ) e ^ { - i \omega x }, & x \to + \infty, \end{cases}$$

and

 in spacetimes that are asymptotically flat, remaining bounded as t → + ∞ for /Ifractur ω &lt; 0. Herein, A in , A out , B in , and B out represent the reflection and transmission coe ffi cients, which, although their exact expressions might not be specified in our derivations, are well-defined for a specific metric. These waveform amplitudes adhere to the relations [64]

$$g ( \omega, x ) & \sim \begin{cases} B _ { \text{out} } ( \omega ) e ^ { + i \omega x } + B _ { \text{in} } ( \omega ) e ^ { - i \omega x }, & x \to + \infty, \\ e ^ { + i \omega x }, & x \to + \infty, \end{cases} \\ \text{nes that are asymptotically flat, remaining bounded as $t \to +\infty$ for $\Im\omega<0$. Here,} \\ \text{$i_{in}$, and $B_{out}$ represent the reflection and transmission coefficients, which, although their} \\ \cdot \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{cases}$$

$$B _ { \text{out} } \, & = \, A _ { \text{in} }, \\ B _ { \text{in} } \, & = \, - A _ { \text{out} } ^ { * },$$

as dictated by principles of completeness and conservation of flux. Moreover, the reflection and transmission coe ffi cients for black holes are specified as for waves originating from x →-∞ .

$$\begin{array} { c c c }, \\ \widetilde { \mathcal { R } } _ { B H } ( \omega ) & = & \frac { B _ { \text{in} } } { B _ { \text{out} } }, \\ \widetilde { \mathcal { T } } _ { B H } ( \omega ) & = & \frac { 1 } { B _ { \text{out} } }, \\ \to - \infty. \\ \text{identified by the poles of Green's function. While the inherent pole} \end{array}$$

Echoes, akin to QNMs, are identified by the poles of Green's function. While the inherent pole structure primarily emerges from the zeros of the Wronskian Eq. (8), it may undergo alterations due to phenomena like pole skipping [65, 66] or additional distortions introduced by an external source [67]. Echoes are categorized within this context. The discussion outlines two scenarios in which particular e ff ective potentials lead to a unique set of poles associated with echoes in compact objects. For the first scenario, the Wronskian is determined by segregating two potential barriers by a specified distance using a spatial displacement operator, interpreting the resultant echoes as repeated GW reflections within a potential well, with the echo period defined by twice the distance separating the two local maxima. The second scenario incorporates discontinuity at a certain order of the e ff ective potential, where the Wronskian is derived from the explicit forms of transmission amplitudes at the discontinuity point. In both instances, it is demonstrated that the resulting mathematical expressions share similarities, and the echoes correspond to a new series of quasinormal modes characterized by minor real parts, with the echo period correlating to the spacing between successive modes in the frequency domain.

## II.1. The first type of echo

In this subsection, we explore Damour-Solodukhin wormholes as a paradigm for the first type of echo. The waves described by Eqs. (9) and (10) can be associated with two distinct e ff ective potentials, V (1) and V (2) , respectively. On the right-hand side of V (1) and left-hand side of V (2) , the asymptotic forms of f and g are as follows:

$$f ( \omega, x ) \sim A _ { \text{out} } ^ { ( 1 ) } ( \omega ) e ^ { + i \omega x } + A _ { \text{in} } ^ { ( 1 ) } ( \omega ) e ^ { - i \omega x },$$

and

$$g ( \omega, x ) \sim B _ { \text{out} } ^ { ( 2 ) } ( \omega ) e ^ { + i \omega x } + B _ { \text{in} } ^ { ( 2 ) } ( \omega ) e ^ { - i \omega x },$$

in an asymptotically flat spacetime.

Accounting for the potential well involves shifting the e ff ective potential V (2) rightward by a distance L , leading to

$$\tilde { g } ( \omega, x ) & \sim \tilde { B } _ { \text{out} } ^ { ( 2 ) } ( \omega ) e ^ { + i \omega x } + \tilde { B } _ { \text{in} } ^ { ( 2 ) } ( \omega ) e ^ { - i \omega x },$$

where the transformed outgoing wave results from appropriately applying the spatial displacement operator in the frequency domain:

$$\hat { P } ( L ) = e ^ { i \omega L },$$

to the equation above. Specifically, it follows that

$$\tilde { B } _ { \text{in} } ^ { ( 2 ) } ( \omega ) & \ = \ B _ { \text{in} } ^ { ( 2 ) } ( \omega ), \\ \tilde { B } _ { \text{out} } ^ { ( 2 ) } ( \omega ) & \ = \ e ^ { 2 i \omega L } B _ { \text{out} } ^ { ( 2 ) } ( \omega ).$$

Utilizing these equations, the Wronskian is found to be

$$\text{use equivalently, use with formulas, is to unit to use} \\ W ( \omega ) = - i \omega \left [ A _ { \text{out} } ^ { ( 1 ) } ( \omega ) B _ { \text{in} } ^ { ( 2 ) } ( \omega ) - e ^ { 2 i \omega L } A _ { \text{in} } ^ { ( 1 ) } ( \omega ) B _ { \text{out} } ^ { ( 2 ) } ( \omega ) \right ]. \\ \text{mergence of echoes is inferred from the zeros of the obtained Wronskian, which are poles} \\ \text{s.e. to use} \quad \text{s.e. to use} \quad \text{s.e. to use} \quad \text{s.e. to use} \quad \text{s.e. to use}.$$

The emergence of echoes is inferred from the zeros of the obtained Wronskian, which are poles of the Green's function. Assuming the original e ff ective potentials V (1) and V (2) do not inherently produce echoes, the frequency ω 's amplitudes A (1) in ( ω ), A (1) out ( ω ), B (2) in ( ω ), and B (2) out ( ω ) are moderate. In regions far from their zeros and poles, they are considered constant. This premise leads to a new sequence of zeros for the Wronskian, corresponding to poles in the Green's function, distributed evenly and parallel to the real frequency axis. Under the above apprximation, if ω is a root, then ω + n π/ L (for any integer n ) is also a root, forming a branch of the QNM spectrum. The impact of these poles on the time-domain Green's function results in echoes. The spacing between successive poles set at 2 L , and the proximity of this branch to the real axis define, respectively, the echo period and observational viability.

## II.2. The second type of echo

Now, we turn our attention to the second type of echo. In this case, the emergence of echoes is closely associated with the presence of a discontinuity in the e ff ective potential. The discontinuity that triggers these echoes may be positioned near what would traditionally be considered the horizon or at a more distant location relative to the compact object. The core principles underlying the formulation of these echoes remain consistent across di ff erent physical systems. As discussed previously, the mathematical framework for deriving echoes from exotic compact objects has been established in [51, 53].

For simplicity, let us consider the discontinuity to be near the horizon. For such a case, the ingoing wave is slightly altered by a fraction of the outgoing wave, while the latter remains una ff ected. Specifically, we describe this interaction as

˜ ( f ω, x ) = f ( ω, x ) + ˜ C ( ω ) g ( ω, x ) , (19) representing a mixture of the ingoing and outgoing waves, as characterized by the asymptotic behaviors delineated in Eqs. (9) and (10).

Assuming the vicinity of the discontinuity, denoted by xc , features a significantly suppressed e ff ective potential V = Vc ∼ 0, the ingoing wave f can be approximated as a superposition of two plane waves

$$\tilde { f } ( \omega, x ) \in e ^ { - i \omega ( x - x _ { c } ) } + \widetilde { \mathcal { R } } ( \omega ) e ^ { i \omega ( x - x _ { c } ) },$$

where the reflection amplitude ˜ R is primarily determined by the unique characteristics of the compact object.

The relationship between ˜ C ( ω ) and R ( ω ) is given by

˜ ˜ with the reflection and transmission amplitudes defined by Eq. (12). The derivation of Eq. (21) involves aligning the coe ffi cients of the wave functions at an arbitrary asymptotical point where their expressions converge to plane waves, detailed further in Appx. A.

$$\text{between $\widetilde{C}(\omega)$ and $\widetilde{R}(\omega)$ is given by} \\ \widetilde { C } ( \omega ) & = \frac { e ^ { - 2 i \omega x _ { c } } \widetilde { \mathcal { T } } _ { \text{BH} } ( \omega ) \widetilde { \mathcal { R } } ( \omega ) } { 1 - e ^ { - 2 i \omega x _ { c } } \widetilde { \mathcal { R } } _ { \text{BH} } ( \omega ) \widetilde { \mathcal { R } } ( \omega ) }, \\ \text{and transmission amplitudes defined by $Eq.$ (12). The derivation of $Eq.$ (21)$} \\ \text{e coefficients of the wave functions at an arbitrary asymptotical point where}$$

Consequently, the frequency domain Green's function, as given by Eq. (7), is computed to be

$$\text{|\text{ently, the frequency domain Green}'s function, as given by $Eq.$(7), is computed to be} \\ \widetilde { G } ( \omega, x, y ) & \, = \, \frac { \tilde { f } ( \omega, y ) g ( \omega, x ) } { W ( g, \tilde { f } ) } \\ & \, = \, \frac { f ( \omega, y ) g ( \omega, x ) } { W ( g, f ) } + \widetilde { C } ( \omega ) \frac { g ( \omega, y ) g ( \omega, x ) } { W ( g, f ) } \\ & \, \equiv \, \widetilde { G } _ { B H } + \widetilde { C } ( \omega ) \frac { g ( \omega, y ) g ( \omega, x ) } { W _ { B H } }, \\ \intertext { i } = W ( g, f ) \, \text{ denotes the Wronskian, and } \widetilde { G } _ { B H } \, \text{ is the Green}'s function corresponding to} \text{|\text{ black hole metric. This formulation suggests that the QNMs inherent to the black hole} }. \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \.$$

Once more, the echoes can be analyzed by assessing the poles of the above Green's function. In this framework, the manifestation of echoes is linked to the phase shift embedded in ˜ C ( ω ) [51, 53, 68], which explicitly depends on the location of discontinuity xc . Again, as the original black hole metric does not imply any echo phenomenon, the quantities ˜ T BH and ˜ R BH are moderate functions of the frequency ω and will be treated as constants when analyzing the echo modes. For a simplified scenario where the discontinuity is implemented by truncating the e ff ective potential at xc from the inside, the reflection amplitude is found to possess the following form where W BH = W g f ( , ) denotes the Wronskian, and ˜ G BH is the Green's function corresponding to the original black hole metric. This formulation suggests that the QNMs inherent to the black hole largely persist within the wormhole metric, with echoes attributed to the poles in ˜ C ( ω ) as defined by Eq. (21).

$$\text{ampute is toum} & \text{to possess one following form} \\ \widetilde { \mathcal { R } } ( \omega ) = \frac { \sqrt { \omega ^ { 2 } - V _ { c } } - \omega } { \sqrt { \omega ^ { 2 } - V _ { c } } + \omega }, \\ \text{nation, whose derivation is also relegated to Appx. A. Using Eq. (21), it} \\ r & \text{ asvmntotal modes} \, \omega \gg 1. \, \text{we have}$$

using the WKB approximation, whose derivation is also relegated to Appx. A. Using Eq. (21), it is readily verified that for asymptotical modes ω /greatermuch 1, we have

$$\widetilde { C } ( \omega ) \sim \frac { 1 } { 2 \omega } e ^ { - 2 i \omega x _ { c } } \widetilde { \mathcal { T } } _ { \text{BH} } ( \omega ) V _ { c } \ll 1. & & ( 2 4 ) \\ \text{games are mainly reminiscent of the preceding subsection. Specifically, if} \\ ( 2 1 ) \text{ and therefore the Green's function Fa. (22). } \omega + n \pi / x. \text{ (where $n$ is an}$$

The succeeding arguments are mainly reminiscent of the preceding subsection. Specifically, if ω is a pole of Eq. (21) and therefore the Green's function Eq. (22), ω + n π/ xc (where n is an arbitrary integer) is also a pole. The union of such poles gives rise to an additional branch of QNMs, contributing to echoes of a period of 2 xc 1 .

1 Rather than representing an absolute coordinate, we note that xc essentially measures a relative distance. It quantifies the o ff set between the location of the e ff ective potential, relative to which the transmission coe ffi cients are defined and evaluated, and the discontinuity. This has been implied in the preceding derivations.

Before closing this section, we observe a nuanced distinction between the two scenarios. The first scenario does not inherently involve any discontinuity within the metric. Here, the echo period is solely determined by the separation of two local maxima within an essentially continuous e ff ective potential. In our derivation, we have assumed that the waveforms within the region spanning these maxima are asymptotic in the sense that they can be expressed as linear combinations of plane waves. This assumption is crucial for the mathematical assessment of the Wronskian. Conversely, the second scenario does not necessitate the presence of a secondary maximum within the e ff ective potential. Instead, it posits that the discontinuity is minor and sits far from the potential's peak, permitting only a negligible portion of the plane wave to move in reverse to merge into the initially outgoing wave. In this case, the echoes hinge critically on the singularity of the coe ffi cient of the added waveform, whose assessment generally requires some form of model-based assumption or simplification. Nevertheless, in both scenarios, the emergent QNM spectrum stems from a term of the form e 2 i ω L , under the presumption that the other terms are moderate frequency functions, notably lacking any exponential dependency on ω . These scenarios are elucidated by examining the zeros in the Wronskian Eq. (18) or the denominator of Eq. (21). Specifically, in the first scenario, the pertinent contribution derives from altering the transit amplitude as indicated in the second line of Eq. (17), as a result of the translation operation. This scenario suggests echoes result from the GWs ricocheting between the two local maxima of the e ff ective potential. In contrast, the second scenario's contribution is linked to the singularity within a minor fraction of the outgoing wave, as delineated by Eqs. (19) and (21), underscoring the distinct physical implications of the discontinuity introduced. As indicated in Sec. IV, echoes emanating from such conditions are usually more substantially attenuated over time.

## III. THE MASTER EQUATION IN THE COMPACT STAR OF UNIFORM DENSITY

The above two types of echoes might be present simultaneously in the context of compact objects. This is because the e ff ective potential of a compact star may possess two local maxima, while the star's surface usually introduces some discontinuity. In the literature, the well-known w-modes [69-72] are readily recognized as echoes in star quasinormal modes. However, to our knowledge, an explicit discussion regarding the nuance di ff erence elaborated in the foregoing section has yet to be carried out.

In what follows, we derive the master equation for axial gravitational perturbations in stars of uniform density. To enhance the reader's understanding, we adopt a minimalist approach that focuses on a spherically symmetric compact star with uniform density, making the physical content more accessible. The spacetime metric possess the form

$$d s ^ { 2 } = - h ( r ) d t ^ { 2 } + \frac { d r ^ { 2 } } { f ( r ) } + r ^ { 2 } \left ( d \theta ^ { 2 } + \sin ^ { 2 } \theta d \varphi ^ { 2 } \right ). \\ \text{outside of the star $r>r_{\lambda}$ where $r_{\lambda}$ is the radial coordinate of the star's surface. the}$$

On the outside of the star r ≥ rb , where rb is the radial coordinate of the star's surface, the metric is essentially Schwarzschild, which reads

$$f ( r ) = h ( r ) = 1 - \frac { 2 M _ { S } } { r },$$

where MS is the mass of the compact object. The tortoise coordinate reads x ≡ r ∗ = ∫ dr / √ f h = r + rS log ( r rS -1 . )

Beneath the surface of a star r ≤ rb with uniform density ρ = 3 MS 4 π r 3 b , one has [73]

$$& \text{the surface of a star } r \leq r _ { b } \text{ with uniform density } \rho = \frac { 2 M s } { 4 \pi r _ { b } }, \text{ one has } [ 7 3 ] \\ & \text{f} ( r ) \, = \, 1 - \frac { 2 M ( r ) } { r }, \\ & \text{h} ( r ) \, = \, \frac { 1 } { 4 } \left ( 3 \, \sqrt { 1 - \frac { 2 M ( r ) } { r } } - \sqrt { 1 - \frac { 2 M ( r ) r ^ { 2 } } { r _ { b } ^ { 3 } } } \right ) ^ { 2 }, \\ & M ( r ) \, = \, \int _ { 0 } ^ { r } 4 \pi r ^ { \prime 2 } \rho d r ^ { \prime } = M _ { S } \frac { r ^ { 3 } } { r _ { b } ^ { 3 } }, \\ & \text{p} ( r ) \, = \, - \rho \left ( \frac { \sqrt { 1 - \frac { 2 M ( r ) } { r } } - \sqrt { 1 - \frac { 2 M ( r ) r ^ { 2 } } { r _ { b } ^ { 3 } } } } { 3 \sqrt { 1 - \frac { 2 M ( r ) } { r } } - \sqrt { 1 - \frac { 2 M ( r ) r ^ { 2 } } { r _ { b } ^ { 3 } } } } \right ). \\ & \text{ise coordinate, the constant of integration is chosen so that } x ( r = 0 ) = 0. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

For the tortoise coordinate, the constant of integration is chosen so that x r ( = 0) = 0.

$$\ = \ - \rho \left ( \frac { \sqrt { 1 - \frac { 2 M ( r ) } { r } } - \sqrt { 1 - \frac { 2 M ( r ) r ^ { 2 } } { r _ { b } ^ { 3 } } } } { 3 \sqrt { 1 - \frac { 2 M ( r ) } { r } } - \sqrt { 1 - \frac { 2 M ( r ) r ^ { 2 } } { r _ { b } ^ { 3 } } } } \right ). \\ \text{constant of integration is chosen so that } \lambda \\ \text{aration of variables, the axial gravitation}$$

By using the method of separation of variables, the axial gravitational perturbations are governed by [74]

$$\frac { \partial ^ { 2 } \Psi } { \partial r ^ { 2 } } - \frac { \partial ^ { 2 } \Psi } { \partial t ^ { 2 } } - V ( r ) \Psi = 0, \\ V _ { \text{star} } ( r ) & = \frac { h } { r ^ { 3 } } \left ( \ell ( \ell + 1 ) r + 4 \pi ( \rho - p ) r ^ { 3 } - 6 M \right ), \\ V _ { \text{BH} } ( r ) & = \frac { h } { r ^ { 3 } } \left ( \ell ( \ell + 1 ) r - 6 M _ { S } \right ), \\ \text{anular momentum} \text{ } At \text{cial infinity} \text{ } the boundary condition dictates that the}$$

where L is the angular momentum. At spatial infinity, the boundary condition dictates that the waveform Ψ is an asymptotic outgoing wave. At the star's center, it must be regular with vanishing flux. Moreover, the junction conditions at the star's surface reads

$$\Psi _ { \text{inside} } ( r _ { b } ) \ & = \ \Psi _ { \text{outside} } ( r _ { b } ), \\ \partial _ { r } \Psi _ { \text{inside} } ( r _ { b } ) \ & = \ \partial _ { r } \Psi _ { \text{outside} } ( r _ { b } ),$$

or

$$\partial _ { x } \Psi _ { i n s i d e } ( r _ { b } ) = \partial _ { x } \Psi _ { o u s i d e } ( r _ { b } ).$$

## IV. NUMERICAL RESULTSABOUT THE EMERGENCE OF TWO TYPES OF ECHOES

The numerical calculations will be performed using the finite di ff erence method [75-77]. We will adopt the scheme recently proposed by some of us [78]. In the literature, most implementations are performed in Eddington-Finkelstein coordinates ( u = t -x , v = t + x ), where one implements the initial condition on v = 0 and the boundary condition on u = 0, In such circumstances, in the far region where the e ff ective potential asymptotically approaches zero, the initial perturbations will never attain the u axis due to causality constraints. However, in the region where the e ff ective potential does not vanish, there is a chance that the speed of signal propagation exceeds the unit, and initial perturbations placed on the u axis might traverse the v axis. In other words, the boundary condition enforced for u = 0 might interfere with the initial perturbations' free propagation.

Figure 1. The grid layouts of finite di ff erence method employed in the present study. As described in the text, the calculations are implemented in ( r ∗ ≡ x t , ) coordinates. The black points indicate the grids to which the initial and boundary conditions are assigned. The finite di ff erence method uses Eqs. (33) and (35) to evaluate the remaining grids in red.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00971v2-figure-2c5ffd79-81a9-45c0-b62c-aad83988810d.jpg)

Conversely, in our proposed scheme, the calculations are implemented in the ( x t , ) coordinates, and the initial condition is placed on the t = 0 slice. As a result, the concerns regarding causality will no longer be relevant. In particular, we adopt the initial conditions Ψ ( t = 0) of a Gaussian form and ∂ t Ψ ( t = 0). One proceeds to discretize the spacetime coordinates and approximates the partial derivatives by first-order finite di ff erences. To be specific, we denote

$$t _ { i } \ & = \ t _ { 0 } + i \Delta t, \\ x _ { j } \ & = \ x _ { 0 } + j \Delta x$$

and

$$\Psi _ { j } ^ { i } \, & = \, \Psi ( t = t _ { i }, x = x _ { j } ), \\ V _ { j } \, & = \, V ( x = x _ { j } ),$$

and therefore, the master equation Eq. (28) becomes

$$\cdots \dots \dots \dots \dots \dots \cdots \cdots \cdots \cdots \\ \Psi _ { j } ^ { i + 1 } = - \Psi _ { j } ^ { i - 1 } + \frac { \Delta t ^ { 2 } } { \Delta x ^ { 2 } } \left ( \Psi _ { j - 1 } ^ { i } + \Psi _ { j } ^ { i - 1 } \right ) + \left ( 2 - 2 \frac { \Delta t ^ { 2 } } { \Delta x ^ { 2 } } - \Delta t ^ { 2 } V _ { j } \right ) \Psi _ { j } ^ { i }, \\ \text{The junction condition Ea. (29)} \, \text{reads}$$

The junction condition Eq. (29) reads

$$\frac { \Psi _ { j _ { b } } ^ { i } - \Psi _ { j _ { b } - 1 } ^ { i } } { \Delta x } = \frac { \Psi _ { j _ { b } + 1 } ^ { i } - \Psi _ { j _ { b } } ^ { i } } { \Delta x },$$

where one uses the subscript b to denote the grid on the star's surface. Eq. (34) can be simplified to give

$$\Psi _ { j _ { b } } ^ { i } & = \frac { 1 } { 2 } \left ( \Psi _ { j _ { b } + 1 } ^ { i } + \Psi _ { j _ { b } - 1 } ^ { i } \right ), \\ \text{dicates that the grid is on the star's surface.}$$

where the subscript b indicates that the grid is on the star's surface.

As shown in Fig. 1, the numerical iterative process can be performed by using Eqs. (33) and (35). Given the black grids on the boundary, the temporal evolution is implemented by inferring the values of the red grids. For most grids, Eq. (33) is utilized to determine the grid values

♣

r

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00971v2-figure-b289da09-00ad-4279-903f-70ed9e96a5df.jpg)

r

Figure 2. The metric functions f ( r ) and h r ( ) (top-left), the mass function M r ( ) (top-right), density ρ ( r ) (bottom-left), and pressure p r ( ) (bottom-right) of the compact star in question shown as a function of the radial coordinate r . The plots are obtained using the metric parameters MS = 1 2 and rb = 1 13. .

✵/a0✁

✂✄☎

❱

✆

✝

✡

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00971v2-figure-bc1a2300-5287-40ad-9adf-48158dafbafd.jpg)

✎

✏

Figure 3. The e ff ective potentials of the compact star as a function of the radial coordinate r (left column) and tortoise coordinate x (right column). The top row is obtained using the metric parameters MS = 1 2 , rb = 1 13, and . /lscript = 2, while the bottom row is evaluated using MS = 1 2 , rb = 1 14, and . /lscript = 2. For both cases, the e ff ective potential features two local maxima at r = 0 and r = 1 65 while a discontinuity is also present . at rb . In the tortoise coordinates, the distances between the two maxima for the two sets of parameters are L = 40 (top) and 32 (bottom).

▲

✑

✒

✔

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00971v2-figure-4f35dcdb-09d5-4b22-ad44-6ce4f8155718.jpg)

✕

✖

Figure 4. The echoes for the two sets of metric parameters whose e ff ective potentials are given in Fig. 3. The spatial-temporal evoluations are shown in the left column, while the temporal evolutions are presented in the right column. The echo periods are found to be T /similarequal 80 (top) and 64 (bottom), in accordance with twice the distance between the two maxima, while irrelevant to the position of the discontinuity.

for the next time step, except for those on the star's surface. The latter is determined by employing Eq. (35). To avoid the von Neumann instability [79, 80], we choose

$$\frac { \Delta t ^ { 2 } } { \Delta x ^ { 2 } } + \frac { \Delta t ^ { 2 } } { 4 } V _ { \max } < 1.$$

Now, we proceed to present the numerical results. For the first types of echoes, we elaborate two sets of parameters: MS = 1 2 , rb = 1 13, . /lscript = 2, and MS = 1 2 , rb = 1 14, . /lscript = 2. The metric functions are presented in Fig. 2, and the corresponding e ff ective potentials are shown in Fig. 3. When presented in radial coordinate r , the forms of the e ff ective potentials are primarily identical, particularly the positions of the two maxima. However, the di ff erence is apparent in terms of the tortoise coordinate x . For both parameters, echoes are demonstrated, as shown in Fig. 4. The spatial-temporal evolution indicates that initial perturbations propagate at v = 1 and are reflected at both maxima. The temporal evolution clearly shows the occurrence of echoes. The distance between the two maxima governs the resulting echoes' period. Even though a discontinuity is present, it is irrelevant to the observed echoes.

We consider the metric parameters MS = 1 2 , rb = 100, and /lscript = 2 to demonstrate the second type of echo. As shown in Fig. 5, since the discontinuity is placed beyond the would-be maximum of Regge-Wheeler e ff ective potential, the resulting e ff ective potential does not feature a second maximum. Moreover, the discontinuity is visually insignificant. The resulting echoes are presented in Fig. 6. The echoes' period is found to be T /similarequal 200, in accordance with twice the distance between

✗

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00971v2-figure-f2f7eaa5-dd5c-4b57-ad97-e0ed12e9676b.jpg)

✙

✚

Figure 5. The e ff ective potentials of the compact star as a function of the radial coordinate r (left) and tortoise coordinate x (right). The results are obtained using the metric parameters MS = 1 2 , rb = 100, and /lscript = 2. The e ff ective potential does not possess any local maximum but features a discontinuity at the star's surface. In both coordinates, it is apparent that the discontinuity is visually insignificant.

✛

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00971v2-figure-02af1a73-e3f1-44c6-9d5b-246d59712dfa.jpg)

✜

✢

Figure 6. The echoes for the third set of metric parameters whose e ff ective potentials are given in Fig. 5. The spatial-temporal and temporal evoluations are shown in the left and right panels. The echoes' period is found to be T /similarequal 200, in accordance with twice the distance between the maximum and the discontinuity of the e ff ective potential.

the maximum and the discontinuity of the e ff ective potential. In contrast to Fig. 4, the magnitude of the echoes is suppressed in time. This explains why such echoes were not observed in Fig. 4, as it is likely buried inside the echoes formed by the first type.

## V. CONCLUDINGREMARKS

Echoes, emerging as a novel observable in the GW spectrum, o ff er a unique window into the physics of compact objects and their environments. By exploring the nuances of echo mechanisms through a unified model, this study contributes to a deeper understanding of the complex interplay between GWs and the structure of compact objects. Through numerical simulations, we demonstrate the possibility of independent triggering of distinct echo types, highlighting the diversity of physical processes that can manifest in the GW signals from compact stars. In particular, for the simplified model adopted in the present study, the more attenuated echoes might actually occur earlier. Therefore, it might still be observed if the first type of echo also exists and will be more persistent during evolution. Regarding the ongoing spaceborne GW detection programs, the find-

ings presented herein not only shed light on the theoretical underpinnings of echo phenomena but also pave the way for future observational strategies to detect and interpret echoes in GW data.

For a given total mass, the results presented in Sec. IV indicate that the second type of echo becomes relevant only when the star's size is significant. Specifically, the criterion for the emergence of the second type of echo is that the star's surface lies outside the would-have-been local maximum of the Regge-Wheeler e ff ective potential, i.e., rb /greaterorapproxeql r max , where r max = 3 M in the eikonal limit. At the same time, one might pursue other observational channels. The above criterion ensures that such an object will likely not possess a photon sphere, specifically rp = 3 M for the Regge-Wheeler potential, resulting in a distinct optical image compared to objects emitting the first type of echo. From an empirical perspective, we argue that this condition is feasible in realistic scenarios, and various compact objects may exhibit the second type of echo. For instance, neutron stars typically have a radius of 10-15 kilometers, while their Schwarzschild radius is about 3 kilometers [81]. This implies that neutron stars usually do not have a photon sphere and are subject to the second type of echo. On the other hand, as the density of the compact object increases further, it might become optically 'dark' due to the presence of a photon sphere while exhibiting the first type of echo. Furthermore, beyond the specific quasinormal oscillations of stars discussed in this work, we speculate that the second type of echo signal could have broader relevance in a universal context. For instance, on a larger scale, discontinuities arising from phenomena like Bondi accretion around compact objects could give rise to such echoes, which might be used to probe their presence. Similarly, these echoes might manifest in the gravitational wave spectrum of compact binaries, as surfaces and resulting discontinuities are inherent features of any compact object. In conclusion, our findings regarding the two types of echoes in compact stars enrich the theoretical framework for investigating the universe's most extreme gravitational environments. This potentially opens new pathways for testing the predictions of general relativity and exploring alternative theories of gravity, contributing to our ongoing quest to understand the cosmos.

## ACKNOWLEDGMENTS

We also gratefully acknowledge the financial support from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), A part of this work was developed under the project Institutos Nacionais de Ciências e Tecnologia - Física Nuclear e Aplicações (INCT FNA) Proc. No. 464898 2014-5. This research is also supported by the Center for Scien-/ / tific Computing (NCC GridUNESP) of São Paulo State University (UNESP). T. Zhu is supported / by the National Key Research and Development Program of China Grant No. 2020YFC2201503, and the Zhejiang Provincial Natural Science Foundation of China under Grant No. LR21A050001 and No. LY20A050002, the National Natural Science Foundation of China under Grant No. 12275238, and the Fundamental Research Funds for the Provincial Universities of Zhejiang in China under Grant No. RF-A2019015.

## Appendix A: The reflection amplitude at the point discontinuity

In this appendix, we derive the specific forms of the transit amplitudes utilized in Eqs. (21) and (23) in the main text.

To derive Eq. (21), one substitutes the asymptotic forms given by the first rows of Eqs. (9)

and (10) into Eq. (19). Subsequently, one matches the coe ffi cients of the wave functions e i ω x and e -i ω x between Eqs. (19) and (20) and finds which can be simplified to give

$$\frac { 1 + \widetilde { C } ( \omega ) B _ { \text{in} } } { e ^ { i \omega x _ { c } } } & = \frac { \widetilde { C } ( \omega ) B _ { \text{out} } } { \widetilde { \mathcal { R } } ( \omega ) e ^ { - i \omega x _ { c } } }, \quad & ( \text{A1} ) \\ \intertext { d to give } \quad \quad \quad$$

where the reflection and transmission amplitudes are defined by Eq. (12).

To derive Eq. (23), one utilizes the WKB approximation and the waveform on the r.h.s. of the truncation point x ≥ xc is given by

$$\widetilde { C } ( \omega ) = \frac { e ^ { - 2 i \omega _ { c } } \widetilde { \mathcal { T } } \text{BH} ( \omega ) \widetilde { \mathcal { R } } ( \omega ) } { 1 - e ^ { - 2 i \omega _ { c } } \widetilde { \mathcal { R } } \text{BH} ( \omega ) \widetilde { \mathcal { R } } ( \omega ) }, &$$

$$\tilde { f } ( \omega, x ) = C e ^ { i S ( x _ { 0 }, x ) } + D e ^ { - i S ( x _ { 0 }, x ) },$$

where, at the lowest-order approximation that su ffi ces for our case, we have

$$\lambda ^ { \nu r } \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \$$

where k x ( ) = √ ω 2 -V x ( ), and x 0 can be taken somewhat arbitrarily in the region where the WKB formula is relevant. Without loss of generosity, we assume x 0 = 0. While for x ≤ xc , we have

$$f & = e ^ { - i \omega x }.$$

The Wronskian vanishes at the point of discontinuity x = xc , which gives

By taking into account S ′ = k x ( ) and V x ( ) = Vc ∼ 0, we can approximate

$$\intertext { u o c o u n } \varrho = \kappa ( x ) \, \text{and} \, v \, ( x ) = v _ { c } \sim \mathfrak { u }, \, \text{we can app} \, \text{may} \, \text{made} \\ \frac { C } { D } = \frac { k ( x ) - \omega } { k ( x ) + \omega } e ^ { - 2 i S ( 0, x _ { c } ) } = \frac { \sqrt { \omega ^ { 2 } - V _ { c } } - \omega } { \sqrt { \omega ^ { 2 } - V _ { c } } + \omega } e ^ { - 2 i \omega x _ { c } }. \\ \text{ing against the form of Eq. (19), one finds}$$

$$\tan vanishes & \text{at the point of discontinuity } x = x _ { c }, \text{ which gives} \\ & i S ^ { \prime } ( x _ { 0 }, x ) \frac { C e ^ { i S ( 0, x ) } - D e ^ { - i S ( 0, x ) } } { C e ^ { i S ( 0, x ) } + D e ^ { - i S ( 0, x ) } } \Big | _ { x = x _ { c } } = - i \omega, & ( A 6 ) \\ \text{account} \, S ^ { \prime } = k ( x ) \text{ and } V ( x ) = V _ { c } \sim 0, \, \text{ we can approximate} \\ & \sim \quad. \, \cdots \, \sqrt { \frac { 1 } { 2 } }.$$

By comparing against the form of Eq. (19), one finds

$$\text{inst use norm or equal. } \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varproptor \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varproptore \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varpropto \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptodelta \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptide \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptore \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptcode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptoded \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptideo \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptude \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptede \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptoda \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptade \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptoder \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptodes \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproproptode \varproptode \varproproptode \varproptode \varproptode \varproproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \ varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \vproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varproptode \varpropro \varproptode \varproptode \varproptode \varproptode$$

which is the result used in Eq. (23).

[1] Virgo, LIGO Scientific, B. P. Abbott et al. , Phys. Rev. Lett. 116 , 061102 (2016), arXiv:1602.03837.

[2] Virgo, LIGO Scientific, B. P. Abbott et al. , Phys. Rev. Lett. 116 , 221101 (2016), arXiv:1602.03841, [Erratum: Phys. Rev. Lett.121,no.12,129902(2018)].

- [3] Virgo, LIGO Scientific, B. P. Abbott et al. , Phys. Rev. Lett. 116 , 241103 (2016), arXiv:1606.04855.
- [4] Virgo, LIGO Scientific, B. P. Abbott et al. , Phys. Rev. Lett. 119 , 141101 (2017), arXiv:1709.09660.
- [5] LISA, P. Amaro-Seoane et al. , (2017), arXiv:1702.00786.
- [6] TianQin, J. Luo et al. , Class. Quant. Grav. 33 , 035010 (2016), arXiv:1512.02076.
- [7] W.-R. Hu and Y.-L. Wu, Natl. Sci. Rev. 4 , 685 (2017).
- [8] M. Tinto and S. V. Dhurandhar, Living Rev. Rel. 17 , 6 (2014).
- [9] M. Tinto and S. V. Dhurandhar, Living Rev. Rel. 24 , 1 (2021).
- [10] P.-P. Wang, Y.-J. Tan, W.-L. Qian, and C.-G. Shao, Phys. Rev. D 104 , 082002 (2021), arXiv:2106.02236.
- [11] P.-P. Wang, Y.-J. Tan, W.-L. Qian, and C.-G. Shao, Phys. Rev. D 103 , 063021 (2021).
- [12] P.-P. Wang, Y.-J. Tan, W.-L. Qian, and C.-G. Shao, Phys. Rev. D 104 , 023002 (2021).
- [13] O. Dreyer et al. , Class. Quant. Grav. 21 , 787 (2004), arXiv:gr-qc 0309007. /
- [14] E. Berti, V. Cardoso, and C. M. Will, Phys. Rev. D73 , 064030 (2006), arXiv:gr-qc 0512160. /
- [15] M. Giesler, M. Isi, M. A. Scheel, and S. Teukolsky, Phys. Rev. X 9 , 041060 (2019), arXiv:1903.08284.
- [16] M. Cabero et al. , Phys. Rev. D 101 , 064044 (2020), arXiv:1911.01361.
- [17] A. Dhani, Phys. Rev. D 103 , 104048 (2021), arXiv:2010.08602.
- [18] H. Liu, C. Zhang, Y. Gong, B. Wang, and A. Wang, Phys. Rev. D102 , 124011 (2020), arXiv:2002.06360.
- [19] M. Visser, Phys. Rev. D46 , 2445 (1992), arXiv:hep-th 9203057. /
- [20] P. T. Leung, Y. T. Liu, W. M. Suen, C. Y. Tam, and K. Young, Phys. Rev. Lett. 78 , 2894 (1997), arXiv:gr-qc 9903031. /
- [21] P. T. Leung, Y. T. Liu, W. M. Suen, C. Y. Tam, and K. Young, Phys. Rev. D59 , 044034 (1999), arXiv:gr-qc 9903032. /
- [22] E. Barausse, V. Cardoso, and P. Pani, Phys. Rev. D89 , 104059 (2014), arXiv:1404.7149.
- [23] H.-P. Nollert, Class. Quant. Grav. 16 , R159 (1999).
- [24] E. Berti, V. Cardoso, and A. O. Starinets, Class. Quant. Grav. 26 , 163001 (2009), arXiv:0905.2975.
- [25] B. Wang, Braz. J. Phys. 35 , 1029 (2005), arXiv:gr-qc 0511133. /
- [26] J. D. Bekenstein, Phys. Rev. Lett. 28 , 452 (1972).
- [27] J. D. Bekenstein, Phys. Rev. D51 , R6608 (1995).
- [28] H.-P. Nollert, Phys. Rev. D53 , 4397 (1996), arXiv:gr-qc 9602032. /
- [29] H.-P. Nollert and R. H. Price, J. Math. Phys. 40 , 980 (1999), arXiv:gr-qc 9810074. /
- [30] W.-L. Qian, K. Lin, C.-Y. Shao, B. Wang, and R.-H. Yue, Phys. Rev. D103 , 024019 (2021), arXiv:2009.11627.
- [31] H.-P. Nollert, Phys. Rev. D47 , 5253 (1993).
- [32] L. Motl, Adv. Theor. Math. Phys. 6 , 1135 (2003), arXiv:gr-qc 0212096. /
- [33] J. L. Jaramillo, R. Panosso Macedo, and L. Al Sheikh, Phys. Rev. X 11 , 031003 (2021), arXiv:2004.06434.
- [34] J. L. Jaramillo, R. Panosso Macedo, and L. A. Sheikh, Phys. Rev. Lett. 128 , 211102 (2022), arXiv:2105.03451.
- [35] K. Destounis, R. P. Macedo, E. Berti, V. Cardoso, and J. L. Jaramillo, Phys. Rev. D 104 , 084091 (2021), arXiv:2107.09673.
- [36] A. Zengino˘ glu, Phys. Rev. D83 , 127502 (2011), arXiv:1102.2451.
- [37] M. H.-Y. Cheung, K. Destounis, R. P. Macedo, E. Berti, and V. Cardoso, Phys. Rev. Lett. 128 , 111103 (2022), arXiv:2111.05415.
- [38] V. Cardoso, E. Franzin, and P. Pani, Phys. Rev. Lett. 116 , 171101 (2016), arXiv:1602.07309, [Erratum: Phys. Rev. Lett.117,no.8,089902(2016)].

- [39] R. Dong and D. Stojkovic, Phys. Rev. D 103 , 024058 (2021), arXiv:2011.04032.
- [40] K. Chakravarti, R. Ghosh, and S. Sarkar, Phys. Rev. D 104 , 084049 (2021), arXiv:2108.02444.
- [41] K. Chakravarti, R. Ghosh, and S. Sarkar, Phys. Rev. D 105 , 044046 (2022), arXiv:2112.10109.
- [42] P. O. Mazur and E. Mottola, Proc. Nat. Acad. Sci. 101 , 9545 (2004), arXiv:gr-qc 0407075. /
- [43] M. Visser and D. L. Wiltshire, Class. Quant. Grav. 21 , 1135 (2004), arXiv:gr-qc 0310107. /
- [44] M. S. Morris and K. S. Thorne, Am. J. Phys. 56 , 395 (1988).
- [45] M. S. Morris, K. S. Thorne, and U. Yurtsever, Phys. Rev. Lett. 61 , 1446 (1988).
- [46] M. Visser, Phys. Rev. D39 , 3182 (1989), arXiv:0809.0907.
- [47] M. Visser, Nucl. Phys. B328 , 203 (1989), arXiv:0809.0927.
- [48] T. Damour and S. N. Solodukhin, Phys. Rev. D76 , 024016 (2007), arXiv:0704.2667.
- [49] S. Biswas, M. Rahman, and S. Chakraborty, Phys. Rev. D 106 , 124003 (2022), arXiv:2205.14743.
- [50] S. Biswas, C. Singha, and S. Chakraborty, Phys. Rev. D 109 , 064043 (2024), arXiv:2307.04836.
- [51] Z. Mark, A. Zimmerman, S. M. Du, and Y. Chen, Phys. Rev. D96 , 084002 (2017), arXiv:1706.06155.
- [52] P. Bueno, P. A. Cano, F. Goelen, T. Hertog, and B. Vercnocke, Phys. Rev. D97 , 024040 (2018), arXiv:1711.00391.
- [53] H. Liu et al. , Phys. Rev. D104 , 044012 (2021), arXiv:2104.11912.
- [54] S.-F. Shen, W.-L. Qian, K. Lin, C.-G. Shao, and Y. Pan, Class. Quant. Grav. 39 , 225004 (2022), arXiv:2203.14320.
- [55] L. Hui, D. Kabat, and S. S. C. Wong, JCAP 12 , 020 (2019), arXiv:1909.10382.
- [56] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astrophys. J. 462 , 563 (1996), arXiv:astro-ph 9508025. /
- [57] B. Moore, F. Governato, T. R. Quinn, J. Stadel, and G. Lake, Astrophys. J. Lett. 499 , L5 (1998), arXiv:astro-ph 9709051. /
- [58] O. Valenzuela et al. , Astrophys. J. 657 , 773 (2007), arXiv:astro-ph 0509644. /
- [59] B. Diemer and A. V. Kravtsov, Astrophys. J. 789 , 1 (2014), arXiv:1401.1216.
- [60] S. Adhikari, N. Dalal, and R. T. Chamberlain, JCAP 1411 , 019 (2014), arXiv:1409.4482.
- [61] E. W. Leaver, Phys. Rev. D34 , 384 (1986).
- [62] A. Anderson and R. H. Price, Phys. Rev. D 43 , 3147 (1991).
- [63] H.-P. Nollert and B. G. Schmidt, Phys. Rev. D45 , 2617 (1992).
- [64] V. P. Frolov and I. D. Novikov, Black Hole Physics: Basic Concepts and New Developments (Kluwer Academic, 1998).
- [65] S. Grozdanov, K. Schalm, and V. Scopelliti, Phys. Rev. Lett. 120 , 231601 (2018), arXiv:1710.00921.
- [66] M. Blake, R. A. Davison, and D. Vegh, JHEP 01 , 077 (2020), arXiv:1904.12883.
- [67] W.-L. Qian, K. Lin, J.-P. Wu, B. Wang, and R.-H. Yue, Eur. Phys. J. C80 , 959 (2020), arXiv:2006.07122.
- [68] W.-L. Qian, K. Lin, X.-M. Kuang, B. Wang, and R.-H. Yue, Eur. Phys. J. C 82 , 188 (2022), arXiv:2109.02844.
- [69] K. D. Kokkotas and B. F. Schutz, Mon. Not. Roy. Astron. Soc. 255 , 119 (1992).
- [70] K. D. Kokkotas, Pulsating relativistic stars, in Relativistic gravitation and gravitational radiation. Proceedings, School of Physics, Les Houches, France, September 26-October 6, 1995 , pp. 89-102, 1995, gr-qc 9603024. /
- [71] K. Tominaga, M. Saijo, and K.-i. Maeda, Phys. Rev. D60 , 024004 (1999), arXiv:gr-qc 9901040. /
- [72] V. Ferrari and K. D. Kokkotas, Phys. Rev. D62 , 107504 (2000), arXiv:gr-qc 0008057. /
- [73] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity (John Wiley and Sons, New York, 1972).
- [74] L. K. Tsui and P. T. Leung, Mon. Not. Roy. Astron. Soc. 357 , 1029 (2005), arXiv:gr-qc 0412024. /
- [75] C. Gundlach, R. H. Price, and J. Pullin, Phys. Rev. D49 , 883 (1994), arXiv:gr-qc 9307009. /

- [76] C. Gundlach, R. H. Price, and J. Pullin, Phys. Rev. D49 , 890 (1994), arXiv:gr-qc 9307010. /
- [77] J. S. F. Chan and R. B. Mann, Phys. Rev. D55 , 7546 (1997), arXiv:gr-qc 9612026. /
- [78] K. Lin and W.-L. Qian, Chin. Phys. C 47 , 085101 (2023), arXiv:2204.09531.
- [79] B. Cuadros-Melgar, J. de Oliveira, and C. E. Pellicer, Phys. Rev. D85 , 024014 (2012), arXiv:1110.4856.
- [80] K. Lin, W.-L. Qian, and A. B. Pavan, Phys. Rev. D94 , 064050 (2016), arXiv:1609.05963.
- [81] A. W. Steiner, J. M. Lattimer, and E. F. Brown, Astrophys. J. Lett. 765 , L5 (2013), arXiv:1205.6871.
Shui-Fa Shen Kai Lin Tao Zhu Yu-Peng Yan Cheng-Gang Shao Wei-Liang Qian
Two distinct types of echoes in compact objects
In the black hole perturbation theory framework, two different physical
pictures for echoes in compact objects have been proposed. The first mechanism
interprets echoes as repeated reflections of gravitational waves within a
potential well, where the echo period is defined by twice the distance related
to the spatial displacement operator that separates two local maxima of the
effective potential. The second mechanism associates echoes with a
discontinuity in the effective potential, potentially associated with specific
accretion processes, without necessarily introducing a second local maximum in
the effective potential. This discontinuity leads to echo signals that are
typically attenuated over time more quickly, with their period dictated by the
characteristics of the transfer amplitudes. In both scenarios, the echoes
correspond to a new category of quasinormal modes with minor real parts, with
their period connected to the spacing between successive modes in the frequency
domain. This work elaborates on a unified framework in compact stars that
encompasses both echo mechanisms. It suggests that these two types of echoes
derive from different physical origins and can be independently triggered. The
occurrence and interplay between these two types of echoes are demonstrated
through numerical simulations. %The observational relevance of this study is
also addressed.