Oxygen Evolution Reaction Kinetic Barriers on Nitrogen-Doped Carbon Nanotubes.

We investigate kinetic barriers for the oxygen evolution reaction (OER) on singly and doubly nitrogen-doped single-walled carbon nanotubes (NCNTs) using the climbing image nudged elastic band method with solvent effects represented by a 45-water-molecule droplet. The studied sites were chosen based on a previous study of the same systems utilizing a thermodynamic model which ignored both solvent effects and kinetic barriers. According to that model, the two studied sites, one on a singly nitrogen-doped CNT and the other on a doubly doped CNT, were approximately equally suitable for OER. For the four-step OER process, however, our reaction barrier calculations showed a clear difference in the rate-determining *OOH formation step between the two systems, with barrier heights differing by more than 0.4 eV. Thus, the simple thermodynamic model may alone be insufficient for identifying optimal OER sites. Of the remaining three reaction steps, the two H2O forming ones were found to be barrierless in all cases. We also performed solvent-free barrier calculations on NCNTs and undoped CNTs. Substantial differences were observed in the energies of the intermediates when the solvent was present. In general, the observed low activation energy barriers for these reactions corroborate both experimental and theoretical findings of the utility of NCNTs for OER catalysis.

Oxygen Evolution Reaction on Nitrogen-Doped Defective Carbon Nanotubes and Graphene.

The realization of a hydrogen economy would be facilitated by the discovery of a water-splitting electrocatalyst that is efficient, stable under operating conditions, and composed of earth-abundant elements. Density functional theory simulations within a simple thermodynamic model of the more difficult half-reaction, the anodic oxygen evolution reaction (OER), with a single-walled carbon nanotube as a model catalyst, show that the presence of 0.3-1% nitrogen reduces the required OER overpotential significantly compared to the pristine nanotube. We performed an extensive exploration of systems and active sites with various nitrogen functionalities (graphitic, pyridinic, or pyrrolic) obtained by introducing nitrogen and simple lattice defects (atomic substitutions, vacancies, or Stone-Wales rotations). A number of nitrogen functionalities (graphitic, oxidized pyridinic, and Stone-Wales pyrrolic nitrogen systems) yielded similar low overpotentials near the top of the OER volcano predicted by the scaling relation, which was seen to be closely observed by these systems. The OER mechanism considered was the four-step single-site water nucleophilic attack mechanism. In the active systems, the second or third step, the formation of attached oxo or peroxo moieties, was the potential-determining step of the reaction. The nanotube radius and chirality effects were examined by considering OER in the limit of large radius by studying the analogous graphene-based model systems. They exhibited trends similar to those of the nanotube-based systems but often with reduced reactivity due to weaker attachment of the OER intermediate moieties.

Ab initio molecular dynamics studies of formic acid dimer colliding with liquid water.

Ab initio molecular dynamics simulations of formic acid (FA) dimer colliding with liquid water at 300 K have been performed using density functional theory. The two energetically lowest FA dimer isomers were collided with a water slab at thermal and high kinetic energies up to 68kBT. Our simulations agree with recent experimental observations of nearly a complete uptake of gas-phase FA dimer: the calculated average kinetic energy of the dimers immediately after collision is 5 ± 4% of the incoming kinetic energy, which compares well with the experimental value of 10%. Simulations support the experimental observation of no delayed desorption of FA dimers following initial adsorption. Our analysis shows that the FA dimer forms hydrogen bonds with surface water molecules, where the hydrogen bond order depends on the dimer structure, such that the most stable isomer possesses fewer FA-water hydrogen bonds than the higher energy isomer. Nevertheless, even the most stable isomer can attach to the surface through one hydrogen bond despite its reduced hydrophilicity. Our simulations further show that the probability of FA dimer dissociation is increased by high collision energies, the dimer undergoes isomerization from the higher energy to the lowest energy isomer, and concerted double-proton transfer occurs between the FA monomers. Interestingly, proton transfer appears to be driven by the release of energy arising from such isomerization, which stimulates those internal vibrational degrees of freedom that overcome the barrier of a proton transfer.

