Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
\LetLtxMacro\oldsqrt

Rotational influence on fermions within negative curvature wormholes

Abdullah Guvendi 1111 E-mail: abdullah.guvendi@erzurum.edu.tr   Semra Gurtas Dogan 2222 E-mail: semragurtasdogan@hakkari.edu.tr (Corr. Auth.)   R. L. L. Vitória 3333 E-mail: ricardo-luis91@hotmail.com
1 Department of Basic Sciences, Erzurum Technical University, 25050, Erzurum, Türkiye
2 Department of Medical Imaging Techniques, Hakkari University, 30000, Hakkari, Türkiye
3 Faculdade de Física, Universidade Federal do Pará, Av. Augusto Corrêa, Guamá, Belém, PA 66075-110, Brazil
Abstract

In this research, we examine relativistic fermions within the rotating frame of negative curvature wormholes. Initially, as is typical in our context, we introduce the wormholes by embedding a curved surface into a higher dimensional flat Minkowski spacetime. Subsequently, we derive the spacetime metric that characterizes the rotating frame of these wormholes. We then investigate analytical solutions of the generalized Dirac equation within this framework. Through exploring a second-order non-perturbative wave equation, we seek exact solutions for fermions within the rotating frame of hyperbolic and elliptic wormholes, also known as negative curvature wormholes. Our analysis provides closed-form energy expressions, and we generalize our findings to Weyl fermions. By considering the impact of the rotating frame and curvature radius of wormholes, we discuss how these factors affect the evolution of fermionic fields, offering valuable insights into their behavior.

Keywords: Fermions; Non-inertial effects; Wormhole; Rotating reference frame; Graphene; Weyl fermions


1 Introduction

Landau and Lifshitz introduced rotation to Minkowski spacetime, analyzing two specific effects: the impact on clocks within a rotating body and the singular behavior at significant distances within a uniformly rotating frame [1]. Notably, the latter effect imposes a limitation on the radial coordinate. Within the domain of relativistic quantum theory, investigating the impact of a rotating frame involves understanding the principles of special relativity, where the quantized nature of fields and particle behaviors are crucial considerations. This exploration unveils intriguing and unconventional effects absent in classical Sagnac-type descriptions, representing the coupling between angular momentum and rotation angular frequency of a uniformly rotating frame, emphasized across various references [2, 3, 4, 5]. In the non-relativistic context, this effect is known as the Page-Werner et al. term [6, 7, 8]. The incorporation of relativistic corrections within this framework underscores the potential for phase shifts induced by relative motion between an observer and the rotating frame. Analyzing relativistic quantum systems within such a rotating reference frame unveils a realm of interesting results. Additionally, for quantum fields within curved spaces, the modified effective gravitational field, influenced by the centrifugal force, alters the energy levels of quantum entities, especially in confined or bound states [9]. In this context, several studies explore the impact of nontrivial rotating spacetimes on quantum mechanical systems. These investigations encompass various scenarios such as scalar bosons within cosmic string spacetime [10], a scalar field within a Kaluza-Klein theory [11], and a scalar field within time-dislocation spacetime [12]. Predicting the effects of a rotating frame on relativistic quantum particles is theoretically plausible; however, practically observing these minute alterations demands precision, control, and advanced techniques in quantum optics, atomic physics, or condensed matter physics, alongside cutting-edge technology [13, 14, 15]. Directly observing these effects necessitates trapping and controlling particles within the rotating frame, potentially detectable through variations in energy profiles, trajectories, or alterations in total angular momentum or spin. While extensive theoretical studies have scrutinized the effects of rotating reference frames on physical systems within non-relativistic and relativistic frameworks [16, 17, 2, 3, 4, 5, 9, 18, 19, 20, 21, 22, 7, 23, 24, 25, 26, 27], there is a noticeable gap in announced results regarding the evolution of relativistic fermions within the rotating frame of graphene or a graphene wormhole characterized by a two-dimensional curved surface of constant negative curvature. This unexplored domain requires further investigation to comprehend the intricate interplay between relativistic fermionic fields and specific curved geometries such as graphene structures.

In astrophysics, wormholes emerge from solutions in general relativity, serving as connections between two distinct regions [28, 29, 30]. These phenomena may manifest in both astrophysical contexts and condensed matter systems (by embedding curved surfaces into higher-dimensional spacetime (see also [31, 32])). However, the latter depiction differs because it involves two-dimensional curved surfaces, with time remaining unaffected by the background curvature [33]. Wormholes are actively researched to understand the fundamental properties of spacetime and whether they can exist within our current understanding of physics. Current research on wormholes spans from astrophysics to condensed matter contexts [34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47]. Research examining quantum systems evolving within curved spaces has garnered considerable interest [33, 49, 48, 50, 51, 52, 53, 54, 55, 56, 57], exemplified by investigations into Weyl particles in optical backgrounds [50], fermion-antifermion pairs in magnetized elliptic wormholes [52], and more [58, 59]. No results were found on fermionic fields in the rotating frame of negative curvature wormholes. This research addresses a gap in the literature by offering non-perturbative insights into fermions within rotating graphene and wormhole-like graphene structures characterized by a rotating surface with negative Gaussian curvature. Pursuing exact solutions holds promise in enabling the control or observation of compelling physical properties, potentially fostering novel phenomena, validating fundamental theories, and facilitating technological advancements. Our focus is on investigating relativistic fermion evolution within the rotating framework of negative curvature wormholes by seeking exact solutions to the corresponding Dirac equation.

Our examination of Dirac particles within the rotating frame of negative curvature wormholes is structured as follows: In section 2, we unveil a general and non-perturbative wave equation, exploring the intricate realm of relativistic fermions in the rotating frame of negative curvature wormholes. Through a thorough analysis, we present the main wave equation governing their behavior in these environments. In section 3, we uncover precise and exact outcomes concerning the evolution of relativistic fermions within two distinct rotating frames: the hyperbolic wormhole and the elliptic wormhole. By investigating these specific scenarios, we reveal compelling insights into the behavior and dynamics of particles under these conditions. Finally, in section 4, we summarize our findings, encapsulating the essence of our results. Additionally, we provide a comprehensive discussion that explores the implications of our results across various physically plausible scenarios and limits. This reflective section aims to highlight the significance and broader ramifications of our work within the realm of theoretical physics.

2 Wave equation

Refer to caption
Figure 1: Geometric structure of the wormhole surface.

In this section, we formulate a radial equation system for relativistic fermions in the rotating frame of negative curvature wormholes. This type of wormhole is characterized geometrically by points on its surface, which are expressed parametrically as shown below (also referenced in [31, 32]):

r(u,v)=x(u,v)i^+y(u,v)j^+zk^,𝑟𝑢𝑣𝑥𝑢𝑣^𝑖𝑦𝑢𝑣^𝑗𝑧^𝑘\vec{r}\left(u,v\right)=x\left(u,v\right)\hat{i}+y\left(u,v\right)\hat{j}+z% \hat{k},over→ start_ARG italic_r end_ARG ( italic_u , italic_v ) = italic_x ( italic_u , italic_v ) over^ start_ARG italic_i end_ARG + italic_y ( italic_u , italic_v ) over^ start_ARG italic_j end_ARG + italic_z over^ start_ARG italic_k end_ARG ,

