Dispersive interactions in water bilayers at metallic surfaces: A comparison of the PBE and RPBE functional including semiempirical dispersion corrections

Orientational Water Bonding on Pt(111): Beyond the Frontier Orbital Principle

For decades, our understanding of water-metal bonding has been dominated by the frontier orbital principle in which globally stable water-metal interactions are ruled by HOMO interacting with metal surfaces. Using density functional theory calculations, herein, we have revealed that the frontier orbital principle cannot be applied to metastable water bonding on Pt(111), where the decisive role of HOMO is replaced by HOMO-1 in terms of the greatest orbital shifts and depopulations as the two different bonding indicators. Unlike the stable water configuration in which both HOMO-1 and HOMO prefer to overlap with metal states through σ-like orbital interactions, metastable configurations exhibit delicate competition or balance between σ-like and π-like orbital interactions exerted by HOMO-1 and HOMO, respectively. These findings have significantly deepened our understanding of orbital roles in water-metal bonding interactions and bridged the gap between theoretical understanding of electrified waters at electrochemical interfaces and water science on metal surfaces.

A Review of Trends in Corrosion-Resistant Structural Steels Research-From Theoretical Simulation to Data-Driven Directions

This paper provides a review of models commonly used over the years in the study of microscopic models of material corrosion mechanisms, data mining methods and the corrosion-resistant performance control of structural steels. The virtual process of material corrosion is combined with experimental data to reflect the microscopic mechanism of material corrosion from a nano-scale to macro-scale, respectively. Data mining methods focus on predicting and modeling the corrosion rate and corrosion life of materials. Data-driven control of the corrosion resistance of structural steels is achieved through micro-alloying and organization structure control technology. Corrosion modeling has been used to assess the effects of alloying elements, grain size and organization purity on corrosion resistance, and to determine the contents of alloying elements.

Electrochemical hydrogen evolution on Pt-based catalysts from a theoretical perspective

Hydrogen evolution reaction (HER) by splitting water is a key technology toward a clean energy society, where Pt-based catalysts were long known to have the highest activity under acidic electrochemical conditions but suffer from high cost and poor stability. Here, we overview the current status of Pt-catalyzed HER from a theoretical perspective, focusing on the methodology development of electrochemistry simulation, catalytic mechanism, and catalyst stability. Recent developments in theoretical methods for studying electrochemistry are introduced, elaborating on how they describe solid-liquid interface reactions under electrochemical potentials. The HER mechanism, the reaction kinetics, and the reaction sites on Pt are then summarized, which provides an atomic-level picture of Pt catalyst surface dynamics under reaction conditions. Finally, state-of-the-art experimental solutions to improve catalyst stability are also introduced, which illustrates the significance of fundamental understandings in the new catalyst design.

Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes

Despite its simple molecular formula, obtaining an accurate in silico description of water is far from straightforward. Many of its very peculiar properties are quite elusive, and in particular, obtaining good estimations of the diffusion coefficients of the solvated proton and hydroxide at a reasonable computational cost has been an unsolved challenge until now. Here, I present extensive results of several unusually long ab initio molecular dynamics (MD) simulations employing different combinations of the Born-Oppenheimer and second-generation Car-Parrinello MD propagation methods with different ensembles (NVE and NVT) and thermostats, which show that these methods together with the RPBE-D3 functional provide a very accurate estimation of the diffusion coefficients of the solvated H₃O⁺ and OH^- ions, together with an extremely accurate description of several properties of neutral water (such as the structure of the liquid and its diffusion and shear viscosity coefficients). In addition, I show that the estimations of D_H3O⁺ and D_OH^- depend dramatically on the simulation length, being necessary to reach timescales in the order of hundreds of picoseconds to obtain reliable results.

Nanosecond solvation dynamics of the hematite/liquid water interface at hybrid DFT accuracy using committee neural network potentials

