We study the Dirichlet problem for the stationary Schrödinger fractional Laplacian equation (−∆) s u + V u = f posed in bounded domain Ω ⊂ R n with zero outside conditions. We consider general nonnegative potentials V ∈ L 1 loc (Ω) and prove well-posedness of very weak solutions when the data are chosen in an optimal class of weighted integrable functions f . Important properties of the solutions, such as its boundary behaviour, are derived. The case of super singular potentials that blow up near the boundary is given special consideration since it leads to so-called flat solutions. We comment on related literature.