We investigate quasi-Monte Carlo rules for the numerical integration of multivariate periodic functions from Besov spaces S r p,q B(T d ) with dominating mixed smoothness 1/ p < r < 2. We show that order 2 digital nets achieve the optimal rate of convergence N −r (log N ) (d−1)(1−1/q) . The logarithmic term does not depend on r and hence improves the known bound of Dick (SIAM J Numer Anal 45: 2007) for the special case of Sobolev spaces H r mix (T d ). Secondly, the rate of convergence is independent of the integrability p of the Besov space, which allows for sacrificing integrability in order to gain Besov regularity. Our method combines characterizations of periodic Besov spaces with dominating mixed smoothness via Faber bases with sharp estimates of Haar coefficients for the discrepancy function of order 2 digital nets. Moreover, we provide numerical computations which indicate that this bound also holds for the case r = 2.Mathematics Subject Classification 11K06 · 11J71 · 42C10 · 46E35 · 65C05 · 65D30 · 65D32 · 91G60