Metal oxide/water interfaces play an important role in biology, catalysis, energy storage and photocatalytic water splitting. The atomistic structure at these interfaces is often difficult to characterize by experimental techniques, whilst results from ab initio molecular dynamics simulations tend to be uncertain due to the limited length and time scales accessible. In this work, we train a committee neural network potential to simulate the hematite/water interface at the hybrid DFT level of theory to reach the nanosecond timescale and systems containing more than 3000 atoms. The NNP enables us to converge dynamical properties, not possible with brute-force ab initio molecular dynamics. Our simulations uncover a rich solvation dynamics at the hematite/water interface spanning three different time scales: picosecond H-bond dynamics between surface hydroxyls and the first water layer, in-plane/out-of-plane tilt motion of surface hydroxyls on the 10 ps time scale, and diffusion of water molecules from the oxide surface characterized by a mean residence lifetime of about 60 ps. Calculation of vibrational spectra confirm that H-bonds between surface hydroxyls and first layer water molecules are stronger than H-bonds in bulk water. Our study showcases how state of the art machine learning approaches can routinely be utilized to explore the structural dynamics at transition metal oxide interfaces with complex electronic structure. It foreshadows that c-NNPs are a promising tool to tackle the sampling problem in ab initio electrochemistry with explicit solvent molecules.

Structures of liquid and aqueous water isotopologues at ambient temperature from ab initio path integral simulations

The heavy hydrogen isotopes D and T are found in trace amounts in water, and when their concentration increases they can play an intricate role in modulating the physical properties of the liquid. We present an analysis of the microscopic structures of ambient light water (H₂O(l)), heavy water (D₂O(l)), T₂O(l), HDO(aq) and HTO(aq) studied by ab initio path integral molecular dynamics (PIMD). Unlike previous ab initio PIMD investigations of H₂O(l) and D₂O(l) [Chen et al., Phys. Rev. Lett., 2003, 91, 215503] [Machida et al., J. Chem. Phys., 2017, 148, 102324] we find that D₂O(l) is more structured than H₂O(l), as is predicted by the experiment. The agreement between the experiment and our simulation for H₂O(l) and D₂O(l) allows us to accurately predict the intra- and intermolecular structures of T₂O(l) HDO(aq) and HTO(aq). T₂O(l) is found to have a similar intermolecular structure to that of D₂O(l), while the intramolecular structure is more compact, giving rise to a smaller dipole moment than those of H₂O(l) and D₂O(l). For the mixed isotope species, HDO(aq) and HTO(aq), we find smaller dipole moments and fewer hydrogen bonds when compared with the pure species H₂O and D₂O. We can attribute this effect to the relative compactness of the mixed isotope species, which results in a lower dipole moment than that of the pure species.

Ab Initio Simulations of Water/Metal Interfaces

Structures and processes at water/metal interfaces play an important technological role in electrochemical energy conversion and storage, photoconversion, sensors, and corrosion, just to name a few. However, they are also of fundamental significance as a model system for the study of solid-liquid interfaces, which requires combining concepts from the chemistry and physics of crystalline materials and liquids. Particularly interesting is the fact that the water-water and water-metal interactions are of similar strength so that the structures at water/metal interfaces result from a competition between these comparable interactions. Because water is a polar molecule and water and metal surfaces are both polarizable, explicit consideration of the electronic degrees of freedom at water/metal interfaces is mandatory. In principle, ab initio molecular dynamics simulations are thus the method of choice to model water/metal interfaces, but they are computationally still rather demanding. Here, ab initio simulations of water/metal interfaces will be reviewed, starting from static systems such as the adsorption of single water molecules, water clusters, and icelike layers, followed by the properties of liquid water layers at metal surfaces. Technical issues such as the appropriate first-principles description of the water-water and water-metal interactions will be discussed, and electrochemical aspects will be addressed. Finally, more approximate but numerically less demanding approaches to treat water at metal surfaces from first-principles will be briefly discussed.

A first principle study of water adsorbed on flat and stepped silver surfaces

The structural, electronic and vibrational properties of a water layer on Ag(100) and Ag(511) have been studied by first principles calculations and ab initio molecular dynamics simulations. The most stable...

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

In situ Raman spectroscopy reveals the structure and dissociation of interfacial water

