2025-12-29 | | Total: 31
We investigate the observational tests of generalized mass-to-horizon entropic cosmology by incorporating large-scale structure growth data in addition to purely geometric probes. The theoretical framework is constructed from a generalized mass-to-horizon scaling relation, $M \propto L^n$, which implies a corresponding generalized entropic functional $S_n \propto L^{n+1}$. Within this setting, cosmic acceleration arises as an emergent phenomenon driven by an entropic force acting on the cosmological horizon. While earlier studies demonstrated that these entropic cosmologies can reproduce the background expansion history of the standard $Λ$CDM model, here we present a comprehensive observational analysis that jointly employs Pantheon+ Type Ia supernova data with SH0ES calibration, DESI DR2 baryon acoustic oscillation measurements, cosmic microwave background (CMB) distance priors, and a suite of cosmological structure growth observations. A Bayesian model comparison indicates that the entropic models are statistically preferred over the conventional $Λ$CDM scenario, thereby providing strong support for an entropic origin of the observed late-time cosmic acceleration in place of a fundamental cosmological constant.
A classical result by Penrose establishes that null geodesics generating a black hole event horizon can only intersect at their entrance to the horizon in ``crossover'' points. This points together with limit points of this set, namely caustics, form the so-called "crease set". Light rays enter into the horizon through the crease set, characterizing the latter as the birth of the horizon. A natural question in this context refers to the classification and stability of the structural possibilities of black hole crease sets. In this work we revisit the strategy adopted by Gadioux & Reall for such a classification in the setting of singularity theory in contact geometry. Specifically, in such contact geometry setting, the event horizon is identified as a component (not connected to null infinity) of a so-called ``BigFront''. The characterization of BigFronts as Legendrian projections of Legendrian submanifolds permits to classify the crease sets and ``cuspidal sets'' (or caustics in Penrose's terminology) by applying classical results established by V.I. Arnol'd. Here we refine the stability discussion presented by Gadioux & Reall of that connected component of the crease set that is not causally connected to null infinity and that constitutes the event horizon birth. In addition, we identify the existence of other components of the crease set that lie in the part of the BigFront that is causally connected to null infinity.
We study thermodynamic properties and the entanglement island of a black hole in asymptotically safe quantum gravity, analyzing key thermodynamic quantities such as the Hawking temperature, heat capacity, and entropy, as well as the mass-horizon radius relation. Unlike Schwarzschild black holes, the temperature decreases with mass near the evaporation endpoint, signaling a phase transition and possible stable remnant. The entanglement entropy of Hawking radiation is obtained both with and without island contributions. Without islands, the radiation entropy grows linearly indefinitely, leading to the information paradox. By including island contributions and extremizing the generalized entropy functional, we resolve this paradox. At late times, the radiation entropy saturates at the Bekenstein-Hawking entropy, confirming unitary evolution. From this, we derive the Page time and scrambling time by equating early- and late-time entanglement entropies. The result of this study establishes the finiteness of the radiation entropy and consistency with quantum mechanics.
In this article, we are devoted to discuss different compact stars admitting anisotropic interiors in a particular modified theory of gravity. For this purpose, a spherically symmetric metric is adopted to formulate the field equations corresponding to two different $f(\mathcal{R},\mathcal{T},\mathcal{Q})$ models, where $\mathcal{Q}=\mathcal{R}_{αγ}\mathcal{T}^{αγ}$. Since the field equations contain extra degrees of freedom, we choose Finch-Skea metric and MIT bag model equation of state to make them solvable. We also use matching conditions to calculate a constant triplet in the chosen ansatz. The resulting solutions are then graphically analyzed for particular values of the bag constant and model parameter in the interior of 4U 1820-30 compact star. The viability and stability of the modified models are also checked through certain tests. Further, we calculate the values of model parameter through the vanishing radial pressure constraint that correspond to the observed data (radii and masses) of eight different star candidates. Finally, we conclude that our models I and II are in well-agreement with the conditions needed for physically relevant interiors to exist.
A tertiary body near a coalescing binary can imprint its influence on the gravitational waves (GWs) emitted by that binary in the form of center-of-mass (CoM) acceleration. An example of such a scenario is a binary black hole (BBH) merging near a supermassive black hole, which is touted to occur frequently. The limited low-frequency sensitivity of current GW detectors makes it challenging to detect these effects, as the associated waveform phase remains elusive. However, next-generation (3G) detectors such as Cosmic Explorer (CE) and Einstein Telescope (ET), with improved sensitivity at lower frequencies, are expected to be capable of capturing such signatures. In our study, we focus on the stellar-mass BBHs and explore the parameter space where the CoM acceleration will play a dominant role affecting parameter inference of the binary. We demonstrate that an unaccounted CoM acceleration of a low-mass binary with a total mass of $5\, \rm{M}_{\odot} $ can lead to significant systematic biases, exceeding statistical errors in the estimation of the chirp mass and symmetric mass ratio when the CoM parameter $α$ is as small as $\sim 10^{-9}$ and $10^{-10}$ $\rm{s}^{-1}$ for CE and ET, respectively. We also find that asymmetric binaries are more susceptible to systematic bias when CoM acceleration is neglected. When the effect of CoM acceleration is included in the GW phase, then $α= 10^{-7} \rm s^{-1}$ can be constrained with a bound of $10^{-9} (10^{-11})\, \rm s^{-1}$ for CE (ET). Our study thus highlights the crucial implications of considering the presence of a tertiary body in the GW emitted by a stellar-mass BBH when observed in 3G detectors.
Mini-Extreme-Mass-Ratio Inspirals (mini-EMRIs), comprising a sub-solar exotic compact object (such as a primordial black hole or boson star) orbiting a much heavier stellar-origin or exotic compact object, represent key targets for ground-based gravitational-wave detectors to probe the early universe and the nature of dark matter. However, detecting such systems, which could spend hours to years in LIGO, Virgo and KAGRA data, poses a computational challenge to standard matched-filtering methods. However, semi-coherent methods are constrained by the quasi-monochromatic assumption, which restricts the coherence time to avoid spectral leakage caused by frequency evolution. In this work, we extend the development of our method, $Σ$Track, to the regime in which the quasi-monochromatic approximation is relaxed, in two ways. First, we establish an analytical model for the spectral leakage, extending the validity of conventional analyses beyond the quasi-monochromatic regime. Second, we propose the $ΣR$ statistic -- a novel detection metric formed by a weighted summation of power ratios -- which effectively recovers the signal energy dispersed across adjacent frequency bins. Building on this framework, we further introduce an innovative frequency-layered search strategy that dynamically optimizes the coherence time across the observation band. We benchmark our method against a globally optimized Hough transform pipeline using a fiducial mini-EMRI signal from a binary with masses $(1.5, 10^{-5})\,M_\odot$. The results demonstrate that our framework achieves an order-of-magnitude enhancement in the effective detection volume, significantly expanding the horizon for discovering mini-EMRIs and sub-solar exotic compact objects with ground-based gravitational wave detectors. This approach can be similarly applied to EMRI searches for future space-based gravitational wave detectors.
Mini-extreme-mass-ratio inspirals (mini-EMRIs), composed of a stellar-mass compact object and a much lighter companion, are promising sources of continuous gravitational waves in the frequency band of ground-based interferometers such as LIGO-Virgo-KAGRA. Such systems, consisting of sub-solar-mass compact objects, provide a unique probe of exotic compact objects, including primordial black holes. Detecting such long-lived signals, however, remains challenging. Here, we adapt standard methods used in searches for quasi-monochromatic signals to search for mini-EMRIs, and derive a statistical framework that explicitly handles spectral leakage. In particular, we introduce a new method that sums along the tracks in the time-frequency plane carved out by possible mini-EMRI signals, which we call $Σ$Track. This refinement establishes a general basis for analyzing long-duration transient signals with rapid frequency evolutions, regardless of the underlying mechanism for gravitational-wave emission. We also compute a new semi-analytic sensitivity estimate within our new statistical framework, which is valid under the assumption that the signal is weak with respect to the noise level. We then establish a statistic that quantifies how to discretize the search parameter space for our method, which works for mini-EMRIs, as well as arbitrary signal types. Our results provide a foundation for mini-EMRI searches and demonstrate the potential of current ground-based detectors to probe the existence of sub-solar-mass compact objects.
A unified equation is employed to analytically investigate the scattering of massless spin particles by a Schwarzschild-type medium black hole. It is found that for spin particles, curved spacetime induces an effective complex potential analogous to a Coulomb field. While the real part of this potential contributes a real logarithmic term to the phase, the imaginary part gives rise to a corresponding imaginary logarithmic term. Crucially, this imaginary term is precisely responsible for generating the correct asymptotic decay of the wave function. From this framework, a unified analytical expression for the differential cross section is derived, applicable to all particle types considered. Given the successful fabrication of a Schwarzschild-equivalent medium via transformation optics, our theoretical scattering predictions can be tested experimentally by transmitting plane electromagnetic waves through such a structure. Insights gained from these experiments could, in turn, shed light on the scattering of other massless fields (e.g., gravitational waves) by actual black holes.
This study investigates the quasinormal modes of Kerr-Newman-de Sitter black holes for massless spin particles using the unified equation. We derive analytical expressions for both the quasinormal mode frequencies and the radial wave functions. The frequencies are determined exclusively by the black hole parameters and the quantum numbers $n$ and $m$, while the radial wave functions also depend on the quantum number $k$, indicating a degeneracy in frequency. For identical quantum numbers, the frequency expression and the degree of degeneracy are the same for all massless spin particles, regardless of their specific properties. This implies that, through the observation of quasinormal modes, one can not only determine the black hole's parameters but also observe the phenomenon in which one type of particle reproduces the quasinormal mode of another. Our work thus provides a theoretical foundation for understanding this mimicking behavior as well.
We present a detailed evaluation of constrained minisuperspace path integrals in Jackiw-Teitelboim (JT) gravity and in biaxial Bianchi IX quantum cosmology, employing the Gelfand-Yaglom theorem to compute the relevant functional determinants. In both settings, integrating out the dilaton or a minisuperspace variable produces a functional delta that enforces the classical constraint equation, thereby localizing the remaining path integral onto classical configurations. The associated Jacobian, equivalently, the functional determinant of the operator obtained by linearizing the constraint about the classical solution, fixes the semiclassical prefactor and the correct measure. We evaluate this determinant exactly via the Gelfand-Yaglom method and obtain the fully normalized fixed-lapse propagators. We further extend the JT analysis to a quadratic dilaton potential $U(φ)=m^{2}φ^{2}$ and comment on the corresponding saddle-point structure of the lapse integral. Finally, we apply the same approach to Bianchi IX quantum cosmology and derive the fixed-lapse propagator, including its prefactor. Our results provide a systematic and broadly applicable prescription for treating constraint structures in gravitational path integrals using functional determinant techniques, with potential applications to a wider class of minisuperspace quantum cosmology and quantum gravity.
This short report is dedicated to the 40th anniversary of International Journal of Modern Physics A (IJMPA) and Modern Physics Letters A (MPLA). While the report is based on a series of papers[1-8], its content reflects my personal viewpoints. Therefore I am solely responsible for all the statements in the report. In this report we discuss conservation of energies in a curved spacetime including general relativity. We argue that the matter energy is not necessarily conserved in a curved space time due to a lack of time translational invariance, and adding energy of gravitational fields to recover the conservation law of energies fails due to Noether's 2nd theorem. We show that there exists conserved quantities associated with the matter
Fundamental physics today is best defined operationally: it is the program of identifying the microscopic degrees of freedom, symmetries, and dynamical laws that (i) reproduce the Standard Model (SM) of particle physics, General Relativity (GR), and the $Λ$CDM cosmological model in their regimes of validity, and (ii) explain the observed phenomena that these baseline theories do not account for (dark matter, neutrino masses, baryogenesis, dark energy), while resolving conceptual inconsistencies (quantum gravity, naturalness, the cosmological constant problem, the measurement problem in quantum theory, information in black holes) and providing predictive unification. This review first lays out the SM+GR+$Λ$CDM baseline, the best current evidence for its parameters, and the concrete anomalies and missing ingredients. It then surveys the most relevant theoretical directions (effective field theories; amplitude/positivity programs; lattice and many-body methods; symmetry-based model building; cosmological EFTs; quantum information approaches to QFT/gravitation) and the experimental/observational landscape, including ground and space platforms, astronomical messengers, and in-situ tests. Throughout we emphasize: (a) how each observable maps to energy scales and couplings; (b) the dominant statistical and systematic limitations; (c) the sensitivity required for decisive progress. A staged roadmap is given only after the technical review, organized by decision points and cross-checks rather than by specific projects.
We introduce GraviBERT, a novel deep learning framework for inference of gravitational-wave time series, which relies on an Inception-inspired multi-scale convolutional feature extractor combined with a transformer encoder and a suitable regression head. GraviBERT is trained in two stages: a BERT-style pretraining phase, in which the model learns to predict masked segments in feature space to capture universal patterns and physics, followed by supervised fine-tuning for accurate parameter estimation. This approach demonstrates impressive improvement across multiple metrics compared to training from scratch. On in-domain data, it reduces the mean absolute error for point-estimate parameter inference by up to $30\%$, and training convergence accelerates by up to a factor of six. Moreover, at low signal-to-noise ratio, the mean relative precision of the inferred masses and distances reaches the few-percent level, while the mean absolute error in the effective spin is about $10^{-3}$. For domain adaptation to new detector noise profiles, the pretrained model demonstrates remarkable efficiency, converging up to $15\times$ faster on small target datasets and reducing estimation errors by up to approximately $45\%$, indicating that it learns sufficient detector-agnostic representations. Cross-approximant transfer demonstrates comparable performance, achieving up to $44\%$ reductions in mean absolute error across all parameters and up to $15\times$ training speedups, with $R^2$ scores consistently exceeding 0.9 for mass parameters at signal-to-noise ratio 10, compared to 0.74 - 0.87 when training from scratch. Notably, GraviBERT works directly with noisy waveforms. The final regression head of the model can be adapted for a range of downstream tasks after pretraining, positioning it as a step towards foundation-style models in gravitational-wave and multi-messenger astronomy.
Motivated by the work of Scully \textit{et al.} [ \textcolor{blue}{Proc. Nat. Acad. Sci. 115, 8131 (2018)}] and Camblong \textit{et al.}[ \textcolor{blue}{Phys. Rev. D 102, 085010 (2020)}], we investigate horizon-brightened acceleration radiation (HBAR) for freely falling two-level atoms in the geometry of a Bardeen regular black hole. Building on the quantum-optics approach to acceleration radiation and its near-horizon conformal quantum mechanics (CQM) structure, we show that the dominant physics is again governed by an inverse-square potential in the radial Klein-Gordon equation, with an effective coupling fixed by the Bardeen surface gravity. Using geodesic expansions and a near-horizon CQM reduction of the scalar field, we derive the excitation probability for atoms falling through a Boulware-like vacuum in the presence of a stretched-horizon mirror. The resulting spectrum is Planckian in the mode frequency, with a temperature determined by the Bardeen Hawking temperature. We analyze how the regular core parameter controls the strength of the radiation and demonstrate that the excitation probability is strongly suppressed as the geometry approaches the extremal (cold remnant) limit. Numerical results illustrate the dependence of the spectrum on the Bardeen parameter and on the atomic transition frequency.
Signatures of massive particles during inflation are highly informative targets for cosmological experiments. With recent progress on both theoretical and observational frontiers, we have reached the point where these novel signals of primordial non-Gaussianities (PNG) can be systematically tested with increasingly precise data. In this paper, we present the results of improved CMB data analysis for cosmological collider signals using Planck CMB data. To set the stage, we first construct a set of simplified but characteristic collider templates which are accurate over a broad range of particle masses, spins and sound speeds. In order to break degeneracies with single-field PNG, we propose an orthogonalization scheme such that the collider templates are uncorrelated with the highly constrained equilateral and orthogonal shapes. On this basis, we deploy the Modal bispectrum estimator for the Planck analysis and perform a systematic scan of parameters to search for the most significant collider signal. The maximum signal-to-noise ratio is found to be $2.35σ$ for massive spin-0 exchange after taking into account the look-elsewhere effect. In addition, we cross-validate the Modal analysis with the CMB-BEST pipeline, which demonstrates the consistency of results across the benchmark examples of collider templates. Given the low signal-to-noise ratio regime we find at the current stage of PNG observations, we believe the orthogonalization procedure provides an optimized strategy for future tests of the cosmological collider with the ability to rule out single field inflation.
The excursion-set formalism enables us to infer the mass distribution of collapsed objects, such as primordial black holes (PBHs), by the language of stochastic processes. Within the framework, this article investigates how a smooth coarse-graining procedure affects the resulting PBH mass function. As a demonstrative example, we employ a Gaussian window function, for which the stochastic noise becomes fully correlated across scales. It is found that these correlated noises result in a mass function of PBHs, whose maximum and its neighbourhood are predominantly determined by the probability that the density contrast exceeds a given threshold at each mass scale. Our results clarify the role of noise correlations induced by smooth coarse-graining and highlight their importance in predicting the abundance of PBHs.
Orbital eccentricity remains one of the least accessible parameters in observations of binary black hole (BBH) systems, largely erased by gravitational radiation long before detection. We introduce a new method to recover this lost parameter by using a more accessible and routinely measurable quantity: spin-orbit misalignment. In isolated binary evolution, a natal kick from the second supernova both tilts the orbital plane and injects orbital eccentricity, forging a direct and quantifiable connection between spin-tilt and post-supernova eccentricity. By measuring this spin-tilt using gravitational waves, we can not only constrain the natal kick, but we can also reconstruct the binary's formation eccentricity. We apply this method to GW190412 and GW241011, assuming an isolated formation channel, and show how their eccentricity at formation can be constrained even in the absence of direct eccentricity measurements. As more advanced detectors come online, improved signal-to-noise ratios will tighten spin-tilt constraints, allowing more precise and reliable estimates of BBH formation eccentricity. Combining this method with multiband observations from LISA and next-generation (XG) detectors will allow us to recover not only eccentricity but also the binary's orbital separation and redshift at formation, offering a clearer picture of the birth environments of BBH systems and processes that drive their merger.
Auxiliary field techniques have recently gained interest in four-dimensional non-linear electrodynamics and two-dimensional integrable sigma models. In these settings, coupling a suitable ``seed'' theory to auxiliary fields provides a powerful mechanism to generate infinite families of models while preserving key dynamical properties, such as electromagnetic duality invariance in four dimensions and classical integrability in two dimensions. Deformations induced through auxiliary fields are closely related to $T\bar{T}$-like deformations and, in two dimensions, also to their higher-spin generalisations. In this paper, we analyse and clarify the relations between different auxiliary field formulations in two and four dimensions, showing how they are governed by Legendre transformations of the interaction functions combined with appropriate field redefinitions. In four-dimensional electrodynamics, we establish a correspondence between the auxiliary field model of Russo and Townsend and the Ivanov--Zupnik formalism. In two dimensions, we develop the analogue of the Ivanov--Zupnik $μ$-frame to deform Principal Chiral, symmetric-space, non-Abelian T-dual, and (bi-)Yang-Baxter sigma models. We discuss how integrability is preserved and use properties of the $μ$-frame to further extend known families of integrable deformations.
Neutron stars such as pulsars and magnetars lose angular momentum primarily through electromagnetic dipole radiation, gravitational waves, $r$-mode oscillation, and also affected by fallback accretion processes. However, anomalous spin variations, particularly sudden enhanced spin-down rates, indicate additional spin-down mechanisms. We propose superradiant spin-down as a potential explanation for these events. By modelling the interplay between conventional and superradiant spin-down channels, we evaluate their impact on neutron star rotational evolution. We also discuss gravitational-wave emission produced by quadrupole deformation, $r$-mode oscillations, and axion-induced bosonic clouds around an isolated neutron star, highlighting their potential as distinct multimessenger probes in upcoming detectors.
We study the partition function associated with the democratic formulation of M-theory, focusing on its global definition and quantization properties. Using a path integral representation that makes manifest the underlying cohomological structure, we analyze the coupled system of form fields $C_4$ and $C_7$ and its associated gauge transformations. We show that the resulting description is naturally captured by a non-linear differential cocycle, reflecting the presence of a quadratic coupling between electric and magnetic degrees of freedom. This framework provides a transparent characterization of the global structure of the theory, clarifies the role of higher-form gauge symmetry, and allows for a consistent definition of the partition function in terms of higher-dimensional auxiliary manifolds.
Gravitational waves (GWs) from compact binary coalescences are standard sirens that provide a direct measure of the source's luminosity distance, enabling an independent measurement of the Hubble constant (H0). While a bright siren -- a GW event with an identified electromagnetic (EM) counterpart -- provided the first such constraint, most detections, currently dominated by black hole mergers, lack EM signatures. A measurement of H0 is still possible with these dark sirens by statistically associating GW events with galaxies in existing catalogs based on the sky localization. In this work, we explore the potential of two subsets of sirens: rare golden dark sirens, for which a single galaxy dominates the H0 posterior, and silver dark sirens, which are far more common but have a larger set of plausible host galaxies. Using the fifth internal data release of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), we assess the suitability of the Visible Integral-field Replicable Unit Spectrograph (VIRUS) for spectroscopic follow-up of dark sirens. VIRUS exposures of the standard HETDEX depth provide precise redshifts and exquisite completeness within z = 0.2. After a single year of observations with the upgraded LIGO-A# network, the combined sample of golden and silver dark sirens with z < 0.2 and follow-up VIRUS observations can potentially yield a few-percent constraint on H0. Our predictions suggest that spectroscopic redshift surveys such as HETDEX can play a key role in realizing high-precision cosmology with dark sirens in the near future. Standard-siren distance measurements offer a critical, fully independent path to the local value of H0 to resolve the Hubble tension.
The Moon has been regarded as a natural Weber bar capable of amplifying gravitational waves (GWs) for detecting events across a wide range of frequencies. However, accurately determining the amplification effects remains challenging due to the absence of 3D numerical simulation methods. In this study, we develop a high-order 3D finite element method (spectral-element method, SEM) to numerically simulate the lunar response to GWs below 20 mHz. We verify the accuracy of our method by comparing the resonant peaks of our results with those from semi-analytical solutions and find that the frequency deviation is less than 3% for the first peak at about 1 mHz and less than 0.8% for the subsequent peaks up to 10 mHz. Using this method, we evaluate the amplification of GW signals due to 3D topographic effects of the Moon, and we find enhancements at a series of specific frequency components. These results highlight the non-negligible effect of surface topography on the lunar response to GWs, as a fundamental factor that holds significant implications across both global and regional analyses. Our work paves the way for a comprehensive evaluation of the Moon's resonant response to GWs, helpful for the strategic planning of lunar GW detections.
Reheating is essential for transforming the cold, vacuum dominated Universe at the end of inflation into the hot thermal bath required by the Standard Model. In many well motivated inflationary models, however, the inflaton has no direct couplings to other fields, raising the question of how the Universe becomes repopulated with particles. We address this question within the framework of geometric reheating, where energy transfer occurs purely through gravitational effects. Focusing on a Higgs inflationary scenario with a non-minimal curvature coupling $ξφ^2 R$, we derive the post-inflationary dynamics and compute particle production using the Bogoliubov formalism. We show that the rapid, oscillatory evolution of the curvature scalar after inflaton acts as a time dependent gravitational pump, creating scalar spectator particles even in the absence of explicit interactions. This curvature driven production mechanism provides a natural and efficient route to reheating, demonstrating that gravity alone can initiate the standard thermal history and bridge inflation with radiation domination in minimal, coupling free models of the early Universe.
We establish an exact duality between the extended thermodynamics of five-dimensional charged Gauss-Bonnet AdS black holes and the thermodynamic framework of the dual boundary conformal field theory (CFT). The thermodynamics of the dual CFT involves two central charges originating from the trace anomaly. We demonstrate a precise correspondence between the extended first laws on the bulk and boundary sides. Moreover, the topological charges of the CFT thermodynamics, associated with the phase transition and critical point, coincide with those of the corresponding bulk black hole.
The region near a black hole horizon may be modified by quantum gravity effects that resolve the singularity. Such geometry may be represented by an exotic compact object. Because the horizon is enclosed by a photon sphere, it is difficult to probe this region directly. In this paper, we develop a method to study the region inside the photon sphere by applying the AdS/CFT correspondence. We extract signatures of the modified geometry from the retarded Green functions of the dual conformal field theory. The retarded Green functions can be computed from bulk wave functions of scalar field. We show that exotic compact objects leave two characteristic imprints: bulk-cone singularities and echoes. The bulk-cone singularities correspond to null geodesics in the bulk, allowing us to detect null trajectories that are specific to exotic compact objects. The echoes arise from wave modes trapped inside the photon sphere, and thus signal the absence of a horizon. As concrete examples, we study AdS gravastar and AdS wormhole. We compute the corresponding bulk wave functions both via the WKB approximation and through numerical analysis and observe the bulk-cone singularities and echoes explicitly.