Theoretical investigations of a platinum-water interface using quantum-mechanics-molecular-mechanics based molecular dynamics simulations.

Pt-water interfaces have been of immense interest in the field of energy storage and conversion. Studying this interface using both experimental and theoretical tools is challenging. On the theoretical front, typically one uses classical molecular dynamics (MD) simulations to handle large system sizes or time scales while for a more accurate quantum mechanical description Born Oppenheimer MD (BOMD) is typically used. The latter is limited to smaller system sizes and time-scales. In this study using quantum-mechanics-molecular-mechanics (QMMM), we have performed atomistic MD simulations to have a microscopic understanding of the structure of the Pt-water interface using a system size that is much larger than that accessible when using BOMD simulations. In contrast to recent reports using BOMD simulations, our study reveals that the water molecules typically form two distinct layers above the Pt-surface before they form bulk like structures. Further, we also find that a significant fraction of the water molecules at the interface are pointed towards the surface thereby disrupting the H-bond network. Consistent with this observation, the layer resolved oxygen-oxygen radial distribution function for the water molecules belonging to the solvating water layer shows a high density liquid like behaviour even though the overall water behaves like a low density liquid. A charge transfer analysis reveals that this solvating water layer donates electrons to the Pt atoms in contact with it thereby resulting in the formation of an interface dipole that is pointing towards the surface. Our results suggest that, using QMMM-MD, on one hand it is possible to study more realistic models of solid-liquid interfaces that are inaccessible with BOMD, while on the other hand one also has access to information about such systems that are not obtained from conventional classical MD simulations.

Intermediate phase of the one dimensional half-filled Hubbard-Holstein model.

We present a numerical study of the Hubbard-Holstein model in one dimension at half filling, including finite-frequency quantum phonons. At half filling, the effects of the electron-phonon and electron-electron interactions compete with the Holstein phonon coupling acting as an effective negative Hubbard on-site interaction U that promotes on-site electron pairs and a Peierls charge-density wave state. Most previous work on this model has assumed that only Peierls or Mott phases are possible at half filling. However, there has been speculation that a third metallic phase exists between the Peierls and Mott phases. We confirm the intermediate phase, and show that the Luttinger liquid correlation exponent K(rho) >1 in this region, indicating dominant superconducting pair correlations. We explore the full phase diagram as a function of Hubbard U, phonon coupling constant, and phonon frequency.