Deprotonation of formic acid in collisions with a liquid water surface studied by molecular dynamics and metadynamics simulations.

Deprotonation of organic acids at aqueous surfaces has important implications in atmospheric chemistry and other disciplines, yet it is not well-characterized or understood. This article explores the interactions of formic acid (FA), including ionization, in collisions at the air-water interface. Ab initio molecular dynamics simulations with dispersion-corrected density functional theory were used. The 8-50 picosecond duration trajectories all resulted in the adsorption of FA within the interfacial region, with no scattering, absorption into the bulk or desorption into the vapor. Despite the known weak acidity of FA, spontaneous deprotonation of the acid was observed at the interface on a broad picosecond timescale, ranging from a few picoseconds typical for stronger acids to tens of picoseconds. Deprotonation occurred in 4% of the trajectories, and was followed by Grotthuss proton transfer through adjacent water molecules. Both sequential and ultrafast concerted proton transfer were observed. The formation of contact ion pairs and solvent-separated ion pairs, and finally the reformation of neutral FA, both trans and cis conformers, occurred in different stages of the dynamics. To better understand the deprotonation mechanisms at the interface compared with the process in bulk water, we used well-tempered metadynamics to obtain deprotonation free energy profiles. While in bulk water FA deprotonation has a free energy barrier of 14.8 kJ mol-1, in fair agreement with the earlier work, the barrier at the interface is only 7.5 kJ mol-1. Thus, at the air-water interface, FA may dissociate more rapidly than in the bulk. This finding can be rationalized with reference to the dissimilar aqueous solvation and hydrogen-bonding environments in the interface compared to those in bulk liquid water.

Temperature and collision energy effects on dissociation of hydrochloric acid on water surfaces.

Collisions of HCl at the air-water interface modelled by a 72 molecule water slab are studied for a range of various impact energies and temperatures using ab initio molecular dynamics with density functional theory. A range of short-timescale events can follow the collision, from direct scattering to nondissociative trapping on the surface. In most cases, HCl dissociation occurs within a few picoseconds, followed by the formation of a solvent-separated ion pair, or rarely, the reformation of HCl. With increasing impact energy and/or system temperature, dissociation occurs more rapidly, with Cl(-) tending to diffuse deeper into the slab. At temperatures corresponding to the frozen water regime, dissociation is seen only once out of the five thermal collisions, but with the addition of a total of 4kT or more of kinetic energy to HCl, it occurs in all our trajectories within a few ps.

Ab initio water pair potential with flexible monomers.

A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects.

Computational studies of atmospherically-relevant chemical reactions in water clusters and on liquid water and ice surfaces.

