On the Range of Water Structure Models Compatible with X-ray and Neutron Diffraction Data

Revisiting the activity origin of the PtAu24(SR)18 nanocluster for enhanced electrocatalytic hydrogen evolution by combining first-principles simulations with the experimental in situ FTIR technique

Thiolate-protected metal nanoclusters (NCs) have been widely used in various electrocatalytic reactions, yet the dynamic evolution of metal NCs during electrocatalysis has been rarely explored and the activity origin remains largely ambiguous. Herein, using a PtAu24(SCH3)18 NC as a prototype model, we combined advanced first-principles calculations and attenuated total reflection surface-enhanced infrared spectroscopy (ATR-SEIRAS) to re-examine its active site and reaction dynamics in the hydrogen evolution reaction (HER). It has been previously assumed that the central Pt is the only catalytic center. However, differently, we observed the spontaneous desorption of thiolate ligands under moderate potential, and the dethiolated PtAu24 exhibits excellent HER activity, which is contributed not only by the central Pt atom but also by the exposed bridged Au sites. Particularly, the exposed Au exhibits high activity even comparable to Pt, and the synergistic effect between them makes dethiolated PtAu24 an extraordinary HER electrocatalyst, even surpassing the commercial Pt/C catalyst. Our predictions are further verified by electrochemical activation experiments and in situ FTIR (ATR-SEIRAS) characterization, where evident adsorption of Au-H* and Pt-H* bonds is monitored. This work detected, for the first time, the Au-S interfacial dynamics of the PtAu24 nanocluster in electrocatalytic processes, and quantitatively evaluated the essential catalytic role of the exposed Au sites that has been largely overlooked in previous studies.

-SR removal or -R removal? A mechanistic revisit on the puzzle of ligand etching of Au25(SR)18 nanoclusters during electrocatalysis

Accurate identification of active sites is highly desirable for elucidation of the reaction mechanism and development of efficient catalysts. Despite the promising catalytic performance of thiolated metal nanoclusters (NCs), their actual catalytic sites remain elusive. Traditional first-principles calculations and experimental observations suggested dealkylated S and dethiolated metal, respectively, to be the active centers. However, the real kinetic origin of thiolate etching during the electrocatalysis of NCs is still puzzling. Herein, we conducted advanced first-principles calculations and electrochemical/spectroscopic experiments to unravel the electrochemical etching kinetics of thiolate ligands in prototype Au25(SCH3)18 NC. The electrochemical processes are revealed to be spontaneously facilitated by dethiolation (i.e., desorption of -SCH3), forming the free HSCH3 molecule after explicitly including the solvent effect and electrode potential. Thus, exposed under-coordinated Au atoms, rather than the S atoms, serve as the real catalytic sites. The thermodynamically preferred Au-S bond cleavage arises from the selective attack of H from proton/H2O on the S atom under suitable electrochemical bias due to the spatial accessibility and the presence of S lone pair electrons. Decrease of reduction potential promotes the proton attack on S and significantly accelerates the kinetics of Au-S bond breakage irrespective of the pH of the medium. Our theoretical results are further verified by the experimental electrochemical and spectroscopic data. At more negative electrode potentials, the number of -SR ligands decreased with concomitant increase of the vibrational intensity of S-H bonds. These findings together clarify the atomic-level activation mechanism on the surface of Au25(SR)18 NCs.

Thermodynamic calculations using reverse Monte Carlo: Simultaneously tuning multiple short-range order parameters for 2D lattice adsorption problem