where x(u,v)=χ(u)cos(v)𝑥𝑢𝑣𝜒𝑢𝑣x\left(u,v\right)=\chi\left(u\right)\cos(v)italic_x ( italic_u , italic_v ) = italic_χ ( italic_u ) roman_cos ( italic_v ), y(u,v)=χ(u)sin(v)𝑦𝑢𝑣𝜒𝑢𝑣y\left(u,v\right)=\chi\left(u\right)\sin(v)italic_y ( italic_u , italic_v ) = italic_χ ( italic_u ) roman_sin ( italic_v ), and z(u)=\oldsqrt[]1χ,u2duz\left(u\right)=\int{\oldsqrt[\ ]{1-\chi_{,u}^{2}}}\,duitalic_z ( italic_u ) = ∫ [ ] 1 - italic_χ start_POSTSUBSCRIPT , italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_u, with ,u indicating the derivative with respect to xusuperscript𝑥𝑢x^{u}italic_x start_POSTSUPERSCRIPT italic_u end_POSTSUPERSCRIPT. The shape of the wormhole is generally described by the function χ(u)𝜒𝑢\chi\left(u\right)italic_χ ( italic_u ). The constraints on z𝑧zitalic_z arise from the relationship:

ds2=dx2+dy2+dz2=du2+χ2(u)dv2,𝑑superscript𝑠2𝑑superscript𝑥2𝑑superscript𝑦2𝑑superscript𝑧2𝑑superscript𝑢2superscript𝜒2𝑢𝑑superscript𝑣2ds^{2}=dx^{2}+dy^{2}+dz^{2}=du^{2}+\chi^{2}\left(u\right)dv^{2},italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_d italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_d italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_d italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_d italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_u ) italic_d italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ,

which results in the Hilbert horizon at χ,u=1\chi_{,u}=1italic_χ start_POSTSUBSCRIPT , italic_u end_POSTSUBSCRIPT = 1 [33, 48, 49]. Embedding a three-dimensional wormhole into a (3+1)-dimensional spacetime generates effective gravity or curvature within the altered spacetime, leading to the experience of spacetime curvature by any particle or quasi-particle present in this modified background spacetime. Now, let us define the (2+1)-dimensional wormhole spacetime coordinates as dxμ={cdT,du,dv}𝑑superscript𝑥𝜇𝑐𝑑𝑇𝑑𝑢𝑑𝑣dx^{\mu}=\{cdT,du,dv\}italic_d italic_x start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = { italic_c italic_d italic_T , italic_d italic_u , italic_d italic_v } and the (3+1)-dimensional Minkowski spacetime coordinates as dxμ={cdT,dx,dy,dz}𝑑superscript𝑥superscript𝜇𝑐𝑑𝑇𝑑𝑥𝑑𝑦𝑑𝑧dx^{\mu^{{}^{\prime}}}=\{cdT,dx,dy,dz\}italic_d italic_x start_POSTSUPERSCRIPT italic_μ start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT = { italic_c italic_d italic_T , italic_d italic_x , italic_d italic_y , italic_d italic_z }, where Greek indices denote the coordinates within the curved space. The transformation between these coordinates can be accomplished via the following matrix:

xμxν=(10000χ,ucos(v)χ,usin(v)\oldsqrt[]1χ,u20χ,usin(v)χ,ucos(v)0).\frac{\partial x^{\mu^{{}^{\prime}}}}{\partial x^{\nu}}=\left(\begin{array}[]{% cccc}1&0&0&0\\ 0&\chi_{,u}\cos(v)&\chi_{,u}\sin(v)&\oldsqrt[\ ]{1-\chi_{,u}^{2}}\\ 0&-\chi_{,u}\sin(v)&\chi_{,u}\cos(v)&0\end{array}\right).divide start_ARG ∂ italic_x start_POSTSUPERSCRIPT italic_μ start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_x start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT end_ARG = ( start_ARRAY start_ROW start_CELL 1 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL italic_χ start_POSTSUBSCRIPT , italic_u end_POSTSUBSCRIPT roman_cos ( italic_v ) end_CELL start_CELL italic_χ start_POSTSUBSCRIPT , italic_u end_POSTSUBSCRIPT roman_sin ( italic_v ) end_CELL start_CELL [ ] 1 - italic_χ start_POSTSUBSCRIPT , italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL - italic_χ start_POSTSUBSCRIPT , italic_u end_POSTSUBSCRIPT roman_sin ( italic_v ) end_CELL start_CELL italic_χ start_POSTSUBSCRIPT , italic_u end_POSTSUBSCRIPT roman_cos ( italic_v ) end_CELL start_CELL 0 end_CELL end_ROW end_ARRAY ) .

Consequently, the metric of the considered wormhole spacetime takes the following form (with a negative signature) [33, 50, 51, 52]:

ds2=c2dT2du2χ2(u)dv2.𝑑superscript𝑠2superscript𝑐2𝑑superscript𝑇2𝑑superscript𝑢2superscript𝜒2𝑢𝑑superscript𝑣2\displaystyle ds^{2}=c^{2}dT^{2}-du^{2}-\chi^{2}\left(u\right)dv^{2}.italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_T start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_d italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_u ) italic_d italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (2.1)

A wormhole of this nature is geometrically depicted as shown in Fig. 1. In Eq. (2.1), c𝑐citalic_c symbolizes the speed of light, and the function χ(u)𝜒𝑢\chi(u)italic_χ ( italic_u ) takes the form χ(u)=acosh(ur0)𝜒𝑢𝑎𝑢subscript𝑟0\chi\left(u\right)=a\,\cosh\left(\frac{u}{r_{0}}\right)italic_χ ( italic_u ) = italic_a roman_cosh ( divide start_ARG italic_u end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) for a hyperbolic wormhole and χ(u)=bsinh(ur0)𝜒𝑢𝑏𝑢subscript𝑟0\chi\left(u\right)=b\ \sinh\left(\frac{u}{r_{0}}\right)italic_χ ( italic_u ) = italic_b roman_sinh ( divide start_ARG italic_u end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) for an elliptic wormhole. Here, a𝑎aitalic_a (b𝑏bitalic_b) denotes the radius of the associated wormhole at the midpoint (u=0𝑢0u=0italic_u = 0) between two ends, while r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT represents the radius of curvature of the wormhole surface u^^𝑢\hat{u}over^ start_ARG italic_u end_ARG. It is noteworthy that the Gaussian curvature for the spacetime background as expressed in Eq. (2.1) is K=χ,uuχK=-\frac{\chi_{,uu}}{\chi}italic_K = - divide start_ARG italic_χ start_POSTSUBSCRIPT , italic_u italic_u end_POSTSUBSCRIPT end_ARG start_ARG italic_χ end_ARG. Furthermore, Eq. (2.1) describes a curved surface exhibiting negative Gaussian curvature, making such curved surfaces suitable for describing a monolayer graphene sheet and graphene wormholes. These surfaces can also be derived using the optical metric of the Banados-Teitelboim-Zanelli (BTZ) black hole [48], referred to as negative curvature wormholes (generally) or hyperbolic/elliptic wormholes depending on the shape of the radius function χ(u)𝜒𝑢\chi\left(u\right)italic_χ ( italic_u ) (for further details, see [33, 48, 49]). In this scenario, we examine the rotating frame of spacetime given by Eq. (2.1). This can be achieved through transformations provided by T=t𝑇𝑡T=titalic_T = italic_t, u=ρ𝑢𝜌u=\rhoitalic_u = italic_ρ, and v=ϕ+ωrfT𝑣italic-ϕsubscript𝜔𝑟𝑓𝑇v=\phi+\omega_{rf}Titalic_v = italic_ϕ + italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT italic_T, where ωrfsubscript𝜔𝑟𝑓\omega_{rf}italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT denotes the angular frequency of the rotating frame. These transformations allow us to derive the spacetime interval describing a uniformly rotating frame of negative curvature wormholes (also see [60]).

