1 Introduction

Braneworld scenarios [1, 2], show a new viewpoint of spacetime and enables a new approaches to explain a large number of outstanding issues such as the hierarchy problem [2], the cosmological problem [3], the nature of dark matter [4] and dark energy. Furthermore, by assuming a warped geometry, the propagation of the gravitational field [5, 6] and the gauge field, [7], as well as fermionic fields [8], are governed by the bulk curvature in general relativity (GR). An equivalent theory is the known teleparallel equivalent of general relativity (TEGR) [9,10,11,12,13,14,15,16,17], that is constructed in the Weitzenböck manifold, which has vanishing curvature but nonvanishing torsion. Equivalence with GR is provided by the ratio of the torsion T scalar and the Ricci R scalar, which is the boundary term B, making TEGR have the same field equations as GR.

The localization mechanism employed for matter fields living in a 5D braneworld scenarios has been the subject of many studies. The study of fermion localization on branes is rich and interesting, yet the most popular method for the localization of fermions is formulated in a rather speculative way [8, 18,19,20,21,22,23,24,25]. The same is true in 6D [26,27,28]. This is because of the freedom one has to propose the Yukawa coupling term. The location of the fermion was studied in several modified gravity models, such as gravity f(R) [29,30,31] and gravity f(T) [32, 33].

A new teleparallel gravity model is the f(TB) gravity, where B is the boundary term [34,35,36], that has attracted a lot of attention due this model features have, as well as good agreement with observational data to describe the accelerated expansion of the universe [37, 38], and their significant results in cosmological perturbations and thermodynamics, and dark energy, and gravitational waves [39,40,41,42,43,44,45]. Furthermore, the gravity f(TB) was studied in a brane scenario, where it was possible to observe that the additional term B induces changes on the energy density causing a split in the brane, also changing the gravitational perturbations [46].

Inspired on the results obtained in [32, 46], we investigate the issue of fermion localization in f(TB) gravity. In Sect. 2 we review the main definitions of the teleparallel f(TB) theory and build the respective braneworld. In the Sect. 3, we obtain the solutions of the system and we examined the energy density components of the brane. The Sect. 4 deals with the fermionic sector of the model using the Yukawa coupling. Finally, additional comments and results are discussed in Sect. 5.

2 Metric equations

In teleparallel gravity, the vielbein, \(h^{{\overline{M}}}\ _M\) (rather than the metric) are the actual gravitational dynamic variables. We will use the Latin letters (\(M , N , Q ,\ldots =0, 1, 2, 3, 4\)) for the indices related to the bulk, and barred indices (\({\overline{M}}, {\overline{N}}, {\overline{Q}},\ldots =0, 1, 2, 3, 4\)) for indices related to tangent space. We assume a mostly plus metric signature, i.e., \(diag(-1, 1, 1, 1,1)\).

The relevant connection for TEGR is the so-called Weitzenböck connection. An important feature of this connection is that the corresponding spin connection vanishes identically. Thus, the Weitzenböck connection is represented as \({\widetilde{\varGamma }}^P\ _{NM}=h_{{\overline{M}}}\ ^P\partial _M h^{{\overline{N}}}\ _M\), that is, \(\nabla _Q h^{{\overline{M}}}\ _M\equiv \partial _Q h^{{\overline{M}}}\ _M- {\widetilde{\varGamma }}^P\ _{QM} h^{{\overline{M}}}\ _P =0\), which is called the condition of absolute parallelism [14]. The Weitzenböck and Levi-Civita connections are related by

$$\begin{aligned} {\widetilde{\varGamma }}^P\ _{NM}= \varGamma ^P\ _{NM} + K^P\ _{NM}, \end{aligned}$$
(1)

where \(\varGamma ^P\ _{NM}\) is the Levi-Civita connection of general relativity, and

$$\begin{aligned} K^P\ _{NM}=( T_N\ ^P\ _M +T_M\ ^P\ _N - T^P\ _{NM})/2 \end{aligned}$$
(2)