Lattice simulations are an important class of problems in crystalline solids, surface science, alloys, adsorption, absorption, separation, catalysis, to name a few. We describe a fast computational method for performing lattice thermodynamic calculations that is based on the use of the reverse Monte Carlo (RMC) technique and multiple short-range order (SRO) parameters. The approach is comparable in accuracy to the Metropolis Monte Carlo (MC) method. The equilibrium configuration is determined in 5-10 Newton-Raphson iterations by solving a system of coupled nonlinear algebraic SRO growth rate equations. This makes the RMC-based method computationally more efficient than MC, given that MC typically requires sampling of millions of configurations. The technique is applied to the interacting 2D adsorption problem. Unlike grand canonical MC, RMC is found to be adept at tackling geometric frustration, as it is able to quickly and correctly provide the ordered c(2 × 2) adlayer configuration for Cl adsorbed on a Cu (100) surface.

A probabilistic microkinetic modeling framework for catalytic surface reactions

We present a probabilistic microkinetic modeling (MKM) framework that incorporates the short-ranged order (SRO) evolution for adsorbed species (adspecies) on a catalyst surface. The resulting model consists of a system of ordinary differential equations. Adsorbate-adsorbate interactions, surface diffusion, adsorption, desorption, and catalytic reaction processes are included. Assuming that the adspecies ordering/arrangement is accurately described by the SRO parameters, we employ the reverse Monte Carlo (RMC) method to extract the relevant local environment probability distributions and pass them to the MKM. The reaction kinetics is faithfully captured as accurately as the kinetic Monte Carlo (KMC) method but with a computational time requirement of few seconds on a standard desktop computer. KMC, on the other hand, can require several days for the examples discussed. The framework presented here is expected to provide the basis for wider application of the RMC-MKM approach to problems in computational catalysis, electrocatalysis, and material science.

Hydration in aqueous NaCl

Atomistic details about the hydration of ions in aqueous solutions are still debated due to the disordered and statistical nature of the hydration process. However, many processes from biology, physical chemistry to materials sciences rely on the complex interplay between solute and solvent. Oxygen K-edge X-ray excitation spectra provide a sensitive probe of the local atomic and electronic surrounding of the excited sites. We used ab initio molecular dynamics simulations together with extensive spectrum calculations to relate the features found in experimental oxygen K-edge spectra of a concentration series of aqueous NaCl with the induced structural changes upon solvation of the salt and distill the spectral fingerprints of the first hydration shells around the Na⁺- and Cl^--ions. By this combined experimental and theoretical approach, we find the strongest spectral changes to indeed result from the first hydration shells of both ions and relate the observed shift of spectral weight from the post- to the main-edge to the origin of the post-edge as a shape resonance.

Pseudo-adsorption and long-range redox coupling during oxygen reduction reaction on single atom electrocatalyst

Fundamental understanding of the dynamic behaviors at the electrochemical interface is crucial for electrocatalyst design and optimization. Here, we revisit the oxygen reduction reaction mechanism on a series of transition metal (M = Fe, Co, Ni, Cu) single atom sites embedded in N-doped nanocarbon by ab initio molecular dynamics simulations with explicit solvation. We have identified the dissociative pathways and the thereby emerged solvated hydroxide species for all the proton-coupled electron transfer (PCET) steps at the electrochemical interface. Such hydroxide species can be dynamically confined in a "pseudo-adsorption" state at a few water layers away from the active site and respond to the redox event at the catalytic center in a coupled manner within timescale less than 1 ps. In the PCET steps, the proton species (in form of hydronium in neutral/acidic media or water in alkaline medium) can protonate the pseudo-adsorbed hydroxide without needing to travel to the direct catalyst surface. This, therefore, expands the reactive region beyond the direct catalyst surface, boosting the reaction kinetics via alleviating mass transfer limits. Our work implies that in catalysis the reaction species may not necessarily bind to the catalyst surface but be confined in an active region.

Structural Analysis of Molecular Materials Using the Pair Distribution Function

