Adsorption of Glycine on TiO2 in Water from On-the-fly Free-Energy Calculations and In Situ Electrochemical Impedance Spectroscopy.

We report here an experimental-computational study of hydrated TiO2 anatase nanoparticles interacting with glycine, where we obtain quantitative agreement of the measured adsorption free energies. Ab initio simulations are performed within the tight binding and density functional theory in combination with enhanced free-energy sampling techniques, which exploit the thermodynamic integration of the unbiased mean forces collected on-the-fly along the molecular dynamics trajectories. The experiments adopt a new and efficient setup for electrochemical impedance spectroscopy measurements based on portable screen-printed gold electrodes, which allows fast and in situ signal assessment. The measured adsorption free energy is -30 kJ/mol (both from experiment and calculation), with preferential interaction of the charged NH3+ group which strongly adsorbs on the TiO2 bridging oxygens. This highlights the importance of the terminal amino groups in the adsorption mechanism of amino acids on hydrated metal oxides. The excellent agreement between computation and experiment for this amino acid opens the doors to the exploration of the interaction free energies for other moderately complex bionano systems.

A Simple Electron-Density Based Force Field Model for High-Energy Interactions between Atoms and Molecules.

In high-energy molecular dynamics or Monte Carlo simulations, standard force fields optimized for simulations at ambient temperatures are inadequate. This is largely because their repulsive parts have been regarded as not very significant, even well below zero interaction energies. It is, therefore, not obvious which force fields to resort to for simulating hot gases or plasmas. A force field model that uses the electronic densities of noninteracting atoms or molecules within the pair approximation is introduced. We start by deriving a naïve model that neglects any exchange and correlation effects between the electronic clouds and then correct this model by adding a term calibrated from ab initio calculations using the CCSD(T)/cc-pVTZ level of theory. The resulting expression for this term can be regarded as a simple exchange-correlation function. We compare the results for the repulsive part of the potential energy hypersurfaces with the force fields commonly used on some dimers of small molecules.

Evolutionary Monte Carlo of QM Properties in Chemical Space: Electrolyte Design.

Optimizing a target function over the space of organic molecules is an important problem appearing in many fields of applied science but also a very difficult one due to the vast number of possible molecular systems. We propose an evolutionary Monte Carlo algorithm for solving such problems which is capable of straightforwardly tuning both exploration and exploitation characteristics of an optimization procedure while retaining favorable properties of genetic algorithms. The method, dubbed MOSAiCS (Metropolis Optimization by Sampling Adaptively in Chemical Space), is tested on problems related to optimizing components of battery electrolytes, namely, minimizing solvation energy in water or maximizing dipole moment while enforcing a lower bound on the HOMO-LUMO gap; optimization was carried out over sets of molecular graphs inspired by QM9 and Electrolyte Genome Project (EGP) data sets. MOSAiCS reliably generated molecular candidates with good target quantity values, which were in most cases better than the ones found in QM9 or EGP. While the optimization results presented in this work sometimes required up to 106 QM calculations and were thus feasible only thanks to computationally efficient ab initio approximations of properties of interest, we discuss possible strategies for accelerating MOSAiCS using machine learning approaches.

H2O2(s) and H2O2·2H2O(s) crystals compared with ices: DFT functional assessment and D3 analysis.

The H2O and H2O2 molecules resemble each other in a multitude of ways as has been noted in the literature. Here, we present density functional theory (DFT) calculations for the H2O2(s) and H2O2·2H2O(s) crystals and make selected comparisons with ice polymorphs. The performance of a number of dispersion-corrected density functionals-both self-consistent and a posteriori ones-are assessed, and we give special attention to the D3 correction and its effects. The D3 correction to the lattice energies is large: for H2O2(s) the D3 correction constitutes about 25% of the lattice energy using PBE, much more for RPBE, much less for SCAN, and it primarily arises from non-H-bonded interactions out to about 5 Å.The large D3 corrections to the lattice energies are likely a consequence of several effects: correction for missing dispersion interaction, the ability of D3 to capture and correct various other kinds of limitations built into the underlying DFT functionals, and finally some degree of cell-contraction-induced polarization enhancement. We find that the overall best-performing functionals of the twelve examined are optPBEvdW and RPBE-D3. Comparisons with DFT assessments for ices in the literature show that where the same methods have been used, the assessments largely agree.