Understanding the structure and dynamic process of water at the solid-liquid interface is an extremely important topic in surface science, energy science and catalysis^1-3. As model catalysts, atomically flat single-crystal electrodes exhibit well-defined surface and electric field properties, and therefore may be used to elucidate the relationship between structure and electrocatalytic activity at the atomic level^4,5. Hence, studying interfacial water behaviour on single-crystal surfaces provides a framework for understanding electrocatalysis^6,7. However, interfacial water is notoriously difficult to probe owing to interference from bulk water and the complexity of interfacial environments⁸. Here, we use electrochemical, in situ Raman spectroscopic and computational techniques to investigate the interfacial water on atomically flat Pd single-crystal surfaces. Direct spectral evidence reveals that interfacial water consists of hydrogen-bonded and hydrated Na⁺ ion water. At hydrogen evolution reaction (HER) potentials, dynamic changes in the structure of interfacial water were observed from a random distribution to an ordered structure due to bias potential and Na⁺ ion cooperation. Structurally ordered interfacial water facilitated high-efficiency electron transfer across the interface, resulting in higher HER rates. The electrolytes and electrode surface effects on interfacial water were also probed and found to affect water structure. Therefore, through local cation tuning strategies, we anticipate that these results may be generalized to enable ordered interfacial water to improve electrocatalytic reaction rates.

Ab initio study of nuclear quantum effects on sub- and supercritical water

The structures of water in the ambient, subcritical, and supercritical conditions at various densities were studied systematically by ab initio path integral molecular dynamics simulations. It was found that the nuclear quantum effects (NQEs) have a significant impact on the structure of hydrogen bonds in close contact, not only in the ambient condition but also in the sub- and supercritical conditions. The NQEs on the structure beyond the hydrogen bond contact are important in ambient water, but not much for water in the sub- and supercritical conditions. The NQEs are furthermore important for determining the number of hydrogen bonds in the ambient conditions, and this role is, however, diminished in the sub- and supercritical conditions. The NQEs do, nevertheless, show their importance in determining the intramolecular structure of water and the close contact structures of the hydrogen bonds, even at sub- and supercritical conditions. Using the RPBE-D3 functional, the computed radial distribution functions for ambient water are in excellent agreement with experimental data, upgrading our previous results using the BLYP-D2 functional [Machida et al., J. Chem. Phys. 148, 102324 (2018)]. The computed radial distribution functions for water in the sub- and supercritical conditions were carefully compared with experiment. In particular, we found that the first peak in hydrogen pair distribution functions matches only when the NQEs are taken into account.

Static and dynamic water structures at interfaces: A case study with focus on Pt(111)

An accurate atomistic treatment of aqueous solid-liquid interfaces necessitates the explicit description of interfacial water ideally via ab initio molecular dynamics simulations. Many applications, however, still rely on static interfacial water models, e.g., for the computation of (electro)chemical reaction barriers and focus on a single, prototypical structure. In this work, we systematically study the relation between density functional theory-derived static and dynamic interfacial water models with specific focus on the water-Pt(111) interface. We first introduce a general construction protocol for static 2D water layers on any substrate, which we apply to the low index surfaces of Pt. Subsequently, we compare these with structures from a broad selection of reference works based on the Smooth Overlap of Atomic Positions descriptor. The analysis reveals some structural overlap between static and dynamic water ensembles; however, static structures tend to overemphasize the in-plane hydrogen bonding network. This feature is especially pronounced for the widely used low-temperature hexagonal ice-like structure. In addition, a complex relation between structure, work function, and adsorption energy is observed, which suggests that the concentration on single, static water models might introduce systematic biases that are likely reduced by averaging over consistently created structural ensembles, as introduced here.

Properties of the Pt(111)/electrolyte electrochemical interface studied with a hybrid DFT-solvation approach

Self-consistent modeling of the interface between solid metal electrode and liquid electrolyte is a crucial challenge in computational electrochemistry. In this contribution, we adopt the effective screening medium reference interaction site method (ESM-RISM) to study the charged interface between a Pt(111) surface that is partially covered with chemisorbed oxygen and an aqueous acidic electrolyte. This method proves to be well suited to describe the chemisorption and charging state of the interface at controlled electrode potential. We present an in-depth assessment of the ESM-RISM parameterization and of the importance of computing near-surface water molecules explicitly at the quantum mechanical level. We found that ESM-RISM is able to reproduce some key interface properties, including the peculiar, non-monotonic charging relation of the Pt(111)/electrolyte interface. The comparison with independent theoretical models and explicit simulations of the interface reveals strengths and limitations of ESM-RISM for modeling electrochemical interfaces.

Structure and Energetics of Dye-Sensitized NiO Interfaces in Water from Ab Initio MD and Large-Scale GW Calculations

