The hydrogen plasma is studied in the very high density (atomic and metallic) regime by extensive ab initio molecular dynamics simulations. Protons are treated classically, and electrons in the Born-Oppenheimer framework, within the local density approximation to density functional theory. Densities and temperatures studied fall within the strong coupling regime of the protons. We address the question of the validity of linear screening, and we find it yields a reasonably good description up to ${\mathit{r}}_{\mathit{s}}$\ensuremath{\approxeq}0.5, but already too crude for ${\mathit{r}}_{\mathit{s}}$=1 (with ${\mathit{r}}_{\mathit{s}}$=(3/4\ensuremath{\pi}\ensuremath{\rho}${)}^{1/3}$ the ion sphere radius). These values are typical of Jovian planets interiors. Finite-size and Brillouin zone sampling effects in metallic systems are studied and shown to be very delicate also in the fluid (liquid metal) phase. We analyze the low-temperature phase diagram and the melting transition. A remarkably fast decrease of the melting temperature with decreasing density is found, up to a point when it becomes comparable to the Fermi temperature of the protons. The possible vicinity of a triple point bcc-hcp(fcc)-liquid is discussed in the region of ${\mathit{r}}_{\mathit{s}}$\ensuremath{\approxeq}1.1 and T\ensuremath{\approxeq}100-200 K. The fluid phase is studied in detail for several temperatures. The structure of the fluid is found to be reminiscent of the underlying bcc (solid) phase. Proton-electron correlations show a weak temperature dependence, and proton-proton correlations exhibit a well-defined first coordination shell, thus characterizing fluid H in this regime as an atomic liquid. Diffusion coefficients are computed and compared to the values for the one-component plasma. Vibrational densities of states (VDOS) show a plasmon renormalization due to electron screening, and the presence of a plasmon-coupled single-particle mode up to very high temperatures. Collective modes are studied through dynamical structure factors. In close relationship with the VDOS, the simulations reveal the remarkable persistence of a weakly damped high-frequency ion-acoustic mode, even under conditions of strong electron screening. The possibility of using this observation as a diagnostic for the plasma phase transition to the fluid molecular phase at lower densities is discussed. \textcopyright{} 1996 The American Physical Society.