This is a review of atomic pair distribution function (PDF) analysis as applied to the study of molecular materials. The PDF method is a powerful approach to study short- and intermediate-range order in materials on the nanoscale. It may be obtained from total scattering measurements using X-rays, neutrons, or electrons, and it provides structural details when defects, disorder, or structural ambiguities obscure their elucidation directly in reciprocal space. While its uses in the study of inorganic crystals, glasses, and nanomaterials have been recently highlighted, significant progress has also been made in its application to molecular materials such as carbons, pharmaceuticals, polymers, liquids, coordination compounds, composites, and more. Here, an overview of applications toward a wide variety of molecular compounds (organic and inorganic) and systems with molecular components is presented. We then present pedagogical descriptions and tips for further implementation. Successful utilization of the method requires an interdisciplinary consolidation of material preparation, high quality scattering experimentation, data processing, model formulation, and attentive scrutiny of the results. It is hoped that this article will provide a useful reference to practitioners for PDF applications in a wide realm of molecular sciences, and help new practitioners to get started with this technique.

Comparison of molecular dynamics simulations of water with neutron and X-ray scattering experiments

The atomistic structure and dynamics obtained from molecular dynamics (MD) simulations with the example of TIP3P (rigid and flexible) and TIP4P/2005 (rigid) water is compared to neutron and X-ray scattering data at ambient conditions. Neutron and X-ray diffractograms are calculated from the simulations for four isotopic substitutions as well as the incoherent intermediate scattering function for neutrons. The resulting curves are compared to each other and to published experimental data. Differences between simulated and measured intermediate scattering functions are quantified by fitting an analytic model to the computed values. The sensitivity of the scattering curves to the parameters of the MD simulations is demonstrated on the example of two parameters, bond length and angle.

Simulation of Liquids with the Tight-Binding Density-Functional Approach and Improved Atomic Charges

Theoretical description of liquids, especially liquid water, is an ongoing subject with important implications in various domains such as homogeneous catalysis; solvation of molecular, ionic, and biomolecular species; and reactivity. Various formalisms exist to describe liquids, each one displaying its own balance between accuracy and computational cost that defines its range of applications. The present article revisits the ability of the density-functional-based tight-binding (SCC-DFTB) approach to model liquids by focusing on liquid water and liquid benzene under ambient conditions. To do so, we benchmark a recent correction for the SCC-DFTB atomic charges that allows for a drastic improvement of the pair radial distribution functions of liquid water as compared to both experimental data and density-functional theory results performed in the generalized-gradient approximation. We also report the coupling of the deMonNano and i-PI codes to perform path-integral molecular dynamics. This allows us to rationalize the impact of nuclear quantum effects on the SCC-DFTB description of liquid water. This study evidences the rather good ability of SCC-DFTB to describe liquid water and liquid benzene. As the first example of application, we also present results for a benzene molecule solvated in water with the perspectives of further studies devoted to solvent/water interfaces.

Decoding signatures of structure, bulk thermodynamics, and solvation in three-body angle distributions of rigid water models

A tetrahedral structure resulting from hydrogen bonding is a hallmark of liquid water and plays a significant role in determining its unique thermophysical properties. This water feature has helped understand anomalous properties and physically interpret and model hydrophobic solvation thermodynamics. Tetrahedrality is well described by the geometric relationship of any central water molecule with two of its nearest neighbors in the first coordination shell, as defined by the corresponding "three-body" angle. While order parameters and even full water models have been developed using specific or average features of the three-body angle distribution, here we examine the distribution holistically, tracking its response to changes in temperature, density, and the presence of model solutes. Surprisingly, we find that the three-body distribution responds by varying primarily along a single degree of freedom, suggesting a remarkably simplified view of water structure. We characterize three-body angle distributions across temperature and density space and identify principal components of the variations with state conditions. We show that these principal components embed physical significance and trace out transitions between tetrahedral and simple-fluid-like behavior. Moreover, we find that the ways three-body angles vary within the hydration shells of model colloids of different types and sizes are nearly identical to the variations seen in bulk water across density and temperature. Importantly, through the principal directions of these variations, we find that perturbations to the hydration-water distributions well predict the thermodynamics associated with colloid solvation, in particular, the relative entropy of this process that captures indirect, solvent-mediated contributions to the hydration free energy.