CONSPECTUS: Reactions on water and ice surfaces and in other aqueous media are ubiquitous in the atmosphere, but the microscopic mechanisms of most of these processes are as yet unknown. This Account examines recent progress in atomistic simulations of such reactions and the insights provided into mechanisms and interpretation of experiments. Illustrative examples are discussed. The main computational approaches employed are classical trajectory simulations using interaction potentials derived from quantum chemical methods. This comprises both ab initio molecular dynamics (AIMD) and semiempirical molecular dynamics (SEMD), the latter referring to semiempirical quantum chemical methods. Presented examples are as follows: (i) Reaction of the (NO(+))(NO3(-)) ion pair with a water cluster to produce the atmospherically important HONO and HNO3. The simulations show that a cluster with four water molecules describes the reaction. This provides a hydrogen-bonding network supporting the transition state. The reaction is triggered by thermal structural fluctuations, and ultrafast changes in atomic partial charges play a key role. This is an example where a reaction in a small cluster can provide a model for a corresponding bulk process. The results support the proposed mechanism for production of HONO by hydrolysis of NO2 (N2O4). (ii) The reactions of gaseous HCl with N2O4 and N2O5 on liquid water surfaces. Ionization of HCl at the water/air interface is followed by nucleophilic attack of Cl(-) on N2O4 or N2O5. Both reactions proceed by an SN2 mechanism. The products are ClNO and ClNO2, precursors of atmospheric atomic chlorine. Because this mechanism cannot result from a cluster too small for HCl ionization, an extended water film model was simulated. The results explain ClNO formation experiments. Predicted ClNO2 formation is less efficient. (iii) Ionization of acids at ice surfaces. No ionization is found on ideal crystalline surfaces, but the process is efficient on isolated defects where it involves formation of H3O(+)-acid anion contact ion pairs. This behavior is found in simulations of a model of the ice quasi-liquid layer corresponding to large defect concentrations in crystalline ice. The results are in accord with experiments. (iv) Ionization of acids on wet quartz. A monolayer of water on hydroxylated silica is ordered even at room temperature, but the surface lattice constant differs significantly from that of crystalline ice. The ionization processes of HCl and H2SO4 are of high yield and occur in a few picoseconds. The results are in accord with experimental spectroscopy. (v) Photochemical reactions on water and ice. These simulations require excited state quantum chemical methods. The electronic absorption spectrum of methyl hydroperoxide adsorbed on a large ice cluster is strongly blue-shifted relative to the isolated molecule. The measured and calculated adsorption band low-frequency tails are in agreement. A simple model of photodynamics assumes prompt electronic relaxation of the excited peroxide due to the ice surface. SEMD simulations support this, with the important finding that the photochemistry takes place mainly on the ground state. In conclusion, dynamics simulations using quantum chemical potentials are a useful tool in atmospheric chemistry of water media, capable of comparison with experiment.

First and second deprotonation of H2SO4 on wet hydroxylated (0001) α-quartz.

We present an ab initio molecular dynamics study of deprotonation of sulfuric acid on wet quartz, a topic of atmospheric interest. The process is preferred, with 65% of our trajectories at 250 K showing deprotonation. The time distribution of the deprotonation events shows an exponential behavior and predicts an average deprotonation time of a few picoseconds. The process is exoergic, with most of the temperature increase being due to formation of hydrogen bonds prior to deprotonation. In agreement with existing studies of H2SO4 in water clusters, in liquid water, and at the air-water interface, the main determinant of deprotonation is the degree of solvation of H2SO4 by neighboring water molecules. However, we find that if both hydrogens of H2SO4 are simultaneously donated to water oxygens, deprotonation is disfavored. Predicted spectroscopic signatures showing the presence of solvated hydronium and bisulfate are presented. Increasing the temperature up to 330 K accelerates the process but does not change the main features of the deprotonation mechanisms or the spectroscopic signatures. The second deprotonation of H2SO4, studied only at 250 K, occurs provided there is sufficient solvation of the bisulfate by additional water molecules. In comparison to HCl deprotonation on the identical surface examined in our previous work, the first deprotonation of H2SO4 occurs more readily and releases more energy.

Raman spectroscopy of solutions and interfaces containing nitrogen dioxide, water, and 1,4 dioxane: evidence for repulsion of surface water by NO2 gas.

The interaction of water, 1,4 dioxane, and gaseous nitrogen dioxide, has been studied as a function of distance measured through the liquid-vapour interface by Raman spectroscopy with a narrow (<0.1 mm) laser beam directed parallel to the interface. The Raman spectra show that water is present at the surface of a dioxane-water mixture when gaseous NO2 is absent, but is virtually absent from the surface of a dioxane-water mixture when gaseous NO2 is present. This is consistent with recent theoretical calculations that show NO2 to be mildly hydrophobic.

Nitrogen dioxide at the air-water interface: trapping, absorption, and solvation in the bulk and at the surface.

