We obtain regularity results in weighted Sobolev spaces for the solution of the obstacle problem for the integral fractional Laplacian (−∆) s . The weight is a power of the distance to the boundary that accounts for the singular boundary behavior of the solution for any 0 < s < 1. These bounds then serve us as a guide in the design and analysis of an optimal finite element scheme over graded meshes.definitions of the operator (−∆) s , motivated by applications, here we choose the so-called integral one; that is, for a sufficiently smooth function v : R n → R we setOur choice of definition is justified by the fact that, unlike the regional or the spectral ones, the integral fractional Laplacian of order s is the infinitesimal generator of a 2s-stable Lévy process. These processes have been widely employed for modeling market fluctuations, both for risk management and option pricing purposes. It is in this context that, as mentioned above, the fractional obstacle problem arises as a pricing model for American options. More precisely, if u represents the rational price of a perpetual American option, modeling the assets prices by a Lévy process X t and denoting by χ the payoff function, then u solves (1.3). We refer the reader to [16] for an overview of the use of jump processes in financial modeling.Taking into account their applications in finance, it is not surprising that numerical schemes for integrodifferential inequalities have been proposed and analyzed in the literature; we refer the reader to [26] for a survey on these methods. These applications aim to approximate the price of a number of assets; therefore, the consideration of a logarithmic price leads to problems posed in the whole space R n . For the numerical solution, it is usual to perform computations on a sufficiently large tensor-product domain. Among the schemes based on Galerkin discretizations, reference [46] utilizes piecewise linear Lagrangian finite elements, while [29] proposes the use of wavelet bases in space. As for approximations of variational inequalities involving integral operators on arbitrary bounded domains, an a posteriori error analysis is performed in [36].Since the seminal work of Silvestre [45], the fractional obstacle problem started to draw the attention of the mathematical community. Using potential theoretic methods, reference [45] shows that if the obstacle is of class C 1,s , then the solution to the fractional obstacle problem is of class C 1,α for all α ∈ (0, s); optimal C 1,s regularity of solutions was derived assuming convexity of the contact set. The pursuit of the optimal regularity of solutions without a convexity hypothesis, in turn, motivated the celebrated extension by Caffarelli and Silvestre [11] for the fractional Laplacian in R n . Using this extension technique, Caffarelli, Salsa and Silvestre proved, in [10], the optimal regularity of solutions (cf. Proposition 3.4 below). It is important to notice, however, that this is only an interior regularity result. Nothing is said about the boundary behavior o...