Intermediate range O-O correlations in supercooled water down to 235 K

Wide angle x-ray scattering of supercooled water down to 234.8 K was studied using high energy x rays at the European Synchrotron Radiation Facility. The oxygen-oxygen pair distribution function (PDF) was calculated from the scattering pattern out to the 5th peak at an intermolecular distance, r ≈ 11 Å. We observe that the 4th peak and the 5th peak in the PDF increase in height upon supercooling. We also observe that the 4th peak position (r₄) shifts to shorter distances upon supercooling consistent with previous studies, but we see a more rapid change at the lowest temperature. The running oxygen-oxygen coordination number is calculated for 5 different temperatures, and an isosbestic point at r_iso = 3.31 ± 0.05 Å was found corresponding to a coordination number of 4.39 ± 0.15. The comparison of the PDF of the coldest water with that of amorphous ice shows distinct differences. We propose that there are 5-member pentamer rings in low density liquid-like structures giving rise to the sharp correlations at r ≈ 9 Å and r ≈ 11 Å.

Understanding the microscopic origin behind heterogeneous properties of water confined in and around Aβ17-42 protofilaments

Aggregation of amyloid beta (Aβ) peptides in the brain is responsible for one of the most devastating neurodegenerative diseases, namely, Alzheimer's disease. In this study, we have carried out atomistic molecular dynamics simulations to explore the effects of non-uniform structural distortions of Aβ_17-42 pre-fibrillar aggregates of different sizes on the microscopic structure and ordering of water molecules confined within their amphiphilic nanocores. The calculations revealed non-uniform peptide-water interactions resulting in simultaneous existence of both highly ordered and disordered water molecules within the spatially heterogeneous confined environment of the protofilament cores. It is found that the high degree of ordering originates from a sizable fraction of doubly coordinated core water molecules, while the randomly oriented ones are those that are coordinated with three neighbors in their first coordination shells. Furthermore, it is quantitatively demonstrated that relative fractions of these two types of water molecules are correlated with the protofilament core topology and the degree of confinement within that. It is proposed that the ordered core waters are likely to stabilize the Aβ protofilaments by screening the residue charges and favoring water-mediated salt bridge formations, while the randomly oriented ones can drive further growth of the protofilaments by being displaced easily during the docking of additional peptides. In that way, both types of core water molecules can play equally important roles in controlling the growth and stability of the Aβ-aggregates.

Magnetic properties and core electron binding energies of liquid water

The magnetic properties and the core and inner valence electron binding energies of liquid water are investigated. The adopted methodology relies on the combination of molecular dynamics and electronic structure calculations. Born-Oppenheimer molecular dynamics with the Becke and Lee-Yang-Parr functionals for exchange and correlation, respectively, and includes an empirical correction (BLYP-D3) functional and classical molecular dynamics with the TIP4P/2005-F model were carried out. The Keal-Tozer functional was applied for predicting magnetic shielding and spin-spin coupling constants. Core and inner valence electron binding energies in liquid water were calculated with symmetry adapted cluster-configuration interaction. The relationship between the magnetic shielding constant σ(¹⁷O), the role played by the oxygen atom as a proton acceptor and donor, and the tetrahedral organisation of liquid water are investigated. The results indicate that the deshielding of the oxygen atom in water is very dependent on the order parameter (q) describing the tetrahedral organisation of the hydrogen bond network. The strong sensitivity of magnetic properties on changes of the electronic density in the nuclei environment is illustrated by a correlation between σ(¹⁷O) and the energy gap between the 1a₁[O1s] (core) and the 2a₁ (inner valence) orbitals of water. Although several studies discussed the eventual connection between magnetic properties and core electron binding energies, such a correlation could not be clearly established. Here, we demonstrate that for liquid water this correlation exists although involving the gap between electron binding energies of core and inner valence orbitals.