The energy-level alignment across solvated molecule/semiconductor interfaces is a crucial property for the correct functioning of dye-sensitized photoelectrodes, where, following the absorption of solar light, a cascade of interfacial hole/electron transfer processes has to efficiently take place. In light of the difficulty of performing X-ray photoelectron spectroscopy measurements at the molecule/solvent/metal-oxide interface, being able to accurately predict the level alignment by first-principles calculations on realistic structural models would represent an important step toward the optimization of the device. In this respect, dye/NiO surfaces, employed in p-type dye-sensitized solar cells, are undoubtedly challenging for ab initio methods and, also for this reason, much less investigated than the n-type dye/TiO₂ counterpart. Here, we consider the C343-sensitized NiO surface in water and combine ab initio molecular dynamics (AIMD) simulations with GW (G₀W₀) calculations, performed along the MD trajectory to reliably describe the structure and energetics of the interface when explicit solvation and finite temperature effects are accounted for. We show that the differential perturbative correction on the NiO and molecule states obtained at the GW level is mandatory to recover the correct (physical) interfacial energetics, allowing hole transfer from the semiconductor valence band to the highest occupied molecular orbital (HOMO) of the dye. Moreover, the calculated average driving force quantitatively agrees with the experimental estimate.

Pairing double hybrid functionals with a tailored basis set for an accurate thermochemistry of hydrocarbons

A collection of five challenging datasets, including noncovalent interactions, reaction barriers and electronic rearrangements of medium-sized hydrocarbons, has been selected to verify the robustness of double-hybrid functionals used in conjunction with the small DH-SVPD basis set, especially developed for noncovalent interactions. The analysis is completed by other, more standard functionals, for a total of 17 models, including also empirical corrections for dispersion. The obtained results show that the chemical accuracy threshold, that is an error lower than 1.0 kcal mol⁻¹, can be obtained by pairing the nonempirical PBE-QIDH functional with the DH-SVPD basis set, as well as by other semi-empirical functionals, such as DSD-PBEP86, using larger basis sets and empirical corrections. More in general, a significant improvement can be obtained using the DH-SVPD basis set with DHs, without resorting to any empirical corrections. This choice leads to a fast computational protocol that, avoiding any empirical potential, remains on a fully quantum ground.

The pairing of the PBE-QIDH double-hybrid functional with a tailored split-valence basis set leads to a fast computational protocol for the accurate evaluation of hydrocarbon thermochemistry, without resorting to any empirical correction.

Aqueous TMAO solution under high hydrostatic pressure

Trimethylamine N-oxide (TMAO) is a well known osmolyte in nature, which is used by deep sea fish to stabilize proteins against High Hydrostatic Pressure (HHP). We present a combined ab initio molecular dynamics, force field molecular dynamics, and THz absorption study of TMAO in water up to 12 kbar to decipher its solvation properties upon extreme compression. On the hydrophilic oxygen side of TMAO, AIMD simulations at 1 bar and 10 kbar predict a change of the coordination number from a dominating TMAO·(H₂O)₃ complex at ambient conditions towards an increased population of a TMAO·(H₂O)₄ complex at HHP conditions. This increase of the TMAO-oxygen coordination number goes in line with a weakening of the local hydrogen bond network, spectroscopic shifts and intensity changes of the corresponding intermolecular THz bands. Using a pressure-dependent HHP force field, FFMD simulations predict a significant increase of hydrophobic hydration from 1 bar up to 4-5 kbar, which levels off at higher pressures up to 10 kbar. THz spectroscopic data reveal two important pressure regimes with spectroscopic inflection points of the dominant intermolecular modes: The first regime (1.5-2 kbar) is barely recognizable in the simulation data. However, it relates well with the observation that the apparent molar volume of solvated TMAO is nearly constant in the biologically relevant pressure range up to 1 kbar as found in the deepest habitats on Earth in the ocean. The second inflection point around 4-5 kbar is related to the amount of hydrophobic hydration as predicted by the FFMD simulations. In particular, the blueshift of the intramolecular CNC bending mode of TMAO at about 390 cm^-1 is the spectroscopic signature of increasingly pronounced pressure-induced changes in the solvation shell of TMAO. Thus, the CNC bend can serve as local pressure sensor in the multi-kbar pressure regime.

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C₂H₄) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousands of C₂H₄ spectra, computed from DFT, which regress surface binding-type and generalized coordination number, outperform those trained using CO and NO spectra. A python package, pDOS_overlap, for implementing the electron density-based analysis on any combination of adsorbates and materials, is also made available.

