Implications of a Simpson–Visser solution in Verlinde’s framework

This study focuses on investigating a regular black hole within the framework of Verlinde’s emergent gravity. In particular, we explore the main aspects of the modified Simpson–Visser solution. Our analysis reveals the presence of a unique physical event horizon under certain conditions. More-over, we study the thermodynamic properties, including the Hawking temperature, the entropy, and the heat capacity. Based on these quantities, our results indicate several phase transitions. Geodesic trajectories for photon–like particles, encompassing photon spheres and the formation of black hole shadows, are also calculated to comprehend the behavior of light in the vicinity of the black hole. Additionally, we also provide the calculation of the time delay and the deflection angle. Corroborating our results, we include an additional application in the context of high—energy astrophysical phenomena: neutrino energy deposition. Finally, we investigate the quasinormal modes using third–order WKB approximation.


I. INTRODUCTION
General Relativity, while a remarkable theory, is acknowledged to be incomplete when it comes to describing the behavior of the universe at both the classical and quantum levels.
One significant issue arises from the existence of singularities, which are problematic within this framework.Classical Einstein's theory of gravity encounters challenges when dealing with unavoidable singularities found in solutions like Schwarzschild, Reisner-Nordström, and Kerr metrics, which exhibit these peculiar features within their interiors.The scientific consensus recognizes the need for modifications to general relativity in regions where spacetime curvature becomes significantly high.These modifications are crucial for a more comprehensive understanding of gravity in extreme conditions.
In addition to addressing the presence of singularities, the modification of general relativity is essential for achieving a theory that is ultraviolet (UV) complete.Various proposals have emerged to accomplish this necessary modification.Extensive research has shown that incorporating higher-order curvature terms and terms involving higher derivatives can significantly improve the UV properties of Einstein's gravity [1][2][3][4].
In the absence of a specific theory, exploring potential modifications that could arise when gravity achieves UV completeness provides valuable insights.Such investigations are particularly informative when certain "natural" assumptions about the properties of a comprehensive theory are taken into account.In this context, our focus centers on the exploration of regular (non-singular) models of black holes, aiming to investigate black hole metrics that lack curvature singularities.The pioneering work of Bardeen [25] introduced the concept of a non-singular black hole, where the singularity was replaced by a charged matter core resulting from the collapse of charged matter.Additionally, a variety of models depicting such a feature have been proposed and discussed, including neutral, charged, and rotating configurations .
Nevertheless, a comprehensive grasp of gravitational waves and their properties is indispensable in exploring a myriad of physical phenomena, ranging from cosmological events in the primordial universe to astrophysical processes such as the evolution of stellar oscillations [47][48][49] and binary systems [50][51][52][53].These waves exhibit a diverse spectrum of intensities and distinct characteristic modes, with their spectral traits profoundly influenced by the underlying phenomena that engender them [54].When matter undergoes gravitational collapse, giving rise to the formation of a black hole, it enters a perturbed state, emitting radiation that encompasses an array of discrete frequencies unrelated to the collapse process itself [55].Termed quasinormal modes, these perturbations exhibit distinct frequencies that define their unique nature [56,57].
Significant progress has been made in the field of gravitational wave detection, enabling the identification of waves emitted by various physical phenomena [83][84][85][86].Ground-based interferometers, including VIRGO, LIGO, TAMA-300, and EO-600, have played a pivotal role in these detections [87][88][89][90].Over time, these detectors have substantially improved their precision, approaching a level of genuine sensitivity [91].These aspects gained from these advancements have provided valuable knowledge about the nature of astrophysical entities, encompassing boson and neutron stars for instance.
The detection of gravitational waves has far-reaching consequences for the study of black holes.By observing the emitted gravitational radiation, it becomes possible to directly confirm the existence of perturbed black holes [92].The pioneering work of Regge and Wheeler focused on exploring the stability of Schwarzschild black holes, laying the foundation for investigations into black hole perturbations [93].Subsequently, Zerilli made seminal contributions to the study of perturbations, significantly advancing our understanding in this field [94,95].
In recent years, there has been a significant surge of interest in the study of gravitational solutions involving scalar fields due to their remarkable characteristics.Notably, the behavior of black holes with nontrivial scalar fields has challenged the well-established "no-hair theorem" [96].This departure from the expected behavior has opened up exciting avenues of research, delving into the existence of long-lived scalar field patterns [97], the exploration of exotic astrophysical scenarios such as gravastars [98][99][100], and the formation of boson stars [101][102][103].Moreover, considering Klein-Gordon scalar fields on curved backgrounds has revealed a plethora of fascinating phenomena, including the intriguing concept of black hole bombs [104][105][106], where the scalar field can enhance the extraction of energy from the black hole.Additionally, the phenomenon of superradiance [107] arises when scalar fields interact with rotating black holes, leading to amplification effects.In Verlinde's theory [108], dark matter is proposed to emerge as a consequence of gravity arising from the distribution of baryonic matter.According to this theory, an additional gravitational effect is postulated due to the volume law contribution to entropy associated with positive dark energy.The hypothesis posits that the distribution of baryonic matter decreases the universe's overall entropy, triggering an elastic reaction in the fundamental microscopic system.This reaction gives rise to an extra gravitational force, commonly labeled as the dark matter effect, which is inherent to gravity's essence.Verlinde's theory offers a comprehensive model to interpret numerous observational phenomena, such as the rotation curves of galaxies.It proposes that the interplay between the distribution of baryonic matter and the perceived dark matter can elucidate the consistently observed flat rotation curves in galaxies.
In this study, our focus lies in examining a regular black hole within the framework of Verlinde's emergent gravity.Specifically, we place emphasis on investigating the modified Simpson-Visser solution.Through our analysis, we reveal the presence of a single physical event horizon.By studying the Hawking temperature and heat capacity, we unveil the existence of phase transitions.Furthermore, we provide the calculation of geodesic trajectories for photon-like particles, including critical orbits known as photon spheres.Additionally, our study includes calculations for both time delay and deflection angle.Moreover, to substantiate our findings, we introduce an additional application within the realm of high-energy astrophysical phenomena: the deposition of neutrino energy.investigate the quasinormal modes employing third-order WKB approximations.

