  Interfacial behaviours in multiphase systems containing H2 are crucial to
underground H2 storage but are not well understood. Molecular dynamics
simulations were conducted to study interfacial properties of the H2O+H2 and
H2O+H2+silica/kerogen systems over a wide range of temperatures (298 - 523 K)
and pressures (1 - 160 MPa). The combination of the H2 model with the INTERFACE
force field and TIP4P/2005 H2O model can accurately predict the interfacial
tensions (IFTs) from the experiment. The IFTs from simulations are also in good
agreement with those from the density gradient theory coupled to the PC-SAFT
equation of state. Generally, the IFTs decrease with pressure and temperature.
However, at relatively high temperatures and pressures, the IFTs increase with
pressure. The opposite pressure effect on IFTs can be explained by the
inversion of the sign of the relative adsorption of H2. The enrichment of H2 in
the interfacial regions was observed in density profiles. Meanwhile, the
behaviours of contact angles (CAs) in the H2O+H2+silica system are noticeably
different from those in the H2O+H2+kerogen system. The H2O CAs for the
H2O+H2+silica and H2O+H2+kerogen systems increase with pressure and decrease
with temperature. However, the effect of temperature and pressure on these CAs
is less pronounced for the H2O+H2+silica system at low temperatures. The
behaviours of CAs were understood based on the variations of IFTs in the H2O+H2
system (fluid-fluid interaction) and adhesion tensions (fluid-solid
interaction). Furthermore, the analysis of the atomic density profiles shows
that the presence of H2 in between the H2O droplet and the silica/kerogen
surface is almost negligible. Nevertheless, the adsorption of H2O on the silica
surface outside the H2O droplet is strong, while less H2O adsorption is seen on
the kerogen surface.
