turn out to be at the TeV scale.\
Using the solution of $\zeta^n_{L,R}(y)$ (in eqn. (\[solution3\])), we determine the coupling between massive KK fermion modes and the radion field as, $$\begin{aligned}
    &\lambda&^{(n)} = \sqrt{\frac{k}{24M^3}}e^{A(\pi r_c)} e^{\frac{\kappa}{2\sqrt{3}}\big(<\Phi>+\xi(\pi r_c)\big)}
    \exp{\bigg[-\frac{p}{4k}e^{(2+\nu)k\pi r_c}\bigg]}\nonumber\\
   &\exp&{\big[-\frac{1}{4k}\sqrt{q^2e^{4(2+\nu)k\pi r_c}+k^2+4m_{n}^2e^{2k\pi r_c}+4m_5^2}\big]}*
   \bigg[HypergeometricU\bigg[\frac{5}{8}\nonumber\\
   &+&\frac{qe^{2(2+\nu)k\pi r_c}}{8k} + \frac{\sqrt{q^2e^{4(2+\nu)k\pi r_c}+k^2+4m_{n}^2e^{2k\pi r_c}+4m_5^2}}{8k},\nonumber\\
   &1&+ \frac{\sqrt{q^2e^{4(2+\nu)k\pi r_c}+k^2+4m_{n}^2e^{2k\pi r_c}+4m_5^2}}{4k},\frac{p}{4k}e^{(2+\nu)k\pi r_c}\bigg]\bigg]^2
  \label{coupling}
 \end{aligned}$$ where $\lambda^{(n)}$ is the coupling between $n$ th KK fermion mode and the radion field and $p$, $q$ are given in eqn.(\[p\]), eqn.(\[q\]) respectively. With increasing $\alpha$, the third argument of Hypergeometric function increases in the above expression and as a result $\lambda^{(n)}$ decreases.\
Eqn.(\[solution3\]) indicates the relation between $\zeta_{L,R}$ and $y$ for various values of $l$ from which one can find the dependence of localization for massive KK fermion modes on the backreaction parameter. From this, the behaviour of the first KK mode ($n=1$) is described in figure (\[plot massive mode\]).

![$\zeta^{(1)}_{L,R}(=\zeta^{(1)}_{L,R}(F(R))$ in the diagram) vs $y$ for $M=1$, $\Lambda=-1$ and $m_5=0.5$[]{data-label="plot massive mode"}](massive.eps){width="3.0in" height="2.0in"}

Figure (\[plot massive mode\]) clearly depicts that the wave function for first massive KK mode gets more and more localized near Planck brane with increasing value of higher curvature parameter. Consequently, the coupling parameter decreases near the visible 