Water on ceria{111}: Comparison between 23 experimental vibrational studies in the literature and new modeling.

Theoretical and experimental vibrational signatures of H2O and OH- (dissociated water) adsorbed on stoichiometric ceria{111} surfaces are compared. The experimental ones were collected from low-coverage experiments in the literature, and the theoretical anharmonic frequencies were generated using density functional theory calculations employing the optPBE-vdW functional for coverages from 0.5 to a few monolayers. It is found that (i) the experiments and our calculations overall agree well, lending credibility to both; (ii) the calculations manage to resolve the large class of H-bonded motifs into frequency classes that can guide experimental assignments; (iii) it is possible to find a geometrical H-bond definition that also captures the OH vibrational frequency downshifts well: R(H⋯O) ≤2.5 Å and the O-H⋯O angle θ ≥ 100°; and (iv) the frequency vs electric field relations for water and hydroxides (i.e., dissociated water) follow different and well-separated curves.

Origin of the Hydrophobic Behaviour of Hydrophilic CeO2.

The nature of the hydrophobicity found in rare-earth oxides is intriguing. The CeO2 (100) surface, despite its strongly hydrophilic nature, exhibits hydrophobic behaviour when immersed in water. In order to understand this puzzling and counter-intuitive effect we performed a detailed analysis of the confined water structure and dynamics. We report here an ab-initio molecular dynamics simulation (AIMD) study which demonstrates that the first adsorbed water layer, in immediate contact with the hydroxylated CeO2 surface, generates a hydrophobic interface with respect to the rest of the liquid water. The hydrophobicity is manifested in several ways: a considerable diffusion enhancement of the confined liquid water as compared with bulk water at the same thermodynamic condition, a weak adhesion energy and few H-bonds above the hydrophobic water layer, which may also sustain a water droplet. These findings introduce a new concept in water/rare-earth oxide interfaces: hydrophobicity mediated by specific water patterns on a hydrophilic surface.

Toward an efficient f-in-core/f-in-valence switchable description for DFTB calculations of Ce 4f states in ceria.

A computational protocol is developed for efficient studies of partially reduced redox-active oxides using the self-consistent charge density functional tight-binding method. The protocol is demonstrated for ceria, which is a prototypical reducible oxide material. The underlying idea is to achieve a consistent (and harmonized) set of Slater-Koster (SK) tables with connected repulsive potentials that enable switching on and off the in-valence description of the Ce 4f states without serious loss of accuracy in structure and energetics. The implicit treatment of the Ce 4f states, with the use of f-in-core SK-tables, is found to lead to a significant decrease in computational time. More importantly, it allows for explicit control of the oxidation states of individual Ce atoms. This makes it possible to "freeze" the electronic configuration, thereby allowing the exploration of the energetics for various meta-stable configurations. We anticipate that the outlined strategy can help to shed light on the interplay between the size, shape, and redox activity for nanoceria and other related materials.

Predicting Frequency from the External Chemical Environment: OH Vibrations on Hydrated and Hydroxylated Surfaces.

Robust correlation curves are essential to decipher structural information from IR-vibrational spectra. However, for surface-adsorbed water and hydroxides, few such correlations have been presented in the literature. In this paper, OH vibrational frequencies are correlated against 12 structural descriptors representing the quantum mechanical or geometrical environment, focusing on those external to the vibrating molecule. A nonbiased fitting procedure based on Gaussian process regression (GPR) was used alongside simple analytical functional forms. The training data consist of 217 structurally unique OH groups from 38 water/metal oxide interface systems for MgO, CaO and CeO2, all optimized at the DFT level, and the fully anharmonic and uncoupled OH vibrational signatures were calculated. Among our results, we find the following: (i) The intermolecular R(H···O) hydrogen bond distance is particularly strong, indicating the primary cause of the frequency shift. (ii) Similarly, the electric field along the H-bond vector is also a good descriptor. (iii) Highly detailed machine learning descriptors (ACSF, SOAP) are less intuitive but were found to be more capable descriptors. (iv) Combinations of geometric and QM descriptors give the best predictions, supplying complementary information.