The interaction of NO(2) with water surfaces in the troposphere is of major interest in atmospheric chemistry. We examined an initial step in this process, the uptake of NO(2) by water through the use of molecular dynamics simulations. An NO(2)-H(2)O intermolecular potential was obtained by fitting to high-level ab initio calculations. We determined the binding of NO(2)-H(2)O to be about two times stronger than that previously calculated. From scattering simulations of an NO(2) molecule interacting with a water slab we observed that the majority of the scattering events resulted in outcomes in which the NO(2) molecule became trapped at the surface or in the interior of the water slab. Typical surface-trapped/adsorbed and bulk-solvated/absorbed trajectories were analyzed to obtain radial distribution functions and the orientational propensity of NO(2) with respect to the water surface. We observed an affinity of the nitrogen atom for the oxygen in water, rather than hydrogen-bonding which was rare. The water solvation shell was less tight for the bulk-absorbed NO(2) than for the surface-adsorbed NO(2). Adsorbed NO(2) demonstrated a marked orientational preference, with the oxygens pointing into the vacuum. Such behavior is expected for a mildly hydrophobic and surfactant molecule like NO(2). Estimates based on our results suggest that at high NO(2) concentrations encountered, for example, in some sampling systems, adsorption and reaction of NO(2) at the surface may contribute to the formation of gas-phase HONO.

Study of ion specific interactions of alkali cations with dicarboxylate dianions.

Alkali metal cations often show pronounced ion-specific interactions and selectivity with macromolecules in biological processes, colloids, and interfacial sciences, but a fundamental understanding about the underlying microscopic mechanism is still very limited. Here we report a direct probe of interactions between alkali metal cations (M(+)) and dicarboxylate dianions, (-)O(2)C(CH(2))(n)CO(2)(-) (D(n)(2-)) in the gas phase by combined photoelectron spectroscopy (PES) and ab initio electronic structure calculations on nine M(+)-D(n)(2-) complexes (M = Li, Na, K; n = 2, 4, 6). PES spectra show that the electron binding energy (EBE) decreases from Li(+) to Na(+) to K(+) for complexes of M(+)-D(2)(2-), whereas the order is Li(+) < Na(+) ≈ K(+) when M(+) interacts with a more flexible D(6)(2-) dianion. Theoretical modeling suggests that M(+) prefers to interact with both ends of the carboxylate -COO(-) groups by bending the flexible aliphatic backbone, and the local binding environments are found to depend upon backbone length n, carboxylate orientation, and the specific cation M(+). The observed variance of EBEs reflects how well each specific dicarboxylate dianion accommodates each M(+). This work demonstrates the delicate interplay among several factors (electrostatic interaction, size matching, and strain energy) that play critical roles in determining the structures and energetics of gaseous clusters as well as ion specificity and selectivity in solutions and biological systems.

Semiempirical self-consistent polarization description of bulk water, the liquid-vapor interface, and cubic ice.

We have applied an efficient electronic structure approach, the semiempirical self-consistent polarization neglect of diatomic differential overlap (SCP-NDDO) method, previously parametrized to reproduce properties of water clusters by Chang, Schenter, and Garrett [ J. Chem. Phys. 2008 , 128 , 164111 ] and now implemented in the CP2K package, to model ambient liquid water at 300 K (both the bulk and the liquid-vapor interface) and cubic ice at 15 and 250 K. The SCP-NDDO potential retains its transferability and good performance across the full range of conditions encountered in the clusters and the bulk phases of water. In particular, we obtain good results for the density, radial distribution functions, enthalpy of vaporization, self-diffusion coefficient, molecular dipole moment distribution, and hydrogen bond populations, in comparison to experimental measurements.

Improving the density functional theory description of water with self-consistent polarization.

We applied the self-consistent polarization density functional theory (SCP-DFT) to water. SCP-DFT requires only minimal parametrization, self-consistently includes the dispersion interaction neglected by standard DFT functionals, and has a cost similar to standard DFT despite its improved performance. Compared to the DFT functionals BLYP and BLYP-D (where the latter contains a simple dispersion correction), SCP-DFT yields interaction energies per molecule and harmonic frequencies of clusters in better agreement with experiment, with errors in the former of only a few tenths of a kcal/mol. BLYP and BLYP-D underbind and overbind the clusters, respectively, by up to about 1 kcal/mol. For liquid water, both BLYP and SCP-DFT predict radial distribution functions that are similar and overstructured compared to experiment. However, SCP-DFT improves over BLYP in predicting the experimental enthalpy of vaporization. A decomposition of the dimer interaction energy attempts to rationalize the performance of SCP-DFT. The SCP-DFT approach holds promise as an efficient and accurate method for describing large hydrogen-bonded systems, and has the potential to model complex systems with minimal parametrization.