ds2=(c2ωrf2χ2)dt22ωrfχ2dtdϕdρ2χ2dϕ2,𝑑superscript𝑠2superscript𝑐2superscriptsubscript𝜔𝑟𝑓2superscript𝜒2𝑑superscript𝑡22subscript𝜔𝑟𝑓superscript𝜒2𝑑𝑡𝑑italic-ϕ𝑑superscript𝜌2superscript𝜒2𝑑superscriptitalic-ϕ2\displaystyle ds^{2}=(c^{2}-\omega_{rf}^{2}\chi^{2})dt^{2}-2\omega_{rf}\chi^{2% }dtd\phi-d\rho^{2}-\chi^{2}d\phi^{2},italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_t italic_d italic_ϕ - italic_d italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_ϕ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (2.2)

for which the Gaussian curvature remains unchanged but for uρ𝑢𝜌u\rightarrow\rhoitalic_u → italic_ρ. In this context, it is crucial to observe that the coordinate ρ𝜌\rhoitalic_ρ needs to adhere to a condition derived from the relationship χ(ρ)<c/ωrf𝜒𝜌𝑐subscript𝜔𝑟𝑓\chi\left(\rho\right)<c/\omega_{rf}italic_χ ( italic_ρ ) < italic_c / italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT for the fermions to reside within the light-cone. This introduces an additional requirement, indicating that the corresponding wave function must vanish at the point determined by the condition χ(ρ)=c/ωrf𝜒𝜌𝑐subscript𝜔𝑟𝑓\chi\left(\rho\right)=c/\omega_{rf}italic_χ ( italic_ρ ) = italic_c / italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT. Referring to the line element Eq. (2.2), one can derive the covariant (gμν)subscript𝑔𝜇𝜈(g_{\mu\nu})( italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ) and contravariant (gμν)superscript𝑔𝜇𝜈(g^{\mu\nu})( italic_g start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT ) metric tensors

gμν=(c2ωrf2χ20ωrfχ2010ωrfχ20χ2),subscript𝑔𝜇𝜈superscript𝑐2superscriptsubscript𝜔𝑟𝑓2superscript𝜒20subscript𝜔𝑟𝑓superscript𝜒2010subscript𝜔𝑟𝑓superscript𝜒20superscript𝜒2\displaystyle g_{\mu\nu}=\left(\begin{array}[]{ccc}c^{2}-\omega_{rf}^{2}\chi^{% 2}&0&-\omega_{rf}\chi^{2}\\ 0&-1&0\\ -\omega_{rf}\chi^{2}&0&-\chi^{2}\end{array}\right),italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = ( start_ARRAY start_ROW start_CELL italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL start_CELL 0 end_CELL start_CELL - italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL - 1 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL - italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL start_CELL 0 end_CELL start_CELL - italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL end_ROW end_ARRAY ) ,
gμν=(1/c20ωrf/c2010ωrf/c20ωrf2/c2χ2).superscript𝑔𝜇𝜈1superscript𝑐20subscript𝜔𝑟𝑓superscript𝑐2010subscript𝜔𝑟𝑓superscript𝑐20superscriptsubscript𝜔𝑟𝑓2superscript𝑐2superscript𝜒2\displaystyle g^{\mu\nu}=\left(\begin{array}[]{ccc}1/c^{2}&0&-\omega_{rf}/c^{2% }\\ 0&-1&0\\ -\omega_{rf}/c^{2}&0&\omega_{rf}^{2}/c^{2}-\chi^{-2}\end{array}\right).italic_g start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT = ( start_ARRAY start_ROW start_CELL 1 / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL start_CELL 0 end_CELL start_CELL - italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL - 1 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL - italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL start_CELL 0 end_CELL start_CELL italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_χ start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT end_CELL end_ROW end_ARRAY ) .

Now, let us express the Dirac equation in its generalized form within this specific background [33, 50]

[γμμ+im~2]Ψ(xμ)=0,μ=t,ρ,ϕ,formulae-sequencedelimited-[]superscript𝛾𝜇subscriptscript-D̸𝜇𝑖~𝑚subscript2Ψsuperscript𝑥𝜇0𝜇𝑡𝜌italic-ϕ\displaystyle\left[\gamma^{\mu}\not{\mathcal{D}}_{\mu}+i\tilde{m}\mathcal{I}_{% 2}\right]\Psi(x^{\mu})=0,\quad\mu=t,\rho,\phi,[ italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT caligraphic_D̸ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT + italic_i over~ start_ARG italic_m end_ARG caligraphic_I start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ] roman_Ψ ( italic_x start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ) = 0 , italic_μ = italic_t , italic_ρ , italic_ϕ ,
μ=μΓμ,m~=mc.formulae-sequencesubscriptscript-D̸𝜇subscript𝜇subscriptΓ𝜇~𝑚𝑚𝑐Planck-constant-over-2-pi\displaystyle\not{\mathcal{D}}_{\mu}=\partial_{\mu}-\Gamma_{\mu},\quad\tilde{m% }=\frac{mc}{\hbar}.caligraphic_D̸ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT = ∂ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT - roman_Γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT , over~ start_ARG italic_m end_ARG = divide start_ARG italic_m italic_c end_ARG start_ARG roman_ℏ end_ARG . (2.3)

In this context, m𝑚mitalic_m denotes the fermion’s rest mass, Planck-constant-over-2-pi\hbarroman_ℏ represents the conventional Planck constant, c𝑐citalic_c stands for the speed of light, 2subscript2\mathcal{I}_{2}caligraphic_I start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT signifies the 2-dimensional identity matrix, Ψ(xμ)Ψsuperscript𝑥𝜇\Psi(x^{\mu})roman_Ψ ( italic_x start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ) symbolizes the Dirac field dependent on the spacetime position vector xμsuperscript𝑥𝜇x^{\mu}italic_x start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT, and γμsuperscript𝛾𝜇\gamma^{\mu}italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT denotes the space-dependent Dirac matrices, derived from the relation: γμ=ekμγksuperscript𝛾𝜇subscriptsuperscript𝑒𝜇𝑘superscript𝛾𝑘\gamma^{\mu}=e^{\mu}_{k}\gamma^{k}italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = italic_e start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT, where ekμsubscriptsuperscript𝑒𝜇𝑘e^{\mu}_{k}italic_e start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT are the inverse tetrad fields (k=0,1,2.𝑘012k=0,1,2.italic_k = 0 , 1 , 2 .) and γksuperscript𝛾𝑘\gamma^{k}italic_γ start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT denote the free Dirac matrices. The space-independent (free) Dirac matrices are represented in terms of the Pauli spin matrices (σx,σy,σzsuperscript𝜎𝑥superscript𝜎𝑦superscript𝜎𝑧\sigma^{x},\sigma^{y},\sigma^{z}italic_σ start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT , italic_σ start_POSTSUPERSCRIPT italic_y end_POSTSUPERSCRIPT , italic_σ start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT) as follows: γ0=σz,γ1=iσx,γ2=iσyformulae-sequencesuperscript𝛾0superscript𝜎𝑧formulae-sequencesuperscript𝛾1𝑖superscript𝜎𝑥superscript𝛾2𝑖superscript𝜎𝑦\gamma^{0}=\sigma^{z},\gamma^{1}=i\sigma^{x},\gamma^{2}=i\sigma^{y}italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT = italic_σ start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT , italic_γ start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT = italic_i italic_σ start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT , italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_i italic_σ start_POSTSUPERSCRIPT italic_y end_POSTSUPERSCRIPT, in accordance with the signature of the line element [61, 62, 63]. In Equation (2.3), ΓμsubscriptΓ𝜇\Gamma_{\mu}roman_Γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT refers to the spinorial affine connections for the Dirac field, attainable through the relation: Γλ=14[eν,λkekτΓνλτ]𝒮μν\Gamma_{\lambda}=\frac{1}{4}\left[e^{k}_{\nu_{,\lambda}}e^{\tau}_{k}-\Gamma_{% \nu\lambda}^{\tau}\right]\mathcal{S}^{\mu\nu}roman_Γ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 4 end_ARG [ italic_e start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT , italic_λ end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_τ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT - roman_Γ start_POSTSUBSCRIPT italic_ν italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ end_POSTSUPERSCRIPT ] caligraphic_S start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT, where indicates the derivative with respect to xλsuperscript𝑥𝜆x^{\lambda}italic_x start_POSTSUPERSCRIPT italic_λ end_POSTSUPERSCRIPT. Here, ΓνλτsuperscriptsubscriptΓ𝜈𝜆𝜏\Gamma_{\nu\lambda}^{\tau}roman_Γ start_POSTSUBSCRIPT italic_ν italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ end_POSTSUPERSCRIPT represents the Christoffel symbols Γνλτ=12gτϵ[νgλτ+λgϵνϵgνλ]superscriptsubscriptΓ𝜈𝜆𝜏12superscript𝑔𝜏italic-ϵdelimited-[]subscript𝜈subscript𝑔𝜆𝜏subscript𝜆subscript𝑔italic-ϵ𝜈subscriptitalic-ϵsubscript𝑔𝜈𝜆\Gamma_{\nu\lambda}^{\tau}=\frac{1}{2}g^{\tau\epsilon}\left[\partial_{\nu}g_{% \lambda\tau}+\partial_{\lambda}g_{\epsilon\nu}-\partial_{\epsilon}g_{\nu% \lambda}\right]roman_Γ start_POSTSUBSCRIPT italic_ν italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUPERSCRIPT italic_τ italic_ϵ end_POSTSUPERSCRIPT [ ∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_λ italic_τ end_POSTSUBSCRIPT + ∂ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_ϵ italic_ν end_POSTSUBSCRIPT - ∂ start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_ν italic_λ end_POSTSUBSCRIPT ], ekτsubscriptsuperscript𝑒𝜏𝑘e^{\tau}_{k}italic_e start_POSTSUPERSCRIPT italic_τ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT denotes the inverse tetrad fields, and 𝒮μνsuperscript𝒮𝜇𝜈\mathcal{S}^{\mu\nu}caligraphic_S start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT symbolizes the spin operators, defined as 𝒮μν=12[γμ,γν]superscript𝒮𝜇𝜈12superscript𝛾𝜇superscript𝛾𝜈\mathcal{S}^{\mu\nu}=\frac{1}{2}\left[\gamma^{\mu},\gamma^{\nu}\right]caligraphic_S start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT , italic_γ start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ] [61, 62, 63]. Now, the tetrads (and their inverses) can be obtained as follows (also refer to [33, 50]):