Space-Resolved OH Vibrational Spectra of the Hydration Shell around CO2.

The CO2 molecule is weakly bound in water. Here we analyze the influence of a dissolved CO2 molecule on the structure and OH vibrational spectra of the surrounding water. From the analysis of ab initio molecular dynamics simulations (BLYP-D3) we present static (structure, coordination, H-bonding, tetrahedrality) and dynamical (OH vibrational spectra) properties of the water molecules as a function of distance from the solute. We find a weakly oscillatory variation ("ABBA") in the 'solution minus bulk water' spectrum. The origin of these features can largely be traced back to solvent-solute hard-core interactions which lead to variations in density and tetrahedrality when moving from the solute's vicinity out to the bulk region. The high-frequency peak in the solute-affected spectra is specifically analyzed and found to originate from both water OH groups that fulfill the geometric H-bond criteria, and from those that do not (dangling ones). Effectively, neither is hydrogen-bonded.

Vertical-Flow Paper Sensor for On-Site and Prompt Evaluation of Chloride Contamination in Concrete Structures.

Corrosion occurring in reinforced concrete has turned into a primary concern of the current century, concrete being the most ubiquitous and predominant material used in the construction industry. Among the many interrelated processes that trigger corrosion of metallic reinforcements, the penetration of chloride ions into the concrete matrix is the most insidious threat. Herein, we developed the first electrochemical device entirely made of paper that allows for the direct, prompt, and noninvasive evaluation of free chloride ion contamination in concrete-based constructions. Our device is based on a three-layer wax-modified filter paper, consisting of two Ag/AgCl screen-printed electrodes that are interfaced by a junction pad in a sandwich-like configuration. Filter paper allows for generating a vertical-flow potentiometric device capable of measuring the electrochemical potential between two solutions containing different concentrations of chloride ions, which are separately drop-cast on the top and bottom layers. After demonstrating the analytical performance of the device, the same principle was applied to the evaluation of the chloride contents in different concrete samples, exploiting paper as a suitable interfacing material for potentiometric measurements on the cement solid surface. Laboratory-prepared concrete samples with known chloride contents were first assessed, and then, the paper-based vertical-flow device was applied to real concrete structures at the Giacomo Manzù Museum (Ardea, Italy) for the evaluation of chloride contamination caused by the proximity to the seaside. The capability of our device to provide timely warning of the risk conditions of concrete-based artifacts was demonstrated.

Oxygen chemistry of halogen-doped CeO2(111).

We study substitutional fluorine, chlorine and bromine impurities at CeO2(111), and their effects on the oxygen chemistry of the surface, using density functional theory. We find that impurity formation results in a halide ion and one Ce3+ ion for all three halogens, although the formation energy depends strongly on the identity of the halogen; however, once formed, all three halogens exhibit a similar propensity to form impurity-impurity pairs. Furthermore, while the effects of halogen impurities on oxygen vacancy formation are marginal, they are more significant for oxygen molecule adsorption, due to electron transfer from the Ce3+ ion which results in an adsorbed superoxide molecule. We also consider the displacement of a halide ion on to the surface by half of an oxygen molecule, and find that the energy required to do so depends strongly not only on the identity of the halogen, but also on whether or not a second halogen impurity, with its associated Ce3+ ion, is present; if it is, then the process is greatly facilitated. Overall, our results demonstrate the existence of a rich variety of ways in which the oxygen chemistry of CeO2(111) may be modified by the presence of halogen dopants.

Supercooled liquid-like dynamics in water near a fully hydrated titania surface: Decoupling of rotational and translational diffusion.