Supercritical Water is not Hydrogen Bonded

Thinking about water is inextricably linked to hydrogen bonds, which are highly directional in character and determine the unique structure of water, in particular its tetrahedral H-bond network. Here, we assess if this common connotation also holds for supercritical water. We employ extensive ab initio molecular dynamics simulations to systematically monitor the evolution of the H-bond network mode of water from room temperature, where it is the hallmark of its fluctuating three-dimensional network structure, to supercritical conditions. Our simulations reveal that the oscillation period required for H-bond vibrations to occur exceeds the lifetime of H-bonds in supercritical water by far. Instead, the corresponding low-frequency intermolecular vibrations of water pairs as seen in supercritical water are found to be well represented by isotropic van-der-Waals interactions only. Based on these findings, we conclude that water in its supercritical phase is not a H-bonded fluid.

Sulfate, Bisulfate, and Hydrogen Co-adsorption on Pt(111) and Au(111) in an Electrochemical Environment

The co-adsorption of sulfate, bisulfate and hydrogen on Pt(111) and Au(111) electrodes was studied based on periodic density functional calculations with the aqueous electrolyte represented by both explicit and implicit solvent models. The influence of the electrochemical control parameters such as the electrode potential and pH was taken into account in a grand-canonical approach. Thus, phase diagrams of the stable coadsorption phases as a function of the electrochemical potential and Pourbaix diagrams have been derived which well reproduce experimental findings. We demonstrate that it is necessary to include explicit water molecules in order to determine the stable adsorbate phases as the (bi)sulfate adsorbates rows become significantly stabilized by bridging water molecules.

Origins of fast diffusion of water dimers on surfaces

The diffusion of water molecules and clusters across the surfaces of materials is important to a wide range of processes. Interestingly, experiments have shown that on certain substrates, water dimers can diffuse more rapidly than water monomers. Whilst explanations for anomalously fast diffusion have been presented for specific systems, the general underlying physical principles are not yet established. We investigate this through a systematic ab initio study of water monomer and dimer diffusion on a range of surfaces. Calculations reveal different mechanisms for fast water dimer diffusion, which is found to be more widespread than previously anticipated. The key factors affecting diffusion are the balance of water-water versus water-surface bonding and the ease with which hydrogen-bond exchange can occur (either through a classical over-the-barrier process or through quantum-mechanical tunnelling). We anticipate that the insights gained will be useful for understanding future experiments on the diffusion and clustering of hydrogen-bonded adsorbates.

The experimental observation that water dimers diffuse more rapidly than monomers across materials’ surfaces is yet to be clarified. Here the authors show by ab initio calculations classical and quantum mechanical mechanisms for faster water dimer diffusion on a broad range of metal and non-metal surfaces.

Infrared spectroscopy data- and physics-driven machine learning for characterizing surface microstructure of complex materials

There is a need to characterize complex materials and their dynamics under reaction conditions to accelerate materials design. Adsorbate vibrational excitations are selective to adsorbate/surface interactions and infrared (IR) spectra associated with activating adsorbate vibrational modes are accurate, capture details of most modes, and can be obtained operando. Current interpretation depends on heuristic peak assignments for simple spectra, precluding the possibility of obtaining detailed structural information. Here, we combine data-based approaches with chemistry-dependent problem formulation to develop physics-driven surrogate models that generate synthetic IR spectra from first-principles calculations. Using synthetic IR spectra of carbon monoxide on platinum, we implement multinomial regression via neural network ensembles to learn probability distributions functions (pdfs) that describe adsorption sites and quantify uncertainty. We use these pdfs to infer detailed surface microstructure from experimental spectra and extend this methodology to other systems as a first step towards characterizing complex interfaces and closing the materials gap.

A theoretical study on the binding and electrolytic splitting of hydrogen fluoride on Ni(111) and Ni(211)

Anhydrous hydrogen fluoride (HF) is electrochemically activated on Ni electrodes in the Simons process. We study the electrode/electrolyte interface using DFT methods and find the splitting of HF to be exothermic even at low cell potentials.