eμk=(c00010ωrfχ0χ),ekμ=(1/c00010ωrf/c0χ1)formulae-sequencesubscriptsuperscript𝑒𝑘𝜇𝑐00010subscript𝜔𝑟𝑓𝜒0𝜒subscriptsuperscript𝑒𝜇𝑘1𝑐00010subscript𝜔𝑟𝑓𝑐0superscript𝜒1\displaystyle e^{k}_{\mu}=\left(\begin{array}[]{ccc}c&0&0\\ 0&1&0\\ \omega_{rf}\chi&0&\chi\end{array}\right),e^{\mu}_{k}=\left(\begin{array}[]{ccc% }1/c&0&0\\ 0&1&0\\ -\omega_{rf}/c&0&\chi^{-1}\end{array}\right)italic_e start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT = ( start_ARRAY start_ROW start_CELL italic_c end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 1 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT italic_χ end_CELL start_CELL 0 end_CELL start_CELL italic_χ end_CELL end_ROW end_ARRAY ) , italic_e start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = ( start_ARRAY start_ROW start_CELL 1 / italic_c end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 1 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL - italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT / italic_c end_CELL start_CELL 0 end_CELL start_CELL italic_χ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT end_CELL end_ROW end_ARRAY )

since gμν=eμkeνlηklsubscript𝑔𝜇𝜈subscriptsuperscript𝑒𝑘𝜇subscriptsuperscript𝑒𝑙𝜈subscript𝜂𝑘𝑙g_{\mu\nu}=e^{k}_{\mu}e^{l}_{\nu}\eta_{kl}italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_e start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_l end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_η start_POSTSUBSCRIPT italic_k italic_l end_POSTSUBSCRIPT and ekμ=gμνeνlηklsubscriptsuperscript𝑒𝜇𝑘superscript𝑔𝜇𝜈subscriptsuperscript𝑒𝑙𝜈subscript𝜂𝑘𝑙e^{\mu}_{k}=g^{\mu\nu}e^{l}_{\nu}\eta_{kl}italic_e start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = italic_g start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_l end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_η start_POSTSUBSCRIPT italic_k italic_l end_POSTSUBSCRIPT where ηklsubscript𝜂𝑘𝑙\eta_{kl}italic_η start_POSTSUBSCRIPT italic_k italic_l end_POSTSUBSCRIPT is the Minkowski tensor, ηkl=diag(1,1,1)subscript𝜂𝑘𝑙diag111\eta_{kl}=\textrm{diag}\left(1,-1,-1\right)italic_η start_POSTSUBSCRIPT italic_k italic_l end_POSTSUBSCRIPT = diag ( 1 , - 1 , - 1 ). Accordingly, we have γt=σz/c,γρ=iσx,γϕ=ωrf/cσz+iχ1σyformulae-sequencesuperscript𝛾𝑡superscript𝜎𝑧𝑐formulae-sequencesuperscript𝛾𝜌𝑖superscript𝜎𝑥superscript𝛾italic-ϕsubscript𝜔𝑟𝑓𝑐superscript𝜎𝑧𝑖superscript𝜒1superscript𝜎𝑦\gamma^{t}=\sigma^{z}/c,\gamma^{\rho}=i\sigma^{x},\gamma^{\phi}=-\omega_{rf}/c% \,\sigma^{z}+i\,\chi^{-1}\sigma^{y}italic_γ start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT = italic_σ start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT / italic_c , italic_γ start_POSTSUPERSCRIPT italic_ρ end_POSTSUPERSCRIPT = italic_i italic_σ start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT , italic_γ start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT = - italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT / italic_c italic_σ start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + italic_i italic_χ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_σ start_POSTSUPERSCRIPT italic_y end_POSTSUPERSCRIPT [33, 50]. Also, non-vanishing components of the the spinorial affine connections are Γt=i2ωrfχ,ρσz\Gamma_{t}=\frac{i}{2}\omega_{rf}\chi_{,\rho}\sigma^{z}roman_Γ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = divide start_ARG italic_i end_ARG start_ARG 2 end_ARG italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT italic_χ start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT italic_σ start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT and Γϕ=i2χ,ρσz\Gamma_{\phi}=\frac{i}{2}\chi_{,\rho}\sigma^{z}roman_Γ start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT = divide start_ARG italic_i end_ARG start_ARG 2 end_ARG italic_χ start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT italic_σ start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT. Hence, one can verify that [33, 50]

γμΓμ=i2χ,ρχσx.\displaystyle\gamma^{\mu}\Gamma_{\mu}=-\frac{i}{2}\frac{\chi_{,\rho}}{\chi}% \sigma^{x}.italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT = - divide start_ARG italic_i end_ARG start_ARG 2 end_ARG divide start_ARG italic_χ start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT end_ARG start_ARG italic_χ end_ARG italic_σ start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT .

