The question of uniqueness of harmonic functions in a domain Ω ⊂ R d \Omega \subset \mathbb {R}^d with boundary ∂ Ω \partial \Omega , satisfying Dirichlet boundary conditions and with normal derivatives vanishing on a subset of the boundary ω \omega of positive d − 1 d-1 dimensional measure has attracted a lot of attention (see e.g. V. Adolfsson and L. Escauriaza [Comm. Pure Appl. Math. 50 (1997), pp. 935–969]; F.-H. Lin [Comm. Pure Appl. Math. 44 (1991), pp. 287–308]; X. Tolsa [Comm. Pure Appl. Math. 76 (2023), pp. 305–336]). It is essentially a consequence of Carleman estimates when ω \omega contains an open subset of the boundary (uniqueness for second order elliptic operators across an hypersurface). The main open questions (about uniqueness) concern now Lipschitz domains and variable coefficients. Here, using results by Logunov and Malinnikova [Nodal sets of Laplace eigenfunctions: estimates of the Hausdorff measure in dimensions two and three, Birkhäuser/Springer, Cham, 2018, pp. 333-344; Proceedings of the International Congress of Mathematicians–Rio de Janeiro 2018, World Scientific Publishing, Hackensack, NJ, 2018, pp. 2391–2411], we consider the simpler case of W 2 , ∞ W^{2, \infty } domains but prove quantitative uniqueness both for Dirichlet and Neumann boundary conditions. As an application, we deduce quantitative estimates for the Dirichlet and Neumann Laplace eigenfunctions on a W 2 , ∞ W^{2, \infty } domain with boundary.