Assessing the properties of supercritical water in terms of structural dynamics and electronic polarization effects

Evolution of water's structural dynamics from ambient liquid to supercritical dense liquid-like and dilute gas-like conditions.

Water structures on a Pt(111) electrode from ab initio molecular dynamic simulations for a variety of electrochemical conditions

A structural analysis of solvating water layers on a Pt(111) electrode has been performed based on extensive ab initio molecular dynamics simulations. We have emulated different electrochemical conditions by varying the concentration of hydrogen ions in the water layers, which effectively corresponds to a variation in the electrode potential. We present a detailed analysis of the arrangement and orientation of the water molecules and also address their mobility in the solvation layer.

Structure and thermodynamics of aqueous urea solutions from ambient to kilobar pressures: From thermodynamic modeling, experiments, and first principles simulations to an accurate force field description

Molecular simulations based on classical force fields are a powerful method for shedding light on the complex behavior of biomolecules in solution. When cosolutes are present in addition to water and biomolecules, subtle balances of weak intermolecular forces have to be accounted for. This imposes high demands on the quality of the underlying force fields, and therefore force field development for small cosolutes is still an active field. Here, we present the development of a new urea force field from studies of urea solutions at ambient and elevated hydrostatic pressures based on a combination of experimental and theoretical approaches. Experimental densities and solvation shell properties from ab initio molecular dynamics simulations at ambient conditions served as the target properties for the force field optimization. Since urea is present in many marine life forms, elevated hydrostatic pressure was rigorously addressed: densities at high pressure were measured by vibrating tube densitometry up to 500 bar and by X-ray absorption up to 5 kbar. Densities were determined by the perturbed-chain statistical associating fluid theory equation of state. Solvation properties were determined by embedded cluster integral equation theory and ab initio molecular dynamics. Our new force field is able to capture the properties of urea solutions at high pressures without further high-pressure adaption, unlike trimethylamine-N-oxide, for which a high-pressure adaption is necessary.

Initial Steps in Forming the Electrode-Electrolyte Interface: H2O Adsorption and Complex Formation on the Ag(111) Surface from Combining Quantum Mechanics Calculations and Ambient Pressure X-ray Photoelectron Spectroscopy

The interaction of water with metal surfaces is at the heart of electrocatalysis. But there remain enormous uncertainties about the atomistic interactions at the electrode-electrolyte interface (EEI). As the first step toward an understanding of the EEI, we report here the details of the initial steps of H₂O adsorption and complex formation on a Ag(111) surface, based on coupling quantum mechanics (QM) and ambient-pressure X-ray photoelectron spectroscopy (APXPS) experiments. We find a close and direct comparison between simulation and experiment, validated under various isotherm and isobar conditions. We identify five observable oxygen-containing species whose concentrations depend sensitively on temperature and pressure: chemisorbed O* and OH*, H₂O* stabilized by hydrogen bond interactions with OH* or O*, and multilayer H₂O*. We identify the species experimentally by their O 1s core-level shift that we calculate with QM along with the structures and free energies as a function of temperature and pressure. This leads to a chemical reaction network (CRN) that we use to predict the time evolution of their concentrations over a wide range of temperature (298-798 K) and pressure conditions (10^-6-1 Torr), which agree well with the populations determined from APXPS. This multistep simulation CRN protocol should be useful for other heterogeneous catalytic systems.

The crucial role of electron transfer from interfacial molecules in the negative potential shift of Au electrode immersed in ionic liquids

Potential of zero charge (PZC) is essential in electrochemistry to understand physical and chemical phenomena at the interface between an electrode and a solution. A negative potential shift from the work function to the PZC has been experimentally observed in a metal/ionic liquid (IL) system, but the mechanism remains unclear and controversial. In this paper we provide valuable insight into the mechanism on the potential shift in the Au/IL (1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)amide: [BMIM][TFSA]) system using a computational approach combining classical molecular dynamics simulations and first-principles calculations. By separately estimating some contributions to the potential shift, the shift is calculated in an easy-to-understand manner. The resultant PZC is shown to be in good agreement with the experimental one. Among the contributions, the electron redistribution at the Au/IL interface is found to provide the largest negative potential change. This indicates that the redistribution plays a crucial role in determining the potential shift of the Au electrode immersed in the IL. Detailed analyses suggest that the redistribution corresponds to the electron transfer not only from the anionic TFSA but also from the cationic BMIM molecules to the Au electrode surface. This unique observation is understood to originate from the interfacial structure where the IL molecules are in very close proximity to the electrode surface via dispersion interaction.

