[@faraoni; @felice; @barrow; @marino; @bahamonde; @catena; @sumanta_higher; @sumanta_lhc]. In such a scenario, several models have been studied to explore their phenomenological implications. Some of these are formulated by placing the standard model fields inside the bulk . For these models the localization property of a bulk fermion field on brane [@bajc; @smyth; @tp_fermion; @chang; @koley; @ssg3; @grossman; @jm] has been a subject of great interest to explore the chiral nature of massless fermion and also the hierarchial masses of fermions among different generations [@grossman]. In this context, it is observed that the overlap of the bulk fermion wave function on our visible brane plays crucial role in determining the effective radion-fermion coupling which in turn determines phenomenology of the radion with brane matter fields.\
In view of above, it is worthwhile to address the role of higher curvature terms on fermion localization. We aim to address this in the present work.\
Our paper is organized as follows : The mapping between higher curvature and scalar degrees of freedom is discussed in section II. The description of warped geometry in F(R) model and its solutions are described in section III and section IV. Section V addresses the localization property of bulk 