  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 $m\alpha^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 m\alpha^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/R\ll m\alpha^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/R\sim m\alpha^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.