The metric Eq. (2.2) allows us to factorize the Dirac field as the following Ψ(xμ)=eiωteisϕ(ψ1(ρ)ψ2(ρ))TΨsuperscript𝑥𝜇superscripte𝑖𝜔𝑡superscripte𝑖𝑠italic-ϕsuperscriptsubscript𝜓1𝜌subscript𝜓2𝜌T\Psi\left(x^{\mu}\right)=\textrm{e}^{-i\omega t}\textrm{e}^{is\phi}\left(\psi_% {1}\left(\rho\right)\psi_{2}\left(\rho\right)\right)^{\textbf{T}}roman_Ψ ( italic_x start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ) = e start_POSTSUPERSCRIPT - italic_i italic_ω italic_t end_POSTSUPERSCRIPT e start_POSTSUPERSCRIPT italic_i italic_s italic_ϕ end_POSTSUPERSCRIPT ( italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_ρ ) italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_ρ ) ) start_POSTSUPERSCRIPT T end_POSTSUPERSCRIPT where ω𝜔\omegaitalic_ω is the relativistic frequency, s𝑠sitalic_s is the spin and T means transpose of the ρ𝜌\rhoitalic_ρ-dependent spinor. Finally, a system of coupled equations can be derived

[ϖm~]ψ1[ρ+12χ,ρχ+sχ]ψ2=0,\displaystyle\left[\varpi-\tilde{m}\right]\psi_{1}-\left[\partial_{\rho}+\frac% {1}{2}\frac{\chi_{,\rho}}{\chi}+\frac{s}{\chi}\right]\psi_{2}=0,[ italic_ϖ - over~ start_ARG italic_m end_ARG ] italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - [ ∂ start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_χ start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT end_ARG start_ARG italic_χ end_ARG + divide start_ARG italic_s end_ARG start_ARG italic_χ end_ARG ] italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0 ,
[ϖ+m~]ψ2+[ρ+12χ,ρχsχ]ψ1=0.\displaystyle\left[\varpi+\tilde{m}\right]\psi_{2}+\left[\partial_{\rho}+\frac% {1}{2}\frac{\chi_{,\rho}}{\chi}-\frac{s}{\chi}\right]\psi_{1}=0.[ italic_ϖ + over~ start_ARG italic_m end_ARG ] italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + [ ∂ start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_χ start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT end_ARG start_ARG italic_χ end_ARG - divide start_ARG italic_s end_ARG start_ARG italic_χ end_ARG ] italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0 . (2.4)

in which, ϖ=ω/c+sωrf/citalic-ϖ𝜔𝑐𝑠subscript𝜔𝑟𝑓𝑐\varpi=\omega/c+s\omega_{rf}/citalic_ϖ = italic_ω / italic_c + italic_s italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT / italic_c, describing fermions within the rotating frame of negatively curved wormholes (see also [48]). In the subsequent section, we aim to uncover analytical solutions to this wave equation within two specific contexts: (i) hyperbolic wormholes and (ii) elliptic wormholes.

3 Exact results

By decoupling the equations presented in Eq. (2.4), a non-perturbative second-order wave equation for the component ψ2subscript𝜓2\psi_{2}italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is obtained in the subsequent format

[ρ2+χ,ρχρ+Q^+(ϖ2m~2)]ψ2(ρ)=0,\displaystyle\left[\partial_{\rho}^{2}+\frac{\chi_{,\rho}}{\chi}\partial_{\rho% }+\hat{Q}+\left(\varpi^{2}-\tilde{m}^{2}\right)\right]\psi_{2}(\rho)=0,[ ∂ start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_χ start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT end_ARG start_ARG italic_χ end_ARG ∂ start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT + over^ start_ARG italic_Q end_ARG + ( italic_ϖ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - over~ start_ARG italic_m end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ] italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_ρ ) = 0 , (3.1)

where

Q^=12χ,ρρχs2χ2sχ,ρχ214χ,ρ2χ2.\displaystyle\hat{Q}=\frac{1}{2}\frac{\chi_{,\rho\rho}}{\chi}-\frac{s^{2}}{% \chi^{2}}-s\frac{\chi_{,\rho}}{\chi^{2}}-\frac{1}{4}\frac{\chi_{,\rho}^{2}}{% \chi^{2}}.over^ start_ARG italic_Q end_ARG = divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_χ start_POSTSUBSCRIPT , italic_ρ italic_ρ end_POSTSUBSCRIPT end_ARG start_ARG italic_χ end_ARG - divide start_ARG italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - italic_s divide start_ARG italic_χ start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT end_ARG start_ARG italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 1 end_ARG start_ARG 4 end_ARG divide start_ARG italic_χ start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (3.2)

Now, we have the ability to explore the precise solution of this wave equation under two distinct scenarios.

3.1 Relativistic fermions in the rotating frame of hyperbolic wormhole

If we examine a fermionic field within the relativistic framework of a rotating hyperbolic wormhole, the equation (3.1) transforms into the following form

ψ2,ρρ+tanh(ρr0)r0ψ2,ρ+[λ~+Q^1]ψ2=0,\displaystyle\psi_{2_{,\rho\rho}}+\frac{\tanh\left(\frac{\rho}{r_{0}}\right)}{% r_{0}}\,\psi_{2_{,\rho}}+\left[\tilde{\lambda}+\hat{Q}_{1}\right]\,\psi_{2}=0,italic_ψ start_POSTSUBSCRIPT 2 start_POSTSUBSCRIPT , italic_ρ italic_ρ end_POSTSUBSCRIPT end_POSTSUBSCRIPT + divide start_ARG roman_tanh ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG italic_ψ start_POSTSUBSCRIPT 2 start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT end_POSTSUBSCRIPT + [ over~ start_ARG italic_λ end_ARG + over^ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ] italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0 ,
Q^1=12r02s2/a2cosh2(ρr0)s/ar0tanh(ρr0)cosh(ρr0)tanh2(ρr0)4r02,subscript^𝑄112superscriptsubscript𝑟02superscript𝑠2superscript𝑎2superscript2𝜌subscript𝑟0𝑠𝑎subscript𝑟0𝜌subscript𝑟0𝜌subscript𝑟0superscript2𝜌subscript𝑟04subscriptsuperscript𝑟20\displaystyle\hat{Q}_{1}=\frac{1}{2r_{0}^{2}}-\frac{s^{2}/a^{2}}{\cosh^{2}% \left(\frac{\rho}{r_{0}}\right)}-\frac{s/a}{r_{0}}\,\frac{\tanh\left(\frac{% \rho}{r_{0}}\right)}{\cosh\left(\frac{\rho}{r_{0}}\right)}-\frac{\tanh^{2}% \left(\frac{\rho}{r_{0}}\right)}{4\,r^{2}_{0}},over^ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_cosh start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG - divide start_ARG italic_s / italic_a end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG divide start_ARG roman_tanh ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG start_ARG roman_cosh ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG - divide start_ARG roman_tanh start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG start_ARG 4 italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ,

where λ~=ϖ2m~2~𝜆superscriptitalic-ϖ2superscript~𝑚2\tilde{\lambda}=\varpi^{2}-\tilde{m}^{2}over~ start_ARG italic_λ end_ARG = italic_ϖ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - over~ start_ARG italic_m end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. To transform this unfamiliar equation into a more recognizable form, we can employ an assumed function, ψ2(ρ)=(sinh(ρ/r0)+i)p(sinh(ρ/r0)i)p*ψ(ρ)subscript𝜓2𝜌superscript𝜌subscript𝑟0𝑖𝑝superscript𝜌subscript𝑟0𝑖superscript𝑝*𝜓𝜌\psi_{2}\left(\rho\right)=(\sinh(\rho/r_{0})+i)^{p}(\sinh(\rho/r_{0})-i)^{p^{% \textbf{*}}}\psi(\rho)italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_ρ ) = ( roman_sinh ( italic_ρ / italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) + italic_i ) start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT ( roman_sinh ( italic_ρ / italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) - italic_i ) start_POSTSUPERSCRIPT italic_p start_POSTSUPERSCRIPT * end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT italic_ψ ( italic_ρ ), where p=32isr02a𝑝32𝑖𝑠subscript𝑟02𝑎p=\frac{3}{2}-i\frac{sr_{0}}{2a}italic_p = divide start_ARG 3 end_ARG start_ARG 2 end_ARG - italic_i divide start_ARG italic_s italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_a end_ARG and **{}^{\textbf{*}}start_FLOATSUPERSCRIPT * end_FLOATSUPERSCRIPT denotes the complex conjugate. Then, by introducing a new variable, ξ=12[1isinh(ρ/r0)]𝜉12delimited-[]1𝑖𝜌subscript𝑟0\xi=\frac{1}{2}[1-i\ \sinh(\rho/r_{0})]italic_ξ = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ 1 - italic_i roman_sinh ( italic_ρ / italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ], we arrive at the following wave equation