is defined as the contortion tensor of the Weitzenböck connection [14]. We take the torsion as \(T^{P}\ _{NM}= {\widetilde{\varGamma }}^P\ _{MN}-{\widetilde{\varGamma }}^P\ _{NM}\), and we also define a dual torsion tensor, known as a superpotential \(S_{P}\ ^{NM}=( K^{NM}\ _{P}-\delta ^M_P T^{QN}\ _Q+\delta ^N_P T^{QM}\ _Q)/2\) [14]. Therefore, the Lagrangian of TEGR reads

$$\begin{aligned} {\mathscr {L}}=- hT/4 , \end{aligned}$$
(3)

where \(T=T_{PMN}S^{PMN}\) is the torsional scalar, \(h=\sqrt{-g}\) and \(c^4/4\pi G=1\) for simplicity [14]. On the other hand, the Riemann tensor in the Levi-Civita connection

$$\begin{aligned} R^P\ _{MQN}= & {} \partial _Q\varGamma ^P\ _{MN}-\partial _N\varGamma ^P\ _{MQ}\nonumber \\&+\varGamma ^P\ _{S Q}\varGamma ^S\ _{MN}-\varGamma ^P\ _{SN}\varGamma ^S\ _{MQ}. \end{aligned}$$
(4)

From the relation between the Weitzenböck connection and the Levi-Civita connection given by Eq. (1), one can write the Riemann tensor in the form

$$\begin{aligned} R^P\ _{MQN}= & {} \nabla _N K^P\ _{MQ}-\nabla _Q K^P\ _{MN}\nonumber \\&+K^P\ _{SN}K^S\ _{MQ}-K^P\ _{SQ}K^Q\ _{MN}, \end{aligned}$$
(5)

whose associated Ricci tensor can then be written as

$$\begin{aligned} R_{MN}= & {} \nabla _N K^P\ _{MP}-\nabla _P K^P\ _{MN}\nonumber \\&+K^P\ _{SN}K^S\ _{MP}-K^P\ _{SP}K^S\ _{MN}. \end{aligned}$$
(6)

Using the relations \(K^{(MN)S}=T^{M(NS)}=S^{M(NS)}\) and considering that \(S^M\ _{PM}=2K^M\ _{PM}=-2T^M\ _{PM}\) along with Eq. (2), one can get [47,48,49]

$$\begin{aligned} R _{MN}=-\nabla ^P S _{NPM}-g_{MN}\nabla ^P T^S\ _{PS}-S^{PS}\ _{N}\ K_{SPN}. \end{aligned}$$
(7)

In turn, Ricci scalar is

$$\begin{aligned} R=-T-2\nabla ^{M}T^{N}\ _{MN}. \end{aligned}$$
(8)

Thus we can identify the boundary term

$$\begin{aligned} B\equiv -2\nabla ^{M}T^{N}\ _{MN}=\frac{2}{h}\partial _M(h T^M), \end{aligned}$$
(9)

in which \(T^M\) is the torsion tensor that can be defined as \(T_M=T^N\ _{MN}\). We can easily see that GR and TEGR will lead to exactly the same equations [14]. However, this will not be the case if one uses f(R) or f(T) as the Lagrangian of the theory, which therefore corresponds to different gravitational modifications [45]. However, when we consider f(TB) as the Lagrangian of the theory, we have that \(f(T,B)=f(-T+B)\) is the teleparallel equivalent of f(R).

We can consider a modified gravity theory where the gravitational Lagrangian depends on T and B [45]. Therefore, we have a gravitational action to f(TB), namely

$$\begin{aligned} {\mathscr {S}}=-\frac{1}{4}\int h f(T,B)d^5x+\int {\mathscr {L}}_m d^5x, \end{aligned}$$
(10)

where \({\mathscr {L}}_m\) is the matter Lagrangian. We can get the field equations by varying the action with respect to the vielbein [41, 45]

$$\begin{aligned}&\frac{1}{h}f_T\Big [\partial _Q(h S_N\ ^{MQ})-h{\widetilde{\varGamma }}^R\ _{SN}S_R\ ^{MS}\Big ]+\frac{1}{4}\Big [f-Bf_B\Big ]\delta _N^M \nonumber \\&\quad +\Big [(\partial _Qf_T)+(\partial _Qf_B) \Big ]S_N\ ^{MQ}\nonumber \\&\quad +\frac{1}{2}\Big [\nabla ^M\nabla _N f_B-\delta ^M_N\Box f_B\Big ]=-{\mathscr {T}}_N\ ^M, \end{aligned}$$
(11)