Structure formation and surface chemistry of ionic liquids on model electrode surfaces-Model studies for the electrode | electrolyte interface in Li-ion batteries

Ionic liquids (ILs) are considered as attractive electrolyte solvents in modern battery concepts such as Li-ion batteries. Here we present a comprehensive review of the results of previous model studies on the interaction of the battery relevant IL 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide ([BMP]⁺[TFSI]^-) with a series of structurally and chemically well-defined model electrode surfaces, which are increasingly complex and relevant for battery applications [Ag(111), Au(111), Cu(111), pristine and lithiated highly oriented pyrolytic graphite (HOPG), and rutile TiO₂(110)]. Combining surface science techniques such as high resolution scanning tunneling microscopy and X-ray photoelectron spectroscopy for characterizing surface structure and chemical composition in deposited (sub-)monolayer adlayers with dispersion corrected density functional theory based calculations, this work aims at a molecular scale understanding of the fundamental processes at the electrode | electrolyte interface, which are crucial for the development of the so-called solid electrolyte interphase (SEI) layer in batteries. Performed under idealized conditions, in an ultrahigh vacuum environment, these model studies provide detailed insights on the structure formation in the adlayer, the substrate-adsorbate and adsorbate-adsorbate interactions responsible for this, and the tendency for chemically induced decomposition of the IL. To mimic the situation in an electrolyte, we also investigated the interaction of adsorbed IL (sub-)monolayers with coadsorbed lithium. Even at 80 K, postdeposited Li is found to react with the IL, leading to decomposition products such as LiF, Li₃N, Li₂S, Li_xSO_y, and Li₂O. In the absence of a [BMP]⁺[TFSI]^- adlayer, it tends to adsorb, dissolve, or intercalate into the substrate (metals, HOPG) or to react with the substrate (TiO₂) above a critical temperature, forming LiO_x and Ti³⁺ species in the latter case. Finally, the formation of stable decomposition products was found to sensitively change the equilibrium between surface Li and Li⁺ intercalated in the bulk, leading to a deintercalation from lithiated HOPG in the presence of an adsorbed IL adlayer at >230 K. Overall, these results provide detailed insights into the surface chemistry at the solid | electrolyte interface and the initial stages of SEI formation at electrode surfaces in the absence of an applied potential, which is essential for the further improvement of future Li-ion batteries.

Aqueous TMAO solutions as seen by theoretical THz spectroscopy: hydrophilic versus hydrophobic water

All THz resonances of aqueous TMAO solutions are computed and assigned based on ab initio molecular dynamics simulations.

Toward a More Rational Design of the Direct Synthesis of Aniline: A Density Functional Theory Study

Plane-wave density functional theory (PW-DFT) calculations have been used to investigate the direct amination of benzene catalyzed by a Ni(111) surface to explore the reaction intermediates and to understand the role of nickel in this reaction. Adsorption structures, sites, energetics, and proposed reaction pathways relevant to the amination of benzene on the Ni(111) surface were investigated using the spin-polarized slab model with the Perdew-Burke-Ernzerhof functional. The dispersion-corrected DFT-D3 was used to examine the effect of van der Waals interactions on the adsorption energy. Detailed discussion of the adsorption behaviors of NH, NH₂, C₆H₅, C₆H₅NH₂, and C₆H₅NH on the Ni(111) surface is provided. Imide and benzene were predicted to be the most predominant adsorbed species on the Ni(111) surface, and a reaction process involving a surface-bound anilide as an intermediate was predicted to be more thermodynamically favorable than other reaction pathways. The electronic interactions and vibrational frequencies of isolated and adsorbed molecules were also investigated.

Surface configuration and wettability of nickel(oxy)hydroxides: a first-principles investigation

This article explores the wetting behavior of β-type nickel hydroxide, β-Ni(OH)₂, and nickel oxyhydroxide, β-NiOOH, by means of first-principles calculations.