ξ(ξ1)ψ,ξξ[γ(α+β+1)ξ]ψ,ξ+αβψ=0,\displaystyle\xi\left(\xi-1\right)\psi_{,\xi\xi}-[\gamma-(\alpha+\beta+1)\xi]% \psi_{,\xi}+\alpha\beta\psi=0,italic_ξ ( italic_ξ - 1 ) italic_ψ start_POSTSUBSCRIPT , italic_ξ italic_ξ end_POSTSUBSCRIPT - [ italic_γ - ( italic_α + italic_β + 1 ) italic_ξ ] italic_ψ start_POSTSUBSCRIPT , italic_ξ end_POSTSUBSCRIPT + italic_α italic_β italic_ψ = 0 ,
α=1ir0\oldsqrt[]λ~,β=α,γ=32+isr0a,formulae-sequence𝛼1𝑖subscript𝑟0\oldsqrt~𝜆formulae-sequence𝛽superscript𝛼𝛾32𝑖𝑠subscript𝑟0𝑎\displaystyle\alpha=1-ir_{0}\oldsqrt[\ ]{\tilde{\lambda}},\quad\beta=\alpha^{*% },\quad\gamma=\frac{3}{2}+i\frac{sr_{0}}{a},italic_α = 1 - italic_i italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [ ] over~ start_ARG italic_λ end_ARG , italic_β = italic_α start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_γ = divide start_ARG 3 end_ARG start_ARG 2 end_ARG + italic_i divide start_ARG italic_s italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_a end_ARG , (3.3)

in which α𝛼\alphaitalic_α, β𝛽\betaitalic_β and γ𝛾\gammaitalic_γ are constants. This equation is the renowned Gauss’ hypergeometric equation, and its solution in the vicinity of the origin can be expressed as ψ(ξ)=𝒩2F1(α,β;γ;ξ)𝜓𝜉subscript𝒩2subscript𝐹1𝛼𝛽𝛾𝜉\psi\left(\xi\right)=\mathcal{N}\ _{2}F_{1}\left(\alpha,\beta;\gamma;\xi\right)italic_ψ ( italic_ξ ) = caligraphic_N start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_α , italic_β ; italic_γ ; italic_ξ ) [64]. Here, 𝒩𝒩\mathcal{N}caligraphic_N denotes an arbitrary constant. The polynomial condition associated with this function is α=n𝛼𝑛\alpha=-nitalic_α = - italic_n [50, 64], where n𝑛nitalic_n signifies the principal quantum number (n=0,1,2..)(n=0,1,2..)( italic_n = 0 , 1 , 2 . . ). This criterion leads to the quantization principle governing the formation of the system under consideration, consequently yielding the subsequent expression

ns=sωrf±c\oldsqrt[]m2c22(n+1)2r02,subscript𝑛𝑠plus-or-minus𝑠Planck-constant-over-2-pisubscript𝜔𝑟𝑓Planck-constant-over-2-pi𝑐\oldsqrtsuperscript𝑚2superscript𝑐2superscriptPlanck-constant-over-2-pi2superscript𝑛12superscriptsubscript𝑟02\displaystyle\mathcal{E}_{ns}=-s\hbar\omega_{rf}\pm\hbar c\oldsqrt[\ ]{\frac{m% ^{2}c^{2}}{\hbar^{2}}-\frac{(n+1)^{2}}{r_{0}^{2}}},caligraphic_E start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT = - italic_s roman_ℏ italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT ± roman_ℏ italic_c [ ] divide start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG ( italic_n + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (3.4)

for the relativistic energy of the system under scrutiny. In this context, it becomes evident that the energy of such a system relies on the properties of particles, including spin (assuming ωrf0subscript𝜔𝑟𝑓0\omega_{rf}\neq 0italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT ≠ 0), particle’s rest mass, and the curvature radius of the hyperbolic wormhole, in addition to fundamental constants like Planck-constant-over-2-pi\hbarroman_ℏ and c𝑐citalic_c. The background under consideration is a spacetime with constant negative Gaussian curvature, characteristic of structures like monolayer graphene sheets or graphene wormholes since K=1r02𝐾1superscriptsubscript𝑟02K=-\frac{1}{r_{0}^{2}}italic_K = - divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (refer to [48]). Consequently, there is potential utility in adjusting the outcomes for Weyl fermions by setting m=0𝑚0m=0italic_m = 0 and substituting c𝑐citalic_c with νfsubscript𝜈𝑓\nu_{f}italic_ν start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT, where νfsubscript𝜈𝑓\nu_{f}italic_ν start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT denotes the Fermi velocity, approximately equal to c/300𝑐300c/300italic_c / 300. Notably, equation (3.4) transforms accordingly

ns=sωrf±iνf(n+1)r0,subscript𝑛𝑠plus-or-minus𝑠Planck-constant-over-2-pisubscript𝜔𝑟𝑓𝑖Planck-constant-over-2-pisubscript𝜈𝑓𝑛1subscript𝑟0\displaystyle\mathcal{E}_{ns}=-s\hbar\omega_{rf}\pm i\hbar\nu_{f}\frac{(n+1)}{% r_{0}},caligraphic_E start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT = - italic_s roman_ℏ italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT ± italic_i roman_ℏ italic_ν start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT divide start_ARG ( italic_n + 1 ) end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG , (3.5)

if m=0𝑚0m=0italic_m = 0. The Eq. (3.5) indicates that the energy of this system possesses a complex nature. Consequently, these states cannot remain stable, and the associated modes either decay or grow over time since Ψeinstproportional-toΨsuperscripte𝑖subscript𝑛𝑠Planck-constant-over-2-pi𝑡\Psi\propto\textrm{e}^{-i\frac{\mathcal{E}_{ns}}{\hbar}t}roman_Ψ ∝ e start_POSTSUPERSCRIPT - italic_i divide start_ARG caligraphic_E start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT end_ARG start_ARG roman_ℏ end_ARG italic_t end_POSTSUPERSCRIPT. Additionally, the real oscillation of a Weyl particle is solely influenced by the coupling between the particle’s spin and the rotational frequency of the frame. Our findings suggest that the stability of the background spacetime is compromised under perturbations caused by massless particles, while it remains stable under perturbations from Weyl antiparticles (also see [51, 52]) when ωrf=0subscript𝜔𝑟𝑓0\omega_{rf}=0italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT = 0. Notably, the coupling sωrf𝑠subscript𝜔𝑟𝑓s\omega_{rf}italic_s italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT might induce symmetry breaking around the Dirac point when m0𝑚0m\neq 0italic_m ≠ 0. Furthermore, it appears conceivable to measure, in principle, only the energy stemming from the sωrf𝑠subscript𝜔𝑟𝑓s\omega_{rf}italic_s italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT coupling if the Dirac particle possesses a critical (or effective) mass value, mc=(n+1)cr0subscript𝑚𝑐Planck-constant-over-2-pi𝑛1𝑐subscript𝑟0m_{c}=\frac{\hbar\ (n+1)}{c\ r_{0}}italic_m start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG roman_ℏ ( italic_n + 1 ) end_ARG start_ARG italic_c italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG. Such a scenario could arise particularly in condensed matter systems, especially when r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is on the order of a few nanometers (also see [50, 51]).

3.2 Relativistic fermions in the rotating frame of elliptic wormhole

In this scenario, we derive precise outcomes for relativistic fermions within the rotating frame of an elliptical wormhole. Consequently, the wave equation presented in Eq. (3.1) manifests as follows:

ψ2,ρρ+coth(ρr0)r0ψ2,ρ+[λ+Q^2]ψ2=0,\displaystyle\psi_{2_{,\rho\rho}}+\frac{\coth\left(\frac{\rho}{r_{0}}\right)}{% r_{0}}\,\psi_{2_{,\rho}}+\left[\lambda+\hat{Q}_{2}\right]\,\psi_{2}=0,italic_ψ start_POSTSUBSCRIPT 2 start_POSTSUBSCRIPT , italic_ρ italic_ρ end_POSTSUBSCRIPT end_POSTSUBSCRIPT + divide start_ARG roman_coth ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG italic_ψ start_POSTSUBSCRIPT 2 start_POSTSUBSCRIPT , italic_ρ end_POSTSUBSCRIPT end_POSTSUBSCRIPT + [ italic_λ + over^ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ] italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0 ,
Q^2=12r02s2/b2sinh2(ρr0)s/br0coth(ρr0)sinh(ρr0)coth2(ρr0)4r02.subscript^𝑄212superscriptsubscript𝑟02superscript𝑠2superscript𝑏2superscript2𝜌subscript𝑟0𝑠𝑏subscript𝑟0hyperbolic-cotangent𝜌subscript𝑟0𝜌subscript𝑟0superscripthyperbolic-cotangent2𝜌subscript𝑟04subscriptsuperscript𝑟20\displaystyle\hat{Q}_{2}=\frac{1}{2r_{0}^{2}}-\frac{s^{2}/b^{2}}{\sinh^{2}% \left(\frac{\rho}{r_{0}}\right)}-\frac{s/b}{r_{0}}\,\frac{\coth\left(\frac{% \rho}{r_{0}}\right)}{\sinh\left(\frac{\rho}{r_{0}}\right)}-\frac{\coth^{2}% \left(\frac{\rho}{r_{0}}\right)}{4\,r^{2}_{0}}.over^ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_sinh start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG - divide start_ARG italic_s / italic_b end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG divide start_ARG roman_coth ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG start_ARG roman_sinh ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG - divide start_ARG roman_coth start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_ρ end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_ARG start_ARG 4 italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG .

Let us attempt to reformulate this equation using a new variable substitution, z=12[1+cosh(ρ/r0)]𝑧12delimited-[]1𝜌subscript𝑟0z=\frac{1}{2}[1+\cosh(\rho/r_{0})]italic_z = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ 1 + roman_cosh ( italic_ρ / italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ], aiming to eliminate the hyperbolic functions. By expressing Eq. (3.2) in terms of the variable z𝑧zitalic_z, we can render it free of hyperbolic functions, facilitating an exploration of its asymptotic behaviors. Consequently, employing an ansatz function, ψ2(z)=z2sr0b4b(z1)2sr0+b4bψ(z)subscript𝜓2𝑧superscript𝑧2𝑠subscript𝑟0𝑏4𝑏superscript𝑧12𝑠subscript𝑟0𝑏4𝑏𝜓𝑧\psi_{2}\left(z\right)=z^{\frac{2sr_{0}-b}{4b}}\left(z-1\right)^{-\frac{2sr_{0% }+b}{4b}}\psi\left(z\right)italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_z ) = italic_z start_POSTSUPERSCRIPT divide start_ARG 2 italic_s italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_b end_ARG start_ARG 4 italic_b end_ARG end_POSTSUPERSCRIPT ( italic_z - 1 ) start_POSTSUPERSCRIPT - divide start_ARG 2 italic_s italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_b end_ARG start_ARG 4 italic_b end_ARG end_POSTSUPERSCRIPT italic_ψ ( italic_z ), the resulting equation can be transformed into the following wave equation