We report an ab initio molecular dynamics (MD) simulation investigating the effect of a fully hydrated surface of TiO2 on the water dynamics. It is found that the universal relation between the rotational and translational diffusion characteristics of bulk water is broken in the water layers near the surface with the rotational diffusion demonstrating progressive retardation relative to the translational diffusion when approaching the surface. This kind of rotation-translation decoupling has so far only been observed in the supercooled liquids approaching glass transition, and its observation in water at a normal liquid temperature is of conceptual interest. This finding is also of interest for the application-significant studies of the water interaction with fully hydrated nanoparticles. We note that this is the first observation of rotation-translation decoupling in an ab initio MD simulation of water.

Lignin Intermediates on Palladium: Insights into Keto-Enol Tautomerization from Theoretical Modelling.

It has been suggested in the literature that keto-to-enol tautomerization plays a vital role for lignin fragmentation under mild conditions. On the other hand, previous modelling has shown that the adsorbed keto form is more stable than enol on the Pd(111) catalyst. The current density functional theory study of lignin model molecules shows that, in the gas-phase, keto is more stable than enol, but on the Pd surface, we find enol conformers that are at least as stable as keto. This supports the experimental result that the favourable reaction pathway for lignin depolymerization involves keto-enol tautomerization. An energy decomposition analysis gives insights concerning the origin of the fine energy balance between the keto and enol forms, where the molecule-surface interaction (-7 eV) and the molecular strain energy (+3 eV) are the main contributors to the adsorption energy.

Tuning the ORR activity of Pt-based Ti2CO2 MXenes by varying the atomic cluster size and doping with metals.

The rational design of ideal catalysts for the oxygen reduction reaction (ORR) is of great significance for solving the electrocatalytic potential problems in proton exchange membrane fuel cells (PEMFCs). Ptn (n = 1-4) and Pt3Au alloy subnanoclusters supported on a defective Ti2CO2 monolayer with oxygen vacancies (denoted as v-Ti2CO2) are simulated by using density functional theory to investigate their ORR performance. The geometries, energetics, and electronic properties of the different systems are analyzed. It is found that the supported Pt3Au alloy subnanocluster possesses the best ORR activity. The underlying mechanisms of the improved ORR activity originates from the moderate hybridization between the O 2p and the 5d orbitals of Au and Pt according to the density of states analysis. Our study suggests a facile route for designing low-cost MXene-based electrocatalysts by alloying transition metals with Pt catalysts, which may stimulate realization of suitable alternative catalysts for ORR catalysis.

The water/ceria(111) interface: Computational overview and new structures.

Thin film structures of water on the CeO2(111) surface for coverages between 0.5 and 2.0 water monolayers have been optimized and analyzed using density functional theory (optPBE-vdW functional). We present a new 1.0 ML structure that is both the lowest in energy published and features a hydrogen-bond network extending the surface in one-dimension, contrary to what has been found in the literature, and contrary to what has been expected due to the large bulk ceria cell dimension. The adsorption energies for the monolayer and multilayered water structures agree well with experimental temperature programmed desorption results from the literature, and we discuss the stability window of CeO2(111) surfaces covered with 0.5-2.0 ML of water.

Temperature effects on the ionic conductivity in concentrated alkaline electrolyte solutions.

Alkaline electrolyte solutions are important components in rechargeable batteries and alkaline fuel cells. As the ionic conductivity is thought to be a limiting factor in the performance of these devices, which are often operated at elevated temperatures, its temperature dependence is of significant interest. Here we use NaOH as a prototypical example of alkaline electrolytes, and for this system we have carried out reactive molecular dynamics simulations with an experimentally verified high-dimensional neural network potential derived from density-functional theory calculations. It is found that in concentrated NaOH solutions elevated temperatures enhance both the contributions of proton transfer to the ionic conductivity and deviations from the Nernst-Einstein relation. These findings are expected to be of practical relevance for electrochemical devices based on alkaline electrolyte solutions.

Anion-mediated electronic effects in reducible oxides: Tuning the valence band of ceria via fluorine doping.

