Massive U(1) gauge fields and their accompanying scalars in brane world

Figures(2) / Tables(5)

Get Citation
Ye-Hao Yang, Wen-Xuan Ma and Chun-E Fu. Massive U(1) gauge fields and their accompanying scalars in brane world[J]. Chinese Physics C. doi: 10.1088/1674-1137/ae4579
Ye-Hao Yang, Wen-Xuan Ma and Chun-E Fu. Massive U(1) gauge fields and their accompanying scalars in brane world[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ae4579 shu
Milestone
Received: 2025-12-24
Article Metric

Article Views(49)
PDF Downloads(2)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Massive U(1) gauge fields and their accompanying scalars in brane world

    Corresponding author: Chun-E Fu, fuche13@mail.xjtu.edu.cn
  • Institute of Modern Physics, Department of Modern Physics and Astronomy, School of Physics, Xi'an Jiaotong University, Xi'an ‌ 710049, China

Abstract: In brane-world scenarios, the effective action of a massless bulk $U(1)$ gauge field preserves gauge invariance via couplings between massive vector Kaluza-Klein (KK) and scalar KK modes. In this study, we extend this framework by introducing the term $(\nabla^M X_M)^2$ into the massless bulk $U(1)$ gauge action. This modification explicitly breaks full gauge redundancy while preserving residual gauge symmetry both in the bulk and on the brane. In this setup, scalar KK modes can acquire masses from the background geometry. We find that, on the five-dimensional brane, these scalar KK modes are lighter than the vector KK modes. On the six-dimensional brane, two types of scalar modes emerge, and mixed interactions between them give rise to oscillations among these scalar modes.

    HTML

    I.   INTRODUCTION
    • The idea that our Universe may contain additional spatial dimensions beyond the four-dimensional spacetime has drawn increasing attention over the past decades. The earliest attempt in this direction can be traced back to Kaluza–Klein (KK) theory, which proposed that introducing a fifth spatial dimension enables a unified description of gravity and electromagnetism [13]. This idea was later extended dramatically in string theory, where the presence of multiple extra dimensions is not optional but required for mathematical consistency. Building on these developments, modern brane-world models have provided further insights into the hierarchy problem and the cosmological constant problem. Representative examples include the ADD model [4], which assumes large but compact extra dimensions, and the Randall–Sundrum (RS) scenarios [5, 6], where warped extra dimensions give rise to novel gravitational dynamics.

      Within the brane-world framework [716], one particularly compelling line of research concerns the behavior of higher-dimensional matter fields and their associated Kaluza–Klein (KK) modes [1727]. Investigating these KK spectra is crucial not only for ensuring the theoretical consistency of higher-dimensional models but also for probing potential signatures of extra dimensions in observable physics. A growing body of literature has explored how KK excitations may lead to distinctive phenomenology, including modifications to gravitational interactions [2830], corrections to Standard Model processes [31, 32], and potential resonant modes localized near the brane [3335]. These works highlight the importance of KK modes as a window into the possible structure of higher-dimensional spacetime [3639].

      It is well established that massless fields propagating in the bulk give rise to massive KK modes in the effective four-dimensional theory on the brane, with the masses induced by the geometry of the extra dimensions. For a bulk U(1) gauge field, the KK tower on the brane includes both vector and scalar modes. Notably, the effective four-dimensional action preserves gauge invariance due to vector–scalar and scalar–scalar couplings among the KK modes [4042]. As a consequence, sufficient gauge freedom remains in the brane theory to eliminate the unphysical polarization degrees of freedom from the massive vector KK modes. This mechanism further reveals that the massive scalar KK modes play a role analogous to Stückelberg or Goldstone fields.

      Previous studies have indicated that certain mixed vector–scalar and scalar–scalar couplings may exist. However, since the basis functions for the scalar KK modes cannot be uniquely determined, the precise form of these couplings remains uncertain. What is important is that, due to such mixings, the scalar KK modes cannot be fully gauged away—implying that they may carry physical degrees of freedom [43]. To properly define these scalar modes, and motivated by the Stückelberg mechanism, we propose adding a covariant term $(\nabla^M X_M)^2$ to the bulk action. This term introduces a kinetic contribution for the scalar involving $\partial_z X_z$, thereby clarifying how the full geometry influences the scalar KK modes together with the existing term. At the same time, this addition preserves gauge symmetry both in the bulk and on the brane just like the Stückelberg mechanism.

      It is noteworthy that in a brane model with d extra dimensions ($d \gt 1$), there exist d distinct types of scalar KK modes. Interactions among the different types of scalar KK modes induce mixing in the mass matrix of the KK spectrum, leading to oscillation phenomena among these modes. Consequently, although the scalars acquire masses from the geometry of the extra dimensions, the physical states observed on the brane are mass eigenstates arising from this mixing. By diagonalizing the mass matrix, we compute the physical masses of these mixed modes in two illustrative 6D brane models.

      The paper is organized as follows: In Sec. 2, serving as a demonstrative example of the influence of the new term $(\nabla^M X_M)^2$, we investigate the gauge invariance and the mass spectrum of the KK modes in a 5D brane model. Then, in Sec. 3, we extend the study to a 6D brane, where two types of scalar KK modes emerge, leading to interesting oscillation phenomena. Numerical results related to these oscillations are presented in Sec. 4. Finally, a brief summary and discussion are provided in Sec. 5.

    • I.   INTRODUCTION
      • The idea that our Universe may contain additional spatial dimensions beyond the four-dimensional spacetime has gained increasing attention over the past decades. The earliest attempt in this direction can be traced back to the Kaluza–Klein (KK) theory, which proposes that introducing a fifth spatial dimension enables a unified description of gravity and electromagnetism [13]. Subsequently, this idea was extended dramatically in string theory, wherein the presence of multiple extra dimensions is not optional but required for mathematical consistency. Building on these developments, modern brane-world models provided further insights into the hierarchy and cosmological constant problems. Representative examples include the Arkani-Hamed–Dimopoulos–Dvali model [4], which assumes large but compact extra dimensions, and the Randall–Sundrum scenarios [5, 6], where warped extra dimensions induce novel gravitational dynamics.

        Within the brane-world framework [716], one compelling line of research concerns the behavior of higher-dimensional matter fields and their associated KK modes [1727]. Investigating these KK spectra is crucial not only for ensuring the theoretical consistency of higher-dimensional models but also for probing the potential signatures of extra dimensions in observable physics. A growing body of literature has explored how KK excitations can lead to distinctive phenomenology, including modifications to gravitational interactions [2830], corrections to standard model processes [31, 32], and potential resonant modes localized near the brane [3335]. These studies highlight the importance of KK modes as a window into the possible structure of higher-dimensional spacetime [3639].

        Massless fields propagating in the bulk give rise to massive KK modes in the effective four-dimensional (4D) theory on the brane, with masses induced by the geometry of extra dimensions. For a bulk U(1) gauge field, the KK tower on the brane includes both vector and scalar modes. The effective 4D action preserves gauge invariance caused by vector–scalar and scalar–scalar couplings among the KK modes [4042]. Consequently, sufficient gauge freedom remains in the brane theory to eliminate the unphysical polarization degrees of freedom from the massive vector KK modes. This mechanism reveals that massive scalar KK modes play a role analogous to Stückelberg or Goldstone fields.

        Previous studies indicated that certain mixed vector–scalar and scalar–scalar couplings may exist. However, the precise form of these couplings remains uncertain because the basis functions for the scalar KK modes cannot be determined uniquely. Owing to such mixings, scalar KK modes cannot be fully gauged, which implies that they may carry physical degrees of freedom [43]. To properly define these scalar modes, and motivated by the Stückelberg mechanism, we propose adding a covariant term $(\nabla^M X_M)^2$ to the bulk action. This term introduces a kinetic contribution for the scalar involving $\partial_z X_z$, which helps clarify how the full geometry affects the scalar KK modes with the existing term. Meanwhile, this addition preserves gauge symmetry in the bulk and on the brane, similar to the Stückelberg mechanism.

        In a brane model with d extra dimensions ($d \gt 1$), there exist d distinct types of scalar KK modes. Interactions among the different types of scalar KK modes induce mixing in the mass matrix of the KK spectrum, which leads to oscillation phenomena among these modes. Consequently, physical states observed on the brane are mass eigenstates arising from this mixing although the scalars acquire masses from the geometry of the extra dimensions. We compute the physical masses of these mixed modes in two illustrative six-dimensional (6D) brane models by diagonalizing the mass matrix.

        The remainder of this manuscript is organized as follows: In Sec. II, serving as a demonstrative example of the effect of the new term $(\nabla^M X_M)^2$, we investigate the gauge invariance and mass spectrum of the KK modes in a five-dimensional (5D) brane model. In Sec. III, we extend the study to a 6D brane, where two types of scalar KK modes emerge, leading to interesting oscillation phenomena. Numerical results related to these oscillations are presented in Sec. IV. Finally, a brief summary and discussion are provided in Sec. V.

      • I.   INTRODUCTION
        • The idea that our Universe may contain additional spatial dimensions beyond the four-dimensional spacetime has gained increasing attention over the past decades. The earliest attempt in this direction can be traced back to the Kaluza–Klein (KK) theory, which proposes that introducing a fifth spatial dimension enables a unified description of gravity and electromagnetism [13]. Subsequently, this idea was extended dramatically in string theory, wherein the presence of multiple extra dimensions is not optional but required for mathematical consistency. Building on these developments, modern brane-world models provided further insights into the hierarchy and cosmological constant problems. Representative examples include the Arkani-Hamed–Dimopoulos–Dvali model [4], which assumes large but compact extra dimensions, and the Randall–Sundrum scenarios [5, 6], where warped extra dimensions induce novel gravitational dynamics.

          Within the brane-world framework [716], one compelling line of research concerns the behavior of higher-dimensional matter fields and their associated KK modes [1727]. Investigating these KK spectra is crucial not only for ensuring the theoretical consistency of higher-dimensional models but also for probing the potential signatures of extra dimensions in observable physics. A growing body of literature has explored how KK excitations can lead to distinctive phenomenology, including modifications to gravitational interactions [2830], corrections to standard model processes [31, 32], and potential resonant modes localized near the brane [3335]. These studies highlight the importance of KK modes as a window into the possible structure of higher-dimensional spacetime [3639].

          Massless fields propagating in the bulk give rise to massive KK modes in the effective four-dimensional (4D) theory on the brane, with masses induced by the geometry of extra dimensions. For a bulk U(1) gauge field, the KK tower on the brane includes both vector and scalar modes. The effective 4D action preserves gauge invariance caused by vector–scalar and scalar–scalar couplings among the KK modes [4042]. Consequently, sufficient gauge freedom remains in the brane theory to eliminate the unphysical polarization degrees of freedom from the massive vector KK modes. This mechanism reveals that massive scalar KK modes play a role analogous to Stückelberg or Goldstone fields.

          Previous studies indicated that certain mixed vector–scalar and scalar–scalar couplings may exist. However, the precise form of these couplings remains uncertain because the basis functions for the scalar KK modes cannot be determined uniquely. Owing to such mixings, scalar KK modes cannot be fully gauged, which implies that they may carry physical degrees of freedom [43]. To properly define these scalar modes, and motivated by the Stückelberg mechanism, we propose adding a covariant term $(\nabla^M X_M)^2$ to the bulk action. This term introduces a kinetic contribution for the scalar involving $\partial_z X_z$, which helps clarify how the full geometry affects the scalar KK modes with the existing term. Meanwhile, this addition preserves gauge symmetry in the bulk and on the brane, similar to the Stückelberg mechanism.

          In a brane model with d extra dimensions ($d \gt 1$), there exist d distinct types of scalar KK modes. Interactions among the different types of scalar KK modes induce mixing in the mass matrix of the KK spectrum, which leads to oscillation phenomena among these modes. Consequently, physical states observed on the brane are mass eigenstates arising from this mixing although the scalars acquire masses from the geometry of the extra dimensions. We compute the physical masses of these mixed modes in two illustrative six-dimensional (6D) brane models by diagonalizing the mass matrix.

          The remainder of this manuscript is organized as follows: In Sec. II, serving as a demonstrative example of the effect of the new term $(\nabla^M X_M)^2$, we investigate the gauge invariance and mass spectrum of the KK modes in a five-dimensional (5D) brane model. In Sec. III, we extend the study to a 6D brane, where two types of scalar KK modes emerge, leading to interesting oscillation phenomena. Numerical results related to these oscillations are presented in Sec. IV. Finally, a brief summary and discussion are provided in Sec. V.

        II.   GAUGE-INVARIANT MASSIVE VECTORS AND THEIR ACCOMPANYING SCALARS IN THE 5D BRANE MODEL
        • We consider a U(1) gauge field coupling with a dilaton field π in a $ (4+1) $-dimensional braneworld scenario described by a conformal metric with the line element ${\rm d}s^2 = {\rm e}^{2A (z)}(\hat{g}_{\mu\nu}{\rm d}x^{\mu}{\rm d}x^{\nu} + {\rm d}z^2)$, where z represents the extra-dimensional coordinate and the warp factor ${\rm e}^{2A(z)}$depends only on z. The action for the bulk U(1) gauge field is given by

          $ S_1 =- \frac{1}{4} \int {\rm d}^5 x \sqrt{-g}\; {\rm e}^{\tau \pi}\big(F_{\mu\nu}F^{\mu\nu}+2F_{\mu z}F^{\mu z}+(\nabla_M X^M)^2\big), $

          (1)

          where τ represents the coupling constant. This action is invariant under the gauge transformations

          $ X_\mu\rightarrow X_\mu+\partial_\mu \Xi,\; \; X_z\rightarrow X_z+\partial_z \Xi, $

          (2)

          where Ξ represents a bulk scalar field satisfying the condition

          $ \partial_\mu \partial^\mu\Xi+\partial_z^2 \Xi+3(\partial_z \Xi)(\partial_z A)=0. $

          (3)

          Eq. (3) is a linear differential equation, and therefore, nontrivial solutions exist generically for regular warp factors and appropriate boundary conditions. Although the term $ (\nabla_M X^M)^2 $ is formally reminiscent of a gauge-fixing structure, in our framework, it is introduced at the level of the bulk action and leads to genuine physical effects after KK reduction instead of serving as a mere gauge choice.

          To derive the effective brane action, we perform a KK decomposition for the bulk field.

          $ X_\mu = \sum\limits_n \hat{X}^{(n)}_\mu(x^\nu) f^{(n)}(z), $

          (4a)

          $ X_z = \sum\limits_m \phi^{(m)}(x^\nu) \rho^{(m)}(z). $

          (4b)

          Here, $\hat{X}^{(n)}_\mu$ represents the n-level vector KK mode on the brane, and $\phi^{(m)}$ represents the m-level scalar KK modes on the brane. It is assumed that the wavefunction profiles $ \{f^{(n)}(z)\} $ and $ \{\rho^{(m)}(z)\} $ form a complete set along the extra dimension. We highlight three key consequences within the framework of the new bulk action (1).

          ● Gauge invariance of massive vector KK mode

          Substituting the KK decompositions (4) into the bulk action, we arrive at

          $ \begin{aligned}[b] S_{\text{eff}} =\;& -\frac{1}{4} \int {\rm d}^4 x \sqrt{|g|} \sum\limits_{n,m}\big( N^{(nm)}\;\hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\&+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I^{(nm)} + \widetilde{N}^{(nm)}\;\partial_\nu \phi^{(n)}\partial^\nu \phi^{(m)} \\ &-2\widetilde{I}^{(nm)} \hat{X}_\nu^{(n)}\partial^\nu \phi^{(m)} \big],+\big[ (\partial^\mu \hat{X}_\mu^{(n)})(\partial^\mu \hat{X}_\mu^{(m)})N^{(nm)} \\ &+ \phi^{(n)}\phi^{(m)}C^{(nm)} +2\widetilde{C}^{(nm)} \partial^\mu \hat{X}^{(n)}_\mu \phi^{(m)}\big]\big), \end{aligned} $

          (5)

          where

          $\begin{aligned}[b]& N^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; f^{(n)} f^{(m)},\\& \widetilde{N}^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; \rho^{(n)} \rho^{(m)},\end{aligned} $

          (6)

          $\begin{aligned}[b]& I^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; \partial_z f^{(n)} \partial_z f^{(m)},\\&\widetilde{I}^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; \partial_z f^{(n)} \rho^{(m)},\end{aligned} $

          (7)

          $ \begin{aligned}[b]C^{(nm)} =\;&\int {\rm e}^{A+\tau\pi} {\rm d}z\;[\partial_z\rho^{(n)}+3(\partial_z A)\rho^{(n)}]\\&\times[\partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}],\end{aligned} $

          (8)

          $ \widetilde{C}^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; f^{(n)}[\partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}]. $

          (9)

          Orthonormality conditions $ N^{(nn)}=1 $ and $ \widetilde{N}^{(mm)}=1 $ serve as the definitions of the inner product for the basis function sets $ \{ f^{(n)}(z) \} $ and $ \{ \rho^{(m)}(z) \} $, respectively. Consequently, from the definitions of $ \widetilde{I}^{(nm)} $ and $ \widetilde{C}^{(nm)} $, we obtain that

          $ \partial_z f^{(n)}=\sum\limits_k \widetilde{I}^{(nk)} \rho^{(k)},\; \; \; \; \; \partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}=\sum\limits_l \widetilde{C}^{(ml)} f^{(l)}, $

          (10)

          which leads to

          $ I^{(nm)} = \sum\limits_{k,l}\widetilde{I}^{(nk)} \widetilde{I}^{(ml)}\widetilde{N}^{(kl)},\; \; \; \; \; \widetilde{I}^{(nm)}=\sum\limits_k \widetilde{I}^{(nk)} \widetilde{N}^{(km)}, $

          (11)

          $ C^{(nm)} = \sum\limits_{k,l}\widetilde{C}^{(kn)} \widetilde{C}^{(lm)}N^{(km)},\; \; \; \; \widetilde{C}^{(nm)}=\sum\limits_k \widetilde{C}^{(nk)} N^{(km)}. $

          (12)

          Then, we can rewrite the effective action as

          $ \begin{aligned}[b] S_{\text{eff}} = \;&-\frac{1}{4} \sum\limits_{n,m} \int {\rm d}^4 x \sqrt{|g|} \big( \widetilde{N}_1^{(nm)} \hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\& +\widetilde{N}_1^{(nm)} \big[\partial_\nu \phi^{(n)}-\sum\limits_k \widetilde{I}^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \phi^{(m)}-\sum\limits_l \widetilde{I}^{(lm)} \hat{X}^{\nu(l)}\big],\\ & + N_1^{(nm)} \big[\partial^\mu \hat{X}_\mu^{(n)}+\sum\limits_{k} \widetilde{C}^{(nk)} \phi^{(k)}\big] \big[\partial^\mu \hat{X}_\mu^{(m)}+ \sum\limits_{l} \widetilde{C}^{(ml)} \phi^{(l)}\big]. \end{aligned} $

          (13)

          This is gauge invariant under gauge transformations:

          $ \hat{X}_\nu^{(k)} \rightarrow \hat{X}_\nu^{(k)} + \partial_\nu \Lambda^{(k)}, $

          (14)

          $ \phi^{(n)} \rightarrow \phi^{(n)} + \sum\limits_k \widetilde{I}^{(kn)} \Lambda^{(k)}, $

          (15)

          where $ \Lambda^{(k)} $ represents the brane-localized gauge freedom and satisfies

          $ \partial^\nu\partial_\nu \Lambda^{(n)}+\sum\limits_{k,l}(\widetilde{C}^{(nk)}\widetilde{I}^{(lk)} )\Lambda^{(l)}=0. $

          (16)

          This condition originates from the bulk constraint Eq. (3). By applying the separation of variables $ \Xi(x,z) = \sum\limits_n \Lambda^{(n)}(x)\Theta^{(n)}(z) $, the bulk constraint decouples into a linear ordinary differential equation for the extra-dimensional profile $ \Theta^{(n)}(z) $ and effective 4D equation (16). The ODE admits well-defined solutions for regular backgrounds, and therefore, the existence of nontrivial solutions for the bulk parameter $ \Xi(x,z) $ is guaranteed if Eq. (16) is solvable. In the mass eigenbasis, Eq. (16) reduces to a standard massive Klein-Gordon equation, which naturally admits propagating wave solutions. This confirms that the residual gauge symmetry is physically realizable.

          Effective action (13) exhibits a structure reminiscent of the Stueckelberg mechanism; however, it possesses distinct physical origins and implications. Unlike the standard Stueckelberg formalism where the gauge parameter is an arbitrary function, the brane-localized gauge freedom $ \Lambda^{(n)} $ is constrained by the differential condition in Eq. (16). This restriction implies that a residual gauge symmetry is preserved while the full gauge redundancy is broken.

          ● Eigenstate of vector and scalar KK modes

          Starting from the original definitions of $ I^{(nn)} $ and $ C^{(nn)} $, and imposing either Dirichlet or periodic boundary conditions on the basis functions $ f^{(n)} $ and $ \rho^{(m)} $, we derive two constraint equations for the massive vector and scalar KK modes:

          $ -{\rm e}^{-(A+\tau\pi)} \partial_z ({\rm e}^{(A+\tau\pi)}\partial_z f^{(m)}) = m_v^{2} f^{(m)}, $

          (17)

          $ -{\rm e}^{2A-\tau\pi}\partial_z [{\rm e}^{\tau\pi-5A}\partial_z({\rm e}^{3A}\rho^{(m)})] = m_{\phi}^{2} \rho^{(m)}. $

          (18)

          After applying the transformation $f^{(n)}={\rm e}^{-\frac{\sqrt{3}\tau+1}{2}A} \bar{f}^{(n)}$ and $\rho^{(m)}={\rm e}^{-\frac{3}{2} A} \bar{\rho}^{(m)}$, the equations reduce to the following Schrödinger-like forms:

          $ (-\partial^2_z+V_{\text{vec}})\bar{f}^{(n)}=m_v^{2} \bar{f}^{(n)}, $

          (19)

          $ (-\partial^2_z+V_{\text{sca}})\bar{\rho}^{(m)}=m_\phi^{2} \bar{\rho}^{(m)}, $

          (20)

          with the effective potentials given by

          $ V_{\text{vec}}=\frac{(\sqrt{3}\tau+1)^2}{4}A'^2+\frac{\sqrt{3}\tau+1}{2}A'' , $

          (21)

          $ V_{\text{sca}}=\frac{(\sqrt{3} \tau-5)^2 }{4}A^{\prime 2}+\frac{\sqrt{3} \tau-5}{2} A^{\prime \prime}. $

          (22)

          Here, we identify the dilaton field with $ \pi=\sqrt{3}A $ following the literature [44], where the brane solution is given as

          $ A(z)=-\frac{v^2}{9}\big(\ln \cosh ^2(a z)+\frac{1}{2} \tanh ^2(a z)\big). $

          (23)

          We analyze the behavior of the potentials for both vector and scalar KK modes in the limits $ z\rightarrow 0 $ and $ z \rightarrow \infty $, obtaining:

          $\begin{aligned}[b]& V_{\text{vec}}(z\rightarrow 0) \rightarrow -\frac{1}{6} a^2 v^2(1+\sqrt{3} \tau),\\& V_{\text{vec}}(z\rightarrow \infty)\rightarrow \frac{1}{81} a^2v^4(1+\sqrt{3} \tau)^2,\end{aligned} $

          (24)

          $\begin{aligned}[b]& V_{\text{sca}}(z\rightarrow 0) \rightarrow -\frac{1}{6} a^2 v^2(\sqrt{3} \tau-5),\\& V_{\text{sca}}(z\rightarrow \infty)\rightarrow \frac{1}{81} a^2v^4(\sqrt{3} \tau-5)^2. \end{aligned}$

          (25)

          It is evident that, for $ \tau>-1/\sqrt{3} $, a zero vector KK mode is always localized, while a massless scalar mode only appears when $ \tau>5/\sqrt{3} $. When τ is sufficiently large, both massive vector and scalar modes exist. The corresponding effective potentials are illustrated in Fig. 1. Further, as shown in Table 1, the masses of the scalar modes are smaller than those of the vector modes.

          Figure 1.  (color online) Effective potentials for the vector and scalar KK modes with $ a=1,v=1 $.

          $ n_2 $ $ \tau=1.15 $ $ \tau=3.5 $ $ \tau=30 $ $ \tau=40 $
          $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $
          0 0 0 0 0 0 0 0
          1 15.59 13.59 21.38 19.38
          2 26.89 22.82 38.56 34.53
          3 33.41 27.09 51.24 45.11
          4 58.96 50.54

          Table 1.  Squared mass eigenvalues for vector and scalar KK modes.

        II.   GAUGE-INVARIANT MASSIVE VECTORS AND THEIR ACCOMPANYING SCALARS IN 5D BRANE MODEL
        • We consider a U(1) gauge field coupling with a dilaton field π in a $ (4+1) $-dimensional braneworld scenario described by a conformal metric with the line element $ ds^2 = e^{2A (z)}(\hat{g}_{\mu\nu}dx^{\mu}dx^{\nu} + dz^2) $, where z denotes the extra-dimensional coordinate and the warp factor $ e^{2A(z)} $depends only on z. The action for the bulk U(1) gauge field is taken to be

          $ S_1 =- \frac{1}{4} \int d^5 x \sqrt{-g}\; e^{\tau \pi}\big(F_{\mu\nu}F^{\mu\nu}+2F_{\mu z}F^{\mu z}+(\nabla_M X^M)^2\big), $

          (1)

          where τ is the coupling constant. This action is invariant under the gauge transformations

          $ X_\mu\rightarrow X_\mu+\partial_\mu \Xi,\; \; X_z\rightarrow X_z+\partial_z \Xi, $

          (2)

          where Ξ is a bulk scalar field satisfying the condition

          $ \partial_\mu \partial^\mu\Xi+\partial_z^2 \Xi+3(\partial_z \Xi)(\partial_z A)=0. $

          (3)

          Since Eq.(3) is a linear differential equation, nontrivial solutions generically exist for regular warp factors and appropriate boundary conditions. Although the term $ (\nabla_M X^M)^2 $ is formally reminiscent of a gauge-fixing structure, in our framework it is introduced at the level of the bulk action and leads to genuine physical effects after KK reduction, rather than serving as a mere gauge choice.

          To derive the effective brane action, we perform a KK decomposition for the bulk field:

          $ X_\mu = \sum\limits_n \hat{X}^{(n)}_\mu(x^\nu) f^{(n)}(z), $

          (4a)

          $ X_z = \sum\limits_m \phi^{(m)}(x^\nu) \rho^{(m)}(z). $

          (4b)

          Here, $\hat{X}^{(n)}_\mu$ represents the n-level vector KK mode on the brane, and $\phi^{(m)}$ represent the m-level scalar KK modes on the brane, respectively. It is assumed that the wavefunction profiles $ \{f^{(n)}(z)\} $ and $ \{\rho^{(m)}(z)\} $ each form a complete set along the extra dimension. Here we highlight three key consequences within the framework of the new bulk action (1).

          ● The gauge invariance of massive vector KK mode

          Substituting the KK decompositions (4) into the bulk action, we arrive at the following effective action:

          $ \begin{aligned}[b] S_{\text{eff}} =\;& -\frac{1}{4} \int d^4 x \sqrt{|g|} \sum\limits_{n,m}\big( N^{(nm)}\;\hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\&+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I^{(nm)} + \widetilde{N}^{(nm)}\;\partial_\nu \phi^{(n)}\partial^\nu \phi^{(m)} \\ &-2\widetilde{I}^{(nm)} \hat{X}_\nu^{(n)}\partial^\nu \phi^{(m)} \big],+\big[ (\partial^\mu \hat{X}_\mu^{(n)})(\partial^\mu \hat{X}_\mu^{(m)})N^{(nm)} \\ &+ \phi^{(n)}\phi^{(m)}C^{(nm)} +2\widetilde{C}^{(nm)} \partial^\mu \hat{X}^{(n)}_\mu \phi^{(m)}\big]\big), \end{aligned} $

          (5)

          where:

          $\begin{aligned}[b]& N^{(nm)} = \int e^{A+\tau\pi} dz\; f^{(n)} f^{(m)},\\& \widetilde{N}^{(nm)} = \int e^{A+\tau\pi} dz\; \rho^{(n)} \rho^{(m)},\end{aligned} $

          (6)

          $\begin{aligned}[b]& I^{(nm)} = \int e^{A+\tau\pi} dz\; \partial_z f^{(n)} \partial_z f^{(m)},\\&\widetilde{I}^{(nm)} = \int e^{A+\tau\pi} dz\; \partial_z f^{(n)} \rho^{(m)},\end{aligned} $

          (7)

          $ \begin{aligned}[b]C^{(nm)} =\;&\int e^{A+\tau\pi} dz\;[\partial_z\rho^{(n)}+3(\partial_z A)\rho^{(n)}]\\&\times[\partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}],\end{aligned} $

          (8)

          $ \widetilde{C}^{(nm)} = \int e^{A+\tau\pi} dz\; f^{(n)}[\partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}]. $

          (9)

          The orthonormality conditions, $ N^{(nn)}=1 $ and $ \widetilde{N}^{(mm)}=1 $, serve as the definitions of the inner product for the basis function sets $ \{ f^{(n)}(z) \} $ and $ \{ \rho^{(m)}(z) \} $, respectively. Consequently, from the definitions of $ \widetilde{I}^{(nm)} $ and $ \widetilde{C}^{(nm)} $, we obtain that

          $ \partial_z f^{(n)}=\sum\limits_k \widetilde{I}^{(nk)} \rho^{(k)},\; \; \; \; \; \partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}=\sum\limits_l \widetilde{C}^{(ml)} f^{(l)}, $

          (10)

          which lead to

          $ I^{(nm)} = \sum\limits_{k,l}\widetilde{I}^{(nk)} \widetilde{I}^{(ml)}\widetilde{N}^{(kl)},\; \; \; \; \; \widetilde{I}^{(nm)}=\sum\limits_k \widetilde{I}^{(nk)} \widetilde{N}^{(km)}, $

          (11)

          $ C^{(nm)} = \sum\limits_{k,l}\widetilde{C}^{(kn)} \widetilde{C}^{(lm)}N^{(km)},\; \; \; \; \widetilde{C}^{(nm)}=\sum\limits_k \widetilde{C}^{(nk)} N^{(km)}. $

          (12)

          Then we can rewrite the effective action as

          $ \begin{aligned}[b] S_{\text{eff}} = \;&-\frac{1}{4} \sum\limits_{n,m} \int d^4 x \sqrt{|g|} \big( \widetilde{N}_1^{(nm)} \hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\& +\widetilde{N}_1^{(nm)} \big[\partial_\nu \phi^{(n)}-\sum\limits_k \widetilde{I}^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \phi^{(m)}-\sum\limits_l \widetilde{I}^{(lm)} \hat{X}^{\nu(l)}\big],\\ & + N_1^{(nm)} \big[\partial^\mu \hat{X}_\mu^{(n)}+\sum\limits_{k} \widetilde{C}^{(nk)} \phi^{(k)}\big] \big[\partial^\mu \hat{X}_\mu^{(m)}+ \sum\limits_{l} \widetilde{C}^{(ml)} \phi^{(l)}\big]. \end{aligned} $

          (13)

          which is gauge invariant under the gauge transformations:

          $ \hat{X}_\nu^{(k)} \rightarrow \hat{X}_\nu^{(k)} + \partial_\nu \Lambda^{(k)}, $

          (14)

          $ \phi^{(n)} \rightarrow \phi^{(n)} + \sum\limits_k \widetilde{I}^{(kn)} \Lambda^{(k)}, $

          (15)

          Here $ \Lambda^{(k)} $ the brane-localized gauge freedom satisfies

          $ \partial^\nu\partial_\nu \Lambda^{(n)}+\sum\limits_{k,l}(\widetilde{C}^{(nk)}\widetilde{I}^{(lk)} )\Lambda^{(l)}=0. $

          (16)

          This condition originates from the bulk constraint Eq. (3). By applying the separation of variables $ \Xi(x,z) = \sum\limits_n \Lambda^{(n)}(x)\Theta^{(n)}(z) $, the bulk constraint decouples into a linear ordinary differential equation for the extra-dimensional profile $ \Theta^{(n)}(z) $ and the effective 4D equation (16). Since the ODE admits well-defined solutions for regular backgrounds, the existence of nontrivial solutions for the bulk parameter $ \Xi(x,z) $ is guaranteed if Eq. (16) is solvable. In the mass eigenbasis, as we will see, Eq. (16) reduces to a standard massive Klein-Gordon equation, which naturally admits propagating wave solutions. This confirms that the residual gauge symmetry is physically realizable.

          The effective action (13) exhibits a structure reminiscent of the Stueckelberg mechanism, yet it possesses distinct physical origins and implications. Unlike the standard Stueckelberg formalism where the gauge parameter is an arbitrary function, the brane-localized gauge freedom $ \Lambda^{(n)} $ is constrained by the differential condition in Eq. (16). This restriction implies that while the full gauge redundancy is broken, a residual gauge symmetry is preserved.

          ● Eigenstate of vector and scalar KK modes

          Starting from the original definitions of $ I^{(nn)} $ and $ C^{(nn)} $, and imposing either Dirichlet or periodic boundary conditions on the basis functions $ f^{(n)} $ and $ \rho^{(m)} $, we derive two constraint equations for the massive vector and scalar KK modes:

          $ -e^{-(A+\tau\pi)} \partial_z (e^{(A+\tau\pi)}\partial_z f^{(m)}) = m_v^{2} f^{(m)}, $

          (17)

          $ -e^{2A-\tau\pi}\partial_z [e^{\tau\pi-5A}\partial_z(e^{3A}\rho^{(m)})] = m_{\phi}^{2} \rho^{(m)}. $

          (18)

          After applying the transformation $ f^{(n)}=e^{-\frac{\sqrt{3}\tau+1}{2}A} \bar{f}^{(n)} $ and $ \rho^{(m)}=e^{-\frac{3}{2} A} \bar{\rho}^{(m)} $, the equations reduce to the following Schrödinger-like forms:

          $ (-\partial^2_z+V_{\text{vec}})\bar{f}^{(n)}=m_v^{2} \bar{f}^{(n)}, $

          (19)

          $ (-\partial^2_z+V_{\text{sca}})\bar{\rho}^{(m)}=m_\phi^{2} \bar{\rho}^{(m)}, $

          (20)

          with the effective potentials given by

          $ V_{\text{vec}}=\frac{(\sqrt{3}\tau+1)^2}{4}A'^2+\frac{\sqrt{3}\tau+1}{2}A'' , $

          (21)

          $ V_{\text{sca}}=\frac{(\sqrt{3} \tau-5)^2 }{4}A^{\prime 2}+\frac{\sqrt{3} \tau-5}{2} A^{\prime \prime}. $

          (22)

          Here we identify the dilaton field with $ \pi=\sqrt{3}A $ following the literature [47], where the brane solution is given as

          $ A(z)=-\frac{v^2}{9}\big(\ln \cosh ^2(a z)+\frac{1}{2} \tanh ^2(a z)\big). $

          (23)

          We analyze the behavior of the potentials for both vector and scalar KK modes in the limits $ z\rightarrow 0 $ and $ z \rightarrow \infty $, obtaining:

          $\begin{aligned}[b]& V_{\text{vec}}(z\rightarrow 0) \rightarrow -\frac{1}{6} a^2 v^2(1+\sqrt{3} \tau),\\& V_{\text{vec}}(z\rightarrow \infty)\rightarrow \frac{1}{81} a^2v^4(1+\sqrt{3} \tau)^2,\end{aligned} $

          (24)

          $\begin{aligned}[b]& V_{\text{sca}}(z\rightarrow 0) \rightarrow -\frac{1}{6} a^2 v^2(\sqrt{3} \tau-5),\\& V_{\text{sca}}(z\rightarrow \infty)\rightarrow \frac{1}{81} a^2v^4(\sqrt{3} \tau-5)^2. \end{aligned}$

          (25)

          It is evident that for $ \tau>-1/\sqrt{3} $, a zero vector KK mode is always localized, while a massless scalar mode only appears when $ \tau>5/\sqrt{3} $. When τ is sufficiently large, both massive vector and scalar modes exist. The corresponding effective potentials are illustrated in Fig. 1. Moreover, as shown in Table 1, the masses of the scalar modes are smaller than those of the vector modes.

          Figure 1.  (color online) The effective potentials for the vector and scalar KK modes with $ a=1,v=1 $.

          $ n_2 $ $ \tau=1.15 $ $ \tau=3.5 $ $ \tau=30 $ $ \tau=40 $
          $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $
          0 0 0 0 0 0 0 0
          1 15.59 13.59 21.38 19.38
          2 26.89 22.82 38.56 34.53
          3 33.41 27.09 51.24 45.11
          4 58.96 50.54

          Table 1.  Squared mass eigenvalues for vector and scalar KK modes.

        II.   GAUGE-INVARIANT MASSIVE VECTORS AND THEIR ACCOMPANYING SCALARS IN THE 5D BRANE MODEL
        • We consider a U(1) gauge field coupling with a dilaton field π in a $ (4+1) $-dimensional braneworld scenario described by a conformal metric with the line element ${\rm d}s^2 = {\rm e}^{2A (z)}(\hat{g}_{\mu\nu}{\rm d}x^{\mu}{\rm d}x^{\nu} + {\rm d}z^2)$, where z represents the extra-dimensional coordinate and the warp factor ${\rm e}^{2A(z)}$depends only on z. The action for the bulk U(1) gauge field is given by

          $ S_1 =- \frac{1}{4} \int {\rm d}^5 x \sqrt{-g}\; {\rm e}^{\tau \pi}\big(F_{\mu\nu}F^{\mu\nu}+2F_{\mu z}F^{\mu z}+(\nabla_M X^M)^2\big), $

          (1)

          where τ represents the coupling constant. This action is invariant under the gauge transformations

          $ X_\mu\rightarrow X_\mu+\partial_\mu \Xi,\; \; X_z\rightarrow X_z+\partial_z \Xi, $

          (2)

          where Ξ represents a bulk scalar field satisfying the condition

          $ \partial_\mu \partial^\mu\Xi+\partial_z^2 \Xi+3(\partial_z \Xi)(\partial_z A)=0. $

          (3)

          Eq. (3) is a linear differential equation, and therefore, nontrivial solutions exist generically for regular warp factors and appropriate boundary conditions. Although the term $ (\nabla_M X^M)^2 $ is formally reminiscent of a gauge-fixing structure, in our framework, it is introduced at the level of the bulk action and leads to genuine physical effects after KK reduction instead of serving as a mere gauge choice.

          To derive the effective brane action, we perform a KK decomposition for the bulk field.

          $ X_\mu = \sum\limits_n \hat{X}^{(n)}_\mu(x^\nu) f^{(n)}(z), $

          (4a)

          $ X_z = \sum\limits_m \phi^{(m)}(x^\nu) \rho^{(m)}(z). $

          (4b)

          Here, $\hat{X}^{(n)}_\mu$ represents the n-level vector KK mode on the brane, and $\phi^{(m)}$ represents the m-level scalar KK modes on the brane. It is assumed that the wavefunction profiles $ \{f^{(n)}(z)\} $ and $ \{\rho^{(m)}(z)\} $ form a complete set along the extra dimension. We highlight three key consequences within the framework of the new bulk action (1).

          ● Gauge invariance of massive vector KK mode

          Substituting the KK decompositions (4) into the bulk action, we arrive at

          $ \begin{aligned}[b] S_{\text{eff}} =\;& -\frac{1}{4} \int {\rm d}^4 x \sqrt{|g|} \sum\limits_{n,m}\big( N^{(nm)}\;\hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\&+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I^{(nm)} + \widetilde{N}^{(nm)}\;\partial_\nu \phi^{(n)}\partial^\nu \phi^{(m)} \\ &-2\widetilde{I}^{(nm)} \hat{X}_\nu^{(n)}\partial^\nu \phi^{(m)} \big],+\big[ (\partial^\mu \hat{X}_\mu^{(n)})(\partial^\mu \hat{X}_\mu^{(m)})N^{(nm)} \\ &+ \phi^{(n)}\phi^{(m)}C^{(nm)} +2\widetilde{C}^{(nm)} \partial^\mu \hat{X}^{(n)}_\mu \phi^{(m)}\big]\big), \end{aligned} $

          (5)

          where

          $\begin{aligned}[b]& N^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; f^{(n)} f^{(m)},\\& \widetilde{N}^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; \rho^{(n)} \rho^{(m)},\end{aligned} $

          (6)

          $\begin{aligned}[b]& I^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; \partial_z f^{(n)} \partial_z f^{(m)},\\&\widetilde{I}^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; \partial_z f^{(n)} \rho^{(m)},\end{aligned} $

          (7)

          $ \begin{aligned}[b]C^{(nm)} =\;&\int {\rm e}^{A+\tau\pi} {\rm d}z\;[\partial_z\rho^{(n)}+3(\partial_z A)\rho^{(n)}]\\&\times[\partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}],\end{aligned} $

          (8)

          $ \widetilde{C}^{(nm)} = \int {\rm e}^{A+\tau\pi} {\rm d}z\; f^{(n)}[\partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}]. $

          (9)

          Orthonormality conditions $ N^{(nn)}=1 $ and $ \widetilde{N}^{(mm)}=1 $ serve as the definitions of the inner product for the basis function sets $ \{ f^{(n)}(z) \} $ and $ \{ \rho^{(m)}(z) \} $, respectively. Consequently, from the definitions of $ \widetilde{I}^{(nm)} $ and $ \widetilde{C}^{(nm)} $, we obtain that

          $ \partial_z f^{(n)}=\sum\limits_k \widetilde{I}^{(nk)} \rho^{(k)},\; \; \; \; \; \partial_z\rho^{(m)}+3(\partial_z A)\rho^{(m)}=\sum\limits_l \widetilde{C}^{(ml)} f^{(l)}, $

          (10)

          which leads to

          $ I^{(nm)} = \sum\limits_{k,l}\widetilde{I}^{(nk)} \widetilde{I}^{(ml)}\widetilde{N}^{(kl)},\; \; \; \; \; \widetilde{I}^{(nm)}=\sum\limits_k \widetilde{I}^{(nk)} \widetilde{N}^{(km)}, $

          (11)

          $ C^{(nm)} = \sum\limits_{k,l}\widetilde{C}^{(kn)} \widetilde{C}^{(lm)}N^{(km)},\; \; \; \; \widetilde{C}^{(nm)}=\sum\limits_k \widetilde{C}^{(nk)} N^{(km)}. $

          (12)

          Then, we can rewrite the effective action as

          $ \begin{aligned}[b] S_{\text{eff}} = \;&-\frac{1}{4} \sum\limits_{n,m} \int {\rm d}^4 x \sqrt{|g|} \big( \widetilde{N}_1^{(nm)} \hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\& +\widetilde{N}_1^{(nm)} \big[\partial_\nu \phi^{(n)}-\sum\limits_k \widetilde{I}^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \phi^{(m)}-\sum\limits_l \widetilde{I}^{(lm)} \hat{X}^{\nu(l)}\big],\\ & + N_1^{(nm)} \big[\partial^\mu \hat{X}_\mu^{(n)}+\sum\limits_{k} \widetilde{C}^{(nk)} \phi^{(k)}\big] \big[\partial^\mu \hat{X}_\mu^{(m)}+ \sum\limits_{l} \widetilde{C}^{(ml)} \phi^{(l)}\big]. \end{aligned} $

          (13)

          This is gauge invariant under gauge transformations:

          $ \hat{X}_\nu^{(k)} \rightarrow \hat{X}_\nu^{(k)} + \partial_\nu \Lambda^{(k)}, $

          (14)

          $ \phi^{(n)} \rightarrow \phi^{(n)} + \sum\limits_k \widetilde{I}^{(kn)} \Lambda^{(k)}, $

          (15)

          where $ \Lambda^{(k)} $ represents the brane-localized gauge freedom and satisfies

          $ \partial^\nu\partial_\nu \Lambda^{(n)}+\sum\limits_{k,l}(\widetilde{C}^{(nk)}\widetilde{I}^{(lk)} )\Lambda^{(l)}=0. $

          (16)

          This condition originates from the bulk constraint Eq. (3). By applying the separation of variables $ \Xi(x,z) = \sum\limits_n \Lambda^{(n)}(x)\Theta^{(n)}(z) $, the bulk constraint decouples into a linear ordinary differential equation for the extra-dimensional profile $ \Theta^{(n)}(z) $ and effective 4D equation (16). The ODE admits well-defined solutions for regular backgrounds, and therefore, the existence of nontrivial solutions for the bulk parameter $ \Xi(x,z) $ is guaranteed if Eq. (16) is solvable. In the mass eigenbasis, Eq. (16) reduces to a standard massive Klein-Gordon equation, which naturally admits propagating wave solutions. This confirms that the residual gauge symmetry is physically realizable.

          Effective action (13) exhibits a structure reminiscent of the Stueckelberg mechanism; however, it possesses distinct physical origins and implications. Unlike the standard Stueckelberg formalism where the gauge parameter is an arbitrary function, the brane-localized gauge freedom $ \Lambda^{(n)} $ is constrained by the differential condition in Eq. (16). This restriction implies that a residual gauge symmetry is preserved while the full gauge redundancy is broken.

          ● Eigenstate of vector and scalar KK modes

          Starting from the original definitions of $ I^{(nn)} $ and $ C^{(nn)} $, and imposing either Dirichlet or periodic boundary conditions on the basis functions $ f^{(n)} $ and $ \rho^{(m)} $, we derive two constraint equations for the massive vector and scalar KK modes:

          $ -{\rm e}^{-(A+\tau\pi)} \partial_z ({\rm e}^{(A+\tau\pi)}\partial_z f^{(m)}) = m_v^{2} f^{(m)}, $

          (17)

          $ -{\rm e}^{2A-\tau\pi}\partial_z [{\rm e}^{\tau\pi-5A}\partial_z({\rm e}^{3A}\rho^{(m)})] = m_{\phi}^{2} \rho^{(m)}. $

          (18)

          After applying the transformation $f^{(n)}={\rm e}^{-\frac{\sqrt{3}\tau+1}{2}A} \bar{f}^{(n)}$ and $\rho^{(m)}={\rm e}^{-\frac{3}{2} A} \bar{\rho}^{(m)}$, the equations reduce to the following Schrödinger-like forms:

          $ (-\partial^2_z+V_{\text{vec}})\bar{f}^{(n)}=m_v^{2} \bar{f}^{(n)}, $

          (19)

          $ (-\partial^2_z+V_{\text{sca}})\bar{\rho}^{(m)}=m_\phi^{2} \bar{\rho}^{(m)}, $

          (20)

          with the effective potentials given by

          $ V_{\text{vec}}=\frac{(\sqrt{3}\tau+1)^2}{4}A'^2+\frac{\sqrt{3}\tau+1}{2}A'' , $

          (21)

          $ V_{\text{sca}}=\frac{(\sqrt{3} \tau-5)^2 }{4}A^{\prime 2}+\frac{\sqrt{3} \tau-5}{2} A^{\prime \prime}. $

          (22)

          Here, we identify the dilaton field with $ \pi=\sqrt{3}A $ following the literature [44], where the brane solution is given as

          $ A(z)=-\frac{v^2}{9}\big(\ln \cosh ^2(a z)+\frac{1}{2} \tanh ^2(a z)\big). $

          (23)

          We analyze the behavior of the potentials for both vector and scalar KK modes in the limits $ z\rightarrow 0 $ and $ z \rightarrow \infty $, obtaining:

          $\begin{aligned}[b]& V_{\text{vec}}(z\rightarrow 0) \rightarrow -\frac{1}{6} a^2 v^2(1+\sqrt{3} \tau),\\& V_{\text{vec}}(z\rightarrow \infty)\rightarrow \frac{1}{81} a^2v^4(1+\sqrt{3} \tau)^2,\end{aligned} $

          (24)

          $\begin{aligned}[b]& V_{\text{sca}}(z\rightarrow 0) \rightarrow -\frac{1}{6} a^2 v^2(\sqrt{3} \tau-5),\\& V_{\text{sca}}(z\rightarrow \infty)\rightarrow \frac{1}{81} a^2v^4(\sqrt{3} \tau-5)^2. \end{aligned}$

          (25)

          It is evident that, for $ \tau>-1/\sqrt{3} $, a zero vector KK mode is always localized, while a massless scalar mode only appears when $ \tau>5/\sqrt{3} $. When τ is sufficiently large, both massive vector and scalar modes exist. The corresponding effective potentials are illustrated in Fig. 1. Further, as shown in Table 1, the masses of the scalar modes are smaller than those of the vector modes.

          Figure 1.  (color online) Effective potentials for the vector and scalar KK modes with $ a=1,v=1 $.

          $ n_2 $ $ \tau=1.15 $ $ \tau=3.5 $ $ \tau=30 $ $ \tau=40 $
          $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $
          0 0 0 0 0 0 0 0
          1 15.59 13.59 21.38 19.38
          2 26.89 22.82 38.56 34.53
          3 33.41 27.09 51.24 45.11
          4 58.96 50.54

          Table 1.  Squared mass eigenvalues for vector and scalar KK modes.

        III.   GAUGE-INVARIANT MASSIVE VECTORS AND THEIR ACCOMPANYING SCALARS IN 6D BRANE MODEL
        • We consider a massless U(1) gauge field in a $ (4+2) $-dimensional brane world. The line element for the brane is given by

          $ {\rm d}s^2 = {\rm e}^{2A (z,y)}\hat{g}_{\mu\nu} {\rm d}x^{\mu} {\rm d}x^{\nu} + {\rm e}^{2B_1(z,y)} {\rm d}z^2 +{\rm e}^{2B_2(z,y)} {\rm d}y^2, $

          (26)

          where the warp factors ${\rm e}^{2A(z,y)}$ and${\rm e}^{2B_2(z,y)},\; {\rm e}^{2B_2(z,y)}$ are functions of the two extra dimensions.

        III.   GAUGE-INVARIANT MASSIVE VECTORS AND THEIR ACCOMPANYING SCALARS IN 6D BRANE MODEL
        • We consider a massless U(1) gauge field in a $ (4+2) $-dimensional brane world. The line element for the brane is given by:

          $ ds^2 = e^{2A (z,y)}\hat{g}_{\mu\nu}dx^{\mu}dx^{\nu} + e^{2B_1(z,y)} dz^2 +e^{2B_2(z,y)} dy^2, $

          (26)

          where the warp factors $ e^{2A(z,y)} $ and $ e^{2B_2(z,y)},e^{2B_2(z,y)} $ are functions of the two extra dimensions.

        III.   GAUGE-INVARIANT MASSIVE VECTORS AND THEIR ACCOMPANYING SCALARS IN 6D BRANE MODEL
        • We consider a massless U(1) gauge field in a $ (4+2) $-dimensional brane world. The line element for the brane is given by

          $ {\rm d}s^2 = {\rm e}^{2A (z,y)}\hat{g}_{\mu\nu} {\rm d}x^{\mu} {\rm d}x^{\nu} + {\rm e}^{2B_1(z,y)} {\rm d}z^2 +{\rm e}^{2B_2(z,y)} {\rm d}y^2, $

          (26)

          where the warp factors ${\rm e}^{2A(z,y)}$ and${\rm e}^{2B_2(z,y)},\; {\rm e}^{2B_2(z,y)}$ are functions of the two extra dimensions.

        • A.   Review the gauge invariance of massive vector KK modes

        • For a free bulk U(1) gauge field in the 6D brane world, after performing the KK decomposition

          $\begin{aligned}[b] X_\mu &= \sum\limits_n \hat{X}^{(n)}_\mu(x^\nu) f^{(n)}(z,y), X_z \\&= \sum\limits_n \phi^{(n)}(x^\nu) \rho^{(n)}(z,y), X_y = \sum\limits_n \varphi^{(n)}(x^\nu) \chi^{(n)}(z,y),\end{aligned} $

          (27)

          where $\hat{X}^{(n)}_\mu$ represents the n-level vector KK mode and $\phi^{(m)}, \varphi^{(m)}$ represent the two types of scalar ones. The effective action takes the form

          $ \begin{aligned}[b] S_{\text{eff1}} =\;& -\frac{1}{4} \sum\limits_{n,m}\int {\rm d}^4 x \sqrt{|g|} \big( N^{(nm)}\;\hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\ &+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I_1^{(nm)} -2\hat{X}_\nu^{(n)}\partial^\nu \phi^{(m)}\widetilde{I}_1^{(nm)} \\ & + \widetilde{N}_1^{(nm)}\;\partial_\nu \phi^{(n)}\partial^\nu \phi^{(m)}\big]+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I_2^{(nm)} \\ & -2\hat{X}_\nu^{(n)}\partial^\nu \varphi^{(m)}\widetilde{I}_2^{(nm)} + \widetilde{N}_2^{(nm)}\;\partial_\nu \varphi^{(n)}\partial^\nu \varphi^{(m)}\big]\\ & +2\big[C_1^{(nm)}\phi^{(n)}\phi^{(m)}+C_2^{(nm)}\varphi^{(n)}\varphi^{(m)}\\ & -2\widetilde{C}^{(nm)}\phi^{(n)}\varphi^{(m)}\big] \big) \end{aligned} $

          (28)

          where

          $ N^{(nm)} = \int {\rm e}^{B_1+B_2} {\rm d}z {\rm d}y\; f^{(n)} f^{(m)}, $

          (29)

          $\begin{aligned}[b] & \widetilde{N}_1^{(nm)} = \int {\rm e}^{2A-B_1+B_2} {\rm d}z {\rm d}y\; \rho^{(n)} \rho^{(m)},\\&\widetilde{N}_2^{(nm)} = \int {\rm e}^{2A+B_1-B_2} {\rm d}z {\rm d}y\; \chi^{(n)} \chi^{(m)}, \end{aligned}$

          (30)

          $\begin{aligned}[b]& I_1^{(nm)} = \int {\rm e}^{2A-B_1+B_2} {\rm d}z{\rm d}y\; \partial_z f^{(n)} \partial_z f^{(m)},\\& \widetilde{I}_1^{(nm)} = \int {\rm e}^{2A-B_1+B_2} {\rm d}z{\rm d}y\; \partial_z f^{(n)} \rho^{(m)},\end{aligned} $

          (31)

          $\begin{aligned}[b]& I_2^{(nm)} = \int {\rm e}^{2A+B_1-B_2} {\rm d}z{\rm d}y\; \partial_y f^{(n)} \partial_y f^{(m)},\\& \widetilde{I}_2^{(nm)} = \int {\rm e}^{2A+B_1-B_2} {\rm d}z{\rm d}y\; \partial_y f^{(n)} \chi^{(m)},\end{aligned} $

          (32)

          $\begin{aligned}[b]& C_1^{(mm)} = \int {\rm e}^{4A-B_1-B_2} {\rm d}z{\rm d}y\; \partial_y \rho^{(m)} \partial_y \rho^{(m)} ,\\& C_2^{(kk)}= \int {\rm e}^{4A-B_1-B_2} {\rm d}z{\rm d}y\; \partial_z \chi^{(k)} \partial_z \chi^{(k)},\end{aligned} $

          (33)

          $ \widetilde{C}^{(mk)} = \int {\rm e}^{4A-B_1-B_2} {\rm d}z{\rm d}y\; \partial_y \rho^{(m)} \partial_z \chi^{(k)}. $

          (34)

          The scalar KK modes $ \phi^{(m)} $ and $ \varphi^{(k)} $ mix with each other through the term $ F_{yz}F^{yz} $ in the 6D bulk action. At first sight, the presence of these scalar mass and mixing terms make gauge invariance difficult to maintain. However, as shown below, the effective action remains gauge invariant because of nontrivial relations among the vector-scalar and scalar-scalar couplings.

          Although there are two types of scalar KK modes, they couple to the vector KK modes in the same manner. Therefore, an analogous relation holds between $ I_1^{(nn)} $ ($ I_2^{(nm)} $) and $ \widetilde{I}_1^{(nm)} $ ($ \widetilde{I}_2^{(nm)} $), as specified in (11). By defining two auxiliary matrices

          $ T_1^{(nm)} = \int {\rm e}^{B_1+B_2}{\rm d}z{\rm d}y \;{\rm e}^{2A-B_1-B_2}\partial_y \rho^{(m)} f^{(n)}, $

          (35)

          $ T_2^{(nk)} = \int {\rm e}^{B_1+B_2}{\rm d}z{\rm d}y \; {\rm e}^{2A-B_1-B_2}\partial_z \chi^{(k)} f^{(n)}, $

          (36)

          with the completeness of the basis functions, we obtain the other relationships.

          $\begin{aligned}[b] C_1^{(nm)}&=\sum\limits_{k,l} T_1^{(kn)}T_1^{(lm)}N^{(kl)}, C_2^{(nm)}\\&=\sum\limits_{k,l} T_2^{(kn)}T_2^{(lm)}N^{(kl)}, \widetilde{C}^{(nm)} \\&= \sum\limits_{k,l} T_1^{(kn)}T_2^{(lm)}N^{(kl)}.\end{aligned} $

          With these identities, the effective action can be rewritten as

          $ \begin{aligned}[b] S^{(4)} =\;& -\frac{1}{4} \sum\limits_{n,m}\int {\rm d}^4 x \sqrt{|g|} \big( \hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(n)}\\& + \widetilde{N}_1^{(nm)} \big[\partial_\nu \phi^{(n)}-\sum\limits_k \widetilde{I}_1^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \phi^{(m)}-\sum\limits_l \widetilde{I}_1^{(lm)} \hat{X}^{\nu(l)}\big]\\ & +\widetilde{N}_2^{(nm)} \big[\partial_\nu \varphi^{(n)}-\sum\limits_k \widetilde{I}_2^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \varphi^{(m)}-\sum\limits_l \widetilde{I}_2^{(lm)} \hat{X}^{\nu(l)}\big]\\ & +\sum\limits_{k,l} N^{(kl)} \big[T_1^{(kn)}\phi^{(n)}-T_2^{(kn)}\varphi^{(n)}\big] \big[T_1^{(lm)}\phi^{(m)}-T_2^{(lm)}\varphi^{(m)}\big] . \end{aligned}$

          (37)

          We observe that, under the gauge transformations,

          $\begin{aligned}[b]& \hat{X}_\nu^{(k)} \rightarrow \hat{X}_\nu^{(k)} + \partial_\nu \Lambda^{(k)}, \quad \phi^{(n)} \rightarrow \phi^{(n)} + \sum\limits_k\widetilde{I}_1^{(kn)} \Lambda^{(k)}, \\&\varphi^{(n)} \rightarrow \varphi^{(n)} + \sum\limits_l\widetilde{I}_2^{(ln)} \Lambda^{(l)} , \end{aligned}$

          (38)

          where the first three terms remain gauge invariant. However, an additional term of the form $ \sum\nolimits_{k,l} (T_1^{(kl)} \widetilde{I}_1^{(nl)} - T_2^{(kl)} \widetilde{I}_2^{(nl)})\Lambda^{(n)} $ emerges in the final expression. Utilizing the symmetry of mixed partial derivatives $ \partial_y \partial_z f^{(n)} = \partial_z \partial_y f^{(n)} $, we obtain

          $ \sum\limits_{l}\left( T_1^{(kl)}\widetilde{I}_1^{(nl)} - T_2^{(kl)}\widetilde{I}_2^{(nl)} \right)=0. $

          (39)

          Thus, this term then vanishes. Consequently, the effective action is gauge-invariant, which is consistent with the result established in Ref. [43].

          However, when we examine the equations that govern the KK modes under Dirichlet or periodic boundary conditions on the basis functions,

          $\begin{aligned}[b]& -{\rm e}^{-(B_1+B_2)} \big(\partial_z ({\rm e}^{2A-B_1+B_2}\partial_z f^{(m)}) + \partial_y( {\rm e}^{2A+B_1-B_2}\partial_y f^{(m)})\big) \\ & = m_v^{2} f^{(m)},\end{aligned} $

          (40)

          $ -{\rm e}^{-(2A-B_1+B_2)} \partial_y ({\rm e}^{4A-B_1-B_2} \partial_y \rho^{(m)} ) = m_\phi^2 \rho^{(m)}, $

          (41)

          $ -{\rm e}^{-(2A+B_1-B_2)} \partial_z ({\rm e}^{4A-B_1-B_2} \partial_z \chi^{(k)} ) = m_\varphi^2 \chi^{(k)}, $

          (42)

          we note that the vector fields acquire masses from the full extra spacetime, while the scalar fields receive mass contributions from only one of the extra dimensions, either y or z. Both the warp factor and wave functions depend on y and z; therefore, the equations for the scalar modes cannot be solved independently. As in the 5D case, it is necessary to modify the bulk action by including the term $ (\nabla_M X^M)^2 $.

        • A.   Review the gauge invariance of massive vector KK modes

        • For a free bulk U(1) gauge field in the 6D brane world, after performing the KK decomposition

          $\begin{aligned}[b] X_\mu &= \sum\limits_n \hat{X}^{(n)}_\mu(x^\nu) f^{(n)}(z,y), X_z \\&= \sum\limits_n \phi^{(n)}(x^\nu) \rho^{(n)}(z,y), X_y = \sum\limits_n \varphi^{(n)}(x^\nu) \chi^{(n)}(z,y),\end{aligned} $

          (27)

          where $\hat{X}^{(n)}_\mu$ represents the n-level vector KK mode and $\phi^{(m)}, \varphi^{(m)}$ represent the two types of scalar ones. The effective action takes the form

          $ \begin{aligned}[b] S_{\text{eff1}} =\;& -\frac{1}{4} \sum\limits_{n,m}\int {\rm d}^4 x \sqrt{|g|} \big( N^{(nm)}\;\hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\ &+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I_1^{(nm)} -2\hat{X}_\nu^{(n)}\partial^\nu \phi^{(m)}\widetilde{I}_1^{(nm)} \\ & + \widetilde{N}_1^{(nm)}\;\partial_\nu \phi^{(n)}\partial^\nu \phi^{(m)}\big]+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I_2^{(nm)} \\ & -2\hat{X}_\nu^{(n)}\partial^\nu \varphi^{(m)}\widetilde{I}_2^{(nm)} + \widetilde{N}_2^{(nm)}\;\partial_\nu \varphi^{(n)}\partial^\nu \varphi^{(m)}\big]\\ & +2\big[C_1^{(nm)}\phi^{(n)}\phi^{(m)}+C_2^{(nm)}\varphi^{(n)}\varphi^{(m)}\\ & -2\widetilde{C}^{(nm)}\phi^{(n)}\varphi^{(m)}\big] \big) \end{aligned} $

          (28)

          where

          $ N^{(nm)} = \int {\rm e}^{B_1+B_2} {\rm d}z {\rm d}y\; f^{(n)} f^{(m)}, $

          (29)

          $\begin{aligned}[b] & \widetilde{N}_1^{(nm)} = \int {\rm e}^{2A-B_1+B_2} {\rm d}z {\rm d}y\; \rho^{(n)} \rho^{(m)},\\&\widetilde{N}_2^{(nm)} = \int {\rm e}^{2A+B_1-B_2} {\rm d}z {\rm d}y\; \chi^{(n)} \chi^{(m)}, \end{aligned}$

          (30)

          $\begin{aligned}[b]& I_1^{(nm)} = \int {\rm e}^{2A-B_1+B_2} {\rm d}z{\rm d}y\; \partial_z f^{(n)} \partial_z f^{(m)},\\& \widetilde{I}_1^{(nm)} = \int {\rm e}^{2A-B_1+B_2} {\rm d}z{\rm d}y\; \partial_z f^{(n)} \rho^{(m)},\end{aligned} $

          (31)

          $\begin{aligned}[b]& I_2^{(nm)} = \int {\rm e}^{2A+B_1-B_2} {\rm d}z{\rm d}y\; \partial_y f^{(n)} \partial_y f^{(m)},\\& \widetilde{I}_2^{(nm)} = \int {\rm e}^{2A+B_1-B_2} {\rm d}z{\rm d}y\; \partial_y f^{(n)} \chi^{(m)},\end{aligned} $

          (32)

          $\begin{aligned}[b]& C_1^{(mm)} = \int {\rm e}^{4A-B_1-B_2} {\rm d}z{\rm d}y\; \partial_y \rho^{(m)} \partial_y \rho^{(m)} ,\\& C_2^{(kk)}= \int {\rm e}^{4A-B_1-B_2} {\rm d}z{\rm d}y\; \partial_z \chi^{(k)} \partial_z \chi^{(k)},\end{aligned} $

          (33)

          $ \widetilde{C}^{(mk)} = \int {\rm e}^{4A-B_1-B_2} {\rm d}z{\rm d}y\; \partial_y \rho^{(m)} \partial_z \chi^{(k)}. $

          (34)

          The scalar KK modes $ \phi^{(m)} $ and $ \varphi^{(k)} $ mix with each other through the term $ F_{yz}F^{yz} $ in the 6D bulk action. At first sight, the presence of these scalar mass and mixing terms make gauge invariance difficult to maintain. However, as shown below, the effective action remains gauge invariant because of nontrivial relations among the vector-scalar and scalar-scalar couplings.

          Although there are two types of scalar KK modes, they couple to the vector KK modes in the same manner. Therefore, an analogous relation holds between $ I_1^{(nn)} $ ($ I_2^{(nm)} $) and $ \widetilde{I}_1^{(nm)} $ ($ \widetilde{I}_2^{(nm)} $), as specified in (11). By defining two auxiliary matrices

          $ T_1^{(nm)} = \int {\rm e}^{B_1+B_2}{\rm d}z{\rm d}y \;{\rm e}^{2A-B_1-B_2}\partial_y \rho^{(m)} f^{(n)}, $

          (35)

          $ T_2^{(nk)} = \int {\rm e}^{B_1+B_2}{\rm d}z{\rm d}y \; {\rm e}^{2A-B_1-B_2}\partial_z \chi^{(k)} f^{(n)}, $

          (36)

          with the completeness of the basis functions, we obtain the other relationships.

          $\begin{aligned}[b] C_1^{(nm)}&=\sum\limits_{k,l} T_1^{(kn)}T_1^{(lm)}N^{(kl)}, C_2^{(nm)}\\&=\sum\limits_{k,l} T_2^{(kn)}T_2^{(lm)}N^{(kl)}, \widetilde{C}^{(nm)} \\&= \sum\limits_{k,l} T_1^{(kn)}T_2^{(lm)}N^{(kl)}.\end{aligned} $

          With these identities, the effective action can be rewritten as

          $ \begin{aligned}[b] S^{(4)} =\;& -\frac{1}{4} \sum\limits_{n,m}\int {\rm d}^4 x \sqrt{|g|} \big( \hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(n)}\\& + \widetilde{N}_1^{(nm)} \big[\partial_\nu \phi^{(n)}-\sum\limits_k \widetilde{I}_1^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \phi^{(m)}-\sum\limits_l \widetilde{I}_1^{(lm)} \hat{X}^{\nu(l)}\big]\\ & +\widetilde{N}_2^{(nm)} \big[\partial_\nu \varphi^{(n)}-\sum\limits_k \widetilde{I}_2^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \varphi^{(m)}-\sum\limits_l \widetilde{I}_2^{(lm)} \hat{X}^{\nu(l)}\big]\\ & +\sum\limits_{k,l} N^{(kl)} \big[T_1^{(kn)}\phi^{(n)}-T_2^{(kn)}\varphi^{(n)}\big] \big[T_1^{(lm)}\phi^{(m)}-T_2^{(lm)}\varphi^{(m)}\big] . \end{aligned}$

          (37)

          We observe that, under the gauge transformations,

          $\begin{aligned}[b]& \hat{X}_\nu^{(k)} \rightarrow \hat{X}_\nu^{(k)} + \partial_\nu \Lambda^{(k)}, \quad \phi^{(n)} \rightarrow \phi^{(n)} + \sum\limits_k\widetilde{I}_1^{(kn)} \Lambda^{(k)}, \\&\varphi^{(n)} \rightarrow \varphi^{(n)} + \sum\limits_l\widetilde{I}_2^{(ln)} \Lambda^{(l)} , \end{aligned}$

          (38)

          where the first three terms remain gauge invariant. However, an additional term of the form $ \sum\nolimits_{k,l} (T_1^{(kl)} \widetilde{I}_1^{(nl)} - T_2^{(kl)} \widetilde{I}_2^{(nl)})\Lambda^{(n)} $ emerges in the final expression. Utilizing the symmetry of mixed partial derivatives $ \partial_y \partial_z f^{(n)} = \partial_z \partial_y f^{(n)} $, we obtain

          $ \sum\limits_{l}\left( T_1^{(kl)}\widetilde{I}_1^{(nl)} - T_2^{(kl)}\widetilde{I}_2^{(nl)} \right)=0. $

          (39)

          Thus, this term then vanishes. Consequently, the effective action is gauge-invariant, which is consistent with the result established in Ref. [43].

          However, when we examine the equations that govern the KK modes under Dirichlet or periodic boundary conditions on the basis functions,

          $\begin{aligned}[b]& -{\rm e}^{-(B_1+B_2)} \big(\partial_z ({\rm e}^{2A-B_1+B_2}\partial_z f^{(m)}) + \partial_y( {\rm e}^{2A+B_1-B_2}\partial_y f^{(m)})\big) \\ & = m_v^{2} f^{(m)},\end{aligned} $

          (40)

          $ -{\rm e}^{-(2A-B_1+B_2)} \partial_y ({\rm e}^{4A-B_1-B_2} \partial_y \rho^{(m)} ) = m_\phi^2 \rho^{(m)}, $

          (41)

          $ -{\rm e}^{-(2A+B_1-B_2)} \partial_z ({\rm e}^{4A-B_1-B_2} \partial_z \chi^{(k)} ) = m_\varphi^2 \chi^{(k)}, $

          (42)

          we note that the vector fields acquire masses from the full extra spacetime, while the scalar fields receive mass contributions from only one of the extra dimensions, either y or z. Both the warp factor and wave functions depend on y and z; therefore, the equations for the scalar modes cannot be solved independently. As in the 5D case, it is necessary to modify the bulk action by including the term $ (\nabla_M X^M)^2 $.

        • A.   Review the gauge invariance of massive vector KK modes

        • For a free bulk U(1) gauge field in the 6D brane world, after performing the KK decomposition

          $\begin{aligned}[b] X_\mu &= \sum\limits_n \hat{X}^{(n)}_\mu(x^\nu) f^{(n)}(z,y), X_z \\&= \sum\limits_n \phi^{(n)}(x^\nu) \rho^{(n)}(z,y), X_y = \sum\limits_n \varphi^{(n)}(x^\nu) \chi^{(n)}(z,y),\end{aligned} $

          (27)

          where $\hat{X}^{(n)}_\mu$ denotes n-level vector KK mode, and $\phi^{(m)}, \varphi^{(m)}$ are two types of scalar ones, the effective action takes the form

          $ \begin{aligned}[b] S_{\text{eff1}} =\;& -\frac{1}{4} \sum\limits_{n,m}\int d^4 x \sqrt{|g|} \big( N^{(nm)}\;\hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(m)}\\ &+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I_1^{(nm)} -2\hat{X}_\nu^{(n)}\partial^\nu \phi^{(m)}\widetilde{I}_1^{(nm)} \\ & + \widetilde{N}_1^{(nm)}\;\partial_\nu \phi^{(n)}\partial^\nu \phi^{(m)}\big]+ 2\big[ \hat{X}^{(n)}_\nu \hat{X}^{\nu(m)} I_2^{(nm)} \\ & -2\hat{X}_\nu^{(n)}\partial^\nu \varphi^{(m)}\widetilde{I}_2^{(nm)} + \widetilde{N}_2^{(nm)}\;\partial_\nu \varphi^{(n)}\partial^\nu \varphi^{(m)}\big]\\ & +2\big[C_1^{(nm)}\phi^{(n)}\phi^{(m)}+C_2^{(nm)}\varphi^{(n)}\varphi^{(m)}\\ & -2\widetilde{C}^{(nm)}\phi^{(n)}\varphi^{(m)}\big] \big) \end{aligned} $

          (28)

          where

          $ N^{(nm)} = \int e^{B_1+B_2} dzdy\; f^{(n)} f^{(m)}, $

          (29)

          $\begin{aligned}[b] & \widetilde{N}_1^{(nm)} = \int e^{2A-B_1+B_2} dzdy\; \rho^{(n)} \rho^{(m)},\\&\widetilde{N}_2^{(nm)} = \int e^{2A+B_1-B_2} dzdy\; \chi^{(n)} \chi^{(m)}, \end{aligned}$

          (30)

          $\begin{aligned}[b]& I_1^{(nm)} = \int e^{2A-B_1+B_2} dzdy\; \partial_z f^{(n)} \partial_z f^{(m)},\\& \widetilde{I}_1^{(nm)} = \int e^{2A-B_1+B_2} dzdy\; \partial_z f^{(n)} \rho^{(m)},\end{aligned} $

          (31)

          $\begin{aligned}[b]& I_2^{(nm)} = \int e^{2A+B_1-B_2} dzdy\; \partial_y f^{(n)} \partial_y f^{(m)},\\& \widetilde{I}_2^{(nm)} = \int e^{2A+B_1-B_2} dzdy\; \partial_y f^{(n)} \chi^{(m)},\end{aligned} $

          (32)

          $\begin{aligned}[b]& C_1^{(mm)} = \int e^{4A-B_1-B_2} dzdy\; \partial_y \rho^{(m)} \partial_y \rho^{(m)} ,\\& C_2^{(kk)}= \int e^{4A-B_1-B_2} dzdy\; \partial_z \chi^{(k)} \partial_z \chi^{(k)},\end{aligned} $

          (33)

          $ \widetilde{C}^{(mk)} = \int e^{4A-B_1-B_2} dzdy\; \partial_y \rho^{(m)} \partial_z \chi^{(k)}. $

          (34)

          We note that the scalar KK modes $ \phi^{(m)} $ and $ \varphi^{(k)} $ mix with each other through the term $ F_{yz}F^{yz} $ in the 6D bulk action. At the first sight the presence of these scalar mass and mixing terms appears to make gauge invariance difficult to maintain. However, as we show below, the effective action remains gauge invariant due to nontrivial relations among the vector-scalar and scalar-scalar couplings.

          Although there are two types of scalar KK modes, they couple to the vector KK modes in the same manner. Therefore, an analogous relation holds between $ I_1^{(nn)} $ ($ I_2^{(nm)} $) and $ \widetilde{I}_1^{(nm)} $ ($ \widetilde{I}_2^{(nm)} $) as specified in (11). While by defining two auxiliary matrices

          $ T_1^{(nm)} = \int e^{B_1+B_2}dzdy \;e^{2A-B_1-B_2}\partial_y \rho^{(m)} f^{(n)}, $

          (35)

          $ T_2^{(nk)} = \int e^{B_1+B_2}dzdy \; e^{2A-B_1-B_2}\partial_z \chi^{(k)} f^{(n)}, $

          (36)

          together with the completeness of the basis functions, we obtain other relations:

          $\begin{aligned}[b] C_1^{(nm)}&=\sum\limits_{k,l} T_1^{(kn)}T_1^{(lm)}N^{(kl)}, C_2^{(nm)}\\&=\sum\limits_{k,l} T_2^{(kn)}T_2^{(lm)}N^{(kl)}, \widetilde{C}^{(nm)} \\&= \sum\limits_{k,l} T_1^{(kn)}T_2^{(lm)}N^{(kl)}.\end{aligned} $

          With these identities, the effective action can be rewritten as

          $ \begin{aligned}[b] S^{(4)} =\;& -\frac{1}{4} \sum\limits_{n,m}\int d^4 x \sqrt{|g|} \big( \hat{F}_{\mu\nu}^{(n)} \hat{F}^{\mu\nu(n)}\\& + \widetilde{N}_1^{(nm)} \big[\partial_\nu \phi^{(n)}-\sum\limits_k \widetilde{I}_1^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \phi^{(m)}-\sum\limits_l \widetilde{I}_1^{(lm)} \hat{X}^{\nu(l)}\big]\\ & +\widetilde{N}_2^{(nm)} \big[\partial_\nu \varphi^{(n)}-\sum\limits_k \widetilde{I}_2^{(kn)} \hat{X}^{(k)}_\nu\big] \big[\partial^\nu \varphi^{(m)}-\sum\limits_l \widetilde{I}_2^{(lm)} \hat{X}^{\nu(l)}\big]\\ & +\sum\limits_{k,l} N^{(kl)} \big[T_1^{(kn)}\phi^{(n)}-T_2^{(kn)}\varphi^{(n)}\big] \big[T_1^{(lm)}\phi^{(m)}-T_2^{(lm)}\varphi^{(m)}\big] . \end{aligned}$

          (37)

          We observe that under the gauge transformations

          $\begin{aligned}[b]& \hat{X}_\nu^{(k)} \rightarrow \hat{X}_\nu^{(k)} + \partial_\nu \Lambda^{(k)}, \phi^{(n)} \rightarrow \phi^{(n)} + \sum\limits_k\widetilde{I}_1^{(kn)} \Lambda^{(k)}, \\&\varphi^{(n)} \rightarrow \varphi^{(n)} + \sum\limits_l\widetilde{I}_2^{(ln)} \Lambda^{(l)} \end{aligned}$

          (38)

          the first three terms remain gauge invariant. But an additional term of the form $ \sum\nolimits_{k,l} (T_1^{(kl)} \widetilde{I}_1^{(nl)} - T_2^{(kl)} \widetilde{I}_2^{(nl)})\Lambda^{(n)} $ emerges in the final expression. Utilizing the symmetry of mixed partial derivatives $ \partial_y \partial_z f^{(n)} = \partial_z \partial_y f^{(n)} $, we obtain

          $ \sum\limits_{l}\left( T_1^{(kl)}\widetilde{I}_1^{(nl)} - T_2^{(kl)}\widetilde{I}_2^{(nl)} \right)=0. $

          (39)

          This term thus vanishes. Consequently, the effective action is gauge-invariant, consistent with the result established in Ref.[43].

          However, when we examine the equations that govern the KK modes under Dirichlet or periodic boundary conditions on the basis functions

          $\begin{aligned}[b]& -e^{-(B_1+B_2)} \big(\partial_z (e^{2A-B_1+B_2}\partial_z f^{(m)})+ \partial_y( e^{2A+B_1-B_2}\partial_y f^{(m)})\big) \\=\;& m_v^{2} f^{(m)},\end{aligned} $

          (40)

          $ -e^{-(2A-B_1+B_2)} \partial_y (e^{4A-B_1-B_2} \partial_y \rho^{(m)} ) = m_\phi^2 \rho^{(m)}, $

          (41)

          $ -e^{-(2A+B_1-B_2)} \partial_z (e^{4A-B_1-B_2} \partial_z \chi^{(k)} ) = m_\varphi^2 \chi^{(k)}, $

          (42)

          we note that the vector fields acquire masses from the full extra spacetime, whereas the scalar fields receive mass contributions from only one of the extra dimensions, either y or z. Because both the warp factor and the wave functions depend on y and z, the equations for the scalar modes cannot be solved independently. Therefore, as in the five-dimensional case, it is necessary to modify the bulk action by including the term $ (\nabla_M X^M)^2 $.

        • B.   New bulk action for the U(1) gauge field and their KK modes

        • By the KK decomposition, the additional term $ (\nabla_M X^M)^2 $ becomes

          $ \begin{aligned}[b]& \int dzdy\sqrt{-g}(\nabla_M X^M)^2=\sum\limits_{n,m}\sqrt{-\hat{g}}\big( N^{(nm)}(\partial^\mu \hat{X}_\mu^{(n)}) (\partial^\mu \hat{X}_\mu^{(m)})\\&\quad +C_1^{'(nm)}\phi^{(n)} \phi^{(m)} +C_2^{'(nm)} \varphi^{(n)} \varphi^{(m)}\\&\quad +2\widetilde{I}_1^{'(nm)}\partial^\mu \hat{X}^{(n)}_\mu \phi^{(m)} +2\widetilde{I}_2^{'(nm}\partial^\mu \hat{X}^{(n)}_\mu \varphi^{(m)} +2\widetilde{I}^{'(nm)} \phi^{(n)} \varphi^{(m)} \big), \end{aligned} $

          (43)

          where

          $ \begin{aligned}[b]C_1^{'(nm)} =\;& \int {e}^{-(4 A+B_1+B_2)} dzdy \;(\partial_z({e}^{4 A-B_1+B_2} \rho^{(n)}))\\&\times(\partial_z({e}^{4 A-B_1+B_2} \rho^{(m)})), \end{aligned}$

          (44)

          $\begin{aligned}[b] C_2^{'(nm)} =\;& \int {e}^{-(4 A+B_1+B_2)} dzdy \;(\partial_y({e}^{4 A+B_1-B_2}\chi^{(n)}))\\&\times(\partial_y({e}^{4 A+B_1-B_2}\chi^{(m)})),\end{aligned} $

          (45)

          $ \widetilde{I}_1^{'(nm)} = \int e^{-2A}dzdy \;f^{(n)}(\partial_z({e}^{4 A-B_1+B_2} \rho^{(m)})), $

          (46)

          $ \widetilde{I}_2^{'(nm)} = \int e^{-2A}dzdy \;f^{(n)}(\partial_y({e}^{4 A+B_1-B_2}\chi^{(m)})), $

          (47)

          $\begin{aligned}[b] \widetilde{C}^{'(nm)} =\;& \int {e}^{-(4 A+B_1+B_2)} dzdy \;(\partial_z({e}^{4 A-B_1+B_2} \rho^{(n)}))\\&\times(\partial_y({e}^{4 A+B_1-B_2}\chi^{(m)})).\end{aligned} $

          (48)

          Considering

          $ {e}^{-2A-B_1-B_2}(\partial_z({e}^{4 A-B_1+B_2} \rho^{(m)})) = \sum\limits_n \widetilde{I}_1^{'(nm)} f^{(n)} , $

          (49)

          $ {e}^{-2A-B_1-B_2}(\partial_y({e}^{4 A+B_1-B_2} \chi^{(k)})) = \sum\limits_n \widetilde{I}_2^{'(nk)} f^{(n)} , $

          (50)

          we get

          $\begin{aligned}[b]& C_1^{'(nm)}=\sum\limits_{k,l}\widetilde{I}_1^{'(kn)}\widetilde{I}_1^{'(lm)}N^{(kl)},\\& C_2^{'(nm)}=\sum\limits_{k,l}\widetilde{I}_2^{'(kn)}\widetilde{I}_2^{'(lm)}N^{(kl)},\\& \widetilde{I}^{'(nm)}= \sum\limits_{k,l}\widetilde{I}_1^{'(kn)}\widetilde{I}_2^{'(lm)}N^{(kl)}.\end{aligned} $

          This makes the effective term (43) for $ (\nabla_M X^M)^2 $ at last turn to

          $ \sum\limits_{n}\sqrt{-\hat{g}}\big[\partial^\mu \hat{X}_\mu^{(n)} +\sum\limits_k \widetilde{I}_1^{'(nk)}\;\phi^{(k)} +\sum\limits_k \widetilde{I}_2^{'(nk)}\;\varphi^{(k)}\big]^2, $

          (51)

          which is gauge invariant under the gauge transformations (38) with the condition

          $ \partial^\mu\partial_\mu \Lambda^{(n)} +\sum\limits_{k,l}(\widetilde{I}_1^{'(nk)}\widetilde{I}_1^{(lk)}+\widetilde{I}_2^{'(nk)}\widetilde{I}_2^{(lk)})\Lambda^{(l)}=0. $

          (52)

          Now for the scalars, the masses of them are $ C_1^{(nm)}+C_1^{'(nm)}=m_{\phi}^{2}\delta^{nm} $ and $ C_2^{(nm)}+C_2^{'(nm)}=m_{\varphi}^{2}\delta^{nm} $, then there are

          $\begin{aligned}[b]& -e^{-(2A-B_1+B_2)} \partial_y (e^{4A-B_1-B_2} \partial_y \rho^{(m)} ) \\&-{e}^{2 A}\partial_z\big( {e}^{-(4 A+B_1+B_2)}\partial_z({e}^{4 A-B_1+B_2} \rho^{(n)})\big)=m_\phi^2 \rho^{(m)}, \end{aligned} $

          (53)

          $\begin{aligned}[b]& -e^{-(2A+B_1-B_2)} \partial_z (e^{4A-B_1-B_2} \partial_z \chi^{(k)} ) \\&-{e}^{2 A}\partial_y\big( {e}^{-(4 A+B_1+B_2)}\partial_y({e}^{4 A+B_1-B_2}\chi^{(k)})\big) =m_\varphi^2 \chi^{(k)},\end{aligned} $

          (54)

          The eigenfunctions for these three equations (17), (53) and (54) are chosen as the basis functions $ f^{(n)} $, $ \rho^{(n)} $ and $ \chi^{(n)} $. It is noted that these two types of scalars interact with each other and can exhibit mixing.

          We note that in most 6D brane models, the warp factors are assumed to depend only on a single extra dimension y. By separating the three basis functions as $ f^n = R_1^{(n,l)}(r) \Theta^{(l)}(\theta) $, $ \rho^m = R_2^{(m,l)}(r) \Theta^{(l)}(\theta) $, and $ \chi^k = R_3^{(k,l)}(r) \Theta^{(l)}(\theta) $, assuming $ \partial_\theta^2 \Theta^{(l)} + l^2 \Theta^{(l)} = 0 $, and applying the transformations $ d\bar{r}=e^{-A+B_1}dr $, $ \bar{R}_1^{(n,l)}=e^{\frac{1}{2}(A+B_2)}R_1^{(n,l)} $, $ \bar{R}_2^{(m,l)}=e^{\frac{1}{2}(3A-2B_1+B_2)}R_2^{(m,l)} $, $ \bar{R}_3^{(k,l)}=e^{\frac{1}{2}(3A-B_2)}R_3^{(k,l)} $, we find that the three equations (17), (53), and (54) reduce to Schrödinger-like equations:

          $ -\partial_{\bar{r}}^2 \bar{R}_1^{(n,l)} + V_{\text{eff1}} \bar{R}_1^{(n,l)} = m_v^{2} \bar{R}_1^{(n,l)}, $

          (55)

          $ -\partial_{\bar{r}}^2 \bar{R}_2^{(m,l)} + V_{\text{eff2}}\bar{R}_2^{(m,l)} = m_\phi^2 \bar{R}_2^{(m,l)}, $

          (56)

          $ -\partial_{\bar{r}}^2 \bar{R}_3^{(k,l)} + V_{\text{eff3}}\bar{R}_3^{(k,l)} =m_\varphi^2 \bar{R}_3^{(k,l)}. $

          (57)

          Here, the effective potentials are given by

          $ V_{\text{eff1}} = \frac{1}{4}(\partial{\bar{r}} (A+B_2))^2+\frac{1}{2}\partial{\bar{r}}^2(A+B_2)+ e^{2A-2B_2}l^2, $

          (58)

          $\begin{aligned}[b] V_{\text{eff2}} =\;& \frac{1}{4}(\partial_{\bar{r}} (3A-2B_1+B_2))^2+\frac{1}{2}\partial_{\bar{r}}^2(3A-2B_1+B_2)\\&+ e^{2A-2B_2}l^2 -\big[\partial_{\bar{r}}^2( 4A-B_1+B_2)\\&+\partial_{\bar{r}}(-A-B_1)\partial_{\bar{r}}( 4A-B_1+B_2)\big],\end{aligned} $

          (59)

          $ V_{\text{eff3}} = \frac{1}{4}(\partial_{\bar{r}} (3A-B_2))^2+\frac{1}{2}\partial_{\bar{r}}^2(3A-B_2)+ e^{2A-2B_2}l^2. $

          (60)

          These effective potentials can then be analyzed to study the mass spectra of the Kaluza–Klein modes in specific brane models.

        • B.   New bulk action for the U(1) gauge field and their KK modes

        • By the KK decomposition, the additional term $ (\nabla_M X^M)^2 $ becomes

          $ \begin{aligned}[b]& \int {\rm d}z{\rm d}y\sqrt{-g}(\nabla_M X^M)^2=\sum\limits_{n,m}\sqrt{-\hat{g}}\big( N^{(nm)}(\partial^\mu \hat{X}_\mu^{(n)}) (\partial^\mu \hat{X}_\mu^{(m)})\\&\quad +C_1^{'(nm)}\phi^{(n)} \phi^{(m)} +C_2^{'(nm)} \varphi^{(n)} \varphi^{(m)}\\&\quad +2\widetilde{I}_1^{'(nm)}\partial^\mu \hat{X}^{(n)}_\mu \phi^{(m)} +2\widetilde{I}_2^{'(nm}\partial^\mu \hat{X}^{(n)}_\mu \varphi^{(m)} +2\widetilde{I}^{'(nm)} \phi^{(n)} \varphi^{(m)} \big), \end{aligned} $

          (43)

          where

          $ \begin{aligned}[b]C_1^{'(nm)} =\;& \int {\rm e}^{-(4 A+B_1+B_2)} {\rm d}z{\rm d}y \;(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(n)}))\\&\times(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(m)})), \end{aligned}$

          (44)

          $\begin{aligned}[b] C_2^{'(nm)} =\;& \int {\rm e}^{-(4 A+B_1+B_2)} {\rm d}z{\rm d}y \;(\partial_y({\rm e}^{4 A+B_1-B_2}\chi^{(n)}))\\&\times(\partial_y({\rm e}^{4 A+B_1-B_2}\chi^{(m)})),\end{aligned} $

          (45)

          $ \widetilde{I}_1^{'(nm)} = \int {\rm e}^{-2A}{\rm d}z{\rm d}y \;f^{(n)}(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(m)})), $

          (46)

          $ \widetilde{I}_2^{'(nm)} = \int {\rm e}^{-2A}{\rm d}z{\rm d}y \;f^{(n)}(\partial_y({e}^{4 A+B_1-B_2}\chi^{(m)})), $

          (47)

          $\begin{aligned}[b] \widetilde{C}^{'(nm)} =\;& \int {\rm e}^{-(4 A+B_1+B_2)} {\rm d}z{\rm d}y \;(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(n)}))\\&\times(\partial_y({\rm e}^{4 A+B_1-B_2}\chi^{(m)})).\end{aligned} $

          (48)

          Considering

          $ {\rm e}^{-2A-B_1-B_2}(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(m)})) = \sum\limits_n \widetilde{I}_1^{'(nm)} f^{(n)} , $

          (49)

          $ {\rm e}^{-2A-B_1-B_2}(\partial_y({\rm e}^{4 A+B_1-B_2} \chi^{(k)})) = \sum\limits_n \widetilde{I}_2^{'(nk)} f^{(n)} , $

          (50)

          we get

          $\begin{aligned}[b]& C_1^{'(nm)}=\sum\limits_{k,l}\widetilde{I}_1^{'(kn)}\widetilde{I}_1^{'(lm)}N^{(kl)},\\& C_2^{'(nm)}=\sum\limits_{k,l}\widetilde{I}_2^{'(kn)}\widetilde{I}_2^{'(lm)}N^{(kl)},\\& \widetilde{I}^{'(nm)}= \sum\limits_{k,l}\widetilde{I}_1^{'(kn)}\widetilde{I}_2^{'(lm)}N^{(kl)}.\end{aligned} $

          This makes the effective term (43) for $ (\nabla_M X^M)^2 $ at last turn to

          $ \sum\limits_{n}\sqrt{-\hat{g}}\big[\partial^\mu \hat{X}_\mu^{(n)} +\sum\limits_k \widetilde{I}_1^{'(nk)}\;\phi^{(k)} +\sum\limits_k \widetilde{I}_2^{'(nk)}\;\varphi^{(k)}\big]^2, $

          (51)

          which is gauge invariant under gauge transformations (38) with the condition

          $ \partial^\mu\partial_\mu \Lambda^{(n)} +\sum\limits_{k,l}(\widetilde{I}_1^{'(nk)}\widetilde{I}_1^{(lk)}+\widetilde{I}_2^{'(nk)}\widetilde{I}_2^{(lk)})\Lambda^{(l)}=0. $

          (52)

          For the scalars, the masses are $ C_1^{(nm)}+C_1^{'(nm)}=m_{\phi}^{2}\delta^{nm} $ and $ C_2^{(nm)}+C_2^{'(nm)}=m_{\varphi}^{2}\delta^{nm} $. Then, there are

          $\begin{aligned}[b]& -{\rm e}^{-(2A-B_1+B_2)} \partial_y ({\rm e}^{4A-B_1-B_2} \partial_y \rho^{(m)} ) \\&-{\rm e}^{2 A}\partial_z\big( {\rm e}^{-(4 A+B_1+B_2)}\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(n)})\big)=m_\phi^2 \rho^{(m)}, \end{aligned} $

          (53)

          $\begin{aligned}[b]& -{\rm e}^{-(2A+B_1-B_2)} \partial_z ({\rm e}^{4A-B_1-B_2} \partial_z \chi^{(k)} ) \\&-{\rm e}^{2 A}\partial_y\big( {\rm e}^{-(4 A+B_1+B_2)}\partial_y({\rm e}^{4 A+B_1-B_2}\chi^{(k)})\big) =m_\varphi^2 \chi^{(k)},\end{aligned} $

          (54)

          The eigenfunctions for (17), (53), and (54) are selected as the basis functions $ f^{(n)} $, $ \rho^{(n)} $, and $ \chi^{(n)} $. These two types of scalars interact with each other and can exhibit mixing.

          In most 6D brane models, warp factors are assumed to depend only on a single extra dimension y. By separating three basis functions as $ f^n = R_1^{(n,l)}(r) \Theta^{(l)}(\theta) $, $ \rho^m = R_2^{(m,l)}(r) \Theta^{(l)}(\theta) $, and $ \chi^k = R_3^{(k,l)}(r) \Theta^{(l)}(\theta) $, assuming $ \partial_\theta^2 \Theta^{(l)} + l^2 \Theta^{(l)} = 0 $ and applying the transformations ${\rm d}\bar{r}={\rm e}^{-A+B_1}{\rm d}r$, $\bar{R}_1^{(n,l)}={\rm e}^{\frac{1}{2}(A+B_2)}R_1^{(n,l)}$, $\bar{R}_2^{(m,l)}={\rm e}^{\frac{1}{2}(3A-2B_1+B_2)}R_2^{(m,l)}$, and $\bar{R}_3^{(k,l)}={\rm e}^{\frac{1}{2}(3A-B_2)}R_3^{(k,l)}$, we find that (17), (53), and (54) reduce to Schrödinger-like equations:

          $ -\partial_{\bar{r}}^2 \bar{R}_1^{(n,l)} + V_{\text{eff1}} \bar{R}_1^{(n,l)} = m_v^{2} \bar{R}_1^{(n,l)}, $

          (55)

          $ -\partial_{\bar{r}}^2 \bar{R}_2^{(m,l)} + V_{\text{eff2}}\bar{R}_2^{(m,l)} = m_\phi^2 \bar{R}_2^{(m,l)}, $

          (56)

          $ -\partial_{\bar{r}}^2 \bar{R}_3^{(k,l)} + V_{\text{eff3}}\bar{R}_3^{(k,l)} =m_\varphi^2 \bar{R}_3^{(k,l)}. $

          (57)

          Here, the effective potentials are given by

          $ V_{\text{eff1}} = \frac{1}{4}(\partial_{\bar{r}} (A+B_2))^2+\frac{1}{2}\partial_{\bar{r}}^2(A+B_2)+ {\rm e}^{2A-2B_2}l^2, $

          (58)

          $\begin{aligned}[b] V_{\text{eff2}} =\;& \frac{1}{4}(\partial_{\bar{r}} (3A-2B_1+B_2))^2+\frac{1}{2}\partial_{\bar{r}}^2(3A-2B_1+B_2)\\&+ {\rm e}^{2A-2B_2}l^2 -\big[\partial_{\bar{r}}^2( 4A-B_1+B_2)\\&+\partial_{\bar{r}}(-A-B_1)\partial_{\bar{r}}( 4A-B_1+B_2)\big],\end{aligned} $

          (59)

          $ V_{\text{eff3}} = \frac{1}{4}(\partial_{\bar{r}} (3A-B_2))^2+\frac{1}{2}\partial_{\bar{r}}^2(3A-B_2)+{\rm e}^{2A-2B_2}l^2. $

          (60)

          These effective potentials can be analyzed to study the mass spectra of the KK modes in specific brane models.

        • B.   New bulk action for the U(1) gauge field and their KK modes

        • By the KK decomposition, the additional term $ (\nabla_M X^M)^2 $ becomes

          $ \begin{aligned}[b]& \int {\rm d}z{\rm d}y\sqrt{-g}(\nabla_M X^M)^2=\sum\limits_{n,m}\sqrt{-\hat{g}}\big( N^{(nm)}(\partial^\mu \hat{X}_\mu^{(n)}) (\partial^\mu \hat{X}_\mu^{(m)})\\&\quad +C_1^{'(nm)}\phi^{(n)} \phi^{(m)} +C_2^{'(nm)} \varphi^{(n)} \varphi^{(m)}\\&\quad +2\widetilde{I}_1^{'(nm)}\partial^\mu \hat{X}^{(n)}_\mu \phi^{(m)} +2\widetilde{I}_2^{'(nm}\partial^\mu \hat{X}^{(n)}_\mu \varphi^{(m)} +2\widetilde{I}^{'(nm)} \phi^{(n)} \varphi^{(m)} \big), \end{aligned} $

          (43)

          where

          $ \begin{aligned}[b]C_1^{'(nm)} =\;& \int {\rm e}^{-(4 A+B_1+B_2)} {\rm d}z{\rm d}y \;(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(n)}))\\&\times(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(m)})), \end{aligned}$

          (44)

          $\begin{aligned}[b] C_2^{'(nm)} =\;& \int {\rm e}^{-(4 A+B_1+B_2)} {\rm d}z{\rm d}y \;(\partial_y({\rm e}^{4 A+B_1-B_2}\chi^{(n)}))\\&\times(\partial_y({\rm e}^{4 A+B_1-B_2}\chi^{(m)})),\end{aligned} $

          (45)

          $ \widetilde{I}_1^{'(nm)} = \int {\rm e}^{-2A}{\rm d}z{\rm d}y \;f^{(n)}(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(m)})), $

          (46)

          $ \widetilde{I}_2^{'(nm)} = \int {\rm e}^{-2A}{\rm d}z{\rm d}y \;f^{(n)}(\partial_y({e}^{4 A+B_1-B_2}\chi^{(m)})), $

          (47)

          $\begin{aligned}[b] \widetilde{C}^{'(nm)} =\;& \int {\rm e}^{-(4 A+B_1+B_2)} {\rm d}z{\rm d}y \;(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(n)}))\\&\times(\partial_y({\rm e}^{4 A+B_1-B_2}\chi^{(m)})).\end{aligned} $

          (48)

          Considering

          $ {\rm e}^{-2A-B_1-B_2}(\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(m)})) = \sum\limits_n \widetilde{I}_1^{'(nm)} f^{(n)} , $

          (49)

          $ {\rm e}^{-2A-B_1-B_2}(\partial_y({\rm e}^{4 A+B_1-B_2} \chi^{(k)})) = \sum\limits_n \widetilde{I}_2^{'(nk)} f^{(n)} , $

          (50)

          we get

          $\begin{aligned}[b]& C_1^{'(nm)}=\sum\limits_{k,l}\widetilde{I}_1^{'(kn)}\widetilde{I}_1^{'(lm)}N^{(kl)},\\& C_2^{'(nm)}=\sum\limits_{k,l}\widetilde{I}_2^{'(kn)}\widetilde{I}_2^{'(lm)}N^{(kl)},\\& \widetilde{I}^{'(nm)}= \sum\limits_{k,l}\widetilde{I}_1^{'(kn)}\widetilde{I}_2^{'(lm)}N^{(kl)}.\end{aligned} $

          This makes the effective term (43) for $ (\nabla_M X^M)^2 $ at last turn to

          $ \sum\limits_{n}\sqrt{-\hat{g}}\big[\partial^\mu \hat{X}_\mu^{(n)} +\sum\limits_k \widetilde{I}_1^{'(nk)}\;\phi^{(k)} +\sum\limits_k \widetilde{I}_2^{'(nk)}\;\varphi^{(k)}\big]^2, $

          (51)

          which is gauge invariant under gauge transformations (38) with the condition

          $ \partial^\mu\partial_\mu \Lambda^{(n)} +\sum\limits_{k,l}(\widetilde{I}_1^{'(nk)}\widetilde{I}_1^{(lk)}+\widetilde{I}_2^{'(nk)}\widetilde{I}_2^{(lk)})\Lambda^{(l)}=0. $

          (52)

          For the scalars, the masses are $ C_1^{(nm)}+C_1^{'(nm)}=m_{\phi}^{2}\delta^{nm} $ and $ C_2^{(nm)}+C_2^{'(nm)}=m_{\varphi}^{2}\delta^{nm} $. Then, there are

          $\begin{aligned}[b]& -{\rm e}^{-(2A-B_1+B_2)} \partial_y ({\rm e}^{4A-B_1-B_2} \partial_y \rho^{(m)} ) \\&-{\rm e}^{2 A}\partial_z\big( {\rm e}^{-(4 A+B_1+B_2)}\partial_z({\rm e}^{4 A-B_1+B_2} \rho^{(n)})\big)=m_\phi^2 \rho^{(m)}, \end{aligned} $

          (53)

          $\begin{aligned}[b]& -{\rm e}^{-(2A+B_1-B_2)} \partial_z ({\rm e}^{4A-B_1-B_2} \partial_z \chi^{(k)} ) \\&-{\rm e}^{2 A}\partial_y\big( {\rm e}^{-(4 A+B_1+B_2)}\partial_y({\rm e}^{4 A+B_1-B_2}\chi^{(k)})\big) =m_\varphi^2 \chi^{(k)},\end{aligned} $

          (54)

          The eigenfunctions for (17), (53), and (54) are selected as the basis functions $ f^{(n)} $, $ \rho^{(n)} $, and $ \chi^{(n)} $. These two types of scalars interact with each other and can exhibit mixing.

          In most 6D brane models, warp factors are assumed to depend only on a single extra dimension y. By separating three basis functions as $ f^n = R_1^{(n,l)}(r) \Theta^{(l)}(\theta) $, $ \rho^m = R_2^{(m,l)}(r) \Theta^{(l)}(\theta) $, and $ \chi^k = R_3^{(k,l)}(r) \Theta^{(l)}(\theta) $, assuming $ \partial_\theta^2 \Theta^{(l)} + l^2 \Theta^{(l)} = 0 $ and applying the transformations ${\rm d}\bar{r}={\rm e}^{-A+B_1}{\rm d}r$, $\bar{R}_1^{(n,l)}={\rm e}^{\frac{1}{2}(A+B_2)}R_1^{(n,l)}$, $\bar{R}_2^{(m,l)}={\rm e}^{\frac{1}{2}(3A-2B_1+B_2)}R_2^{(m,l)}$, and $\bar{R}_3^{(k,l)}={\rm e}^{\frac{1}{2}(3A-B_2)}R_3^{(k,l)}$, we find that (17), (53), and (54) reduce to Schrödinger-like equations:

          $ -\partial_{\bar{r}}^2 \bar{R}_1^{(n,l)} + V_{\text{eff1}} \bar{R}_1^{(n,l)} = m_v^{2} \bar{R}_1^{(n,l)}, $

          (55)

          $ -\partial_{\bar{r}}^2 \bar{R}_2^{(m,l)} + V_{\text{eff2}}\bar{R}_2^{(m,l)} = m_\phi^2 \bar{R}_2^{(m,l)}, $

          (56)

          $ -\partial_{\bar{r}}^2 \bar{R}_3^{(k,l)} + V_{\text{eff3}}\bar{R}_3^{(k,l)} =m_\varphi^2 \bar{R}_3^{(k,l)}. $

          (57)

          Here, the effective potentials are given by

          $ V_{\text{eff1}} = \frac{1}{4}(\partial_{\bar{r}} (A+B_2))^2+\frac{1}{2}\partial_{\bar{r}}^2(A+B_2)+ {\rm e}^{2A-2B_2}l^2, $

          (58)

          $\begin{aligned}[b] V_{\text{eff2}} =\;& \frac{1}{4}(\partial_{\bar{r}} (3A-2B_1+B_2))^2+\frac{1}{2}\partial_{\bar{r}}^2(3A-2B_1+B_2)\\&+ {\rm e}^{2A-2B_2}l^2 -\big[\partial_{\bar{r}}^2( 4A-B_1+B_2)\\&+\partial_{\bar{r}}(-A-B_1)\partial_{\bar{r}}( 4A-B_1+B_2)\big],\end{aligned} $

          (59)

          $ V_{\text{eff3}} = \frac{1}{4}(\partial_{\bar{r}} (3A-B_2))^2+\frac{1}{2}\partial_{\bar{r}}^2(3A-B_2)+{\rm e}^{2A-2B_2}l^2. $

          (60)

          These effective potentials can be analyzed to study the mass spectra of the KK modes in specific brane models.

        IV.   THE MIXING AND OSCILLATIONS OF THE TWO SCALAR KK MODES
        • Note that the two types of scalars couple together; the coupling coefficients $ C^{(mk)} $ read

          $ C^{(mk)}=\widetilde{C}^{(mk)}+\widetilde{C}^{'(mk)}=\int d\bar{r}\; e^{-2B_2}\partial_{\bar{r}}(e^{A+B_2} \bar{R}_2^{(m,l)}\bar{R}_3^{(k,l)}). $

          (61)

          These couplings lead to the oscillation of the corresponding scalar KK modes. Typically, the physical masses $ {\cal{M}}^{2(n)}_A $ of these oscillations are obtained by diagonalizing the mass matrix $ {\cal{M}}^2 $, which is given by:

          $ {\cal{M}}^2=\left(\begin{array}{ll} {\left[m_\phi^2\right]_{n \times n}} & {[C]_{n \times n}} \\ {\left[C^T\right]_{n \times n}} & {\left[m_{\varphi}^2\right]_{n \times n}} \end{array}\right) $

          (62)

          Here, the diagonal blocks correspond to the KK mass-squared of the two types of scalar fields, while the off-diagonal blocks characterize the coupling between them.

          In the following, we will calculate the masses of these oscillations in two different 6D brane models.

        IV.   MIXING AND OSCILLATIONS OF THE TWO SCALAR KK MODES
        • The two types of scalars couple together. The coupling coefficients $ C^{(mk)} $ read

          $ C^{(mk)}=\widetilde{C}^{(mk)}+\widetilde{C}^{'(mk)}=\int {\rm d}\bar{r}\; {\rm e}^{-2B_2}\partial_{\bar{r}}({\rm e}^{A+B_2} \bar{R}_2^{(m,l)}\bar{R}_3^{(k,l)}). $

          (61)

          These couplings lead to the oscillation of the corresponding scalar KK modes. The physical masses $ {\cal{M}}^{2(n)}_A $ of these oscillations are obtained by diagonalizing the mass matrix $ {\cal{M}}^2 $, which is given by

          $ {\cal{M}}^2=\left(\begin{array}{ll} {\left[m_\phi^2\right]_{n \times n}} & {[C]_{n \times n}} \\ {\left[C^T\right]_{n \times n}} & {\left[m_{\varphi}^2\right]_{n \times n}} \end{array}\right). $

          (62)

          The diagonal blocks correspond to the KK mass-squared of the two types of scalar fields, while the off-diagonal blocks characterize the coupling between them.

          In the following, we calculate the masses of these oscillations in two different 6D brane models.

        IV.   MIXING AND OSCILLATIONS OF THE TWO SCALAR KK MODES
        • The two types of scalars couple together. The coupling coefficients $ C^{(mk)} $ read

          $ C^{(mk)}=\widetilde{C}^{(mk)}+\widetilde{C}^{'(mk)}=\int {\rm d}\bar{r}\; {\rm e}^{-2B_2}\partial_{\bar{r}}({\rm e}^{A+B_2} \bar{R}_2^{(m,l)}\bar{R}_3^{(k,l)}). $

          (61)

          These couplings lead to the oscillation of the corresponding scalar KK modes. The physical masses $ {\cal{M}}^{2(n)}_A $ of these oscillations are obtained by diagonalizing the mass matrix $ {\cal{M}}^2 $, which is given by

          $ {\cal{M}}^2=\left(\begin{array}{ll} {\left[m_\phi^2\right]_{n \times n}} & {[C]_{n \times n}} \\ {\left[C^T\right]_{n \times n}} & {\left[m_{\varphi}^2\right]_{n \times n}} \end{array}\right). $

          (62)

          The diagonal blocks correspond to the KK mass-squared of the two types of scalar fields, while the off-diagonal blocks characterize the coupling between them.

          In the following, we calculate the masses of these oscillations in two different 6D brane models.

        • A.   Case 1

        • In Ref. [45], some 6D thick-brane models are considered from a cosmological perspective, where the metric is assumed as

          $ {\rm d} s^2=M^2(r) g_{\mu \nu} {\rm d} x^\mu {\rm d} x^\nu+ {\rm d} r^2+L(r)^2 {\rm d} \theta^2. $

          (63)

          The metric functions $ M(r) $ and $ L(r) $ depend only on the radial coordinate r, with a periodic angular coordinate $ \theta \in[0,2 \pi] $. To solve equations for the vector and scalar KK modes, we consider one of the solutions:

          $ M(r)=\cosh (a r), \quad L(r)=\Re a \sinh (a r), $

          (64)

          where a is related to the Hubble constant $ H_0 $, and $ \Re $ depends on the brane tension. For convenience, we set $ \Re=1 $.

          Comparing the notation of the metric in (26), we obtain

          $ {\rm e}^{2 A}=M^2, \quad {\rm e}^{2 B_1}=1, \quad {\rm e}^{2 B_2}=L^2, $

          (65)

          where $ L(r)>0 $ for $ r>0 $; otherwise, $ {\rm e}^B=|L| $. The corresponding effective potentials for the KK modes are plotted in Fig. 2. Further, we compute the mass spectra for these modes, as summarized in Table 2.

          Figure 2.  (color online) Effective potentials for the KK modes in a 6D brane with the line element (63).

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $
          0 13.90 31.77 22.73 19.21 33.70 24.24 25.49 38.88 28.33 32.75 46.31 34.27
          1 45.96 81.81 63.67 54.99 84.59 66.02 65.37 91.95 72.31 76.67 102.35 81.23
          2 96.23 150.1 122.81 109.06 153.71 126.02 123.44 163.24 134.48 138.78 176.59 146.37
          3 164.73 236.62 200.18 181.34 241.07 204.23 199.7 252.73 214.85 219.07 269.02 229.70

          Table 2.  Mass spectra for the KK modes in a 6D brane with the line element (63).

          In this brane model, the corresponding coupling matrix for different l is provided in Table 3. We can calculate the masses of these modes by diagonalizing the mass matrix $ {\cal{M}}^2 $, which is shown in Table 4. The vector KK modes are heavier than the scalars except for $ l=0, n=0 $.

          $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ C^{(nm)} $ $\left( \begin{array}{cccc}9.54&7.20&6.22&5.58\\4.09&14.03&11.93&10.85\\5.11&7.90&18.74&16.70\\3.71&9.64&12.27&23.83 \end{array} \right)$ $\left( \begin{array}{cccc} 8.95&6.40&5.21&4.47\\3.30&12.73&10.34&8.87\\4.33&6.20&16.39&14.04\\2.73&7.94&9.25&19.97 \end{array} \right) $ $ \left( \begin{array}{cccc} 6.39&3.07&1.89&1.28\\0.29&8.25&5.05&3.46\\1.90&1.26&10.10&6.92\\0.19&3.49&2.50&11.90 \end{array} \right)$ $ \left( \begin{array}{cccc}7.27&4.17&2.94&2.23\\1.35&9.88&6.81&5.27\\2.58&3.14&12.51&9.35\\1.01&4.83&5.26&15.19 \end{array} \right)$

          Table 3.  Coupling coefficients of scalar-scalar couplings in a 6D brane with the line element (63).

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $
          0 13.90 15.33 19.21 18.00 25.49 25.24 32.75 30.67
          1 45.96 36.00 54.99 37.70 65.37 41.59 76.67 49.18
          2 96.23 55.83 109.06 59.61 123.44 69.45 138.78 77.48
          3 164.73 86.00 181.34 88.48 199.7 94.31 219.07 105.13

          Table 4.  Mass spectra for the vector and physical scalar KK modes in a 6D brane with the line element (63).

        • A.   Case 1

        • In Ref. [48] some 6D thick-brane models are considered from a cosmological perspective, where the metric is assumed as:

          $ d s^2=M^2(r) g_{\mu \nu} d x^\mu d x^\nu+d r^2+L(r)^2 d \theta^2 $

          (63)

          The metric functions $ M(r) $ and $ L(r) $ depend only on the radial coordinate r, with a periodic angular coordinate $ \theta \in[0,2 \pi] $. To solve equations for the vector and scalar KK mode, we consider one of the solutions :

          $ M(r)=\cosh (a r), \quad L(r)=\Re a \sinh (a r), $

          (64)

          where a is related to the Hubble constant $ H_0 $, and $ \Re $ depends on the brane tension. For convenience, we set $ \Re=1 $.

          Comparing the notation of the metric in (26), we obtain:

          $ e^{2 A}=M^2, \quad e^{2 B_1}=1, \quad e^{2 B_2}=L^2 $

          (65)

          where $ L(r)>0 $ for $ r>0 $; otherwise, $ e^B=|L| $. The corresponding effective potentials for the Kaluza–Klein (KK) modes are plotted in Figure 2. We also compute the mass spectra for these modes, as summarized in Table 2.

          Figure 2.  (color online) Effective potentials for the KK modes in a 6D brane with line element (63).

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $
          0 13.90 31.77 22.73 19.21 33.70 24.24 25.49 38.88 28.33 32.75 46.31 34.27
          1 45.96 81.81 63.67 54.99 84.59 66.02 65.37 91.95 72.31 76.67 102.35 81.23
          2 96.23 150.1 122.81 109.06 153.71 126.02 123.44 163.24 134.48 138.78 176.59 146.37
          3 164.73 236.62 200.18 181.34 241.07 204.23 199.7 252.73 214.85 219.07 269.02 229.70

          Table 2.  Mass spectra for the KK modes in a 6D brane with line element (63).

          In this brane model, the corresponding coupling matrix for different l is provided in Table 3. We can calculate the masses of these modes by diagonalizing the mass matrix $ {\cal{M}}^2 $, which is shown in Table 4. It suggests that the vector KK modes are heavier than the scalars except for $ l=0, n=0 $.

          $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ C^{(nm)} $ $\left( \begin{array}{cccc}9.54&7.20&6.22&5.58\\4.09&14.03&11.93&10.85\\5.11&7.90&18.74&16.70\\3.71&9.64&12.27&23.83 \end{array} \right)$ $\left( \begin{array}{cccc} 8.95&6.40&5.21&4.47\\3.30&12.73&10.34&8.87\\4.33&6.20&16.39&14.04\\2.73&7.94&9.25&19.97 \end{array} \right) $ $ \left( \begin{array}{cccc} 6.39&3.07&1.89&1.28\\0.29&8.25&5.05&3.46\\1.90&1.26&10.10&6.92\\0.19&3.49&2.50&11.90 \end{array} \right)$ $ \left( \begin{array}{cccc}7.27&4.17&2.94&2.23\\1.35&9.88&6.81&5.27\\2.58&3.14&12.51&9.35\\1.01&4.83&5.26&15.19 \end{array} \right)$

          Table 3.  Coupling coefficients of scalar-scalar couplings in a 6D brane with line element (63).

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $
          0 13.90 15.33 19.21 18.00 25.49 25.24 32.75 30.67
          1 45.96 36.00 54.99 37.70 65.37 41.59 76.67 49.18
          2 96.23 55.83 109.06 59.61 123.44 69.45 138.78 77.48
          3 164.73 86.00 181.34 88.48 199.7 94.31 219.07 105.13

          Table 4.  Mass spectra for the vector and physical scalar KK modes in a 6D brane with line element (63).

        • A.   Case 1

        • In Ref. [45], some 6D thick-brane models are considered from a cosmological perspective, where the metric is assumed as

          $ {\rm d} s^2=M^2(r) g_{\mu \nu} {\rm d} x^\mu {\rm d} x^\nu+ {\rm d} r^2+L(r)^2 {\rm d} \theta^2. $

          (63)

          The metric functions $ M(r) $ and $ L(r) $ depend only on the radial coordinate r, with a periodic angular coordinate $ \theta \in[0,2 \pi] $. To solve equations for the vector and scalar KK modes, we consider one of the solutions:

          $ M(r)=\cosh (a r), \quad L(r)=\Re a \sinh (a r), $

          (64)

          where a is related to the Hubble constant $ H_0 $, and $ \Re $ depends on the brane tension. For convenience, we set $ \Re=1 $.

          Comparing the notation of the metric in (26), we obtain

          $ {\rm e}^{2 A}=M^2, \quad {\rm e}^{2 B_1}=1, \quad {\rm e}^{2 B_2}=L^2, $

          (65)

          where $ L(r)>0 $ for $ r>0 $; otherwise, $ {\rm e}^B=|L| $. The corresponding effective potentials for the KK modes are plotted in Fig. 2. Further, we compute the mass spectra for these modes, as summarized in Table 2.

          Figure 2.  (color online) Effective potentials for the KK modes in a 6D brane with the line element (63).

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $ $ m_v^{2} $ $ m_{\phi}^{2} $ $ m_{\varphi}^{2} $
          0 13.90 31.77 22.73 19.21 33.70 24.24 25.49 38.88 28.33 32.75 46.31 34.27
          1 45.96 81.81 63.67 54.99 84.59 66.02 65.37 91.95 72.31 76.67 102.35 81.23
          2 96.23 150.1 122.81 109.06 153.71 126.02 123.44 163.24 134.48 138.78 176.59 146.37
          3 164.73 236.62 200.18 181.34 241.07 204.23 199.7 252.73 214.85 219.07 269.02 229.70

          Table 2.  Mass spectra for the KK modes in a 6D brane with the line element (63).

          In this brane model, the corresponding coupling matrix for different l is provided in Table 3. We can calculate the masses of these modes by diagonalizing the mass matrix $ {\cal{M}}^2 $, which is shown in Table 4. The vector KK modes are heavier than the scalars except for $ l=0, n=0 $.

          $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ C^{(nm)} $ $\left( \begin{array}{cccc}9.54&7.20&6.22&5.58\\4.09&14.03&11.93&10.85\\5.11&7.90&18.74&16.70\\3.71&9.64&12.27&23.83 \end{array} \right)$ $\left( \begin{array}{cccc} 8.95&6.40&5.21&4.47\\3.30&12.73&10.34&8.87\\4.33&6.20&16.39&14.04\\2.73&7.94&9.25&19.97 \end{array} \right) $ $ \left( \begin{array}{cccc} 6.39&3.07&1.89&1.28\\0.29&8.25&5.05&3.46\\1.90&1.26&10.10&6.92\\0.19&3.49&2.50&11.90 \end{array} \right)$ $ \left( \begin{array}{cccc}7.27&4.17&2.94&2.23\\1.35&9.88&6.81&5.27\\2.58&3.14&12.51&9.35\\1.01&4.83&5.26&15.19 \end{array} \right)$

          Table 3.  Coupling coefficients of scalar-scalar couplings in a 6D brane with the line element (63).

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $
          0 13.90 15.33 19.21 18.00 25.49 25.24 32.75 30.67
          1 45.96 36.00 54.99 37.70 65.37 41.59 76.67 49.18
          2 96.23 55.83 109.06 59.61 123.44 69.45 138.78 77.48
          3 164.73 86.00 181.34 88.48 199.7 94.31 219.07 105.13

          Table 4.  Mass spectra for the vector and physical scalar KK modes in a 6D brane with the line element (63).

        • B.   Case 2

        • Let us consider another 6D brane model presented in Ref. [46], which features football-shaped extra dimensions. The line element is assumed to be

          $ {\rm d} s^2=\eta_{\mu\nu} {\rm d}x^\mu {\rm d}x^\nu+\psi(r)({\rm d}r^2+r^2{\rm d}\phi^2). $

          (66)

          The solution for $ \psi(r) $ is $ \dfrac{4 \alpha^2 a_0^2}{r^2\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2} $, where $ r_0 $ is an arbitrary parameter, while α and $ a_0 $ are observable quantities related to the brane tension.

          Comparing this with the line element used in the present work, we obtain

          $\begin{aligned}[b]& A=0,\; B_1=\log\left(\frac{4 \alpha^2 a_0^2}{r^2\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2}\right), \\ & B_2=\log\left(\frac{4 \alpha^2 a_0^2}{\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2}\right).\end{aligned} $

          (67)

          Interestingly, effective potentials $ V_{\text{eff2}} $ and $ V_{\text{eff3}} $ coincide. Consequently, mass spectra for the two types of scalar KK modes are identical. Following the same procedure as in case 1, we calculate the physical masses of the scalar KK modes arising from their interactions; the results are presented in Table 5.

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $
          0 1.12 1.05 8.25 10.85 24.06 28.03 48.06 52.28
          1 12.36 4.21 25.28 11.33 48.33 28.04 80.25 52.28
          2 32.95 19.25 50.98 30.34 80.88 53.94 120.54 85.70
          3 62.80 22.84 85.79 31.45 121.90 54.01 168.96 85.71

          Table 5.  Mass spectra for the vector and physical scalar KK modes in the 6D brane with the line element (66).

          We observe that for $ n\neq 0 $, the vector modes remain heavier than the scalar modes. For larger values of l, some degenerate states of the scalars appear.

        • B.   Case 2

        • Let us consider another 6D brane model presented in Ref. [49], which features football-shaped extra dimensions. The line element is assumed to be

          $ ds^2=\eta_{\mu\nu}dx^\mu dx^\nu+\psi(r)(dr^2+r^2d\phi^2). $

          (66)

          The solution for $ \psi(r) $ is $ \dfrac{4 \alpha^2 a_0^2}{r^2\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2} $, where $ r_0 $ is an arbitrary parameter, while α and $ a_0 $ are observable quantities related to the brane tension.

          Comparing this with the line element used in the present work, we obtain

          $\begin{aligned}[b]& A=0,\; B_1=\log\left(\frac{4 \alpha^2 a_0^2}{r^2\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2}\right), \\& B_2=\log\left(\frac{4 \alpha^2 a_0^2}{\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2}\right).\end{aligned} $

          (67)

          Interestingly, in this case the effective potentials $ V_{\text{eff2}} $ and $ V_{\text{eff3}} $ coincide. Consequently, the mass spectra for the two types of scalar KK modes are identical. Following the same procedure with case 1, we can calculate the physical masses of the scalar KK modes arising from their interactions; the results are presented in Table 5.

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $
          0 1.12 1.05 8.25 10.85 24.06 28.03 48.06 52.28
          1 12.36 4.21 25.28 11.33 48.33 28.04 80.25 52.28
          2 32.95 19.25 50.98 30.34 80.88 53.94 120.54 85.70
          3 62.80 22.84 85.79 31.45 121.90 54.01 168.96 85.71

          Table 5.  Mass spectra for the vector and physical scalar KK modes in the 6D brane with the line element (66).

          We observe that for $ n\neq 0 $, the vector modes remain heavier than the scalar modes. For larger values of l, some degenerate states of the scalars appear.

        • B.   Case 2

        • Let us consider another 6D brane model presented in Ref. [46], which features football-shaped extra dimensions. The line element is assumed to be

          $ {\rm d} s^2=\eta_{\mu\nu} {\rm d}x^\mu {\rm d}x^\nu+\psi(r)({\rm d}r^2+r^2{\rm d}\phi^2). $

          (66)

          The solution for $ \psi(r) $ is $ \dfrac{4 \alpha^2 a_0^2}{r^2\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2} $, where $ r_0 $ is an arbitrary parameter, while α and $ a_0 $ are observable quantities related to the brane tension.

          Comparing this with the line element used in the present work, we obtain

          $\begin{aligned}[b]& A=0,\; B_1=\log\left(\frac{4 \alpha^2 a_0^2}{r^2\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2}\right), \\ & B_2=\log\left(\frac{4 \alpha^2 a_0^2}{\left[(r / r_0)^\alpha+(r / r_0)^{-\alpha}\right]^2}\right).\end{aligned} $

          (67)

          Interestingly, effective potentials $ V_{\text{eff2}} $ and $ V_{\text{eff3}} $ coincide. Consequently, mass spectra for the two types of scalar KK modes are identical. Following the same procedure as in case 1, we calculate the physical masses of the scalar KK modes arising from their interactions; the results are presented in Table 5.

          n $ l=0 $ $ l=1 $ $ l=2 $ $ l=3 $
          $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $ $ m_v^{2} $ $ {\cal{M}}^{2(n)}_A $
          0 1.12 1.05 8.25 10.85 24.06 28.03 48.06 52.28
          1 12.36 4.21 25.28 11.33 48.33 28.04 80.25 52.28
          2 32.95 19.25 50.98 30.34 80.88 53.94 120.54 85.70
          3 62.80 22.84 85.79 31.45 121.90 54.01 168.96 85.71

          Table 5.  Mass spectra for the vector and physical scalar KK modes in the 6D brane with the line element (66).

          We observe that for $ n\neq 0 $, the vector modes remain heavier than the scalar modes. For larger values of l, some degenerate states of the scalars appear.

        V.   SUMMARY AND DISCUSSION
        • In this study, we derived a Stückelberg-like effective action for the massive KK modes of the gauge field by incorporating a $(\nabla^M X_M)^2$ term for a bulk $U(1)$ gauge field. The gauge invariance of the effective action is ensured by an imposed condition on the brane-localized gauge freedom. This construction deviates from the conventional Stückelberg mechanism in two critical aspects: (1) The scalar KK modes are not introduced ad hoc but emerge naturally from the extra-dimensional components of the bulk $U(1)$ gauge field. (2) The scalar KK modes can acquire masses from the entire geometric structure. Additional scalar modes arise when the brane possesses more than one extra dimension.

          We first examined the scenario within a five-dimensional brane-world framework. Numerical analyses of this model revealed that the massive vector KK modes were heavier than those of the scalar KK modes. Extending this analysis to 6D brane models, two distinct types of scalar KK modes emerge. We focused specifically on the scalar–scalar mixed couplings induced by the geometric background. These mixed couplings give rise to mass oscillations among the scalar modes. We obtained the physical masses of the scalar sector by diagonalizing the mass matrix.

          In summary, the discussed scalar KK modes are distinguishable from conventional Stückelberg scalars and Goldstone bosons, which can be attributed to their unique coupling structure with vector modes and among themselves. For future work, it would be intriguing to further explore the properties of these scalar modes and their mixed couplings within the framework of quantum field theory, with specific emphasis on their phenomenological implications.

        V.   SUMMARY AND DISCUSSION
        • In this study, we derived a Stückelberg-like effective action for the massive KK modes of the gauge field by incorporating a $(\nabla^M X_M)^2$ term for a bulk $U(1)$ gauge field. The gauge invariance of the effective action is ensured by an imposed condition on the brane-localized gauge freedom. This construction deviates from the conventional Stückelberg mechanism in two critical aspects: (1) The scalar KK modes are not introduced ad hoc but emerge naturally from the extra-dimensional components of the bulk $U(1)$ gauge field. (2) The scalar KK modes can acquire masses from the entire geometric structure. Additional scalar modes arise when the brane possesses more than one extra dimension.

          We first examined the scenario within a five-dimensional brane-world framework. Numerical analyses of this model revealed that the massive vector KK modes were heavier than those of the scalar KK modes. Extending this analysis to 6D brane models, two distinct types of scalar KK modes emerge. We focused specifically on the scalar–scalar mixed couplings induced by the geometric background. These mixed couplings give rise to mass oscillations among the scalar modes. We obtained the physical masses of the scalar sector by diagonalizing the mass matrix.

          In summary, the discussed scalar KK modes are distinguishable from conventional Stückelberg scalars and Goldstone bosons, which can be attributed to their unique coupling structure with vector modes and among themselves. For future work, it would be intriguing to further explore the properties of these scalar modes and their mixed couplings within the framework of quantum field theory, with specific emphasis on their phenomenological implications.

        V.   SUMMARY AND DISCUSSION
        • In this work, by incorporating a $(\nabla^M X_M)^2$ term for a bulk $U(1)$ gauge field, we derived a Stückelberg-like effective action for the massive Kaluza-Klein (KK) modes of the gauge field. Gauge invariance of the effective action is ensured by an imposed condition on the brane-localized gauge freedom. This construction deviates from the conventional Stückelberg mechanism in two critical aspects: first, the scalar KK modes are not introduced ad hoc but emerge naturally from the extra-dimensional components of the bulk $U(1)$ gauge field; second, the scalar KK modes can acquire masses from the entire geometric structure. Notably, when the brane possesses more than one extra dimension, additional scalar modes arise.

          We first examined the scenario within a five-dimensional brane-world framework. Numerical analyses of this model reveal that the massive vector KK modes are heavier than the scalar KK modes. Extending this analysis to six-dimensional brane models, two distinct types of scalar KK modes emerge. We focused specifically on the scalar–scalar mixed couplings induced by the geometric background. These mixed couplings give rise to mass oscillations among the scalar modes; by diagonalizing the mass matrix, we obtained the physical masses of the scalar sector.

          In summary, the scalar KK modes discussed herein are distinguishable from both conventional Stückelberg scalars and Goldstone bosons, attributable to their unique coupling structure with vector modes and among themselves. For future work, it would be intriguing to further explore the properties of these scalar modes and their mixed couplings within the framework of quantum field theory, with particular emphasis on their phenomenological implications.

      Reference (46)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return