We consider a class of nonhomogeneous elliptic equations with fractional Laplacian and nonlinear gradient terms, namely (−Δ)α2u=V(x)u+g(u,∇u)+f in ℝn, where 0<α<n, g is the nonlinearity, V the potential and f is a forcing term. Some examples of nonlinearities dealt with are u|u|ρ−1, |∇u|ρ and |u|ρ1|∇u|ρ2, covering large values of ρ,ρ1,ρ2, and particularly variational supercritical powers for u and super-α ones for |∇u| (superquadratic if α=2). Moreover, we are able to consider some exponential growths, g belonging to certain classes of power series, or g satisfying some conditions in the Lipschitz spirit. We obtain results on existence, uniqueness, symmetry, and other qualitative properties in a new framework, namely modulation-type spaces based on Lorentz spaces. For that, we need to develop properties and estimates in those spaces such as complex interpolation, Hölder-type inequality, estimates for product, convolution and Riesz potential operators, among others. In order to handle the nonlinearity, other ingredients are estimates for composition operators in our setting.