Van der Waals interactions between two neutral but polarizable systems at a separation $R$ much larger than the typical size of the systems are at the core of a broad sweep of contemporary problems in settings ranging from atomic, molecular and condensed matter physics to strong interactions and gravity. We reexamine the dispersive van der Waals interactions between two hydrogen atoms. The novelty of the analysis resides in the usage of nonrelativistic EFTs of QED. In this framework, the van der Waals potential acquires the meaning of a matching coefficient in an EFT suited to describe the low energy dynamics of an atom pair. It may be computed systematically as a series in $R$ times some typical atomic scale and in the fine structure constant $alpha$. The van der Waals potential gets short range contributions and radiative corrections, which we compute in dimensional regularization and renormalize here for the first time. Results are given in $d$ spacetime dimensions. One can distinguish among different regimes depending on the relative size between $1/R$ and the typical atomic bound state energy $malpha^2$. Each regime is characterized by a specific hierarchy of scales and a corresponding tower of EFTs. The short distance regime is characterized by $1/R gg malpha^2$ and the LO van der Waals potential is the London potential. We compute also NNNLO corrections. In the long distance regime we have $1/Rll malpha^2$. In this regime, the van der Waals potential contains contact terms, which are parametrically larger than the Casimir-Polder potential that describes the potential at large distances. In the EFT the Casimir-Polder potential counts as a NNNLO effect. In the intermediate distance regime, $1/Rsim malpha^2$, a significantly more complex potential is obtained which we compare with the two previous limiting cases. We conclude commenting on the hadronic van der Waals case.