Self-consistent polarization density functional theory: application to argon.

We present a comprehensive set of results for argon, a case study in weak interactions, using the self-consistent polarization density functional theory (SCP-DFT). With minimal parametrization, SCP-DFT is found to give excellent results for the dimer interaction energy, the second virial coefficient, the liquid structure, and the lattice constant and cohesion energy of the face-centered cubic crystal compared to both accurate theoretical and experimental benchmarks. Thus, SCP-DFT holds promise as a fast, efficient, and accurate method for performing ab initio dynamics that include additional polarization and dispersion interactions for large, complex systems involving solvation and bond breaking.

Toward an accurate and efficient theory of physisorption. I. Development of an augmented density-functional theory model.

Currently available density functionals cannot describe the dispersion component of the interaction energy present in weakly bound complexes. Moreover, the exchange energy as obtained from the density-functional theory is often incorrect. Examples of problematic cases include clusters of van der Waals-bound rare-gas atoms and most hydrogen-bonded molecular systems. Thus, accurate ab initio methods to treat intermolecular forces should be used in such systems. These methods are, however, too slow to be applicable to the large systems needed to model adsorption. This is why DFT continues to be used, where, in addition, a quite common compensation of errors sometimes produces some sort of agreement with the corresponding experimental data. In this paper, we analyze in detail the inadequacy of standard DFT for describing the weak binding present in a few rare gas-rare gas, metal atom-rare gas, and metal atom-metal atom dimers.Inspired by the success of the Hartree-Fock plus (damped) dispersion (HFD) method, we test the use of an improved hybrid model in which to a density-functional interaction energy (with corrected exchange and avoidance of double-counting of dispersion), a (damped) dispersion expansion is added in the usual way.Comparisons with accurate theoretical or experimental benchmarks show that our DFdD method using the revPBEx or revPBEx+VWNc functionals and accurate dispersion coefficients is found to recover the interaction energy curves very well for many of the tested systems. The sec and paper in this series will describe the use of the DFdD method for physisorption for the previously well-studied (but not solved) case of Xe/Cu(111).

(HCl)2 and (HF)2 in small helium clusters: quantum solvation of hydrogen-bonded dimers.

We present a rigorous theoretical study of the solvation of (HCl)(2) and (HF)(2) by small ((4)He)(n) clusters, with n=1-14 and 30. Pairwise-additive potential-energy surfaces of He(n)(HX)(2) (X=Cl and F) clusters are constructed from highly accurate four-dimensional (rigid monomer) HX-HX and two-dimensional (rigid monomer) He-HX potentials and a one-dimensional He-He potential. The minimum-energy geometries of these clusters, for n=1-6 in the case of (HCl)(2) and n=1-5 for (HF)(2), correspond to the He atoms in a ring perpendicular to and bisecting the HX-HX axis. The quantum-mechanical ground-state energies and vibrationally averaged structures of He(n)(HCl)(2) (n=1-14 and 30) and He(n)(HF)(2) (n=1-10) clusters are calculated exactly using the diffusion Monte Carlo (DMC) method. In addition, the interchange-tunneling splittings of He(n)(HCl)(2) clusters with n=1-14 are determined using the fixed-node DMC approach, which was employed by us previously to calculate the tunneling splittings for He(n)(HF)(2) clusters, n=1-10 [A. Sarsa et al., Phys. Rev. Lett. 88, 123401 (2002)]. The vibrationally averaged structures of He(n)(HX)(2) clusters with n=1-6 for (HCl)(2) and n=1-5 for (HF)(2) have the helium density localized in an effectively one-dimensional ring, or doughnut, perpendicular to and at the midpoint of the HX-HX axis. The rigidity of the solvent ring varies with n and reaches its maximum for the cluster size at which the ring is filled, n=6 and n=5 for (HCl)(2) and (HF)(2), respectively. Once the equatorial ring is full, the helium density spreads along the HX-HX axis, eventually solvating the entire HX dimer. The interchange-tunneling splitting of He(n)(HCl)(2) clusters hardly varies at all over the cluster size range considered, n=1-14, and is virtually identical to that of the free HCl dimer. This absence of the solvent effect is in sharp contrast with our earlier results for He(n)(HF)(2) clusters, which show a approximately 30% reduction of the tunneling splitting for n=4. A tentative explanation for this difference is proposed. The implications of our results for the interchange-tunneling dynamics of (HCl)(2) in helium nanodroplets are discussed.