II. REGULAR BLACK HOLE IN VERLINDE'S GRAVITY
In our manuscript, our main focus is to correlate baryonic matter and apparent dark matter in order to develop black hole solutions within the framework of a theory VEG (Verlinde's emergent gravity) [109].In other words, our goal is to investigate the effects arising from the existence of apparent dark matter on the spacetime geometry in the subsequent sections.
Verlinde says that in the context of spherical symmetry, a connection exists between the quantities of apparent dark matter M D (r) and baryonic matter M B (r) [109].This relationship can be mathematically written as follows: In this context, we have introduced a constant denoted by a 0 .In our investigation, we focus on the simplest scenario involving spherically symmetric black hole solutions.These ones are described by the line element: where f (r) = −g 00 = g −1 11 = 1 − 2m(r)/r, g 22 = r 2 , g 33 = r 2 sin 2 θ, and m(r), being given by [110] Building upon the recent findings presented in Ref. [109], we aim to derive a regular black hole solution.Subsequently, we will explore the impact of dark matter on the geometry of the spacetime in the upcoming sections.It is important to highlight that in Ref. [111] the authors revisited the general properties of regular models with a dS core in light of Verlinde's idea.Additionally, we shall analyze the discernible features arising from such a spacetime, which include the event horizon, thermodynamics, geodesics, shadows, and the characteristics of quasinormal modes.

III. THERMODYNAMICS
We proceed by following the approach presented in reference [112], where the metric describing the geometry of the spacetime is given by: with If we consider that Q → 0, the following solution gives rise to [109]: Above solution represents a generalization of the black-bounce spacetime geometry examined in reference [112].However, a notable distinction arises from the presence of apparent dark matter, which renders the spacetime non-asymptotically flat.
Despite the presence of six roots in the aforementioned expression, it is important to note that only one of these roots corresponds to a physical horizon (labeled as r + ): where, and O(M 5 ) and O(M 6 ) are higher-order mass terms, i.e., fifth-and sixth-orders respectively.It is important to mention that all terms (even the higher order ones) will be taken into account in our calculations.
In order to get a better comprehension of Eq. ( 6), we provide Fig. 1.To the left, we showcase f (r) across different values of a with M = 0.1, and to the right, we present f (r) for varying M with a = 0.1.Here, the horizon r + has a particularity, i.e., for possessing real positive defined values, the mass should should satisfy the following constraint: M > 0.14699 (when a = 1).Notice that such a condition will limit the whole thermodynamic properties of system from below.More so, a point of singularity is encountered here when M = 0.25 (for a = 1).All these features as well as the comparison with the Schwarzschild case are displayed in Fig. 2.
Upon obtaining the event horizon, we can now accurately examine the modified Hawking temperature as follows Figure 3 illustrates the behavior of the modified Hawking temperature as a function of mass.Remarkably, our results suggest the occurrence of a phase transition when M ≈ 0.1568.These remarks can be directly attributed to the influence of dark matter effects on For our purpose, we shall consider M + only.Among other reasons, when we consider a = 1 and r runs, the unique solution which gives rise to a real positive defined values of mass is M + .In this sense, we can perform the calculation of the Hawking temperature via first law of thermodynamics, T , as follows Notice that Eqs. ( 14) and ( 12) did not match to each other.In other words, it is clearly seen that the first law of thermodynamics did not give us the correct expression.In general lines, this is because regular black holes do not obey the so-called Bekenstein-Hawking area law.For the sake of overcoming this situation, a correction must be implemented in Eq. ( 14) [113].Instead, the corrected temperature is written as [113] Υ(r + , a)dM In this context,

≁
T represents the corrected version of the Hawking temperature obtained through the application of the first law of thermodynamics, with S denoting the entropy.
Notably, the function Υ(r + , a), intricately dependent on the terms of the mass function, not only governs the first law for regular black holes within the specific scenario under consideration but also assumes a crucial role in formulating the first law for diverse classes of regular black holes, as argued in [41].
As highlighted in Ref. [113], the general expression for Υ(r + , a) is articulated as follows: Here, the notation T 0 0 refers to the stress-energy component corresponding to energy density.
IV. GEODESIC TRAJECTORIES The study of particle motion within the realm of emergent gravity has garnered significant interest owing to its profound implications from a theoretical viewpoint [82,108,136].
Of particular significance is comprehending the geodesic characteristics of Simpson-Visser black holes, as it holds crucial relevance in understanding various astrophysical phenomena associated with these entities, including the nature of accretion disks and shadows for instance.In essence, our attention is directed towards comprehensively investigating the behavior dictated by the geodesic equation.To achieve this goal, we write where s is an arbitrary affine parameter.This investigation leads to the four coupled partial differential equations, which can be expressed as follows: As evident from the derived geodesic equations, we encounter four coupled partial differential equations.These ones serve as the foundation for determining the behavior of photon particles.Consequently, we perform ray tracing simulations to visualize the paths of light, as depicted in Fig. 6.

V. CRITICAL ORBITS AND SHADOWS
A comprehensive understanding of particle dynamics and the behavior of light rays in the vicinity of black hole structures necessitates a profound knowledge of critical orbits.These orbits hold a pivotal significance in elucidating the properties of spacetime influenced by dark matter effects in our specific context.
To achieve a more profound comprehension of the photon sphere's influence, often referred to as the critical orbit, in our black hole scenario, we will utilize the Lagrangian method to calculate null geodesics.This approach provides a clearer and more accessible understanding for readers compared to the utilization of the previously presented geodesic equation.Through this analysis, we aim to explore how the black hole mass impacts the photon sphere, shedding light on the gravitational effects inherent in the Simpson-Visser solution within Verlinde's emergent gravity framework.In this sense, we write: Upon considering a fixed angle of θ = π/2, the aforementioned expression undergoes simplification, resulting in: with L is the angular momentum and E being the energy.Next, Eq. ( 29) reads, with is the effective potential.In order to ascertain the critical radius, we must solve the equation ∂ ≁ V /∂r = 0. Here, considering a small, there exist two physical solutions for this equation (when the following range is considered for the values of M : [0.1, 0.2]), i.e., two photon spheres.These ones are represented by r c − (the inner photon sphere) and r c + (the outer photon sphere).It is worth mentioning that, although we have considered small values for a, all outputs will numerically be displayed.
In order to gain a deeper understanding of the behavior of r c − and r c + , we present Tables II and III.Examining the left side of Table II, we observe that the inner photon sphere expands with increasing a (for M = 0.1).On the right side of the same table, we note that there exists an augmented radius for the critical orbit r c − when the mass M increases (for a = 0.1).Turning our attention to Table III, it illustrates the values of r c + for different parameters M and a.Here, the increase of a and M corresponds to a higher value of r c + .
Recently in the literature, Ref. [82] addressed a similar study also in this context of Verlinde's emergent gravity.Additionally, it is important to mention that the appearance of two photon spheres instead was recently reported within the context Simpson-Visser solution [137,138] and others [139].
Furthermore, the study of shadows in the context of emergent gravity and black hole structures holds immense significance as it offers a remarkable perspective of understanding where K is the Carter constant.After some algebraic manipulations, we get In our quest to determine the radius of the shadow, we will employ the celestial coordinates α and β [57,140,141] in the following manner: α = −ξ and β = ± √ η.Utilizing these coordinates, we obtain the shadow radius: In Fig. 7, the illustration features circles that represent the shadows cast by varying mass M and a associated with the photon sphere r c − .To elaborate, on the left side, the sequence begins with the outer radius at M = 0.10, followed by M = 0.11, M = 0.12, and culminating with the inner radius at M = 0.13, all under the condition a = 0.10.Here, it can be observed that as the mass M increases, the shadow radius decreases.This interesting phenomenon illustrates the gravitational effect of the dark matter on the surrounding space, causing a noticeable reduction in the apparent size of the shadow when M increases.Moving to the right side, we examine different values of a for M = 0.10, ranging from the innermost at a = 0.10 to the outermost at a = 0.13.In other words, the considered values for this latter scenario encompass a = 0.10, a = 0.11, a = 0.12, and a = 0.13.Naturally, it is inherently verifiable that with the augmentation of the parameter a, there is a concomitant increase in the magnitude of the cast shadows.
In Fig. 8, the diagram depicts circles symbolizing shadows cast by varying mass M and a associated with the photon sphere r c + .On the left side, the sequence initiates with the inner radius at M = 0.10, followed by M = 0.11, and concluding with the outer radius at M = 0.12, all within the confines of the condition a = 0.1.This noteworthy phenomenon underscores the gravitational impact of dark matter on the surrounding space, inducing a perceptible increment in the apparent size of the shadow as M increases for a fixed value of a, i.e., a = 0.10, and for the following range of M : [0.10, 0.12].
Transitioning to the right side, we scrutinize various values of a for M = 1, spanning from the outermost at a = 0.10 to the innermost at a = 0.12.Specifically, the values considered in this latter scenario encompass a = 0.10, a = 0.11, and a = 0.12.Inevitably, it is demonstrably evident that an increase in the parameter a is intricately associated with a simultaneous reduction in the magnitude of the cast shadows for R + .These results are in contrast to those ones encountered in the analysis of R − .
Furthermore, the quantitative values of R ± are displayed in Tabs.IV and V for a better comprehension to the reader.

VI. TIME DELAY AND DEFLECTION ANGLE
The examination of time delay in dark matter scenarios holds profound significance within theoretical physics.This exploration provides profound insights into the fundamental aspects of a Simpson-Visser solution in Verlinde gravity.To determine the time delay, we employ the Lagrangian for null geodesics, enabling an accurate quantification of the temporal delay experienced by particles.This sheds light on the intricate dynamics of such a spacetime, as expressed by the equation: Also, knowing that dr/dt = 0 at r = r min , we derive: Thereby, the time delay can be computed as Furthermore, the bending of light as it traverses the curves of spacetime stands as a fundamental and intriguing phenomenon, serving as a crucial tool in the field of scientific inquiry.Light, exhibiting an elegant conformity to the contours of spacetime, reveals a captivating interplay that imparts invaluable understanding into the intricate physics governing gravitational sources [142][143][144][145][146][147].
To quantitatively describe this deflection angle, we employ a well-established formula [148,149] expressed as: to the innermost at a = 0.12.More explicitly, the values considered for this latter scenario encompass a = 0.10, a = 0.11, and a = 0.12.
In our analysis, we use the symbol α to denote the deflection angle.Our approach entails the initial determination of the minimum radius, r min , along the trajectory corresponding to each specific a and M .Particularly, in Tab.VI, we display the deflection angle α, considering different values of a with a fixed M .Here, we notice a very small difference between them.Such an aspect is expected, since we have been dealing with small a.

VII. NEUTRINO ENERGY DEPOSITION IN HIGH-ENERGY ASTROPHYSI-CAL PHENOMENA
The study is focused on examining the energy deposition rate resulting from the νν − → e + e − process, with the aim of elucidating its relevance to gamma ray burst emissions.The investigated scenario pertains to the concluding phase of neutron star mergers, conceptualized as a black hole accompanied by an accretion disk.
Salmonson and Wilson, as highlighted in Ref. [150,151], pioneered the exploration of the impacts within strong gravitational field regimes.Their groundbreaking work revealed that, within a Schwarzschild spacetime and for neutrinos emitted from the central core, the efficiency of the annihilation process νν − → e + e − undergoes substantial amplification-approximately 30 times greater than its Newtonian counterpart-especially in the context of collapsing neutron stars.
Building upon this foundation, subsequent studies [152,153] further examined the ramifications of general relativity on neutrino pair annihilation.This exploration extended to regions near the neutrinosphere and in proximity to a thin accretion disk (assuming an isothermal profile), with the gravitational background characterized by both Schwarzschild and Kerr geometries.
Our analysis centers around a black hole (BH) encircled by a slender accretion disk that emits neutrinos, as elaborated in [153].We concentrate on an idealized model, free from reliance on the specifics of disk formation and excluding considerations of self-gravitational effects.This disk is characterized by well-defined inner and outer edges, represented by radii denoted as R in and R out , respectively [154].
The Hamiltonian proves instrumental in analyzing the trajectory of a test particle within spacetime.It facilitates the computation of crucial parameters, including the energy and angular momentum of the test particle, and allows for the derivation of its equations of motion.In the scenario of a test particle traversing a curved background, the Hamiltonian takes the form Our emphasis is on determining the rate of energy deposition in close proximity to the axis, which is perpendicular to the disk, particularly at θ = 0 o .To evaluate the energy emitted within a half cone with an angular extent of approximately ∆θ ≈ 10 o , we must consider the scalar product of the momenta of a neutrino and an antineutrino at θ = 0 o .This scalar product can be expressed as Here, the term E ν is defined as the energy of the neutrino, calculated as E ν 0 / √ g 00 .Also, E ν 0 is defined as the observer energy of the neutrino measured at infinity and In this context, ρ ν is precisely defined as the ratio of the angular momentum L ν to the observed energy E 0ν .Moreover, considering geometric factors, there exist both a minimum and maximum value, denoted as θ m and θ M respectively, for a neutrino originating from , where r c denotes the radius of the photosphere.Furthermore, it can be shown that the following relationship holds [153]: where, r 0 represents the closest distance between the particle and the center before reaching θ = 0.The ultimate component is the trajectory equation, written as [154] ˆdr The equation above considers the emission of neutrinos from the position (R, π/2), where R ranges between R in and R out .These neutrinos then travel to the position (r, 0).As a result, the energy deposition rate resulting from neutrino pair annihilation is elaborated upon in [153] dE where k stands for the Boltzmann constant, G F represents the Fermi constant, T eff (2r c ) denotes the effective temperature at a radius of 2r c , and ζ(s) is the Riemann zeta function as follows After an algebraic manipulation, we get where T eff denotes the effective temperature as measured by a local observer and T 0 represents the temperature observed at infinity, having its general form is given by we do not account for the reabsorption of deposited energy by the black hole.Consequently, our focus centers on a scenario featuring a straightforward temperature gradient [153] T Moreover, the assumptions pertaining to temperature values and the configuration of the gradient model align with recent discoveries in neutrino-cooled accretion disk models, as shown in [155][156][157].
Anticipated within the theoretical framework is an effective maximum temperature, denoted as T eff , typically falling in the order of O(10 MeV).This magnitude assumes pivotal significance in achieving the observed neutrino disk luminosity, thereby rendering the disk luminosity comparatively invariant across diverse models.Given our non-engagement in numerical simulations, we adopt the assumption T eff ∼ O(10 MeV) to ensure a standardized basis for comparing the impacts of distinct gravitational models under identical conditions.
It is imperative to underscore that despite these theoretical suppositions, the exact temperature profile can only be ascertained through a comprehensive disk simulation originating from the merging of neutron stars with a meticulously defined geometry.Here, we set Furthermore, G(r) plays a pivotal role in computing the energy deposition rate and, consequently, in determining the energy deposition for a gamma ray burst event.
We calculate such an energy deposition rate within an infinitesimal angle dθ, considering a characteristic angle of 10 • degrees and a temperature of 10 MeV [153,154] dE 0 dt ≃ 4.41 × 10 48 ∆θ 10 VIII.THE QUASINORMAL MODES Throughout the ringdown phase, a remarkable phenomenon known as quasinormal modes emerges, displaying distinct oscillation patterns that remain not affected by the initial perturbations.Indeed, these modes manifest the intrinsic characteristics of the system and originate from the innate oscillations of spacetime, independent of specific initial conditions.
Unlike normal modes, which pertain to closed systems, quasinormal modes are associated with open systems.As a result, these modes dissipate energy gradually through the emission of gravitational waves.Mathematically, they can be characterized as poles of the complex Green function.
For determining their frequencies, one needs to find solutions to the wave equation within a system governed by a background metric g µν .However, acquiring analytical solutions for these modes is often a challenging procedure.
In the scientific literature, various techniques have been addressed to obtain solutions for these modes.Among them, the WKB (Wentzel-Kramers-Brillouin) approach stands out as one of the most widely used methods.Its development can be traced back to the groundbreaking work of Will and Iyer [158,159], and subsequent advancements up to the sixth order were made by Konoplya [160].For our specific calculations, we focus on analyzing perturbations using the scalar field, which involves considering the Klein-Gordon equation within the context of a curved spacetime While the exploration of backreaction effects in this particular scenario is intriguing, this manuscript does not provide this aspect and instead places its emphasis on other aspects.
Specifically, our primary focus revolves around studying the scalar field as a small perturbation.Furthermore, the presence of spherical symmetry allows us to take advantage of this opportunity to decompose the scalar field in a specific manner, as elaborated below: where we express the spherical harmonics as Y lm (θ, φ).Also, we can substitute the decomposition of the scalar field, as shown in Eq. ( 52), into Eq.( 51).This substitution transforms the equation into a Schrödinger-like form, endowed with wave-like properties, making it highly suitable for our analysis The potential V ef f is widely recognized as the Regge-Wheeler potential or the effective potential, carrying crucial information about the black hole's geometry.Additionally, we introduce the tortoise coordinate r * , which spans the entire spacetime as r * → ±∞.It is defined as dr * = [1/f (r) 2 ]dr.After some algebraic manipulations, the effective potential reads: Figure 9 illustrates the effective potential V ef f as a function of the tortoise coordinate r * for a particular value of l.Here, our primary objective is to derive stationary solutions for the system, achieved by assuming that Ψ(t, r) can be expressed as Ψ(t, r) = e −iωt ψ(r), where ω denotes the frequency.
This assumption allows us to conveniently separate the time-independent component of Eq. ( 53) by employing the following approach: To successfully solve Eq. ( 55), careful consideration of the appropriate boundary conditions becomes crucial.In our specific case, the solutions that satisfy the necessary conditions are characterized by exhibiting purely ingoing behavior near the horizon l (ω)e +iωr * (r * → +∞).
In our analysis, the complex constants C l (ω), A Also, it is worthy to be mentioned that the spectrum of quasinormal modes is determined by the eigenvalues of Eq. (55).To analyze these frequencies, we adopt the WKB method, i.e., a semi-analytical technique that draws parallels with quantum mechanics.This approach empowers us to gain valuable insights into the behavior of quasinormal modes near the event horizon and at infinity, unraveling essential aspects of black hole dynamics and gravitational wave phenomena.
The WKB approximation, first introduced by Schutz and Will [161], has become a valuable tool for computing quasinormal modes in the context of particle scattering around black holes.Over time, this method has seen further refinements, notably by Konoplya [160,162].However, it is essential to note that the validity of this approach hinges on the potential exhibiting a barrier-like shape, approaching constant values as r * → ±∞.By fitting the power series of the solution near the turning points of the maximum potential, the quasinormal modes can be reliably obtained [59].Then, the Konoplya formula reads: Konoplya's formula for the quasinormal modes, as mentioned earlier, incorporates several elements.The term V ′′ 0 represents the second derivative of the potential evaluated at its maximum point r 0 , and Λ j are the constants that depend on the effective potential as well as of its derivatives at the maximum.It is worth noting that recent advancements in the field have introduced a 13th-order WKB approximation, proposed by Matyjasek and Opala [163], which significantly improve the accuracy for the calculation of the quasinormal frequencies.
Presented in Tables VII, VIII, and IX, we find a comprehensive compilation of quasinormal frequencies obtained using the third-order WKB method.These tables are organized based on the multipole number l and mass M .Notably, ω 0 , ω 1 , and ω 2 for M = 0.1, turn out to be unstable.Such behaviors may be attributed to the influence of dark matter, which is introduced in the initial conditions to yield a Simpson-Visser-like black hole solution.
It is crucial to emphasize that the quasinormal modes linked to the scalar field exhibit a negative imaginary part.This significant characteristic implies that these modes experience exponential decay over time, signifying the dissipation of energy through scalar waves.This finding aligns with earlier investigations examining scalar, electromagnetic, and gravitational perturbations in spherically symmetric geometries [55,57,61,164].
In a broader context, it becomes evident that as M varies, there is an increase in the real part, coupled with a decrease in the imaginary part of the quasinormal modes.This observation highlights the crucial role played by the mass parameter in governing the damping behavior of the scalar waves.Additionally, for all values of l under consideration, the quasinormal modes demonstrate a trend of exhibiting damper frequencies as ω n increases.
It is important to mention that similar studies have been recently proposed in the literature, considering Hayward-like regular black holes [82,165].Furthermore, Ref. [166] analyzed a regular black-bounce model similar to the Simpson-Visser solution, deriving the thermodynamics, phase transition, Hawking radiation, light ring, and the quasinormal modes.

IX. CONCLUSION
In this study, we have focused on exploring the properties of a regular black hole within the framework of Verlinde's emergent gravity, with particular emphasis on the modified Simpson-Visser solution.Our analysis has revealed the existence of a single physical event horizon under certain conditions.
we have incorporated an extra application within the context of high-energy astrophysical phenomena, specifically focusing on the deposition of neutrino energy.
Furthermore, we have employed third-order WKB approximations to study the behavior of quasinormal modes, offering valuable information about the damping and decay properties of scalar waves in the context of the black hole surrounded by dark matter.Within the context of the dark matter scenario, another aspect worthy of investigation is the influence of charge Q on various aspects, such as the quasinormal modes, thermodynamic properties, shadows, and the photon sphere.
In the realm of future research, a fascinating direction for further investigation pertains to the analysis of quantum tunneling radiation and its associated backreaction effects, similar to the study encountered in Ref. [125].These and other ideas are now under development.

Figure 1 :
Figure 1: On the left side, we depict f (r) for various values of a (with M = 0.1), while on the right side, we illustrate f (r) for different values of M (with a = 0.1).

Figure 2 :
Figure 2: The modified Simpson-Visser horizon, r + , is displayed.A discontinuity occurs when M = 0.25 and the mass must obey the following constraint: M > 0.14699 (when a = 1).It is also shown the comparison between the Simpson-Visser horizon and the Schwarzschild one (the bottom plot).

Figure 3 :
Figure 3: The modified Hawking temperature is displayed, emphasizing the occurrence of a phase transition at M ≈ 0.1568.

Fig. 4 .Figure 4 :
Fig. 4. The entropy reveals at least two more points where the phase transition occurs, i.e., M = 0.25 and M = 0.5.Now, the last remaining thermodynamic quantity under consideration is the heat capacity

Figure 5
Figure 5 depicts the heat capacity for different values of mass M .Consistent with our analysis of the Hawking temperature, we observe a phase transition occurring around M ≈ 0.1568.This feature also occurs for other configurations of the system, specifically M = 0.25, M = 0.5, and M = 1.2.The first two values of mass are supported by the entropy behavior.

Figure 5 :
Figure 5: The modified heat capacity is displayed, highlighting multiple phase transition points.

5 Figure 6 :
Figure 6: The depicted trajectories illustrate the paths of light, with the photon sphere indicated by red dashed lines.Additionally, the event horizon r + is denoted by a black dot.

Figure 7 :
Figure 7: The circles in the illustration depict the shadows cast by varying values of massM and a, corresponding to the photon sphere r c − .Specifically, on the left side, the outer radius corresponds to M = 0.10, followed by M = 0.11, M = 0.12, and the inner radius at M = 0.13, all for a = 0.1.Furthermore, on the right side, we explore values for a at M = 0.10, ranging from the innermost at a = 0.01 to the outermost at a = 0.13, i.e., the values regarded to this latter case are a = 0.01, a = 0.11, a = 0.12, and a = 0.13.

Figure 8 :
Figure 8: The circles symbolize shadows cast by varying masses M corresponding to the photon sphere r c + .To elaborate, on the left side, the sequence progresses from the inner radius at M = 0.10, followed by M = 0.11, culminating with the outer radius at M = 0.12, all with a fixed parameter a = 0.1.Conversely, on the right side, we systematically explore different values of a for a constant mass M = 0.10, ranging from the outermost at a = 0.10

Figure 9 :
Figure 9: The effective potential V ef f is depicted as a function of the tortoise coordinate r * , specifically considering a particular value of

l
(ω) hold paramount significance.These ones are fundamental to explore the quasinormal modes of a black hole, represented by frequencies ω nl that satisfy the condition A (−) l (ω nl ) = 0.Such modes exhibit a unique behavior, with purely outgoing waves at spatial infinity and purely ingoing waves at the event horizon.The integers n and l represent the overtone and multipole numbers.
The author would like to thank Fundação de Apoio à Pesquisa do Estado da Paraíba (FAPESQ) and Conselho Nacional de Desenvolvimento Cientíıfico e Tecnológico (CNPq) -[200486/2022-5] and [150891/2023-7] for the financial support.Most of the calculations were performed by using the Mathematica software.Furthermore, the author extends gratitude to N. Heidari for invaluable suggestions and for providing the code used to compute the revised version of this manuscript.

Table I :
The thermodynamic properties of the Schwarzschild and Simpson-Visser like black holes.

Table VI :
Variation of r min and the deflection angle for various a and impact parameter (L/E) = 10 (when M = 0.1).