Density functional molecular dynamics method is applied to hydrogen plasma in the domain of liquid metallic hydrogen. Tentative results for the proton-proton pair distribution function and the electron-proton pair distribution function are obtained. It is shown that with the increase of the parameter γ(s), we have increasingly strong screening of proton charge by electrons and the decrease of electron density in the domain between protons.