Probing the triplet correlation function in liquid water by experiments and molecular simulations

Three-body information of liquid water is extracted using X-ray diffraction experiment as well as in molecular simulations via isothermal pressure derivative of structure factor term.

Dowser++, a new method of hydrating protein structures

A new method of hydrating protein structures, which we call Dowser++, is presented. The method is based on a semi-empirical modification of a popular program for protein hydration Dowser, and the usage of protocols AutoDock Vina, and WaterDock. The positions of water molecules predicted by Dowser++ were compared with experimental data for a set of 14 high-resolution crystal structures of oligopeptide-binding protein (OppA) containing a large number of resolved internal water molecules, as well as for the D- and K-channels of cytochrome c oxidase, and the recent data on PSII. Comparison is also made with the predictions of the original Dowser, and its improved version, Dowser+, described in our previous publication. We also present a model for quantitative estimation of the quality of water molecules placement made by a program, which includes an assumption of possible false negative data from the crystallographic analysis. The comparison of predictions made by Dowser++, Dowser and Dowser+ demonstrates significant improvement of predictive power of the new method. Proteins 2016; 84:1347-1357. © 2016 Wiley Periodicals, Inc.

Potential Paths for the Hydrogen-Bond Relaxing With (H2O)N Cluster Size

Relaxation of the inter- and intra-molecular interactions for the hydrogen bond (O:H-O) between undercoordinated molecules determines the unusual behavior of water nanodroplets and nanobubbles. However, probing such potentials remains unreality. Here we show that the Lagrangian solution [Huang et al., J. Phys. Chem. B, 2013. 117: 13639] transforms the observed H-O bond (x = H) and O:H nonbond (x = L) lengths and phonon frequencies (dx, x) [Sun et al., J. Phys. Chem. Lett., 2013. 4: 2565] into the respective force constants and bond energies (kx, Ex) and hence enables the mapping of the potential paths for the O:H-O bond relaxing with water cluster size. Results show that molecular undercoordination not only reduces the molecular size (dH) with enhanced H-O energy from the bulk value of 3.97 to 5.10 eV for a H2O monomer, but also enlarges the molecular separation (dL) with reduced O:H energy from 95 to 35 meV for a dimer. The H-O energy gain raises the melting point from bulk value 273 to 310 K for the skin and the O:H energy loss lowers the freezing temperature from bulk value 258 to 202 K for 1.4 nm sized droplet, by dispersing the quasisolid phase boundaries.

Triplet correlation functions in liquid water

Triplet correlations have been shown to play a crucial role in the transformation of simple liquids to anomalous tetrahedral fluids [M. Singh, D. Dhabal, A. H. Nguyen, V. Molinero, and C. Chakravarty, Phys. Rev. Lett. 112, 147801 (2014)]. Here we examine triplet correlation functions for water, arguably the most important tetrahedral liquid, under ambient conditions, using configurational ensembles derived from molecular dynamics (MD) simulations and reverse Monte Carlo (RMC) datasets fitted to experimental scattering data. Four different RMC data sets with widely varying hydrogen-bond topologies fitted to neutron and x-ray scattering data are considered [K. T. Wikfeldt, M. Leetmaa, M. P. Ljungberg, A. Nilsson, and L. G. M. Pettersson, J. Phys. Chem. B 113, 6246 (2009)]. Molecular dynamics simulations are performed for two rigid-body effective pair potentials (SPC/E and TIP4P/2005) and the monatomic water (mW) model. Triplet correlation functions are compared with other structural measures for tetrahedrality, such as the O-O-O angular distribution function and the local tetrahedral order distributions. In contrast to the pair correlation functions, which are identical for all the RMC ensembles, the O-O-O triplet correlation function can discriminate between ensembles with different degrees of tetrahedral network formation with the maximally symmetric, tetrahedral SYM dataset displaying distinct signatures of tetrahedrality similar to those obtained from atomistic simulations of the SPC/E model. Triplet correlations from the RMC datasets conform closely to the Kirkwood superposition approximation, while those from MD simulations show deviations within the first two neighbour shells. The possibilities for experimental estimation of triplet correlations of water and other tetrahedral liquids are discussed.

