So far the basic physical properties of matter forming a thin accretion disc in the static and spherically symmetric space–time metric of the vacuum f(R) modified gravity models (Pun et al. in Phys. Rev. D 78:024043, 2008) and building radiative models of thin accretion disks for both Schwarzschild and Kerr black holes in f(R) gravity (Perez et al. in Astron. Astrophys. 551:4, 2013) were addressed properly. Also von Zeipel surfaces and convective instabilities in f(R) -Schwarzschild(Kerr) background have been investigated recently (Alipour et al. in Mon. Not. R. Astron. Soc. 454:1992, 2015). In this streamline, here we study the effects of radial and angular pressure gradients on thick accretion disks in Schwarzschild geometries in a constant curvature f(R) modified gravity. Since thick accretion disks have high accretion rate, we study configuration and structure of thick disks by focusing on the effect of pressure gradient on formation of the disks. We clarify our study by assuming two types of equation of state: polytropic and Clapeyron equation of states.
