Abstract. We introduce a new Petrov-Galerkin multiscale method for the numerical approximation of the Helmholtz equation with large wave number κ in bounded domains in R d . The discrete trial and test spaces are generated from standard mesh-based finite elements by local subscale correction in the spirit of numerical homogenization. The precomputation of the correction involves the solution of coercive cell problems on localized subdomains of size ℓH; H being the mesh size and ℓ being the oversampling parameter. If the mesh size and the oversampling parameter are such that Hκ and log(κ)/ℓ fall below some generic constants and if the cell problems are solved sufficiently accurate on some finer scale of discretization, then the method is stable and its error is proportional to H; pollution effects are eliminated in this regime.