z(z1)ψ,zz[γ~(α~+β~+1)z]ψ,z+α~β~ψ=0,,\displaystyle z\left(z-1\right)\psi_{,zz}-[\tilde{\gamma}-(\tilde{\alpha}+% \tilde{\beta}+1)z]\psi_{,z}+\tilde{\alpha}\tilde{\beta}\psi=0,\quad,italic_z ( italic_z - 1 ) italic_ψ start_POSTSUBSCRIPT , italic_z italic_z end_POSTSUBSCRIPT - [ over~ start_ARG italic_γ end_ARG - ( over~ start_ARG italic_α end_ARG + over~ start_ARG italic_β end_ARG + 1 ) italic_z ] italic_ψ start_POSTSUBSCRIPT , italic_z end_POSTSUBSCRIPT + over~ start_ARG italic_α end_ARG over~ start_ARG italic_β end_ARG italic_ψ = 0 , ,
α~=ir0\oldsqrt[]λ~,β~=α~,γ~=12+sr0b.formulae-sequence~𝛼𝑖subscript𝑟0\oldsqrt~𝜆formulae-sequence~𝛽superscript~𝛼~𝛾12𝑠subscript𝑟0𝑏\displaystyle\tilde{\alpha}=-ir_{0}\oldsqrt[\ ]{\tilde{\lambda}},\quad\tilde{% \beta}=\tilde{\alpha}^{*},\quad\tilde{\gamma}=\frac{1}{2}+\frac{sr_{0}}{b}.over~ start_ARG italic_α end_ARG = - italic_i italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [ ] over~ start_ARG italic_λ end_ARG , over~ start_ARG italic_β end_ARG = over~ start_ARG italic_α end_ARG start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , over~ start_ARG italic_γ end_ARG = divide start_ARG 1 end_ARG start_ARG 2 end_ARG + divide start_ARG italic_s italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_b end_ARG .

In this context, α~~𝛼\tilde{\alpha}over~ start_ARG italic_α end_ARG, β~~𝛽\tilde{\beta}over~ start_ARG italic_β end_ARG, and γ~~𝛾\tilde{\gamma}over~ start_ARG italic_γ end_ARG represent fundamental constants within Gauss hypergeometric wave equation, as extensively documented in references [50, 64]. The regular solution of this equation in the vicinity of the origin can be elegantly expressed as ψ(z)=𝒞2F1(α~,β~;γ~;z)𝜓𝑧subscript𝒞2subscript𝐹1~𝛼~𝛽~𝛾𝑧\psi\left(z\right)=\mathcal{C}\ _{2}F_{1}\left(\tilde{\alpha},\tilde{\beta};% \tilde{\gamma};z\right)italic_ψ ( italic_z ) = caligraphic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( over~ start_ARG italic_α end_ARG , over~ start_ARG italic_β end_ARG ; over~ start_ARG italic_γ end_ARG ; italic_z ), where 𝒞𝒞\mathcal{C}caligraphic_C stands for an arbitrary constant. Remarkably, when the parameter α~~𝛼\tilde{\alpha}over~ start_ARG italic_α end_ARG assumes the value n𝑛-n- italic_n, with n𝑛nitalic_n representing the radial quantum number, the solution function ψ(z)𝜓𝑧\psi(z)italic_ψ ( italic_z ) transforms into a polynomial of degree n𝑛nitalic_n concerning the variable z𝑧zitalic_z. This establishes a link between the polynomial nature of the solution function and the quantum characteristics of the system, offering valuable insights into the system’s energy behavior at different quantum states. This condition establishes an intriguing connection to the system’s energy spectrum, unveiling the following set of energy levels