where \(\Box \equiv \nabla ^M\nabla _M\), \(f\equiv f(T,B)\), \(f_T\equiv \partial f(T,B)/\partial T\) e \(f_{B}\equiv \partial f(T,B)/\partial B\) and \({\mathscr {T}}_N\ ^M\) is the stress-energy tensor, which in terms of the matter Lagrangian is given by \({\mathscr {T}}_a\ ^M=-\delta {\mathscr {L}}_m/\delta h^a\ _M\). The matter Lagrangian density is taken as

$$\begin{aligned} {\mathscr {L}}_m=-h\left[ \frac{1}{2}\partial ^M\phi \partial _M\phi +V(\phi )\right] , \end{aligned}$$
(12)

where \(\phi \equiv \phi (y)\) is a background scalar field that generates the brane.

In our work, we consider the static codimension one braneworld scenario whose metric can be written as

$$\begin{aligned} ds^2=e^{2A(y)}\eta _{\mu \nu }dx^\mu dx^\nu +dy^2, \end{aligned}$$
(13)

where \(e^{A(y)}\) is the warped factor. We can choose the vielbein in the form

$$\begin{aligned} h_{{\overline{M}}}\ ^M=diag(e^A, e^A, e^A, e^A, 1). \end{aligned}$$
(14)

Using the Weitzenböck connection, the torsional scalar and the boundary term are given by

$$\begin{aligned} T=-12A'^2&,&B=-8(A''+4A'^2), \end{aligned}$$
(15)

