obtain Figure (\[plot massless1\]) between $\zeta_{L,R}$ and $y$ for various values of the higher curvature parameter $\alpha$. We focus into the region near the TeV brane (see Figure (\[plot massless1\])) to depict the localization properties of left and right chiral modes.\

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

Figure (\[plot massless1\]) clearly demonstrates that for $m_5=0$, the two chiral modes get more and more localized on TeV brane as the value of the parameter $\alpha$ increases. On the other hand, for small values of $\alpha$, the fermions are clearly localized deep inside the bulk spacetime. Thus without any bulk mass term, the fermions can be localized at different regions inside the bulk by adjusting the value of higher curvature parameter.\
From eqn.(\[solution2\]), we obtain the plots (Figure (\[plot massless2\]) and Figure (\[plot massless3\])) of left and right chiral massless modes for various values of the parameter $\alpha$ in presence of non zero bulk fermionic mass.\

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

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

Figure (\[plot massless2\]) and Figure (\[plot massless3\]) reveal that once again as the 