ns=sωrf±c\oldsqrt[]m2c22n2r02.subscript𝑛𝑠plus-or-minus𝑠Planck-constant-over-2-pisubscript𝜔𝑟𝑓Planck-constant-over-2-pi𝑐\oldsqrtsuperscript𝑚2superscript𝑐2superscriptPlanck-constant-over-2-pi2superscript𝑛2superscriptsubscript𝑟02\displaystyle\mathcal{E}_{ns}=-s\hbar\omega_{rf}\pm\hbar c\oldsqrt[\ ]{\frac{m% ^{2}c^{2}}{\hbar^{2}}-\frac{n^{2}}{r_{0}^{2}}}.caligraphic_E start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT = - italic_s roman_ℏ italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT ± roman_ℏ italic_c [ ] divide start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (3.6)

While bearing resemblance to the expression presented in Eq. (3.4), the relativistic energy pertaining to a fermion precisely equals the rest mass energy ±mc2plus-or-minus𝑚superscript𝑐2\pm mc^{2}± italic_m italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT under specific conditions where ωrf=0subscript𝜔𝑟𝑓0\omega_{rf}=0italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT = 0 and n=0𝑛0n=0italic_n = 0. This distinction stands as the fundamental dissimilarity between the outcomes delineated in Eqs. (3.4) and (3.6). Our investigations substantiate that within the rotating frame of an elliptic wormhole, the energy of a relativistic fermion manifests an intriguing independence from the curvature radius of the wormhole as the system converges toward its ground state. Notably, our observations assert that analogous deliberations, explicated in sec. 3.1, remain viable and generalizable when transitioning solely through n+1n𝑛1𝑛n+1\longrightarrow nitalic_n + 1 ⟶ italic_n. The inference arises from the energy profiles derived in both situations, suggesting that the baseline energy level for relativistic fermions in the rotating frame of a hyperbolic wormhole appears akin to the first excited state for those in the rotating frame of an elliptic wormhole. Moreover, this correlation remains intact even when ωrf=0subscript𝜔𝑟𝑓0\omega_{rf}=0italic_ω start_POSTSUBSCRIPT italic_r italic_f end_POSTSUBSCRIPT = 0.

4 Summary and discussions

In this study, we analyse the influence of non-inertial effects arising from the uniformly rotating reference frame of negative curvature wormholes, such as the hyperbolic wormhole and elliptic wormhole, on relativistic fermions. Our approach involves the investigation of analytically permissible solutions of the associated Dirac equation. Initially, we derive a comprehensive non-perturbative second-order wave equation. Subsequently, we proceed to ascertain the exact solutions of this wave equation within two distinct scenarios. In the first scenario, we focus on a background that characterizes the hyperbolic wormhole. Within this context, we determine the solution function expressed in terms of Gauss hypergeometric function. Consequently, we derive a non-perturbative energy spectrum in a closed mathematical form. The expression for this energy spectrum is provided by equation (3.4).

The findings reveal an intriguing interdependence between the energy states of relativistic fermions, the curvature radius of a wormhole, and the angular frequency of the uniformly rotating frame. Alongside the intrinsic quantum attributes of the fermionic field, our investigation underscores the direct influence of the particle’s spin on the energy manifestation. Notably, this interaction may induce symmetry breaking centered around the Dirac point (zero energy), with energy levels diminishing until they reach a critical curvature radius of the wormhole, precisely at r0=(n+1)mcsubscript𝑟0Planck-constant-over-2-pi𝑛1𝑚𝑐r_{0}=\frac{\hbar(n+1)}{mc}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = divide start_ARG roman_ℏ ( italic_n + 1 ) end_ARG start_ARG italic_m italic_c end_ARG, where n𝑛nitalic_n signifies the principle quantum number. Importantly, the system’s stability is contingent upon the relationship between r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and (n+1)mcPlanck-constant-over-2-pi𝑛1𝑚𝑐\frac{\hbar(n+1)}{mc}divide start_ARG roman_ℏ ( italic_n + 1 ) end_ARG start_ARG italic_m italic_c end_ARG; instability ensues if r0subscript𝑟0r_{0}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT surpasses (n+1)mcPlanck-constant-over-2-pi𝑛1𝑚𝑐\frac{\hbar(n+1)}{mc}divide start_ARG roman_ℏ ( italic_n + 1 ) end_ARG start_ARG italic_m italic_c end_ARG, irrespective of whether the angular frequency of the rotating frame is null or not. Our investigation further reveals intriguing insights into Weyl fermions (m=0𝑚0m=0italic_m = 0), suggesting that states associated with these fermions may not exhibit stability due to potential damped modes characterized by a decay time τn=/|nsIm|subscript𝜏𝑛Planck-constant-over-2-pisubscript𝑛subscript𝑠𝐼𝑚\tau_{n}=\hbar/|\mathcal{E}_{ns_{Im}}|italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = roman_ℏ / | caligraphic_E start_POSTSUBSCRIPT italic_n italic_s start_POSTSUBSCRIPT italic_I italic_m end_POSTSUBSCRIPT end_POSTSUBSCRIPT |, where nsImsubscript𝑛subscript𝑠𝐼𝑚\mathcal{E}_{ns_{Im}}caligraphic_E start_POSTSUBSCRIPT italic_n italic_s start_POSTSUBSCRIPT italic_I italic_m end_POSTSUBSCRIPT end_POSTSUBSCRIPT denotes the imaginary component of the energy. In an alternative scenario involving an elliptic wormhole background in spacetime, our analysis of the primary equation yields precise solutions for the relativistic energy within the system under consideration. Surprisingly, the ground state energy identified in the earlier scenario appears to correspond to the first excited state of relativistic fermions within the rotating frame of an elliptic wormhole. Intriguingly, when the angular velocity of the rotating frame reaches zero, our results demonstrate that the relativistic energy of the system aligns with the rest mass energy of the fermion upon reaching the ground state (n=0𝑛0n=0italic_n = 0). This alignment reflects an inevitable outcome given our understanding that the system will ultimately settle into its ground state. Furthermore, our outcomes suggest the possibility of measuring energy solely derived from the interplay between the particle’s spin and the angular frequency of the uniformly rotating frame. This measurement becomes feasible if the particle possesses a critical (or effective) mass value represented by mc=(n+1)cr0subscript𝑚𝑐Planck-constant-over-2-pi𝑛1𝑐subscript𝑟0m_{c}=\frac{\hbar(n+1)}{cr_{0}}italic_m start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG roman_ℏ ( italic_n + 1 ) end_ARG start_ARG italic_c italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG and mc=ncr0subscript𝑚𝑐Planck-constant-over-2-pi𝑛𝑐subscript𝑟0m_{c}=\frac{\hbar n}{cr_{0}}italic_m start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG roman_ℏ italic_n end_ARG start_ARG italic_c italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG for the first and second scenarios, respectively. Particularly intriguing is the potential manifestation of this specialized case in condensed matter systems, where c𝑐citalic_c can be equated to νfsubscript𝜈𝑓\nu_{f}italic_ν start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT (Fermi velocity). Notably, this scenario appears plausible when the curvature radius of the wormhole approximates a few nanometers, assuming the particle’s mass is the typical electron mass.

Recently, there has been a notable focus on performing thermodynamic evaluations using energy levels from quantum systems in complex environments. Several studies in the field have investigated this, setting the stage for future exploration of thermodynamic properties [65, 66, 67, 68, 69, 70, 71]. Our precise findings could be valuable for analyzing the thermal properties of graphene structures, presenting an intriguing avenue for future research.

Data availability

This manuscript does not contain any associated data.

Conflicts of interest statement

The authors have disclosed no conflicts of interest.

Funding

This research has not received any funding.

References