Combining experimental spectroscopy and hybrid density functional theory calculations, we show that the incorporation of fluoride ions into a prototypical reducible oxide surface, namely, ceria(111), can induce a variety of nontrivial changes to the local electronic structure, beyond the expected increase in the number of Ce3+ ions. Our resonant photoemission spectroscopy results reveal new states above, within, and below the valence band, which are unique to the presence of fluoride ions at the surface. With the help of hybrid density functional calculations, we show that the different states arise from fluoride ions in different atomic layers in the near surface region. In particular, we identify a structure in which a fluoride ion substitutes for an oxygen ion at the surface, with a second fluoride ion on top of a surface Ce4+ ion giving rise to F 2p states which overlap the top of the O 2p band. The nature of this adsorbate F--Ce4+ resonant enhancement feature suggests that this bond is at least partially covalent. Our results demonstrate the versatility of anion doping as a potential means of tuning the valence band electronic structure of ceria.

Multiscale Modeling of Agglomerated Ceria Nanoparticles: Interface Stability and Oxygen Vacancy Formation.

The interface formation and its effect on redox processes in agglomerated ceria nanoparticles (NPs) have been investigated using a multiscale simulation approach with standard density functional theory (DFT), the self-consistent-charge density functional tight binding (SCC-DFTB) method, and a DFT-parameterized reactive force-field (ReaxFF). In particular, we have modeled Ce40O80 NP pairs, using SCC-DFTB and DFT, and longer chains and networks formed by Ce40O80 or Ce132O264 NPs, using ReaxFF molecular dynamics simulations. We find that the most stable {111}/{111} interface structure is coherent whereas the stable {100}/{100} structures can be either coherent or incoherent. The formation of {111}/{111} interfaces is found to have only a very small effect on the oxygen vacancy formation energy, E vac. The opposite holds true for {100}/{100} interfaces, which exhibit significantly lower E vac values than the bare surfaces, despite the fact that the interface formation eliminates reactive {100} facets. Our results pave the way for an increased understanding of ceria NP agglomeration.

STM Images of Anionic Defects at CeO2(111)-A Theoretical Perspective.

We present a theoretically oriented analysis of the appearance and properties of plausible candidates for the anionic defects observed in scanning tunneling microscopy (STM) experiments on CeO2(111). The simulations are based on density functional theory (DFT) and cover oxygen vacancies, fluorine impurities and hydroxyl groups in the surface and sub-surface layers. In the surface layer, all three appear as missing spots in the oxygen sublattice in filled state simulated STM images, but they are distinguishable in empty state images, where surface oxygen vacancies and hydroxyls appear as, respectively, diffuse and sharp bright features at oxygen sites, while fluorine defects appear as triangles of darkened Ce ions. In the sub-surface layer, all three defects present more complex patterns, with different combinations of brightened oxygen ion triangles and/or darkened Ce ion triangles, so we provide image maps to support experimental identification. We also discuss other properties that could be used to distinguish the defects, namely their diffusion rates and distributions.

Fifty Shades of Water: Benchmarking DFT Functionals against Experimental Data for Ionic Crystalline Hydrates.

We propose that crystalline ionic hydrates constitute a valuable resource for benchmarking theoretical methods for aqueous ionic systems. Many such structures are known from the experimental literature, and they contain a large variety of water-water and ion-water structural motifs. Here we have collected a data set (CRYSTALWATER50) of 50 structurally unique "in-crystal" water molecules, involved in close to 100 nonequivalent O-H···O hydrogen bonds. A dozen well-known DFT functionals were benchmarked with respect to their ability to describe these experimental structures and their OH vibrational frequencies. We find that the PBE, RPBE-D3, and optPBE-vdW methods give the best H-bond distances and that anharmonic OH frequencies generated from B3LYP//optPBE-vdW energy scans outperform the other methods, i.e., here we performed B3LYP energy scans along the OH stretching coordinate while the rest of the structure was kept fixed at the optPBE-vdW-optimized positions.