A continuous mixture of two different dimers in liquid water

Liquid water is formed by a continuous mixture of two different dimers (cis and trans) with distinct energies related to different relative water molecule orientations.

Natural polarizability and flexibility via explicit valency: the case of water

As the dominant physiological solvent, water drives the folding of biological macromolecules, influences conformational changes, determines the ionization states of surface groups, actively participates in catalytic events, and provides "wires" for long-range proton transfer. Elucidation of all these roles calls for atomistic simulations. However, currently available methods do not lend themselves to efficient simulation of proton transfer events, or even polarizability and flexibility. Here, we report that an explicit account of valency can provide a unified description for the polarizability, flexibility, and dissociability of water in one intuitive and efficient setting. We call this approach LEWIS, after the chemical theory that inspires the use of valence electron pairs. In this paper, we provide details of the method, the choice of the training set, and predictions for the neat ambient liquid, with emphasis on structure, dynamics, and polarization. LEWIS water provides a good description of bulk properties, and dipolar and quadrupolar responses.

The Structure of Hydrated Poly (D, L - Lactic Acid) Studied With X-Ray Diffraction and Molecular Simulation Methods

The effect of hydration on the molecular structure of amorphous poly (D, L-lactic acid) (PDLLA) with 50:50 L-to-D ratio has been studied by combining experiments with molecular simulations. X-ray diffraction measurements revealed significant changes upon hydration in the structure functions of the copolymer. Large changes in the structure functions at ~ 10 days of incubation coincided with the large increase in the water uptake from ~1 to ~40% and the formation of voids in the film. Computer modeling based on the recently developed TIGER2/TIGER3 mixed sampling scheme was used to interpret these changes by efficiently equilibrating both dry and hydrated models of PDLLA. Realistic models of bulk amorphous PDLLA structure were generated as demonstrated by close agreement between the calculated and the experimental structure functions. These molecular simulations were used to identify the interactions between water and the polymer at the atomic level including the change of positional order between atoms in the polymer due to hydration. Changes in the partial O-O structure functions, about 95% of which were due to water-polymer interactions, were apparent in the radial distribution functions. These changes, and somewhat smaller changes in the C-C and C-O partial structure functions, clearly demonstrated the ability of the model to capture the hydrogen bonding interactions between water and the polymer, with the probability of water forming hydrogen bonds with the carbonyl oxygen of the ester group being about four times higher than with its ether oxygen.

Wide-angle X-ray diffraction and molecular dynamics study of medium-range order in ambient and hot water

We have developed wide-angle X-ray diffraction measurements with high energy-resolution and accuracy to study water structure at three different temperatures (7, 25 and 66 °C) under normal pressure. Using a spherically curved Ge crystal an energy resolution better than 15 eV has been achieved which eliminates influence from Compton scattering. The high quality of the data allows for a reliable Fourier transform of the experimental data resolving shell structure out to ~12 Å, i.e. 5 hydration shells. Large-scale molecular dynamics (MD) simulations using the TIP4P/2005 force-field reproduce excellently the experimental shell-structure in the range 4-12 Å although less agreement is seen for the first peak in the intermolecular pair-correlation function (PCF). The Shiratani-Sasai Local Structure Index [J. Chem. Phys. 104, 7671 (1996)] identifies a tetrahedral minority giving the intermediate-range oscillations in the O-O PCF and a disordered majority providing a more featureless background in this range. The current study supports the proposal that the structure of liquid water, even at high temperatures, can be described in terms of a two-state fluctuation model involving local structures related to the high-density and low-density forms of liquid water postulated in the liquid-liquid phase transition hypothesis.