Intermolecular potential energy surface and spectra of He-HCl with generalization to other rare gas-hydrogen halide complexes.

A two-dimensional (rigid monomer) intermolecular potential energy surface (PES) of the He-HCl complex has been obtained from ab initio calculations utilizing the symmetry-adapted perturbation theory (SAPT) and an spdfg basis set including midbond functions. The bond length in HCl was chosen to be equal to the expectation value in the ground vibrational state of isolated HCl. The rigid-monomer potential should be a very good approximation to the complete (three-dimensional) potential for H-Cl distances corresponding to the lowest vibrational levels of the monomer since the He-HCl interaction energy was found to be only weakly dependent on the HCl bond length in this region, at least as compared to systems such as Ar-HF. The calculated points were fitted using an analytic function with ab initio computed asymptotic coefficients. As expected, the complex is loosely bound, with the dispersion energy providing the majority of the attraction. Our SAPT PES agrees with the semiempirical PES of Willey et al. [J. Chem. Phys. 96, 898 (1992)], in finding that, atypically for rare gas-hydrogen halide complexes including the lighter halide atoms, the global minimum is on the Cl side (with intermonomer separation 3.35 A and depth of 32.8 cm(-1)), rather than on the H side, where there is only a local minimum (3.85 A, 30.8 cm(-1)). The ordering of the minima was confirmed by single-point calculations in larger basis sets and complete basis set extrapolations, and also using higher levels of theory. We show that the opposite findings in the recent calculations of Zhang and Shi [J. Mol. Struct: THEOCHEM 589, 89 (2002)] are due to the lack of midbond functions in their basis set. Despite the closeness in depth of the two linear minima, the existence of a relatively high barrier between them invalidates the assumption of isotropy, a feature of some literature potentials. The trends concerning the locations of minima within the family of rare gas-hydrogen halide complexes are rationalized in terms of the physical components of the intermolecular forces and related to monomer properties. The accuracy of the SAPT PES was tested by performing calculations of rovibrational levels. The transition frequencies obtained were found to be in excellent agreement (to within 0.02 cm(-1)) with the measurements of Lovejoy and Nesbitt [J. Chem. Phys. 93, 5387 (1990)]. The SAPT PES predicts a dissociation energy for the complex of 7.74 cm(-1) which is probably more accurate than the experimental value of 10.1+/-1.2 cm(-1). Our analysis of the ground-state rovibrational wave function shows that the He-HCl configuration is favored over the He-ClH configuration despite the ordering of minima. This is due to the greater volume of the well in the former case. We have also determined positions and widths of three low-lying resonance states through scattering calculations. These predictions are expected to be more accurate than values derived from experiment.

Efficient generation of flexible-monomer intermolecular potential energy surfaces.

A new method of generating flexible-monomer intermolecular interaction potentials has been proposed. The method, based on symmetry-adapted perturbation theory, extends a rigid-monomer potential into a flexible-monomer one at a cost negligible compared to performing calculations on a full-dimensional grid (i.e., including internal degrees of freedom of monomers). The non-rigidity effects are accounted for by density-overlap integrals and by asymptotic expansion coefficients. Results for a model system (Ar-HF) demonstrate that the method recovers a substantial portion of these effects.