Optimal boundary regularity and Green function estimates for nonlocal equations in divergence form

In this article we prove for the first time the Cˢ boundary regularity for solutions to nonlocal elliptic equations with H\"older continuous coefficients in divergence form in C^1, domains. So far, it was only known that solutions are H\"older continuous up to the boundary, and establishing their optimal regularity has remained an open problem in the field. Our proof is based on a delicate higher order Campanato-type iteration at the boundary, which we develop in the context of nonlocal equations and which is quite different from the local theory. As an application of our results, we establish sharp two-sided Green functions estimates in C^1, domains for the same class of operators. Previously, this was only known under additional structural assumptions on the coefficients and in more regular domains.