where the prime \(( ' )\) denotes differentiation with respect to y.

Thus, the gravitational field equations are given as

$$\begin{aligned}&\phi ''+4A'\phi =\frac{d V}{d\phi }, \end{aligned}$$
(16)
$$\begin{aligned}&\frac{1}{4}\Big [f+8(A''+4A'^2)f_B\Big ]+6A'^2f_T= \frac{1}{2}\phi ^2-V, \end{aligned}$$
(17)
$$\begin{aligned}&12\Big [A'(A'''+8A'A'')(f_{BB}+f_{TB})+3A''A'^2(f_{TT}+f_{BT})\Big ]\nonumber \\&\quad -\frac{1}{2}(A''+4A'^2)(4f_{B}+3f_{T})+\frac{1}{4}f= \frac{1}{2}\phi ^2+V. \end{aligned}$$
(18)

We can rewrite Eqs. (17) and (18) as

$$\begin{aligned} 6A'^2= & {} -\frac{1}{f_T}(P+P_{TB} ), \end{aligned}$$
(19)
$$\begin{aligned} 3A''+12A'^2= & {} -\frac{2}{f_T}(\rho +\rho _{TB}), \end{aligned}$$
(20)

where

$$\begin{aligned} P_{TB}= & {} \frac{1}{4}\Big [f+8(A''+4A'^2)f_B\Big ], \end{aligned}$$
(21)
$$\begin{aligned} \rho _{TB}= & {} -12\Big [A'(A'''+8A'A'')(f_{BB}+f_{TB})\nonumber \\&+3A''A'^2(f_{TT}+f_{BT})\Big ]\nonumber \\&+2(A''+4A'^2)f_{B}+\frac{1}{4}f. \end{aligned}$$
(22)

Note that the left side of Eqs. (19) and (20) is equivalent to that obtained in TEGR. So we can state that modified gravity equations of the motion of the f(TB) gravities are similar to an inclusion of an additional source with \(\rho _{TB}\) and \(P_{TB}\).

The diagonal tetrad (14) represents a good choice among all the possible vielbeins giving metric (13). In fact, the gravitational field equations do not involve any additional constraints on the function f(TB) or the scalars T and B. Thus, the choice in Eq. (14) can be regarded as a “good vielbein”. Similarly, in the FRW cosmological models the f(TB) gravitational dynamics preserves the form of the usual Friedmann equations (two equations) [37,38,39,40,41,42]. The choice of vielbein is a rather important issue, as it fixes the number of degrees of freedom of the theory, as seen particularly in a gravitational wave analysis of f(TB) gravity [50].

3 Thick brane solutions

Since the Eqs. (17) and (18) form a second-order derivative system, it is difficult to provide an analytical solution for this case. For simplicity, we can take an ansatz [51]

$$\begin{aligned} e^{2A(y)}=\cosh ^{-2p}(\lambda y), \end{aligned}$$
(23)

where the p parameter modifies the warp variation within the brane core, and \(\lambda \) determines the brane width.

We will then propose the cases where \(f_1(T,B)=T+k_1B^{n_1}\) and \(f_2(T,B)=B+k_2T^{n_2}\) , where \(k_{1,2}\) and \(n_{1,2}\) are parameters controlling the deviation of the usual teleparallel theory [46].

We follow the approach carried out in Ref. [46], where by manipulating the (17) and (18) equations, an equation relating the metric components and the scalar field is obtained. In this case, for \(f_1(T,B)\) with Eqs. (17), (18) and (23), we get [46]

$$\begin{aligned} \phi '^2(y)= & {} \frac{3}{2}p\lambda ^2\mathrm {sech}^2(\lambda y)\nonumber \\&-\frac{\alpha ^{n_1}}{\beta ^2}\Big [8^{n_1-1}k_1n_1(n_1-1)(1+4p)\sinh ^2(\lambda y)\Big ],\nonumber \\ \end{aligned}$$
(24)
$$\begin{aligned} V(\phi (y))= & {} \frac{3}{4}\alpha +2^{3n_1-4}k_1(n_1-1)\nonumber \\&\quad (p\lambda ^2)^2\alpha ^{n_1-2}\Big \{4 \mathrm {sech}^4(\lambda y)+64p^2\tanh ^4(\lambda y)\nonumber \\&-[3n_1+4(8+3n_1)]\mathrm {sech}^2(\lambda y)\tanh ^2(\lambda y)\Big \}.\nonumber \\ \end{aligned}$$
(25)

We can solve Eq. (24) to find a function \(\phi =g(y)\) that may be inverted to give \(y=g^{-1}(\phi )\), which allows us to write the potential in the usual way \(V=V(\phi )\). The thick brane solution for \(n_1=1\) is the same obtained in Refs. [51, 52], which does not depend on the \(k_1\) parameter. For \(n_1=2\) we get as a solution the first and second kind elliptic integrals, which depends on parameter \(k_1\) [46].

In Fig. 1, we plotted the \(\phi (y)\) field for \(f_1(T,B)\). The \(n_1=1\) configuration (Fig. 1a) is a kink solution. For \(n_1=2\) configuration, for a decreasing value of \(k_1\), the solution goes from kink to double-kink, as shown in Fig. 1b. This feature reflects the brane internal structure, which tends to split the brane. A similar result was obtained in Ref. [32].

Fig. 1
figure 1

The shape of the scalar \(\phi (y)\) for \(f_1(T,B)\), where \(p=\lambda =1\). a For \(n_1=1\). b For \(n_1=2\)

For \(f_2(T,B)\), we get

$$\begin{aligned} \phi '^2(y)= & {} -2^{2n_2-3}(-3)^{n_2}p^{-1}k_2n_2\nonumber \\&\times (2n_2-1){\mathrm {csch}}^2(\lambda y)[p\lambda \tanh (\lambda y)]^{2n_2}, \end{aligned}$$
(26)
$$\begin{aligned} V(\phi (y))= & {} 4^{n_2-2}(-3)^{n_2}p^{-1}k_2\nonumber \\&\times (2n_2-1)[4p-n_2{\mathrm {csch}}^2(\lambda y)][p\lambda \tanh (\lambda y)]^{2n_2}.\nonumber \\ \end{aligned}$$
(27)

The setting \( n_2 = 4, 6, ... \) (even numbers) does not present a pleasant solution.

In Fig. 2, we plotted the \(\phi (y)\) field for \(f_2(T,B)\) varying the parameter \(k_2\). For \(n_2=1\) configuration (Fig. 2a) we have a kink solution, whereas for \(n_2=3\) configuration we have a double-kink solution (Fig. 2b). Again, this feature reflects the brane internal structure, which tends to split the brane.

Fig. 2
figure 2

The shape of the scalar \(\phi (y)\) for \(f_2(T,B)\), where \(p=\lambda =1\). a For \(n_2=1\). b For \(n_2=3\)

The energy densities for \(f_1(T,B)\) are [46]

$$\begin{aligned} \rho _1(y)= & {} \alpha \Bigg [2^{3n_1-2}(n_1-1)k_1\nonumber \\&-\frac{8^{n_1-1}(n_1-1)3n_1k_1(1+4p)\sinh ^2(\lambda y)}{\beta ^2}\Bigg ]\nonumber \\&-3\Big [(p\lambda )^2-\frac{1}{2}(1+2p)p\lambda ^2\mathrm {sech}^2(\lambda y)\Big ], \end{aligned}$$
(28)

where we defined the functions \(\alpha \equiv p \lambda ^2[\mathrm {sech}^2(\lambda y)-4p\tanh ^2(\lambda y)]\), and \(\beta \equiv 1+2[1-\cosh (2\lambda y)]p\). In Fig. 3, we plotted the energy densities \(\rho _1(y)\) for \(f_1(T,B)\), with \(n_1=1\) (Fig. 3a) and \(n_1=2\) (Fig. 3b), that includes a new peak varying the parameter \(k_1\).

Fig. 3
figure 3

Energy density in the brane for \(f_1(T,B)\), where \(p=\lambda =1\). a For \(n_1=1\). b For \(n_1=2\)

The energy densities for \(f_2(T,B)\) are [46]

$$\begin{aligned}&\rho _2(y)=-2^{2n_2-3}(-3)^{n_2}p^{-1}k_2(2n_2-1)\beta \mathrm {csch}^2(\lambda y)\nonumber \\&\quad \times [p\lambda \tanh (\lambda y)]^{2n_2}. \end{aligned}$$
(29)

In Fig. 4, we plotted the energy densities \(\rho _2(y)\) for \(f_2(T,B)\) varying the parameter \(k_2\), for \(n_2=1\) (Fig. 4a) and \(n_2=3\) (Fig. 4b), which has two peaks.

Fig. 4
figure 4

Energy density in the brane for \(f_2(T,B)\), where \(p=\lambda =1\). a For \(n_2=1\). b For \(n_2=3\)

4 Spin 1/2 fermions

In this section, we explore the effects of the modified teleparallel f(TB) on the matter (fermionic) sector. We changed the variable from y to z in the metric (13), and so, \(dz=e^{-A(y)}dy\) and the metric turns to \(ds^2=e^{2A}(\eta ^{\mu \nu }dx^\mu dx^\nu +dz^2)\). Considering a well-known Yukawa coupling between the fermion and the scalar field \(\phi \), the 5-dimensional Dirac action of a spin 1/2 fermion minimally coupled to the gravity and to the background scalar \(\phi \) is

$$\begin{aligned} {\mathscr {S}}_{1/2}=\int h {\overline{\varPsi }} \Big (\varGamma ^M D_M\varPsi -\xi \phi \varPsi \Big )d^5x, \end{aligned}$$
(30)

where \(\varGamma ^M=h_{{\overline{M}}}\ ^M \varGamma ^{{\overline{M}}}\) are the Dirac curved matrices defined from the Dirac flat matrices \(\varGamma ^{{\overline{M}}}\) through the vielbeins. These matrices obey Clifford algebra \(\{\varGamma ^M,\varGamma ^N\}=2g^{MN}\). \(D_M\) is the covariant derivative given by \( D_M=\partial _M +\varOmega _M\), where [24, 25]

$$\begin{aligned} \varOmega _M=\frac{1}{4}\Big (K_M\ ^{{{\overline{N}}}{{\overline{Q}}}}\Big )\ \varGamma _{{\overline{N}}}\varGamma _{{\overline{Q}}}, \end{aligned}$$
(31)

is the spin connection, which for our case is such that \(\varOmega _\mu =\frac{1}{4}(-\partial _z A)\varGamma _{\mu }\varGamma ^{z}\) and \(\varOmega _z=\partial _z A\). By choosing the spinor representation [8, 28, 53]

$$\begin{aligned}&\varPsi \equiv \varPsi (x,z)=\left( \begin{array}{cccccc} \psi \\ 0\\ \end{array}\right) ,\quad \varGamma ^{{\overline{\mu }}}=\left( \begin{array}{cccccc} 0&{}\gamma ^{{\overline{\mu }}}\\ \gamma ^{{\overline{\mu }}}&{}0\\ \end{array}\right) ,\nonumber \\&\varGamma ^{{\overline{z}}}=\left( \begin{array}{cccccc} 0&{}\gamma ^4\\ \gamma ^4&{}0\\ \end{array}\right) , \end{aligned}$$
(32)

the Dirac equation takes the form

$$\begin{aligned} \Big [\gamma ^{\mu }\partial _\mu +\gamma ^4\partial _z-\xi e^A\phi \Big ]\psi =0. \end{aligned}$$
(33)

We apply a decomposition to the spinor \(\psi =\sum _n[\psi _{L,n}(x)\varphi _{L,n}(z)+\psi _{R,n}(x)\varphi _{R,n}(z)]\), being \(\gamma ^4\psi _{R,L}=\pm \psi _{R,L}\) e \(\gamma ^\mu \partial _\mu \psi _{R,L}=m\psi _{L,R}\). So, we have the coupled equations

$$\begin{aligned}&\Big [\partial _z+\xi e^A \phi \Big ]\varphi _{L}(z)=m \varphi _{R}(z),\nonumber \\&\Big [\partial _z-\xi e^A\phi \Big ]\varphi _{R}(z)=m \varphi _{L}(z). \end{aligned}$$
(34)

These equations can be decoupled and reduced to Schrödinger-like equations

$$\begin{aligned}&\Big [-\partial ^2_z+V_L(z)\Big ]\varphi _{L}(z)=m^2 \varphi _{L}(z),\nonumber \\&\Big [-\partial ^2_z+V_R(z)\Big ]\varphi _{R}(z)=m^2 \varphi _{R}(z), \end{aligned}$$
(35)

where

$$\begin{aligned} V_L(z)= & {} U^2 -\partial _{z}U,\nonumber \\ V_R(z)= & {} U^2 +\partial _{z}U, \end{aligned}$$
(36)

and \(U=\xi e^A \phi \) is the so-called superpotential. The supersymmetric structure of the potentials (36) leads to a massless mode in the form

$$\begin{aligned} \varphi _{R0,L0}(z)\propto \exp {\Bigg [\pm \int \xi \phi e^{A}dz\Bigg ]}, \end{aligned}$$
(37)

where \(\phi e^{A}|_{z\rightarrow \pm \infty }\rightarrow 0\). In other words, the zero mode for fermions can be localized on the brane for positive \(\xi \) [32].

For \(f_1(T,B)\) with \(n_1=1\) only left-chiral fermions can be localized on the brane. For \(n_1=2\) only right-chiral fermions can be localized on the brane. In this case note that the smaller the \(k_1\) parameter, the more localized the mode becomes, as can be seen in the Fig. 5. For \(f_2(T,B)\) with \(n_2=1\) only left-chiral fermions can be localized on the brane. Now, the higher the parameter \( k_1 \), the more localized the mode becomes. On the other hand, for \(n_1=3\) only right-chiral fermions can be localized on the brane, and so, the smaller the \(k_1\) parameter, the more localized the mode becomes, as can be seen in the Fig. 6.

Fig. 5
figure 5

\(V_L\) (a) and \(\varphi _{L}\) (b) for \(f_1(T,B)\) with \(n_1=1\). \(V_R\) (c) and \(\varphi _{R}\) (d) for \(f_1(T,B)\) with \(n_1=2\) (\(p=\lambda =\xi =1\))

Fig. 6
figure 6

\(V_L\) (a) and \(\varphi _{L}\) (b) for \(f_2(T,B)\) with \(n_2=1\). \(V_R\) (c) and \(\varphi _{R}\) (d) for \(f_2(T,B)\) with \(n_2=3\) (\(p=\lambda =\xi =1\))

Note that the effective potentials are even functions, so the wave functions will be either even or odd. We can analyze numerically the Eq. (35). For that we impose the following conditions: \(\varphi _{even}(0)=c\), \(\partial _z\varphi _{even}(0)=0\), \(\varphi _{odd}(0)=c\), and \(\partial _z\varphi _{odd}(0)=0\), where c is a constant [8, 19, 20]. Here \(\varphi _{even}\) and \(\varphi _{odd}\) denote the even and odd parity modes of \(\varphi _{R,L}(z)\), respectively.

Fig. 7
figure 7

Non-normalized massive fermionic modes for \(f_1(T,B)\), with \(n_1=2\) and \(k_1=-0.5\). For \(\varphi _{even}\) (a) and \(\varphi _{odd}\) (b). By varying \(k_1\), \(\varphi _{even}\) with \(m=3.492\) (c) and \(\varphi _{odd}\) with \(m=4.233\) (d). For \(f_2(T,B)\) with \(n_2=3\), \(\varphi _{even}\) with \(m=6.72\) (e) and \(\varphi _{odd}\) with \(m=6.711\) (f) (\(p=\lambda =\xi =1\))

As depicted in Fig. 7, the asymptotic divergence of the massive modes shows that they form non-localized states, which is a behavior typical of plane wave oscillations, characteristic of a free mode. This shows that these modes represent massive fermions that certainly will be leaked from the brane.

For \(f_1(T,B)\), both for \(n_1=1\) and \(n_1=2\), the greater the mass, the more oscillations we obtain, as can be seen in the Fig. 7a, b for \(n_1=2\). In Fig. 7c, d, we observe that when decreasing the value of \(k_1\), the greater the amplitude of the oscillation, mainly near the brane. For \(f_2(T,B)\), both for \(n_2=1\) and \(n_2=3\), the greater the mass, the more oscillations we obtain. Increasing the value of \(k_2\), greater the amplitude of the oscillation, mainly near the brane as we can see in the Fig. 7e, f for \(n_2=3\).

5 Final remarks

In this work we considered a braneworld in the context of the f(TB) modified teleparallel gravity constructed with one scalar field. We propose two particular cases for f(TB), namely \(f_1(T,B)=T+k_1B^{n_1}\) and \(f_2(T,B)=B+k_2T^{n_2}\). In both cases the torsion and boundary term produce an inner brane structure tending to split the brane. We also find that the \(n_{1,2}\) and \(k_{1,2}\) parameters determines whether the domain wall solution is a kink or double-kink. For \(f_1(T,B)\) where \(n_1=2\), with the decreasing of the contribution of \(k_1\), the configuration of the solution changes from a kink to double-kink. The same is true for \(f_2(T,B)\) where \(n_2=3\), when we increase the value of parameter \(k_2\). The thick brane undergoes a phase transition evinced by the energy density components. Similar behavior was found for f(T) in Ref. [32].

We considered a simple Yukawa coupling between the scalar and the spinor field. We notice that potentials feel the division of brane when we vary \(n_{1,2}\) and \(k_{1,2}\), the same happens with the zero modes, which become more localized. We note that for \(f_1(T,B)\) where \(n_1=1\), only left-chiral fermions are located, the same is true for \(f_2(T,B)\) with \(n_2=1\). For \(f_1(T,B)\) where \(n_1=2\), only right-chiral fermions are located, the same is true for \(f_2(T,B)\) with \(n_2=3\). The massive fermionic modes are dependent on the parameters that control torsion and boundary term. This is well evidenced for \(f_1(T,B)\) with \( n_1 = 3 \) since that decreasing the value of \(k_1\), increases the amplitude of the ripples making them more intense and presenting ripples within the brane. The same goes for \(f_2(T,B)\) when increasing the value of \(k_2\), which is very evident for \(n_2=3\). Therefore, the brane splitting process leads to modifications of the massive fermionic modes inside the thick brane. The interaction of the massive modes with the torsion and boundary term is more intense in the brane core where the amplitude and the rate of growth depend on the parameters \(n_{1,2}\) and \(k_{1,2}\).

Although only one massless chiral mode was found for each configuration \(n_{1,2}\) and \(k_{1,2}\), only left-handed spin 1/2 fermions were detected so far. The configurations where the right-handed massless mode is localized on the brane are beyond the standard model states. The absence of left-handed massless mode can be used to rule out those configurations where only right-handed massless mode are trapped.

In addition, it is worthwhile to mention the role played by the parameter \(k_{1,2}\) on the brane internal structure. As \(k_{1,2}\) grows the brane undergoes a transition from a single into a two-brane. Therefore, the parameter \(k_{1,2}\) can be regarded as phase transition parameter controlling the brane splitting process.